Tổng quan luận án

Đề tài "Nghiên cứu đặc điểm sinh thái học và di truyền của Thằn lằn bóng đốm Eutropis macularius (Blyth, 1853) ở vùng Cao nguyên Buôn Ma Thuột – Buôn Hồ" do nghiên cứu sinh Trương Bá Phong thực hiện dưới sự hướng dẫn khoa học của GS. Ngô Đắc Chứng và cố PGS. Ngô Văn Bình, thuộc chuyên ngành Sinh thái học (mã số: 9420120), bảo vệ tại Trường Đại học Đà Lạt năm 2023.

Tính cấp thiết và khoảng trống nghiên cứu

Lớp Bò sát trên thế giới ghi nhận 10.940 loài tính đến tháng 12 năm 2022. Theo đánh giá của Liên minh Bảo tồn Thiên nhiên Quốc tế (IUCN) và công bố của Cox cùng cộng sự, khoảng 21% số loài bò sát toàn cầu đang đối mặt với nguy cơ tuyệt chủng do suy thoái môi trường sống, biến đổi khí hậu và áp lực khai thác tài nguyên. Việt Nam sở hữu mức độ đa dạng sinh học cao với khoảng 650 loài bò sát được ghi nhận đến năm 2016. Khu vực Tây Nguyên, đặc biệt là vùng Cao nguyên Buôn Ma Thuột – Buôn Hồ (tỉnh Đắk Lắk), nằm trong vùng khí hậu nhiệt đới gió mùa mưa hè với mùa khô kéo dài 3–4 tháng, có sinh cảnh đặc trưng là rừng khộp rụng lá theo mùa và các đồn điền cây công nghiệp. Tuy nhiên, hệ sinh thái này đang chịu áp lực lớn từ nạn phá rừng, canh tác nông nghiệp, xây dựng thủy điện và săn bắt động vật hoang dã.

Thằn lằn bóng đốm Eutropis macularius (Blyth, 1853) thuộc họ Scincidae, bộ Squamata, là loài bò sát ăn côn trùng, đóng vai trò quan trọng trong việc cân bằng chuỗi thức ăn tự nhiên và hỗ trợ kiểm soát sinh vật gây hại cho nông, lâm nghiệp. Mặc dù vậy, các nghiên cứu trước đây tại Tây Nguyên chủ yếu dừng lại ở mức độ điều tra danh lục thành phần loài. Trước thời điểm luận án thực hiện, chưa có công trình nào nghiên cứu toàn diện và có hệ thống về đặc điểm hình thái, sai khác giới tính, sinh thái học quần thể (mật độ, sử dụng vi môi trường, mô hình chiếm cứ điểm), sinh thái dinh dưỡng, sinh thái sinh sản kết hợp với cấu trúc di truyền phân tử của loài này tại vùng Cao nguyên Buôn Ma Thuột – Buôn Hồ cũng như toàn vùng Tây Nguyên.

Mục tiêu và nội dung nghiên cứu

Mục tiêu tổng quát của luận án là bổ sung và hệ thống hóa cơ sở dữ liệu khoa học về đặc điểm hình thái, sinh thái học và đa dạng di truyền của loài Thằn lằn bóng đốm Eutropis macularius tại vùng Cao nguyên Buôn Ma Thuột – Buôn Hồ, làm cơ sở khoa học cho công tác bảo tồn và quản lý loài.

