Tổng quan về luận án

Hệ sinh thái núi đá vôi (karst) tại khu vực Đông Nam Á từ lâu đã được giới sinh học quốc tế thừa nhận là "phòng thí nghiệm thiên nhiên lý tưởng cho các nghiên cứu về phân loại học, sinh thái học, tiến hóa và địa lý động vật học" (Viossanges, 2017). Mặc dù lớp Bò sát (Reptilia) trên toàn cầu đã chạm mốc 11.570 loài vào tháng 5/2021 (Uetz et al., 2021), song "ước tính có khoảng 20% số loài bò sát trên toàn cầu bị đe dọa tuyệt chủng; bên cạnh đó công tác bảo tồn gặp nhiều khó khăn do hiểu biết về sự đa dạng, đặc điểm sinh học - sinh thái học còn hạn chế" (Boehm et al., 2013). Nước Cộng hòa Dân chủ Nhân dân Lào sở hữu diện tích núi đá vôi khoảng 11.344 km² (tương đương 1,134 triệu ha, chiếm 5% diện tích cả nước), tập trung chủ yếu ở miền Trung và miền Bắc, nhưng lịch sử nghiên cứu bò sát tại đây từng ghi nhận những khoảng trống dữ liệu kéo dài: từ cá thể Gekko gecko đầu tiên được Linnaeus ghi nhận năm 1758, phải đến năm 1899 loài thứ hai là Dixonius siamensis mới được Boulenger mô tả.

+----------------------------------------------------------------------------------------------------+
|                                    KHUNG TỔNG QUAN LUẬN ÁN TIẾN SĨ                                |
|                                                                                                    |
|  Tác giả: Saly Sitthivong (2022) | Đơn vị: Trường Đại học Lâm nghiệp Việt Nam (VNUF)               |
|  Chuyên ngành: Quản lý Tài nguyên Rừng (Mã số: 9 62 02 11)                                          |
|  Người hướng dẫn: PGS. Lưu Quang Vinh, GS. Nguyễn Thế Nhã                                          |
|  Đối tượng & Địa bàn: Họ Tắc kè (Gekkonidae) tại 6 tỉnh núi đá vôi miền Bắc & Trung Lào           |
|  Quy mô thực địa: 112 ngày, 7 đợt, 13 điểm, 38 tuyến | Mẫu vật: 138 cá thể, 34 mẫu phân tử DNA    |
+----------------------------------------------------------------------------------------------------+

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng: trong khi các quốc gia láng giềng như Việt Nam và Thái Lan đã công bố tới 94 loài thuộc họ Tắc kè (Gekkonidae) với hệ thống cơ sở dữ liệu phân tử hoàn chỉnh, thì tài liệu tại Lào phân tán, thiếu vắng các khảo sát định lượng đa tầng và chưa tích hợp phát sinh chủng loại phân tử cho toàn bộ các giống phân bố ở các khối karst biệt lập (Stuart et al., 1999; Teynié et al., 2004; Luu et al., 2016).

Luận án của nghiên cứu sinh Saly Sitthivong giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm định 4 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Thành phần loài và mức độ phong phú của họ Gekkonidae tại các khối núi đá vôi miền Bắc và Trung Lào thực tế đạt quy mô nào?
  • RQ2: Mối quan hệ phát sinh chủng loại phân tử giữa các quần thể thuộc giống Thằn lằn ngón (Cyrtodactylus), Thằn lằn chân lá (Dixonius) và Tắc kè (Gekko) được cấu trúc ra sao dựa trên chỉ thị gen ty thể?
  • RQ3: Mức độ tương đồng thành phần loài biến thiên như thế nào theo không gian địa lý (Bắc - Trung Lào) và các lát cắt sinh cảnh, đai độ cao, vi sinh cảnh?
  • RQ4: Những áp lực nhân sinh nào đang đe dọa các loài tắc kè karst đặc hữu và đâu là ma trận giải pháp bảo tồn tối ưu?
  • H1: Tính đa dạng loài họ Gekkonidae ở Lào bị đánh giá thấp nghiêm trọng do hiện tượng loài ẩn (cryptic species) và các khối karst hoạt động như những "đảo lục địa" thúc đẩy phân hóa loài cục bộ.
  • H2: Sự tương đồng thành phần loài giữa các khối karst suy giảm tỷ lệ nghịch với khoảng cách địa lý và bị chia cắt bởi các rào cản thủy văn lớn (như lưu vực sông Mê Kông).
  • H3: Các loài tắc kè karst có sự phân tầng sinh thái nghiêm ngặt theo đai độ cao và vi sinh cảnh (vách đá, thân cây, nền đất) nhằm giảm thiểu cạnh tranh cùng loài và khác loài.
  • H4: Các can thiệp nhân sinh từ khai thác đá xi măng và canh tác nương rẫy là nhân tố chi phối chính làm suy thoái các điểm nóng đa dạng sinh học tắc kè.

Dựa trên Khung phân loại học tích hợp (Integrative Taxonomy Framework) và Thuyết địa lý sinh vật đảo (Island Biogeography Theory), công trình mang lại tác động đột phá khi định loại 28 loài thuộc 6 giống từ 138 mẫu vật thu thập trong 112 ngày thực địa (2018–2020) tại 13 điểm thuộc 6 tỉnh (Viêng Chăn, Luông Pha Bang, Khăm Muôn, U Đôm Xay, Húa Phăn, Xiêng Khoảng). Đóng góp mang tính lịch sử là việc phát hiện và công bố 06 loài mới cho khoa học: Cyrtodactylus houaphanensis, Cyrtodactylus muangfuangensis, Cyrtodactylus ngoiensis, Dixonius lao, Dixonius somchanhae, Gekko khunkhamensis, cùng 02 ghi nhận phân bố mới cho tỉnh U Đôm Xay (Cyrtodactylus wayakonei, Hemiphyllodactylus kiziriani).

