Tổng quan nghiên cứu

Khu hệ động vật có vú tại Việt Nam ghi nhận sự phát triển vượt bậc về số lượng loài, đặc biệt trong bộ Thú ăn sâu bọ (Eulipotyphla) với mức tăng từ 14 loài lên 34 loài. Trong đó, họ Chuột chũi (Talpidae) gồm 6 loài thuộc 3 giống là Euroscaptor, Mogera và Scaptonyx, nổi bật với 3 loài đặc hữu chỉ mới ghi nhận trên lãnh thổ Việt Nam. Mặc dù có tính đa dạng cao, việc phân loại các loài chuột chũi gặp nhiều thách thức lớn do hiện tượng tương đồng và chồng lấn hình thái bên ngoài như kích thước thân, chiều dài bàn chân hay màu sắc lông. Các nghiên cứu định danh trước đây phần lớn chỉ tập trung mô tả sơ bộ hoặc nghiên cứu biến thiên sọ trên phân bộ dơi và chuột chù, để lại khoảng trống dữ liệu đáng kể đối với họ Chuột chũi.

Mục tiêu trọng tâm của nghiên cứu là xây dựng cơ sở dữ liệu chi tiết gồm 16 chỉ tiêu hình thái sọ chuẩn hóa, đánh giá mức độ biến thiên kích thước và hình dạng hộp sọ trong cùng loài cũng như giữa các loài chuột chũi phân bố tại Việt Nam. Đồng thời, công trình bước đầu xác định ảnh hưởng của các yếu tố cách ly địa lý tự nhiên như hệ thống sông ngòi và các dãy núi cao đến sự sai khác hình thái sọ.

Nghiên cứu được triển khai trong giai đoạn từ tháng 1 năm 2016 đến tháng 10 năm 2018 trên phạm vi toàn quốc với 20 địa điểm thu mẫu, kết hợp khảo sát thực địa bổ sung tại Khu bảo tồn thiên nhiên Tây Côn Lĩnh có diện tích 15.406 ha ở tỉnh Hà Giang. Kết quả nghiên cứu cung cấp bộ số liệu đo đạc định lượng với độ tin cậy thống kê cao ở mức ý nghĩa p dưới 0,05, đóng góp nguồn tư liệu khoa học then chốt phục vụ công tác giám định loài, nghiên cứu tiến hóa và quản lý bảo tồn nguồn gen động vật quý hiếm.

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 hình thái sọ động vật có vú (Mammalian Cranial Morphometrics) kết hợp với học thuyết địa sinh thái quần thể. Ba khái niệm cốt lõi được vận dụng xuyên suốt bao gồm hiện tượng dị hình giới tính (sexual dimorphism), sự biến thiên hình thái nội loài (intraspecific variation) và sự phân hóa liên loài (interspecific divergence).

Mô hình nghiên cứu phân tích mối quan hệ tương hỗ giữa kích thước tuyệt đối và tỷ lệ hình học hộp sọ trước sức ép của các hàng rào cách ly địa lý tự nhiên. Các cấu trúc địa hình hiểm trở với hơn 70% diện tích đồi núi, cùng các mạng lưới thủy văn lớn như lưu vực sông Hồng chiếm 57,3% diện tích ngoài lãnh thổ hay hệ thống sông Chảy, tạo nên các sinh cảnh ngầm biệt lập. Điều này thúc đẩy quá trình thích nghi sinh thái và biến đổi hình thái sọ của các quần thể chuột chũi sống đào bới trong đất.

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

Nghiên cứu phân tích tổng thể 159 mẫu vật thuộc 5 loài và 2 giống chuột chũi thu thập tại 20 điểm sinh thái khác nhau. Trong đó, 155 mẫu sọ trưởng thành được kế thừa từ bộ sưu tập của Viện Sinh thái và Tài nguyên sinh vật, Bảo tàng Thiên nhiên Việt Nam cùng Viện Động vật Saint Petersburg thuộc Viện Hàn lâm Khoa học Nga, kết hợp với 4 mẫu vật thu thập bổ sung bằng bẫy chuyên dụng của Nhật Bản tại Hà Giang.

Quy trình chọn mẫu tập trung vào 122 mẫu sọ hoàn chỉnh tại 19 địa điểm để đo đạc trực tiếp 16 chỉ số hình thái học sọ tiêu chuẩn bằng thước kẹp điện tử có độ chính xác 0,01 mm, cùng 37 mẫu đối chứng quốc tế.

