Tổng quan nghiên cứu

Bách xanh (Calocedrus macrolepis Kurz) là loài cây gỗ lá kim bản địa quý hiếm thuộc họ Hoàng đàn (Cupressaceae), mang giá trị kinh tế và khoa học đặc biệt cao. Hiện nay, theo ước tính điều tra thực địa, trên phạm vi toàn quốc không còn quá 500 cây Bách xanh có đường kính thân trên 10 cm. Trước áp lực khai thác gỗ thương mạị, làm đồ mỹ nghệ, sản xuất bột hương cùng với tình trạng phá rừng làm nương rẫy, sinh cảnh sống tự nhiên của loài đang bị thu hẹp và chia cắt nghiêm trọng. Tổ chức Bảo tồn Thiên nhiên Quốc tế (IUCN, 2014) đã xếp Bách xanh vào cấp độ Sẽ nguy cấp (VU A2cd), đồng thời pháp luật Việt Nam đưa loài vào Nhóm IIA trong Danh mục thực vật rừng nguy cấp, quý, hiếm nhằm hạn chế khai thác.

Vấn đề đặt ra là dù Bách xanh có khả năng tái sinh bằng hạt tốt, nhưng tỷ lệ cây con sống sót phát triển thành cây trưởng thành trong tự nhiên rất thấp. Trong khi đó, các nghiên cứu trước đây chủ yếu tập trung vào đặc điểm hình thái và vùng phân bố địa lý, còn dữ liệu về đa dạng di truyền nguồn gen ở cấp độ phân tử tại khu vực Tây Nguyên gần như chưa được khai phá.

Mục tiêu trọng tâm của nghiên cứu là đánh giá mức độ đa dạng di truyền và cấu trúc nguồn gen của 3 quần thể Bách xanh tự nhiên tại Lâm Đồng, Đắk Lắk và Gia Lai bằng hai hệ thống chỉ thị phân tử ISSR và SSR. Nghiên cứu thuộc khuôn khổ Đề tài TN3/T15 thuộc Chương trình Tây Nguyên 3, thực hiện tại Bảo tàng Thiên nhiên Việt Nam và Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội. Ý nghĩa của đề tài thể hiện qua việc định lượng chính xác các chỉ số di truyền quần thể (tỷ lệ phân đoạn đa hình đạt từ 50,39% đến 95,24%, hệ số gen dị hợp tử mong đợi He đạt từ 0,168 đến 0,367), tạo tiền đề khoa học vững chắc cho công tác quy hoạch bảo tồn nguyên vị và phục hồi nguồn gen loài cây lá kim đặc hữu này.

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ể (Population Genetics Theory) và sinh học bảo tồn thực vật. Theo nguyên lý cân bằng Hardy-Weinberg và động lực học quần thể ngẫu phối, sự đa dạng di truyền là thước đo khả năng thích nghi và tồn tại lâu dài của loài trước biến đổi môi trường. Quần thể giao phối chéo đóng vai trò là đơn vị sinh sản cơ bản; khi quy mô quần thể bị thu hẹp và cô lập địa lý, hiện tượng giao phối cận noãn (inbreeding) cùng sự trôi dạt di truyền (genetic drift) sẽ làm gia tăng tần số đồng hợp tử, dẫn đến suy thoái nguồn gen.

Các khái niệm chính được vận dụng trong nghiên cứu bao gồm:

  • Đa dạng di truyền (Genetic Diversity): Mức độ biến dị của các tổ hợp nucleotide trong hệ gen giữa các cá thể và giữa các quần thể.
  • Hệ số dị hợp tử quan sát (Ho) và mong đợi (He): Chỉ số phản ánh mức độ phong phú alen và cấu trúc kiểu gen của quần thể theo lý thuyết Nei (1973).
  • Hàm lượng thông tin đa hình (Polymorphism Information Content - PIC): Thước đo năng lực phân biệt kiểu gen của các locus chỉ thị phân tử.
  • Phân tích biến đổi phân tử (AMOVA): Phương pháp thống kê phi tham số phân tách phương sai di truyền giữa các quần thể và trong nội bộ quần thể.

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