Literature Review và Positioning

Các công trình nghiên cứu bò sát tại Lào phát triển qua nhiều giai đoạn nhưng không đồng đều. Stuart (1999) mở đầu với danh lục 109 loài, tiếp nối bởi Teynié et al. (2004) bổ sung 29 loài, Stuart & Heatwole (2008) ghi nhận thêm 13 loài, và Uetz et al. (2021) nâng tổng số bò sát lên 225 loài. Riêng họ Gekkonidae ghi nhận bước nhảy vọt: "từ năm 2004 tới năm 2020 chỉ khoảng 16 năm, các nhà nghiên cứu về bò sát đã phát hiện 30 loài trong họ Tắc kè tại Lào, trong đó có 20 loài là nằm trong giống Thằn lằn ngón Cyrtodactylus" (trích luận án).

Tiến trình ghi nhận số loài Bò sát tại CHDCND Lào:
1999 (Stuart et al.)         : 109 loài
2004 (Teynié et al.)          : 138 loài (+29)
2008 (Stuart & Heatwole)      : 151 loài (+13)
2010 (Uetz et al.)            : 157 loài (+6)
2018 (Uetz et al.)            : 181 loài (+24)
2020 (Uetz et al.)            : 212 loài (+31)
2021 (Uetz et al.)            : 225 loài (+13)

Tổng hợp các dòng nghiên cứu phân tử chính trong khu vực:

  1. Dòng nghiên cứu phân loại giống Cyrtodactylus: Nazarov et al. (2014) phân tích 673 bp gen COI xác lập 3 nhóm loài (nhóm A: C. phongnhakebangensis; nhóm B: C. wayakonei; nhóm C: C. irregularis). Luu et al. (2016) mở rộng tập mẫu với 668 bp gen COI, phân định các nhánh Bắc Lào - Tây Bắc Việt Nam và Trung Lào - Bắc Trung Bộ Việt Nam. Ngô Thị Hạnh et al. (2017) với 226 mẫu vật Đông Dương đã phân chia giống này thành 5 nhóm đơn ngành (A, B, C, D, E).
  2. Dòng nghiên cứu giống Gekko: Luu et al. (2015) kết hợp hình thái và sinh học phân tử xác lập các tổ hợp loài nhóm Gekko japonicus (G. bonkowskii, G. sengchanthavongi) và nhóm Gekko petricolus (G. boehmei) tại tỉnh Khăm Muôn.
  3. Dòng nghiên cứu đa gen khu vực: Brennan et al. (2017) tích hợp hai đoạn gen ty thể COI và ND2 trên các mẫu thu từ Ấn Độ, Myanmar, Malaysia, Thái Lan và Đông Dương, khẳng định sự chia cắt sâu sắc giữa các nhánh Nam - Trung - Bắc Đông Dương.

Trong y văn tồn tại hai cuộc tranh luận lý thuyết lớn:

  • Tranh luận 1 (Đồng danh hình thái vs. Đơn vị tiến hóa độc lập): Ngô Thị Hạnh et al. (2017) chứng minh C. paradoxus là đồng danh của C. condorensis, và C. thuongae là đồng danh của C. roesleri. Ngược lại, nhiều tác giả quốc tế cảnh báo việc dựa thuần túy vào tương đồng hình thái ngoài dễ bỏ sót các loài ẩn thích nghi cao độ trên từng tháp đá vôi cô lập.
  • Tranh luận 2 (Độ bao phủ của chỉ thị đơn gen ty thể): Một số ý kiến cho rằng chỉ thị COI có thể phản ánh sai lệch do hiện tượng bắt giữ dòng gen ty thể (mitochondrial introgression), trong khi trường phái phân loại học tích hợp (Padial et al., 2010) khẳng định COI hoàn toàn đủ độ nhạy để phân tách loài khi kết hợp chặt chẽ với các sai khác hình thái gián đoạn (discrete morphological gaps).

So sánh với các nghiên cứu quốc tế: Tại Việt Nam, hệ thống bò sát karst được khảo sát bài bản tại Phong Nha - Kẻ Bàng (Ziegler et al., 2013, 2016), ghi nhận 48 loài Cyrtodactylus và 17 loài Gekko. Tại Thái Lan, công trình của Taylor (1963) và các cập nhật của Uetz et al. (2021) xác nhận 38 loài Cyrtodactylus và 15 loài Gekko. Luận án của Saly Sitthivong định vị chính xác vào "khoảng trống xám" địa lý nối liền giữa khối karst Trường Sơn (Việt Nam) và các cao nguyên đá vôi phía Bắc Thái Lan, lấp đầy mắt xích tiến hóa quan trọng của khu hệ Gekkonidae Đông Dương.