Luận án xác định 4 nội dung nghiên cứu cụ thể được đánh số theo văn bản:

  1. Xác định đặc điểm hình thái và phân tích tương quan giữa những sai khác về hình thái theo giới tính.
  2. Phân tích các đặc điểm sinh thái học: mật độ quần thể; sử dụng vi môi trường sống; xác suất phát hiện loài và tỉ suất chiếm cứ điểm.
  3. Phân tích đặc điểm sinh học về dinh dưỡng và sinh sản.
  4. Đánh giá mức độ đa dạng di truyền và so sánh với các vùng khác ở khu vực Tây Nguyên.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Loài Thằn lằn bóng đốm Eutropis macularius (Blyth, 1853).
  • Phạm vi không gian: Vùng Cao nguyên Buôn Ma Thuột – Buôn Hồ (tỉnh Đắk Lắk), tập trung khảo sát sinh thái tại Vườn quốc gia (VQG) Yok Don, Thành phố Buôn Ma Thuột, huyện Buôn Đôn và huyện Krông Búk. Phạm vi nghiên cứu di truyền được mở rộng đối sánh trên 4 tỉnh Tây Nguyên gồm Đắk Lắk, Đắk Nông, Gia Lai và Kon Tum.
  • Phạm vi thời gian: Toàn bộ dữ liệu thực địa và phân tích thực nghiệm được tiến hành liên tục trong giai đoạn từ năm 2017 đến năm 2021.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước đây

  1. Phân loại học và tên gọi:

    • Trên thế giới, Blyth (1853) mô tả loài lần đầu tiên dưới tên Euprepis macularius. Günther (1864) mô tả Euprepes rufescens và năm 1875 đổi thành Eumeces brevis. Anderson (1871) dùng tên Euprepes macularius, Boulenger (1887) đổi thành Mabuia macularia. Mausfeld và cộng sự (2002), Mausfeld & Schmitz (2003) phân tích quan hệ di truyền và chuyển các loài thằn lằn bóng châu Á sang giống Eutropis. Murthy (2010) chuẩn hóa danh pháp Eutropis macularius theo Luật Danh pháp Quốc tế (ICZN).
    • Tại Việt Nam, Bourret (1934), Trần Kiên, Nguyễn Văn Sáng, Hồ Thu Cúc (1981), Ngô Đắc Chứng (1998, 2007, 2008), Lê Nguyên Ngật (2001, 2002) từng sử dụng tên Mabuya macularia. Từ năm 2009, các tác giả Trương Thị Vinh Hương, Lê Nguyên Ngật (2009), Hoàng Xuân Quang, Ngô Đắc Chứng (2009, 2012) bắt đầu sử dụng tên Eutropis macularia và chuẩn hóa thành Eutropis macularius.

  2. Đặc điểm hình thái và sai khác giới tính (SSD):

    • Nghiên cứu hình thái quốc tế do Blyth (1853) và Boulenger (1887) đặt nền móng. Anderson (1994), Brana (1996), Schwarzkopf (2005) nghiên cứu cơ chế chọn lọc giới tính. Olsson và cộng sự (2002) nghiên cứu trên Niveoscincus microlepidotus, Huang (2006a, 2006b) nghiên cứu trên Eutropis longicaudata chỉ ra sự khác biệt kích thước đầu và ổ bụng giữa con đực và con cái. Cox, Skelly, Alder (2003) xây dựng công thức tính chỉ số sai khác kích thước giới tính (SSD) dựa trên chiều dài thân (SVL).
    • Tại Việt Nam, Hoàng Ngọc Thảo, Hoàng Xuân Quang, Ngô Đắc Chứng (2009, 2012, 2013) mô tả hình thái loài tại VQG Bạch Mã và Bắc Trung Bộ. Ngô Đắc Chứng, Ngô Văn Bình (2014) khảo sát sai khác giới tính trên loài Eutropis multifasciatus.

  3. Phân bố, sinh thái học, dinh dưỡng và sinh sản:

    • Phân bố loài được ghi nhận tại Nam Á và Đông Nam Á (Smith, 1935; Taylor, 1963; Grismer et al., 2008; Adil et al., 2020). Tại Việt Nam, loài được ghi nhận ở Cát Tiên (Nguyễn Văn Sáng, Hồ Thu Cúc, 2002), Khánh Hòa (Ngô Đắc Chứng, Trương Tấn Mỹ, 2007), Đắk Nông (Trương Thị Vinh Hương, Lê Nguyên Ngật, 2009), Thừa Thiên Huế (Lê Thắng Lợi, Ngô Đắc Chứng, 2009; Ngô Đắc Chứng et al., 2020), Quảng Trị, An Giang, Đồng Tháp (Ngô Đắc Chứng et al., 2012).
    • Sinh thái học mật độ và vi môi trường được nghiên cứu bởi Bullock & Evans (1990), Vitt et al. (1997), Van Schingen et al. (2014). Phương pháp đánh giá xác suất phát hiện loài và mô hình chiếm cứ điểm phát triển bởi MacKenzie và cộng sự (2002, 2003, 2004, 2006), Cao Thị Thanh Nguyên, Ngô Văn Bình, Ngô Đắc Chứng (2018, 2019).
    • Về dinh dưỡng và sinh sản, Huang và cộng sự (2006a, 2006b, 2011, 2012) khảo sát tập tính trên Eutropis longicaudata; Ngô Đắc Chứng và cộng sự (2007, 2014, 2015, 2020) nghiên cứu dinh dưỡng và chu kỳ sinh sản của Eutropis multifasciatus và Eutropis macularius tại miền Trung.

  4. Nghiên cứu di truyền:

    • Ota và cộng sự (1996, 2001) khảo sát bộ nhiễm sắc thể (2n = 32 hoặc 2n = 38). Kỹ thuật RAPD được dùng trên bò sát (Williams et al., 1990; Lougheed et al., 2000; Hoang, 2018). Phân tích trình tự gen ty thể 12S và 16S rRNA được thực hiện bởi Honda et al. (1999, 2000), Mausfeld (2000), Whiting et al. (2006), Datta-Roy et al. (2012).

