Tổng quan luận án

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

Thạch tùng răng (Huperzia serrata (Thunb.) Trevis) thuộc họ Thông đất (Lycopodiaceae), là loài thực vật bậc cao có mạch cổ xưa phân bố tại các vùng rừng nhiệt đới và ôn đới núi cao. Loài này có giá trị dược liệu quan trọng nhờ chứa các hợp chất alkaloid, tiêu biểu là Huperzine A có hoạt tính sinh học cao dùng trong hỗ trợ điều trị các bệnh thoái hóa thần kinh, đặc biệt là bệnh Alzheimer. Do giá trị kinh tế và dược liệu lớn, các quần thể Thạch tùng răng trong tự nhiên bị khai thác thương mại quá mức dẫn đến suy giảm nghiêm trọng. Tại Việt Nam, loài đã được đưa vào danh mục thực vật rừng nguy cấp, quý, hiếm thuộc nhóm IIA theo Nghị định số 06/2019/NĐ-CP.

Bên cạnh áp lực khai thác, các quần thể Thạch tùng răng đang chịu tác động tiêu cực từ biến đổi khí hậu toàn cầu, sự gia tăng nhiệt độ, hạn hán kéo dài và tình trạng phân mảnh sinh cảnh. Đối với thực vật có đời sống cố định, khả năng duy trì quần thể trước biến động môi trường phụ thuộc chủ yếu vào tính mềm dẻo kiểu hình và tiềm năng thích nghi dựa trên nền tảng đa dạng di truyền. Mặc dù các nghiên cứu quốc tế đã tiếp cận một số khía cạnh sinh học, sinh thái và hóa dược của loài, việc đánh giá định lượng và đồng bộ ảnh hưởng của các nhân tố khí hậu đến các đặc điểm hình thái ngoài, giải phẫu hiển vi và cấu trúc di truyền của các quần thể Thạch tùng răng phân bố dọc theo chiều dài địa lý Việt Nam vẫn chưa được thực hiện đầy đủ. Luận án được triển khai nhằm giải quyết khoảng trống nghiên cứu này.

Mục tiêu nghiên cứu

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

  1. Cung cấp dữ liệu về sự biến đổi một số đặc điểm hình thái bên ngoài cũng như đặc điểm hình thái bên trong và đặc điểm di truyền trong mối tương quan với sự thay đổi của điều kiện khí hậu;
  2. Dự đoán chiều hướng biến đổi cũng như vùng sinh thái có khả năng thích nghi của loài trong bối cảnh biến đổi khí hậu toàn cầu;
  3. Đề xuất định hướng quy hoạch và bảo tồn vốn gen của loài.

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

  • Đối tượng nghiên cứu: Loài Thạch tùng răng (Huperzia serrata (Thunb.) Trevis), họ Thông đất (Lycopodiaceae).
  • Phạm vi không gian: Nghiên cứu được tiến hành trên 4 quần thể tự nhiên phân bố từ Bắc vào Nam tại Việt Nam:
    1. Vườn Quốc gia Hoàng Liên (tỉnh Lào Cai);
    2. Vườn Quốc gia Bạch Mã (tỉnh Thừa Thiên Huế);
    3. Khu Bảo tồn Thiên nhiên Ngọc Linh (tỉnh Quảng Nam);
    4. Vườn Quốc gia Bidoup – Núi Bà (tỉnh Lâm Đồng). Trong đó, đặc điểm hình thái và giải phẫu được nghiên cứu chi tiết tại 3 quần thể (Hoàng Liên, Ngọc Linh, Bidoup); đặc điểm đa dạng di truyền được đánh giá trên cả 4 quần thể (Hoàng Liên, Bạch Mã, Ngọc Linh, Bidoup).
  • Phạm vi thời gian: Luận án thực hiện từ tháng 11 năm 2016 đến tháng 11 năm 2020. Dữ liệu khí tượng phục vụ phân tích tương quan sinh thái được thu thập trong giai đoạn 2013 – 2016.
  • Giới hạn nội dung: Tập trung vào 4 nhân tố khí hậu chính: nhiệt độ không khí trung bình năm (MAT, °C), tổng lượng mưa trung bình năm (MAP, mm), tổng số giờ nắng trung bình năm (SuH, giờ), và độ ẩm không khí tương đối trung bình năm (Hu, %), cùng các chỉ số khí hậu phân chia theo mùa (mùa hè, mùa đông). Các chỉ tiêu hình thái bao gồm lá, túi bào tử, bào tử, khí khổng, tế bào biểu bì, quản bào xylem và sinh trưởng chiều cao hàng năm.

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

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

