Tổng quan nghiên cứu

Rừng nhiệt đới vùng Đông Nam Bộ từng sở hữu hệ sinh thái cây họ Dầu vô cùng phong phú, nhưng đang chịu sức ép nghiêm trọng từ quá trình suy thoái sinh cảnh. Theo số liệu thống kê lâm nghiệp, diện tích rừng họ Dầu tại khu vực này đã sụt giảm mạnh mẽ từ 1.146.275 ha (chiếm khoảng 49% diện tích toàn vùng vào năm 1959) xuống chỉ còn 183.081 ha (tương đương khoảng 8% vào năm 1992). Trong số các loài cây gỗ lớn bản địa, Dầu song nàng (Dipterocarpus dyeri Pierre) là loài cây gỗ ưu thế tầng vượt tán có giá trị kinh tế và sinh thái đặc biệt cao. Tuy nhiên, do nạn khai thác gỗ quá mức cùng với sự phân mảnh môi trường sống, loài này đã bị suy giảm nghiêm trọng về số lượng cá thể và hiện được xếp vào cấp Sẽ nguy cấp (VU A1c,d + 2c,d) theo Sách đỏ Việt Nam.

Trước thực trạng nguồn gen bản địa đứng trước nguy cơ xói mòn, nghiên cứu được triển khai nhằm mục tiêu đánh giá toàn diện mức độ đa dạng di truyền và các thông số thụ phấn chéo của loài Dầu song nàng tại 3 quần thể tự nhiên trọng điểm thuộc vùng Đông Nam Bộ, bao gồm Rừng phòng hộ Tân Phú (21 mẫu), Khu Bảo tồn Thiên nhiên và Văn hóa Đồng Nai tại Mã Đà (25 mẫu), và Vườn Quốc gia Lò Gò – Xa Mát (27 mẫu), với tổng quy mô 73 cá thể trưởng thành.

Nghiên cứu ứng dụng kỹ thuật chỉ thị phân tử vi vệ tinh (SSR) với 9 locus đặc hiệu nhằm xác định cấu trúc di truyền và mức độ biến dị phân tử. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi cung cấp hệ thống dữ liệu sinh học phân tử chuẩn xác với chỉ số thông tin đa hình PIC trung bình đạt 0,491. Đây là căn cứ khoa học vững chắc giúp các nhà quản lý lâm nghiệp và bảo tồn thiên nhiên xây dựng chiến lược phục hồi, quy hoạch vùng thu hái giống và bảo tồn nguồn gen bền vững cho loài cây gỗ quý này trong giai đoạn tới.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng lý thuyết di truyền học quần thể kinh điển và lý thuyết sinh thái học bảo tồn phân tử. Trong cấu trúc rừng nhiệt đới, họ Dầu (Dipterocarpaceae) bao gồm khoảng 580 đến 680 loài thân gỗ chiếm ưu thế, trong đó có 33 loài đang bị đe dọa ở mức độ toàn cầu và 11 loài được ghi nhận trong Sách đỏ Việt Nam. Quá trình tiến hóa và khả năng thích nghi của quần thể thực vật phụ thuộc chặt chẽ vào mức độ đa dạng di truyền được duy trì qua các thế hệ.

Khung lý thuyết của luận văn tập trung vào 3 mô hình cốt lõi:

  • Lý thuyết trôi dạt di truyền và phân mảnh sinh cảnh (Genetic Drift & Habitat Fragmentation): Khi nơi sống bị thu hẹp và chia cắt, kích thước hữu hiệu của quần thể giảm mạnh, dẫn đến hiện tượng giao phối cận noãn và suy giảm tính đa dạng alen.
  • Mô hình cân bằng Hardy-Weinberg: Đánh giá sự biến động tần số alen và mức độ lệch pha dị hợp tử trong các quần thể tự nhiên chịu tác động nhân vi.
  • Lý thuyết hệ thống sinh sản và thụ phấn chéo (Mating System Theory): Phân tích tương quan giữa thụ phấn chéo đa locus và thụ phấn cận noãn trong việc định hình cấu trúc di truyền thế hệ con.