Khoảng trống luận án giải quyết

Luận án giải quyết khoảng trống dữ liệu thực chứng tại vùng Cao nguyên Buôn Ma Thuột – Buôn Hồ: cung cấp bộ số liệu đo đạc hình thái chi tiết kèm kiểm định sai khác giới tính; lượng hóa mật độ quần thể và xác định các loại vi môi trường sống đặc trưng; áp dụng mô hình toán học chiếm cứ điểm theo không gian và mùa; phân tích toàn diện quang phổ thức ăn (chỉ số IRI) và chu kỳ sinh sản thông qua giải phẫu lát cắt mô học; đồng thời giải trình tự và giải đoán cấu trúc di truyền quần thể dựa trên gen ty thể 16S rDNA tại 4 tỉnh Tây Nguyên.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và khung phân tích

  • Lý thuyết sai khác kích thước giới tính (Sexual Size Dimorphism - SSD): Đánh giá mức độ phân hóa kích thước giữa hai giới theo chỉ số của Cox, Skelly và Alder (2003): $$\text{SSD} = \left(\frac{\text{kích thước giới lớn hơn}}{\text{kích thước giới nhỏ hơn}}\right) - 1$$ (quy ước mang dấu âm nếu con đực lớn hơn con cái và dấu dương nếu con cái lớn hơn con đực).
  • Lý thuyết chiếm cứ điểm (Site Occupancy Modeling): Sử dụng khung phân tích xác suất của MacKenzie và cộng sự (2002, 2006), đánh giá tỉ suất chiếm cứ điểm ($\psi$) và xác suất phát hiện loài ($p$) thông qua tiêu chuẩn thông tin Akaike (AIC) nhằm chọn lọc mô hình tối ưu có tính đến tác động của sinh cảnh và điều kiện vi khí hậu.
  • Lý thuyết sinh thái dinh dưỡng: Lượng hóa vai trò con mồi qua Chỉ số quan trọng tương đối (Index of Relative Importance - IRI) của Pinkas et al. (1971): $$\text{IRI} = (%N + %V) \times %F$$ kết hợp đánh giá độ phong phú và độ đồng đều thành phần thức ăn thông qua chỉ số đa dạng Simpson ($1/D$).
  • Lý thuyết di truyền quần thể và kiểm định tính trung lập: Đo lường các chỉ số đa dạng haplotype ($Hd$), đa dạng nucleotide ($\pi$), khoảng cách di truyền, hệ số phân hóa quần thể ($F_{st}$) và các chỉ số kiểm định trung lập (Tajima's $D$, Fu & Li's $D^$ và $F^$).

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

  • Khảo sát thực địa và thu mẫu: Thiết lập 72 ô tiêu chuẩn tại VQG Yok Don (chia thành vùng lõi và vùng đệm, sinh cảnh rừng khộp và rừng trồng). Bắt mẫu bằng tay hoặc bẫy hố, ghi nhận tọa độ GPS, đo đạc nhiệt độ, độ ẩm môi trường bằng máy đo chuyên dụng.
  • Phân tích hình thái: Đo 16 chỉ tiêu kích thước cơ thể bằng thước kẹp điện tử (độ chính xác 0,01 mm) và cân điện tử (độ chính xác 0,01 g), bao gồm: chiều dài thân (SVL), chiều dài đuôi (TaL), chiều dài đầu (HL), chiều rộng đầu (HW), chiều rộng miệng (MW), khối lượng cơ thể (BM)...
  • Phân tích dinh dưỡng: Sử dụng phương pháp rửa dạ dày không gây chết hoặc mổ dạ dày các mẫu thu thập; định danh con mồi dưới kính hiển vi soi nổi đến cấp bộ hoặc họ; đếm số lượng ($N$), đo thể tích ($V$) bằng phương pháp dìm nước và tính tần suất xuất hiện ($F$).
  • Phân tích sinh sản và cấu trúc mô học: Đo kích thước, khối lượng tuyến sinh dục (tinh hoàn, buồng trứng) và gan; phân loại 3 giai đoạn phát triển của trứng; làm tiêu bản mô học tinh hoàn và buồng trứng bằng phương pháp cắt lát nhuộm Hematoxylin - Eosin (H&E) để quan sát các giai đoạn phát triển của ống sinh tinh và noãn hoàng.
  • Kỹ thuật sinh học phân tử: Tách chiết DNA tổng số từ mô cơ/đuôi; nhân bản đoạn gen ty thể 16S rDNA bằng kỹ thuật PCR; giải trình tự nucleotide trực tiếp; xử lý và căn chỉnh chuỗi bằng phần mềm BioEdit, so sánh đối chiếu với chuỗi chuẩn mã số AB057394 trên GenBank; xây dựng cây phả hệ di truyền và tính toán các chỉ số di truyền bằng phần mềm MEGA và DnaSP.
  • Xử lý số liệu thống kê: Dữ liệu được phân tích trên phần mềm SPSS, R và PRESENCE; sử dụng phân tích phương sai (ANOVA), kiểm định $t$-test, phân tích hồi quy tương quan tuyến tính và đa biến.

