Tổng quan về luận án

Nghiên cứu sinh thái học và di truyền quần thể bò sát tại các hệ sinh thái rừng khộp và cao nguyên nhiệt đới gió mùa giữ vai trò then chốt trong đánh giá sức khỏe sinh thái trước áp lực nhân sinh và biến đổi khí hậu toàn cầu (Cox et al., 2022). Luận án tiến sĩ sinh học chuyên ngành Sinh thái học (mã số: 9420120) của 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à PGS. Ngô Văn Bình tại Trường Đại học Đà Lạt (2023), mang tiêu đề: "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ồ". Công trình tiên phong thiết lập cơ sở dữ liệu tích hợp đa tầng từ hình thái, vi sinh cảnh, mô hình chiếm cứ điểm, dinh dưỡng, sinh sản mô học cho đến cấu trúc phả hệ phân tử của loài Thằn lằn bóng đốm Eutropis macularius tại khu vực Tây Nguyên.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định rõ: phần lớn các công trình trước đây chỉ dừng lại ở mức điều tra danh lục thành phần loài đơn thuần (Nguyễn Văn Sáng & Hồ Thu Cúc, 1996; Nguyen et al., 2009) hoặc khảo sát dinh dưỡng cục bộ với dung lượng mẫu hạn chế (Hoàng Xuân Quang et al., 2009; Ngô Đắc Chứng et al., 2020), hoàn toàn khuyết thiếu các mô hình toán học đánh giá xác suất chiếm cứ điểm kết hợp đa biến sinh cảnh, cũng như chưa từng có phân tích di truyền quần thể bằng chỉ thị phân tử 16S rDNA cho loài này tại vùng chuyển tiếp sinh thái Buôn Ma Thuột – Buôn Hồ.

Luận án giải quyết 4 câu hỏi nghiên cứu và giả thuyết tương ứng:

  1. Câu hỏi 1: Mức độ dị hình giới tính kích thước (SSD) biểu hiện như thế nào và các biến số hình thái tương quan ra sao? Giả thuyết H1: Tồn tại sự phân hóa kích thước đầu và chiều dài thân giữa hai giới nhằm giảm thiểu cạnh tranh nội loài.
  2. Câu hỏi 2: Mật độ cá thể và mô hình chiếm cứ điểm ($\psi$) bị chi phối bởi những yếu tố sinh cảnh và thời tiết nào? Giả thuyết H2: Sinh cảnh rừng khộp cùng tổ hợp nhiệt - ẩm tối ưu hóa xác suất phát hiện ($p$) và nâng cao tỉ suất chiếm cứ thực tế so với tỉ suất chiếm cứ thuần túy (naïve occupancy).
  3. Câu hỏi 3: Chiến lược sử dụng tài nguyên dinh dưỡng và chu kỳ sinh sản vận hành đồng bộ ra sao theo mùa? Giả thuyết H3: Loài thích nghi với chế độ ăn phổ rộng, kích thước con mồi tỷ lệ thuận với độ rộng miệng, và chu kỳ phát triển giao tử đực - cái đạt đỉnh đồng pha trong mùa mưa ẩm.
  4. Câu hỏi 4: Cấu trúc di truyền quần thể tại 4 tỉnh Tây Nguyên có phân hóa sâu sắc không? Giả thuyết H4: Các quần thể chịu sự chi phối của tiến hóa trung tính ngẫu nhiên và phản ánh hiện tượng mở rộng quần thể sau thắt cổ chai địa lý.

Khung lý thuyết của luận án được xây dựng vững chắc trên nền tảng: Thuyết Dị hình kích thước giới tính (Sexual Size Dimorphism Theory - Fairbairn, 1997; Cox et al., 2003), Lý thuyết Sử dụng vi sinh cảnh và Ổ sinh thái (Niche Partitioning Theory - Pianka, 1973), Khung mô hình Toán học Chiếm cứ điểm (Site Occupancy Modeling - MacKenzie et al., 2002, 2006), Thuyết Tối ưu hóa Dinh dưỡng (Optimal Foraging Theory - MacArthur & Pianka, 1966), và Thuyết Tiến hóa Trung tính Phân tử (Neutral Theory of Molecular Evolution - Kimura, 1983; Tajima, 1989).