Tổng quan tài liệu trong luận án điểm lại các hướng nghiên cứu nền tảng:

  • Tác động của nhân tố sinh thái và tính mềm dẻo kiểu hình: Mayr (1970/1981) khẳng định vai trò của biến dị di truyền trong việc đảm bảo khả năng thích ứng của quần thể trước các biến động cực đoan của môi trường. Woodward (1975, 1979), Körner và Woodward (1987), Marcysiak (2012), Atkin và cộng sự đã chứng minh ảnh hưởng của nhiệt độ và độ cao đến tốc độ mở rộng phiến lá, chỉ số độ dày lá (SLA) và phản ứng thích nghi ở thực vật núi cao. Hovenden và Schoor (2006) phân tích sự thay đổi hình thái lá và mật độ khí khổng của loài sồi Nothofagus cunninghamii dưới các điều kiện chiếu sáng và đai cao khác nhau. Về phương diện tiến hóa, Pigliucci (2001), Sultan (1995, 2000, 2001), Matesanz và cộng sự (2010), Schlichting (1986, 2008) phân tích cơ chế chọn lọc tự nhiên tác động lên tính mềm dẻo kiểu hình và các "mức phản ứng ẩn" (hidden reaction norms) khi sinh vật đối mặt với sinh cảnh mới. Dudley và Schmitt (1996), Donohue và cộng sự (2000, 2001) cung cấp bằng chứng thực nghiệm trên loài Impatiens capensis về giá trị thích nghi của phản ứng kéo dài thân và thời gian ra hoa.
  • Nghiên cứu về xylem, bào tử và khí khổng: Sanio (1872), Dinwoodie (1963), White (1963), Carlquist (1975), Graaff và Baas (1974) khảo sát biến dị cấu trúc quản bào và mạch gỗ theo đai cao và vĩ độ. Ejsmond và cộng sự (2011), Pacini (1990), Pacini và Franchi (2020) nghiên cứu sự thay đổi kích thước bào tử để chống mất nước. Franks và Farquhar (2007) mô tả cơ chế đóng mở khí khổng đặc thù ở loài Huperzia prolifera.
  • Nghiên cứu về phân loại, hình thái và sinh sản chi Huperzia: Holloway (1919) và Wilce (1972) thiết lập hệ thống phân loại nhóm Selago dựa trên cấu trúc hình thái và bào tử. Wang và cộng sự (2011) mô tả cấu trúc chén truyền thể và truyền thể sinh dưỡng, đồng thời đề xuất phương pháp dùng các lớp chén truyền thể hàng năm để xác định tuổi cây. Long và cộng sự (2016) nghiên cứu chu kỳ phát triển kéo dài gần một năm của túi bào tử Huperzia serrata. Soltis và Soltis (1988, 1990) phân tích hệ thống giao phối và dòng chảy gen ở thực vật bào tử đồng hình.
  • Nghiên cứu về sinh thái, hóa dược và đa dạng di truyền Thạch tùng răng: Ai và Zhang (2005), Li và cộng sự (2005), Huang và cộng sự (2011) xác định các đặc trưng sinh cảnh chuyên hóa của loài (rừng thường xanh ẩm, độ che bóng cao, đất chua giàu mùn, phân bố theo cụm). Về dược liệu, Xie và cộng sự (2013), Gao và cộng sự (2000), Toribio và cộng sự (2007), Wu và Gu (2006), Wei và cộng sự (2008), Ma và cộng sự nghiên cứu cấu trúc, phương pháp chiết xuất và sự biến động hàm lượng Huperzine A theo mùa và sinh cảnh. Về di truyền phân tử, Huang và He (2010) sử dụng chỉ thị AFLP trên 10 quần thể H. serrata tại Trung Quốc (ghi nhận PPB = 86,5%, He = 0,20, I = 0,33); Zhang và cộng sự (2011) đánh giá 7 quần thể H. serrata var. longipetiolata (PPB = 69,38%, He = 0,272, I = 0,392, GST = 0,168); Ji và cộng sự (2014) khảo sát 4 quần thể tại vùng núi Wuling bằng AFLP.

