Tổng quan nghiên cứu

Việt Nam nằm trọn trong 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 tại nước ta liên tục gia tăng nhanh chóng qua các mốc nghiên cứu, từ 258 loài năm 1996 lên 296 loài năm 2005, 368 loài năm 2009 và đạt khoảng 430 loài vào năm 2018. Riêng phân bộ rắn (Serpentes) đã ghi nhận sự mở rộng vượt bậc từ 171 loài năm 2005 lên 244 loài vào tháng 5 năm 2018. Trong bức tranh đa dạng sinh học đó, giống rắn khuyết (Lycodon Fitzinger, 1826) thuộc họ Rắn nước (Colubridae) đại diện cho một trong những nhóm đa dạng nhất với 58 loài trên thế giới, phân bố trải rộng từ Đông Iran đến Nam Trung Quốc và Đông Nam Á.

Tuy nhiên, các loài thuộc giống rắn khuyết sở hữu đặc điểm hình thái tương đồng cao, cấu trúc vảy phức tạp và hiện tượng biến dị nội loài lớn, gây ra nhiều nhầm lẫn trong phân loại học truyền thống. Nhằm giải quyết triệt để vấn đề này, luận văn thạc sĩ chuyên ngành Động vật học tập trung thực hiện hai mục tiêu cốt lõi: cập nhật thành phần loài, phân tích hình thái chi tiết và giải mã mối quan hệ di truyền phát sinh loài của giống rắn khuyết tại Việt Nam. Phạm vi nghiên cứu bao quát 145 mẫu vật hình thái và 13 mẫu phân tử thu thập tại 27 tỉnh thành trên cả nước, kết hợp 5 đợt khảo sát thực địa bổ sung với 42 ngày nghiên cứu thực tế tại Nghệ An, Thừa Thiên Huế và Cà Mau. Công trình cung cấp cơ sở dữ liệu vảy học chuẩn mực và phân tử đáng tin cậy phục vụ công tác quy hoạch, quản lý và bảo tồn bền vững nguồn tài nguyên sinh vật 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

Luận văn vận dụng nền tảng của lý thuyết phân loại học tích hợp (Integrative Taxonomy), kết hợp chặt chẽ giữa hình thái học so sánh truyền thống và phát sinh chủng loại học phân tử hiện đại. Khung nghiên cứu áp dụng lý thuyết tiến hóa phân tử và đồng hồ sinh học để đánh giá mức độ sai khác di truyền trên chuỗi gen ty thể Cytochrome b (cyt-b), khắc phục hoàn toàn những giới hạn mà phương pháp phân loại hình thái đơn thuần gặp phải đối với các loài có hình thái tương đồng (cryptic species).

Mô hình nghiên cứu tập trung vào 4 khái niệm khoa học trọng tâm:

  1. Giống rắn khuyết (Lycodon Fitzinger, 1826): Nhóm rắn nước đặc trưng bởi xương hàm trên cong về phía trước, răng hàm trên chia thành 3 nhóm kích thước riêng biệt với khoảng trống ở giữa, mắt tròn có con ngươi hình bầu dục đứng và vảy thân từ 15 đến 19 hàng.
  2. Quan hệ di truyền phát sinh loài (Phylogenetic relationships): Cấu trúc phả hệ thể hiện mức độ phân nhánh tiến hóa và khoảng cách di truyền giữa các đơn vị phân loại.
  3. Loài đặc hữu (Endemic species): Nhóm sinh vật chỉ phân bố tự nhiên trong phạm vi địa lý giới hạn tại Việt Nam.
  4. Mã vạch di truyền ty thể (Mitochondrial DNA Barcoding): Sử dụng đoạn gen cyt-b có kích thước 1117 cặp bazơ làm chỉ thị phân tử nhận diện loài.

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

Nghiên cứu sử dụng tổng cộng 145 mẫu vật hình thái lưu trữ tại 8 bảo tàng khoa học uy tín trong nước và quốc tế (Viện Sinh thái và Tài nguyên sinh vật, Bảo tàng Thiên nhiên Việt Nam, Bảo tàng Alexander Koenig Đức) cùng 13 mẫu mô phục vụ phân tích sinh học phân tử từ 27 tỉnh thành trải dài từ Bắc vào Nam. Phương pháp chọn mẫu thực địa là chọn mẫu có chủ đích dọc theo các sinh cảnh đặc trưng như ven suối, vách đá, đường mòn rừng và gốc cây mục. Khảo sát được thực hiện chủ yếu vào ban đêm từ 18:00 đến 24:00 – khung giờ hoạt động mạnh nhất của các loài rắn khuyết, kết hợp định vị tọa độ bằng máy GPS Garmin Dakota 20.

