Tổng quan luận án

Luận án tiến sĩ sinh học của tác giả Nguyễn Thị Định (chuyên ngành Sinh thái học, mã số: 9 42 01 20) nghiên cứu cấu trúc quần xã bọ cánh cứng ăn lá và mối quan hệ dinh dưỡng của chúng với thực vật tại Vườn quốc gia (VQG) Núi Chúa, tỉnh Ninh Thuận. Đề tài được thực hiện dưới sự hướng dẫn khoa học của PGS. Nguyễn Văn Sinh và TS. Jesús Gómez-Zurita, hoàn thành năm 2018 tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, phối hợp cùng Viện Nghiên cứu Sinh học Tiến hóa (Institut de Biologia Evolutiva - IBE, Tây Ban Nha).

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

Trong các hệ sinh thái trên cạn, các quần xã sinh vật liên kết và phụ thuộc lẫn nhau; trong đó, quần xã thực vật chịu sự chi phối mạnh mẽ bởi các yếu tố vật lý của môi trường (nhiệt độ, độ ẩm, bức xạ mặt trời) thay đổi theo độ cao địa hình. Quần xã thực vật là nguồn cung cấp thức ăn trực tiếp, quyết định đến tính đa dạng và sự phân bố không gian của các nhóm côn trùng ăn thực vật. Họ Bọ cánh cứng ăn lá (Chrysomelidae, thuộc bộ Cánh cứng Coleoptera) là một trong những họ côn trùng phong phú nhất với tập tính ăn thực vật chuyên hóa cao và liên kết chặt chẽ với cây chủ trong suốt vòng đời.

VQG Núi Chúa sở hữu hệ sinh thái rừng bán khô hạn đặc thù với sự phân hóa môi trường liên tục từ rừng khô hạn vùng đất thấp, qua đới chuyển tiếp bán ẩm, lên đến rừng ẩm thường xanh trên núi cao. Mặc dù vậy, các nghiên cứu trước đây tại Việt Nam chủ yếu tiếp cận theo phương pháp hình thái truyền thống, ghi nhận thực vật chủ bằng quan sát thực địa vốn chứa đựng nhiều sai số, đồng thời chưa từng có nghiên cứu nào ứng dụng công cụ sinh học phân tử để đánh giá đồng thời đa dạng loài và quan hệ dinh dưỡng của Chrysomelidae dọc theo dải biến đổi môi trường tại khu vực này.

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

Luận án xác định 02 mục tiêu nghiên cứu cụ thể:

  1. Đánh giá sự đa dạng loài và sự biến động theo không gian của họ Chrysomelidae ở VQG Núi Chúa.
  2. Xác định thức ăn của các loài họ Chrysomelidae và đánh giá sự thay đổi thức ăn của chúng theo không gian.

Để đạt được các mục tiêu trên, luận án triển khai 03 nội dung nghiên cứu chính:

  1. Sử dụng bộ ba mã vạch ADN để đánh giá sự đa dạng loài của họ Chrysomelidae ở VQG Núi Chúa.
  2. Sử dụng bộ ba mã vạch ADN để xác định thức ăn của các loài họ Chrysomelidae ở VQG Núi Chúa.
  3. Xác định sự biến đổi thành phần của họ Chrysomelidae và thức ăn của chúng theo không gian ở VQG Núi Chúa.

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

  • Đối tượng nghiên cứu: Các loài côn trùng thuộc họ Bọ cánh cứng ăn lá (Chrysomelidae, Coleoptera) và nguồn thức ăn thực vật chủ trong ruột của chúng.
  • Phạm vi không gian: Khu vực nghiên cứu có diện tích khoảng 10 km² thuộc VQG Núi Chúa (nằm trên địa bàn 02 huyện Ninh Hải và Thuận Bắc, tỉnh Ninh Thuận), trải dài từ dải rừng khô ven biển ở độ cao sát mực nước biển lên đến các sinh cảnh rừng ẩm ở độ cao trên 300 m (điểm thu mẫu cao nhất đạt 486 m trên các tuyến chính và 714 m tại tuyến khảo sát).
  • Phạm vi thời gian: Mẫu vật được thu thập qua các đợt điều tra thực địa từ tháng 5/2012 đến tháng 4/2013 (khảo sát tháng 5/2012; thu mẫu định kỳ 1 lần/tháng vào các tháng 6, 7, 9/2012 và tháng 1/2013; thu mẫu 2 lần/tháng từ tháng 2 đến tháng 4/2013).

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