Nghiên cứu tiến hành thu thập tổng cộng 70 mẫu lá và mô gỗ từ 3 quần thể Bách xanh tự nhiên đại diện tại khu vực Tây Nguyên:

  • Quần thể Đatanla (thành phố Đà Lạt, tỉnh Lâm Đồng): 32 mẫu (ký hiệu Cm1 - Cm32), thu ở tọa độ 11°54'2.8" N và độ cao 1315 m.
  • Quần thể Hòa Sơn (huyện Krông Bông, tỉnh Đắk Lắk): 34 mẫu (ký hiệu Cm33 - Cm66), thu ở tọa độ 12°25'05" N - 108°22'17" E và độ cao 1200 m.
  • Quần thể Kon Chư Răng (huyện K'Bang, tỉnh Gia Lai): 4 mẫu (ký hiệu Cm67 - Cm70), thu ở tọa độ 14°30'52" N - 108°33'21" E và độ cao từ 1040 m đến 1057 m.

Phương pháp chọn mẫu đảm bảo tính đại diện sinh cảnh tự nhiên theo từng tầng tán rừng. Toàn bộ DNA tổng số được tách chiết theo quy trình CTAB cải tiến của Porebski và cộng sự (1997), đạt nồng độ tinh sạch cao từ 640 đến 1000 ng/µl với tỷ lệ hấp thụ quang phổ OD260/OD280 dao động trong khoảng 1,732 đến 2,000.

Lý do lựa chọn kết hợp hai hệ thống chỉ thị phân tử ISSR (Inter Simple Sequence Repeat) và SSR (Simple Sequence Repeat): Chỉ thị ISSR với 30 mồi đơn giúp quét rộng toàn bộ hệ gen mà không cần biết trước trình tự nucleotide; trong khi 17 cặp mồi SSR là chỉ thị đồng trội locus đặc hiệu, cho phép phát hiện chính xác trạng thái dị hợp tử và các alen hiếm.

Số liệu điện di gel agarose 1,5% (ISSR) và polyacrylamide 5% (SSR) được mã hóa nhị phân và phân tích ma trận tương đồng di truyền qua các phần mềm chuyên dụng: NTSYS 2.0 (phân tích UPGMA, PCA), GENALEX 6.5 (tính toán các chỉ số Na, Ne, Ho, He, Fst, Nm, AMOVA) và kiểm định Bootstrap 1000 lần lặp lại qua phần mềm Win-Boot.

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

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

Quá trình phân tích thực nghiệm trên 70 mẫu Bách xanh mang lại các phát hiện cốt lõi:

Thứ nhất, tính đa hình DNA thể hiện rõ rệt qua cả hai loại chỉ thị phân tử. Với 30 mồi ISSR, có 25 mồi cho phản ứng đa hình, nhân bản được 129 phân đoạn DNA (kích thước 250 - 2000 bp), trong đó có 65 phân đoạn đa hình, đạt tỷ lệ 50,39% với chỉ số PIC trung bình 0,101. Với 17 cặp mồi SSR, có 15 cặp mồi đa hình, nhân bản 42 phân đoạn (kích thước 100 - 430 bp), trong đó có 40 phân đoạn đa hình, đạt tỷ lệ vượt trội 95,24% cùng giá trị PIC trung bình 0,272.

Thứ hai, mức độ đa dạng di truyền có sự phân hóa mạnh mẽ giữa các quần thể. Theo chỉ thị ISSR, quần thể Đatanla đạt mức đa dạng cao nhất với tỷ lệ phân đoạn đa hình PPB = 35,66%, chỉ số Shannon I = 0,192 và He = 0,130; kế tiếp là Hòa Sơn với PPB = 34,88%, I = 0,175 và He = 0,117; thấp nhất là Kon Chư Răng với PPB chỉ đạt 3,88%, I = 0,022 và He = 0,015. Kết quả phân tích SSR khẳng định xu hướng tương tự: Đatanla dẫn đầu với PPB = 88,24%, I = 0,496 và He = 0,321; Hòa Sơn đạt PPB = 76,47%, I = 0,466 và He = 0,320; Kon Chư Răng thấp nhất với PPB = 47,06%, I = 0,313 và He = 0,211. Ở cấp độ loài, các chỉ số đạt Na = 2,471, Ne = 1,733, Ho = 0,339 và He = 0,367.

Thứ ba, cấu trúc biến dị phân tử duy trì chủ yếu trong nội bộ quần thể. Phân tích AMOVA với chỉ thị ISSR cho thấy 63,67% tổng phương sai biến đổi di truyền xuất phát từ sự khác biệt giữa các cá thể trong cùng một quần thể, trong khi biến đổi giữa các quần thể chỉ chiếm 36,33% (giá trị thống kê p < 0,001).

Thứ tư, mối quan hệ di truyền phản ánh chính xác ranh giới địa lý. Hệ số tương đồng di truyền Nei (1973) giữa các quần thể đạt mức cao từ 0,861 (Đatanla - Kon Chư Răng) đến 0,955 (Đatanla - Hòa Sơn), trung bình là 0,900; khoảng cách di truyền dao động từ 0,046 đến 0,150.

Thảo luận kết quả

Sự suy giảm đa dạng di truyền nghiêm trọng tại quần thể Kon Chư Răng bắt nguồn trực tiếp từ kích thước quần thể quá nhỏ (chỉ tìm thấy 4 cá thể trưởng thành cao từ 15 đến 20 m), dẫn đến hiện tượng nghẽn cổ chai di truyền và thiếu hụt dòng gen. Ngược lại, quần thể Đatanla sở hữu hơn 10 cây trưởng thành cao trên 10 m cùng lớp cây tái sinh phong phú, tạo điều kiện thuận lợi cho quá trình thụ phấn chéo tự nhiên, duy trì hệ số giao phấn cận noãn âm (Fis = -0,025 ở Đatanla và -0,124 tính trung bình toàn loài).

Mối quan hệ phát sinh giữa 70 cá thể được thể hiện trực quan qua biểu đồ hình cây UPGMA và biểu đồ tọa độ không gian 2D/3D (PCA). Trên đồ thị phân nhóm, 70 mẫu tách biệt thành 2 nhánh chính: Nhánh I quy tụ toàn bộ 4 mẫu của Kon Chư Răng (Cm67 - Cm70) với độ tương đồng nội nhóm từ 96,5% đến 98,0%. Nhánh II tập hợp 66 mẫu của Đatanla và Hòa Sơn với độ tương đồng từ 85,5% đến 99,1%. Biểu đồ điện di gel polyacrylamide của các mồi SSR như Cm3 và Cm7 cho thấy rõ sự hiện diện của các băng alen đặc thù (vị trí 415 bp ở mồi Cm3 chỉ xuất hiện tại Gia Lai và một số ít cá thể Lâm Đồng), giải thích vì sao quần thể Gia Lai lại phân nhánh tách biệt.

Khi đối chiếu với các nghiên cứu quốc tế, mức độ đa dạng di truyền của Bách xanh Tây Nguyên (ISSR: He = 0,168; SSR: He = 0,367) cao hơn quần thể Bách xanh tại Tây Nam Trung Quốc (He = 0,1114), nhưng thấp hơn đáng kể so với các loài thông phân bố rộng tại châu Âu như Pinus sylvestris (He = 0,262; PPB = 99,76%). Điều này chứng minh nguồn gen Bách xanh tự nhiên tại Tây Nguyên đang đứng trước nguy cơ xói mòn di truyền cao nếu không có biện pháp can thiệp kịp thời.

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

Dựa trên các chỉ số di truyền định lượng được, đề tài đưa ra 4 nhóm giải pháp cấp bách nhằm bảo tồn và phát triển nguồn gen Bách xanh:

  • Thiết lập phân khu bảo tồn nghiêm ngặt tại chỗ (In situ): Ban Quản lý Vườn Quốc gia Chư Yang Sin và Khu Bảo tồn Thiên nhiên Kon Chư Răng cần tiến hành cắm mốc ranh giới, lập hồ sơ quản lý định danh cho 70 cá thể Bách xanh đã nghiên cứu, đặc biệt là bảo vệ tuyệt đối 4 cây mẹ tại Gia Lai và hơn 50 cây tại khu vực suối Đatanla, ngăn chặn triệt để 100% các hành vi chặt hạ, khai thác trong giai đoạn 2026 - 2028.
  • Thực hiện nhân giống và xây dựng vườn bảo tồn chuyển vị (Ex situ): Sở Nông nghiệp và Phát triển Nông thôn các tỉnh Tây Nguyên phối hợp với Xí nghiệp Giống Lâm nghiệp vùng Tây Nguyên thu hái hạt giống từ 32 cây mẹ ưu trội tại Đatanla và 34 cây tại Hòa Sơn để gieo ươm 5.000 cây giống đạt tiêu chuẩn xuất vườn trong thời hạn 24 tháng, phục vụ trồng bổ sung vào rừng nguyên sinh.
  • Thúc đẩy dòng gen nhân tạo giữa các quần thể cô lập: Viện Hàn lâm Khoa học và Công nghệ Việt Nam cùng các đơn vị lâm nghiệp địa phương triển khai kỹ thuật hỗ trợ thụ phấn chéo nhân tạo, chuyển giao nguồn hạt phấn giữa quần thể Đatanla và Kon Chư Răng nhằm gia tăng hệ số di nhập gen (hiện đạt Nm = 4,176), phá vỡ tình trạng suy thoái do nội phối trong chu kỳ 3 - 5 năm tới.
  • Ứng dụng chỉ thị phân tử giám sát định kỳ chất lượng cây con: Trung tâm Giống lâm nghiệp ứng dụng bộ 15 cặp mồi SSR đa hình cao (như Cm3, Cm7, Cm12, RPS1b) để sàng lọc di truyền định kỳ 3 năm/lần cho các lô cây con nhân tạo, cam kết duy trì hệ số dị hợp tử mong đợi He đạt trên mức 0,350 trước khi đưa vào trồng tái thiết rừng.

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

Nội dung và kết quả của luận văn mang giá trị ứng dụng thiết thực cho nhiều nhóm đối tượng chuyên môn:

  • Cán bộ quản lý vườn quốc gia và lực lượng kiểm lâm: Sử dụng dữ liệu định vị tọa độ, độ cao và mức độ nguy cấp di truyền của 3 quần thể để thiết kế phương án tuần tra, bảo vệ nghiêm ngặt các cá thể cây mẹ hạt trần quý hiếm tại Lâm Đồng, Đắk Lắk và Gia Lai.
  • Nhà khoa học, giảng viên và học viên chuyên ngành Di truyền học, Lâm nghiệp, Công nghệ sinh học: Kế thừa quy trình tách chiết DNA CTAB cải tiến và bộ thông số mồi chuẩn (30 mồi ISSR, 17 cặp mồi SSR) phục vụ công tác giảng dạy, nghiên cứu đa dạng sinh học phân tử trên các đối tượng thực vật có hàm lượng hợp chất thứ sinh cao.
  • Kỹ sư chọn tạo giống và chuyên viên lâm nghiệp: Khai thác dữ liệu về hệ số tương đồng di truyền (0,861 - 0,955) để tuyển chọn các cặp bố mẹ có khoảng cách di truyền xa, tối ưu hóa hiệu quả nhân giống sinh dưỡng và tránh thoái hóa giống do giao phối cận huyết.
  • Cơ quan quản lý nhà nước và tổ chức bảo tồn đa dạng sinh học: Tham khảo các số liệu định lượng (He = 0,367; PPB = 95,24%) để cập nhật Danh lục Đỏ, Sách Đỏ Việt Nam và xây dựng đề án phân bổ ngân sách cho các chương trình bảo tồn nguồn gen quốc gia.

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

Luận văn sử dụng những hệ thống chỉ thị phân tử nào để đánh giá đa dạng di truyền Bách xanh?
Nghiên cứu kết hợp hai loại chỉ thị phân tử là 30 mồi ISSR và 17 cặp mồi SSR. Sự kết hợp này mang tính toàn diện vì ISSR giúp quét đa hình trên toàn hệ gen mà không cần biết trước trình tự, còn SSR là chỉ thị đồng trội giúp định lượng chính xác các alen hiếm và hệ số dị hợp tử.

Mức độ đa dạng di truyền của loài Bách xanh ở Tây Nguyên được đánh giá như thế nào?
Mức độ đa dạng di truyền của loài Bách xanh ở Tây Nguyên nằm ở mức trung bình thấp (với chỉ thị ISSR đạt He = 0,168; với chỉ thị SSR đạt He = 0,367 và I = 0,602). Mức này cao hơn Bách xanh ở Tây Nam Trung Quốc nhưng thấp hơn nhiều so với các loài lá kim phân bố tự nhiên ở châu Âu.

Tại sao quần thể Bách xanh ở Kon Chư Răng lại có chỉ số đa dạng di truyền thấp nhất?
Quần thể Kon Chư Răng chỉ thu thập được 4 cá thể trưởng thành, sinh cảnh bị chia cắt mạnh khiến việc giao phấn ngẫu nhiên bị hạn chế nghiêm trọng. Điều này dẫn tới tỷ lệ phân đoạn đa hình ISSR chỉ đạt 3,88% và SSR đạt 47,06%, phản ánh hiện tượng xói mòn nguồn gen do quy mô quần thể quá nhỏ.

Quy trình tách chiết DNA trong nghiên cứu đạt tiêu chuẩn chất lượng ra sao?
DNA tổng số tách chiết từ lá và gỗ theo phương pháp Porebski cải tiến cho các băng điện di đậm, sắc nét, không bị đứt gãy. Nồng độ DNA đo quang phổ dao động từ 640 đến 1000 ng/µl với tỷ lệ OD260/OD280 đạt từ 1,732 đến 2,000, hoàn toàn đáp ứng độ tinh sạch chuẩn cho phản ứng PCR.

Ý nghĩa của kết quả phân tích biến đổi phân tử (AMOVA) trong nghiên cứu là gì?
Phân tích AMOVA với chỉ thị ISSR chỉ ra rằng 63,67% biến đổi di truyền tập trung trong nội bộ quần thể và 36,33% diễn ra giữa các quần thể (p < 0,001). Kết quả này khẳng định chiến lược bảo tồn cần ưu tiên giữ gìn sự đa dạng giữa từng cá thể trong cùng một vùng sinh cảnh.

Kết luận

  • Hệ thống chỉ thị phân tử chứng minh hiệu quả cao với 65/129 phân đoạn ISSR (50,39%) và 40/42 phân đoạn SSR (95,24%) thể hiện tính đa hình rõ rệt trên 70 mẫu nghiên cứu.
  • Xác định mức độ đa dạng di truyền tổng thể của loài Bách xanh tại Tây Nguyên ở mức trung bình thấp (ISSR: He = 0,168; SSR: He = 0,367; PPB = 51,16% - 88,24%).
  • Khẳng định sự phân hóa di truyền sâu sắc theo vùng: Quần thể Đatanla sở hữu độ đa dạng cao nhất (SSR: I = 0,496), trong khi quần thể Kon Chư Răng suy thoái nghiêm trọng (SSR: I = 0,313; ISSR: PPB = 3,88%).
  • Chứng minh phương sai biến đổi di truyền chủ yếu duy trì bên trong nội bộ quần thể (chiếm 63,67%) với hệ số tương đồng di truyền trung bình giữa các quần thể đạt 0,900.
  • Xây dựng thành công sơ đồ phát sinh chủng loại UPGMA phân lập 70 cá thể thành 2 nhóm cấu trúc di truyền tương ứng với khoảng cách phân bố địa lý.

Đóng góp lớn nhất của luận văn là cung cấp bộ dữ liệu sinh học phân tử có hệ thống đầu tiên về quần thể Bách xanh Tây Nguyên thuộc Chương trình Tây Nguyên 3. Về kế hoạch tiếp theo trong giai đoạn 3 - 5 năm tới, các cơ quan lâm nghiệp cần khẩn trương mở rộng giải trình tự hệ gen lục lạp và áp dụng kỹ thuật thụ phấn nhân tạo để phục hồi nguồn gen. Hãy liên hệ với các đơn vị bảo tồn và viện nghiên cứu lâm nghiệp chuyên trách để chung tay hành động bảo vệ loài Bách xanh quý hiếm trước nguy cơ tuyệt chủng.