Phạm vi không gian nghiên cứu bao quát hệ thống 72 ô tiêu chuẩn tại Vườn quốc gia (VQG) Yok Don (vùng lõi và vùng đệm thuộc Đắk Lắk) cùng các điểm thu mẫu mở rộng tại Đắk Nông, Gia Lai và Kon Tum trong khung thời gian thực địa liên tục từ năm 2017 đến năm 2021. Đóng góp đột phá của luận án mang tính định lượng rõ rệt: xác định mật độ quần thể đạt 14 cá thể/ha, định lượng 6 kiểu vi môi trường với $83,00%$ tập trung ở trảng bụi và thảm lá khô, phân lập 17 bộ thức ăn với chỉ số quan trọng (IRI) đạt $77,43%$ ở 7 nhóm mồi chính, và giải trình tự thành công 8 haplotype từ 16 mẫu đại diện phục vụ quản trị bảo tồn sinh thái.

Literature Review và Positioning

Các công trình nghiên cứu trên thế giới về giống Eutropis (trước đây thuộc giống Mabuya) đã hình thành những dòng nghiên cứu lớn. Về mặt phân loại học và danh pháp, Mausfeld & Schmitz (2003) cùng Murthy (2010) đã tái cấu trúc chi và xác lập tên khoa học chuẩn hóa là Eutropis macularius (Blyth, 1853), thay thế cho các danh pháp đồng vật lịch sử như Euprepis macularius hay Mabuya macularia (Boulenger, 1887; Smith, 1935). Về sinh thái học dinh dưỡng và hành vi sinh sản, các nghiên cứu kinh điển của Huang (2006a, 2006b, 2011) trên đảo Lan-yu (Đài Loan) đã chứng minh mối liên kết chặt chẽ giữa tính sẵn có của con mồi với kích thước khoang miệng, cũng như sự đánh đổi sinh học trong tập tính bảo vệ tổ trứng trước áp lực sinh vật săn mồi (Shine, 1980; Shine & Harlow, 1996).

Tuy nhiên, tồn tại những tranh luận học thuật sâu sắc:

  • Tranh luận 1: Động lực chi phối dị hình giới tính (SSD). Nhóm quan điểm chọn lọc sinh sản (Olsson, 2002; Warner & Shine, 2005, 2007) cho rằng con cái phát triển chiều dài thân (SVL) lớn hơn để tối đa hóa không gian mang trứng (fecundity selection), trong khi nhóm quan điểm chọn lọc giới tính đối kháng (Anderson, 1994; Brana, 1996) khẳng định con đực phải phát triển đầu to hơn (HW, HL) phục vụ giao tranh lãnh thổ.
  • Tranh luận 2: Phương pháp ước tính phân bố và mật độ. Phương pháp đánh dấu bắt lại cổ điển (Pollock, 1982) bộc lộ nhược điểm làm xáo trộn sinh vật, đối lập với trường phái mô hình hóa xác suất phát hiện không hoàn hảo dựa trên lịch sử quan sát lặp lại (MacKenzie et al., 2002, 2006), vốn chưa được ứng dụng hệ thống trên bò sát tại Việt Nam ngoài các thử nghiệm sơ khởi của Cao Thị Thanh Nguyên et al. (2018, 2019) trên loài nhông cát.

Luận án định vị chính xác ở điểm giao thoa giữa sinh thái học thực địa và di truyền học bảo tồn. Khi so sánh với nghiên cứu của Bullock & Evans (1990) tại Dominica (mật độ bò sát đạt ngưỡng cực đại 751-2148 cá thể/ha ở rừng ven biển) và công trình của Vitt et al. (1997) tại Amazon về loài Mabuya nigropunctata, luận án đã chứng minh mật độ 14 cá thể/ha ở rừng khộp Tây Nguyên chịu sự điều tiết khắt khe của tính phân mùa nhiệt đới gió mùa. Đồng thời, so với các dữ liệu di truyền học tế bào của Ota et al. (1996, 2001) tại Thái Lan ghi nhận biến dị nhiễm sắc thể giữa các vùng Mae Yom và Phu Wua, luận án của Trương Bá Phong đã tiến xa hơn khi cung cấp dữ liệu trình tự nucleotide gen ty thể 16S rDNA, lấp đầy lỗ hổng cấu trúc di truyền phân tử của loài ở hạ lưu lưu vực sông Mê Kông và cao nguyên Tây Nguyên.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và kiểm chứng Thuyết Dị hình Kích thước Giới tính của Fairbairn (1997) và công thức định lượng SSD của Cox, Skelly & Alder (2003). Kết quả xác định chỉ số $\text{SSD} = 0,012$, chứng minh mức độ sai khác hình thái giới tính ở E. macularius ở mức độ tinh vi, trong đó kích thước chiều dài thân (SVL) có tương quan tương hỗ với khối lượng thân (BM), chiều dài đuôi (TaL) và kích thước đầu (HW, MW), phù hợp với mô hình tối ưu hóa vận động cơ học trong vi sinh cảnh thảm lá rụng.

