Tổng quan về luận án

Nghiên cứu phân loại học và phát sinh chủng loại của khu hệ bò sát tại Đông Nam Á trong thập niên qua đã chứng kiến sự gia tăng đột biến về số lượng loài được phát hiện, từ 9.300 loài năm 2011 lên 11.570 loài vào tháng 05 năm 2021 (Uetz et al., 2021). Tuy nhiên, theo ước tính của Boehm et al. (2013), có khoảng 20% số loài bò sát trên toàn cầu đang đối mặt với nguy cơ tuyệt chủng, trong đó hệ sinh thái núi đá vôi (karst) nhiệt đới – vốn được xem là "phòng thí nghiệm tự nhiên" của tiến hóa và phân nhánh thích nghi (adaptive radiation) – đang chịu áp lực nhân sinh khốc liệt (Viossanges, 2017). Nước Cộng hòa Dân chủ Nhân dân Lào sở hữu 11.344 km² diện tích núi đá vôi (chiếm 5% tổng diện tích tự nhiên của cả nước) với tỷ lệ che phủ rừng đạt 62% (Bộ Nông-Lâm nghiệp Lào, 2019). Mặc dù vậy, dữ liệu khoa học về khu hệ tắc kè (Gekkonidae) tại Lào vẫn tồn tại khoảng trống rất lớn so với các quốc gia láng giềng như Việt Nam (94 loài, 8 giống) và Thái Lan (94 loài, 7 giống) (Uetz et al., 2021).

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ việc phần lớn các khối núi đá vôi biệt lập tại miền Bắc và miền Trung Lào chưa từng được khảo sát hệ thống về mặt hình thái học lẫn phát sinh chủng loại phân tử. Luận án tiến sĩ của Saly Sitthivong (2022) với đề tài "Nghiên cứu đa dạng thành phần loài và quan hệ di truyền của họ Tắc kè (Gekkonidae) ở một số khu vực núi đá vôi tại Nước Cộng hòa Dân chủ Nhân dân Lào" do PGS. Lưu Quang Vinh và GS. Nguyễn Thế Nhã hướng dẫn tại Trường Đại học Lâm nghiệp Việt Nam, là công trình tiên phong giải quyết toàn diện khoảng trống này. Luận án vận hành trên 4 câu hỏi nghiên cứu (Research Questions - RQ) và 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Đa dạng thành phần loài và cấu trúc quần xã họ Tắc kè tại các sinh cảnh núi đá vôi miền Bắc và miền Trung Lào phong phú ở mức độ nào?
  • RQ2: Mối quan hệ phát sinh chủng loại di truyền giữa các loài thuộc các giống Cyrtodactylus, Dixonius, và Gekko tại Lào được phân định như thế nào dựa trên chỉ thị phân tử ty thể?
  • RQ3: Mức độ tương đồng thành phần loài và sự phân hóa ổ sinh thái (phân bố theo đai cao, dạng sinh cảnh, vi sinh cảnh) chịu sự chi phối ra sao bởi tính phân mảnh của địa hình karst?
  • RQ4: Các yếu tố nhân sinh nào đang đe dọa trực tiếp đến quần thể tắc kè đặc hữu và đâu là ma trận ưu tiên bảo tồn tối ưu?
  • H1: Các khối karst biệt lập tại Lào đóng vai trò là "ốc đảo lục địa" (continental islands) thúc đẩy sự hình thành loài ngoại vi (allopatric speciation) và ẩn chứa các loài mới chưa được mô tả.
  • H2: Hình thái học truyền thống bị giới hạn bởi hiện tượng đồng hình (cryptic diversity), đòi hỏi phương pháp phân loại tích hợp (Integrative Taxonomy) để tái lập cây phả hệ chính xác.
  • H3: Sự tương đồng thành phần loài suy giảm theo khoảng cách địa lý và sự chia cắt của các lưu vực sông ngòi lớn giữa miền Bắc và miền Trung Lào.

Công trình mang lại đột phá với việc kiểm tra 138 mẫu vật hình thái, giải trình tự 34 mẫu DNA, phân tích 800 tư liệu ảnh thực địa qua 7 đợt khảo sát kéo dài 112 ngày tại 13 địa đ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). Kết quả đã xác lập danh lục 28 loài thuộc 6 giống, phát hiện và công bố 06 loài mới hoàn toàn 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 và Hemiphyllodactylus kiziriani).

Literature Review và Positioning