Các hướng nghiên cứu trên thế giới

Họ Chrysomelidae có quy mô phân loại rất lớn với hơn 37.000 loài đã được mô tả thuộc hơn 2.000 giống (Jolivet, Petitpierre & Hsiao 1988; Bouchard et al. 2011). Mối quan hệ giữa bọ cánh cứng ăn lá và thực vật chủ đã được tổng kết sâu rộng trong các công trình kinh điển của Jolivet & Hawkeswood (1995), Wagner (1999), Erwin (1982), Novotny & Miller (2014). Tại các vùng nhiệt đới, độ phong phú loài của họ này rất cao, được ghi nhận qua các nghiên cứu của Mohamedsaid (2004) tại Malaysia (1.073 loài), Thormann et al. (2016) tại Ecuador (266 loài), hay Sánchez-Reyes et al. (2014, 2016) tại Mexico (157 loài).

Về mặt phương pháp luận, việc ứng dụng mã vạch ADN (DNA barcoding) dựa trên đoạn gen ty thể Cytochrome c oxidase subunit I (cox1 / COI) dài khoảng 600–800 bp đã trở thành tiêu chuẩn quốc tế trong phân định và kiểm kê nhanh đa dạng loài động vật (Hebert et al. 2003, 2004, 2010; Savolainen et al. 2005; Kress & Erickson 2008). Trong nghiên cứu tương tác sinh thái, thay vì các phương pháp quan sát thực địa, nuôi thử nghiệm (Barone 1998) hoặc phân tích đồng vị bền (Otte & Joern 1976; Joern & Parker 1978) vốn gặp nhiều hạn chế ở địa hình đồi núi hiểm trở, các nhà khoa học đã phát triển kỹ thuật giải mã trình tự ADN lục lạp (psbA-trnH, rbcLa, matK) từ mô thực vật bị tiêu hóa trong ruột côn trùng để định danh thức ăn (Jurado-Rivera et al. 2009; Kishimoto-Yamada et al. 2013; García-Robledo et al. 2013; Liu J et al. 2015; Papadopoulou et al. 2013, 2015).

Về mặt sinh thái học quần xã và độ cao, các nghiên cứu của Whittaker (1972), Baselga (2010), Colwell & Lees (2000), Haller et al. (2013), Pausas & Austin (2001), Jill et al. (2009), Yang et al. (2015), Arnan et al. (2006) và Fahrig (2003) đã thiết lập các khung lý thuyết về phân rã đa dạng sinh học (alpha, beta, gamma), hiệu ứng trung tâm lãnh thổ (mid-domain effect) và phản ứng của sinh vật với sự phân mảnh, biến đổi sinh cảnh. Khung phân tích mạng lưới dinh dưỡng nhị phân (bipartite network) và phân rã đa dạng beta tương tác được chuẩn hóa bởi Lewinsohn et al. (2006), Dormann et al. (2008, 2009) và Poisot et al. (2012).

Tình hình nghiên cứu tại Việt Nam

