Tổng quan nghiên cứu

Việt Nam sở hữu mức độ đa dạng sinh học cao với 1.181 loài bướm ngày thuộc 6 họ chính đã được ghi nhận. Tuy nhiên, dữ liệu di truyền học phục vụ công tác định danh còn rất hạn chế khi Hệ thống dữ liệu mã vạch của sự sống (BOLD) tính đến tháng 7 năm 2021 mới chỉ lưu trữ 80 trình tự của 45 loài, tương đương 3,81% tổng số loài trên toàn quốc. Tại Cao nguyên đá Đồng Văn (tỉnh Hà Giang) – khu vực có diện tích 2.356,8 km² với hơn 186 loài bướm và tỷ lệ đặc hữu vượt trên 16% – nguồn cơ sở dữ liệu mã vạch trước đây chỉ vỏn vẹn 7 trình tự thuộc 6 loài.

Phương pháp phân loại học hình thái truyền thống thường gặp nhiều rào cản khi xử lý các mẫu vật hư hại, các cặp loài đồng hình hoặc các cá thể biến đổi hình thái theo độ cao địa lý. Nhằm giải quyết triệt để vấn đề này, luận văn thạc sĩ chuyên ngành Sinh học thực nghiệm tập trung thực hiện ba mục tiêu cụ thể:

  • Giải mã thành công trình tự nucleotide vùng gen Cytochrome c oxidase 1 (CO1) cho 28 mẫu bướm đại diện 14 loài thuộc 5 họ thu thập tại Lũng Cú, Đồng Văn, Hà Giang ở độ cao 1.600 mét.
  • Đăng ký thành công toàn bộ 28 trình tự lên Ngân hàng dữ liệu gen quốc tế GenBank, mở rộng độ phủ dữ liệu di truyền bướm Việt Nam.
  • Đánh giá mức độ tương đồng di truyền, kiểm tra tính chính xác của phương pháp định danh hình thái và làm rõ quan hệ phát sinh chủng loại đối với các loài có hình thái phức tạp.

Nghiên cứu được triển khai thực địa từ ngày 20 đến ngày 25 tháng 5 năm 2019 và hoàn thiện phân tích di truyền vào năm 2021. Kết quả đóng góp trực tiếp 28 mã vạch chuẩn, tăng độ bao phủ mã vạch DNA cho khu hệ bướm ngày Hà Giang lên 400%, tạo nền tảng vững chắc cho công tác giám định loài và bảo tồn sinh học vùng núi cao Đông Bắc.

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 mã vạch DNA (DNA Barcoding) làm nền tảng cốt lõi để nhận diện và phân loại sinh vật. Trong giới động vật, đoạn gen chuẩn Cytochrome c oxidase 1 (CO1) thuộc hệ gen ti thể có chiều dài tiêu chuẩn 658 cặp bazơ (bp) được quốc tế công nhận là công cụ tối ưu nhờ tính bảo thủ cao nhưng vẫn tích lũy đủ các biến dị nucleotide để phân biệt loài.

Nghiên cứu tích hợp mô hình phân tích phát sinh chủng loại phân tử (Molecular Phylogenetics) dựa trên 4 khái niệm trung tâm:

  • Mã vạch DNA (DNA Barcode): Trình tự nucleotide ngắn, chuẩn hóa dùng để định danh loài nhanh chóng.
  • Ngưỡng sai khác liên loài (Interspecific Divergence): Mức độ khác biệt nucleotide giữa các loài khác nhau trong cùng một giống, thường được xác lập ở ngưỡng từ 3,0% trở lên đối với nhóm côn trùng cánh vẩy.
  • Ngưỡng biến dị nội loài (Intraspecific Divergence): Tỷ lệ sai khác di truyền giữa các cá thể trong cùng một loài, thường dao động dưới 1,0% (từ 0 đến 2 nucleotide).
  • Phân tích nhánh ngoài (Outgroup Comparison) và giá trị lặp lại Bootstrap: Phương pháp xác thực độ tin cậy của cấu trúc cây tiến hóa, trong đó các nhánh đạt chỉ số bootstrap từ 70% trở lên sau 1.000 lần lặp được xem là có ý nghĩa phân loại vững chắc.

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