Đó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 trực tiếp ba lý thuyết sinh học tiến hóa nền tảng:

  1. Lý thuyết hình thành loài trên đảo karst (Karst Inselberg Speciation Hypothesis - Clements et al., 2006; Grismer et al., 2018): Mở rộng luận điểm cho rằng các khối núi đá vôi nhiệt đới đóng vai trò như các đảo sinh thái biệt lập trên đất liền. Sự phân cắt địa hình ngăn cản dòng gen, dẫn đến hiện tượng biệt hóa ngoại vi (allopatric speciation) và tạo ra mức độ đặc hữu cục bộ cực cao ở các giống Cyrtodactylus và Dixonius.
  2. Khái niệm loài thống nhất (Unified Species Concept - de Queiroz, 2007) và Khung phân loại học tích hợp (Integrative Taxonomy - Dayrat, 2005; Padial et al., 2010): Xác lập ranh giới loài dựa trên sự hội tụ đa chiều giữa độ phân kỳ di truyền vượt ngưỡng (uncorrected p-distance > 6–10% ở gen COI) và các chỉ số hình thái gián đoạn (meristic and morphometric diagnostic traits).
  3. Thuyết phân chia ổ sinh thái (Niche Partitioning Theory - MacArthur & Levins, 1967): Chứng minh sự chung sống hòa hợp (co-existence) của nhiều loài tắc kè sympatric thông qua sự phân hóa kích thước cơ thể, cấu trúc giác bám và tập tính lựa chọn vi sinh cảnh (vách đá đứng, cành cây cao, thảm lá mục chân núi).
SƠ ĐỒ KHUNG MÔ HÌNH PHÂN HÓA LOÀI TRÊN NÚI ĐÁ VÔI LÀO:
+-------------------------------------------------------------------------+
|                  CÁC KHỐI NÚI ĐÁ VÔI BIỆT LẬP (KARST INSELBERGS)        |
|  - Rào cản địa hình thung lũng & mạng lưới sông suối                     |
|  - Vi khí hậu hang động và vách đá đặc thù                              |
+-------------------------------------------------------------------------+
+-------------------------------------------------------------------------+
|                      CƠ CHẾ PHÂN HÓA TIẾN HÓA                           |
|  - Cách ly địa lý (Geographic Isolation)                                |
|  - Chọn lọc thích nghi vi sinh cảnh (Vách đá vs Thân cây vs Mặt đất)   |
|  - Trôi dạt di truyền (Genetic Drift) & Đột biến gen COI                 |
+-------------------------------------------------------------------------+
+-------------------------------------------------------------------------+
|                  KẾT QUẢ ĐA DẠNG HÓA SINH HỌC                           |
|  - 06 loài mới công bố (*C. houaphanensis*, *D. lao*, *G. khunkhamensis*...) |
|  - Phân tầng thích nghi theo đai cao (<200m -> >1200m)                  |
|  - Tính đặc hữu vi mô cục bộ (Micro-endemism)                           |
+-------------------------------------------------------------------------+

Các mệnh đề lý thuyết (Theoretical Propositions) được xác lập:

  • Proposition 1 (P1): Khoảng cách phân kỳ di truyền giữa các đơn vị phân loại Gekkonidae tỷ lệ thuận với mức độ phân mảnh của khối núi đá vôi, trong đó các giống có độ linh động thấp (Dixonius, Hemiphyllodactylus) biểu hiện mức độ đặc hữu vi mô chặt chẽ hơn các giống có khả năng phát tán rộng (Hemidactylus).
  • Proposition 2 (P2): Sự chồng lấn về vùng phân bố địa lý (sympatry) giữa các loài trong cùng một giống bắt buộc phải đi kèm với sự phân hóa triệt để về vi sinh cảnh hoặc kích thước hình thái nhằm triệt tiêu cạnh tranh loại trừ.

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

Khung phân tích của luận án tích hợp một cách đa tầng:

  • Hệ thống 39 chỉ tiêu hình thái học: Gồm 19 chỉ số đo kích thước (SVL, TaL, AG, HL, HW, HH, SE, EyeEar, ForeL, FemurL, CrusL, LD4A, LD4P, OD, EarL, RW, RH, MW, ML) đo bằng thước kẹp điện tử chính xác 0,01 mm và 20 chỉ số đếm vảy (CS, N, I, SL, IL, IO, PO, PM, GP, DTR, GSDT, SBL, SR, V, LF1, LF4, LT1, LT4, PP, PAT) dưới kính hiển vi điện tử.
  • Phân tích phả hệ đa giải thuật: Đối chiếu song song Tiết kiệm tối đa (Maximum Parsimony - MP), Khả năng tối đa (Maximum Likelihood - ML) và Suy luận Bayes (Bayesian Inference - BI).
  • Phân tích cấu trúc sinh thái quần xã: Ứng dụng hệ số tương đồng Sorensen-Dice $d_{jk} = \frac{2M}{2M+N}$ kết hợp phân tích cụm (Cluster Analysis) trên phần mềm PAST để lượng hóa sự biến thiên thành phần loài giữa các vùng địa lý.
  • Điều kiện biên (Boundary Conditions): Khung phân tích được chuẩn hóa cho hệ sinh thái rừng trên núi đá vôi nhiệt đới gió mùa tại Lào (độ cao từ 150 m đến 1.397 m), không áp dụng suy diễn đồng nhất cho các hệ sinh thái rừng trên đất cát hoặc rừng ngập nước.

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

Thiết kế nghiên cứu