Nội dung chính theo từng chương

Chương 1: Tổng quan về vấn đề nghiên cứu

Chương 1 hệ thống hóa toàn diện các tài liệu khoa học trong nước và quốc tế về loài Thằn lằn bóng đốm Eutropis macularius. Nội dung chương làm rõ lịch sử biến động danh pháp phân loại học, các dẫn liệu về hình thái học, phân bố địa lý từ Nam Á đến Đông Nam Á và các tỉnh thành ở Việt Nam. Đồng thời, chương 1 điểm luận chi tiết các cơ sở lý luận về sinh thái vi môi trường, mật độ quần thể, các khung mô hình chiếm cứ điểm hiện đại, sinh thái học dinh dưỡng, tập tính sinh sản và các công trình ứng dụng chỉ thị phân tử (NST, RAPD, giải trình tự 12S/16S rDNA) trong nghiên cứu di truyền thằn lằn bóng. Cuối chương, tác giả khái quát các đặc trưng tự nhiên của vùng Cao nguyên Buôn Ma Thuột – Buôn Hồ: địa hình cao nguyên lượn sóng (độ cao trung bình 400–800 m), nền đất đỏ bazan và đất xám bán bình nguyên Ea Súp; chế độ khí hậu nhiệt đới gió mùa với mùa mưa (tháng V–X) chiếm 84% lượng mưa năm và mùa khô (tháng XI–IV) gây bốc hơi mạnh (cao nhất vào tháng III đạt 230 mm).

Chương 2: Đối tượng, địa điểm và phương pháp nghiên cứu

Chương 2 mô tả chi tiết đối tượng nghiên cứu là Eutropis macularius, địa bàn triển khai chính tại VQG Yok Don, huyện Buôn Đôn, huyện Krông Búk, TP. Buôn Ma Thuột (Đắk Lắk) và các điểm thu mẫu bổ sung tại Đắk Nông, Gia Lai, Kon Tum. Chương này chuẩn hóa quy trình 8 nhóm phương pháp nghiên cứu thực nghiệm:

  1. Phương pháp khảo sát thực địa, đặt ô tiêu chuẩn và thu thập mẫu vật.
  2. Phương pháp đo đếm các chỉ tiêu hình thái và kiểm định dị hình giới tính SSD.
  3. Phương pháp ước tính mật độ quần thể trên ô tiêu chuẩn.
  4. Phương pháp đo đạc thông số vi khí hậu và xác định 6 dạng vi môi trường sống.
  5. Phương pháp xây dựng ma trận lịch sử phát hiện (0/1) để ước lượng tỉ suất chiếm cứ điểm ($\psi$) và xác suất phát hiện ($p$) trên phần mềm PRESENCE.
  6. Phương pháp phân tích sinh thái dinh dưỡng định lượng qua hệ thống chỉ số $%N$, $%V$, $%F$, IRI và chỉ số Simpson.
  7. Phương pháp phân tích sinh sản giải phẫu kết hợp làm tiêu bản mô học lát cắt tinh hoàn và buồng trứng theo các tháng trong năm.
  8. Quy trình tách chiết DNA, phản ứng PCR nhân đoạn gen 16S rDNA và các thuật toán phân tích di truyền học quần thể.

