Tổng quan nghiên cứu

Chi Hoya thuộc họ Thiên lý (Asclepiadaceae) sở hữu khoảng 350 đến 450 loài trên phạm vi toàn cầu, trong đó hệ thực vật Việt Nam đã ghi nhận sự xuất hiện của ít nhất 24 loài tự nhiên. Cây Lan tai cáo (Hoya parasitica (Roxb.) Wight), hay còn gọi là Dây lưỡi lợn hoặc Hồ hoa ký sinh, là loài thảo dược sống phụ sinh lâu năm tại các vùng rừng nhiệt đới ẩm. Loài cây này có đặc điểm hình thái đặc trưng với thân dây leo màu đỏ đậm, lá dày mọng nước và cụm hoa dạng tán chứa khoảng 22 hoa nhỏ màu trắng hồng. Dịch chiết từ lá tươi Lan tai cáo chứa nhiều hoạt chất sinh học quý có tác dụng kháng khuẩn, chống oxy hóa, hỗ trợ làm lành vết thương, giảm sẹo và phòng ngừa suy thận cấp.

Tuy nhiên, sự đa dạng hình thái và mức độ tương đồng kiểu hình rất cao giữa các loài trong chi Hoya gây khó khăn lớn cho công tác phân loại truyền thống. Khi mẫu thực vật ở dạng bột khô, bị biến dạng hoặc chưa phát triển đầy đủ cơ quan sinh sản, việc nhận dạng thuần túy bằng mắt thường rất dễ dẫn đến sai sót. Nghiên cứu này được thực hiện nhằm phân lập, giải mã đặc điểm trình tự ba vùng gene gồm matK, rbcL (từ hệ gene lục lạp) và ITS (từ hệ gene nhân), từ đó phân tích sự phát sinh loài của Lan tai cáo mẫu thu tại tỉnh Thái Nguyên (ký hiệu TN1). Nghiên cứu được hoàn thiện vào tháng 8 năm 2023 trong khuôn khổ đề tài cấp Bộ Giáo dục và Đào tạo mang mã số B2022-TNA-22-CT562. Kết quả mang lại ý nghĩa khoa học và thực tiễn sâu sắc, thiết lập bộ dữ liệu chuẩn với độ phân giải di truyền đạt từ 95,98% đến 99,82%, tạo nền tảng vững chắc cho việc định danh nguồn gốc dược liệu và phục vụ bảo tồn nguồn gene thực vật quý tại Việt Nam.

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 vận dụng lý thuyết Tiến hóa phân tử (Molecular Evolution) kết hợp với Phát sinh chủng loại học phân tử (Molecular Phylogenetics) và mô hình thay thế nucleotide Tamura-Nei. Cơ sở của phương pháp dựa trên nguyên lý: mức độ sai khác trong trình tự nucleotide của các phân tử DNA phản ánh chính xác khoảng cách di truyền và lịch sử tiến hóa giữa các taxon sinh vật.

Khung nghiên cứu tích hợp ba khái niệm chuyên ngành then chốt:

  • Hệ gene lục lạp (chloroplast DNA - cpDNA): Phân tử DNA vòng kép có kích thước khoảng 85 đến 2000 kb, bảo tồn cấu trúc cao, di truyền theo dòng mẹ và mã hóa 30 vùng phức hợp quang hợp. Hai vùng gene lục lạp chủ đạo được phân tích là matK (mã hóa enzyme maturase K) và rbcL (mã hóa tiểu đơn vị lớn của enzyme RuBisCO).
  • Vùng đệm phiên mã nội ribosom (Internal Transcribed Spacer - ITS): Thuộc hệ gene nhân, sở hữu tốc độ tiến hóa nhanh, tạo sự đa hình phong phú để nhận diện loài.
  • Mã vạch DNA (DNA Barcoding): Khung phân loại chuẩn do Hiệp hội Mã vạch Sự sống (CBOL) đề xuất nhằm giám định sinh vật bằng các phân đoạn DNA ngắn tiêu chuẩn.

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

