Tổng quan nghiên cứu

Việt Nam sở hữu hệ thực vật phong phú với gần 12000 loài thực vật bậc cao có mạch thuộc 2256 chi và 305 họ, mở ra tiềm năng to lớn cho ngành công nghệ sinh học và phát triển dược liệu. Trong số các loài cây thuốc quý, cây dừa cạn (Catharanthus roseus (L.) G. Don thuộc họ Trúc đào Apocynaceae) đóng vai trò đặc biệt quan trọng khi chứa hơn 60 loại alkaloid terpenoid indole khác nhau. Các hoạt chất tiêu biểu như vinblastine và vincristine chiếm tỷ lệ khoảng 1 phần vạn trong lá khô nhưng là thành phần thiết yếu trong phác đồ điều trị ung thư bạch cầu và ung thư hạch. Tuy nhiên, việc nhận diện và phân loại chính xác các giống dừa cạn thương mại hiện nay gặp nhiều rào cản khi chỉ dựa vào hình thái truyền thống, nhất là ở giai đoạn cây non hoặc khi dược liệu đã qua sơ chế, băm sấy.

Đề tài nghiên cứu được thực hiện với mục tiêu phân lập, giải trình tự đoạn gen lục lạp rpoC1 và so sánh đặc điểm hình thái giữa ba giống dừa cạn có màu sắc hoa khác nhau (hoa hồng tím, hoa trắng vàng, hoa trắng đỏ). Nghiên cứu được triển khai trong phạm vi thời gian từ tháng 6 năm 2015 đến tháng 4 năm 2016 tại hai địa bàn thu mẫu là thành phố Hà Giang và thành phố Thái Nguyên, thực nghiệm tại Phòng thí nghiệm Công nghệ Gen, Trường Đại học Sư phạm thuộc Đại học Thái Nguyên. Kết quả nghiên cứu cung cấp cơ sở dữ liệu phân tử với độ tương đồng chuỗi đạt 96,6% đến 100%, thiết lập dữ liệu mã vạch DNA chuẩn xác nhằm phục vụ công tác giám định nguồn gen và nâng cao hiệu quả kiểm soát chất lượng dược liệu trong y tế cộng đồng.

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 ứng dụng lý thuyết phân loại học phân tử kết hợp với khái niệm mã vạch DNA (DNA barcoding) do Paul Hebert đề xuất từ năm 2003. Phương pháp này dựa trên nguyên lý sử dụng một đoạn trình tự nucleotide chuẩn, ngắn và bảo thủ trong hệ gen để định danh chính xác loài sinh vật tương tự như cơ chế quét mã vạch hàng hóa. Trong hệ gen thực vật bậc cao, bộ gen lục lạp (cpDNA) có cấu trúc dạng vòng kép với kích thước từ 120000 đến 170000 cặp base (bp), chứa khoảng 100 gen chức năng và có tính di truyền theo dòng mẹ.

Khung lý thuyết của đề tài tập trung vào gen rpoC1, một gen lục lạp có chiều dài dao động từ 409 đến 545 bp, đảm nhiệm chức năng mã hóa cho tiểu đơn vị beta prime của enzyme RNA polymerase. Nhờ tốc độ tiến hóa ổn định và mức độ bảo thủ cao, gen rpoC1 trở thành chỉ thị phân tử lý tưởng để phát hiện các đột biến điểm (SNPs) và đánh giá sai khác di truyền giữa các giống dưới loài. Bên cạnh đó, lý thuyết hóa phân loại học về nhóm alkaloid terpenoid indole cũng được tích hợp để giải thích mối tương quan giữa đặc điểm phân loại học và giá trị dược lý của loài Catharanthus roseus.

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

Nguồn dữ liệu của nghiên cứu bao gồm 4 mẫu thực vật dừa cạn tươi được thu thập trực tiếp và 2 trình tự đối chứng quốc tế mang mã số KC561139 và JN115007 tải từ Ngân hàng Gen Quốc tế NCBI. Phương pháp chọn mẫu là chọn mẫu có chủ đích, thu thập các cá thể trưởng thành khỏe mạnh, đại diện cho 3 kiểu hình màu hoa tại tỉnh Hà Giang (TIM_HG, TV_HG) và tỉnh Thái Nguyên (TIM_TN, TRANG_TN). Cỡ mẫu thực nghiệm gồm 4 mẫu đại diện được thuần hóa và chăm sóc trong cùng điều kiện ánh sáng và chế độ nước để loại trừ biến dị thích nghi môi trường.