Nguồn dữ liệu của luận văn được xây dựng từ cỡ mẫu gồm 28 cá thể bướm trưởng thành đại diện cho 14 loài thuộc 5 họ: Hesperiidae, Lycaenidae, Nymphalidae, Pieridae và Riodinidae (mỗi loài gồm 2 mẫu đối chứng). Phương pháp chọn mẫu có chủ đích được áp dụng nhằm ưu tiên các nhóm loài đặc trưng cho sinh cảnh rừng núi cao 1.600 mét, các loài có dạng hình thái dị hình giữa cá thể đực và cái, các loài chưa từng có dữ liệu mã vạch tại Việt Nam, và các đơn vị phân loại còn tồn tại tranh cãi học thuật.

Quy trình thực nghiệm được chuẩn hóa nghiêm ngặt qua 4 bước:

  • Tách chiết DNA tổng số: Sử dụng bộ kit DNeasy Blood and Tissue (hãng Qiagen) trên mô cơ chân bướm đã được nghiền cơ học để phá vỡ lớp vỏ kitin, đảm bảo hiệu suất thu hồi DNA tinh sạch mà không làm tổn hại cấu trúc tiêu bản cánh.
  • Phản ứng PCR và giải trình tự: Sử dụng cặp mồi phổ quát LCO1490 và HCO2198 để khuếch đại vùng gen CO1 qua 35 chu kỳ nhiệt (nhiệt độ gắn mồi 45°C trong 40 giây, kéo dài ở 72°C trong 60 giây). Sản phẩm PCR sau điện di thạch agarose 0,9% được tinh sạch bằng bộ kit QIAquick và giải trình tự hai chiều trên hệ thống tự động của Macrogen.
  • Xử lý và phân tích tin sinh học: Kiểm tra sắc ký đồ trên phần mềm ChromasPro 2.10 để loại bỏ tín hiệu nhiễu. Sử dụng thuật toán BLAST và ClustalW trên phần mềm MEGA X để căn gióng chuỗi nucleotide với 45 mẫu đối sánh truy xuất từ cơ sở dữ liệu BOLD.
  • Lý do lựa chọn phương pháp phân tích: Cây phát sinh chủng loại được xây dựng theo phương pháp Tiết kiệm tối đa (Maximum Parsimony – MP) với 1.000 lần lặp bootstrap. Phương pháp MP được lựa chọn vì tính khách quan cao, tốc độ xử lý nhanh và nguyên lý tối ưu hóa số lượng biến đổi đột biến ít nhất, rất phù hợp cho việc phân tích các chuỗi gen ti thể có tính bảo thủ cao ở côn trùng. Toàn bộ lộ trình nghiên cứu thực địa và phân tích phòng thí nghiệm được hoàn tất trong mốc thời gian từ tháng 5 năm 2019 đến tháng 7 năm 2021.

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 đã mang lại những kết quả khoa học nổi bật với các chỉ số đo lường chính xác:

  • Hiệu suất khuếch đại và giải trình tự đạt 100%: Toàn bộ 28/28 mẫu vật đều cho sản phẩm PCR chất lượng cao, các đỉnh huỳnh quang sắc nét. Có 22 mẫu đạt độ dài chuẩn 658 bp và 6 mẫu thuộc giống Dodona đạt chiều dài 600 bp. Toàn bộ 28 trình tự đã được GenBank cấp mã định danh quốc tế chính thức từ OK342112 đến OK342274 và OK340914 đến OK340919.
  • Độ đồng nhất di truyền nội loài đạt từ 99,85% đến 100%: Các mẫu trong cùng một loài thể hiện sự tương đồng gần như tuyệt đối. Điển hình như loài Symbrenthia lilaea chỉ xuất hiện 1 đột biến điểm tại vị trí nucleotide 41 (tương đồng 99,85%), loài Pieris canidia chỉ sai khác 1 nucleotide tại vị trí 196 (tương đồng 99,85%), trong khi 12 loài còn lại có mức độ trùng khớp 100% giữa 2 cá thể nghiên cứu.
  • Ngưỡng phân tách liên loài rõ ràng vượt mức 3,0%: Mẫu Sebastonyma dolopia thể hiện mức sai khác 3,04% với 20 vị trí nucleotide biến dị khi so sánh với loài họ hàng gần Sebastonyma medoensis (mã số MN199421 trên BOLD). Đặc biệt, mẫu Dodona dipoea có mức sai biệt lên tới 5,17% với 31 vị trí nucleotide đột biến so với loài đối sánh gần nhất Dodona elvira (mã số KX866692).
  • Bổ sung dữ liệu đầu tiên cho họ Riodinidae tại Việt Nam: Nghiên cứu đã công bố thành công mã vạch cho 3 loài thuộc họ Riodinidae gồm Dodona dipoea, Dodona eugenes và Dodona ouida, xóa bỏ hoàn toàn khoảng trống dữ liệu 0% của họ bướm này trên hệ sinh thái BOLD Việt Nam.