Về nguồn dữ liệu và chọn mẫu: Nghiên cứu sử dụng 100 mg mô lá non tươi thu thập từ mẫu Lan tai cáo TN1 tại Thái Nguyên, lưu giữ tại Vườn thực nghiệm và Phòng thí nghiệm Công nghệ tế bào thực vật, Trường Đại học Sư phạm - Đại học Thái Nguyên. Phương pháp chọn mẫu có chủ đích được áp dụng để đảm bảo mẫu nghiên cứu đạt độ thuần khiết cao, không sâu bệnh và đối sánh toàn diện với 82 trình tự tham chiếu quốc tế trích xuất từ ngân hàng gene GenBank (NCBI).

Về kỹ thuật tách chiết và nhân bản: Tách chiết DNA tổng số theo phương pháp CTAB cải tiến của Doyle. Đo quang phổ hấp thụ ghi nhận nồng độ DNA đạt 356,20 ng/µL với tỷ số tinh sạch A260/A280 đạt 2,07. Phản ứng PCR thể tích 50 µL thực hiện qua 35 chu kỳ với các cặp mồi chuyên biệt: mồi matK (nhiệt độ bắt mồi 58°C, khuếch đại đoạn 0,9 kb), mồi rbcL (nhiệt độ 56°C, khuếch đại đoạn 0,6 kb) và mồi ITS (nhiệt độ 65°C, khuếch đại đoạn 0,9 kb).

Về giải trình tự và phân tích tiến hóa: Sản phẩm PCR được tinh sạch bằng bộ kít GenJET và giải trình tự nucleotide tự động trên máy ABI 3500xL. Trình tự sau đó được xử lý bằng phần mềm BioEdit và đối soát BLAST. Lý do lựa chọn phương pháp Cực đại xác suất (Maximum Likelihood - ML) kết hợp kiểm định bootstrap 1000 lần lặp trên phần mềm MEGA X là vì phương pháp này ước lượng chính xác nhất các tham số thay thế nucleotide, loại bỏ sai lệch do phân kỳ dị thời và đem lại cấu trúc cây phát sinh loài có độ tin cậy thống kê cao nhất. Toàn bộ thực nghiệm hoàn thành trong lộ trình 2022 - 2023.

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

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

Quá trình phân tích thực nghiệm và giải trình tự nucleotide đã mang lại các phát hiện cốt lõi sau:

  • Phát hiện về vùng gene matK: Phân lập thành công đoạn gene matK có chiều dài 729 nucleotide từ mẫu Lan tai cáo TN1. Khi đối sánh với hai trình tự Hoya parasitica trên GenBank (mã số MH598756 và MH598772), tỷ lệ tương đồng đạt 98,61%. Giữa ba mẫu tồn tại 16 vị trí sai khác nucleotide; trong đó mẫu TN1 khác biệt 11 vị trí so với mẫu WT104 và 13 vị trí so với mẫu ZP28. Cây phát sinh loài xây dựng từ 29 trình tự phân chia thành 3 nhóm rõ rệt, mẫu TN1 kết cụm chặt chẽ cùng nhóm với mẫu MH598756 với giá trị bootstrap đạt 98%.
  • Phát hiện về vùng gene rbcL: Thu nhận đoạn gene rbcL có kích thước 547 nucleotide. Mức độ tương đồng giữa Lan tai cáo TN1 và các loài thuộc chi Hoya đạt rất cao, dao động từ 98,90% đến 99,82%. Trình tự này chỉ sai khác 2 nucleotide tại vị trí 529 và 547 so với Hoya carnosa (NC_045868.1), và chỉ khác 1 nucleotide tại vị trí 547 so với Hoya sp. gp-51. Cây phát sinh loài từ 37 trình tự chia làm 2 nhánh chính, mẫu TN1 nằm trong nhóm lớn gồm 32 loài.
  • Phát hiện về vùng gene ITS: Phân lập thành công vùng ITS kích thước 820 nucleotide, cho tỷ lệ tương đồng 95,98% so với 16 trình tự Hoya trên NCBI. Trình tự mẫu TN1 sai khác 14 vị trí so với Hoya monetteae (HE794431.1) và 13 vị trí so với Hoya cf. salweenica (HE794456.1). Trên cây phát sinh chủng loại 17 trình tự, mẫu TN1 thuộc nhánh IV (nhóm 2b) cùng hai loài trên với chỉ số bootstrap đạt 92%.