Các khái niệm then chốt được vận dụng xuyên suốt gồm: Đa dạng di truyền (Genetic Diversity), Chỉ số thông tin đa hình (Polymorphic Information Content - PIC), Chỉ số phân biệt (Power of Discrimination - PD), Chỉ số phân giải (Resolving Power - Rp), Hệ số cận noãn (F) và Mức độ dị hợp tử quan sát cũng như kỳ vọng.

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

Nghiên cứu thu thập dữ liệu sơ cấp từ 73 mẫu vỏ thân cây Dầu song nàng trưởng thành tại 3 khu rừng bảo tồn tiêu biểu ở Đông Nam Bộ: 21 mẫu tại Rừng phòng hộ Tân Phú (độ cao 116 - 127 m), 25 mẫu tại Khu Bảo tồn Thiên nhiên và Văn hóa Đồng Nai (độ cao 98 - 119 m) và 27 mẫu tại Vườn Quốc gia Lò Gò – Xa Mát (độ cao 10 - 20 m). Phương pháp lấy mẫu ngẫu nhiên theo tuyến khảo sát thực địa được áp dụng kết hợp với quy phạm kỹ thuật xây dựng rừng giống và vườn giống lâm nghiệp (QPN 15-93) để chọn 9 cây trội ưu tú phục vụ phân tích thụ phấn chéo.

Quy trình phân tích trong phòng thí nghiệm được tiến hành theo các bước chuẩn hóa:

  • Tách chiết DNA tổng số: Áp dụng phương pháp CTAB có cải tiến theo Doyle và Doyle (1987) nhằm loại bỏ triệt để hợp chất thứ sinh và tannin trong mô vỏ cây gỗ, đạt tỷ lệ quang phổ OD260/OD280 tối ưu từ 1,8 đến 2,0.
  • Phân tích PCR-SSR: Sử dụng 9 cặp mồi microsatellite đặc hiệu (Dipt1, Dipt2, Dipt3, Dipt4, Dipt6, Dipt7, Dipt8, Shc2, Shc7) với chu nhiệt 30 chu kỳ trên hệ thống máy PCR System 9700.
  • Điện di phân tích: Sản phẩm PCR được phân tách trên gel polyacrylamide nồng độ 6%, nhuộm huỳnh quang GelRed và chụp ảnh dưới ánh sáng tử ngoại.
  • Lý do lựa chọn phương pháp SSR: SSR là chỉ thị phân tử đồng trội, có độ lặp lại cao, khả năng phát hiện đa hình vượt trội và nhận diện chính xác trạng thái dị hợp tử ở thế hệ bố mẹ lẫn con cái so với chỉ thị RAPD hay RFLP.
  • Xử lý thống kê: Sử dụng phần mềm chuyên dụng GenAlEx để tính toán tần số alen, PIC, PD, FST; phần mềm Arlequin 3.1 để phân tích biến thiên phân tử AMOVA; POPTREE2 để thiết lập cây quan hệ di truyền UPGMA; và phần mềm MLTR để xác định thông số thụ phấn chéo. Toàn bộ lộ trình nghiên cứu thực địa và phân tích phòng thí nghiệm được hoàn thiện trong giai đoạn 2017 – 2018.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Nghiên cứu đã mang lại các phát hiện khoa học định lượng có giá trị cao về mặt sinh học phân tử của loài Dầu song nàng tại khu vực Đông Nam Bộ:

  • Mức độ đa hình di truyền đạt 100%: Tất cả 9 cặp mồi microsatellite khảo sát đều cho kết quả đa hình rõ rệt trên 73 mẫu nghiên cứu, khuếch đại thành công tổng cộng 35 allele khác nhau với kích thước đoạn DNA dao động từ 114 bp đến 286 bp. Số alen trung bình trên mỗi locus đạt 4,4 alen (biến thiên từ 3 alen ở các locus Dipt1, Dipt4, Dipt6 đến 5 alen ở các locus Dipt2, Dipt7).
  • Hiệu quả phân tích vượt trội của hệ thống chỉ thị: Giá trị PIC trung bình của 9 locus đạt 0,491, trong đó cao nhất là locus Dipt7 đạt 0,745 và thấp nhất là Dipt1 đạt 0,219. Chỉ số phân biệt PD trung bình đạt 0,629, với hai locus Dipt2 (0,837) và Dipt7 (0,815) thể hiện khả năng phân biệt cá thể cực kỳ chuẩn xác. Chỉ số phân giải Rp đạt trung bình 2,614, cao nhất tại locus Dipt4 đạt 3,168.
  • Biến động tần số alen và sự xuất hiện alen hiếm: Nghiên cứu phát hiện sự sai khác rõ rệt về tần số alen giữa 3 quần thể. Tại locus Dipt7, alen 150 xuất hiện dưới dạng alen hiếm với tần số 4,8% ở Mã Đà và 3,7% ở Tân Phú, nhưng hoàn toàn biến mất ở quần thể Lò Gò – Xa Mát. Ngược lại, alen 120 chiếm ưu thế áp đảo tại Tân Phú (61,1%) trong khi ở Mã Đà chỉ đạt 28,6%.
  • Cấu trúc sinh sản và khả năng duy trì dòng gen: Kết quả phân tích thụ phấn trên 9 gia đình cây trội ở Mã Đà chỉ ra rằng loài Dầu song nàng vẫn duy trì tỷ lệ thụ phấn chéo cao, đóng vai trò nền tảng bảo vệ mức độ dị hợp tử tự nhiên trong điều kiện sinh cảnh bị cô lập.