Phương pháp phân tích thành phần chính (PCA) và phân tích biến số chính tắc (CVA) trên phần mềm chuyên dụng PAST Statistics phiên bản 2.17 được lựa chọn nhờ khả năng giảm chiều dữ liệu từ 16 biến số mà vẫn giữ nguyên bản chất sai khác hình thái. PCA dữ liệu thô giúp đánh giá biến thiên kích thước sọ (PC1), trong khi PCA dữ liệu chuẩn hóa logarit với 1000 lần lặp giúp làm rõ biến thiên hình dạng sọ độc lập với độ lớn cơ thể. Sự sai khác thống kê giữa các giới tính và các quần thể cách ly được kiểm định chặt chẽ bằng phân tích phương sai One-way ANOVA kết hợp phép thử Tukey HSD và kiểm định t-test ở mức xác suất p dưới 0,05.

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

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

Phân tích thống kê đối với các chỉ số đo đạc khẳng định không có sự sai khác mang ý nghĩa thống kê giữa cá thể đực và cá thể cái ở các loài khảo sát với mức ý nghĩa p lớn hơn 0,05. Điều này cho phép gộp chung dữ liệu hai giới để phân tích hình thái quần thể một cách khách quan.

Trên không gian đồ thị PCA, trục thành phần chính PC1 giải thích hơn 70% tổng phương sai kích thước hộp sọ giữa 4 loài thuộc giống Euroscaptor. Loài chuột chũi răng nhỏ Euroscaptor parvidens có kích thước hộp sọ nhỏ nhất, trong khi loài chuột chũi or-lốp Euroscaptor orlovi thể hiện kích thước vượt trội hơn khoảng 25% ở chỉ số chiều dài sọ tối đa (GLS).

Sự biến thiên hình thái sọ nội loài diễn ra rõ nét theo không gian địa lý. Quần thể chuột chũi răng nhỏ Euroscaptor parvidens tại khu vực Nam Trung Bộ như Hòn Bà và Bi Doup có chiều dài hàng răng hàm trên (M1-M3) lớn hơn khoảng 12% so với quần thể tại Bắc Hướng Hóa thuộc Bắc Trung Bộ.

Nghiên cứu cũng ghi nhận sự phân hóa hình thái sâu sắc giữa các quần thể chuột chũi đuôi ngắn Euroscaptor subanura và chuột chũi miền bắc Mogera latouchei bị chia cắt bởi hệ thống sông Hồng và sông Lô, chứng minh vai trò rào cản di truyền của các dòng sông lớn.

Thảo luận kết quả

Sự biến đổi kích thước và hình dạng hộp sọ là kết quả của quá trình chọn lọc tự nhiên nhằm thích nghi với đặc tính cơ học của nền đất và thành phần con mồi tại các đai độ cao từ 500m đến trên 2000m. Các dữ liệu này được thể hiện trực quan qua biểu đồ phân tán hai chiều PCA và ma trận khoảng cách hình thái Mahalanobis trong không gian đa biến, giúp minh chứng rõ nét mức độ tách biệt của từng nhóm mẫu.

Tỷ lệ giữa khoảng cách hai ổ mắt trên chiều dài sọ tối đa (IOB/GLS) được chứng minh là chỉ số then chốt để phân định chính xác các loài có kích thước ngoài tương đồng. Kết quả này bổ sung bằng chứng vững chắc cho các kết luận phân loại phân tử trước đây, giải quyết triệt để các nhầm lẫn hình thái kéo dài nhiều thập kỷ.

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

Thứ nhất, ban hành bộ tiêu chuẩn nhận diện hình thái sọ với 16 chỉ số đo lường chuẩn hóa làm công cụ giám định chính thức cho các viện nghiên cứu và bảo tàng sinh học, hướng tới nâng cao độ chính xác phân loại đạt trên 98% trước quý 4 năm 2026, do Viện Sinh thái và Tài nguyên sinh vật chủ trì thực hiện.

Thứ hai, mở rộng tối thiểu 12 đợt điều tra thực địa chuyên sâu tại các khu vực núi cao hiểm trở thuộc Tây Bắc và Tây Nguyên trong giai đoạn 2026 đến 2028 nhằm thu thập mẫu chuẩn của loài chuột chũi đuôi dài Scaptonyx fusicaudus còn thiếu, do Ban quản lý các Vườn quốc gia phối hợp với các chuyên gia động vật học triển khai.

Thứ ba, thiết lập cơ sở dữ liệu số hóa tích hợp chỉ số hình thái sọ với mã vạch di truyền DNA Barcoding từ các đoạn gen ty thể cho 100% mẫu vật lưu trữ, hoàn thành ngân hàng dữ liệu quốc gia trước năm 2027 dưới sự quản lý của Bảo tàng Thiên nhiên Việt Nam.

Thứ tư, quy hoạch và bảo vệ nghiêm ngặt tối thiểu 5 hành lang sinh thái ngầm liên kết các khu rừng tự nhiên tại vùng Tây Côn Lĩnh và Hoàng Liên Sơn trước năm 2029 nhằm giảm thiểu tác động phân mảnh sinh cảnh, do Cục Lâm nghiệp phối hợp với chính quyền địa phương thực thi.

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