Thảo luận kết quả

Sự khác biệt về mức độ đa hình giữa các vùng gene phản ánh rõ nét cơ chế tiến hóa phân tử ở thực vật. Vùng gene rbcL chịu áp lực chọn lọc tự nhiên thanh lọc rất mạnh nhằm duy trì chức năng xúc tác cố định carbon của enzyme RuBisCO trong quang hợp. Do đó, rbcL có tốc độ biến đổi nucleotide cực kỳ chậm, dẫn đến độ tương đồng vượt trên 98,90% giữa các loài Hoya. Đặc tính này khiến rbcL chỉ phù hợp làm chỉ thị định danh ở cấp độ chi hoặc họ, không đủ độ nhạy để phân biệt giữa các loài gần gũi.

Ngược lại, vùng gene lục lạp matK (16 vị trí đột biến trên 729 bp) và vùng rDNA nhân ITS (22 vị trí biến dị trên 820 bp) sở hữu tốc độ thay thế nucleotide nhanh hơn hẳn. Các đột biến điểm này tạo nên tính đa hình đặc hiệu loài, cho phép tách biệt rõ ràng các quần thể Lan tai cáo phân bố tại các vùng địa lý khác nhau.

Dữ liệu nghiên cứu có thể được trực quan hóa hiệu quả thông qua hai dạng thức trình bày:

  • Biểu đồ cây phát sinh chủng loại dạng nhánh ngang (cladogram) thể hiện cấu trúc 3 nhóm phân loại của matK và 4 nhánh tiến hóa của ITS kèm giá trị bootstrap từ 92% đến 98% tại các nút nhánh.
  • Bảng ma trận khoảng cách di truyền đối sánh từng cặp loài nhằm biểu thị trực quan tỷ lệ đột biến thay thế nucleotide và phần trăm tương đồng từ 95,98% đến 99,82%.

Kết quả chứng minh sự phối hợp đa locus giữa matK và ITS là giải pháp tối ưu giúp nâng tỷ lệ định danh chính xác loài Lan tai cáo lên mức xấp xỉ 100%, khắc phục hoàn toàn nhược điểm của phương pháp hình thái học.

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

Nhằm phát huy giá trị của kết quả nghiên cứu trong thực tiễn kiểm nghiệm dược liệu và bảo tồn sinh học, bốn giải pháp trọng tâm được đề xuất như sau:

  • Ứng dụng tổ hợp mã vạch DNA đa locus matK, rbcL và ITS vào kiểm định chất lượng dược liệu. Chủ thể thực hiện là Viện Dược liệu phối hợp với các Trung tâm Kiểm nghiệm Thuốc trên toàn quốc. Mục tiêu hướng tới đạt độ chính xác 100% trong việc xác định nguồn gốc nguyên liệu Lan tai cáo thương phẩm, loại bỏ triệt để nguy cơ pha tạp các loài Hoya kém chất lượng. Thời gian triển khai hoàn thiện quy trình trong 6 tháng đầu năm 2024.
  • Xây dựng cơ sở dữ liệu số hóa mã vạch DNA cho 24 loài thuộc chi Hoya tại Việt Nam. Chủ thể thực hiện gồm các trường đại học sư phạm, viện nghiên cứu công nghệ sinh học và di truyền học. Mục tiêu đăng ký bổ sung từ 50 đến 100 trình tự chuẩn quốc gia lên ngân hàng dữ liệu quốc tế GenBank trong giai đoạn 2024 - 2025.
  • Thiết kế bộ sinh phẩm chẩn đoán nhanh dạng kít PCR đặc hiệu loài dựa trên 16 vị trí nucleotide sai khác của gene matK. Chủ thể thực hiện là các doanh nghiệp công nghệ sinh học và phòng thí nghiệm trọng điểm. Mục tiêu giảm thời gian giám định từ 48 giờ xuống dưới 4 giờ, đồng thời cắt giảm khoảng 60% chi phí phân tích so với quy trình giải trình tự toàn phần. Thời gian thử nghiệm trong vòng 12 tháng.
  • Mở rộng phạm vi thu thập mẫu Lan tai cáo tại các vùng sinh thái trọng điểm như Võ Xá (Hà Nam), Núi Dinh (Bà Rịa - Vũng Tàu) và khu vực Tây Nguyên. Chủ thể thực hiện là nhóm nghiên cứu chuyên ngành phối hợp cùng các ban quản lý rừng đặc dụng. Mục tiêu thu thập tối thiểu 30 mẫu địa phương trong chu kỳ 2024 - 2026 để đánh giá toàn diện mức độ phân hóa di truyền nội loài.

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

