Tổng quan nghiên cứu

Thực trạng suy giảm tài nguyên rừng tại khu vực Đông Nam Bộ đang đặt ra những thách thức nghiêm trọng đối với công tác bảo tồn đa dạng sinh học. Theo số liệu thống kê lâm nghiệp, diện tích rừng tự nhiên có cây họ Dầu tại vùng Đông Nam Bộ đã suy giảm mạnh mẽ từ 1.146.275 ha (chiếm 49% diện tích toàn vùng vào năm 1959) xuống còn 834.050 ha (năm 1968), tiếp tục giảm xuống 416.900 ha (năm 1982) và chỉ còn vỏn vẹn 183.081 ha (tương đương 8% diện tích) vào năm 1992. Sự suy giảm tới hơn 84% diện tích rừng này không chỉ phá vỡ cấu trúc hệ sinh thái mà còn đẩy 33 loài cây họ Dầu trên thế giới và 11 loài tại Việt Nam vào tình trạng bị đe dọa tuyệt chủng theo các tiêu chí của Sách đỏ.

Loài Dầu song nàng (Dipterocarpus dyeri Pierre) thuộc họ Dầu (Dipterocarpaceae) là loài cây gỗ lớn vượt tán có chiều cao vượt trội trên 40 m, đóng vai trò sinh thái chủ chốt và sở hữu giá trị kinh tế đặc biệt cao nhờ nguồn gỗ tốt cùng chất lượng nhựa dầu quý phục vụ ngành công nghiệp sơn và đóng tàu. Do áp lực khai thác thương mại quá mức cùng sự phân mảnh sinh cảnh sống, loài thực vật này đã được xếp vào nhóm nguy cấp bậc VU (Sẽ nguy cấp theo tiêu chuẩn A1c,d + 2c,d) trong Sách đỏ Việt Nam năm 2007.

Luận văn tập trung giải quyết vấn đề cấp bách: Đánh giá hiện trạng di truyền, đo lường sự biến đổi cấu trúc gen và xác định các thông số thụ phấn chéo của 73 cá thể Dầu song nàng trưởng thành tại 3 khu bảo tồn trọng điểm thuộc Đông Nam Bộ bao gồm Rừng phòng hộ Tân Phú, Khu Bảo tồn Thiên nhiên và Văn hóa Đồng Nai (Mã Đà) và Vườn Quốc gia Lò Gò – Xa Mát. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp bộ cơ sở dữ liệu phân tử chuẩn xác nhằm định hình các biện pháp khoanh nuôi phục hồi, chuyển vị giống và phát triển bền vững nguồn gen thực vật bản địa quý hiếm.

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 trên nền tảng vững chắc của thuyết di truyền học quần thể và sinh học phân tử hiện đại. Trong di truyền học quần thể, quần thể được xác định là đơn vị tiến hóa cơ bản; mức độ đa dạng di truyền phản ánh khả năng thích ứng của loài trước những biến động bất lợi của môi trường sống. Hiện tượng phân mảnh sinh cảnh làm suy giảm kích thước quần thể, dẫn đến nguy cơ gia tăng hệ số giao phối cận noãn (Inbreeding coefficient - F), suy thoái giống và làm biến mất các alen hiếm.

Bên cạnh đó, lý thuyết chỉ thị phân tử chứng minh kỹ thuật vi vệ tinh (Simple Sequence Repeat - SSR hay Microsatellite) sở hữu nhiều ưu thế vượt trội so với các chỉ thị phân tử thế hệ cũ như RAPD (Random Amplified Polymorphic DNA) hay RFLP (Restriction Fragment Length Polymorphism). Chỉ thị SSR là marker đồng trội (co-dominant), gồm các đoạn trình tự lặp lại nối tiếp từ 1 đến 6 cặp bazơ (bp), cho phép nhận diện chính xác các kiểu gen dị hợp tử và đồng hợp tử ở thế hệ bố mẹ cũng như con lai F1 với độ tin cậy và khả năng lặp lại rất cao.

