Tổng quan nghiên cứu

Thế giới hiện ghi nhận khoảng 350.000 loài thực vật đã biết, trong đó có hơn 70.000 loài mang hoạt tính dược lý giá trị phục vụ y học cổ truyền và công nghiệp dược phẩm hiện đại. Tuy nhiên, áp lực từ biến đổi khí hậu cùng các hoạt động khai thác quá mức đang đẩy hệ sinh thái vào khủng hoảng suy thoái, với tốc độ tuyệt chủng ước tính từ 1 đến 50 loài mỗi ngày và dự báo khoảng 25% tổng số loài thực vật có thể biến mất trước năm 2050 nếu thiếu các biện pháp bảo tồn đồng bộ. Tại Việt Nam, chi Cẩm lai (Dalbergia) thuộc họ Đậu (Fabaceae) là một trong những nhóm thực vật đặc biệt quý hiếm, sở hữu các hợp chất thứ cấp quan trọng như isoflavonoid, neoflavonoid có khả năng chống dị ứng và ức chế tế bào ung thư. Dù vậy, do giá trị kinh tế vượt trội của phần gỗ lõi bền chắc và hương thơm đặc trưng, nhiều loài Dalbergia đã bị khai thác kiệt quệ đến mức Nguy cấp (EN) theo Danh lục Đỏ.

Thực trạng phân loại truyền thống dựa trên cơ quan sinh dưỡng và sinh sản gặp vô vàn rào cản khi mẫu vật lưu thông trên thị trường thường đã qua sơ chế cắt nhỏ, phơi khô hoặc ở dạng bột. Để giải quyết triệt để vấn đề định danh và bảo tồn, nghiên cứu được thực hiện với mục tiêu phân tích toàn diện các đặc điểm thực vật học kết hợp giải mã trình tự các vùng gen chỉ thị matK, rbcL và ITS2 của loài Dalbergia sp., thường gọi là cây Cựa gà. Mẫu nghiên cứu được thu thập tại xã Hợp Hòa, huyện Sơn Dương, tỉnh Tuyên Quang (tọa độ địa lý 21°37′24″B và 105°26′25″Đ). Đề tài mang ý nghĩa then chốt trong việc xây dựng cơ sở dữ liệu mã vạch DNA quốc gia, cung cấp công cụ giám định phân tử có độ chính xác đạt trên 96% đến 99%, phục vụ trực tiếp cho công tác bảo tồn nguồn gen bản địa và quản lý chất lượng dược liệu.

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 dựa trên lý thuyết mã vạch DNA (DNA Barcoding) do Liên minh Mã vạch Sự sống (CBOL) đề xuất, xem các đoạn trình tự nucleotide tiêu chuẩn ngắn như mã định danh phân tử đặc hiệu cho từng loài sinh vật. Khung lý thuyết kết hợp song song hai hệ thống gen: các gen lục lạp bảo thủ gồm maturase K (matK) và ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit (rbcL), cùng với vùng đệm phiên mã nội bộ trong nhân tế bào (Internal Transcribed Spacer - ITS2 nằm giữa đoạn 5.8S và 28S rRNA).

Mô hình nghiên cứu vận dụng lý thuyết tiến hóa phân tử Tamura-Nei và thuật toán suy luận quan hệ phát sinh chủng loại Maximum Likelihood (ML). Các khái niệm cốt lõi được vận hành trong nghiên cứu bao gồm: đột biến thay thế nucleotide (đồng chuyển và đảo chuyển), hệ số tương đồng di truyền, và giá trị độ tin cậy bootstrap với 1000 lần lặp lại nhằm kiểm chứng mức độ vững chắc của các phân nhánh tiến hóa.

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

Nghiên cứu thực hiện lấy 10 mẫu mô lá non tươi của cây Cựa gà tại xã Hợp Hòa, huyện Sơn Dương, tỉnh Tuyên Quang. Phương pháp chọn mẫu tập trung vào các chồi lá bánh tẻ không bị sâu bệnh, được làm sạch bề mặt bằng cồn 70% ngay tại thực địa trước khi đưa vào bảo quản lạnh sâu. Lý do lựa chọn mô lá non là nhằm hạn chế hàm lượng polyphenol và polysaccharide phức tạp thường tích lũy nhiều ở lá già, vốn là các tác nhân ức chế nghiêm trọng hoạt tính enzyme trong các phản ứng phân tử.