Luận văn là tài liệu chuyên khảo giá trị cho bốn nhóm đối tượng chính:

  • Các nhà nghiên cứu và giảng viên chuyên ngành Di truyền học, Sinh học phân tử, Thực vật học: Có thể khai thác quy trình kỹ thuật tách chiết DNA đạt nồng độ 356,20 ng/µL, chu trình nhiệt PCR 35 chu kỳ và bộ thông số mồi chuẩn để phục vụ công tác giảng dạy, nghiên cứu phân loại học phân tử họ Thiên lý.
  • Chuyên viên kiểm nghiệm dược phẩm và kiểm soát chất lượng tại các doanh nghiệp đông dược: Ứng dụng dữ liệu 729 bp của gene matK và 820 bp của vùng ITS làm chuẩn đối soát, phát hiện nguyên liệu giả mạo và bảo đảm an toàn dược phẩm cho các sản phẩm trị bỏng, mau liền sẹo và phục hồi chức năng thận.
  • Học viên cao học và nghiên cứu sinh ngành Công nghệ sinh học: Tham khảo mô hình phân tích tiến hóa Maximum Likelihood với 1000 lần lặp bootstrap trên MEGA X làm khung phương pháp chuẩn cho các đề tài phân tích di truyền tiến hóa phân tử ở thực vật bậc cao.
  • Cán bộ quản lý tài nguyên thiên nhiên và các tổ chức bảo tồn đa dạng sinh học: Sử dụng cơ sở dữ liệu di truyền để giám định chính xác các loài trong quần thể 350 đến 450 loài thuộc chi Hoya, phục vụ xây dựng kế hoạch bảo tồn nguồn gene bản địa quý hiếm tại Việt Nam.

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

  • Tại sao cần kết hợp cả ba chỉ thị matK, rbcL và ITS thay vì chỉ sử dụng một chỉ thị đơn lẻ?
    Gene rbcL (547 bp) có tính bảo thủ cao với độ tương đồng 98,90% đến 99,82%, rất hữu hiệu để định danh ở cấp độ chi nhưng khó phân biệt giữa các loài gần nhau. Trong khi đó, matK (729 bp) và ITS (820 bp) sở hữu độ biến dị lớn (sai khác từ 16 đến 22 vị trí nucleotide). Việc kết hợp cả chỉ thị lục lạp và chỉ thị nhân giúp kết quả giám định đạt độ chính xác xấp xỉ 100%.

  • Độ tinh sạch của mẫu DNA tách chiết từ lá Lan tai cáo đạt tiêu chuẩn như thế nào?
    Áp dụng quy trình CTAB cải tiến trên 100 mg lá non, DNA tổng số thu được đạt hàm lượng 356,20 ng/µL cùng tỷ số quang phổ A260/A280 bằng 2,07. Sau bước tinh sạch phân đoạn PCR, nồng độ DNA đạt từ 33,6 ng/µL đến 61,9 ng/µL, hoàn toàn đáp ứng các tiêu chuẩn khắt khe để giải trình tự trên hệ thống ABI 3500xL.

  • Giá trị bootstrap 98% trên cây phát sinh chủng loại của gene matK phản ánh điều gì?
    Chỉ số bootstrap 98% qua 1000 lần tái lấy mẫu ngẫu nhiên chứng minh độ tin cậy thống kê cực cao của nhánh tiến hóa. Kết quả này khẳng định mẫu Lan tai cáo Thái Nguyên (TN1) có quan hệ di truyền rất gần gũi với mẫu Hoya parasitica mã số MH598756 trên GenBank, loại trừ khả năng phân nhóm do ngẫu nhiên.

  • Lan tai cáo phân bố ở đâu tại Việt Nam và có công dụng chữa bệnh gì nổi bật?
    Tại Việt Nam, Lan tai cáo mọc phụ sinh tự nhiên tại Hà Nam (Võ Xá), Bà Rịa - Vũng Tàu (Núi Dinh) và Thái Nguyên. Cây mang cụm hoa khoảng 22 hoa màu trắng hồng. Dịch chiết lá tươi có hoạt tính kháng khuẩn, chống oxy hóa, hỗ trợ làm lành vết thương ngoài da nhanh chóng, lợi sữa và hỗ trợ điều trị suy thận cấp.

  • Đóng góp khoa học nổi bật nhất của luận văn đối với ngành sinh học phân tử là gì?
    Luận văn đã công bố bộ trình tự hoàn chỉnh gồm 729 bp gene matK, 547 bp gene rbcL và 820 bp vùng ITS của loài Hoya parasitica tại Thái Nguyên. Đây là công trình đầu tiên tại Việt Nam giải mã đồng thời cả ba locus mã vạch DNA tiêu chuẩn trên đối tượng này, đóng góp dữ liệu quan trọng cho kho tàng di truyền thực vật quốc gia.