Khoảng trống luận án lựa chọn giải quyết

Các công trình trước đây tại Trung Quốc chủ yếu sử dụng chỉ thị AFLP và khảo sát cục bộ từng khía cạnh sinh học riêng lẻ. Luận án của Nguyễn Thị Ái Minh là công trình đầu tiên tại Việt Nam kết hợp đồng thời kỹ thuật chỉ thị phân tử ISSR (Inter Simple Sequence Repeats) và SCoT (Start Codon Targeted polymorphism) cùng các phép đo giải phẫu vi thể và thống kê đa biến (RDA, hồi quy tuyến tính từng bước) để lượng hóa mức độ tác động của các nhân tố khí hậu lên hình thái và đa dạng di truyền của loài Huperzia serrata.

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

Luận án ứng dụng lý thuyết sinh thái học quần thể, lý thuyết di truyền tiến hóa (Mayr, 1970/1981) và lý thuyết về tính mềm dẻo kiểu hình (Pigliucci, 2001; Matesanz et al., 2010). Khung phân tích được xây dựng trên giả định rằng các yếu tố khí hậu địa phương (nhiệt độ, lượng mưa, ánh sáng, độ ẩm) tạo nên áp lực sinh thái chi phối trực tiếp đến các tính trạng hình thái sinh dưỡng, cấu trúc cơ quan sinh sản, giải phẫu mô dẫn, cấu trúc khí khổng và sự biệt hóa di truyền giữa các quần thể.

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

Luận án kết hợp phương pháp nghiên cứu thực địa, giải phẫu thực vật, sinh học phân tử và phân tích số liệu định lượng:

[Mẫu thực địa tại 4 quần thể (Hoàng Liên, Bạch Mã, Ngọc Linh, Bidoup)]
  1. Phương pháp phân tích khí hậu: Sử dụng tiêu chuẩn phân loại khí hậu Köppen và trích xuất dữ liệu khí tượng giai đoạn 2013 – 2016 từ các trạm khí tượng tương ứng với vùng phân bố của 4 quần thể (nhiệt độ, lượng mưa, số giờ nắng, độ ẩm theo năm và theo mùa vụ).
  2. Phương pháp hình thái học và giải phẫu thực vật:
    • Xác định tuổi và sinh trưởng chiều cao hàng năm (G, mm) dựa trên khoảng cách giữa các lớp chén truyền thể kế tiếp nhau dọc theo thân.
    • Đo lường kích thước lá: chiều dài lá ($d_{lá}$), chiều rộng lá ($r_{lá}$), diện tích lá ($S_{lá}$), chỉ số diện tích lá/khối lượng (SLA).
    • Cơ quan sinh sản: chiều dài túi bào tử ($d_{túi\ bào\ tử}$), chiều rộng túi bào tử ($r_{túi\ bào\ tử}$), diện tích túi bào tử ($S_{túi\ bào\ tử}$), kích thước bào tử.
    • Giải phẫu vi thể lá và thân: đo chiều dài ($d_{khí\ khổng}$), chiều rộng ($r_{khí\ khổng}$), diện tích khí khổng ($S_{khí\ khổng}$), chiều dài khe lỗ khí ($d_{khe\ lỗ\ khí}$), chiều rộng khe lỗ khí ($r_{khe\ lỗ\ khí}$), tỷ lệ diện tích khe lỗ khí trên diện tích lá ($S_{khe\ lỗ\ khí}/S_{lá}$), mật độ khí khổng (MK), chỉ số khí khổng (SI), mức độ xẻ thùy của tế bào biểu bì xung quanh khí khổng, và kích thước quản bào xylem.
  3. Phương pháp sinh học phân tử:
    • Tách chiết DNA tổng số từ lá tươi bằng quy trình CTAB/SDS.
    • Phân tích đa hình DNA bằng hai hệ chỉ thị phân tử: ISSR và SCoT.
    • Điện di sản phẩm PCR trên gel agarose, ghi nhận ma trận nhị phân (có band = 1, không có band = 0).
  4. Phương pháp phân tích thống kê và mô hình hóa:
    • Sử dụng phần mềm thống kê chuyên dụng (Statgraphics Centurion XV).
    • Hàm hồi quy tuyến tính đa biến từng bước (Stepwise multivariate linear regression) để đánh giá tác động của tổ hợp các nhân tố khí hậu lên từng chỉ tiêu hình thái riêng lẻ.
    • Phân tích mức độ tương quan RDA (Redundancy Analysis) để xác định tỷ lệ giải thích của ma trận khí hậu đối với tổ hợp các đặc điểm hình thái và di truyền.
    • Tính toán các chỉ số di truyền quần thể: số allele trung bình ($A$), số allele hữu hiệu ($A_e$), mức dị hợp trông đợi ($H_e$), chỉ số đa dạng Shannon ($I$), tỷ lệ band đa hình ($PPB$), chỉ số biệt hóa di truyền ($G_{ST}$), hệ số cố định ($F$), dòng chảy gen ($N_m$), khoảng cách và mức độ tương đồng di truyền theo Nei (Nei’s I), phân tích phương sai phân tử AMOVA và lập sơ đồ phân nhóm di truyền UPGMA.

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

Bố cục cấu trúc luận án

Luận án có độ dài 130 trang (không bao gồm phụ lục và tài liệu tham khảo), được kết cấu thành các phần chính: Mở đầu (4 trang), Chương 1: Tổng quan (44 trang), Chương 2: Nội dung và phương pháp nghiên cứu (12 trang), Chương 3: Kết quả nghiên cứu (68 trang), Kết luận và kiến nghị (2 trang). Đi kèm luận án là 24 bảng số liệu, 34 hình ảnh minh họa (biểu đồ, bản đồ, sơ đồ), 31 phụ lục và 262 tài liệu tham khảo.

Chương 1. Tổng quan

Chương 1 hệ thống hóa cơ sở lý luận sinh thái học và di truyền học thực vật liên quan đến tính mềm dẻo kiểu hình và phản ứng thích nghi trước biến đổi khí hậu; phân tích các cơ chế kiểm soát sự phát triển của mô dẫn xylem, hình thái phiến lá, khí khổng và bào tử; tổng thuật chi tiết lịch sử phân loại, sinh thái học quần thể và tình hình nghiên cứu về thành phần hóa học (Huperzine A) cũng như đa dạng di truyền của chi Huperzia và loài Huperzia serrata trên phạm vi quốc tế và tại Việt Nam.

Chương 2. Nội dung và phương pháp nghiên cứu