Các khái niệm then chốt được áp dụng trong đánh giá bao gồm: Số alen trên mỗi locus (Na), Tần số alen, Độ dị hợp tử quan sát (Ho), Độ dị hợp tử kỳ vọng (He), Hàm lượng thông tin đa hình (Polymorphism Information Content - PIC) và Chỉ số sai khác cá thể (Power of Discrimination - PD).

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

Nguồn dữ liệu thực nghiệm được thu thập từ 73 mẫu vỏ thân cây Dầu song nàng trưởng thành phân bố ngẫu nhiên theo tuyến tại 3 sinh cảnh rừng tự nhiên Đông Nam Bộ: 21 mẫu tại Rừng phòng hộ Tân Phú (độ cao 116 - 127 m so với mực nước biển), 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). Ngoài ra, nghiên cứu thu thập hạt giống từ 9 cây mẹ trội (ký hiệu từ SN01 đến SN09) tại Mã Đà đem gieo ươm tại Trung tâm Lâm nghiệp Biên Hòa để phân tích cây con.

Quy trình phân tích trong phòng thí nghiệm được thực hiện nghiêm ngặt qua các giai đoạn:

  1. Tách chiết DNA tổng số theo quy trình CTAB cải tiến (Doyle & Doyle, 1987) với việc bổ sung PVP 10% và Sodium metabisulfit 0,5% để loại bỏ triệt để các hợp chất phenolic và tạp chất; nồng độ và độ sạch của dịch DNA đạt chuẩn cao với tỷ số quang phổ OD260/OD280 ổn định trong khoảng 1,8 đến 2,0.
  2. Thực hiện phản ứng PCR khuếch đại chuỗi với thể tích 40 µl trên máy PCR System 9700, ứng dụng 9 cặp mồi microsatellite đặc hiệu (Dipt1, Dipt2, Dipt3, Dipt4, Dipt6, Dipt7, Dipt8, Shc2, Shc7) với nhiệt độ gắn mồi từ 54°C đến 56°C qua 30 chu kỳ nhiệt.
  3. Điện di kiểm tra trên gel Agarose 1% và phân tách sản phẩm PCR siêu nhỏ trên gel Polyacrylamide 5% đến 6% nhuộm GelRed Nucleic Acid Stain.

Lý do lựa chọn phương pháp phân tích này là nhằm thu được độ phân giải kích thước đoạn khuếch đại ở mức chênh lệch chỉ 2 bp. Dữ liệu phổ điện di được số hóa và xử lý bằng các phần mềm chuyên dụng: GenAlEx xác định các chỉ số di truyền quần thể, Arlequin 3.1 phân tích biến lượng phân tử AMOVA, POPTREE2 vẽ cây phả hệ di truyền theo phương pháp UPGMA, và MLTR (Ritland) xác định chính xác các hệ số thụ phấn chéo (tm, ts) và tự thụ phấn (s).

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

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

Phân tích 73 mẫu Dầu song nàng trên 9 cặp mồi vi vệ tinh đã ghi nhận những phát hiện nổi bật:

Thứ nhất, phát hiện tổng cộng 35 alen khác nhau trên 9 locus nghiên cứu, với kích thước phân đoạn dao động từ 114 bp đến 286 bp. Tỷ lệ locus đa hình đạt mức tuyệt đối 100% ở cả 3 quần thể khảo sát. Số alen trung bình trên một locus đạt 4,4; dao động từ 3 alen (ở các locus Dipt1, Dipt4, Dipt6) đến 5 alen (ở các locus Dipt2, Dipt7).

Thứ hai, hàm lượng thông tin đa hình (PIC) của các cặp mồi đạt giá trị trung bình là 0,491. Trong đó, locus Dipt7 đạt giá trị PIC cao nhất là 0,745 (thể hiện khả năng phân biệt cực kỳ mạnh mẽ), tiếp theo là Dipt2 đạt 0,665, và thấp nhất là locus Dipt1 với giá trị 0,219.