Lý do lựa chọn phương pháp đo đếm hình thái 59 chỉ tiêu vảy học (chiều dài cơ thể SVL, tỷ lệ đuôi TaL/TL, số hàng vảy bụng, vảy lưng và vảy đầu) bằng thước kẹp điện tử độ chính xác 0,01 mm là nhằm định lượng hóa tối đa các biến dị giải phẫu. Về phân tích phân tử, DNA tổng số được tách chiết bằng bộ kit DNeasy Blood and Tissue của Qiagen, khuếch đại gen cyt-b qua phản ứng PCR với cặp mồi chuyên biệt L14910 và H16064 qua 35 chu kỳ nhiệt. Sản phẩm được tinh sạch bằng GeneJET và giải trình tự hai chiều tại Firstbase (Malaysia). Cây phát sinh chủng loại được xây dựng theo mô hình tiến hóa tối ưu từ ModelTest v3.7 và suy luận Bayesian qua 10.000.000 thế hệ trên phần mềm MrBayes v3.2, đảm bảo độ tin cậy nhánh với xác suất hậu nghiệm trên 95% và tính toán khoảng cách di truyền qua 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 mang lại 4 phát hiện khoa học mang tính đột phá về khu hệ rắn khuyết Việt Nam:

Thứ nhất, xác định và cập nhật danh lục 15 loài rắn khuyết Lycodon phân bố tại Việt Nam, chiếm khoảng 25,9% tổng số 58 loài rắn khuyết được ghi nhận trên toàn thế giới. Đáng chú ý, nghiên cứu phát hiện 2 loài chưa định danh thuộc nhóm Lycodon fasciatus gồm Lycodon sp1 thu tại Quảng Nam, Kon Tum và Lycodon sp2 thu tại Khánh Hòa với các đặc điểm hình thái và di truyền sai khác rõ rệt.

Thứ hai, mở rộng vùng phân bố địa lý mới cho 11 loài rắn khuyết (chiếm tỷ lệ 73,3% tổng số loài của giống tại Việt Nam). Nổi bật là loài Lycodon futsingensis được bổ sung vùng phân bố tại 8 tỉnh (Cao Bằng, Tuyên Quang, Sơn La, Hòa Bình, Nghệ An, Quảng Nam, Kon Tum, Lâm Đồng) và loài Lycodon fasciatus được ghi nhận mới tại 5 tỉnh miền Trung và Tây Nguyên.

Thứ ba, lần đầu tiên phát hiện và công bố dẫn liệu hình thái cá thể cái của loài rắn khuyết ít đốm (Lycodon paucifasciatus) – loài đặc hữu cực kỳ quý hiếm được xếp ở bậc Sẽ nguy cấp (VU) trong Danh lục Đỏ IUCN 2017, đồng thời mở rộng vùng phân bố của loài này tới tỉnh Hà Tĩnh.

Thứ tư, ghi nhận giới hạn phân bố độ cao thấp nhất của loài Lycodon capucinus ở độ cao chỉ từ 20 đến 60 mét tại Phú Yên và đảo Hòn Khoai (Cà Mau), phá vỡ nhận định trước đây cho rằng loài này chỉ phân bố ở độ cao trên 500 mét; đồng thời ghi nhận bằng chứng sinh sản với 5 trứng trong buồng trứng cá thể cái và phổ thức ăn chứa tắc kè đuôi vàng (Cnemaspis psychedelica).

Thảo luận kết quả

Sự phân hóa phong phú của giống rắn khuyết tại Việt Nam bắt nguồn từ điều kiện tự nhiên đặc thù với 7 kiểu khí hậu nhiệt đới gió mùa và mạng lưới thủy văn phức tạp, nơi lưu vực các sông suối lớn chiếm 59,7% dòng chảy từ bên ngoài vào. Địa hình chia cắt sâu sắc bởi dãy Trường Sơn tạo thành các ranh giới địa lý sinh vật học biệt lập, thúc đẩy quá trình phân loài địa phương.

So sánh với các nghiên cứu phân loại trước đây của Smith (1943), Nguyễn Văn Sáng và cộng sự (2005) hay Vogel và cộng sự (2009), kết quả nghiên cứu đã giải quyết thỏa đáng các tranh cãi học thuật kéo dài. Phân tích phân tử gen cyt-b ủng hộ quan điểm sáp nhập các giống Dinodon và Dryocalamus vào giống Lycodon với mức độ sai khác di truyền giữa các nhánh dao động từ 2,5% đến hơn 14,8%. Để tối ưu hóa việc phân tích dữ liệu, các dẫn liệu phân loại được trực quan hóa qua biểu đồ cây phát sinh chủng loại Bayesian với xác suất hậu nghiệm đạt trên 95% và bảng đối sánh 59 chỉ tiêu vảy học, giúp phân định rõ ràng các loài chị em có hình thái tương đồng cao trong giống Lycodon.

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