Chương 2 mô tả chi tiết địa điểm thu mẫu, đặc trưng khí hậu tại 4 khu vực nghiên cứu theo tiêu chuẩn Köppen, khoảng cách từ điểm thu mẫu đến các trạm khí tượng; quy trình giải phẫu, đo đếm các chỉ tiêu hình thái học ngoài và cấu trúc vi thể; quy trình tách chiết DNA, tối ưu hóa phản ứng PCR với các mồi ISSR và SCoT; cùng các công thức toán học và phương pháp thống kê sinh học đa biến (RDA, hồi quy tuyến tính đa biến, UPGMA, AMOVA).

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

Chương 3 trình bày toàn bộ số liệu thực nghiệm và kết quả phân tích thống kê theo 5 cụm nội dung chính:

1. Đặc điểm mẫu Thạch tùng răng thu thập và hình thái sinh sản

Cây có kiểu phân nhánh lưỡng phân đều, mang các lá bé và lá lớn xếp xen kẽ. Cấu trúc sinh sản là túi bào tử thật hình thận, màu vàng tươi, nằm ở nách lá bào tử. Bào tử thuộc kiểu Selago của nhóm "lõm - mặt rỗ" với đường viền tam giác, cạnh lõm, bề mặt có các lỗ và rãnh đặc trưng. Truyền thể sinh dưỡng mọng nước, cấu tạo gồm 1 chồi và 5 lá truyền thể, hình thành trên các chén truyền thể đính dọc theo thân.

2. Biến dị hình thái ngoài và giải phẫu lá ở 3 quần thể

  • Kích thước lá: Tổ hợp các yếu tố nhiệt độ không khí trung bình năm (MAT), MAT mùa hè, MAT mùa đông và lượng mưa năm (MAP) có ảnh hưởng chi phối đến chiều dài ($d_{lá}$) và chiều rộng lá ($r_{lá}$). MAP và MAT là hai nhân tố tác động mạnh nhất lần lượt đến diện tích lá ($S_{lá}$) và chỉ số độ dày lá (SLA). Phân tích RDA chỉ ra rằng tổ hợp các nhân tố khí hậu (MAT, MAP, Hu, SuH) giải thích được 83,76% sự biến thiên của tổ hợp kích thước lá, trong khi 16,24% biến thiên do các yếu tố khác. Kích thước lá tương quan thuận với MAT và MAP (MAP có mức độ tác động mạnh hơn MAT), tương quan nghịch nhưng không đáng kể với SuH và Hu.
  • Kích thước túi bào tử: MAT là nhân tố ảnh hưởng mạnh nhất đến chiều dài, chiều rộng và diện tích túi bào tử. Tổ hợp các nhân tố khí hậu giải thích 82,04% sự biến thiên kích thước túi bào tử (17,96% do các yếu tố khác). Kích thước túi bào tử tương quan nghịch với MAT và tương quan thuận với MAP, SuH và Hu (trong đó MAP có mức độ ảnh hưởng lớn hơn SuH và Hu).
  • Giải phẫu khí khổng và biểu bì lá: Số giờ nắng mùa hè (SuH mùa hè) ảnh hưởng mạnh nhất đến chiều dài và chiều rộng khí khổng. Nhiệt độ mùa hè (MAT mùa hè) tác động mạnh nhất đến hầu hết các đặc tính khí khổng và khe lỗ khí (ngoại trừ chiều dài khí khổng). Số giờ nắng mùa đông (SuH mùa đông) tác động mạnh nhất đến mức độ xẻ thùy của tế bào biểu bì quanh khí khổng. Phân tích RDA cho thấy khí hậu giải thích 62,74% sự biến thiên của tổ hợp đặc điểm giải phẫu lá (37,23% do yếu tố khác). Tổ hợp giải phẫu lá tương quan nghịch với MAT và MAP, tương quan thuận với SuH và Hu (MAT có tác động mạnh nhất).
  • Các tính trạng ổn định: Kích thước bào tử, độ mở khe lỗ khí và kích thước quản bào không thể hiện tương quan rõ ràng với sự khác biệt khí hậu giữa các vùng phân bố.