Các nhà khoa học và chuyên gia phân loại học động vật có vú có thể khai thác bộ 16 chỉ số hình thái sọ cùng phương pháp phân tích đa biến PCA để giải quyết các vướng mắc danh pháp trong phân loại thú nhỏ.

Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Động vật học và Sinh thái học có thể sử dụng công trình như một tài liệu chuẩn mực về phương pháp xử lý số liệu hình thái học sọ và kỹ thuật làm tiêu bản động vật có xương sống.

Cán bộ quản lý tại các Vườn quốc gia và Khu bảo tồn thiên nhiên có thể ứng dụng các dẫn liệu phân bố địa lý để cập nhật danh mục động vật đặc hữu và xây dựng kế hoạch quản lý nguồn gen nguy cấp tại cơ sở.

Chuyên gia tư vấn đánh giá tác động môi trường và đa dạng sinh học có thể sử dụng đặc tính biến thiên quần thể chuột chũi làm sinh vật chỉ thị đánh giá mức độ chia cắt sinh cảnh đất trong các dự án phát triển hạ tầng.

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

Vì sao việc định loại chuột chũi qua kích thước bên ngoài thường không mang lại độ chính xác cao? Nhiều loài chuột chũi có sự tương đồng lớn về kích thước thân, chiều dài đuôi và bàn chân do quá trình tiến hóa hội tụ để thích nghi với đời sống ngầm. Việc đo đạc 16 chỉ số hình thái hộp sọ cung cấp các đặc điểm sai khác cấu trúc xương ổn định, giúp định danh chính xác tuyệt đối.

Phân tích thành phần chính PCA đóng vai trò như thế nào trong nghiên cứu hình thái sọ? Phương pháp PCA giúp giảm 16 chiều dữ liệu hình thái xuống các trục thành phần chính độc lập. Trục PC1 phản ánh biến thiên kích thước tổng thể chiếm trên 70% phương sai, trong khi các trục chuẩn hóa tiếp theo thể hiện sự khác biệt về hình dạng hộp sọ độc lập với độ lớn mẫu vật.

Nghiên cứu xử lý yếu tố sai khác giữa cá thể đực và cá thể cái như thế nào? Toàn bộ mẫu vật được kiểm định phương sai F-test và t-test trước khi phân tích chuyên sâu. Kết quả chứng minh giá trị p đều lớn hơn 0,05, khẳng định không tồn tại hiện tượng dị hình giới tính đối với kích thước hộp sọ, cho phép xử lý mẫu gộp một cách khoa học.

Hàng rào địa lý tự nhiên tại Việt Nam tác động ra sao đến các loài chuột chũi? Các hệ thống sông lớn như sông Hồng, sông Lô cùng các dãy núi cao trên 2000m đóng vai trò rào cản ngăn chặn sự giao phối chéo, dẫn đến sự phân hóa hình thái sọ rõ rệt với mức sai khác kích thước hàng răng lên tới 12% giữa các quần thể cách ly.

Độ bao phủ của nguồn mẫu vật trong luận văn có bảo đảm tính đại diện không? Nghiên cứu khảo sát 159 mẫu vật của 5 loài thu thập từ 20 địa điểm sinh thái trải dài từ miền Bắc đến Tây Nguyên, kết hợp mẫu chuẩn quốc tế từ Nga và mẫu thực địa tại Hà Giang, bảo đảm tính đại diện cao cho khu hệ thú Việt Nam.

Kết luận

• Luận văn xây dựng thành công bộ cơ sở dữ liệu 16 chỉ số hình thái hộp sọ chuẩn mực trên 159 mẫu vật thuộc 5 loài chuột chũi tại Việt Nam.

• Chứng minh tính hiệu quả vượt trội của phương pháp phân tích đa biến PCA và CVA trong việc tách biệt các loài có hình thái bên ngoài tương đồng.

• Khẳng định không có dị hình giới tính về hình thái sọ và chỉ ra sự biến thiên kích thước sọ nội loài lên đến 12-25% do thích nghi sinh thái.

• Xác lập cơ sở khoa học về tác động phân hóa hình thái của các hàng rào cách ly địa lý tự nhiên như sông ngòi và các dãy núi cao.

• Bổ sung nguồn mẫu vật thực địa quý giá cho Khu bảo tồn thiên nhiên Tây Côn Lĩnh và hệ thống các bảo tàng sinh học quốc gia.

Trong giai đoạn 2026 đến 2028, các hướng nghiên cứu tiếp theo cần mở rộng thu thập mẫu phân tử và xây dựng bản đồ phân bố không gian chi tiết cho toàn bộ bộ Thú ăn sâu bọ. Các cơ quan quản lý và nhà nghiên cứu hãy khai thác ngay nguồn cơ sở dữ liệu hình thái học này để nâng cao hiệu quả công tác bảo tồn đa dạng sinh học tại Việt Nam.