Thứ ba, chỉ số sai khác giữa các cặp cá thể (PD) đạt giá trị trung bình là 0,629; với 2 locus đạt giá trị vượt trội trên 0,80 là Dipt2 (0,837) và Dipt7 (0,815). Chỉ số sai khác của mồi (Rp) biến động từ 2,244 (Dipt1) đến 3,168 (Dipt4) với giá trị trung bình toàn bộ đạt 2,614.

Thứ tư, nghiên cứu phát hiện sự hiện diện của nhiều alen hiếm (tần số xuất hiện dưới 5% hoặc dưới 10%). Điển hình như alen 150 tại locus Dipt7 chỉ chiếm 4,8% ở quần thể Mã Đà và 3,7% ở Tân Phú, trong khi hoàn toàn vắng bóng tại Lò Gò – Xa Mát. Ngược lại, tại locus Dipt1, alen 193 chiếm tỷ lệ áp đảo lên đến 90% ở quần thể Lò Gò – Xa Mát nhưng phân bố đồng đều hơn ở Mã Đà và Tân Phú.

Thảo luận kết quả

Khi xem xét cấu trúc di truyền, các dữ liệu có thể được mô hình hóa và trình bày rất trực quan thông qua bảng ma trận khoảng cách di truyền Nei (1978) kết hợp biểu đồ phân nhánh hình cây (Dendrogram) xây dựng bằng thuật toán UPGMA. Đồ thị phân cụm chỉ ra mối quan hệ di truyền chặt chẽ giữa hai quần thể Mã Đà và Tân Phú do cùng nằm trong vùng sinh thái rừng tự nhiên tỉnh Đồng Nai, trong khi quần thể Vườn Quốc gia Lò Gò – Xa Mát (Tây Ninh) có khoảng cách di truyền tách biệt rõ rệt.

So sánh với các nghiên cứu di truyền trên những loài cây hạt trần và cây họ Dầu bản địa khác tại Việt Nam như Thông Xuân Nha (Pinus armandii subsp. xuannhaensis với PIC trung bình chỉ đạt 0,324), loài Dầu song nàng tại Đông Nam Bộ vẫn duy trì được độ đa dạng di truyền ở mức khá cao (PIC trung bình đạt 0,491; PD trung bình đạt 0,629).

Tuy nhiên, sự biến động lớn về tần số alen giữa các quần thể (chẳng hạn alen 120 chiếm tới 61,1% tại Tân Phú nhưng chiếm tỷ lệ thấp ở các vùng khác) phản ánh hệ quả của quá trình phân mảnh sinh cảnh sống kéo dài trong hơn 30 năm qua. Diện tích rừng thu hẹp làm cản trở dòng gen tự nhiên (Gene flow) giữa các địa phương, dẫn tới nguy cơ tích lũy hiện tượng cận noãn cục bộ tại từng tiểu quần thể biệt lập. Kết quả phân tích hệ số giao phối từ 9 cây mẹ trội bước đầu cảnh báo cần có biện pháp trao đổi nguồn phấn để nâng cao sức sống di truyền của các thế hệ tái sinh tự nhiên.

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