Thảo luận kết quả

Mức độ đa dạng di truyền của Dầu song nàng ở mức trung bình khá (PIC = 0,491) cao hơn đáng kể so với loài Thông Xuân Nha (Pinus armandii subsp. xuannhaensis, với PIC = 0,324), và tương đồng với các nghiên cứu trên loài Dầu rái (Dipterocarpus alatus) tại cùng khu vực Đông Nam Bộ. Sự biến đổi tần số alen và nguy cơ thất thoát các alen hiếm tại Vườn Quốc gia Lò Gò – Xa Mát bắt nguồn từ nguyên nhân phân mảnh sinh cảnh kéo dài. Việc diện tích rừng suy giảm khoảng 41% trong giai đoạn 1959 - 1992 đã tạo ra các rào cản địa lý, làm suy yếu sự trao đổi phấn hoa giữa các vùng rừng cách biệt.

Để biểu diễn trực quan các phát hiện này trong các báo cáo chuyên sâu, dữ liệu có thể được trình bày thông qua biểu đồ cột so sánh tần số alen giữa 9 locus của 3 quần thể, kết hợp biểu đồ hình quạt phân tích phương sai phân tử AMOVA (thể hiện tỷ lệ biến thiên trong nội bộ quần thể so với giữa các quần thể) và sơ đồ nhánh cây quan hệ di truyền UPGMA dựa trên hệ số tương đồng Nei (1978). Cấu trúc biểu đồ này giúp làm nổi bật khoảng cách di truyền gần gũi giữa quần thể Mã Đà và Tân Phú so với khoảng cách xa hơn tới quần thể Lò Gò – Xa Mát.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu sinh học phân tử và đặc điểm sinh thái của loài Dầu song nàng, 4 nhóm giải pháp mang tính hành động cao được đề xuất:

  • Thiết lập hành lang sinh thái kết nối dòng gen: Xây dựng các dải rừng liên kết sinh cảnh giữa Khu Bảo tồn Mã Đà và Rừng phòng hộ Tân Phú trong giai đoạn 2024 – 2028 nhằm nâng cao hệ số trao đổi di truyền lên mức trên 1,0; cơ quan chịu trách nhiệm thực hiện là Chi cục Kiểm lâm các tỉnh Đồng Nai và Tây Ninh phối hợp cùng các Ban quản lý rừng đặc dụng.
  • Tăng cường bảo tồn nội vi nghiêm ngặt tại quần thể gốc: Thực hiện khoanh vùng bảo vệ nguyên vẹn 100% các cá thể Dầu song nàng mẹ trưởng thành tại Mã Đà và Tân Phú, duy trì độ tàn che tán rừng tối ưu từ 0,5 đến 0,6 nhằm kích thích tái sinh chồi và hạt tự nhiên trong khung thời gian 12 đến 24 tháng tới; chủ thể thực hiện là Ban quản lý Khu bảo tồn Đồng Nai và Vườn Quốc gia Lò Gò – Xa Mát.
  • Triển khai bảo tồn ngoại vi và xây dựng vườn giống nhân tạo: Thu thập hạt giống từ 9 cây trội ưu tú đã định danh (từ SN01 đến SN09) tại rừng Mã Đà, tiến hành gieo ươm và nhân giống tối thiểu 5.000 cây con đạt chuẩn di truyền tại Vườn ươm Biên Hòa trước quý IV năm 2026; giao Trung tâm Lâm nghiệp Biên Hòa và Viện Công nghệ sinh học Lâm nghiệp chủ trì kỹ thuật.
  • Ứng dụng bộ chỉ thị SSR vào quy chuẩn kiểm định giống lâm nghiệp: Ban hành quy trình ứng dụng 9 locus SSR (đặc biệt là Dipt2, Dipt7, Dipt4) làm công cụ giám sát nguồn gốc và độ đa dạng di truyền (đảm bảo chỉ số PIC đạt từ 0,45 trở lên) cho các lô giống phục vụ chương trình trồng rừng gỗ lớn từ năm 2025; cơ quan chủ quản là Cục Lâm nghiệp - Bộ Nông nghiệp và Phát triển Nông thôn.