Dựa trên các kết quả nghiên cứu thực nghiệm, luận văn đề xuất 4 nhóm giải pháp hành động cụ thể nhằm nâng cao hiệu quả quản lý và bảo tồn:

  1. Khoanh vùng và thiết lập chế độ bảo vệ nghiêm ngặt sinh cảnh cốt lõi:
  • Hành động: Thành lập các khu vực ưu tiên bảo tồn sinh cảnh rừng thường xanh núi đá vôi và rừng ẩm nhiệt đới.
  • Chỉ số mục tiêu: Bảo vệ nguyên vẹn 100% sinh cảnh sống của loài rắn khuyết ít đốm (Lycodon paucifasciatus) và 2 loài chưa định danh tại các khu vực trọng điểm.
  • Thời gian thực hiện: Giai đoạn 2024–2026.
  • Chủ thể thực hiện: Chi cục Kiểm lâm các tỉnh Hà Tĩnh, Thừa Thiên Huế, Quảng Nam, Kon Tum phối hợp với Ban quản lý các Vườn Quốc gia và Khu bảo tồn thiên nhiên.
  1. Mở rộng nghiên cứu di truyền phân tử đa gen để công bố loài mới:
  • Hành động: Tiến hành giải trình tự bổ sung các gen nhân (như c-mos, RAG1) kết hợp gen ty thể ND4 nhằm xác định vị trí phân loại chính xác.
  • Chỉ số mục tiêu: Phân tích 100% mẫu vật của hai dạng loài Lycodon sp1 và Lycodon sp2 để hoàn thiện hồ sơ công bố loài mới cho khoa học.
  • Thời gian thực hiện: Hoàn thành trước năm 2026.
  • Chủ thể thực hiện: Viện Sinh thái và Tài nguyên sinh vật chủ trì, phối hợp với Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội.
  1. Xây dựng sổ tay nhận dạng thực địa và chuẩn hóa quy trình phân loại vảy học:
  • Hành động: Biên tập và phát hành bộ tài liệu hướng dẫn tra cứu nhanh dựa trên 59 chỉ tiêu hình thái và khóa định loại chuẩn.
  • Chỉ số mục tiêu: Tổ chức các lớp tập huấn kỹ thuật cho ít nhất 150 cán bộ kiểm lâm và nhân viên bảo tồn tại 27 tỉnh có loài phân bố.
  • Thời gian thực hiện: Giai đoạn 2025–2027.
  • Chủ thể thực hiện: Bộ Nông nghiệp và Phát triển Nông thôn phối hợp với các viện nghiên cứu chuyên ngành.
  1. Ứng dụng hệ thống thông tin địa lý (GIS) trong quan trắc biến động quần thể:
  • Hành động: Xây dựng cơ sở dữ liệu số hóa toàn diện về tọa độ, sinh cảnh và điều kiện sinh khí hậu của giống Lycodon.
  • Chỉ số mục tiêu: Số hóa 100% dữ liệu phân bố và mô hình hóa nguy cơ suy giảm sinh cảnh dưới tác động của biến đổi khí hậu.
  • Thời gian thực hiện: Tầm nhìn đến năm 2030.
  • Chủ thể thực hiện: Bộ Khoa học và Công nghệ tài trợ kinh phí thông qua Quỹ Phát triển khoa học và công nghệ Quốc gia (NAFOSTED).

Đố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 có giá trị thực tiễn cao dành cho 4 nhóm đối tượng chính:

  1. Nhà nghiên cứu và giảng viên chuyên ngành Động vật học, Sinh học tiến hóa:
  • Lợi ích và ứng dụng: Cung cấp nguồn tài liệu chuẩn mực về 59 chỉ tiêu hình thái và quy trình phân tích phát sinh chủng loại phân tử gen cyt-b, phục vụ trực tiếp cho công tác giảng dạy sau đại học và phát triển các công trình nghiên cứu so sánh khu hệ bò sát Đông Nam Á.
  1. Ban quản lý các Vườn Quốc gia và Khu bảo tồn thiên nhiên:
  • Lợi ích và ứng dụng: Ứng dụng dữ liệu phân bố mới của 11 loài rắn khuyết tại 27 tỉnh thành để rà soát danh mục đa dạng sinh học cơ sở, từ đó xây dựng phương án tuần tra rừng và bảo vệ các sinh cảnh đặc thù có sự xuất hiện của các loài quý hiếm.
  1. Cơ quan quản lý nhà nước và lực lượng Kiểm lâm:
  • Lợi ích và ứng dụng: Sử dụng khóa định loại hình thái chính xác để hỗ trợ công tác giám định loài trong thực thi pháp luật, kiểm soát buôn bán động vật hoang dã và tham mưu cập nhật Sách Đỏ Việt Nam cũng như Danh mục các loài nguy cấp, quý, hiếm.
  1. Học viên cao học, sinh viên và các tổ chức phi chính phủ bảo tồn thiên nhiên:
  • Lợi ích và ứng dụng: Tiếp cận phương pháp luận phân loại học tích hợp (Integrative Taxonomy) chuẩn mực, học hỏi quy trình xử lý dữ liệu Bayesian với 10 triệu thế hệ để ứng dụng cho các đề tài nghiên cứu khoa học và dự án bảo tồn thực địa.

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