Thảo luận kết quả

Cây phát sinh chủng loại tổng thể được thiết lập từ 42 trình tự (bao gồm mẫu nghiên cứu và các chuỗi đối chiếu quốc tế) cùng 1 mẫu ngoài Acherontia lachesis (mã MT080189) cho thấy tỷ lệ phân biệt các loài đạt mức 100% với giá trị ủng hộ bootstrap dao động từ 91% đến 100%. Dữ liệu ma trận khoảng cách di truyền và sơ đồ cây tiến hóa hình quạt được trình bày trực quan qua các bảng so sánh nucleotide và biểu đồ phát sinh, phản ánh chính xác ranh giới phân loại sinh học.

Sự tách biệt di truyền rõ rệt của khu hệ bướm tại Lũng Cú (Hà Giang) được giải thích bởi yếu tố sinh thái và địa hình karst đặc thù. Độ cao 1.600 mét cùng các dãy núi đá vôi dốc đứng đóng vai trò như các hàng rào địa lý tự nhiên, thúc đẩy quá trình cách ly sinh sản và lưu giữ các biến dị di truyền đặc trưng.

Về mặt phân loại học, nghiên cứu đã giải quyết thành công những nghi vấn kéo dài về phân loài:

  • Đối với loài Curetis bulis, mẫu thu tại Hà Giang (độ cao 1.600 mét) có độ sai khác di truyền 0,91% (6 nucleotide biến dị) so với mẫu phân loài Curetis bulis stigmata tại vùng đất thấp Malaysia (mã KF491662). Cây phát sinh chủng loại phân nhánh rõ rệt với độ tin cậy 100%, khẳng định mẫu Hà Giang thuộc phân loài danh định Curetis bulis bulis phân bố tại vùng núi cao Đông Dương.
  • Đối với giống Dodona và Dilipa, các cây phát sinh chủng loại đơn nhánh đã chứng minh tính độc lập di truyền hoàn toàn của Dodona eugenes (trùng khớp 100% với mẫu KX866727) và Dilipa morgiana (trùng khớp 100% với mẫu AB501204), loại bỏ nguy cơ định danh nhầm do biến thái màu sắc cánh.

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

