Tổng quan nghiên cứu

Chi Boehmeria thuộc họ Tầm ma (Urticaceae) hiện diện trên toàn cầu với quy mô khoảng 100 loài, phân bố chủ yếu tại các khu vực nhiệt đới và cận nhiệt đới. Trong đó, cây Gai toàn tơ (Boehmeria holosericea Blume), hay còn gọi là cây gai hoa bạc, là loài cây thân thảo hoặc cây bụi lâu năm với chiều cao trung bình từ 1,5 m đến 2,0 m. Loài thực vật này sở hữu giá trị kép nổi bật: vừa cung cấp sợi tự nhiên bền chắc giàu cellulose cho ngành dệt may sinh thái, vừa là nguồn dược liệu quý trong y học cổ truyền dùng để thanh nhiệt, tiêu viêm, chữa cảm sốt, viêm thận, phù thũng và an thai. Tuy nhiên, rào cản lớn nhất trong thực tiễn là cây Gai toàn tơ có đặc điểm hình thái bên ngoài rất tương đồng với nhiều loài khác trong cùng chi. Khi dược liệu được thu hái, phơi khô, thái lát hoặc nghiền thành dạng bột mịn, việc phân loại bằng mắt thường gần như bất khả thi, dẫn đến nguy cơ nhầm lẫn loài và gian lận thương mại nghiêm trọng.