Quy trình phân tích giải phẫu học áp dụng phương pháp cắt lát vi phẫu tạm thời và nhuộm kép theo tiêu chuẩn thực vật học để quan sát chi tiết cấu trúc mô rễ, thân, lá dưới kính hiển vi quang học ở độ phóng đại 10x và 40x. Quá trình tách chiết DNA tổng số được tiến hành theo quy trình đệm CTAB cải tiến từ Dellaporta và Doyle, bao gồm các bước ủ nhiệt ở 65 độ C trong 60 phút, kết tủa bằng hỗn hợp Chloroform - Isoamyl alcohol (tỷ lệ 24:1) và làm sạch tủa qua cồn 70% cùng cồn 95%. Hàm lượng và độ tinh sạch của DNA được định lượng bằng máy đo quang phổ hấp thụ ở bước sóng 260 nm và 280 nm, với tiêu chuẩn tỷ số OD260/OD280 nằm trong ngưỡng tối ưu từ 1,8 đến 2,0.

Phản ứng chuỗi trùng hợp PCR khuếch đại 3 vùng gen matK, rbcL và ITS2 được thiết lập trong tổng thể tích 15 µL chứa 1,5 µL đệm DreamTaq 10X, 0,075 µL dNTP 40 mM, 0,5 µL mồi xuôi và ngược đặc hiệu, 0,1 µL enzyme Taq polymerase (1 UI/µL) và 1 µL DNA khuôn. Chu trình nhiệt bao gồm 35 chu kỳ với nhiệt độ bắt mồi dao động từ 52 độ C đến 55 độ C và giai đoạn kéo dài ở 72 độ C. Lý do lựa chọn hệ thống giải trình tự tự động ABI 3500xL Genetic Analyzer và phần mềm tin sinh học MEGA X là nhằm đảm bảo tốc độ đọc chuẩn xác từng nucleotide, loại bỏ sai số chủ quan, đồng thời cho phép đối sánh đồng thời với hàng trăm trình tự tương đồng trên ngân hàng gen quốc tế GenBank thông qua thuật toán BLAST.

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

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

Nghiên cứu đã ghi nhận 4 phát hiện quan trọng có giá trị khoa học cao về mặt hình thái học và di truyền phân tử:

Thứ nhất, về đặc điểm thực vật học và cấu trúc giải phẫu, loài Cựa gà tại Tuyên Quang là dạng cây bụi lớn có chiều dài thân phát triển vượt trên 10 m, vỏ cành chuyển từ xanh nâu sang xám trắng, thân mang gai nhọn mọc so le phân nhánh đặc trưng. Lá kép lông chim một lần có từ 3 đến 7 lá chét, cụm hoa mọc ở nách lá dài 4 đến 9 cm với hoa nhỏ màu trắng ngà trổ vào tháng 4 và tháng 5; quả dẹt dài 6 đến 7 cm mang hạt hình thận dài khoảng 12 mm thu hái vào tháng 10 đến tháng 12. Vi phẫu cắt ngang rễ cho thấy tỷ lệ mô gỗ sơ cấp chiếm khoảng 50% cấu trúc, trong khi cấu tạo thân thể hiện các bó mạch dẫn xếp thành vòng tròn liên tục, mang đầy đủ đặc trưng giải phẫu của nhóm thực vật hai lá mầm thuộc họ Fabaceae.

Thứ hai, kết quả tách chiết từ 10 mẫu lá non cho thấy mẫu số 10 có độ tinh sạch vượt trội với nồng độ DNA đạt 248,6 ng/µL và tỷ số hấp thụ OD260/OD280 đạt 1,87 (dao động toàn bộ 10 mẫu từ 1,83 đến 2,08). Vùng gen matK sau khi nhân bản bằng PCR và giải trình tự đạt chiều dài chính xác 700 bp. Đối sánh trên GenBank xác định độ tương đồng của đoạn gen matK với loài Dalbergia rimosa đạt từ 98,58% đến 99,40%, xuất hiện 22 vị trí sai khác nucleotide nằm từ vị trí 53 đến 752.