Tại Việt Nam, các nghiên cứu về họ Chrysomelidae từ năm 1975 đến 2015 tập trung chủ yếu vào phân loại hình thái học và kiểm kê thành phần loài ở miền Bắc, miền Trung và Tây Nguyên (Đặng Thị Đáp & Medvedev 1982, 1983, 1989; Trần Thiếu Dư và cs 2005, 2006, 2007, 2008; Medvedev 1985–2015; Delobel et al. 2011; Kimoto 1981, 1982, 1989, 2000). Toàn quốc đã ghi nhận khoảng 700 loài và ước tính có thể vượt 1.000 loài. Về thực vật chủ, Medvedev & Đặng Thị Đáp (1982, 1983) đã lập danh mục thức ăn cho 212 loài dựa trên quan sát bắt gặp tại thực địa. Về khu hệ côn trùng tại VQG Núi Chúa, công trình điều tra năm 2005 của Tạ Huy Thịnh và cộng sự ghi nhận 10 bộ, 95 họ Côn trùng nhưng chưa khảo sát họ Chrysomelidae.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

  1. Thiếu hụt hoàn toàn dữ liệu về thành phần loài và cấu trúc quần xã bọ cánh cứng ăn lá Chrysomelidae tại VQG Núi Chúa.
  2. Chưa có nghiên cứu nào tại Việt Nam áp dụng công cụ mã vạch ADN ty thể (cox1) để phân định loài Chrysomelidae và mã vạch ADN lục lạp (psbA-trnH) để xác định chính xác thực vật chủ từ thức ăn trong ống tiêu hóa của côn trùng.
  3. Chưa có công trình phân tích cấu trúc mạng lưới tương tác thực vật - côn trùng ăn lá và sự biến đổi đa dạng beta của tương tác dọc theo dải biến đổi môi trường liên tục (gradient độ cao dưới 1.000 m) ở vùng rừng khô hạn nhiệt đới.

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 đa dạng sinh học: Sử dụng các khái niệm đa dạng alpha, đa dạng beta và đa dạng gamma theo Whittaker (1972).
  • Khung phân rã đa dạng beta: Tiếp cận theo mô hình của Baselga (2010), Cardoso et al. (2014) và Carvalho (2012), chia đa dạng beta tổng ($\beta_{total}$) thành hai thành phần độc lập: thành phần thay thế loài ($\beta_{repl}$ hoặc $\beta_{sim}$) và thành phần khác biệt độ giàu loài / tạo ổ ($\beta_{rich}$ hoặc $\beta_{sne}$).
  • Mô hình hiệu ứng trung tâm lãnh thổ (Mid-domain effect): Đánh giá sự tập trung độ giàu loài ở vùng giữa dải độ cao (Colwell & Lees 2000).
  • Mô hình mạng lưới tương tác sinh thái: Phân tích mạng lưới hai phía (bipartite network) giữa thực vật chủ và bọ cánh cứng ăn lá (Lewinsohn et al. 2006; Dormann et al. 2008). Đa dạng beta của tương tác sinh thái được phân rã theo Poisot et al. (2012): $\beta_{WN} = \beta_{OS} + \beta_{ST}$ (trong đó $\beta_{WN}$ là sự khác biệt tương tác tổng thể, $\beta_{OS}$ là khác biệt tương tác giữa các loài chia sẻ chung, và $\beta_{ST}$ là khác biệt tương tác do sự thay thế loài).

Phương pháp thu thập mẫu vật và phân chia sinh cảnh

  • Thiết kế tuyến: Khảo sát 05 tuyến đường rừng cố định lan tỏa theo các hướng tại VQG Núi Chúa: Mái Nhà (MN), Đá Đỏ (DD), Ao Hồ (AH), Núi Ông (NO), Suối Trục (ST) và 01 tuyến khảo sát ban đầu là Đá Hang (DH).
  • Phương pháp thu mẫu: Chuẩn hóa kỹ thuật đập cành/lá (beating method) trong 10 phút tại mỗi điểm cố định dọc theo đai cao, thu mẫu từ thảm cỏ, cây bụi đến tầng cành cao tối đa 2,5 m. Mẫu vật thu được ngâm ngay vào cồn tuyệt đối 96% để bảo quản ADN. Xác mẫu vật sau khi tách chiết ADN được giữ lại làm khô và gắn tiêu bản đối chứng.
  • Phân định sinh cảnh: Kết hợp dữ liệu nhiệt độ, thảm thực vật và phương pháp cửa sổ trượt (sliding window, Barton et al.) theo các khoảng cao 40 m liền kề bằng chỉ số SØrensen-Dice để xác định ranh giới sinh thái.
Tuyến thu mẫu Ký hiệu Tọa độ địa lý (Vĩ độ N - Kinh độ E) Dải độ cao (m) Số điểm thu mẫu
Đá Hang DH $11,6956389^\circ - 11,7221667^\circ \text{N}; 109,162736^\circ - 109,136111^\circ \text{E}$ 109 – 714 5
Ao Hồ AH $11,7279722^\circ - 11,7352222^\circ \text{N}; 109,202944^\circ - 109,219500^\circ \text{E}$ 50 – 265 10
Núi Ông NO $11,7334444^\circ - 11,7456944^\circ \text{N}; 109,184389^\circ - 109,175028^\circ \text{E}$ 122 – 486 14
Mái Nhà MN $11,7279830^\circ - 11,7192500^\circ \text{N}; 109,189704^\circ - 109,170222^\circ \text{E}$ 38 – 364 15
Đá Đỏ DD $11,7226667^\circ - 11,7164167^\circ \text{N}; 109,181306^\circ - 109,177278^\circ \text{E}$ 111 – 243 11
Suối Trục ST $11,7346667^\circ - 11,7534444^\circ \text{N}; 109,190583^\circ - 109,192875^\circ \text{E}$ 101 – 405 10