Xuất phát từ yêu cầu cấp thiết đó, công trình nghiên cứu của tác giả Nguyễn Văn Trình dưới sự hướng dẫn khoa học của Giáo sư Chu Hoàng Mậu tại Trường Đại học Sư phạm - Đại học Thái Nguyên đã tập trung phân tích toàn diện cấu trúc hình thái, vi phẫu giải phẫu thân, rễ, lá; đồng thời phân lập và giải trình tự 5 vùng gene mục tiêu gồm ITS, matK, trnL-trnF, trnH-psbA và rpoC1. Mẫu thực vật mang ký hiệu BOE202405SL được thu thập trực tiếp vào tháng 5 năm 2024 tại bản Áng Ung, xã Chiềng Ban, huyện Mai Sơn, tỉnh Sơn La (tọa độ 21°13'49"N). Đề tài mang ý nghĩa học thuật và thực tiễn to lớn khi bổ sung bộ dữ liệu trình tự phân tử chuẩn mực lên hệ thống Mendeley Data và ngân hàng gen quốc tế GenBank, mở ra giải pháp định danh loài chính xác với hệ số tin cậy thống kê bootstrap đạt từ 94% đến 100%.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết mã vạch DNA (DNA barcoding) do Liên đoàn Mã vạch Sự sống (CBOL Plant Working Group) khởi xướng. Theo các tiêu chuẩn phân tử quốc tế, một mã vạch DNA thực vật lý tưởng phải đáp ứng kích thước chuẩn từ 400 bp đến 800 bp, dễ dàng khuếch đại bằng kỹ thuật PCR từ mô thực vật thô, có tỷ lệ biến đổi nucleotide đủ cao để phân biệt các loài thân thuộc nhưng đủ bảo tồn để nhận diện cấp độ chi và họ. Hai vùng mã vạch thực vật cốt lõi được ứng dụng chủ đạo là vùng đệm phiên mã nội nhân ITS (Internal Transcribed Spacer) thuộc hệ gene nhân và gene matK (maturase K) thuộc hệ gene lục lạp.

Bên cạnh đó, nghiên cứu vận dụng lý thuyết tiến hóa phân tử và khoảng cách di truyền thông qua mô hình thay thế nucleotide Tamura-Nei kết hợp thuật toán Maximum Likelihood (ML) với 1.000 lần lặp lại bootstrap để xây dựng cây phát sinh chủng loại. Hệ gene lục lạp thực vật di truyền theo dòng mẹ, sở hữu cấu trúc bốn phần đặc trưng gồm vùng sao chép đơn lớn LSC, vùng sao chép đơn nhỏ SSC và hai vùng lặp lại ngược IR. Điển hình ở cây Gai Ramie (Boehmeria nivea), bộ gene lục lạp có chiều dài 156.065 bp mang 131 gene với hàm lượng GC đạt 36,33%; trong khi loài Boehmeria umbrosa và Boehmeria spicata có kích thước bộ gene lục lạp lần lượt là 170.920 bp và 170.958 bp. Các khái niệm phân loại học cốt lõi như nhánh đơn ngành (monophyletic) và nhánh cận đơn ngành (paraphyletic) được sử dụng để xác định vị trí phân loại chính xác của cây Gai toàn tơ.

Phương pháp nghiên cứu

Nguồn dữ liệu thực nghiệm bao gồm mẫu tươi thu thập tại thực địa Sơn La kết hợp với 100 trình tự nucleotide tương đồng cao nhất được truy xuất từ cơ sở dữ liệu NCBI GenBank. Phương pháp lấy mẫu thực địa có chủ đích được áp dụng tại sinh cảnh bìa rừng thứ sinh ven suối, chọn lọc các cá thể trưởng thành sinh trưởng tốt, không sâu bệnh và lưu giữ mẫu tiêu bản chuẩn tại Khoa Sinh học, Trường Đại học Sư phạm - Đại học Thái Nguyên.

Phương pháp giải phẫu thực vật sử dụng kỹ thuật cắt lát ngang tự do bằng tay đối với rễ, thân, cuống lá, gân lá và phiến lá, xử lý làm sạch nội chất bằng dung dịch Javen trong 30 phút, ngâm acid acetic 1% trong 10 phút để tăng lực bắt màu, sau đó nhuộm kép bằng dung dịch phẩm màu Carmin (bắt màu đỏ tím trên vách cellulose) và xanh Methylen (bắt màu xanh lục trên vách hóa gỗ lignin). Tiêu bản được quan sát và chụp ảnh vi thể trên kính hiển vi quang học Kruss ở độ phóng đại 10X và 40X.

Đối với phân tích phân tử, DNA tổng số được tách chiết từ 100 mg lá non bằng dung dịch đệm CTAB 2X tại 65 độ C nhằm loại trừ tạp chất polysaccharide và polyphenol. Phản ứng PCR thể tích 50 µL được vận hành trên máy Biorad T100 để khuếch đại 5 phân đoạn gene bằng các cặp mồi chuyên biệt. Sản phẩm PCR tinh sạch từ gel agarose được giải trình tự Sanger tự động trên hệ thống Applied Biosystems 3500xL Genetic Analyzer. Trình tự nucleotide được hiệu chỉnh bằng BioEdit và phân tích phát sinh chủng loại qua phần mềm MEGA 11 trong khoảng thời gian từ tháng 5 năm 2024 đến tháng 5 năm 2025.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, kết quả phân tích giải phẫu vi phẫu chứng minh cây Gai toàn tơ mang cấu trúc điển hình của thực vật hai lá mầm, với dấu hiệu đặc trưng nổi bật là sự phân bố dày đặc của các tinh thể canxi oxalat (calcium oxalate) hình cầu gai và khối đa diện trong tất cả các cơ quan rễ, thân, lá. Tại thân cây, mô mềm ruột phát triển rất rộng, chiếm tới 2/3 bán kính vi phẫu ngang; rễ cây sở hữu trụ bì gồm 2 đến 3 lớp tế bào hóa sợi cứng tập trung thành chùy; phiến lá có mô giậu xếp sít nhau cùng 6 đến 8 cặp gân lá nổi rõ ở mặt dưới.

Thứ hai, nghiên cứu đã phân lập và đọc trình tự thành công 5 phân đoạn DNA từ hệ gen cây Gai toàn tơ với kích thước cụ thể: vùng ITS dài 659 bp; gene matK dài 1.219 bp; vùng đệm liên gene trnL-trnF dài 840 bp; gene rpoC1 dài 523 bp; và vùng đệm liên gene trnH-psbA dài 427 bp. Trong đó, gene matK có chiều dài nucleotide lớn nhất (1.219 bp) và vùng đệm trnH-psbA có kích thước ngắn nhất (427 bp), toàn bộ đều nằm trong giới hạn tối ưu cho các phản ứng khuếch đại phân tử chuẩn.

Thứ ba, so sánh đối chiếu với các trình tự tương đồng cao nhất trên GenBank đã phát hiện các điểm đột biến nucleotide độc bản. Vùng ITS của Gai toàn tơ có 4 vị trí sai khác tại nucleotide thứ 21 (T chuyển thành G), 22 (G chuyển thành T), 63 (A chuyển thành G) và 556 (C chuyển thành A) khi so sánh với loài Boehmeria clidemioides var. diffusa. Gene matK chỉ xuất hiện duy nhất 1 điểm đột biến thay thế tại nucleotide thứ 234 (A chuyển thành G) so với loài Boehmeria macrophylla. Gene rpoC1 có 2 vị trí sai khác tại nucleotide 20 và 26 (đều là G chuyển thành A) so với loài Boehmeria densiflora; trong khi phân đoạn trnL-trnF hoàn toàn tương đồng 100% với Boehmeria macrophylla.

Thứ tư, sơ đồ cây phát sinh chủng loại xây dựng từ gene matK thể hiện Gai toàn tơ tách biệt thành một nhóm đơn ngành rõ rệt với độ tin cậy bootstrap đạt mức tuyệt đối 100%. Vùng đệm trnH-psbA phân bố loài nghiên cứu chung nhánh với Boehmeria clidemioides với chỉ số bootstrap 96% và Boehmeria umbrosa với bootstrap 94%. Gene rpoC1 liên kết Gai toàn tơ cùng nhóm với Boehmeria spicata với độ ủng hộ bootstrap 100%.

Thảo luận kết quả

Mật độ dày đặc của tinh thể canxi oxalat trong mô mềm và thịt lá phản ánh cơ chế tích lũy khoáng chất thích nghi sinh thái tại các vùng núi đá có độ ẩm cao, đóng vai trò như một đặc điểm vi học định lượng quý giá để nhận dạng nhanh bột dược liệu dưới kính hiển vi. Dữ liệu phân tử của 5 vùng gene có thể được trực quan hóa sinh động thông qua biểu đồ cột so sánh tỷ lệ nucleotide tương đồng và bảng ma trận khoảng cách di truyền, kết hợp với các sơ đồ cây phát sinh chủng loại cladogram thể hiện hệ số bootstrap sau 1.000 lần lặp.

Sự hình thành các nhánh đơn ngành vững chắc với giá trị bootstrap đạt từ 94% đến 100% trên cây chủng loại matK, rpoC1 và trnH-psbA chứng tỏ các vùng gene này có tốc độ tiến hóa phân tử thích hợp, vừa đảm bảo tính bảo tồn ở cấp độ chi, vừa sở hữu các đột biến điểm đặc hiệu loài. Do đó, việc kết hợp gene matK cùng vùng đệm trnH-psbA cung cấp giải pháp nhận dạng phân tử tối ưu, giải quyết triệt để tình trạng nhầm lẫn hình thái giữa các loài trong chi Boehmeria khi lưu thông trên thị trường dược liệu.

Đề xuất và khuyến nghị

Thứ nhất, chuẩn hóa quy trình kiểm nghiệm dược liệu bằng kỹ thuật multiplex PCR kết hợp cặp chỉ thị matK và trnH-psbA. Mục tiêu nhằm rút ngắn thời gian định danh loài xuống dưới 24 giờ với độ chính xác đạt 100% đối với các mẫu dược liệu thô hoặc dạng bột nghiền mịn. Cơ quan chịu trách nhiệm thực hiện là Viện Dược liệu Trung ương và hệ thống các Trung tâm Kiểm nghiệm Thuốc trên toàn quốc, triển khai áp dụng trong giai đoạn từ năm 2025 đến năm 2026.

Thứ hai, mở rộng phạm vi khảo sát thực địa và tăng quy mô lấy mẫu lên trên 60 cá thể cây Gai toàn tơ tại các vùng sinh thái đặc trưng như Lào Cai, Hòa Bình, Quảng Ninh và khu vực Tây Nguyên. Mục tiêu nhằm đánh giá toàn diện độ biến thiên nucleotide nội loài và tính đa dạng di truyền với độ tin cậy thống kê trên 99%. Đơn vị thực hiện là các nhóm nghiên cứu Di truyền học phối hợp cùng Chi cục Kiểm lâm các địa phương trong giai đoạn 2026 - 2027.

Thứ ba, triển khai giải trình tự toàn bộ hệ gene lục lạp (Complete Chloroplast Genome) của loài Boehmeria holosericea với kích thước dự kiến khoảng 156.000 bp đến 171.000 bp. Mục tiêu nhằm xác định đầy đủ các trình tự lặp lại vi vệ tinh SSR và tái cấu trúc cây phát sinh chủng loại toàn diện cho chi Boehmeria. Nhiệm vụ này do các phòng thí nghiệm trọng điểm thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam chủ trì hoàn thành trước năm 2028.

Thứ tư, xây dựng mô hình nhân giống vô tính sạch bệnh và bảo tồn nguồn gen Gai toàn tơ phục vụ phát triển vùng nguyên liệu công nghiệp. Mục tiêu cung cấp 20.000 cây giống đạt chuẩn chất lượng mỗi năm, hướng tới chiết xuất các hợp chất tự nhiên quý như ruspolinone và loliolide phục vụ dệt may sinh học và dược phẩm chức năng. Đơn vị chủ trì là Sở Khoa học và Công nghệ tỉnh Sơn La phối hợp cùng các hợp tác xã nông nghiệp công nghệ cao từ năm 2025 đến năm 2030.

Đối tượng nên tham khảo luận văn

Thứ nhất, học viên cao học, nghiên cứu sinh và các giảng viên đại học thuộc chuyên ngành Di truyền học, Thực vật học và Công nghệ sinh học. Luận văn cung cấp case study hoàn chỉnh về việc tích hợp kỹ thuật hình thái vi phẫu với giải trình tự 5 vùng DNA, đồng thời hướng dẫn chi tiết cách vận hành phần mềm MEGA 11 để dựng cây phát sinh phân tử theo thuật toán Maximum Likelihood với 1.000 lần lặp.

Thứ hai, kỹ thuật viên và chuyên viên kiểm định chất lượng tại các Viện Kiểm nghiệm Thuốc, Trung tâm Kiểm định Dược phẩm và Mỹ phẩm. Nghiên cứu cung cấp thông số kích thước chuẩn của 5 phân đoạn DNA (matK 1.219 bp, trnL-trnF 840 bp, ITS 659 bp, rpoC1 523 bp, trnH-psbA 427 bp) và đặc điểm tinh thể canxi oxalat để làm quy chuẩn đối sánh, loại bỏ dược liệu giả mạo.

Thứ ba, các doanh nghiệp chế biến dược phẩm thiên nhiên, sản xuất thực phẩm bảo vệ sức khỏe và các đơn vị dệt may sinh thái. Dữ liệu luận văn giúp doanh nghiệp lựa chọn đúng nguồn gen nguyên liệu chuẩn, kiểm soát chất lượng đầu vào của sợi tự nhiên và các hoạt chất chống oxy hóa, kháng viêm đạt tiêu chuẩn GACP-WHO.

Thứ tư, cán bộ quản lý tài nguyên thiên nhiên tại các Ban quản lý Vườn quốc gia, Khu bảo tồn thiên nhiên và Chi cục Kiểm lâm. Tài liệu là công cụ hỗ trợ hữu ích trong công tác phân loại thực vật, giám sát quần thể loài bản địa và xây dựng cơ sở dữ liệu số hóa tài nguyên thực vật rừng Việt Nam.

Câu hỏi thường gặp

Cây Gai toàn tơ (Boehmeria holosericea) khác biệt gì so với cây Gai làm bánh (Boehmeria nivea)? Về mặt hình thái, cây Gai toàn tơ có phiến lá mọc đối hình mác với 6 đến 8 cặp gân nổi rõ, toàn thân phủ lớp lông tơ màu trắng xám mịn màng ánh bạc. Về mặt phân tử, Gai toàn tơ có trình tự 5 vùng gene đặc trưng với các đột biến nucleotide riêng biệt, khác biệt hoàn toàn với hệ gen lục lạp dài 156.065 bp của cây Gai làm bánh.

Tại sao nghiên cứu phải kết hợp cả phương pháp hình thái vi phẫu và chỉ thị phân tử DNA? Trong thực tế, hình thái bên ngoài của các loài chi Boehmeria rất giống nhau và dễ bị biến đổi hoàn toàn khi phơi khô, nghiền bột. Việc kết hợp quan sát tinh thể canxi oxalat qua vi phẫu và phân tích các đoạn DNA mã vạch giúp quá trình định danh đạt độ chính xác 100%, khắc phục nhược điểm của từng phương pháp đơn lẻ.

Vùng gene nào trong số 5 vùng nghiên cứu có tiềm năng làm mã vạch DNA tốt nhất cho chi Boehmeria? Gene matK (1.219 bp) và vùng đệm liên gene trnH-psbA (427 bp) là hai ứng viên xuất sắc nhất. Trên cây phát sinh chủng loại, gene matK tạo nên nhánh đơn ngành tách biệt với hệ số bootstrap 100%, trong khi trnH-psbA có độ đa dạng nucleotide cao với hệ số bootstrap hỗ trợ lên tới 96%, rất phù hợp để phân biệt các loài cận duyên.

Đặc điểm giải phẫu vi phẫu nào của cây Gai toàn tơ có giá trị chẩn đoán cao nhất? Đặc điểm nhận diện vi học giá trị nhất là sự xuất hiện với mật độ dày đặc của các tinh thể canxi oxalat hình cầu gai ở phiến lá và hình khối đa diện rải rác trong mô mềm vỏ, libe và ruột của cả thân, rễ và cuống lá, kết hợp với vùng mô mềm ruột chiếm tới 2/3 bán kính vi phẫu thân.

Dữ liệu trình tự gene của cây Gai toàn tơ trong luận văn có thể tra cứu ở đâu? Toàn bộ tập dữ liệu trình tự nucleotide của 5 vùng gene ITS, matK, trnL-trnF, trnH-psbA và rpoC1 phân lập từ mẫu BOE202405SL đã được công bố chính thức trên kho dữ liệu khoa học quốc tế Mendeley Data và đăng ký trên ngân hàng gen NCBI GenBank, phục vụ việc tra cứu và đối sánh phân tử trên toàn thế giới.

Kết luận

  • Luận văn đã mô tả toàn diện đặc điểm hình thái và giải phẫu vi phẫu rễ, thân, lá của cây Gai toàn tơ, xác lập dấu hiệu chẩn đoán đặc trưng là mật độ dày đặc của các tinh thể canxi oxalat trong toàn bộ các mô cơ quan.
  • Phân lập, nhân bản và giải trình tự thành công 5 phân đoạn DNA chuẩn với kích thước cụ thể: ITS đạt 659 bp, matK đạt 1.219 bp, trnL-trnF đạt 840 bp, rpoC1 đạt 523 bp và trnH-psbA đạt 427 bp.
  • Xác định chính xác các điểm đột biến thay thế nucleotide đặc hiệu loài trên vùng ITS (4 vị trí), gene matK (1 vị trí) và gene rpoC1 (2 vị trí) so với các loài họ hàng gần nhất trên GenBank.
  • Xây dựng thành công cây phát sinh chủng loại phân tử với thuật toán Maximum Likelihood, khẳng định gene matK, rpoC1 và vùng đệm trnH-psbA là những ứng viên mã vạch DNA tiềm năng cao với độ tin cậy bootstrap lên tới 100%.
  • Công bố toàn bộ tập dữ liệu phân tử mới trên cơ sở dữ liệu quốc tế Mendeley Data, đóng góp nguồn tư liệu chuẩn mực cho phân loại học chi Boehmeria tại Việt Nam.

Lộ trình nghiên cứu tiếp theo sẽ tập trung giải mã toàn bộ hệ gene lục lạp và đánh giá hoạt tính sinh học của các hợp chất chiết xuất trong giai đoạn 2025 - 2027. Các doanh nghiệp và đơn vị kiểm nghiệm dược liệu cần chủ động liên kết với nhóm tác giả để ứng dụng ngay bộ chỉ thị phân tử này vào quy trình kiểm định chất lượng giống cây trồng và dược phẩm thương mại.