Về lý thuyết sinh thái học quần thể, nghiên cứu đóng góp mô hình lý thuyết Occupancy-Detection Probability: $$\psi(\text{Habitat}) \times p(\text{Weather, Survey})$$ chứng minh rằng xác suất quan sát một cá thể bò sát nhỏ không chỉ phản ánh sự hiện diện sinh học thuần túy mà là một hàm xác suất bị chi phối bởi cấu trúc tán che rừng khộp và các biến thời tiết vi khí hậu tức thời.

Luận án củng cố Thuyết Tiến hóa Trung tính của Kimura (1983) và mô hình kiểm định trung tính Tajima's $D$, chứng minh rằng sự phân hóa di truyền giữa các quần thể Tây Nguyên ($F_{ST}$ từ $0,14%$ đến $2,66%$) vận hành chủ yếu theo cơ chế phiêu bạt di truyền (genetic drift) và mở rộng quần thể cục bộ sau chia cắt địa hình tự nhiên, thay vì áp lực chọn lọc định hướng mạnh.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên ngành 3 trục lý thuyết:

  1. Trục Hình thái - Chức năng: Khám phá tương quan giữa kích thước cấu trúc cơ thể (SVL, MW) với kích thước và thể tích con mồi tiêu thụ (Trophic Niche).
  2. Trục Mô hình Toán sinh thái Không gian: Ứng dụng mô hình chiếm cứ điểm đơn mùa (Single-Season Occupancy Model) của MacKenzie et al. (2006) kết hợp tiêu chuẩn chọn lọc mô hình thông tin Akaike (AIC - Akaike Information Criterion) để phân tích tương tác giữa kiểu rừng và vi khí hậu.
  3. Trục Phả hệ - Di truyền Quần thể: Tích hợp phân tích haplotype ty thể với các chỉ số trung lập (Tajima’s $D$, Fu & Li’s $D^, F^$).

Điều kiện biên (boundary conditions) được xác định cụ thể: các mô hình sinh thái áp dụng tối ưu cho sinh cảnh rừng rụng lá theo mùa (rừng khộp), trảng cây bụi nhiệt đới gió mùa tại vùng cao nguyên có độ cao từ 300m đến 800m so với mực nước biển, nơi có sự tương phản sâu sắc giữa mùa khô (tháng XI - IV) và mùa mưa (tháng V - X).

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

Thiết kế nghiên cứu

Luận án vận hành theo hệ hình thực chứng luận (Positivism) với quan điểm nhận thức luận khách quan, kết hợp thiết kế nghiên cứu hỗn hợp định lượng đa cấp độ (Multi-level Quantitative Design):

  • Cấp độ cá thể: Hình thái học, giải phẫu sinh dưỡng, cấu trúc mô học cơ quan sinh sản.
  • Cấp độ quần thể - không gian: Mật độ cá thể, sử dụng vi sinh cảnh, mô hình chiếm cứ 72 ô tiêu chuẩn.
  • Cấp độ phân tử: Phân tích đa dạng di truyền nucleotide và phả hệ haplotype từ 4 tỉnh.

Mẫu nghiên cứu hình thái và sinh thái được thu thập thực địa liên tục, đảm bảo tính đại diện thống kê trên toàn bộ diện tích rừng lõi và rừng đệm VQG Yok Don, không để xảy ra sai số chọn mẫu cục bộ.

Quy trình nghiên cứu rigorous