Chương 3: Kết quả nghiên cứu và thảo luận

1. Đặc điểm hình thái và sai khác hình thái theo giới tính

Số liệu đo đạc chỉ ra loài có kích thước trung bình. Chỉ số sai khác hình thái giới tính tính được là $\text{SSD} = 0,012$, cho thấy con cái có kích thước chiều dài thân (SVL) trung bình hơi nhỉnh hơn con đực. Tuy nhiên, các chỉ số về kích thước đầu như chiều dài đầu (HL), chiều rộng đầu (HW) và chiều rộng miệng (MW) ở cá thể đực lại lớn hơn cá thể cái. Phân tích tương quan cho thấy chiều dài thân (SVL) có mối liên hệ thuận chặt chẽ với khối lượng cơ thể (BM), chiều dài đuôi (TaL), chiều rộng đầu (HW) và chiều rộng miệng (MW) ở cả hai giới.

2. Mật độ quần thể và sử dụng vi môi trường sống

  • Mật độ quần thể: Khảo sát trên hệ thống ô tiêu chuẩn tại VQG Yok Don ghi nhận mật độ trung bình của loài là 14 cá thể/ha. Mật độ tại vùng lõi đạt 15 cá thể/ha, cao hơn vùng đệm (12 cá thể/ha); sự sai khác này có ý nghĩa thống kê ($P = 0,036$).
  • Sử dụng vi môi trường: Loài xuất hiện tại 6 dạng vi môi trường sống khác nhau, trong đó trảng cây bụi và thảm lá khô là 2 sinh cảnh chiếm ưu thế vượt trội:
Loại vi môi trường sống Tỉ lệ sử dụng (%) Nhiệt độ trung bình (°C) Độ ẩm trung bình (%)
Trảng cây bụi 44,53 28,47 ± 0,49 66,36 ± 2,48
Thảm lá khô 38,47 28,65 ± 0,46 66,04 ± 2,25
Gốc cây thân gỗ Có ghi nhận — —
Trên thân cây Có ghi nhận — —
Bụi tre Có ghi nhận — —
Môi trường khác Có ghi nhận — —

Phân tích hồi quy bội chứng minh nhiệt độ không khí và độ ẩm tại vi môi trường tác động có ý nghĩa thống kê đến tần suất sử dụng không gian sống của thằn lằn ($P = 0,037$).

3. Ước lượng xác suất phát hiện và tỉ suất chiếm cứ điểm

Khảo sát 72 ô tiêu chuẩn tại VQG Yok Don cho kết quả:

  • Tỉ suất chiếm cứ điểm thuần túy (naïve occupancy): Mùa mưa đạt 0,4722; mùa khô đạt 0,5417.
  • Chọn lọc mô hình AIC: Mô hình tổng quát nhiều thông số nhất $[\psi(\text{RK}), p(\text{ND}, \text{N}, \text{M}, \text{KXĐ})]$ (kết hợp yếu tố sinh cảnh rừng khộp với biến số nhiệt độ không khí, nắng, mưa và không xác định) là mô hình tối ưu. Khi đưa các biến môi trường vào mô hình, xác suất chiếm cứ điểm ước tính trong mùa mưa là 0,4723 và mùa khô là 0,6054, đều cao hơn tỉ suất thuần túy. Sinh cảnh rừng khộp (RK) có ảnh hưởng tích cực rõ rệt đến xác suất phát hiện loài so với sinh cảnh rừng trồng (RT).

4. Đặc điểm sinh thái học dinh dưỡng