Đối tượng nên tham khảo luận văn

  • Học viên cao học và nghiên cứu sinh ngành Công nghệ sinh học và Lâm học: Tiếp cận quy trình tách chiết DNA bằng CTAB cải tiến từ mô thân gỗ chứa nhiều hợp chất thứ sinh, đồng thời tham khảo phương pháp vận hành 9 chỉ thị SSR trong nghiên cứu di truyền thực vật rừng nhiệt đới.
  • Ban quản lý các Vườn Quốc gia và Khu bảo tồn thiên nhiên: Khai thác trực tiếp cơ sở dữ liệu di truyền của 73 mẫu cá thể để xây dựng kế hoạch quản lý nguồn gen, xác định các quần thể ưu tiên bảo vệ tại Đồng Nai và Tây Ninh.
  • Cán bộ kỹ thuật tại các Trung tâm giống lâm nghiệp và vườn ươm: Ứng dụng nguồn giống thu thập từ 9 cây trội ưu tú (SN01 đến SN09) cùng các thông số sinh thái về độ ẩm đất và độ che bóng (tàn che 0,5 - 0,6) để nâng cao tỷ lệ sống của cây giống khi ươm tạo.
  • Nhà hoạch định chính sách bảo tồn đa dạng sinh học: Sử dụng các bằng chứng thực nghiệm về sự phân hóa alen và hệ số FST để xây dựng quy hoạch các khu bảo tồn liên vùng và ban hành danh mục nguồn gen cây họ Dầu cần ưu tiên tài trợ nghiên cứu.

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