Thứ ba, vùng gen rbcL được khuếch đại thành công với kích thước 583 bp, đạt tỷ lệ tương đồng rất cao từ 99,12% đến 99,30% so với trình tự Dalbergia rimosa mang mã số KM510295.1 trên GenBank. Đoạn gen này ghi nhận 18 vị trí thay thế nucleotide và hiện tượng mất 24 nucleotide ở các vị trí đầu và cuối chuỗi so với mẫu đối chứng.

Thứ tư, phân tích vùng gen chỉ thị nhân ITS2 cho thấy sản phẩm PCR có kích thước 372 bp, đạt độ tương đồng 96,49% với loài Dalbergia rimosa. Dù có sự sai khác ở 37 vị trí nucleotide (trong đó có 35 vị trí thay thế base), việc đối sánh phát sinh chủng loại đã nhóm mẫu Tuyên Quang chung nhánh với Dalbergia rimosa với độ tin cậy bootstrap đạt mức ấn tượng từ 91% đến 95%. Điều này khẳng định mẫu cây Cựa gà nghiên cứu chính là loài Dalbergia rimosa.

Thảo luận kết quả

Sự tương hợp giữa các chỉ số di truyền phân tử và đặc điểm hình thái học phản ánh chính xác nguồn gốc tiến hóa của loài. Trong khi vùng gen lục lạp rbcL và matK có tốc độ tiến hóa chậm, giữ vai trò nhận diện tin cậy ở cấp độ chi, thì vùng đệm phiên mã nội bộ ITS2 trong nhân tế bào lại thể hiện tốc độ đột biến thay thế nucleotide đủ nhanh để phân biệt ranh giới giữa các loài gần gũi. Tỷ lệ tương đồng 96,49% cùng hệ số bootstrap 91% đến 95% ở chỉ thị ITS2 đã khắc phục hoàn toàn nhược điểm hệ số bootstrap thấp của gen matK, giúp khẳng định danh tính phân loại của Dalbergia rimosa một cách tuyệt đối.

Kết quả nghiên cứu này hoàn toàn đồng thuận với nhiều công bố học thuật quốc tế trên các nhóm thực vật thuộc họ Fabaceae và Rosaceae, nơi ITS2 luôn được chứng minh đạt tỷ lệ phân biệt loài thành công vượt trội từ 78% đến 100%. Toàn bộ dữ liệu phân tử của nghiên cứu có thể được trực quan hóa tối ưu thông qua các bảng ma trận đối sánh sai khác từng vị trí nucleotide và mô hình biểu đồ cây phát sinh chủng loại Maximum Likelihood. Trên biểu đồ cây phát sinh, các giá trị bootstrap hiển thị trực tiếp tại từng nút nhánh phân nhánh giúp người đọc dễ dàng nhận thấy sự tách biệt rõ ràng giữa 2 phân nhóm loài chính của chi Dalbergia và xác định vị trí phân loại chính xác của mẫu Tuyên Quang trong nhóm II.

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

Từ những phát hiện thực nghiệm, nghiên cứu đề xuất 4 nhóm giải pháp chiến lược nhằm thúc đẩy bảo tồn và ứng dụng nguồn dược liệu quý:

Thứ nhất, mở rộng ngân hàng mã vạch DNA đa gen bằng cách tiếp tục giải trình tự bổ sung hai vùng gen lạp thể trnL-trnF và psbA-trnH. Mục tiêu là nâng cao độ phân giải định danh loài đạt mức tuyệt đối 99,9%, xây dựng bộ cơ sở dữ liệu số hóa hoàn chỉnh cho toàn bộ chi Dalbergia tại Việt Nam trong khung thời gian 2024 - 2026, do Viện Nghiên cứu Hệ gen phối hợp cùng Đại học Thái Nguyên chủ trì triển khai.