Dạ dày thằn lằn chứa 17 loại thức ăn. Trong đó, 7 nhóm con mồi quan trọng nhất chiếm tổng chỉ số $\text{IRI} = 77,43%$ gồm:

  1. Ấu trùng côn trùng (Insect Larvae)
  2. Bộ Cánh màng (Hymenoptera)
  3. Bộ Cánh thẳng (Orthoptera)
  4. Mối (Isoptera / Blattodea)
  5. Thực vật (Plant)
  6. Bộ Cánh cứng (Coleoptera) / Bộ Chuồn chuồn (Odonata)
  7. Bộ Nhện (Araneae)

Kích thước con mồi tiêu thụ có mối tương quan thuận rõ nét với chiều dài thân (SVL) và độ rộng miệng (MW) của thằn lằn ở cả hai giới.

5. Đặc điểm sinh học sinh sản

  • Cá thể đực: Sinh sản có tính chu kỳ mùa vụ rõ rệt. Quá trình phát triển tuyến sinh dục bắt đầu từ tháng III và kết thúc vào cuối tháng VII. Thể tích tinh hoàn tăng mạnh vào các tháng giữa năm, tương ứng với sự gia tăng kích thước gan và hoạt động sinh tinh bên trong các ống sinh tinh (quan sát rõ cấu trúc phân chia tế bào sinh dục trên lát cắt mô học tháng VI đến tháng VIII).
  • Cá thể cái: Trứng giai đoạn 1 (noãn non) xuất hiện chủ yếu trong dạ con vào các tháng I–II và VIII–XII; giai đoạn từ tháng III đến tháng VII không thấy trứng giai đoạn 1. Trứng giai đoạn 2 bắt đầu tích lũy noãn hoàng từ cuối tháng I đến tháng IV. Trứng giai đoạn 3 (thể tích cực đại, sẵn sàng đẻ) xuất hiện tập trung vào giữa mùa hè. Chu kỳ sinh dưỡng và tích lũy noãn hoàn toàn đồng pha với chu kỳ phát triển tinh hoàn của cá thể đực.

6. Đặc điểm di truyền phân tử

Phân tích đoạn gen ty thể 16S rDNA trên 16 mẫu Thằn lằn bóng đốm thu từ 4 tỉnh Tây Nguyên ghi nhận 8 haplotype khác nhau:

Tỉnh thu mẫu Số lượng mẫu Đa dạng Haplotype ($Hd$) Đa dạng Nucleotide ($\pi$)
Kon Tum 4 0,833 0,00277
Gia Lai 4 0,833 0,00277
Đắk Lắk 4 0,500 0,00277
Đắk Nông 4 0,500 0,02415

Mức độ khác biệt di truyền giữa các quần thể dao động trong khoảng 0,14% – 2,66%. Các chỉ số kiểm định tính trung lập chứng minh quần thể Thằn lằn bóng đốm Tây Nguyên tiến hóa theo hướng chọn lọc ngẫu nhiên, trung tính, quần thể mở rộng sau thời kỳ bị ngăn cách địa lý, với sự hiện diện của nhiều alen hiếm ở tần số cao.

Kết quả và những đóng góp mới

Những đóng góp mới về mặt lý luận và khoa học

  1. Hình thái học: Cung cấp bộ dữ liệu định lượng chuẩn xác về 16 chỉ số kích thước cơ thể và xác lập chỉ số sai khác kích thước giới tính ($\text{SSD} = 0,012$) của loài Eutropis macularius tại Việt Nam; xác định mối quan hệ tương quan giữa kích thước đầu, miệng và thân với cơ chế bắt mồi.
  2. Sinh thái học quần thể: Lần đầu tiên công bố số liệu thực nghiệm về mật độ quần thể (14 cá thể/ha), định lượng tỉ lệ sử dụng 6 dạng vi môi trường và chứng minh sinh cảnh rừng khộp cùng thảm lá khô, trảng bụi là ổ sinh thái cốt lõi của loài.
  3. Mô hình chiếm cứ điểm: Lần đầu tiên ứng dụng thành công khung phân tích mô hình chiếm cứ điểm nhiều thông số của MacKenzie tại Tây Nguyên, lượng hóa được tác động của các biến số thời tiết (nhiệt độ, mưa, nắng) và kiểu rừng đến xác suất bắt gặp loài.
  4. Sinh học sinh sản và dinh dưỡng: Làm sáng tỏ phổ thức ăn gồm 17 loại con mồi (7 nhóm chiếm $\text{IRI} = 77,43%$); chứng minh chu kỳ sinh sản theo mùa (tháng III đến tháng VII) bằng dẫn liệu giải phẫu học kết hợp cấu trúc vi thể mô học ống sinh tinh và noãn hoàng.
  5. Di truyền học phân tử: Thiết lập cơ sở dữ liệu trình tự 16S rDNA cho thằn lằn bóng tại Tây Nguyên; phát hiện 8 haplotype với mức độ đa dạng di truyền cao tại Kon Tum, Gia Lai và mức phân hóa nucleotide cao nhất tại Đắk Nông.