Nhóm tính trạng Tỷ lệ biến thiên do khí hậu giải thích (RDA) Tỷ lệ do yếu tố khác Nhân tố khí hậu tác động mạnh nhất Chiều hướng tương quan chính
Kích thước lá 83,76% 16,24% MAP (diện tích lá), MAT (độ dày lá SLA) Tương quan thuận với MAT, MAP; nghịch với SuH, Hu
Kích thước túi bào tử 82,04% 17,96% MAT (chiều dài, chiều rộng, diện tích) Tương quan nghịch với MAT; thuận với MAP, SuH, Hu
Giải phẫu vi thể lá 62,74% 37,23% MAT mùa hè (khí khổng), SuH mùa hè ($d, r$ khí khổng) Tương quan nghịch với MAT, MAP; thuận với SuH, Hu

3. Tăng trưởng chiều cao hàng năm

Được xác định thông qua việc đo chiều dài các lóng thân giữa các lớp chén truyền thể hàng năm ($G$, mm). Kết quả cho thấy tăng trưởng chiều cao thân hàng năm không có mối liên hệ biến thiên rõ rệt đối với các nhân tố khí hậu giữa các quần thể nghiên cứu.

4. Đa dạng di truyền quần thể Thạch tùng răng tại Việt Nam

Đánh giá thông qua việc kết hợp hai kỹ thuật mồi ISSR và SCoT. Kết quả ghi nhận mức độ đa dạng di truyền tổng thể ở mức loài và mức quần thể, xác định tỷ lệ band đa hình (PPB), mức dị hợp trông đợi ($H_e$), chỉ số Shannon ($I$), chỉ số biệt hóa di truyền ($G_{ST}$), hệ số dòng gen ($N_m$) và khoảng cách di truyền giữa 4 quần thể (Hoàng Liên, Bạch Mã, Ngọc Linh, Bidoup). Cây phân nhóm di truyền UPGMA phản ánh mối quan hệ nguồn gốc di truyền giữa các quần thể Bắc - Trung - Nam.

5. Ảnh hưởng của nhân tố sinh thái đến đa dạng di truyền

Phân tích tương quan giữa các chỉ số di truyền và các thông số khí tượng cho thấy tổng số giờ nắng hàng năm (SuH) là nhân tố khí hậu có ảnh hưởng mạnh nhất đến tỷ lệ band đa hình (PPB) của các quần thể Thạch tùng răng tại Việt Nam.

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