Chiến lược chọn mẫu và quy trình thu thập dữ liệu tuân thủ chuẩn mực quốc tế:

  • Khảo sát hình thái: Sử dụng thước kẹp điện tử Mitutoyo (độ chính xác 0,01 mm) và cân điện tử (độ chính xác 0,01 g) đo đạc các chỉ số: SVL, TaL, HL, HW, MW, BM theo quy chuẩn của Cox et al. (2003).
  • Mô hình Occupancy: Thiết lập 72 ô tiêu chuẩn (diện tích $100,\text{m}^2$/ô) với khoảng cách tối thiểu giữa các ô là 200m để đảm bảo tính độc lập thống kê. Mỗi ô được quan sát lặp lại 3-4 lần trong cả hai mùa. Ghi nhận ma trận phát hiện (1) / không phát hiện (0) đồng thời đo đạc tức thời nhiệt độ không khí, độ ẩm vi môi trường bằng máy đo chuyên dụng.
  • Nghiên cứu dinh dưỡng: Phân tích dạ dày theo phương pháp súc rửa dạ dày không sát thương kết hợp giải phẫu mẫu tiêu bản; định loại con mồi đến cấp bộ/họ dưới kính hiển vi soi nổi; đo thể tích con mồi ($V$) theo công thức hình elip: $$V = \frac{4}{3} \pi \left(\frac{L}{2}\right) \left(\frac{W}{2}\right)^2$$ tính toán Chỉ số Quan trọng Tương đối (Index of Relative Importance - IRI): $$\text{IRI} = (%N + %V) \times %F$$
  • Nghiên cứu sinh thái sinh sản và mô học: Mổ khám giải phẫu xác định kích thước buồng trứng, thể tích tinh hoàn ($V_{\text{testis}}$), đo thể tích gan ($V_{\text{liver}}$) đánh giá dự trữ năng lượng. Cắt lát mô học cơ quan sinh sản cố định trong dung dịch Bouin, đúc parafin, nhuộm Hematoxylin - Eosin (HE) và quan sát cấu trúc ống sinh tinh/noãn bào dưới kính hiển vi quang học Olympus.
  • Nghiên cứu di truyền phân tử: Tách chiết DNA tổng số bằng bộ kit thương mại chuyên dụng; nhân bản đoạn gen ty thể 16S rDNA bằng kỹ thuật PCR với cặp mồi đặc hiệu; giải trình tự nucleotide trực tiếp trên hệ thống tự động ABI 3500; đối chiếu trình tự chuẩn quốc tế mang mã số GenBank AB057394 (Eutropis macularia).

Data và phân tích

Dữ liệu được xử lý bằng hệ thống phần mềm chuyên ngành tiêu chuẩn cao:

  • Phần mềm PRESENCE (version 12.x) và gói lệnh tương đương để ước tính xác suất chiếm cứ ($\psi$) và xác suất phát hiện ($p$), phân tích đối sánh tập hợp mô hình ứng viên thông qua giá trị trọng số Akaike ($w_i$).
  • Phần mềm SPSS 22.0 và R để thực hiện kiểm định t-test độc lập, Phân tích phương sai (One-way ANOVA), hồi quy tuyến tính đa biến, phân tích hệ số tương quan Pearson và tính chỉ số đa dạng Simpson nghịch đảo ($1/D$).
  • Phần mềm BioEdit, MEGA 11 và DnaSP 6.0 để căn chỉnh chuỗi (alignment), tính toán độ đa dạng haplotype ($H_d$), độ đa dạng nucleotide ($\pi$), khoảng cách di truyền Kimura 2-parameter ($F_{ST}$), và kiểm định trung tính tiến hóa Tajima's $D$, Fu & Li's $D^$ và $F^$.

Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án mang lại 5 phát hiện cốt lõi với minh chứng số liệu định lượng chuẩn xác:

"Mật độ quần thể của Thằn lằn bóng đốm trên các ô tiêu chuẩn ở VQG Yok Don là 14 cá thể/ha. Trong đó, mật độ quần thể ở vùng lõi là 15 cá thể/ha và ở vùng đệm là 12 cá thể/ha. Sự sai khác mật độ quần thể ở vùng đệm và vùng lõi có ý nghĩa thống kê ($P = 0,036$)."