Kiến nghị và giải pháp đề xuất

  • Đối với nghiên cứu khoa học tiếp theo: Mở rộng nghiên cứu giải mã thêm các chỉ thị gen nhân và gen ty thể khác (như COI, Cyt-b, 12S rRNA) nhằm làm rõ hơn bức tranh phát sinh chủng loại; tiếp tục theo dõi biến động mật độ quần thể trong mối tương quan lâu dài với biến đổi khí hậu toàn cầu.
  • Đối với công tác quản lý và bảo tồn:
    • Bảo vệ nguyên vẹn cấu trúc thảm thực vật tự nhiên của hệ sinh thái rừng khộp tại VQG Yok Don và các vùng phụ cận, hạn chế tối đa việc đốt rừng, dọn sạch thảm lá khô và phá bỏ trảng bụi trong mùa khô.
    • Quản lý chặt chẽ hoạt động sử dụng hóa chất bảo vệ thực vật tại các vườn cây công nghiệp giáp ranh rừng tự nhiên nhằm bảo tồn nguồn côn trùng làm thức ăn cho thằn lằn.
    • Kiểm soát hoạt động săn bắt trái phép thằn lằn bóng phục vụ mục đích thương mại và thực phẩm tại các địa phương thuộc vùng Tây Nguyên.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế của luận án:
    • Dữ liệu khảo sát chi tiết về mật độ và mô hình chiếm cứ điểm tập trung chủ yếu tại sinh cảnh VQG Yok Don, chưa thể phủ kín toàn bộ các tiểu vùng sinh thái khác trên toàn Cao nguyên Buôn Ma Thuột – Buôn Hồ.
    • Phân tích di truyền mới thực hiện trên một đoạn chỉ thị phân tử ty thể đơn lẻ (16S rDNA) với dung lượng mẫu gồm 16 cá thể đại diện cho 4 tỉnh.
  • Hướng nghiên cứu tiếp theo:
    • Mở rộng phạm vi đánh giá chiếm cứ điểm trên các sinh cảnh nhân tạo (vườn cà phê, cao su, điều, bơ) ở các độ cao khác nhau.
    • Tăng cường cỡ mẫu và kết hợp đồng thời các chỉ thị gen ty thể và gen nhân để đánh giá sâu hơn dòng gen (gene flow) giữa các quần thể bị chia cắt bởi hoạt động phát triển hạ tầng giao thông và nông nghiệp tại Tây Nguyên.

Giá trị tham khảo

  • Đối với nghiên cứu sinh và học viên sau đại học ngành Sinh thái học, Động vật học: Cung cấp tài liệu mẫu mực về phương pháp luận nghiên cứu kết hợp giữa khảo sát sinh thái thực địa định lượng (mô hình chiếm cứ điểm PRESENCE, chỉ số IRI, giải phẫu mô học) với công nghệ sinh học phân tử hiện đại.
  • Đối với giảng viên và các viện nghiên cứu: Bổ sung tư liệu giảng dạy chuyên sâu về khu hệ bò sát Tây Nguyên, sinh thái học hành vi và di truyền học bảo tồn động vật có xương sống.
  • Đối với các nhà quản lý môi trường và ban quản lý vườn quốc gia, khu bảo tồn: Cung cấp các thông số định lượng cụ thể về sinh cảnh sống ưu thích (rừng khộp, thảm lá khô, nhiệt độ 28,47–28,65°C) để xây dựng kế hoạch phân khu bảo tồn và quy trình giám sát đa dạng sinh học định kỳ tại VQG Yok Don và các vùng lân cận.

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