Quy trình phân tích phân tử được thực hiện theo các bước nghiêm ngặt:

  • Tách chiết DNA tổng số từ 200 mg mô lá non bằng phương pháp CTAB 2% với nitơ lỏng ở âm 196 độ C. Phương pháp CTAB được lựa chọn vì khả năng loại bỏ triệt để polysaccharide và polyphenol đặc thù trong lá dừa cạn.
  • Khuếch đại đoạn gen rpoC1 qua phản ứng PCR với cặp mồi đặc hiệu rpoC1-1f và rpoC1-3r, chu trình nhiệt 35 chu kỳ với nhiệt độ gắn mồi 55 độ C trong 40 giây. Kỹ thuật PCR Master Mix 2X tối ưu hóa độ chính xác của phản ứng nhân bản.
  • Tinh sạch sản phẩm PCR qua cột lọc GenJET PCR Purification và giải trình tự nucleotide tự động trên hệ thống máy ABI PRISM 3100 Avant Genetic Analyzer bằng hóa chất BigDye Terminator v3.1.
  • Xử lý số liệu và phân tích phát sinh chủng loại bằng các phần mềm tin sinh học chuyên dụng BioEdit, BLAST và DNAstar.

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

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

Thứ nhất, về đặc điểm thực vật học, cả ba giống dừa cạn khảo sát đều có thân thảo cao tới 100 cm, phân cành nhiều, rễ cái cắm sâu từ 10 đến 20 cm (có thể đạt 35 đến 40 cm), lá đơn mọc đối dài 3 đến 8 cm với 12 đến 14 cặp gân phụ. Sự khác biệt hình thái rõ nét nhất nằm ở màu sắc cánh hoa và tâm hoa: giống hoa hồng tím có thân xanh phớt hồng, giống hoa trắng đỏ có cánh trắng tâm đỏ, và giống hoa trắng vàng có thân xanh nhạt cùng cánh trắng tâm vàng.

Thứ hai, nghiên cứu đã phân lập và nhân bản thành công đoạn gen rpoC1 với kích thước đồng nhất là 505 bp ở cả 4 mẫu nghiên cứu, mã hóa chuỗi peptide gồm 168 amino acid. Thành phần nucleotide ở các mẫu TIM_HG, TIM_TN và TV_HG gồm 139 base loại A, 161 base loại T, 84 base loại C và 121 base loại G; riêng mẫu TRANG_TN có 138 base loại A và 162 base loại T.

Thứ ba, phân tích BLAST khẳng định độ tương đồng chuỗi giữa các mẫu nghiên cứu với dữ liệu chuẩn trên NCBI đạt mức 99%. So sánh nội bộ cho thấy mẫu TIM_HG và TV_HG tương đồng 100% với mã số KC561139, trong khi mẫu TIM_TN đạt 99,6% và TRANG_TN đạt 99,4%.

Thứ tư, nghiên cứu phát hiện các vị trí đột biến thay thế nucleotide đơn lẻ tại các vị trí base 494, 495 và 496. Sự thay đổi này dẫn đến đột biến sai nghĩa tại vị trí amino acid thứ 165: mẫu TIM_HG và TV_HG mang Methionine, mẫu TIM_TN chuyển thành Serine, còn mẫu TRANG_TN chuyển thành Cysteine.

Thảo luận kết quả

Sự sai khác nucleotide tại các vị trí 494, 495 và 496 là minh chứng rõ ràng cho tính đa dạng di truyền ở cấp độ dưới loài của cây dừa cạn tại khu vực miền núi phía Bắc. Dữ liệu thực nghiệm có thể được trực quan hóa tối ưu thông qua bảng đối chiếu ma trận tương đồng nucleotide và sơ đồ cây phát sinh chủng loại Dendrogram xây dựng theo thuật toán Neighbor-Joining.