Những đóng góp mới của luận án

  1. Lượng hóa ảnh hưởng sinh thái đến hình thái học: Luận án cung cấp các dẫn liệu định lượng bằng mô hình toán sinh học, chứng minh tổ hợp khí hậu giải thích phần lớn sự biến thiên kích thước phiến lá (83,76%), kích thước túi bào tử (82,04%) và giải phẫu khí khổng (62,74%) của Thạch tùng răng tại Việt Nam.
  2. Xác định các nhân tố khí hậu chi phối then chốt:
    • Tổng lượng mưa (MAP) và nhiệt độ trung bình năm (MAT) là hai nhân tố kiểm soát chính đối với sự mở rộng diện tích và chỉ số độ dày lá.
    • Nhiệt độ trung bình năm (MAT) có mối tương quan nghịch rõ rệt với kích thước túi bào tử.
    • Số giờ nắng mùa hè (SuH mùa hè) và nhiệt độ mùa hè (MAT mùa hè) chi phối kích thước và hình thái cấu trúc bộ máy khí khổng.
    • Tổng số giờ nắng năm (SuH) là nhân tố tác động mạnh nhất đến tỷ lệ band đa hình di truyền (PPB) của quần thể.
  3. Làm sáng tỏ tính phản ứng thích nghi: Chứng minh rằng các đặc điểm kích thước bào tử, cấu trúc quản bào và tăng trưởng chiều cao thân ít biến đổi trước dao động khí hậu ngắn hạn, trong khi kích thước lá và giải phẫu khí khổng thể hiện tính mềm dẻo kiểu hình cao nhằm điều hòa thoát hơi nước và quang hợp.
  4. Dữ liệu di truyền toàn diện: Luận án là công trình đầu tiên tại Việt Nam thiết lập dữ liệu đa dạng di truyền quần thể Thạch tùng răng bằng sự kết hợp của cả hai hệ chỉ thị phân tử ISSR và SCoT.

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

  • Bảo tồn nguyên vị (In-situ): Ưu tiên bảo vệ nghiêm ngặt sinh cảnh rừng nguyên sinh lá rộng thường xanh đai cao tại 4 khu vực phân bố tự nhiên (Vườn Quốc gia Hoàng Liên, Khu Bảo tồn Thiên nhiên Ngọc Linh, Vườn Quốc gia Bidoup – Núi Bà và Vườn Quốc gia Bạch Mã). Giữ vững độ tàn che của tán rừng để duy trì độ ẩm cao và hạn chế bức xạ trực tiếp.
  • Quy hoạch bảo tồn trước biến đổi khí hậu: Cần chú trọng các khu vực phân bố có lượng mưa dồi dào và nền nhiệt độ mát mẻ, do sự gia tăng nhiệt độ không khí (MAT) có xu hướng làm suy giảm kích thước túi bào tử, ảnh hưởng bất lợi đến tiềm năng sinh sản hữu tính của loài.
  • Quản lý khai thác: Tuyệt đối cấm khai thác kiểu tận thu (nhổ bật gốc cả cụm) vì hệ rễ của loài ăn nông dưới 5cm và chu kỳ tái sinh tự nhiên từ thể giao tử ngầm đòi hỏi thời gian kéo dài (15 – 20 năm để cá thể thành thục).

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

Hạn chế của nghiên cứu

  • Do Thạch tùng răng có đặc tính sinh học khó nhân giống nhân tạo, bào tử khó nảy mầm và thể giao tử sống ngầm cộng sinh với nấm trong nhiều năm, nghiên cứu chưa thể triển khai các nghiệm thức thực nghiệm trong điều kiện vườn chung (common garden experiment) để tách biệt hoàn toàn biến dị di truyền với biến dị kiểu hình.
  • Dữ liệu khí tượng được thu thập gián tiếp từ các trạm khí tượng phụ cận trong khung thời gian 2013 – 2016, chưa phản ánh trọn vẹn các biến đổi vi khí hậu tức thời tại từng tiểu sinh cảnh dưới tán rừng.
  • Phạm vi nghiên cứu giải phẫu vi thể mới thực hiện trên 3 quần thể đại diện.

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

  • Mở rộng quan trắc vi khí hậu tại chỗ bằng các thiết bị cảm biến tự động đặt trực tiếp tại các vi sinh cảnh có Thạch tùng răng phân bố.
  • Ứng dụng công nghệ sinh học phân tử hiện đại (giải trình tự gen thế hệ mới, transcriptomics) để nghiên cứu cơ chế biểu hiện gen liên quan đến tính mềm dẻo kiểu hình và con đường sinh tổng hợp hoạt chất Huperzine A dưới áp lực nhiệt độ và ánh sáng.
  • Nghiên cứu sâu hơn về mối quan hệ cộng sinh giữa nguyên tản Thạch tùng răng và nấm nội sinh phục vụ công tác nhân giống và bảo tồn ngoại vi (ex-situ).

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị học thuật và ứng dụng thực tiễn:

  • Đối với nghiên cứu sinh và giảng viên ngành Sinh thái học, Thực vật học, Di truyền học: Cung cấp phương pháp luận kết hợp phân tích giải phẫu hiển vi thực vật với các mô hình thống kê đa biến (RDA, hồi quy đa biến từng bước) và kỹ thuật sinh học phân tử (ISSR, SCoT) trong nghiên cứu sinh thái học quần thể.
  • Đối với các nhà quản lý lâm nghiệp và bảo tồn thiên nhiên: Cung cấp cơ sở khoa học trực tiếp cho Ban Quản lý các Vườn Quốc gia Hoàng Liên, Bạch Mã, Bidoup – Núi Bà và Khu Bảo tồn Thiên nhiên Ngọc Linh trong việc xây dựng phương án quản lý, khoanh nuôi và bảo vệ loài dược liệu quý nhóm IIA.
  • Đối với lĩnh vực Dược liệu học: Cung cấp các thông tin sinh thái phục vụ việc định vùng trồng, di thực hoặc khai thác hợp lý nguồn gen Huperzia serrata chứa hoạt chất Huperzine A tại Việt Nam.

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 giá đa dạng di truyền của Thạch tùng răng?