Phương pháp sinh học phân tử và định loại

  • Tách chiết ADN: Sử dụng bộ kit DNeasy Blood and Tissue (Qiagen Iberia, Madrid, Tây Ban Nha) trên toàn bộ phần mô mềm của mẫu vật.
  • Khuếch đại và giải trình tự gen cox1 (định danh côn trùng): Khuếch đại đoạn gen ty thể cox1 dài 830 bp bằng cặp mồi TL-N-3014 (Simon et al. 1994) và C1-J-2183 cải biên (Gómez-Zurita et al.). Chu trình nhiệt PCR: 94°C biến tính (30s), 50°C gắn mồi (30s), 72°C kéo dài (60s). Sản phẩm PCR được làm sạch bằng ammonium acetate và isopropanol lạnh, điện di gel agarose 1,5%, giải trình tự 2 chiều bằng kit BigDye Terminator v3. Trình tự được chỉnh sửa trên phần mềm Geneious Pro 5.6 và nộp lên ngân hàng gen EMBL-EBI (mã truy cập LT160095 - LT160588).
  • Phân định loài phân tử: Đánh giá mô hình tiến hóa nucleotide trên jModelTest 2.1.5 (lựa chọn mô hình GTR + I + G). Cây phát sinh loài ML được dựng bằng RAxML 7.2.8 và làm nhẵn thành cây ultrametric bằng PATHd8 hoặc r8s 1.71. Cây Bayesian Inference (BI) được xây dựng bằng BEAST 1.8.2 (chạy 50 triệu thế hệ MCMC). Thực hiện phân định loài bằng mô hình Generalized Mixed Yule-Coalescent (GMYC đơn ngưỡng và đa ngưỡng qua gói splits trong R 3.0.2) và mô hình Poisson Tree Processes (bPTP server).
  • Phân loại hình thái đối chứng: Soi kính lúp nổi phân chia phân họ theo Paul (2009), định loại đến chi/loài theo khóa của Kimoto (1981, 1982, 1989, 2000), kết hợp mổ kiểm tra cơ quan sinh dục khi có biến dị màu sắc.
  • Xác định thực vật chủ (ADN lục lạp): Khuếch đại vùng gen lục lạp psbA-trnH trong ruột bọ cánh cứng bằng cặp mồi PsbA và TrnH (hoặc mồi nội bộ PsbAinT2 và TrnHinT1 cho đoạn ngắn). Giải trình tự và định danh họ thực vật bằng phần mềm BAGpipe (Papadopoulou et al. 2015) đối chiếu dữ liệu GenBank.

Phương pháp phân tích sinh thái học và xử lý số liệu

  • Độ giàu loài tiềm năng: Xây dựng đường cong tích lũy loài trên phần mềm EstimateS 9.1 (Colwell 2013; Hortal et al. 2006).
  • Phân tích biến động không gian và đa dạng beta: Sử dụng các gói betapart, BAT 1.1 và betalink 2.0 trong môi trường R để tính toán các chỉ số Jaccard, SØrensen, $\beta_{sim}$, $\beta_{sne}$, $\beta_{repl}$, $\beta_{rich}$, $\beta_{WN}$, $\beta_{OS}$, $\beta_{ST}$.
  • Phân tích sắp xếp hợp quy chuẩn (CCA): Sử dụng gói vegan 2.0 kiểm định tác động của 4 biến số giải thích: Sinh cảnh (khô/ẩm theo mốc 300 m), Loài bọ cánh cứng (dữ liệu ADN), Tuyến thu mẫu, và Thời gian thu mẫu (trước/sau mùa mưa).
  • Phân tích mạng lưới tương tác: Sử dụng gói Bipartite 2.05 trong R (các hàm visweb, plotweb, networklevel) để phân tích mạng lưới tương tác Galerucinae - thực vật chủ ở hai sinh cảnh khô và ẩm.

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

Chương 1: Tổng quan tài liệu

Chương này tổng kết toàn diện các cơ sở lý luận về đa dạng sinh học và các vấn đề liên quan; phân tích sâu về cuộc khủng hoảng đa dạng sinh học, tốc độ tuyệt chủng và "trở ngại phân loại" (taxonomic impediment) khiến hàng triệu loài chưa được mô tả trước khi biến mất. Tác giả luận chứng sự cần thiết của công cụ sinh học phân tử, đặc biệt là mã vạch ADN ty thể (cox1) ở động vật và ADN lục lạp (psbA-trnH, rbcLa, matK) ở thực vật, giúp đẩy nhanh tiến độ kiểm kê loài và giải mã tương tác dinh dưỡng trong ruột côn trùng mà không phụ thuộc vào việc quan sát trực tiếp ngoài tự nhiên.