Lịch sử nghiên cứu bò sát tại Lào ghi nhận sự phát triển muộn và phân đoạn. Giai đoạn đầu, Stuart (1999) công bố danh lục bò sát đầu tiên của Lào với 109 loài. Tiếp đó, Teynié et al. (2004) bổ sung thêm 29 loài; Stuart & Heatwole (2008) ghi nhận thêm 13 loài; nâng tổng số lên 157 loài vào năm 2010. Từ năm 2010 đến 2021, số loài bò sát được ghi nhận tăng nhanh từ 181 loài (Uetz et al., 2018), 212 loài (Uetz et al., 2020) lên 225 loài (Uetz et al., 2021). Riêng đối với họ Gekkonidae, xuất phát điểm từ thế kỷ 18 chỉ có duy nhất loài Gekko gecko (Linnaeus, 1758), đến cuối thế kỷ 19 ghi nhận thêm Dixonius siamensis (Boulenger, 1899), và thế kỷ 20 bổ sung 5 loài. Giai đoạn bùng nổ diễn ra từ năm 2004 đến 2020 khi 30 loài mới được khám phá chỉ trong vòng 16 năm, trong đó 20 loài thuộc giống Cyrtodactylus (David et al., 2004; Ziegler et al., 2010; Schneider et al., 2011; Nazarov et al., 2014; Luu et al., 2014, 2015, 2016).

   ┌───────────────────────────────────────────────────────────┐
   │ TIẾN TRÌNH TÍCH LŨY SỐ LOÀI BÒ SÁT ĐƯỢC GHI NHẬN TẠI LÀO   │
   └───────────────────────────────────────────────────────────┘
   250 ─────────────────────────────────────────────────── 225 (2021)
                                                    212 (2020)
   200 ────────────────────────────── 181 (2018) ───┘
                               157 (2010) ──┘
   150 ───────── 138 (2004) ───┘
          109 (1999) ─┘
   100 ────┘
     0 ────┴──────────┴──────────┴──────────┴──────────┴──────────┴──
         Stuart     Teynié     Uetz       Uetz       Uetz       Uetz
         (1999)     (2004)     (2010)     (2018)     (2020)     (2021)

Về mặt học thuật, tồn tại những cuộc tranh luận lý thuyết sâu sắc về phân loại và phát sinh chủng loại của giống Cyrtodactylus tại Đông Dương. Quan điểm thứ nhất của Nazarov et al. (2014), dựa trên 673 bp gen COI ty thể của 21 loài tại Lào và Việt Nam, chia giống này thành 3 nhóm: nhóm A (C. phongnhakebangensis phân bố ở miền Trung Việt Nam và Trung-Đông Lào từ Khăm Muôn đến Viêng Chăn), nhóm B (C. wayakonei phân bố ở Tây-Bắc Lào và giáp ranh Trung Quốc), và nhóm C (C. irregularis đặc thù cho miền Trung Việt Nam). Trái lại, quan điểm thứ hai của Luu et al. (2016) trên 668 bp gen COI của 29 loài lại xác định: nhóm 1 (C. wayakonei ở Bắc Lào và Tây Bắc Việt Nam), nhóm 2 (C. phongnhakebangensis ở Nam Lào, Nam và Trung Việt Nam), và nhóm 3 (C. irregularis phân bố ở Trung Lào và Bắc Trung Bộ Việt Nam). Nghiên cứu của Brennan et al. (2017) và Ngô Thị Hạnh (2017) trên 226 mẫu vật Đông Dương tiếp tục mở rộng mô hình phân nhóm khi tích hợp cả chỉ thị COI và ND2, cho thấy các xung đột trong việc phân định ranh giới giữa các đơn vị phân loại tiến hóa độc lập (Evolutionarily Significant Units - ESUs) tại các khối karst chuyển tiếp.

So sánh với các nghiên cứu quốc tế trong khu vực Đông Dương và lân cận, mức độ đa dạng họ Gekkonidae tại Lào (45 loài tính đến 2018) bộc lộ sự thiếu hụt rõ rệt khi đối sánh với Việt Nam (94 loài), Thái Lan (94 loài), Myanmar (77 loài), dù cao hơn Campuchia (25 loài). Luận án của Saly Sitthivong định vị nghiên cứu vào chính tâm điểm giao thoa địa động lực học giữa mảng Đông Bắc Lào và khối karst Trường Sơn Trung Lào, trực tiếp giải quyết sự thiếu hụt mẫu vật chuẩn hóa và cung cấp dữ liệu phát sinh loài đa diện nhằm hợp nhất các luồng quan điểm phân loại học trước đó.

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

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

Luận án đóng góp quan trọng vào việc kiểm chứng và mở rộng Lý thuyết Địa lý sinh vật Đảo (Theory of Island Biogeography - MacArthur & Wilson, 1967) khi áp dụng cho các hệ sinh thái "ốc đảo đá vôi lục địa" (terrestrial karst islands). Kết quả nghiên cứu chứng minh rằng các dãy núi đá vôi bị chia cắt bởi thung lũng nông nghiệp và sông suối đóng vai trò như các hàng rào địa lý tuyệt đối ngăn chặn dòng gen (gene flow), thúc đẩy quá trình phân nhánh tiến hóa cục bộ.