Kết luận

  • Giải mã và xác định thành công kích thước nucleotide chuẩn của ba vùng chỉ thị phân tử: matK (729 bp), rbcL (547 bp) và ITS (820 bp) từ mẫu Lan tai cáo Thái Nguyên (TN1).
  • Chứng minh gene rbcL có mức độ tương đồng rất cao (98,90% đến 99,82%), đóng vai trò là chỉ thị tin cậy để định danh ở cấp độ chi trong họ Thiên lý.
  • Khẳng định vùng gene matK (tương đồng 98,61%, sai khác 16 vị trí) và vùng ITS (tương đồng 95,98%, sai khác 22 vị trí) là hai mã vạch DNA tiềm năng nhất để phân loại chính xác ở cấp độ loài với độ tin cậy bootstrap đạt từ 92% đến 98%.
  • Hoàn thiện quy trình tách chiết DNA đạt nồng độ 356,20 ng/µL (A260/A280 = 2,07) cùng hệ thống thông số nhiệt chu kỳ PCR chuẩn hóa từ 56°C đến 65°C cho từng locus.
  • Kế hoạch triển khai giai đoạn 2024 - 2026: Đẩy mạnh phân tích tổ hợp đa locus (matK + rbcL + ITS) và mở rộng phạm vi thu thập mẫu trên toàn quốc nhằm hoàn thiện kít chẩn đoán thương mại phục vụ kiểm định dược liệu.

Quý cơ quan, viện nghiên cứu và doanh nghiệp có nhu cầu khai thác cơ sở dữ liệu trình tự nucleotide hoặc hợp tác phát triển quy trình giám định dược liệu Lan tai cáo xin vui lòng liên hệ Bộ môn Di truyền học và Công nghệ sinh học, Trường Đại học Sư phạm - Đại học Thái Nguyên để nhận tài liệu hỗ trợ chuyên sâu.