Tổng quan nghiên cứu

Việt Nam nằm trong trung tâm của khu vực sinh thái Đông Dương - Mi-an-ma (Indo-Burma), một trong 25 điểm nóng về đa dạng sinh học toàn cầu với khu hệ bò sát phong phú bậc nhất. Số lượng loài bò sát ghi nhận tại Việt Nam tăng trưởng nhanh chóng từ 258 loài năm 1996 lên 368 loài năm 2009 và vượt mốc 430 loài vào năm 2018. Riêng phân bộ Rắn (Serpentes) đã ghi nhận sự gia tăng mạnh mẽ từ 171 loài năm 2005 lên 244 loài vào năm 2018. Trong đó, giống Rắn khuyết (Lycodon Fitzinger, 1826) thuộc họ Rắn nước (Colubridae) là một trong những nhóm đa dạng nhất với 58 loài trên thế giới, 35 loài ở Đông Nam Á và 16 loài ở Đông Dương. Tuy nhiên, các loài thuộc giống này có hình thái tương đồng cao, cấu trúc vảy phức tạp và lịch sử phân loại biến động khi nhiều giống như Dinodon, Cercaspis, Dryocalamus và Lepturophis được sáp nhập vào Lycodon, dẫn đến nhiều sai sót trong định loại thực địa.

Mục tiêu trọng tâm của nghiên cứu là xác định toàn diện thành phần loài Rắn khuyết ở Việt Nam, cập nhật vùng phân bố, mô tả chi tiết đặc điểm hình thái học và làm sáng tỏ quan hệ di truyền phát sinh chủng loại thông qua chỉ thị sinh học phân tử. Phạm vi nghiên cứu bao quát dữ liệu mẫu vật thu thập tại 27 tỉnh thành từ Bắc vào Nam, kết hợp 5 đợt khảo sát thực địa chuyên sâu với 42 ngày khảo sát tại Nghệ An, Thừa Thiên Huế và Cà Mau trong giai đoạn từ tháng 01/2017 đến tháng 05/2018. Kết quả nghiên cứu có ý nghĩa khoa học và thực tiễn sâu sắc khi hoàn thiện cơ sở dữ liệu phân loại học, cung cấp bằng chứng phục vụ công tác quy hoạch, đánh giá tình trạng bảo tồn các loài bò sát đặc hữu và nguy cấp 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 được xây dựng trên nền tảng phân loại học tích hợp (Integrative Taxonomy), kết hợp chặt chẽ giữa giải phẫu hình thái học truyền thống và phát sinh chủng loại học phân tử (Molecular Phylogenetics). Khung lý thuyết phân loại học bò sát hiện đại áp dụng hệ thống phân loại của họ Rắn nước (Colubridae) và phân họ Colubrinae theo các cập nhật quốc tế mới nhất, tích hợp các giống đồng danh vào giống Rắn khuyết Lycodon.

Bốn khái niệm khoa học cốt lõi được ứng dụng xuyên suốt bao gồm:

  1. Đa dạng hình thái và chỉ số sai khác pholidosis (hệ thống vảy đầu, thân, bụng và huyệt);
  2. Khoảng cách di truyền (Genetic Distance) phản ánh tỷ lệ đột biến nucleotide giữa các bậc phân loại;
  3. Cây quan hệ di truyền theo phương pháp suy luận Bayes (Bayesian Inference) với độ tin cậy được đánh giá qua xác suất hậu nghiệm (Posterior Probability ≥ 95%);
  4. Đơn vị phân loại học tiến hóa độc lập (Evolutionary Significant Units) nhằm phát hiện các loài ẩn kê và quần thể phân hóa địa lý.

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

Nghiên cứu tiến hành phân tích tổng thể trên cỡ mẫu 145 mẫu vật hình thái lưu giữ tại các bảo tàng khoa học lớn như Viện Sinh thái và Tài nguyên sinh vật (IEBR), Bảo tàng Thiên nhiên Việt Nam (VNMN), Bảo tàng Alexander Koenig Bonn (ZFMK) và các trường đại học tại 27 tỉnh thành. Phương pháp chọn mẫu thực địa áp dụng kỹ thuật điều tra theo tuyến sinh cảnh ven suối, vách đá, hang động và thảm mục rừng vào ban đêm (khung giờ từ 18:00 đến 24:00) với tọa độ được định vị bằng máy GPS Garmin Dakota 20.