Dựa trên các kết quả đạt được, luận văn đề xuất 4 nhóm giải pháp mang tính hành động nhằm nâng cao năng lực nghiên cứu và ứng dụng mã vạch DNA tại Việt Nam:

  • Mở rộng phạm vi giải trình tự mã vạch CO1 toàn diện: Tiến hành thu thập và giải mã trình tự cho 100% trong tổng số 186 loài bướm tại Công viên Địa chất Toàn cầu Cao nguyên đá Đồng Văn. Mục tiêu đến năm 2026 hoàn thành ngân hàng mã vạch chuẩn cho ít nhất 150 loài bướm vùng Đông Bắc do Bảo tàng Thiên nhiên Việt Nam và các viện nghiên cứu chuyên ngành chủ trì thực hiện.
  • Ứng dụng kỹ thuật đa chỉ thị phân tử (Multi-locus Barcoding): Kết hợp đoạn gen ti thể CO1 với 2 chỉ thị gen nhân có tốc độ tiến hóa chậm như EF-1a và vùng phiên mã nội bộ ITS. Giải pháp này hướng tới mục tiêu nâng độ phân giải phân loại lên trên 99,5% đối với các phức hợp loài khó định danh, thực hiện định kỳ trong giai đoạn 2025–2028.
  • Chuẩn hóa quy trình tách chiết DNA không xâm lấn: Triển khai áp dụng quy trình trích xuất DNA từ một phần chân hoặc mô tế bào nhỏ cho 100% các bảo tàng sinh học và phòng thí nghiệm trên toàn quốc từ năm 2025. Quy trình này bảo tồn nguyên vẹn hơn 95% giá trị hình thái của các mẫu vật quý hiếm và mẫu vật lưu trữ lịch sử.
  • Xây dựng cổng dữ liệu mã vạch DNA sinh vật quốc gia: Thiết lập nền tảng số liên thông dữ liệu giữa Viện Hàn lâm Khoa học và Công nghệ Việt Nam với các hệ thống GenBank và BOLD quốc tế. Đặt chỉ tiêu số hóa 1.181 loài bướm ngày Việt Nam trong lộ trình 2025–2030 dưới sự điều phối của Bộ Khoa học và Công nghệ.

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

Nội dung luận văn mang giá trị học thuật sâu sắc và tính ứng dụng thực tiễn cao cho 4 nhóm đối tượng:

  • Các nhà nghiên cứu phân loại học và côn trùng học: Cung cấp phương pháp luận kết hợp giữa hình thái học và sinh học phân tử, giúp giải quyết các bài toán về loài đồng hình, xác định chính xác danh tính các loài phụ phức tạp thuộc các giống Dodona, Curetis và Pieris.
  • Ban quản lý các khu bảo tồn thiên nhiên và Vườn quốc gia: Cung cấp công cụ giám sát đa dạng sinh học nhanh chóng thông qua kỹ thuật barcode, hỗ trợ đánh giá tác động của biến đổi khí hậu lên các loài côn trùng chỉ thị tại vùng núi cao.
  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Công nghệ sinh học: Là tài liệu tham khảo chi tiết về quy trình thực nghiệm PCR, phân tích tin sinh học trên phần mềm MEGA X, cách xử lý cây phát sinh chủng loại Maximum Parsimony và đăng ký dữ liệu chuẩn hóa lên GenBank.
  • Cơ quan kiểm định môi trường và cơ quan thực thi bảo tồn: Giúp nhận diện chính xác các loài bướm quý hiếm được liệt kê trong Sách đỏ Việt Nam và Công ước CITES (như Troides aeacus) từ các mẫu vật không còn nguyên vẹn hoặc các sản phẩm thương mại bị phân mảnh.

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

Vì sao đoạn gen CO1 ti thể được chọn làm mã vạch chuẩn cho bướm ngày?

Vùng gen CO1 có số lượng bản sao rất lớn trong mỗi tế bào, cấu trúc bền vững và sở hữu tốc độ đột biến nucleotide phù hợp. Chiều dài tiêu chuẩn 658 bp cung cấp đủ tín hiệu tiến hóa để phân biệt các loài khác nhau với độ chính xác trên 95%, đồng thời các cặp mồi phổ quát như LCO1490/HCO2198 cho tỷ lệ khuếch đại thành công gần như tuyệt đối ở côn trùng.