Về mặt lý thuyết loài, luận án củng cố Khái niệm loài tiến hóa thống nhất (General Lineage Concept / Unified Species Concept - de Queiroz, 1998, 2007) và Khái niệm loài phát sinh chủng loại (Phylogenetic Species Concept - Cracraft, 1983) thông qua phương pháp Phân loại tích hợp (Integrative Taxonomy). Ba mệnh đề lý thuyết (Theoretical Propositions) được xác lập:

  1. Mệnh đề P1 (Karst Micro-endemism): Mỗi hệ thống khối karst cô lập về mặt thủy văn và địa chất tại Lào định hình ít nhất một dòng tiến hóa tắc kè đặc hữu hẹp với sự phân hóa di truyền vượt ngưỡng sai khác loài tiêu chuẩn (>6% phân kỳ gen COI).
  2. Mệnh đề P2 (Microhabitat Niche Partitioning): Cấu trúc phân tầng vi sinh cảnh (vách đá đứng, hang động karst ngầm, thảm thực vật bám đá) tạo ra áp lực chọn lọc phân kỳ hình thái (tỉ lệ giác bám, chiều dài chi, số lượng nốt sần), duy trì sự cùng tồn tại sinh thái của các loài đồng xứ (sympatric species).
  3. Mệnh đề P3 (Biogeographical Suture Zones): Lưu vực sông Mê Kông và đứt gãy địa chất Bắc - Trung Lào đóng vai trò là ranh giới rẽ nhánh phát sinh chủng loại cổ xưa, phân định hai nhóm tiến hóa độc lập của giống Cyrtodactylus và Gekko.

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

Khung phân tích của luận án tích hợp đa tầng giữa Sinh thái học Karst, Sinh học phân tử tiến hóa và Đánh giá không gian bảo tồn.

  ┌─────────────────────────────────────────────────────────────────────────┐
  │                    KHUNG PHÂN TÍCH TÍCH HỢP ĐA TẦNG                     │
  └─────────────────────────────────────────────────────────────────────────┘
                                       │
     ┌─────────────────────────────────┼─────────────────────────────────┐
     ▼                                 ▼                                 ▼
┌──────────────────────────┐ ┌──────────────────────────┐ ┌──────────────────────────┐
│  HÌNH THÁI HỌC SO SÁNH   │ │  DI TRUYỀN HỌC PHÂN TỬ   │ │   SINH THÁI HỌC QUẦN XÃ   │
│  - 39 chỉ tiêu đo/đếm    │ │  - Gen COI (ty thể)      │ │  - Phân tầng sinh cảnh    │
│  - Thước kẹp số 0.1 mm   │ │  - MP, ML, BI algorithms │ │  - Đai cao (mỗi 200 m)    │
│  - Kính hiển vi điện tử  │ │  - Mô hình TVM+I+G       │ │  - Chỉ số Sorensen-Dice   │
└──────────────────────────┘ └──────────────────────────┘ └──────────────────────────┘
     │                                 │                                 │
     └─────────────────────────────────┼─────────────────────────────────┘
                                       ▼
  ┌─────────────────────────────────────────────────────────────────────────┐
  │            ĐỊNH LOẠI TÍCH HỢP & MA TRẬN ƯU TIÊN BẢO TỒN (1-13 ĐIỂM)     │
  └─────────────────────────────────────────────────────────────────────────┘

Cách tiếp cận phân tích bao gồm việc chuẩn hóa 39 chỉ tiêu hình thái đo đếm kết hợp giải trình tự đoạn gen chỉ thị barcode COI, phân tích topo cây phát sinh chủng loại bằng 3 thuật toán độc lập (Tối ưu hóa hình thái - Maximum Parsimony, Khả năng xảy ra tối đa - Maximum Likelihood, và Suy luận Bayes - Bayesian Inference) dưới mô hình tiến hóa nucleotide tối ưu. Đồng thời, cấu trúc quần xã được định lượng thông qua ma trận tương đồng Sorensen-Dice $d_{jk} = 2M / (2M + N)$ trên phần mềm PAST Statistic. Điều kiện biên (boundary conditions) của khung phân tích được xác lập nghiêm ngặt: giới hạn trong bộ Squamata, họ Gekkonidae phân bố tại các sinh cảnh rừng trên núi đá vôi nhiệt đới gió mùa thuộc lãnh thổ Lào.

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 lập trường nhận thức luận Hậu thực chứng (Post-positivism), kết hợp hài hòa giữa quan sát thực nghiệm quy nạp (thu thập mẫu vật, đo đạc hình thái chi tiết) và kiểm định giả thuyết diễn dịch phát sinh loài. Thiết kế nghiên cứu đa tầng (Multi-level design) bao quát 3 cấp độ:

  1. Cấp độ cảnh quan/vùng: So sánh 6 tỉnh đại diện cho 2 miền địa lý sinh vật (Bắc Lào: U Đôm Xay, Luông Pha Bang, Húa Phăn; Trung Lào: Viêng Chăn, Xiêng Khoảng, Khăm Muôn).
  2. Cấp độ hệ sinh thái/địa điểm: 13 huyện trọng điểm với 38 tuyến điều tra thực địa dài từ 2 đến 5 km băng qua các họng hang, vách đá karst hiểm trở.
  3. Cấp độ vi sinh cảnh: Phân loại theo 3 sinh cảnh chính (SC1: Núi đá vôi khu canh tác nông nghiệp; SC2: Rừng thứ sinh trên núi đá vôi nghèo kiệt; SC3: Rừng nguyên sinh giàu trên núi đá vôi), 6 đai độ cao phân chia chi tiết theo bước 200 m (<200 m, 200-400 m, 400-600 m, 600-800 m, 800-1.000 m, >1.000 m), và 3 vị trí bám ngẫu nhiên (vách đá, trên thân cây/dây leo, mặt đất/khe nứt sát chân núi).
   ┌──────────────────────────────────────────────────────────────────────┐
   │ BỐ TRÍ KHẢO SÁT THỰC ĐỊA: 6 TỈNH - 13 ĐỊA ĐIỂM - 112 NGÀY (2018-2020) │
   └──────────────────────────────────────────────────────────────────────┘
   MIỀN BẮC LÀO:
   ├── U Đôm Xay       : Na Mo (620-827 m), Mường Xay (642-782 m)
   ├── Luông Pha Bang  : Mường Ngoi (348-702 m), Luông Pha Bang (298-439 m)
   └── Húa Phăn        : Mường Hiêm (687-1.270 m), Mường Xon (715-926 m), Viêng Xay (698-985 m)
   
   MIỀN TRUNG LÀO:
   ├── Viêng Chăn      : Văng Viêng (230-684 m), Mường Phương (232-675 m)
   ├── Xiêng Khoảng    : Phu Kut (1.397 m), Nong Hẹt
   └── Khăm Muôn       : Khun Khăm (168-524 m), Bua La Pha (170-214 m)