1. Tỉ suất chiếm cứ điểm và xác suất phát hiện loài Thằn lằn bóng đốm tại VQG Yok Don biến động như thế nào giữa hai mùa?
Tỉ suất chiếm cứ điểm thuần túy của loài vào mùa mưa là 0,4722 và mùa khô là 0,5417. Khi áp dụng mô hình tổng quát kết hợp sinh cảnh rừng khộp và các yếu tố thời tiết $[\psi(\text{RK}), p(\text{ND}, \text{N}, \text{M}, \text{KXĐ})]$, xác suất chiếm cứ điểm ước tính tăng lên đạt 0,4723 vào mùa mưa và 0,6054 vào mùa khô, khẳng định điều kiện vi khí hậu và sinh cảnh rừng khộp ảnh hưởng lớn đến khả năng phát hiện loài.

2. Thằn lằn bóng đốm tại vùng nghiên cứu sử dụng chủ yếu những loại vi môi trường nào?
Loài sử dụng 6 dạng vi môi trường sống, nhưng tập trung nhiều nhất ở 2 dạng: trảng cây bụi (chiếm 44,53% với nhiệt độ $28,47 \pm 0,49^\circ\text{C}$, độ ẩm $66,36 \pm 2,48%$) và thảm lá khô (chiếm 38,47% với nhiệt độ $28,65 \pm 0,46^\circ\text{C}$, độ ẩm $66,04 \pm 2,25%$).

3. Thành phần thức ăn quan trọng nhất của loài gồm những nhóm sinh vật nào?
Trong tổng số 17 loại thức ăn được ghi nhận trong dạ dày, 7 nhóm con mồi quan trọng nhất chiếm tới $77,43%$ tổng chỉ số IRI gồm: Ấu trùng côn trùng, bộ Cánh màng (Hymenoptera), bộ Cánh thẳng (Orthoptera), mối (Isoptera/Blattodea), thực vật, bộ Cánh cứng (Coleoptera)/Chuồn chuồn (Odonata) và bộ Nhện (Araneae).

4. Chu kỳ sinh sản của Thằn lằn bóng đốm diễn ra vào thời gian nào trong năm?
Mùa sinh sản của loài bắt đầu từ khoảng tháng III và kết thúc vào cuối tháng VII hàng năm. Trong giai đoạn này, tinh hoàn con đực đạt thể tích tối đa và tích cực sinh tinh, đồng thời buồng trứng con cái phát triển noãn qua các giai đoạn 2 và 3 để chuẩn bị cho việc đẻ trứng.

5. Mức độ đa dạng di truyền dựa trên đoạn gen 16S rDNA giữa các tỉnh Tây Nguyên có đặc điểm gì?
Phân tích 16 mẫu ghi nhận 8 haplotype. Đa dạng haplotype ($Hd$) đạt mức cao ở Kon Tum và Gia Lai ($0,833$) nhưng thấp hơn ở Đắk Lắk và Đắk Nông ($0,500$). Tuy nhiên, độ đa dạng nucleotide ($\pi$) cao nhất lại ghi nhận tại quần thể Đắk Nông ($0,02415$).

Kết luận

Luận án của nghiên cứu sinh Trương Bá Phong đã hệ thống hóa và cung cấp hệ thống dữ liệu thực nghiệm toàn diện về đặc điểm hình thái, sai khác giới tính, sinh thái dinh dưỡng, sinh thái sinh sản, mô hình chiếm cứ điểm và cấu trúc di truyền phân tử của loài Thằn lằn bóng đốm Eutropis macularius tại vùng Cao nguyên Buôn Ma Thuột – Buôn Hồ và Tây Nguyên. Kết quả nghiên cứu khẳng định mối liên hệ chặt chẽ giữa sinh cảnh rừng khộp, thảm lá khô với sự tồn tại của loài, đồng thời chỉ ra đặc điểm thích nghi dinh dưỡng và chu kỳ sinh sản theo mùa vụ rõ rệt. Những phát hiện về 8 haplotype gen ty thể 16S rDNA cùng các kiểm định trung lập cung cấp dẫn liệu khoa học giá trị phục vụ công tác giám sát đa dạng sinh học và đề xuất các giải pháp quản lý, bảo tồn bền vững nguồn lợi bò sát tại khu vực Tây Nguyên.