Luận án quán triệt lập trường triết học Thực chứng sau hiện đại / Hiện thực phê phán (Critical Realism), kết hợp chặt chẽ giữa quan sát thực nghiệm khách quan và suy luận mô hình hóa phân tử. Thiết kế nghiên cứu đa cấp độ (Multi-level design) được triển khai:

  • Cấp độ vĩ mô (Macro-level): Khảo sát so sánh giữa 2 miền sinh địa lý lớn: Miền Bắc Lào (Luông Pha Bang, U Đôm Xay, Húa Phăn, Xiêng Khoảng) và Miền Trung Lào (Viêng Chăn, Khăm Muôn).
  • Cấp độ trung mô (Meso-level): Phân tích qua 3 dạng sinh cảnh thảm thực vật:
    • SC1: Núi đá vôi thuộc khu canh tác nông nghiệp/công nghiệp xen kẽ.
    • SC2: Núi đá vôi thuộc rừng thứ sinh nghèo kiệt, chịu tác động.
    • SC3: Núi đá vôi thuộc rừng nguyên sinh giàu trữ lượng, ít bị tác động. Kèm theo 6 đai độ cao (cách nhau 200 m): <200 m, 200–<400 m, 400–<600 m, 600–<800 m, 800–<1.000 m, 1.000–<1.200 m, và >1.200 m.
  • Cấp độ vi mô (Micro-level): Phân loại chi tiết theo 3 vị trí bắt gặp: Vách đá (khe, lỗ hốc đá), Thân cây (thân, cành, dây leo), và Mặt đất (thảm mục, chân tảng đá).
MA TRẬN ĐỊA BÀN THỰC ĐỊA TẠI 6 TỈNH NÚI ĐÁ VÔI LÀO:

Quy trình nghiên cứu rigorous

  1. Khảo sát thực địa và thu mẫu:
    • Tiến hành theo phương pháp điều tra theo tuyến và điểm có định vị GPS (trackmaker).
    • Khung giờ khảo sát chuẩn hóa ban đêm: 17h00 đến 23h00 (khung giờ hoạt động tối ưu của Gekkonidae).
    • Dụng cụ bắt mẫu: Bắt tay trực tiếp hoặc dùng kẹp chuyên dụng bọc cao su chống tổn thương cơ học.
    • Xử lý tiêu bản: Gây mê nhân đạo bằng bông thấm ethylacetate trong bình kín theo quy chuẩn quốc tế (Simmons, 2002); gắn nhãn mã số trường; cố định trong cồn 90% từ 5–8 giờ; chuyển sang lưu trữ vĩnh viễn trong cồn 70% tại Bộ môn Động vật rừng (VNUF) và Đại học Quốc gia Lào.
  2. Chiết tách DNA và khuếch đại PCR:
    • Mẫu mô cơ/gan/đuôi sâu được bảo quản ở -20°C. Sử dụng kit DNeasy Blood and Tissue (Qiagen, Đức) cho mẫu khó và GeneJET Genomic DNA Purification (ThermoFisher, Lithuania) cho mẫu tươi.
    • Cặp mồi đặc hiệu:
      • Giống Cyrtodactylus: VF1-d (5’-TTCTCAACCAACCACAARGAYATYGG-3’) và VR1-d (5’-TAGACTTCTGGGTGGCCRAARAAYCA-3’) (Ivanova et al., 2006).
      • Giống Dixonius: MetF1 (5’-AAGCTTTCGGGCCCATACC-3’) và COIR1 (5’-AGRGTGCCAATGTCTTTGTGRTT-3’) (Macey et al., 1997).
      • Giống Gekko: GF1 (5'-CAAGCACHATYATYACYATAT-3') và GR1 (5'-CCTATGTGTGCGATTGATGA-3') (Greenbaum et al., 2007).
    • Hỗn hợp phản ứng (21 µl): 1–2 µl DNA khuôn, 2 µl mỗi mồi (10 µM), 5 µl H₂O khử ion, 10 µl HotStarTaq Mastermix (Qiagen) hoặc DreamTaq Mastermix (ThermoFisher). Chu trình nhiệt: Biến tính ban đầu 95°C (5–15 phút); 35 chu kỳ (95°C/30 giây, bắt mồi 48°C–60°C/45 giây, kéo dài 72°C/1 phút); kéo dài cuối 72°C/6 phút. Có đối chứng âm ($ddH_2O$) trong mọi mẻ chạy.
  3. Giải trình tự và kiểm soát chất lượng:
    • Điện di kiểm tra trên gel agarose 1%, đệm TBE 1X, nhuộm ethidium bromide 2 µg/ml ở 90V trong 30 phút, hiển thị dưới tia cực tím trên máy Alphamager MINI (Protein Simple, Mỹ).
    • Tinh sạch sản phẩm PCR bằng kit GeneJET PCR Purification. Giải trình tự Sanger hai chiều tại FirstBase (Malaysia).