Luận văn đã đóng góp bao nhiêu trình tự mới cho hệ thống dữ liệu quốc tế?

Nghiên cứu đã cung cấp 28 trình tự vùng gen CO1 hoàn chỉnh đại diện cho 14 loài thuộc 5 họ bướm ngày thu thập tại Hà Giang. Toàn bộ 28 trình tự này đã được phê duyệt và gắn mã truy cập chính thức trên Ngân hàng dữ liệu gen GenBank (từ OK340914 đến OK342274), bổ sung nguồn dữ liệu quý giá cho hệ thống BOLD toàn cầu.

Ngưỡng sai khác di truyền nào được sử dụng để phân định các loài bướm?

Khoa học phân tử áp dụng ngưỡng sai khác nucleotide liên loài từ 3,0% trở lên để xác định ranh giới loài riêng biệt ở bướm ngày. Luận văn đã chứng minh rõ điều này khi mẫu Sebastonyma dolopia có mức sai khác 3,04% so với Sebastonyma medoensis, và Dodona dipoea sai khác 5,17% so với Dodona elvira, hoàn toàn vượt ngưỡng biến dị phân loài thông thường.

Việc tách chiết DNA có làm hư hại mẫu vật bướm trong bộ sưu tập không?

Nghiên cứu áp dụng quy trình xử lý mô vi lượng, chỉ cắt một đoạn chân rất nhỏ của mẫu bướm để phân rã mô cơ. Phương pháp này đảm bảo thu nhận nồng độ DNA ổn định phục vụ phản ứng PCR, đồng thời bảo tồn nguyên vẹn 100% cấu trúc cánh, thân và các đặc điểm hình thái quan trọng khác của mẫu tiêu bản bảo tàng.

Nghiên cứu đã làm sáng tỏ tình trạng phân loại của loài nào tại Hà Giang?

Nghiên cứu đã làm rõ vị trí phân loài của Curetis bulis bulis thu tại độ cao 1.600 mét ở Hà Giang khi phát hiện mức sai khác di truyền 0,91% so với phân loài Curetis bulis stigmata ở Malaysia. Đồng thời, nghiên cứu cũng phân lập thành công các loài trong giống Dodona và Dilipa morgiana vốn có đặc điểm bên ngoài dễ gây nhầm lẫn.

Kết luận

  • Đã giải mã thành công 100% trình tự vùng gen mã vạch CO1 cho 28 mẫu bướm đại diện 14 loài thuộc 5 họ tại huyện Đồng Văn, tỉnh Hà Giang.
  • Đăng ký thành công 28 mã định danh quốc tế trên GenBank (từ OK340914 đến OK342274), góp phần tăng tỷ lệ bao phủ mã vạch DNA cho bướm ngày Việt Nam.
  • Xác lập độ biến dị nội loài cực thấp (dưới 0,15%) và độ sai khác liên loài rõ rệt (từ 3,04% đến 5,17%), chứng minh tính đặc hiệu vượt trội của chỉ thị CO1.
  • Bổ sung dữ liệu di truyền đầu tiên cho họ Riodinidae tại Việt Nam và giải quyết dứt điểm các vướng mắc về danh pháp phân loài đối với loài Curetis bulis.
  • Hoàn thiện quy trình phân tích phát sinh chủng loại phân tử chuẩn xác, tạo tiền đề khoa học vững chắc để nhân rộng cho các chương trình đánh giá đa dạng sinh học toàn quốc giai đoạn 2025–2030.

Hãy tiếp tục khám phá và tích hợp công nghệ giải mã gen vào công tác bảo tồn đa dạng sinh học để bảo vệ các nguồn gen quý giá của tự nhiên.