"Thằn lằn bóng đốm đã sử dụng 6 loại vi môi trường sống ở khu vực nghiên cứu... Trong đó, 2 vi môi trường sống chủ yếu là trảng cây bụi và thảm lá khô với tỉ lệ lần lượt chiếm $44,53%$ và $38,47%$ với nhiệt độ và độ ẩm trung bình tại vi môi trường là $28,47 \pm 0,49^\circ\text{C}$; $66,36 \pm 2,48%$ và $28,65 \pm 0,46^\circ\text{C}$; $66,04 \pm 2,25%$. Nhiệt độ không khí và độ ẩm đều ảnh hưởng có ý nghĩa đến việc sử dụng vi môi trường sống của loài ($P = 0,037$)."

"Mô hình nhiều thông số nhất $[\psi(\text{RK}), p(\text{ND, N, M, KXĐ})]$ là môi trường sống rừng khộp và có sự kết hợp với các yếu tố như nhiệt độ không khí và tình hình nắng mưa thì xác suất chiếm cứ điểm của loài Thằn lằn bóng đốm là 0,4723 vào mùa mưa và 0,6054 vào mùa khô đều cao hơn tỉ suất chiếm cứ điểm thuần túy."

  1. Phát hiện về Dị hình giới tính và Hình thái: Chỉ số dị hình giới tính đạt $\text{SSD} = 0,012$. Không có sự phân hóa cực đoan về kích thước cơ thể chung giữa hai giới, tuy nhiên các chỉ số đầu (HW, MW) có xu hướng mở rộng ở cá thể đực nhằm phục vụ cạnh tranh giao phối, trong khi tương quan SVL với khối lượng cơ thể (BM) duy trì chặt chẽ ở cả hai giới ($P < 0,01$).
  2. Phát hiện về Vi sinh cảnh và Mô hình Chiếm cứ: Tỉ suất chiếm cứ thuần túy (naïve occupancy) ở mùa mưa là 0,4722 và mùa khô là 0,5417. Khi mô hình hóa bằng PRESENCE, mô hình tổng quát tích hợp sinh cảnh rừng khộp (RK) và các biến thời tiết cho thấy xác suất phát hiện đạt đỉnh vào mùa khô ($p = 0,6054$), khẳng định rừng khộp là sinh cảnh cốt lõi vượt trội so với rừng trồng (RT).
  3. Phát hiện về Sinh thái Dinh dưỡng: Xác định 17 bộ thức ăn trong phổ dinh dưỡng. Nhóm 7 con mồi chủ lực chiếm tổng $\text{IRI} = 77,43%$ bao gồm: ấu trùng côn trùng, Hymenoptera (Cánh màng), Orthoptera (Cánh thẳng), Isoptera (Mối), mô thực vật, Coleoptera (Cánh cứng), và Araneae (Nhện). Kích thước chiều dài thân (SVL) và độ rộng miệng (MW) tương quan thuận có ý nghĩa thống kê với kích thước và thể tích con mồi tiêu thụ ($P < 0,05$).
  4. Phát hiện về Chu kỳ Sinh học Sinh sản và Mô học: Sinh sản mang tính phân mùa nghiêm ngặt từ tháng III đến tháng VII hàng năm. Cấu trúc mô học tinh hoàn chứng minh các ống sinh tinh nở rộng tối đa, chứa đầy tinh trùng hoạt động từ tháng III đến tháng VII và thoái hóa chuyển sang pha nghỉ từ tháng X đến tháng XII. Ở con cái, trứng giai đoạn 1 tích lũy từ tháng I-II và VIII-XII; trứng giai đoạn 2 và 3 (thành thục) phát triển đồng pha từ cuối tháng I đến tháng VII, đồng thời thể tích gan ($V_{\text{liver}}$) biến thiên tỷ lệ thuận biểu thị quá trình tổng hợp noãn hoàng (vitellogenesis).
  5. Phát hiện về Đa dạng Di truyền và Cấu trúc Phân tử: Phân tích 16 mẫu đại diện từ 4 tỉnh Tây Nguyên phát hiện 8 haplotype đặc trưng trên đoạn gen ty thể 16S rDNA. Độ đa dạng haplotype ($H_d$) đạt giá trị cao tại Kon Tum và Gia Lai ($H_d = 0,833$), nhưng thấp hơn tại Đắk Lắk và Đắk Nông ($H_d = 0,500$). Độ đa dạng nucleotide ($\pi$) cao nhất tại Đắk Nông ($\pi = 0,02415$), trong khi 3 tỉnh còn lại đạt trung bình $\pi = 0,00277$. Khoảng cách di truyền ($F_{ST}$) biến thiên từ $0,14%$ đến $2,66%$. Kiểm định trung tính chứng minh các quần thể đang trong pha giãn nở sau tắc nghẽn di truyền với sự ưu thế của các alen hiếm.