Data và phân tích

  • Xử lý chuỗi nucleotide: Kiểm tra tín hiệu chromatogram bằng Sequencher v5.1; truy vấn BLAST trên NCBI GenBank; căn chỉnh đa chuỗi bằng Clustal X v2.0.
  • Mô hình tiến hóa và xây dựng cây phát sinh: Lựa chọn mô hình thay thế nucleotide tối ưu TVM+I+G thông qua jModelTest v2.4.
  • Phân tích thống kê phát sinh loài:
    • Bayesian Inference (BI): Thực hiện trên MrBayes v3.7 với $1 \times 10^7$ thế hệ (10 triệu thế hệ MCMC), điểm cắt bỏ burn-in là 54.000 thế hệ sau khi hàm log-likelihood ($-\ln L$) đạt trạng thái dừng hội tụ ổn định.
    • Maximum Likelihood (ML): Thực hiện trên IQ-TREE v1.1 với 10.000 bản sao bootstrap siêu nhanh (Ultrafast Bootstrap - UFB).
    • Maximum Parsimony (MP): Thực hiện trên PAUP* 4.0b10 với 1.000 bản sao bootstrap chuẩn (BP).
    • Tiêu chuẩn hỗ trợ nhánh: Nút nhánh được coi là có độ tin cậy vững chắc khi đạt đồng thời $BP \ge 70%$, $UFBP \ge 95%$ và Xác suất hậu nghiệm Bayes $PP \ge 0,95$ (Hillis & Bull, 1993; Ronquist et al., 2012). Khoảng cách di truyền p-distance không hiệu chỉnh (uncorrected pairwise divergence) được trích xuất trực tiếp qua PAUP*.

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

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

+----------------------------------------------------------------------------------------------------+
|                       DANH LỤC 06 LOÀI MỚI CÔNG BỐ CHO KHOA HỌC TRONG LUẬN ÁN                      |
+----+-------------------------------+-------------------+-------------------------------------------+
| STT| Tên loài mới                  | Địa bàn phát hiện | Đặc điểm chẩn loại nổi bật                |
+----+-------------------------------+-------------------+-------------------------------------------+
| 1  | Cyrtodactylus houaphanensis   | Hủa Phăn          | Vảy đùi/bụng chuyên biệt, hoa văn lưng độc|
| 2  | Cyrtodactylus muangfuangensis | Viêng Chăn        | Sai khác số lỗ trước huyệt, 39 chỉ số đo  |
| 3  | Cyrtodactylus ngoiensis       | Luông Pha Bang    | Hình thái ngón & đai vảy lưng đặc trưng   |
| 4  | Dixonius lao                  | Viêng Chăn        | Tấm cằm & vảy bụng phân tầng rõ rệt       |
| 5  | Dixonius somchanhae           | Viêng Chăn        | Sai biệt di truyền COI sâu sắc (>8%)      |
| 6  | Gekko khunkhamensis           | Khăm Muôn         | Cỡ cơ thể lớn, nốt sần lưng phân bố dày   |
+----+-------------------------------+-------------------+-------------------------------------------+
  1. Phát hiện và công bố 06 loài mới cho khoa học: Xác lập tên khoa học hợp lệ cho 3 loài Thằn lằn ngón (Cyrtodactylus houaphanensis, Cyrtodactylus muangfuangensis, Cyrtodactylus ngoiensis), 2 loài Thằn lằn chân lá (Dixonius lao, Dixonius somchanhae) và 1 loài Tắc kè (Gekko khunkhamensis).
  2. Mở rộng vùng phân bố địa lý của các loài quý hiếm: Ghi nhận lần đầu tiên sự hiện diện của Cyrtodactylus wayakonei và Hemiphyllodactylus kiziriani tại tỉnh U Đôm Xay, chứng minh cầu nối sinh học liền mạch giữa khu hệ Nam Trung Quốc, Bắc Lào và Tây Bắc Việt Nam.
  3. Cấu trúc đa dạng thành phần loài theo địa phương:
    • Tỉnh Viêng Chăn: Đa dạng cao nhất với 14 loài thuộc 5 giống (4 loài Cyrtodactylus, 3 loài Dixonius, 4 loài Hemidactylus, 2 loài Gehyra, 1 loài Gekko).
    • Tỉnh Luông Pha Bang: Ghi nhận 8 loài thuộc 6 giống.
    • Tỉnh Khăm Muôn: Nổi bật với các loài Gekko và Cyrtodactylus đặc hữu karst thuần túy.
    • Tỉnh U Đôm Xay: Ghi nhận 5 loài thuộc 4 giống.
  4. Quy luật phân tầng sinh thái và vi sinh cảnh:
    • Theo sinh cảnh: Sinh cảnh rừng nguyên sinh trên núi đá vôi (SC3) lưu giữ số lượng loài cao nhất và là nơi duy nhất phát hiện các loài đặc hữu hẹp; sinh cảnh SC2 suy giảm 40% độ phong phú; sinh cảnh SC1 chủ yếu là các loài thạch sùng phổ thông (Hemidactylus frenatus, H. platyurus).
    • Theo đai cao: Mật độ cá thể và số lượng loài tập trung cao nhất ở đai cao từ 200 m đến 600 m. Trên 800 m độ phong phú loài giảm mạnh, và trên 1.200 m chỉ bắt gặp đơn lẻ một số loài chịu lạnh thích nghi cao tại Húa Phăn và Xiêng Khoảng.
    • Theo vi sinh cảnh: Tắc kè karst phân hóa sâu sắc: 64% mẫu bắt gặp trên vách đá (khe nứt đá vôi thẳng đứng); 24% trên thân cây/cành cây trong thung lũng đá vôi; 12% ở mặt đất/chân tảng đá.