Luận án kết hợp đồng thời hai hệ thống chỉ thị phân tử DNA là ISSR (Inter Simple Sequence Repeats) và SCoT (Start Codon Targeted polymorphism) để đánh giá đa dạng di truyền và cấu trúc di truyền của 4 quần thể Thạch tùng răng tại Việt Nam.

2. Các nhân tố khí hậu giải thích bao nhiêu phần trăm sự biến thiên kích thước lá và túi bào tử của loài?

Theo kết quả phân tích mức độ tương quan RDA, tổ hợp các nhân tố khí hậu giải thích 83,76% sự biến thiên của tổ hợp kích thước lá (16,24% do các yếu tố khác) và giải thích 82,04% sự biến thiên kích thước túi bào tử (17,96% do các yếu tố khác).

3. Nhân tố khí hậu nào ảnh hưởng mạnh nhất đến tỷ lệ band đa hình (PPB) của quần thể?

Kết quả phân tích tương quan chỉ ra rằng tổng số giờ nắng hàng năm (SuH) là nhân tố khí hậu có ảnh hưởng mạnh nhất đến tỷ lệ band đa hình (PPB) của các quần thể Thạch tùng răng nghiên cứu.

4. Tăng trưởng chiều cao hàng năm của Thạch tùng răng được xác định bằng phương pháp nào?

Tăng trưởng chiều cao thân hàng năm ($G$, mm) được xác định dựa vào khoảng cách giữa các lớp chén truyền thể kế tiếp nhau trên thân cây, do phần lớn các cá thể Thạch tùng răng trong tự nhiên tạo một lớp chén truyền thể mỗi năm và cấu trúc chén truyền thể này vẫn lưu lại trên thân nhiều năm.

5. Biến đổi kích thước túi bào tử có mối tương quan như thế nào với nhiệt độ không khí?

Nhiệt độ không khí trung bình năm (MAT) là nhân tố có ảnh hưởng mạnh nhất và có mối tương quan nghịch với các chỉ tiêu chiều dài, chiều rộng và diện tích túi bào tử của Thạch tùng răng.

Kết luận

Luận án tiến sĩ của nghiên cứu sinh Nguyễn Thị Ái Minh đã giải quyết trọn vẹn mục tiêu xác định ảnh hưởng của các nhân tố khí hậu đến đặc điểm hình thái ngoài, giải phẫu hiển vi và đa dạng di truyền của loài Thạch tùng răng (Huperzia serrata) tại Việt Nam. Bằng việc kết hợp phương pháp phân tích hình thái số lượng, kỹ thuật sinh học phân tử ISSR và SCoT cùng các mô hình thống kê đa biến RDA và hồi quy tuyến tính, công trình đã định lượng được mức độ chi phối của lượng mưa, nhiệt độ và số giờ nắng lên sự biến dị kiểu hình và cấu trúc di truyền của 4 quần thể tự nhiên. Kết quả nghiên cứu cung cấp cơ sở khoa học xác đáng cho công tác bảo tồn nguyên vị nguồn gen thực vật dược liệu quý hiếm này trong bối cảnh biến đổi khí hậu toàn cầu.