Implications đa chiều

  • Đóng góp lý thuyết: Cung cấp minh chứng thực nghiệm cho Thuyết Tiến hóa Trung tính và Thuyết Dị hình Giới tính tại vùng sinh thái nhiệt đới Đông Dương; chuẩn hóa mô hình Occupancy-AIC trong nghiên cứu bò sát cạn.
  • Đổi mới phương pháp luận: Thiết lập quy trình chuẩn tích hợp giữa mô hình hóa xác suất phát hiện không hoàn hảo với phân tích mô học sinh sản và di truyền nucleotide phân tử, có thể chuyển giao cho nghiên cứu các nhóm bò sát khác.
  • Ứng dụng thực tiễn: Khẳng định vai trò của E. macularius như một tác nhân kiểm soát sinh học tự nhiên đối với các loài côn trùng gây hại nông lâm nghiệp (chiếm trên $70%$ khẩu phần ăn); định hướng bảo tồn trảng bụi và lớp thảm mục lá khô trong canh tác vườn cây công nghiệp (cà phê, cao su, điều).
  • Khuyến nghị chính sách: Cung cấp luận cứ khoa học để VQG Yok Don và Chi cục Kiểm lâm Đắk Lắk thiết lập ranh giới bảo tồn nghiêm ngặt sinh cảnh rừng khộp vùng lõi; kiểm soát việc đốt thực bì vào mùa khô nhằm bảo vệ ổ trứng và vi sinh cảnh lá khô.
  • Khả năng tổng quát hóa: Mô hình sinh thái dinh dưỡng và chiếm cứ điểm có khả năng áp dụng cho toàn bộ dải rừng khộp rụng lá theo mùa thuộc vùng sinh thái hạ lưu sông Mê Kông (Việt Nam, Lào, Campuchia, Thái Lan).

Limitations và Future Research

Luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  1. Giới hạn về chỉ thị phân tử: Nghiên cứu mới dừng lại ở việc phân tích một đoạn gen ty thể 16S rDNA; chưa khai thác các chỉ thị nhân (nuclear genes) như RAG1, C-mos hoặc giải trình tự toàn bộ hệ gen ty thể (mitogenome) để đánh giá sâu hơn dòng gen (gene flow).
  2. Giới hạn cỡ mẫu di truyền: Số lượng 16 trình tự phân tích trên 4 tỉnh Tây Nguyên dù đủ phản ánh cấu trúc haplotype cơ bản nhưng cần mở rộng cỡ mẫu lên hàng trăm cá thể để phân tích mạng lưới quần thể phức hợp.
  3. Giới hạn không gian mô hình chiếm cứ: Dữ liệu chiếm cứ điểm chi tiết mới tập trung tại VQG Yok Don, chưa thực hiện mô hình đa mùa (multi-season dynamic occupancy) kéo dài trên 10 năm để đánh giá tỷ lệ tuyệt chủng và tái chiếm cứ cục bộ trước áp lực biến đổi khí hậu.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Ứng dụng công nghệ giải trình tự thế hệ mới (NGS/RAD-seq) để lập bản đồ cấu trúc di truyền cảnh quan (Landscape Genetics) trên toàn bộ bán đảo Đông Dương.
  • Nghiên cứu sinh lý học nhiệt thích ứng (Thermal Biology & Operative Temperature Models) để dự báo ngưỡng chống chịu nhiệt độ của loài trước các kịch bản ấm lên toàn cầu.
  • Đánh giá độc học sinh thái (Ecotoxicology) về sự tích lũy hóa chất bảo vệ thực vật trong mô mỡ và gan thằn lằn bóng tại các vùng canh tác nông nghiệp chuyên canh ở Tây Nguyên.