Sơ đồ hình cây phân chia 6 mẫu nghiên cứu thành hai nhánh rõ rệt với khoảng cách di truyền nucleotide đạt 5,4%. Nhánh thứ nhất bao gồm các mẫu thu thập tại Việt Nam, trong đó nhóm phụ gồm TIM_HG, TV_HG, TIM_TN liên kết chặt chẽ với nhau, còn TRANG_TN tách thành một phân nhóm độc lập. Đáng chú ý, khi so sánh ở cấp độ chuỗi amino acid suy diễn, khoảng cách di truyền giữa các mẫu giảm xuống chỉ còn 0,4%, phản ánh tính bảo thủ cao trong cấu trúc chức năng của protein RNA polymerase lục lạp. Kết quả này hoàn toàn phù hợp với các nghiên cứu quốc tế về tính ổn định của gen rpoC1, khẳng định giá trị của chỉ thị này trong việc xây dựng hệ thống barcode phân biệt chuẩn xác các giống dừa cạn bản địa.

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

  • Xây dựng ngân hàng dữ liệu mã vạch DNA đa locus: Viện Dược liệu phối hợp cùng các trường đại học thực hiện tích hợp chỉ thị rpoC1 với các vùng gen matK, ITS và rbcL, hướng tới mục tiêu đạt tỷ lệ giám định chính xác 100% cho mọi dạng nguyên liệu dừa cạn thô và bột thuốc trong giai đoạn 2026-2028.
  • Chuẩn hóa quy trình giám định phân tử nhanh: Cơ quan quản lý Dược và các trung tâm kiểm nghiệm cần áp dụng quy trình chiết xuất CTAB cải tiến kết hợp PCR định lượng, rút ngắn thời gian phân tích mẫu xuống dưới 24 giờ kể từ năm 2027.
  • Quy hoạch vùng trồng dược liệu chuẩn GACP-WHO: Sở Nông nghiệp và Phát triển Nông thôn tại các tỉnh Hà Giang, Thái Nguyên và vùng phụ cận cần hướng dẫn người dân chọn lọc giống thuần chủng, đảm bảo kiểm soát độ tinh sạch nguồn gen nhằm duy trì hàm lượng alkaloid toàn phần trong rễ đạt từ 0,7% đến 2,4%.
  • Ứng dụng chỉ thị phân tử trong chọn tạo giống: Các viện công nghệ sinh học cần ứng dụng dấu chuẩn gen rpoC1 để chọn lọc các dòng dừa cạn hoa trắng có hoạt tính vincristine và vinblastine cao vượt trội, nâng cao sản lượng chiết xuất dược chất phục vụ sản xuất thuốc điều trị ung thư trong nước giai đoạn 2027-2030.

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

  • Nhà nghiên cứu và giảng viên chuyên ngành Công nghệ sinh học, Dược học: Sử dụng số liệu giải trình tự gen rpoC1, bảng ma trận phân ly di truyền và quy trình tách chiết DNA làm tài liệu tham khảo trong giảng dạy và phát triển các đề tài phân loại học phân tử thực vật.
  • Chuyên viên kiểm nghiệm tại các trung tâm kiểm nghiệm dược phẩm: Khai thác quy trình phân tích PCR và dữ liệu đột biến nucleotide tại vị trí 494-496 để xây dựng quy chuẩn phát hiện dược liệu giả mạo và kiểm soát chất lượng nguyên liệu dừa cạn đầu vào.
  • Doanh nghiệp sản xuất dược phẩm và đơn vị phát triển vùng trồng: Ứng dụng các đặc điểm thực vật học và dữ liệu phân loại giống hoa trắng, hoa hồng tím để chọn lọc nguồn giống thuần chủng, tối ưu hóa quy trình thu hái theo từng bộ phận thân, lá, rễ.
  • Học viên cao học và sinh viên ngành Sinh học, Nông nghiệp công nghệ cao: Tham khảo phương pháp luận khoa học, cách thức xử lý dữ liệu tin sinh học bằng BioEdit, BLAST và phương pháp trình bày cây phát sinh chủng loại cho các công trình nghiên cứu chuyên sâu.

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

Tại sao gen rpoC1 lại được lựa chọn làm mã vạch DNA cho cây dừa cạn?