Phân tích hình thái được thực hiện thông qua đo đếm 59 chỉ tiêu tiêu chuẩn (11 chỉ số đo kích thước bằng thước kẹp điện tử độ chính xác 0,01 mm, 39 chỉ số đếm vảy và 9 chỉ tiêu hoa văn màu sắc). Đối với phân tích phân tử, 13 mẫu cơ đại diện được tách chiết DNA tổng số bằng bộ kit DNeasy Blood & Tissue (Qiagen), sau đó khuếch đại đoạn gen ty thể Cytochrome b (Cyt-b) bằng phản ứng PCR với cặp mồi chuyên biệt L14910 và H16064 qua 35 chu kỳ nhiệt. Việc lựa chọn gen ty thể Cyt-b là phương pháp tối ưu vì đoạn gen này có tốc độ tiến hóa phù hợp, cung cấp độ phân giải cao để phân định loài và mối quan hệ di truyền giữa các quần thể bò sát. Dữ liệu trình tự được gióng hàng bằng ClustalX, xác định mô hình tiến hóa tối ưu qua Modeltest v3.7, xây dựng cây phân loại bằng MrBayes v3.2 chạy trong 10 triệu thế hệ và tính toán khoảng cách di truyền bằng phần mềm PAUP v4.

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 tổng cộng 15 loài Rắn khuyết thuộc giống Lycodon phân bố tại Việt Nam, bao gồm 13 loài đã được định danh chính thức và 2 đơn vị phân loại chưa định danh thuộc nhóm phức hợp loài Lycodon fasciatus (gồm Lycodon sp1 thu tại Quảng Nam, Kon Tum và Lycodon sp2 thu tại Khánh Hòa). Đáng chú ý, nghiên cứu phát hiện và bổ sung địa điểm phân bố mới cho 11 loài (chiếm 73,3% tổng số loài ghi nhận), mở rộng đáng kể ranh giới địa lý sinh vật học của giống này tại Việt Nam.

Một đóng góp nổi bật là lần đầu tiên cung cấp đầy đủ dữ liệu hình thái học của cá thể cái đối với loài Rắn khuyết ít đốm (Lycodon paucifasciatus) – loài đặc hữu cực kỳ quý hiếm được xếp hạng Sẽ nguy cấp (VU) trong Danh lục Đỏ IUCN 2017. Các ghi nhận phân bố mới tiêu biểu gồm có: loài Lycodon capucinus được phát hiện tại Phú Yên và đảo Hòn Khoai (Cà Mau); loài Lycodon cardamomensis mở rộng tới Quảng Ngãi và đảo Phú Quốc (Kiên Giang); loài Lycodon futsingensis được xác nhận thêm tại 5 tỉnh miền Bắc và miền Trung (Cao Bằng, Tuyên Quang, Sơn La, Hòa Bình, Nghệ An, Quảng Nam); loài Lycodon davisonii mở rộng thêm 4 tỉnh (Hòa Bình, Thừa Thiên Huế, Bình Định, Cà Mau).

Thảo luận kết quả

Kết quả so sánh hình thái đã chỉ ra sự biến thiên địa lý rõ rệt giữa các quần thể. Điển hình ở loài Lycodon fasciatus, quần thể miền Bắc có số lượng khoanh sáng nhiều hơn, bản khoanh hẹp và có đốm đen ở giữa, trong khi quần thể miền Trung sở hữu các khoanh sáng rộng, màu trắng thuần. Về giới hạn sinh thái, nghiên cứu đã mở rộng biên độ độ cao phân bố của nhiều loài: Lycodon capucinus được ghi nhận ở độ cao thấp kỷ lục là 20 m so với mực nước biển, Lycodon futsingensis đạt giới hạn trên tới 996 m, và Lycodon septentrionalis ghi nhận ở vùng núi cao từ 1480 m đến 1730 m.