Dựa trên các dẫn liệu phân tử và đặc điểm sinh thái học của loài, 4 nhóm giải pháp chiến lược được đề xuất triển khai cụ thể:

  1. Thiết lập hành lang trao đổi di truyền và di dời nguồn gen (In-situ & Ex-situ): Chủ động thu thập nguồn hạt giống và cành ghép từ các cá thể mang alen hiếm tại Khu Bảo tồn Thiên nhiên và Văn hóa Đồng Nai và Vườn Quốc gia Lò Gò – Xa Mát để thụ phấn nhân tạo bổ sung hoặc trồng xen dặm. Mục tiêu cụ thể: Nâng cao chỉ số dị hợp tử kỳ vọng (He) toàn quần thể thêm 15% trong vòng 5 năm tới. Chủ thể thực hiện: Chi cục Kiểm lâm phối hợp cùng Ban quản lý Rừng phòng hộ Tân Phú, Khu Bảo tồn Mã Đà và Vườn Quốc gia Lò Gò – Xa Mát.

  2. Xây dựng vườn giống và nhân giống vô tính từ hệ thống cây mẹ trội: Ứng dụng quy trình nhân giống từ 9 cây mẹ trội (SN01 đến SN09) đã được định danh phân tử, kết hợp phương pháp giâm hom cành non (xử lý auxin nồng độ thích hợp) từ cây con 1 năm tuổi để đạt tỷ lệ ra rễ trên 80%. Mục tiêu: Sản xuất 10.000 cây giống sạch bệnh, chuẩn di truyền trong giai đoạn 2024–2028 phục vụ trồng rừng phòng hộ và phục hồi sinh thái. Chủ thể thực hiện: Trung tâm Lâm nghiệp Biên Hòa và các Viện nghiên cứu lâm nghiệp.

  3. Áp dụng kỹ thuật lâm sinh xúc tiến tái sinh tự nhiên dưới tán rừng: Do hạt Dầu song nàng chỉ nảy mầm trong điều kiện râm mát và mất sức nảy mầm nhanh sau 1 - 2 tuần nếu thiếu ẩm, cần điều chỉnh độ tàn che của tán rừng duy trì ở mức tối ưu từ 0,5 đến 0,6. Đồng thời, bảo vệ nghiêm ngặt các chồi gốc tái sinh có đường kính gốc chặt dưới 20 cm. Mục tiêu: Đảm bảo mật độ cây tái sinh đạt trên 600 cây/ha và tỷ lệ sống sót của cây chồi đạt trên 85% sau 3 năm can thiệp. Chủ thể thực hiện: Các đội bảo vệ rừng chuyên trách địa phương.

  4. Định kỳ kiểm toán đa dạng di truyền bằng chỉ thị SSR: Xây dựng quy chế giám sát biến động gen định kỳ 5 năm một lần bằng bộ 9 chỉ thị SSR đặc hiệu (tập trung vào các locus có chỉ số PIC cao như Dipt7, Dipt2, Shc2) nhằm phát hiện sớm các dấu hiệu suy thoái cận noãn ở các thế hệ con cháu. Chủ thể thực hiện: Các đơn vị nghiên cứu Công nghệ sinh học và Quản lý tài nguyên rừng.

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

Nội dung luận văn mang lại giá trị thực tiễn và tính ứng dụng cao cho 4 nhóm đối tượng trọng tâm:

  • Ban quản lý các Khu bảo tồn thiên nhiên và Vườn quốc gia: Cung cấp bản đồ phân bố alen và khoảng cách di truyền chuẩn xác tại các khu vực Mã Đà, Tân Phú, Lò Gò – Xa Mát, hỗ trợ trực tiếp cho việc hoạch định ranh giới vùng bảo tồn nghiêm ngặt và xây dựng kế hoạch cứu hộ loài nguy cấp.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Công nghệ sinh học, Lâm nghiệp: Là tài liệu tham khảo mẫu mực về quy trình tách chiết DNA mẫu mô gỗ/vỏ cây bằng phương pháp CTAB cải tiến, kỹ thuật tối ưu hóa phản ứng PCR-SSR và quy trình tính toán các chỉ số thống kê quần thể bằng GenAlEx và Arlequin.
  • Kỹ sư lâm sinh và cán bộ kỹ thuật tại các Trung tâm giống cây rừng: Nắm bắt chính xác đặc tính sinh lý nảy mầm của hạt giống họ Dầu, kỹ thuật lựa chọn cây trội lấy giống và phương pháp duy trì độ tàn che tán rừng từ 0,5 đến 0,6 trong xúc tiến tái sinh tự nhiên.
  • Các tổ chức phi chính phủ và các quỹ tài trợ bảo tồn đa dạng sinh học (như NAFOSTED, IUCN, WWF): Khai thác nguồn số liệu định lượng xác thực về mức độ đa hình phân tử và thực trạng suy thoái gen để làm căn cứ thẩm định, phê duyệt và cấp vốn cho các dự án bảo tồn thực vật rừng nhiệt đới tại Đông Nam Bộ.

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