PHÂN BỐ TỈ LỆ BẮT GẶP VI SINH CẢNH CỦA HỌ GEKKONIDAE TẠI LÀO:
Vách đá karst (khe nứt, trần hang) : [████████████████████████████████] 64%
Thân cây & Dây leo trong rừng karst : [████████████] 24%
Mặt đất & Chân khối đá vôi         : [██████] 12%
  1. Độ phân kỳ di truyền cao bất thường giữa các quần thể: Cây phát sinh chủng loại BI/ML chỉ ra khoảng cách di truyền không hiệu chỉnh uncorrected p-distance giữa các loài mới và các loài chị em dao động từ 7,8% đến 18,2% trên đoạn gen COI, vượt xa ngưỡng phân định loài tiêu chuẩn (threshold 5–6%) trong phân loại học bò sát Squamata, khẳng định tính độc lập tiến hóa không thể tranh cãi.

Implications đa chiều

  • Về lý thuyết sinh học: Cung cấp bằng chứng thuyết phục chứng minh các dãy núi đá vôi Đông Dương là "bộ máy phát sinh loài" (speciation pump) chứ không đơn thuần là nơi trú ẩn thụ động (refugia).
  • Về phương pháp luận: Chuẩn hóa quy trình tích hợp giữa đo đếm hình thái vi mô 39 chỉ số và quy trình PCR - Sanger hai chiều cho các nghiên cứu lưỡng cư - bò sát tại các quốc gia đang phát triển có nguồn lực hạn chế.
  • Về ứng dụng thực tiễn: Thiết lập cơ sở dữ liệu định danh hình thái kèm mã vạch DNA (DNA barcoding) giúp cơ quan kiểm lâm và hải quan Lào phát hiện các vụ buôn bán trái phép tắc kè qua biên giới.
  • Về chính sách và quản trị: Cung cấp luận cứ khoa học cập nhật Danh lục Đỏ IUCN, bổ sung các loài mới vào Nghị định 08/CP/2021 của Chính phủ Lào và Luật Động vật hoang dã số 07/QH/2007; đồng thời đề xuất phân vùng vùng đệm bảo vệ nghiêm ngặt xung quanh các mỏ khai thác đá vôi xi măng.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu xuất phát từ bối cảnh thực địa:

  1. Độ bao phủ địa lý: Mặc dù đã khảo sát 13 điểm tại 6 tỉnh trọng điểm, nhưng dải karst rộng lớn tại các tỉnh Nam Lào (như Savannakhet, Attapeu, Champasak) và các vùng giáp biên giới hẻo lánh chưa được khảo sát do trở ngại địa hình và bom mìn sót lại sau chiến tranh.
  2. Chỉ thị phân tử đơn locus: Nghiên cứu sử dụng chỉ thị gen ty thể COI. Dù COI phân định loài xuất sắc, việc thiếu các gen nhân (nuclear genes như RAG1, PDC, CMOS) hạn chế khả năng giải quyết triệt để các phân nhánh sâu ở cấp độ giống hoặc phát hiện hiện tượng lai tự nhiên xa xưa.
  3. Một số đơn vị phân loại chưa định danh đến loài: Một số taxon như Cyrtodactylus sp., Gehyra sp., Hemiphyllodactylus sp. thu được số lượng mẫu hạn chế (1–2 cá thể), chưa đủ chuỗi thống kê hình thái để công bố loài mới và tạm thời lưu giữ dưới dạng danh pháp mở.
  4. Thiếu hụt dữ liệu động thái quần thể dài hạn: Dữ liệu sinh thái dựa trên các chuyến thám hiểm cắt ngang (cross-sectional surveys), chưa theo dõi được biến động mật độ cá thể theo chu kỳ nhiều năm dưới tác động của biến đổi khí hậu.

Chương trình nghiên cứu tiếp nối (Future Research Agenda) gồm 5 hướng trọng tâm:

  • Hướng 1: Mở rộng giải trình tự hệ gen thế hệ mới (Next-Generation Sequencing - RAD-seq hoặc Target Capture UCEs) để tái cấu trúc toàn diện cây phát sinh chủng loại các giống Cyrtodactylus và Dixonius toàn Đông Dương.
  • Hướng 2: Nghiên cứu sinh lý học nhiệt thích ứng (Thermal physiology) của các loài tắc kè chuyên biệt hang động karst trước kịch bản nhiệt độ toàn cầu tăng cao.
  • Hướng 3: Khảo sát toàn diện hệ karst miền Nam Lào và các vùng hành lang đa dạng sinh học xuyên biên giới Lào - Việt Nam - Thái Lan.
  • Hướng 4: Đánh giá độc chất học và tác động ô nhiễm bụi đá vôi từ các nhà máy xi măng đến khả năng sinh sản của các loài tắc kè vách đá.
  • Hướng 5: Triển khai nghiên cứu bảo tồn chuyển chỗ (ex-situ) và sinh sản nhân tạo bảo tồn các loài đặc hữu cực kỳ nguy cấp có phạm vi phân bố hẹp dưới 10 km².

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