Quy trình nghiên cứu rigorous

Quy trình khảo sát thực địa được triển khai ban đêm (17h00 - 23h00) bằng phương pháp tìm kiếm tích cực trên các tuyến GPS định vị bằng Trackmaker. Mẫu vật sống được bắt thủ công hoặc dùng kẹp bọc cao su chống tổn thương cơ học. Mẫu sau đó được chụp ảnh sống ở các góc độ giải phẫu, gây mê nhân đạo bằng bông thấm ethylacetate trong bình kín (Simmons, 2002), thu mẫu mô cơ đùi hoặc mô gan/đuôi ngâm riêng trong cồn tuyệt đối Merck 70-96% để phân tích DNA, trước khi cố định mẫu vật tổng thể bằng cồn 90% (5-8 giờ) và lưu giữ vĩnh viễn trong cồn 70%.

Quy trình phân tích hình thái đo đếm (morphometrics & meristics) tuân thủ hệ thống tiêu chuẩn quốc tế nghiêm ngặt theo Nguyen et al. (2013, 2015), Luu et al. (2015, 2016), Ziegler et al. (2013, 2016):

  • 19 chỉ số đo kích thước (mm) bằng thước kẹp điện tử Etopoo (độ chính xác 0,01 mm): Chiều dài mút mõm - hậu môn (SVL), chiều dài đuôi (TaL), khoảng cách nách - háng (AG), chiều dài đầu (HL), chiều rộng đầu (HW), chiều cao đầu (HH), khoảng cách mõm - mắt (SE), khoảng cách mắt - tai (EyeEar), chiều dài chi trước (ForeL), chiều dài đùi (FemurL), chiều dài cẳng chân (CrusL), chiều dài ngón tay IV (LD4A), chiều dài ngón chân IV (LD4P), đường kính mắt (OD), chiều rộng tai (EarL), chiều rộng tấm mõm (RW), chiều cao tấm mõm (RH), chiều rộng tấm cằm (MW), chiều dài tấm cằm (ML).
  • 20 chỉ số đếm vảy giải phẫu dưới kính hiển vi điện tử: Số gai mi mắt (CS), số vảy quanh mũi (N), tấm gian mũi (I), vảy môi trên (SL), vảy môi dưới (IL), vảy góc trước mắt (IO), hàng vảy mõm-mắt (PO), tấm cằm (PM), vảy quanh tấm cằm (GP), hàng củ lồi lưng (DTR), số vảy quanh củ lồi (GSDT), vảy bụng từ sau cằm đến lỗ huyệt (SBL), số hàng vảy quanh giữa thân (SR), số hàng vảy bụng (V), vảy dưới ngón tay I và IV (LF1, LF4), vảy dưới ngón chân I và IV (LT1, LT4), số lỗ trước huyệt/lỗ đùi (PP), và củ lồi gốc đuôi (PAT).
   ┌─────────────────────────────────────────────────────────────┐
   │ QUY TRÌNH PHÂN TÍCH TÁCH CHIẾT DNA & XỬ LÝ PHÂN TỬ CYTOCHROME │
   └─────────────────────────────────────────────────────────────┘
   Mẫu mô (Cơ/Gan/Đuôi bảo quản -4°C)
      │
      ▼
   Tách chiết DNA (Qiagen DNeasy Kit / GeneJET Purification Kit)
      │
      ▼
   Kiểm tra điện di Gel Agarose 1%, Đệm TBE 1X, 70V (30 phút)
      │
      ▼
   Khuếch đại PCR gen COI (Thể tích 21 µl, HotStarTaq/DreamTaq Mastermix)
   - Cyrtodactylus : VF1-d (5'-TTCTCAACCAACCACAARGAYATYGG-3') / VR1-d
   - Dixonius      : MetF1 (5'-AAGCTTTCGGGCCCATACC-3') / COIR1
   - Gekko         : GF1 (5'-CAAGCACHATYATYACYATAT-3') / GR1
      │
      ▼
   Kiểm định PCR & Tinh sạch GeneJET PCR Purification Kit
      │
      ▼
   Giải trình tự Sanger 2 chiều (Công ty FirstBase, Malaysia)
      │
      ▼
   Phân tích phát sinh loài (Sequencher 5.1 -> Clustal X v2 -> jModelTest 2 -> PAUP*/IQ-TREE/MrBayes)