Phân tích di truyền gen Cyt-b củng cố vững chắc quan điểm chuyển các loài thuộc giống Dinodon và Dryocalamus vào giống Lycodon, thể hiện qua các nhánh đơn ngành có giá trị xác suất hậu nghiệm Bayesian đạt 95% đến 100%. Các dữ liệu hình thái chi tiết của 15 loài được chuẩn hóa thành ma trận bảng biểu so sánh 59 chỉ số đo đếm, kết hợp cùng sơ đồ cây phát sinh chủng loại Bayes và bản đồ phân bố GIS trực quan trên 27 tỉnh thành, tạo nên một hệ thống dẫn liệu khoa học hoàn chỉnh và có độ tin cậy vượt trội so với các công bố đơn lẻ trước đây.

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

Dựa trên kết quả nghiên cứu, 4 giải pháp hành động cụ thể được đề xuất nhằm nâng cao hiệu quả nghiên cứu và bảo tồn đa dạng sinh học bò sát:

  1. Thiết lập chương trình giám sát sinh cảnh định kỳ cho loài đặc hữu nguy cấp Lycodon paucifasciatus tại Vườn quốc gia Phong Nha - Kẻ Bàng và Khu bảo tồn thiên nhiên Động Châu - Khe Nước Trong. Mục tiêu duy trì 100% diện tích sinh cảnh tự nhiên nguyên vẹn, thực hiện từ quý 1/2024 đến năm 2026 bởi Ban quản lý các khu bảo tồn phối hợp cùng Chi cục Kiểm lâm địa phương.
  2. Triển khai phân tích sinh học phân tử đa gen đối với 2 taxon tiềm năng chưa định danh (Lycodon sp1 và Lycodon sp2). Tăng cường cỡ mẫu tối thiểu 30 mẫu vật kết hợp giải trình tự các gen nhân như c-mos, RAG1 và gen ty thể ND4 trong giai đoạn 2024 - 2025 do Viện Sinh thái và Tài nguyên sinh vật chủ trì để mô tả loài mới nếu đủ bằng chứng.
  3. Xây dựng bản đồ phân bố không gian và cơ sở dữ liệu số GIS tích hợp chỉ số sinh khí hậu cho 15 loài Rắn khuyết trên phạm vi 27 tỉnh thành. Dự án đặt mục tiêu số hóa 100% điểm tọa độ thực địa, hoàn thành cập nhật hàng năm từ năm 2025 dưới sự quản lý của Cục Bảo tồn thiên nhiên và Đa dạng sinh học.
  4. Tổ chức các khóa tập huấn chuyên sâu về nhận dạng hình thái học thực địa và kỹ thuật cứu hộ bò sát cho hơn 500 cán bộ kiểm lâm tại các vườn quốc gia trọng điểm (Pù Mát, Sao La, Bidoup - Núi Bà) trong giai đoạn 2024 - 2027 nhằm giảm thiểu nhầm lẫn giữa rắn vô hại và rắn độc.

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

Luận văn là tài liệu tham khảo học thuật chuyên sâu và có giá trị ứng dụng cao cho 4 nhóm đối tượng chính:

  1. Các nhà nghiên cứu bò sát học (Herpetologists) và chuyên gia phân loại học động vật: Sử dụng bảng số liệu 59 chỉ tiêu hình thái và trình tự gen chuẩn Cyt-b làm hệ quy chiếu đối sánh khi giám định các mẫu vật phức tạp thuộc họ Colubridae.
  2. Ban quản lý các Vườn quốc gia, Khu bảo tồn thiên nhiên và lực lượng Kiểm lâm: Ứng dụng dữ liệu ghi nhận phân bố mới tại 27 tỉnh để cập nhật danh lục đa dạng sinh học địa phương và xây dựng kế hoạch tuần tra, bảo vệ sinh cảnh rừng.
  3. Giảng viên, học viên cao học và sinh viên các ngành Sinh học, Động vật học, Lâm nghiệp và Quản lý tài nguyên môi trường: Tham khảo cấu trúc phương pháp luận phân loại tích hợp (Integrative Taxonomy) chuẩn mực từ khâu điều tra thực địa đến phân tích số liệu hình thái và phát sinh chủng loại phân tử.
  4. Các tổ chức bảo tồn thiên nhiên trong nước và quốc tế (IUCN, WAR, SVW): Sử dụng dẫn liệu thực tế về loài đặc hữu Lycodon paucifasciatus và 2 loài chưa định danh để xây dựng hồ sơ đánh giá tình trạng bảo tồn và vận động chính sách bảo vệ động vật hoang dã.

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