Chương 1 cũng phân tích các quy luật sinh thái về sự không đồng nhất môi trường dọc theo gradient độ cao và ảnh hưởng của nó đến phân bố loài; chứng minh họ Bọ cánh cứng ăn lá (Chrysomelidae) với hơn 37.000 loài trên thế giới là đối tượng chỉ thị phù hợp nhất để nghiên cứu tính phụ thuộc lẫn nhau trong hệ sinh thái. Phần tổng quan lịch sử nghiên cứu chỉ rõ các công trình tại Việt Nam từ trước đến nay chủ yếu tập trung định loại hình thái, chưa có nghiên cứu sinh thái học phân tử về Chrysomelidae. Cuối cùng, chương 1 cung cấp dữ liệu chi tiết về vị trí địa lý, địa hình lòng chảo, chế độ khí hậu bán khô hạn đặc thù (9 tháng khô, 4 tháng hạn, 2 tháng kiệt, lượng mưa trung bình dưới 700–1.000 mm/năm), mạng lưới thủy văn suối ngắn dốc và 6 kiểu thảm thực vật chính tại VQG Núi Chúa.

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

Chương 2 trình bày chi tiết và chuẩn hóa toàn bộ các quy trình thực nghiệm và phân tích số liệu:

  1. Thiết kế 5 tuyến thu mẫu định kỳ (Mái Nhà, Đá Đỏ, Ao Hồ, Núi Ông, Suối Trục) và 1 tuyến khảo sát (Đá Hang), ghi nhận tọa độ GPS và độ cao chi tiết từ 38 m đến 714 m.
  2. Quy trình thu mẫu chuẩn hóa bằng phương pháp đập 10 phút tại mỗi điểm, bảo quản cồn 96% và lưu trữ mẫu đối chứng.
  3. Kỹ thuật tách chiết ADN bằng kit Qiagen DNeasy, điều kiện phản ứng PCR khuếch đại đoạn gen cox1 (830 bp) và psbA-trnH.
  4. Quy trình tin sinh học phân định loài dựa trên mô hình GMYC (đơn/đa ngưỡng) và bPTP, so sánh đối chiếu với 141 dạng loài hình thái phân loại theo khóa Kimoto.
  5. Quy trình phân tích sinh thái học: phương pháp cửa sổ trượt (sliding window) phân định vùng chuyển tiếp; phần mềm EstimateS 9.1 đánh giá đường cong tích lũy loài; phân tích CCA trong gói vegan; phân rã đa dạng beta quần xã và tương tác trong các gói betapart, BAT, betalink; và mô hình hóa mạng lưới nhị phân bằng gói Bipartite trong môi trường R.

Chương 3: Kết quả nghiên cứu

Chương 3 trình bày các phát hiện thực nghiệm cốt lõi của luận án:

                      ĐIỀU KIỆN MÔI TRƯỜNG VQG NÚI CHÚA
                 (Gradient độ cao & Chế độ nhiệt - ẩm biến đổi)
SINH CẢNH KHÔ (< 300 m)                                   SINH CẢNH ẨM (> 300 m)
[Vùng chuyển tiếp sinh thái: 160 m - 320 m; Ranh giới phân định: 300 m]
QUẦN XÃ CHRYSOMELIDAE                                    THỰC VẬT CHỦ (cpDNA)
- 520 cá thể thu thập                                    - Giải mã locus psbA-trnH
- 141 dạng loài hình thái                                - Xác định qua BAGpipe
- 13 loài mới cho khoa học                               - Phân họ Galerucinae làm thức ăn
  (11 Monolepta, 2 Paleosepharia)                        - Ghi nhận các họ cây chủ chính