Quy trình sinh học phân tử được chuẩn hóa cao độ: Tách chiết DNA tổng số bằng bộ kit DNeasy Blood and Tissue (Qiagen, CHLB Đức) cho mẫu dung lượng thấp và GeneJet Genomic DNA Purification (ThermoFisher Scientific, Lithuania) cho mẫu tươi; kiểm tra nồng độ trên gel agarose 1% đệm TBE 1X ở 70V trong 30 phút, đối chứng kích thước bằng thang chuẩn Marker 1 kb và 100 bp (1st Base, Malaysia), hiện ảnh huỳnh quang UV trên hệ thống Alphamager MINI (Protein Simple, Mỹ). Phản ứng PCR khuếch đại phân đoạn gen ty thể Cytochrome c Oxidase Subunit I (COI) được thực hiện trong tổng thể tích 21 µl (gồm 1-2 µl DNA khuôn, 2 µl mỗi mồi 10 µM, 5 µl $H_2O$, 10 µl HotStarTaq Mastermix hoặc DreamTaq Mastermix). Các cặp mồi chuyên biệt bao gồm:

  • 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). Chu trình nhiệt PCR: Biến tính ban đầu ở 95°C (5-15 phút); 35 chu kỳ gồm 95°C (30 giây), bắt cặp mồi ở 48°C - 60°C (45 giây), kéo dài ở 72°C (1 phút); và bước kéo dài cuối cùng ở 72°C (6 phút). Mẫu tinh sạch bằng GeneJET PCR Purification được giải trình tự hai chiều tại FirstBase (Malaysia).

Data và phân tích

Dữ liệu trình tự thô được hiệu chỉnh và lắp ráp trên Sequencher v5.1 (Gene Codes Corporation), kiểm tra độ tương đồng thông qua BLAST trên NCBI GenBank, và căn chỉnh thẳng hàng đa trình tự bằng Clustal X v2 (Thompson et al., 1997). Lựa chọn mô hình tiến hóa tốt nhất dựa trên tiêu chí AIC thông qua jModelTest v2.1.4 (Darriba et al., 2012), xác định TVM+I+G là mô hình thay thế nucleotide tối ưu cho tập dữ liệu COI.

Phân tích phát sinh chủng loại được xây dựng đồng thời trên 3 nền tảng:

  • Maximum Parsimony (MP): Chạy trên PAUP* 4.0b10 (Swofford, 2001) với 1.000 bản sao bootstrap tái lấy mẫu (Bootstrap Percentages - BP).
  • Maximum Likelihood (ML): Chạy trên IQ-TREE v1.1.1 (Nguyen et al., 2015) với 10.000 bản sao bootstrap siêu nhanh (Ultrafast Bootstrap - UFBP).
  • Bayesian Inference (BI): Chạy trên MrBayes v3.2.7 (Ronquist et al., 2012) với $1 \times 10^7$ thế hệ, lấy mẫu mỗi 1.000 thế hệ qua 4 chuỗi Markov Chain Monte Carlo (MCMC). Điểm hội tụ log-likelihood (-lnL) đạt trạng thái dừng (burn-in) sau 54.000 thế hệ (cut-off 54). Độ tin cậy của các nút nhánh được công nhận là có độ hỗ trợ mạnh mẽ khi BP $\ge 70%$ và Xác suất hậu nghiệm Bayesian (Posterior Probability - PP), UFBP $\ge 95%$ (Hillis & Bull, 1993; Ronquist et al., 2012). Khoảng cách di truyền chưa hiệu chỉnh theo cặp (uncorrected pairwise p-distances) được tính toán trực tiếp trong PAUP*.

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

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

   ┌─────────────────────────────────────────────────────────────┐
   │ TỔNG KẾT THÀNH PHẦN LOÀI VÀ CÁC PHÁT HIỆN ĐỘT PHÁ CỦA ĐỀ TÀI │
   └─────────────────────────────────────────────────────────────┘
   Tổng số ghi nhận: 28 loài thuộc 6 giống họ Gekkonidae
   ├── Cyrtodactylus      : 9 loài (3 loài mới cho khoa học, 1 ghi nhận mới)
   ├── Dixonius           : 4 loài (2 loài mới cho khoa học)
   ├── Gekko              : 4 loài (1 loài mới cho khoa học)
   ├── Hemidactylus       : 4 loài (Các loài phổ biến diện rộng)
   ├── Hemiphyllodactylus : 4 loài (1 ghi nhận mới cho tỉnh)
   └── Gehyra             : 3 loài (Nhiều dòng ẩn danh tiềm năng)

Nghiên cứu mang lại 5 nhóm phát hiện khoa học đột phá:

1. Khám phá 06 loài mới cho khoa học: Luận án đã mô tả chính xác 6 đơn vị phân loại mới với đầy đủ chẩn loại hình thái và bằng chứng phân tử:

  • Cyrtodactylus houaphanensis Sitthivong et al. (thu tại vùng karst Mường Xon, tỉnh Húa Phăn).
  • Cyrtodactylus muangfuangensis Sitthivong et al. (thu tại huyện Mường Phương, tỉnh Viêng Chăn).
  • Cyrtodactylus ngoiensis Sitthivong et al. (thu tại huyện Mường Ngoi, tỉnh Luông Pha Bang).
  • Dixonius lao Sitthivong et al. (loài thằn lằn chân lá mới ghi nhận tại miền Trung Lào).
  • Dixonius somchanhae Sitthivong et al. (loài thằn lằn chân lá mới tại các vách đá vôi Viêng Chăn).
  • Gekko khunkhamensis Sitthivong et al. (thu tại huyện Khun Khăm, tỉnh Khăm Muôn).

2. Ghi nhận phân bố mới và nhận diện các dòng ẩn danh (Cryptic taxa): Bổ sung 02 ghi nhận phân bố mới cho tỉnh U Đôm Xay là Cyrtodactylus wayakonei và Hemiphyllodactylus kiziriani. Đặc biệt, luận án đã định dạng hình thái và phả hệ của 7 dòng tắc kè chưa thể định danh chính xác đến loài gồm: Cyrtodactylus sp. 1, Cyrtodactylus sp. 2, Gehyra sp. 1, Gehyra sp. 2, Dixonius sp., Hemiphyllodactylus sp. 1, và Hemiphyllodactylus sp. 2, mở ra các nghiên cứu mô tả loài mới tiếp theo.

3. Tái cấu trúc quan hệ phát sinh chủng loại 3 giống chủ đạo tại Lào: Cây phát sinh loài gen COI củng cố mạnh mẽ các phân nhóm tiến hóa (PP = 1.0, UFBP = 100%):

  • Giống Cyrtodactylus tại Lào phân lập thành các nhánh tương ứng với địa lý sinh vật: Nhóm C. wayakonei (tập trung ở vùng núi Tây Bắc và Bắc Lào gồm U Đôm Xay, Luông Pha Bang, Húa Phăn); Nhóm C. phongnhakebangensis (chiếm ưu thế ở khối karst Trường Sơn từ Khăm Muôn đến Viêng Chăn).
  • Giống Gekko khẳng định vị trí phân loại của loài mới G. khunkhamensis thuộc nhóm G. japonicus phức tạp, có quan hệ họ hàng gần với G. scientiadventura và G. bonkowskii.
  • Giống Dixonius ghi nhận sự phân nhánh bất ngờ khi D. lao và D. somchanhae hình thành một clade độc lập tách biệt khỏi quần thể Dixonius siamensis truyền thống tại Thái Lan và Campuchia.

4. Phân hóa sinh cảnh và ổ sinh thái đai cao:

  • Sinh cảnh: Rừng nguyên sinh trên núi đá vôi (SC3) duy trì độ giàu thành phần loài cao nhất (chiếm hơn 75% tổng số cá thể bắt gặp), trong khi sinh cảnh nông nghiệp (SC1) chỉ ghi nhận các loài thích nghi rộng như Hemidactylus frenatus, Hemidactylus platyurus.
  • Đai độ cao: Mật độ cá thể và số lượng loài tập trung đậm đặc nhất ở đai cao từ 200 m đến 600 m (chiếm >60% mẫu vật), giảm mạnh ở đai cao trên 1.000 m (chỉ bắt gặp một số loài chịu lạnh thích nghi hẹp tại Húa Phăn và Xiêng Khoảng).
  • Nơi ở: 58% mẫu vật được ghi nhận trực tiếp trên vách đá thẳng đứng và hốc hang karst, 28% trên thân cây bám vách đá, và chỉ 14% bắt gặp trên mặt đất hoặc vụn đá chân núi.

5. Mức độ tương đồng thành phần loài suy giảm mạnh giữa hai miền: Phân tích chỉ số Sorensen-Dice cho thấy độ tương đồng thành phần loài giữa các tỉnh miền Bắc (Húa Phăn, Luông Pha Bang, U Đôm Xay) và miền Trung (Viêng Chăn, Khăm Muôn) đạt mức rất thấp ($d_{jk} < 0,35$). Tỉnh Viêng Chăn có tính đa dạng giống loài cao nhất (14 loài thuộc 5 giống), tiếp đến là Khăm Muôn, phản ánh vai trò của vùng chuyển tiếp sinh thái - địa chất.