Thứ hai, chuẩn hóa quy trình kỹ thuật PCR-barcoding vào hệ thống kiểm định chất lượng dược liệu quốc gia. Khuyến nghị Viện Dược liệu thuộc Bộ Y tế ứng dụng ngay bộ mồi chuẩn matK, rbcL và ITS2 để phát hiện dược liệu giả mạo hoặc pha trộn tạp chất trôi nổi trên thị trường, hướng tới mục tiêu giảm tỷ lệ nhầm lẫn nguồn gốc dược liệu xuống dưới 2% trong vòng 18 tháng áp dụng.

Thứ ba, thiết lập khu vực bảo vệ nguồn gen tại chỗ (in-situ) đối với quần thể Dalbergia rimosa tại huyện Sơn Dương. Sở Nông nghiệp và Phát triển nông thôn tỉnh Tuyên Quang cần phối hợp với chính quyền địa phương khoanh vùng sinh cảnh tự nhiên, đặt biển cảnh báo, nghiêm cấm khai thác gỗ và thân rễ trái phép, phấn đấu tăng mật độ cá thể tái sinh tự nhiên thêm 15% vào năm 2027.

Thứ tư, đẩy mạnh nghiên cứu sàng lọc hoạt tính dược lý chuyên sâu và nhân giống sinh dưỡng. Các trung tâm công nghệ sinh học dược liệu cần tập trung phân lập các hợp chất flavonoid đặc hữu, đồng thời hoàn thiện quy trình nhân giống vô tính in vitro đạt tỷ lệ sống sót của cây con trên 80% trong vòng 2 năm tới, đảm bảo nguồn cung nguyên liệu bền vững cho sản xuất thuốc và thực phẩm bảo vệ sức khỏe.

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

Nội dung và phương pháp tiếp cận của bản luận văn mang lại giá trị thực tiễn to lớn cho 4 nhóm đối tượng trọng tâm:

Thứ nhất, các nhà nghiên cứu, giảng viên, học viên cao học và sinh viên chuyên ngành Thực vật học, Sinh học thực nghiệm và Di truyền học. Tài liệu cung cấp quy trình mẫu chuẩn mực từ khâu lấy mẫu thực địa, giải phẫu hình thái, đến phương pháp chiết tách DNA bằng đệm CTAB cải tiến và phân tích phát sinh loài trên phần mềm MEGA X.

Thứ hai, chuyên viên phân tích tại các trung tâm kiểm nghiệm dược phẩm và cơ quan quản lý thị trường. Luận văn cung cấp bộ dữ liệu trình tự chuẩn (372 bp cho ITS2 và 583 bp cho rbcL) làm căn cứ đối soát phân tử nhanh chóng, giúp nhận diện chính xác các mẫu rễ, thân cây Cựa gà đã qua sơ chế hoặc nghiền bột mà phương pháp hình thái học bất khả thi.

Thứ ba, lực lượng kiểm lâm, cán bộ quản lý tài nguyên rừng và các tổ chức bảo tồn thiên nhiên. Các mô tả chi tiết về đặc điểm hình thái thân gai cựa gà, dạng quả dẹt và chu kỳ trổ hoa tháng 4 - 5 giúp nhân viên thực địa dễ dàng nhận biết, thống kê và lập phương án khoanh nuôi bảo vệ loài Dalbergia rimosa ngoài tự nhiên.

Thứ tư, các doanh nghiệp chế biến dược liệu và sản xuất thực phẩm chức năng. Việc nắm vững quy trình giám định phân tử giúp doanh nghiệp kiểm soát chặt chẽ nguyên liệu đầu vào đạt tiêu chuẩn Thực hành tốt trồng trọt và thu hái dược liệu (GACP-WHO), loại bỏ triệt để nguy cơ pha tạp nguyên liệu kém chất lượng.

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

Cây Cựa gà thu thập tại Tuyên Quang được định danh chính xác là loài gì?

Dựa trên kết quả phân tích giải trình tự vùng gen nhân ITS2 có kích thước 372 bp với độ tương đồng 96,49% và gen lục lạp rbcL đạt 99,30%, cây phát sinh loài suy luận theo mô hình Maximum Likelihood cho hệ số bootstrap đạt từ 91% đến 95%. Kết quả này khẳng định mẫu cây Cựa gà tại xã Hợp Hòa, Sơn Dương chính là loài Dalbergia rimosa.

Vì sao cần kết hợp cả gen lục lạp và gen nhân trong quy trình mã vạch DNA thực vật?