- Mô hình bPTP đồng thuận cao nhất
                    PHÂN TÍCH TƯƠNG TÁC SINH THÁI
      - Phân tích CCA: Sinh cảnh & Loài chi phối tương tác dinh dưỡng
      - Đa dạng Beta: Sự thay thế loài (βrepl, βST) giữ vai trò chủ đạo
      - Mạng lưới Bipartite: Khác biệt cấu trúc giữa rừng khô & rừng ẩm
  • Phân định vùng chuyển tiếp sinh thái: Dựa trên xu hướng giảm nhiệt độ theo độ cao và phân tích chỉ số tương đồng SØrensen-Dice qua phương pháp cửa sổ trượt, luận án xác định đai cao từ 160 m đến 320 m là vùng chuyển tiếp sinh thái giữa rừng khô đất thấp và rừng ẩm núi cao. Điểm phân chia 300 m được xác lập làm ranh giới sinh thái để phân tích so sánh giữa sinh cảnh khô (<300 m) và sinh cảnh ẩm (>300 m).
  • Đa dạng loài Chrysomelidae:
    • Từ 520 mẫu vật thu thập được, phân loại hình thái xác định 141 dạng loài hình thái thuộc nhiều phân họ (Alticinae, Bruchinae, Chlamisinae, Chrysomelinae, Clytrinae, Criocerinae, Cryptocephalinae, Eumolpinae, Galerucinae...).
    • Phát hiện và mô tả 13 loài mới cho khoa học, gồm 11 loài thuộc giống Monolepta Chevrolat (M. ninhthuanensis, M. nuichuaensis, M. trinhdinhcuongi, M. trinhnguyenthuthuyae, M. trinhbinhminhi, M. jesuzi, M. nguyenhansanhi, M. duongnami, M. nguyenquyeti, M. natus, M. thomaswagneri) và 02 loài thuộc giống Paleosepharia Laboissiere (P. nguyenbansini, P. tamdaonus).
    • Về mặt phân tử, phân tích trình tự gen ty thể cox1 bằng mô hình bPTP cho kết quả phân định loài có độ tương thích cao nhất với các dạng loài hình thái so với mô hình GMYC. Bộ dữ liệu trình tự ADN được công bố trên ngân hàng gen quốc tế EMBL-EBI với các mã truy cập LT160095 đến LT160588.
  • Thực vật chủ của phân họ Galerucinae: Ứng dụng thành công kỹ thuật giải trình tự đoạn gen cpDNA psbA-trnH và phân tích qua chương trình BAGpipe đã xác định chính xác danh mục các họ thực vật chủ của các loài thuộc phân họ Galerucinae tại khu vực nghiên cứu.
  • Biến đổi cấu trúc quần xã và tương tác theo độ cao:
    • Kết quả phân tích sắp xếp hợp quy chuẩn (CCA) cho thấy biến số sinh cảnh (độ cao <300 m và >300 m) và biến số phân loại loài bọ cánh cứng có ảnh hưởng mang ý nghĩa thống kê đến cấu trúc tương tác thực vật chủ.
    • Phân tích đa dạng beta cho thấy sự khác biệt thành phần loài giữa sinh cảnh khô và ẩm chủ yếu được thúc đẩy bởi quá trình thay thế loài ($\beta_{repl}$ chiếm ưu thế tuyệt đối so với $\beta_{rich}$; $\beta_{sim}$ chiếm ưu thế so với $\beta_{sne}$).
    • Đa dạng beta của tương tác giữa Galerucinae và thực vật chủ ($\beta_{WN}$) được đóng góp chủ yếu bởi sự thay thế loài ($\beta_{ST}$) thay vì sự thay đổi tương tác giữa các loài chia sẻ chung ($\beta_{OS}$). Mạng lưới dinh dưỡng trong sinh cảnh ẩm thể hiện cấu trúc phức tạp và mức độ liên kết chuyên hóa khác biệt rõ rệt so với sinh cảnh khô.
  • Hoạt động bảo tồn: Luận án cảnh báo các tác động làm suy thoái, phân mảnh sinh cảnh hoặc mất môi trường sống cục bộ tại VQG Núi Chúa sẽ gây đứt gãy mạng lưới liên kết thức ăn thực vật - côn trùng, đòi hỏi chiến lược bảo tồn theo diện rộng bao trùm toàn bộ dải độ cao.

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

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

  1. Dữ liệu phân tử và thành phần loài: Luận án cung cấp bộ dữ liệu mã vạch ADN ty thể gen cox1 đầu tiên tại Việt Nam cho họ Chrysomelidae (gồm các chuỗi nộp trên cơ sở dữ liệu EMBL-EBI mang số hiệu LT160095 - LT160588), khẳng định tính khả thi của mô hình bPTP trong phân định loài bọ cánh cứng ăn lá ở vùng nhiệt đới.
  2. Khám phá phân loại học: Phát hiện và mô tả 13 loài mới cho khoa học thuộc hai giống Monolepta (11 loài) và Paleosepharia (02 loài), làm giàu thêm tri thức phân loại học côn trùng cánh cứng tại Việt Nam và Đông Nam Á.
  3. Phương pháp luận nghiên cứu sinh thái dinh dưỡng: Ứng dụng thành công kỹ thuật giải mã trình tự ADN lục lạp (psbA-trnH) từ chất chứa trong ống tiêu hóa kết hợp công cụ tin sinh học BAGpipe, mở ra hướng tiếp cận chính xác để xác định thực vật chủ của côn trùng mà không cần quan sát thực địa.
  4. Lý luận sinh thái học: Cung cấp dẫn liệu thực nghiệm chứng minh rằng dọc theo dải biến đổi môi trường nhiệt độ - độ ẩm ở VQG Núi Chúa, sự phân hóa quần xã bọ cánh cứng ăn lá và mạng lưới tương tác với thực vật chủ chủ yếu diễn ra theo cơ chế thay thế loài ($\beta_{repl}$, $\beta_{ST}$), với ranh giới sinh thái rõ rệt tại đai cao 300 m.