Giống Rắn khuyết Lycodon có những đặc điểm giải phẫu hình thái đặc trưng nào?
Giống Lycodon đặc trưng bởi xương hàm trên uốn cong vào trong, răng hàm trên chia làm 3 nhóm cách nhau bởi khoảng trống: nhóm đầu 3-6 răng lớn dần, nhóm giữa 7-15 răng nhỏ và nhóm cuối 2-3 răng lớn. Vảy thân gồm 13-19 hàng, con ngươi hình bầu dục đứng, chiều dài toàn bộ cơ thể dao động từ 350 mm đến hơn 1700 mm tùy loài.

Nghiên cứu đã ghi nhận bao nhiêu loài Rắn khuyết tại Việt Nam và có phát hiện gì mới?
Nghiên cứu ghi nhận 15 loài Lycodon (gồm 13 loài định danh và 2 loài chưa định danh thuộc nhóm Lycodon fasciatus). Đề tài đã bổ sung vùng phân bố mới cho 11 loài tại 27 tỉnh thành, đồng thời lần đầu tiên mô tả chi tiết hình thái cá thể cái của loài đặc hữu Lycodon paucifasciatus.

Tại sao cần kết hợp cả đặc điểm hình thái và phân tích chỉ thị gen ty thể Cyt-b?
Do các loài Lycodon có hình thái ngoài rất giống nhau và biến thiên phức tạp giữa các quần thể nên phương pháp đo đếm hình thái đơn thuần dễ dẫn đến định danh sai. Phân tích đoạn gen Cyt-b (dài khoảng 1100 cặp base) giúp cung cấp bằng chứng khách quan, làm sáng tỏ các mối quan hệ phát sinh chủng loại và xác thực các nhóm loài phức hợp.

Quy mô thu thập dữ liệu và khảo sát thực địa của đề tài được tiến hành ra sao?
Đề tài thực hiện 5 đợt điều tra thực địa với 42 ngày khảo sát đêm tại Nghệ An, Thừa Thiên Huế và Cà Mau, đồng thời nghiên cứu trực tiếp 145 mẫu vật lưu trữ từ 27 tỉnh thành trên cả nước cùng 13 mẫu phân tích sinh học phân tử đạt độ chuẩn xác cao.

Loài Rắn khuyết ít đốm Lycodon paucifasciatus có giá trị bảo tồn như thế nào?
Đây là loài rắn đặc hữu của Việt Nam, được xếp ở bậc Sẽ nguy cấp (VU) theo Danh lục Đỏ IUCN 2017. Nghiên cứu đã cung cấp các dẫn liệu bổ sung quan trọng về hình thái học và mở rộng vùng phân bố của loài tại Hà Tĩnh và Quảng Bình, làm căn cứ khoa học thiết yếu cho công tác bảo tồn loài này.

Kết luận

  • Xác định và cập nhật danh lục 15 loài Rắn khuyết thuộc giống Lycodon tại Việt Nam, bao gồm 2 taxon tiềm năng chưa định danh khoa học.
  • Bổ sung địa điểm ghi nhận phân bố địa lý mới cho 11 loài trên phạm vi 27 tỉnh thành và lần đầu tiên công bố dữ liệu hình thái cá thể cái của loài đặc hữu Lycodon paucifasciatus.
  • Phân tích chi tiết 59 chỉ tiêu đo đếm hình thái trên 145 mẫu vật, ghi nhận nhiều dẫn liệu mới về sinh thái học và mở rộng biên độ độ cao sinh cảnh từ 20 m đến 1730 m.
  • Giải trình tự và ứng dụng thành công chỉ thị gen ty thể Cyt-b để xây dựng cây phát sinh chủng loại Bayes, củng cố vị trí phân loại học của giống Lycodon mở rộng.
  • Cung cấp luận cứ khoa học thực tiễn vững chắc phục vụ chiến lược giám sát và quy hoạch bảo tồn đa dạng sinh học bò sát giai đoạn 2024 - 2028.

Công trình là bước tiến quan trọng trong nghiên cứu khu hệ bò sát Việt Nam. Quý độc giả, nhà khoa học và các đơn vị bảo tồn có thể khai thác bộ dữ liệu chuẩn hóa của luận văn để phục vụ công tác nghiên cứu chuyên sâu và quản lý tài nguyên sinh vật hiệu quả.