Các vùng gen lục lạp như matK dài 700 bp và rbcL dài 583 bp có tính bảo thủ cao, dễ khuếch đại, phù hợp để phân loại ở cấp độ chi họ. Trong khi đó, vùng đệm nhân ITS2 sở hữu 37 vị trí nucleotide sai khác, cung cấp độ biến thiên phân tử cần thiết để phân biệt chính xác từng loài thực vật có hình thái tương đồng.

Tiêu chí nào phản ánh chất lượng DNA tách chiết từ lá cây Cựa gà đạt yêu cầu thí nghiệm?

Chất lượng DNA được đánh giá thông qua chỉ số quang phổ hấp thụ OD260/OD280 đạt ngưỡng chuẩn 1,8 đến 2,0 và vạch điện di sắc nét trên gel agarose 1%. Trong nghiên cứu, mẫu số 10 đạt nồng độ 248,6 ng/µL với tỷ số OD260/OD280 là 1,87, đảm bảo loại bỏ sạch protein và tạp chất hữu cơ gây ức chế phản ứng PCR.

Những dấu hiệu hình thái nào giúp nhận biết loài Dalbergia rimosa tại thực địa một cách nhanh nhất?

Loài có thân dạng bụi lớn vươn dài trên 10 m, nổi bật với các gai nhọn phân nhánh mọc so le giống cựa gà. Cây có lá kép lông chim từ 3 đến 7 lá chét bóng nhẵn, chùm hoa nhỏ màu trắng dài 4 đến 9 cm trổ vào tháng 4 đến tháng 5 và quả dẹt dài 6 đến 7 cm chứa 1 hạt hình thận.

Nghiên cứu đóng góp giải pháp gì cho công tác bảo tồn chi Dalbergia đang bị đe dọa tuyệt chủng?

Nhiều loài trong chi Dalbergia đang ở mức Nguy cấp (EN) do nạn chặt phá rừng. Nghiên cứu cung cấp công cụ giám định phân tử nhanh, chính xác ngay cả với mẫu gỗ hay dược liệu khô, thiết lập nền tảng khoa học giúp cơ quan thực thi pháp luật kiểm soát buôn bán trái phép và bảo tồn nguồn gen quý hiếm.

Kết luận

Luận văn đã hoàn thành xuất sắc các mục tiêu nghiên cứu thông qua 5 kết quả cốt lõi:

  • Mô tả đầy đủ và chi tiết đặc điểm hình thái học cùng cấu trúc vi phẫu rễ, thân, lá của cây Cựa gà, xác định gỗ sơ cấp chiếm khoảng 50% cấu tạo rễ đặc trưng cho thực vật hai lá mầm họ Đậu.
  • Phân lập, nhân bản và giải trình tự thành công 3 chỉ thị gen quan trọng gồm matK (700 bp), rbcL (583 bp) và ITS2 (372 bp) với độ tương đồng dao động từ 96,49% đến 99,40%.
  • Xác lập vùng gen ITS2 và rbcL là hai chỉ thị mã vạch DNA tối ưu với hệ số bootstrap hỗ trợ phân nhánh đạt từ 83% đến 95%.
  • Chuẩn hóa quy trình tách chiết DNA tổng số bằng đệm CTAB cải tiến, ghi nhận mẫu số 10 đạt độ tinh sạch lý tưởng với tỷ số OD260/OD280 là 1,87.
  • Định danh chính xác và khẳng định mẫu cây Cựa gà thu thập tại Tuyên Quang thuộc về loài Dalbergia rimosa.

Đóng góp chính của luận văn là đã bổ sung dữ liệu phân tử có giá trị vào ngân hàng mã vạch DNA thực vật Việt Nam. Về kế hoạch tiếp theo trong vòng 12 đến 24 tháng tới, nhóm nghiên cứu cần mở rộng giải mã các vùng liên gen trnL-trnF và psbA-trnH. Các cơ quan quản lý dược liệu và kiểm lâm hãy ứng dụng ngay bộ mã vạch phân tử này vào quy trình kiểm tra thực tế để nâng cao hiệu lực bảo tồn nguồn gen tự nhiên.