Đóng góp về mặt thực tiễn và kiến nghị

  1. Cơ sở dữ liệu bảo tồn: Cung cấp danh mục phân loại chi tiết họ Chrysomelidae và mạng lưới thực vật chủ tương ứng, phục vụ công tác quản lý, giám sát tài nguyên sinh vật tại VQG Núi Chúa.
  2. Kiến nghị bảo tồn hệ sinh thái: Do tính chuyên hóa cao của bọ cánh cứng ăn lá và sự thay thế loài triệt để giữa các đai cao, công tác bảo tồn tại VQG Núi Chúa không thể chỉ tập trung vào vùng lõi rừng ẩm núi cao mà phải bảo vệ tính toàn vẹn, liên tục của toàn bộ cảnh quan từ các dải rừng khô hạn ven biển qua đới chuyển tiếp lên đỉnh núi, ngăn chặn nguy cơ mất môi trường sống và đứt gãy chuỗi dinh dưỡng tự nhiên.

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

Hạn chế của luận án

  • Phạm vi thu mẫu thực địa mới tập trung trong khu vực khoảng 10 km² và dải độ cao khảo sát định kỳ đạt đến 486 m (tuyến khảo sát đạt 714 m), chưa thể bao phủ toàn bộ diện tích 22.513 ha đất liền cũng như độ cao đỉnh núi cao nhất (đỉnh Cô Tuy 1.040 m) của VQG Núi Chúa.
  • Kỹ thuật xác định thực vật chủ bằng mã vạch ADN lục lạp (psbA-trnH) mới được triển khai tập trung và chuyên sâu trên phân họ Galerucinae, chưa thể mở rộng đồng loạt cho toàn bộ các cá thể của các phân họ khác trong mẫu thu thập.

Hướng nghiên cứu tiếp theo

  • Mở rộng phạm vi điều tra lên các đai độ cao trên 500 m đến 1.040 m tại VQG Núi Chúa và các khu vực rừng khô hạn lân cận ở Nam Trung Bộ.
  • Áp dụng kỹ thuật giải trình tự thế hệ mới (Next-Generation Sequencing - NGS/Metabarcoding) để giải mã toàn diện chuỗi thức ăn trong ruột của tất cả các phân họ Chrysomelidae còn lại.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị ứng dụng cao cho nhiều đối tượng:

  • Nghiên cứu sinh và học viên cao học chuyên ngành Sinh thái học, Côn trùng học và Sinh học phân tử: Cung cấp quy trình phương pháp luận mẫu mực về thiết kế tuyến thu mẫu sinh thái, quy trình phân tử tách chiết ADN từ mô côn trùng khô/ngâm cồn, phân định loài bằng thuật toán GMYC/bPTP và các đoạn mã lệnh R ứng dụng trong phân tích đa dạng beta (betapart, BAT, betalink), sắp xếp CCA (vegan) và phân tích mạng lưới (Bipartite).
  • Giảng viên và các nhà nghiên cứu phân loại học Côn trùng: Cung cấp khóa mô tả hình thái đối chiếu dữ liệu phân tử, danh lục các loài mới thuộc chi Monolepta và Paleosepharia, cùng cơ sở dữ liệu mã vạch ADN đối chứng trên EMBL-EBI.
  • Ban Quản lý Vườn quốc gia Núi Chúa và các nhà hoạch định chính sách bảo tồn: Cung cấp luận cứ khoa học về cấu trúc phân đới sinh thái và cảnh báo nguy cơ đứt gãy tương tác dinh dưỡng do suy thoái môi trường sống, hỗ trợ quy hoạch các phân khu bảo vệ nghiêm ngặt gắn liền với dải chuyển tiếp độ cao.

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