Luận văn đã xác định được bao nhiêu loài rắn khuyết thuộc giống Lycodon tại Việt Nam?

Nghiên cứu đã ghi nhận và cập nhật danh lục gồm 15 loài rắn khuyết tại Việt Nam, chiếm khoảng 25,9% trong tổng số 58 loài của giống này trên toàn cầu. Danh lục bao gồm 13 loài đã định danh chính xác và 2 loài chưa định danh thuộc nhóm Lycodon fasciatus tại khu vực miền Trung và Nam Trung Bộ.

Kỹ thuật sinh học phân tử nào được ứng dụng để giải mã quan hệ di truyền của giống Lycodon?

Luận văn sử dụng phương pháp tách chiết DNA tổng số, khuếch đại đoạn gen ty thể Cytochrome b với cặp mồi chuyên biệt L14910 và H16064 qua 35 chu kỳ nhiệt trên 13 mẫu đại diện. Phân tích phát sinh chủng loại được xây dựng theo phương pháp suy luận Bayesian qua 10.000.000 thế hệ trên phần mềm MrBayes v3.2, đạt xác suất hậu nghiệm trên 95%.

Nghiên cứu đã bổ sung vùng phân bố địa lý mới cho những loài nào?

Tác giả đã ghi nhận vùng phân bố mới cho 11 loài rắn khuyết tại 27 tỉnh thành trên cả nước. Tiêu biểu như loài Lycodon futsingensis được phát hiện thêm tại 8 tỉnh từ Bắc vào Nam, loài Lycodon fasciatus tại 5 tỉnh miền Trung, và loài Lycodon capucinus lần đầu được ghi nhận tại tỉnh Phú Yên và đảo Hòn Khoai thuộc tỉnh Cà Mau.

Loài rắn khuyết nào trong nghiên cứu có ý nghĩa bảo tồn đặc biệt cấp thiết?

Loài rắn khuyết ít đốm (Lycodon paucifasciatus) là loài đặc hữu của Việt Nam, được xếp ở bậc Sẽ nguy cấp (VU) trong Danh lục Đỏ IUCN năm 2017. Nghiên cứu đã lần đầu tiên mô tả chi tiết đặc điểm hình thái của cá thể cái và mở rộng vùng phân bố của loài tới tỉnh Hà Tĩnh.

Tác giả đã sử dụng bao nhiêu chỉ tiêu hình thái để định loại các loài rắn khuyết?

Nghiên cứu đã tiến hành đo đếm và phân tích chi tiết 59 chỉ tiêu hình thái (bao gồm kích thước cơ thể, tỷ lệ đuôi, số hàng vảy thân, vảy bụng, vảy dưới đuôi và các tấm vảy đầu) trên 145 mẫu vật bằng thước kẹp điện tử có độ chính xác 0,01 mm nhằm đảm bảo tính chuẩn xác tuyệt đối trong định danh.

Kết luận

  • Xác lập 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, chiếm tỷ lệ 25,9% tổng số loài của giống trên phạm vi toàn cầu.
  • Mở rộng vùng phân bố địa lý cho 11 loài rắn khuyết tại 27 tỉnh thành trên toàn quốc qua 42 ngày khảo sát thực địa chuyên sâu.
  • Bổ sung dẫn liệu hình thái học hoàn chỉnh đầu tiên của cá thể cái loài đặc hữu nguy cấp Lycodon paucifasciatus với hệ thống 59 chỉ tiêu vảy học.
  • Giải mã quan hệ di truyền phát sinh loài của 13 mẫu vật dựa trên đoạn gen ty thể cyt-b qua 10.000.000 thế hệ phân tích Bayesian.
  • Phát hiện 2 loài chưa định danh khoa học đầy đủ, định hình hướng nghiên cứu phân loại chuyên sâu trong giai đoạn 2024–2026.

Lời kêu gọi hành động: Các cơ quan quản lý tài nguyên môi trường, viện nghiên cứu và các tổ chức bảo tồn cần nhanh chóng liên kết nguồn lực, ứng dụng dữ liệu phân loại học tích hợp của luận văn để triển khai các biện pháp bảo tồn sinh cảnh thực địa, góp phần bảo vệ vững chắc nguồn gen bò sát đặc hữu và quý hiếm của Việt Nam.