Implications đa chiều

  • Về mặt học thuật: Cung cấp cơ sở dữ liệu giải phẫu học và barcode di truyền chuẩn mực, giải quyết triệt để các hoài nghi về tình trạng đồng danh (synonymy) và ranh giới loài của khu hệ Squamata Đông Dương.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình tích hợp giữa khảo sát định vị GPS ban đêm, ma trận 39 biến thái học và phân tích phát sinh loài Bayes nhiều dòng gene, áp dụng hiệu quả cho các vùng karst nhiệt đới khó tiếp cận.
  • Về thực tiễn và chính sách: Xây dựng ma trận 5 tiêu chí định lượng (Độ giàu loài, Số loài nguy cấp/đặc hữu, Chất lượng sinh cảnh, Mức độ tác động nhân sinh, Tình trạng bảo vệ hiện tại) với thang điểm 1-13 để xếp hạng ưu tiên bảo tồn cho 13 địa bàn nghiên cứu. Kết quả chỉ rõ: Huyện Khun Khăm (Khăm Muôn), Mường Ngoi (Luông Pha Bang), Mường Phương (Viêng Chăn) và Mường Xon (Húa Phăn) là những "điểm nóng sinh học" (conservation hotspots) cần lập tức thành lập các Khu bảo tồn loài và sinh cảnh đá vôi.
  • Khuyến nghị chính sách thực thi: Đưa 6 loài mới phát hiện vào Danh lục Đỏ IUCN ở thứ hạng Cực kỳ nguy cấp (CR) hoặc Nguy cấp (EN) do phân bố hẹp; đề xuất sửa đổi Luật Bảo vệ Động vật Hoang dã Lào số 07/QH (2007) và Nghị định Chính phủ số 08/CP/2021 nhằm siết chặt quy hoạch các mỏ khai thác đá vôi xi măng, kiểm soát nạn đốt nương làm rẫy và triệt phá đường dây săn bắt tắc kè làm dược liệu/ngâm rượu xuất khẩu.

Limitations và Future Research

Luận án thừa nhận một số hạn chế khách quan:

  1. Chỉ thị phân tử: Nghiên cứu chủ yếu khai thác phân đoạn gen ty thể đơn locus (COI). Mặc dù COI cung cấp độ phân giải loài xuất sắc cho Barcoding, hiện tượng thâm nhập dòng mẹ (mitochondrial introgression) hoặc phân hóa dòng dõi không hoàn toàn (incomplete lineage sorting) có thể gây sai lệch nhẹ ở các nút nhánh sâu.
  2. Độ phủ địa lý: Dù đã bao quát 6 tỉnh, nhiều khối núi đá vôi hiểm trở tại vùng sâu biên giới Lào - Việt (như dải karst giáp Vườn Quốc gia Nam Ét-Phou Louey) chưa thể tiếp cận trọn vẹn do địa hình dốc đứng và hạn chế về an toàn bom mìn sót lại sau chiến tranh.
  3. Mẫu vật một số loài hiếm: Một số taxon chưa xác định định danh (Cyrtodactylus sp., Gehyra sp.) chỉ thu được số lượng cá thể đơn lẻ ($N = 1-2$), chưa đủ đại diện để phân tích biến dị nội loài hoàn chỉnh.

Chương trình nghiên cứu tương lai (Future Research Agenda) cần tập trung vào 4 định hướng:

  • Ứng dụng giải trình tự thế hệ mới (NGS) đa locus nhân (nuclear genes) như UCEs (Ultraconserved Elements) hoặc ddRADseq để làm sáng tỏ các mối quan hệ phát sinh chủng loại ở mức độ dưới giống.
  • Thiết lập mô hình hóa ổ sinh thái (Ecological Niche Modeling - MaxEnt) dưới tác động của kịch bản biến đổi khí hậu toàn cầu đối với các loài tắc kè biến nhiệt phụ thuộc môi trường vi khí hậu hang động.
  • Nghiên cứu sinh học quần thể, tập tính sinh sản và cấu trúc tiếng kêu giao tiếp (Bioacoustics) của các loài tắc kè karst.
  • Mở rộng phạm vi điều tra sang các tỉnh Nam Lào (Attapeu, Champasak, Sekong) để hoàn thiện bản đồ phân bố Gekkonidae toàn quốc.

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

  • Ảnh hưởng học thuật: Các công bố quốc tế bắt nguồn từ luận án trên các tạp chí chuyên ngành phân loại hàng đầu (như Zootaxa, ZooKeys, Amphibian & Reptile Conservation) tạo nên nguồn trích dẫn học thuật quốc tế ổn định, nâng cao vị thế của khoa học lâm nghiệp và sinh thái học Lào trên bản đồ nghiên cứu Đông Dương.
  • Tác động đến bảo tồn thiên nhiên: Cung cấp luận cứ khoa học trực tiếp cho Bộ Nông Lâm nghiệp Lào, Cục Kiểm lâm và các tổ chức quốc tế (IUCN, WWF-Laos, WCS) trong việc điều chỉnh ranh giới Vườn Quốc gia Hin Nam No, VQG Nam Ét-Phou Louey và thành lập các khu bảo tồn vi mô (micro-reserves) trên karst.
  • Đóng góp kinh tế - xã hội: Cảnh báo nguy cơ suy giảm sinh cảnh do khai khoáng đá vôi và nương rẫy, giúp địa phương định hình mô hình du lịch sinh thái khám phá hang động bền vững tại Văng Viêng và Luông Pha Bang, đem lại sinh kế thay thế cho người dân bản địa thay vì khai thác động vật hoang dã.

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

  • Nghiên cứu sinh & Nhà nghiên cứu Bò sát: Tiếp cận bộ dữ liệu hình thái học đo đếm chuẩn hóa (39 biến số), hệ thống primer PCR tối ưu và dữ liệu GenBank làm mẫu chuẩn đối sánh.
  • Giảng viên & Viện nghiên cứu: Giáo trình và tài liệu tham khảo chuyên sâu cho các chương trình đào tạo sau đại học về Quản lý tài nguyên rừng, Đa dạng sinh học và Phân loại học tiến hóa tại Lào và Việt Nam.
  • Cơ quan quản lý và hoạch định chính sách: Cung cấp ma trận điểm ưu tiên bảo tồn (1-13 điểm) làm công cụ ra quyết định phê duyệt các dự án đầu tư hạ tầng giao thông và nhà máy xi măng trên núi đá vôi.
  • Tổ chức phi chính phủ (NGOs) bảo tồn: Dữ liệu nền để xây dựng các đề án tài trợ bảo tồn quốc tế (như CEPF, Rufford Foundation, Rainforest Trust) nhằm cứu trợ các loài nguy cấp có vùng phân bố cực hẹ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? Trả lời: Luận án đã mở rộng Lý thuyết Địa lý sinh vật Đảo (MacArthur & Wilson, 1967) sang hệ thống "đảo đá vôi lục địa" (continental karst islands), chứng minh tính phân mảnh sinh cảnh karst thúc đẩy sự hình thành loài nhanh chóng (allopatric speciation) và tạo ra độ đặc hữu hẹp cực cao ở bò sát biến nhiệt.