+----------------------------------------------------------------------------------------------------+
|                               CÁC TRỤ CỘT TÁC ĐỘNG CỦA LUẬN ÁN                                      |
|                                                                                                    |
|  [HỌC THUẬT QUỐC TẾ]      [CHÍNH SÁCH QUỐC GIA]       [CÔNG NGHIỆP - RỪNG]    [BẢO TỒN TOÀN CẦU]   |
|  Công bố 6 loài mới       Cập nhật Luật ĐVHD Lào      Quy hoạch mỏ đá vôi     Bổ sung RedList IUCN |
|  Cơ sở dữ liệu GenBank    Nghị định 08/CP/2021        ĐTM các dự án xi măng   Bảo tồn cảnh quan    |
|  Trích dẫn Zootaxa/VNUF   Bảo vệ VQG Nam Et-PhouLouey Giám sát đa dạng rừng   karst Đông Dương     |
+----------------------------------------------------------------------------------------------------+
  • Tác động học thuật: Tăng cường vị thế khoa học của Trường Đại học Lâm nghiệp Việt Nam và Đại học Quốc gia Lào trên bản đồ nghiên cứu phân loại học quốc tế; cung cấp 34 trình tự mã vạch chuẩn cho Ngân hàng Gen thế giới (GenBank/NCBI), phục vụ hàng trăm công trình nghiên cứu phát sinh loài tiếp theo.
  • Tác động công nghiệp và tài nguyên: Buộc các tập đoàn công nghiệp khai khoáng và sản xuất xi măng tại Viêng Chăn và Khăm Muôn phải thực hiện Đánh giá tác động môi trường (ĐTM) nghiêm ngặt trước khi nổ mìn phá vỡ các tháp karst nguyên sinh.
  • Ảnh hưởng chính sách quốc gia: Trở thành tài liệu khoa học chính thống giúp Bộ Nông - Lâm nghiệp Lào (MAF) xây dựng Kế hoạch hành động đa dạng sinh học quốc gia; bảo vệ tính nguyên vẹn của các Vườn Quốc gia trọng yếu như Nam Et-Phou Louey, Hin Nam No.
  • Lợi ích xã hội và cộng đồng: Nâng cao nhận thức bảo tồn cho cộng đồng bản địa, giảm thiểu tập quán săn bắt tắc kè bừa bãi làm dược liệu ngâm rượu thương mại, bảo vệ cân bằng sinh thái kiểm soát côn trùng dịch hại nông - lâm nghiệp.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận quy trình thực nghiệm tích hợp chuẩn mực từ thiết kế khảo sát thực địa ban đêm, quy chuẩn đo đếm 39 chỉ số hình thái, đến phân tích tin sinh học phát sinh loài với các dòng lệnh chuyên sâu trên MrBayes và IQ-TREE.
  • Giới nghiên cứu Lưỡng cư - Bò sát học (Herpetologists): Sở hữu bộ dữ liệu so sánh chuẩn xác về các loài Gekkonidae phân bố tại Lào trong mối tương quan với khu hệ Việt Nam, Thái Lan, Campuchia, Myanmar và Trung Quốc.
  • Ban quản lý các Khu Bảo tồn & Vườn Quốc gia tại Lào: Nắm bắt danh mục các loài đặc hữu, tọa độ phân bố chính xác và 5 tiêu chí ưu tiên bảo tồn (với thang điểm 1–13) để phân bổ nguồn lực tuần tra hiệu quả.
  • Các tổ chức Bảo tồn Quốc tế (IUCN, WWF, WCS, Cologne Zoo): Sử dụng dữ liệu thực chứng để đánh giá chính xác tình trạng nguy cấp của các loài tắc kè karst hẹp, từ đó tài trợ các dự án bảo tồn khẩn cấp.

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?

Đóng góp lý thuyết độc đáo nhất là việc chứng minh và hoàn thiện Lý thuyết phân hóa loài trên đảo karst (Karst Inselberg Speciation) đối với họ Gekkonidae tại lục địa Đông Nam Á. Luận án chỉ ra rằng các tháp đá vôi độc lập hoạt động như các hệ thống sinh thái đảo khép kín, nơi áp lực chọn lọc thích nghi vi sinh cảnh (vách đá đứng) kết hợp với rào cản phân cách địa hình thung lũng đã thúc đẩy quá trình phân nhánh tiến hóa (cladogenesis) nhanh chóng, tạo nên hàng loạt loài đặc hữu vi mô có khoảng cách di truyền COI vượt trên 7,8–18,2% so với các tổ tiên chung.

2. Điểm cải tiến đột phá về phương pháp luận so với các nghiên cứu trước đây là gì?

So với các nghiên cứu đơn lẻ trước đây (như Nazarov et al., 2014 chỉ phân tích 21 loài; Brennan et al., 2017 khảo sát mẫu diện rộng phân tán), luận án tạo đột phá bằng cách:

  • Kết hợp đồng thời 39 chỉ tiêu đo đếm hình thái vi mô với phân tích phân tử đa giải thuật (MP, ML, BI) trên mô hình tiến hóa tối ưu TVM+I+G được chọn lọc nghiêm ngặt qua jModelTest.
  • Phân tích định lượng tương đồng sinh thái thông qua chỉ số Sorensen-Dice ($d_{jk}$) xử lý trên phần mềm PAST, kết nối giữa địa lý phân tử và sinh thái học cảnh quan đa tầng (3 dạng sinh cảnh, 6 đai độ cao, 3 vi sinh cảnh).
So sánh phương pháp luận giữa các công trình nghiên cứu:

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện gây ngạc nhiên lớn nhất là sự hiện diện của 02 loài mới thuộc giống Thằn lằn chân lá (Dixonius lao, Dixonius somchanhae) tại các vùng núi đá vôi tỉnh Viêng Chăn. Trước nghiên cứu này, giống Dixonius tại Lào chỉ được biết đến duy nhất qua loài Dixonius siamensis (ghi nhận từ 1899) và thường được xem là nhóm tắc kè phân bố ở rừng đất bằng thấp hoặc savan thưa. Việc phát hiện các loài Dixonius đặc hữu chuyên biệt trên các vách đá vôi karst đảo lộn giả định trước đó về sinh thái học tiến hóa của giống này tại bán đảo Đông Dương.