Tại sao chỉ thị vi vệ tinh SSR lại vượt trội hơn chỉ thị RAPD trong nghiên cứu di truyền cây họ Dầu?
Chỉ thị SSR là marker đồng trội, có độ lặp lại và độ ổn định cao trong phản ứng PCR. Khác với RAPD chỉ phát hiện phân tử dạng trội, SSR phát hiện chính xác 35 allele dị hợp tử trên 9 locus nghiên cứu, cho phép định danh cá thể và xác định tần số alen chuẩn xác.

Vì sao loài Dầu song nàng lại đứng trước nguy cơ suy thoái di truyền nghiêm trọng?
Diện tích rừng họ Dầu ở Đông Nam Bộ đã giảm sút từ 49% xuống 8%, tạo nên sự phân cắt sinh cảnh nặng nề. Thêm vào đó, hạt Dầu song nàng mất sức nảy mầm chỉ sau 1 đến 2 tuần nếu không gặp mưa ẩm, khiến khả năng phát tán tự nhiên bị hạn chế tối đa.

Cỡ mẫu và địa điểm khảo sát trong luận văn được lựa chọn như thế nào?
Luận văn đã phân tích tổng cộng 73 mẫu vỏ thân cây trưởng thành được thu thập ngẫu nhiên tại 3 vùng sinh cảnh đại diện: 21 mẫu tại Rừng phòng hộ Tân Phú, 25 mẫu tại Khu Bảo tồn Thiên nhiên và Văn hóa Đồng Nai (Mã Đà), và 27 mẫu tại Vườn Quốc gia Lò Gò – Xa Mát.

Chỉ số PIC trung bình 0,491 của 9 locus SSR mang lại ý nghĩa lâm sinh gì?
Chỉ số PIC đạt 0,491 phản ánh độ đa dạng di truyền ở mức trung bình khá và khẳng định bộ mồi SSR lựa chọn có tính thông tin cao. Đây là công cụ hữu hiệu để phân biệt các dòng vô tính và kiểm định chất lượng di truyền của các lô hạt giống phục vụ trồng rừng.

Điều kiện sinh thái nào là tối ưu nhất cho sự tái sinh của cây Dầu song nàng?
Cây Dầu song nàng tái sinh chồi tốt nhất trên các gốc chặt có đường kính từ 3 đến 7 cm và chiều cao dưới 20 cm, đặt dưới điều kiện đất ẩm mát và độ tàn che tán rừng duy trì trong khoảng 0,5 đến 0,6.

Kết luận

  • Xác định thành công cấu trúc và mức độ đa dạng di truyền của loài Dầu song nàng trên quy mô 73 cá thể tại 3 quần thể tự nhiên ở Đông Nam Bộ.
  • Chứng minh tính hữu hiệu của 9 cặp chỉ thị SSR khi phát hiện 35 alen với chỉ số PIC trung bình 0,491 và chỉ số phân biệt PD đạt 0,629.
  • Ghi nhận sự sai khác rõ rệt về tần số alen giữa các quần thể và cảnh báo hiện tượng thất thoát alen hiếm do chia cắt sinh cảnh.
  • Đánh giá các thông số sinh sản then chốt và tuyển chọn thành công 9 cây trội ưu tú tại rừng Mã Đà phục vụ công tác nhân giống.
  • Đề xuất khung giải pháp quản lý tích hợp kết hợp giữa bảo tồn nguyên vị, bảo tồn chuyển vị và thiết lập hành lang sinh thái rừng.

Đóng góp quan trọng nhất của luận văn là đã đặt nền móng cơ sở dữ liệu sinh học phân tử đầu tiên cho loài Dầu song nàng tại Việt Nam, mở ra hướng đi ứng dụng chỉ thị DNA vào quản lý giống cây rừng bản địa. Trong giai đoạn 2024 – 2026, các cơ quan lâm nghiệp cần tiếp tục mở rộng quy mô giải trình tự toàn bộ hệ gen và nhân rộng mô hình phục hồi rừng gỗ lớn họ Dầu trên thực địa.