1. Luận án sử dụng những chỉ thị phân tử nào để định danh loài bọ cánh cứng và xác định thực vật chủ?

Luận án sử dụng đoạn gen ty thể Cytochrome c oxidase subunit 1 (cox1 / COI) dài 830 bp để định danh loài bọ cánh cứng ăn lá họ Chrysomelidae (kết hợp mồi TL-N-3014 và C1-J-2183 cải biên). Để xác định thực vật chủ từ thức ăn trong ruột côn trùng, luận án sử dụng đoạn gen lục lạp cpDNA vùng khoảng trống liên gen psbA-trnH (khuếch đại bằng mồi PsbA/TrnH hoặc PsbAinT2/TrnHinT1) và giải mã qua chương trình BAGpipe.

2. Ranh giới chuyển tiếp sinh thái tại VQG Núi Chúa được xác định ở độ cao nào và dựa trên cơ sở gì?

Qua phân tích xu hướng nhiệt độ và phương pháp cửa sổ trượt (sliding window) đánh giá chỉ số tương đồng SØrensen-Dice của các khoảng cao 40 m liền kề, luận án xác định đai độ cao từ 160 m đến 320 m là vùng chuyển tiếp sinh thái giữa rừng khô và rừng ẩm. Trong đó, mốc độ cao 300 m được chọn làm ranh giới phân định rõ nét nhất giữa sinh cảnh khô (<300 m) và sinh cảnh ẩm (>300 m).

3. Nghiên cứu đã phát hiện và ghi nhận được bao nhiêu loài Chrysomelidae tại VQG Núi Chúa?

Từ 520 cá thể thu thập được, phân tích hình thái học ghi nhận 141 dạng loài hình thái thuộc nhiều phân họ khác nhau. Trong đó, tác giả đã phát hiện và mô tả 13 loài mới cho khoa học, gồm 11 loài thuộc giống Monolepta và 02 loài thuộc giống Paleosepharia.

4. Phương pháp phân định loài phân tử nào cho kết quả phù hợp nhất với phân loại hình thái trong luận án?

Mô hình Poisson Tree Processes (bPTP) chạy trực tuyến trên máy chủ bPTP dựa trên cây ML xây dựng từ phần mềm RAxML cho kết quả phân định loài có độ đồng thuận và tương thích cao nhất với các dạng loài hình thái học so với các biến thể đơn ngưỡng và đa ngưỡng của mô hình GMYC.

5. Sự thay đổi cấu trúc quần xã bọ cánh cứng và tương tác dinh dưỡng theo độ cao diễn ra theo cơ chế nào?

Kết quả phân tích đa dạng beta cho thấy sự khác biệt thành phần loài của quần xã bọ cánh cứng ăn lá ($\beta_{total}$) và sự khác biệt mạng lưới tương tác với thực vật chủ ($\beta_{WN}$) giữa sinh cảnh khô và sinh cảnh ẩm chủ yếu do cơ chế thay thế loài ($\beta_{repl}$, $\beta_{ST}$) chi phối, trong khi thành phần tạo ổ/lồng ghép ($\beta_{rich}$, $\beta_{OS}$) chiếm tỷ trọng rất nhỏ.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Thị Định đã hoàn thành các mục tiêu nghiên cứu đề ra thông qua việc kết hợp công cụ mã vạch ADN ty thể (cox1), ADN lục lạp (psbA-trnH) và các phương pháp phân tích sinh thái học quần xã hiện đại. Công trình ghi nhận 141 dạng loài hình thái bọ cánh cứng ăn lá từ 520 mẫu vật thu thập tại VQG Núi Chúa, khám phá 13 loài mới cho khoa học và giải mã thành công mạng lưới thức ăn của phân họ Galerucinae. Kết quả nghiên cứu làm sáng tỏ vai trò chủ đạo của cơ chế thay thế loài dọc theo dải biến đổi độ cao (với ranh giới sinh thái 300 m), cung cấp cơ sở khoa học quan trọng cho công tác bảo tồn đa dạng sinh học tại vùng rừng bán khô hạn nhiệt đới Nam Trung Bộ.