4. Luận án có cung cấp quy trình tái lặp thực nghiệm (Replication Protocol) hoàn chỉnh không?

Hoàn toàn đầy đủ và chi tiết. Luận án cung cấp:

  • Tọa độ GPS chính xác đến từng giây của 13 điểm điều tra tại 6 tỉnh.
  • Trình tự nucleotide chuẩn của 3 cặp mồi PCR (VF1-d/VR1-d, MetF1/COIR1, GF1/GR1).
  • Thông số chu trình nhiệt, nồng độ hóa chất tách chiết, nồng độ mastermix và quy cách kiểm định điện di.
  • Các tham số dòng lệnh cài đặt phân tích phát sinh loài trong PAUP*, IQ-TREE, và MrBayes (số thế hệ $10^7$, burn-in 54.000, ngưỡng hỗ trợ $BP \ge 70$, $UFBP \ge 95%$, $PP \ge 0,95$).

5. Chương trình nghị sự nghiên cứu 10 năm (2022–2032) được phác thảo thế nào?

Lộ trình 10 năm xác định:

  • Giai đoạn 2022–2025: Hoàn thành phân loại học tích hợp cho các loài danh pháp mở (Cyrtodactylus sp., Hemiphyllodactylus sp.); công bố các đánh giá Danh lục Đỏ IUCN cho 6 loài mới.
  • Giai đoạn 2025–2028: Mở rộng điều tra toàn diện các khối karst Nam Lào và dãy Trường Sơn Trung - Nam; ứng dụng công nghệ giải trình tự toàn bộ hệ gen ty thể (Mitogenomics).
  • Giai đoạn 2028–2032: Thiết lập mạng lưới giám sát đa dạng sinh học karst xuyên quốc gia (Lào - Việt Nam - Thái Lan); xây dựng các khu bảo tồn vi mô (micro-reserves) chuyên biệt cho các loài tắc kè karst có nguy cơ tuyệt chủng cao.

Kết luận

Luận án tiến sĩ của Saly Sitthivong xác lập dấu mốc học thuật xuất sắc trong lĩnh vực Quản lý Tài nguyên Rừng và Động vật học tại Đông Dương qua 6 đóng góp cốt lõi:

  1. Phát hiện và công bố 06 loài mới cho khoa học (Cyrtodactylus houaphanensis, C. muangfuangensis, C. ngoiensis, Dixonius lao, D. somchanhae, Gekko khunkhamensis), nâng cao vượt bậc số lượng loài Gekkonidae được ghi nhận tại Lào.
  2. Bổ sung 02 ghi nhận phân bố địa lý mới (Cyrtodactylus wayakonei, Hemiphyllodactylus kiziriani) cho tỉnh U Đôm Xay, góp phần hoàn thiện bức tranh địa lý sinh vật học khu vực Tây Bắc Lào.
  3. Mô tả chi tiết đặc điểm hình thái học 39 chỉ tiêu chuẩn hóa cho toàn bộ 28 loài ghi nhận, thiết lập bộ dữ liệu tra cứu hình thái học chuẩn xác nhất cho bò sát karst Lào.
  4. Giải mã và tái lập cây quan hệ di truyền phân tử cho 3 giống chủ chốt (Cyrtodactylus, Dixonius, Gekko) dựa trên chỉ thị gen COI với các phân tích Bayes và Khả năng tối đa tin cậy cao.
  5. Định lượng hóa cấu trúc phân bố sinh thái của tắc kè karst theo 3 dạng thảm thực vật, 6 đai độ cao và 3 vi sinh cảnh, chỉ rõ mối liên hệ giữa rừng nguyên sinh vách đá (SC3) và tính đặc hữu sinh học.
  6. Xây dựng ma trận đánh giá điểm ưu tiên bảo tồn (thang điểm 1–13) dựa trên 5 tiêu chuẩn thực chứng, nhận diện các điểm nóng bị đe dọa nghiêm trọng bởi công nghiệp khai khoáng và nương rẫy.
+----------------------------------------------------------------------------------------------------+
|                         DI SẢN HỌC THUẬT VÀ TẦM ẢNH HƯỞNG TOÀN CẦU CỦA LUẬN ÁN                      |
|                                                                                                    |
|  1. Mở ra 3 nhánh nghiên cứu mới: Phân loại học bộ gen Karst - Sinh lý học nhiệt thích ứng        |
|     - Độc chất học sinh thái khai khoáng xi măng.                                                  |
|  2. Chuyển dịch mô thức: Khẳng định Karst Lào là trung tâm phát sinh loài độc lập tại Đông Nam Á.  |
|  3. Giá trị ứng dụng trường tồn: Bộ khung pháp lý và dữ liệu DNA phục vụ bảo tồn đa dạng sinh học |
|     quốc gia CHDCND Lào đến năm 2030 và tầm nhìn 2050.                                             |
+----------------------------------------------------------------------------------------------------+

Công trình không chỉ mở ra ba hướng nghiên cứu chuyên sâu về phân loại học tiến hóa, sinh thái cảnh quan karst và bảo tồn gen nguy cấp, mà còn là di sản khoa học nền tảng, định hình chiến lược bảo vệ tài nguyên rừng và thiên nhiên của CHDCND Lào trong nhiều thập kỷ tới.