Tác động và ảnh hưởng

  • Tác động học thuật: Công trình mở ra chuẩn mực nghiên cứu sinh thái học bò sát hiện đại tại Việt Nam, dự kiến thu hút trích dẫn khoa học cao trong các mạng lưới nghiên cứu bò sát học nhiệt đới (Herpetology) và sinh học bảo tồn quốc tế.
  • Tác động ngành và sản xuất: Cung cấp cơ sở khoa học để các nông trường cao su, cà phê, điều tại Đắk Lắk chuyển đổi mô hình nông nghiệp sinh thái bền vững, hạn chế phun thuốc trừ sâu thông qua bảo tồn thiên địch bò sát nhỏ.
  • Ảnh hưởng chính sách: Cung cấp dữ liệu định lượng phục vụ trực tiếp cho Ban Quản lý VQG Yok Don, Sở Nông nghiệp và Phát triển Nông thôn tỉnh Đắk Lắk, Cục Bảo tồn Thiên nhiên và Đa dạng Sinh học (Bộ Tài nguyên và Môi trường) trong việc lập quy hoạch bảo tồn hệ sinh thái rừng rụng lá Tây Nguyên.
  • Lợi ích xã hội và toàn cầu: Góp phần bảo tồn vốn gen bản địa, đóng góp dữ liệu vào Ngân hàng Gen Quốc tế (GenBank) và bổ sung cơ sở dữ liệu cho Sách Đỏ IUCN về các loài bò sát nhiệt đới.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học: Tiếp cận mô hình phương pháp luận tích hợp chuẩn mực giữa PRESENCE, mô học giải phẫu và di truyền phân tử DnaSP/MEGA.
  • Giảng viên, nhà nghiên cứu cao cấp: Kế thừa bộ cơ sở dữ liệu toàn diện về hình thái, sinh thái và trình tự 16S rDNA của giống Eutropis phục vụ các công trình so sánh phát sinh chủng loại quy mô khu vực.
  • Chuyên gia R&D Nông nghiệp Sinh thái: Ứng dụng hiểu biết về phổ dinh dưỡng ăn côn trùng gây hại của loài để xây dựng giải pháp quản lý dịch hại tổng hợp (IPM).
  • Nhà quản lý khu bảo tồn và hoạch định chính sách: Sở hữu các tham số định lượng chính xác (14 cá thể/ha, vi môi trường trảng bụi $44,53%$, mùa sinh sản tháng III-VII) để thiết kế kế hoạch tuần tra, phòng cháy rừng và quản lý vùng đệm tối ưu.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc tích hợp thành công Thuyết Chiếm cứ Điểm Không gian (MacKenzie et al., 2006) với Thuyết Phân hóa Ổ Sinh thái (Pianka, 1973) trên bò sát nhiệt đới gió mùa Việt Nam. Luận án đã mở rộng mô hình chiếm cứ điểm bằng việc chứng minh rằng việc ước tính sự hiện diện loài nếu bỏ qua xác suất phát hiện không hoàn hảo ($p < 1$) sẽ dẫn đến sai lệch lớn, đồng thời lượng hóa được vai trò quyết định của cấu trúc rừng khộp kết hợp với nhiệt - ẩm lên xác suất chiếm cứ thực tế ($\psi = 0,6054$ mùa khô).

2. Đổi mới phương pháp luận của luận án so với ít nhất 2 nghiên cứu đi trước? Trả lời: So với nghiên cứu của Hoàng Xuân Quang et al. (2009) tại Bạch Mã và Ngô Đắc Chứng et al. (2020) tại Thừa Thiên Huế (chỉ thu thập mẫu ngẫu nhiên, mô tả tần số thức ăn giản đơn và không có mô hình toán học), luận án của Trương Bá Phong đã tạo bước đột phá phương pháp luận: (1) Thiết lập lưới 72 ô tiêu chuẩn định vị GPS cố định để chạy mô hình AIC trên phần mềm PRESENCE, loại bỏ sai số quan sát; (2) Kết hợp đồng thời kỹ thuật cắt lát mô học HE kiểm tra chu kỳ ống sinh tinh/noãn bào với giải trình tự gen ty thể 16S rDNA, thay vì chỉ giải phẫu mắt thường hoặc đo vẽ kích thước bên ngoài.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Trả lời: Dữ liệu cho thấy xác suất chiếm cứ điểm và xác suất phát hiện loài vào mùa khô ($\psi = 0,6054$) lại cao hơn đáng kể so với mùa mưa ($\psi = 0,4723$), trái ngược với nhận định thông thường rằng bò sát nhiệt đới hoạt động mạnh nhất trong mùa mưa. Về mặt sinh thái lý thuyết, điều này được giải thích bởi đặc thù của rừng khộp: vào mùa khô cây rụng lá làm tăng độ thông thoáng của tán rừng, buộc thằn lằn phải dồn mật độ về các trảng bụi ẩm và thảm lá khô còn sót lại để tìm kiếm con mồi và phơi nắng điều hòa thân nhiệt, khiến khả năng phát hiện tăng vọt.