Chỉ thị SSR có ưu thế gì vượt trội so với RAPD trong nghiên cứu cây Dầu song nàng?
Chỉ thị SSR là marker đồng trội, có khả năng phát hiện chính xác cả kiểu gen đồng hợp tử và dị hợp tử với độ lặp lại rất cao. Ngược lại, kỹ thuật RAPD là marker át gen, độ ổn định thấp và khó phân biệt trạng thái dị hợp tử ở thế hệ phân ly.

Tỷ lệ 100% locus đa hình và 35 alen phát hiện được phản ánh điều gì?
Kết quả 100% locus đa hình trên cả 9 cặp mồi chứng minh mức độ phong phú di truyền của 73 mẫu Dầu song nàng nghiên cứu còn tương đối tốt. Điều này chứng tỏ tiềm năng thích nghi và khả năng phục hồi của loài vẫn còn nếu được bảo tồn đúng cách.

Quá trình phân mảnh sinh cảnh ảnh hưởng như thế nào đến loài Dầu song nàng tại Đông Nam Bộ?
Sự suy giảm diện tích rừng từ 49% xuống 8% đã cô lập các tiểu quần thể, làm đứt gãy dòng trao đổi phấn hoa giữa các cánh rừng Mã Đà, Tân Phú và Lò Gò – Xa Mát, dẫn tới nguy cơ biến mất các alen hiếm và tăng giao phối cận noãn.

Hạt giống Dầu song nàng có đặc tính sinh học gì cần lưu ý trong nhân giống?
Hạt Dầu song nàng mất sức nảy mầm rất nhanh chỉ sau 1 đến 2 tuần nếu bị khô hạn và chỉ có thể nảy mầm thuận lợi trong môi trường đất đủ ẩm dưới bóng râm với độ tàn che tán rừng đạt từ 0,5 đến 0,6.

Chỉ số PIC trung bình 0,491 của bộ mồi có ý nghĩa lâm sinh như thế nào?
Giá trị PIC đạt 0,491 (trong đó locus Dipt7 đạt 0,745) khẳng định bộ 9 cặp mồi microsatellite được lựa chọn có giá trị thông tin cao, là công cụ phân tử chuẩn xác để kiểm định nguồn gốc giống và đánh giá đa dạng di truyền cho các dự án bảo tồn tiếp theo.

Kết luận

Luận văn đã hoàn thành toàn diện các mục tiêu nghiên cứu với các đóng góp cốt lõi:

  • Đã xây dựng thành công quy trình tách chiết DNA tổng số đạt độ tinh sạch cao (OD260/OD280 từ 1,8 đến 2,0) từ mô vỏ thân cây Dầu song nàng bằng phương pháp CTAB cải tiến.
  • Xác định được 35 alen trên 9 locus vi vệ tinh SSR với tỷ lệ đa hình đạt 100%, số alen trung bình đạt 4,4 alen/locus trên 73 cá thể tại 3 quần thể tự nhiên.
  • Đánh giá định lượng bộ tham số di truyền phân tử với giá trị PIC trung bình đạt 0,491; chỉ số PD đạt 0,629 và chỉ số Rp đạt 2,614.
  • Làm rõ cấu trúc phân hóa di truyền giữa quần thể Lò Gò – Xa Mát với hai quần thể Mã Đà và Tân Phú, phản ánh tác động rõ rệt của hiện tượng phân mảnh sinh cảnh.
  • Đề xuất hệ thống 4 giải pháp lâm sinh và sinh học phân tử đồng bộ nhằm bảo tồn và phát triển bền vững nguồn gen loài Dầu song nàng nguy cấp.

Trong giai đoạn 2024–2030, các cơ quan quản lý và viện nghiên cứu cần tiếp tục mở rộng phân tích di truyền trên toàn bộ vùng phân bố Nam Trung Bộ và Tây Nguyên, kết hợp xây dựng ngân hàng gen ngoại vi (Ex-situ). Hãy chủ động áp dụng ngay bộ chỉ thị phân tử SSR và các khuyến nghị kỹ thuật từ luận văn vào công tác quản lý rừng thực tiễn để bảo vệ vĩnh viễn loài cây quý báu này khỏi nguy cơ tuyệt chủng.