Gen rpoC1 nằm trong hệ gen lục lạp có cấu trúc vòng kép bền vững, kích thước trung bình 505 bp rất thuận lợi cho việc khuếch đại PCR và giải trình tự. Gen này sở hữu các vùng bảo thủ xen kẽ những vị trí đột biến điểm đặc thù, giúp phân biệt rõ ràng các giống cây mà không bị ảnh hưởng bởi môi trường bên ngoài.

Ba giống dừa cạn nghiên cứu có sự khác biệt gì ở cấp độ phân tử?

Phân tích trình tự cho thấy giống hoa hồng tím và hoa trắng vàng tại Hà Giang giống nhau 100%, trong khi mẫu hoa hồng tím Thái Nguyên sai khác ở 2 vị trí nucleotide và mẫu hoa trắng đỏ sai khác ở 3 vị trí nucleotide. Sự biến đổi này dẫn tới thay đổi amino acid thứ 165 giữa Methionine, Serine và Cysteine.

Hoạt chất quý trong cây dừa cạn tập trung nhiều nhất ở bộ phận nào?

Hàm lượng alkaloid toàn phần biến thiên rõ rệt theo cơ quan thực vật. Rễ cây chứa tỷ lệ hoạt chất cao nhất đạt từ 0,7% đến 2,4%, tiếp theo là thân chứa khoảng 0,46% và lá dao động từ 0,37% đến 1,15%, trong đó lá là nguồn chính để chiết tách vinblastine và vincristine.

Phương pháp tách chiết CTAB 2% sử dụng trong nghiên cứu có ưu điểm gì?

Phương pháp CTAB 2% kết hợp nghiền mô trong nitơ lỏng âm 196 độ C giúp phá vỡ hoàn toàn vách tế bào thực vật, đồng thời kết tủa và loại bỏ hiệu quả các hợp chất thứ cấp, tanin và polysaccharide. Nhờ đó, sản phẩm DNA tổng số thu được có độ tinh sạch cao, tạo băng điện di sắc nét.

Khoảng cách di truyền 5,4% ở cấp độ nucleotide phản ánh điều gì?

Khoảng cách di truyền 5,4% chứng minh giữa các giống dừa cạn tại miền Bắc Việt Nam đã có sự phân hóa di truyền nhất định ở mức dưới loài do điều kiện địa lý và chọn lọc tự nhiên. Tuy nhiên, khoảng cách amino acid chỉ 0,4% cho thấy chức năng sinh học của protein vẫn được bảo toàn nguyên vẹn.

Kết luận

  • Đã xác định chi tiết các đặc điểm hình thái thực vật của 3 giống dừa cạn thu tại Hà Giang và Thái Nguyên, chỉ ra sự khác biệt đặc trưng về màu sắc cánh hoa, tâm hoa và sắc tố thân cây.
  • Tách chiết thành công DNA tổng số đạt độ tinh sạch cao và nhân bản đặc hiệu đoạn gen lục lạp rpoC1 kích thước 505 bp bằng cặp mồi rpoC1-1f và rpoC1-3r.
  • Xác định đầy đủ trình tự 505 nucleotide mã hóa 168 amino acid, phát hiện các điểm đột biến nucleotide tại vị trí 494, 495, 496 và sai khác amino acid tại vị trí 165.
  • Đánh giá mức độ tương đồng chuỗi đạt từ 96,6% đến 100%, xác định khoảng cách di truyền 5,4% ở cấp độ nucleotide và 0,4% ở cấp độ amino acid suy diễn.
  • Đóng góp nguồn dữ liệu trình tự chuẩn xác vào hệ thống mã vạch DNA quốc gia, tạo tiền đề khoa học vững chắc cho việc định danh dược liệu và bảo tồn nguồn gen quý.

Để tiếp tục hoàn thiện hệ thống cơ sở dữ liệu mã vạch phân tử, các nghiên cứu tiếp theo cần mở rộng phân tích các vùng gen matK và ITS trên các quần thể dừa cạn tại các vùng sinh thái khác nhau trong giai đoạn tới. Các tổ chức nghiên cứu và doanh nghiệp dược phẩm quan tâm đến quy trình phân tích di truyền hoặc hợp tác phát triển vùng trồng dược liệu chuẩn hóa có thể kết nối với nhóm tác giả qua kênh Học liệu Đại học Thái Nguyên để tiếp nhận chuyển giao kỹ thuật.