4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) hoàn chỉnh không? Trả lời: Luận án cung cấp quy trình tái lặp chi tiết và minh bạch 100%: từ thông số thiết kế ô tiêu chuẩn ($100,\text{m}^2$, giãn cách 200m), công thức đo thể tích elip và tính IRI, quy trình cố định mẫu mô học bằng Bouin và nhuộm HE, cho đến chu trình nhiệt PCR (nhiệt độ gắn mồi, thời gian kéo dài) và mã hiệu primer của gen 16S rDNA đối chiếu GenBank AB057394.

5. Chương trình nghiên cứu 10 năm tiếp theo được phác thảo như thế nào? Trả lời: Chương trình 10 năm tập trung vào 3 trụ cột: (1) Thiết lập trạm quan trắc Occupancy dài hạn đa mùa tại Tây Nguyên để theo dõi biến động quần thể dưới tác động El Niño/La Niña; (2) Mở rộng giải trình tự toàn bộ hệ gen ty thể và hệ gen nhân liên kết cảnh quan trên toàn Đông Dương; (3) Thử nghiệm mô hình bảo tồn ngoại vi (ex-situ) và phát triển các chỉ số bò sát làm sinh vật chỉ thị chất lượng môi trường đất rừng khộp.

Kết luận

  1. Luận án đã hệ thống hóa và cung cấp cơ sở dữ liệu toàn diện, chuẩn xác đầu tiên về hình thái học, sai khác giới tính ($\text{SSD} = 0,012$) 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ồ.
  2. Xác định mật độ quần thể đạt 14 cá thể/ha tại VQG Yok Don (vùng lõi 15 cá thể/ha, vùng đệm 12 cá thể/ha, $P = 0,036$) và định lượng cấu trúc sử dụng 6 kiểu vi môi trường với sự chiếm ưu thế tuyệt đối của trảng cây bụi ($44,53%$) và thảm lá khô ($38,47%$).
  3. Tiên phong áp dụng mô hình chiếm cứ điểm đơn mùa MacKenzie kết hợp tiêu chuẩn AIC, chứng minh sinh cảnh rừng khộp cùng tổ hợp nhiệt - ẩm tối ưu hóa xác suất chiếm cứ điểm thực tế ($\psi = 0,4723$ mùa mưa; $\psi = 0,6054$ mùa khô), vượt trội so với tỉ suất chiếm cứ thuần túy.
  4. Làm rõ đặc điểm sinh thái dinh dưỡng với 17 loại thức ăn, trong đó 7 nhóm mồi then chốt chiếm $\text{IRI} = 77,43%$; chứng minh kích thước con mồi tương quan thuận chặt chẽ với chiều dài thân (SVL) và độ rộng miệng (MW).
  5. Xác định chu kỳ sinh sản phân mùa từ tháng III đến tháng VII thông qua phân tích mô học ống sinh tinh và sự đồng pha của quá trình phát triển buồng trứng với thể tích gan tích lũy vitellogenesis.
  6. Đánh giá đa dạng di truyền phân tử qua trình tự 16S rDNA, phát hiện 8 haplotype/16 mẫu với mức độ phân hóa di truyền từ $0,14%$ đến $2,66%$, chứng minh quần thể Tây Nguyên tiến hóa theo hướng chọn lọc ngẫu nhiên trung tính và đang trong pha mở rộng quần thể.