2. Đổi mới phương pháp luận nổi bật nhất so với các nghiên cứu trước đây là gì? Trả lời: So với các nghiên cứu của Nazarov et al. (2014) và Luu et al. (2016) vốn chỉ phân tích đơn lẻ hoặc hạn chế số lượng mẫu, luận án đã thiết lập quy trình Phân loại tích hợp toàn diện: kết hợp ma trận 39 chỉ tiêu đo đếm hình thái học vi thể với 3 thuật toán phát sinh loài (MP, ML, BI) trên mô hình TVM+I+G, đồng thời phân tích định lượng tương đồng quần xã Sorensen-Dice trên cùng một tập dữ liệu chuẩn hóa 138 mẫu vật.

3. Phát hiện bất ngờ nhất từ tập dữ liệu thực nghiệm là gì? Trả lời: Phát hiện 6 loài mới hoàn toàn cho khoa học và ghi nhận sự tồn tại của 7 dòng ẩn danh chưa thể định danh, khẳng định mức độ đa dạng loài tắc kè tại các khối karst miền Bắc và miền Trung Lào bị đánh giá thấp nghiêm trọng trong các thập kỷ trước.

4. Luận án có cung cấp giao thức nhân bản (Replication protocol) hoàn chỉnh không? Trả lời: Có. Luận án cung cấp chi tiết tuyệt đối: tọa độ GPS và độ cao của 13 điểm nghiên cứu, chu trình nhiệt và trình tự các cặp mồi PCR chuyên biệt (VF1-d/VR1-d, MetF1/COIR1, GF1/GR1), bộ kit tách chiết chuẩn, thông số cài đặt MrBayes ($1 \times 10^7$ thế hệ, burn-in 54) và bảng định nghĩa 39 chỉ số đo đếm hình thái.

5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào? Trả lời: Tập trung vào giải mã hệ gen quy mô lớn (Phylogenomics/RADseq), đánh giá khả năng chống chịu nhiệt độ hang động dưới áp lực biến đổi khí hậu, thiết lập mạng lưới khu bảo tồn vi mô trên karst và hoàn thiện giám định di truyền toàn bộ các loài tắc kè thương mại tại Lào.

Kết luận

  1. Xác lập bộ cơ sở dữ liệu toàn diện đầu tiên về họ Tắc kè (Gekkonidae) tại các hệ sinh thái núi đá vôi miền Bắc và miền Trung Lào với 28 loài thuộc 6 giống được mô tả giải phẫu chi tiết.
  2. Đóng góp đột phá cho khoa học phân loại thế giới thông qua việc công bố 06 loài mới (Cyrtodactylus houaphanensis, C. muangfuangensis, C. ngoiensis, Dixonius lao, D. somchanhae, Gekko khunkhamensis) và bổ sung 02 ghi nhận phân bố mới (C. wayakonei, H. kiziriani).
  3. Giải mã cấu trúc quan hệ phát sinh chủng loại di truyền gen COI của 3 giống chủ đạo (Cyrtodactylus, Dixonius, Gekko), làm sáng tỏ các phân nhóm tiến hóa và chứng minh sự phân kỳ sâu sắc giữa các quần thể karst cô lập.
  4. Định lượng hóa đặc điểm phân bố sinh thái, chỉ ra rằng sinh cảnh rừng nguyên sinh trên núi đá vôi (SC3) và đai độ cao 200–600 m là nơi cư trú trọng yếu nhất của các loài tắc kè đặc hữu.
  5. Xây dựng thành công thang đánh giá và xác định các điểm nóng ưu tiên bảo tồn (Khun Khăm, Mường Ngoi, Mường Phương, Mường Xon), cung cấp cơ sở khoa học cốt lõi cho chiến lược quản lý bền vững tài nguyên rừng và đa dạng sinh học tại CHDCND Lào.