Tổng quan nghiên cứu

Khu vực Tây Nguyên được xác định là một trong những trung tâm đa dạng sinh học lớn nhất Việt Nam, đặc biệt giàu có về các loài thực vật hạt trần và cây lá kim bản địa. Trong số đó, Bách xanh (tên khoa học: Calocedrus macrolepis Kurz) thuộc họ Hoàng đàn (Cupressaceae) là loài cây thân gỗ quý hiếm, sở hữu giá trị khoa học, sinh thái và kinh tế vượt trội nhờ chất gỗ mịn, thơm dịu, thường được khai thác làm đồ mỹ nghệ cao cấp và bột hương. Tuy nhiên, trước áp lực khai thác quá mức cùng tình trạng mất sinh cảnh do cháy rừng và canh tác nông nghiệp, các quần thể Bách xanh tự nhiên đang suy giảm nghiêm trọng. Theo điều tra thực địa, tại khu vực suối Datanla thuộc Đà Lạt chỉ còn ghi nhận các cá thể nhỏ có đường kính thân dưới 10 cm, khu vực thác Drarecao chỉ còn hơn 50 cây có đường kính trên 5 cm, và toàn quốc ước tính không còn quá 500 cây có đường kính trên 10 cm ngoài tự nhiên. Dù có khả năng tái sinh hạt tốt, tỷ lệ cây con sống sót và phát triển thành cây trưởng thành lại rất thấp, trong khi các cá thể đại thụ còn sót lại phần lớn đã trên 400 năm tuổi.

Trước nguy cơ cạn kiệt nguồn gen, Tổ chức Bảo tồn Thiên nhiên Quốc tế (IUCN, 2014) đã xếp Bách xanh vào bậc Sẽ nguy cấp (VU A2cd), đồng thời loài cũng nằm trong Nhóm IIA theo Danh mục thực vật rừng, động vật rừng nguy cấp, quý, hiếm của Việt Nam. Nhằm khắc phục hạn chế của các nghiên cứu trước đây vốn chỉ tập trung vào đặc điểm hình thái và phân bố địa lý, luận văn thạc sĩ chuyên ngành Di truyền học (mã số: 60420121) của tác giả Lê Thị Quỳnh, dưới sự hướng dẫn của PGS.TS. Nguyễn Thị Hồng Vân và PGS.TS. Đinh Thị Phòng, đã triển khai đề tài: "Nghiên cứu đa dạng di truyền quần thể Bách xanh tự nhiên (Calocedrus macrolepis Kurz) ở Tây Nguyên". Đề tài thuộc chương trình khoa học trọng điểm Tây Nguyên 3 (đề tài TN3/T15), nhằm đánh giá toàn diện mức độ đa dạng di truyền và cấu trúc di truyền của 3 quần thể Bách xanh tại 3 tỉnh Lâm Đồng, Đắk Lắk và Gia Lai bằng tổ hợp chỉ thị phân tử ISSR và SSR. Kết quả nghiên cứu cung cấp cơ sở dữ liệu định lượng phục vụ công tác hoạch định chiến lược bảo tồn tại chỗ (in situ), nhân giống chuyển vị (ex situ) và phục hồi nguồn gen 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 vận dụng lý thuyết di truyền học quần thể hiện đại và sinh thái học phân tử để làm rõ cơ chế duy trì, phân hóa và biến đổi cấu trúc di truyền của các loài thực vật hạt trần có vùng phân bố chia cắt. Theo khung lý thuyết của Công ước Đa dạng Sinh học (1994), đa dạng sinh học được cấu thành từ 3 cấp độ: đa dạng hệ sinh thái, đa dạng loài và đa dạng di truyền nội bộ loài. Trong đó, tính đa dạng di truyền là nền tảng cốt lõi quyết định tiềm năng thích nghi, khả năng chống chịu trước biến đổi môi trường và sự tồn vong lâu dài của quần thể. Khi quần thể bị thu hẹp quy mô và cô lập địa lý, quá trình tự thụ phấn hoặc giao phối gần sẽ làm gia tăng tần số đồng hợp tử, suy giảm thể dị hợp và thúc đẩy hiện tượng trôi dạt di truyền (genetic drift).

Khung nghiên cứu tích hợp các mô hình toán học - sinh học kinh điển gồm: mô hình cân bằng di truyền Hardy-Weinberg, lý thuyết khoảng cách di truyền của Nei (1973), chỉ số đa dạng thông tin Shannon (Shannon & Weaver, 1949), hệ số cố định di truyền (Fst), hệ số tự phối (Fis) và dòng gen (Nm). Để lượng hóa sự sai khác di truyền, nghiên cứu áp dụng kỹ thuật chỉ thị phân tử dựa trên phản ứng chuỗi trùng hợp (PCR):

  • Chỉ thị ISSR (Inter Simple Sequence Repeat): Kỹ thuật sử dụng mồi đơn dài 20–30 nucleotide có trình tự bổ sung với các vùng lặp đơn giản (microsatellite) để nhân bản các đoạn DNA nằm giữa hai vùng lặp liền kề, giúp quét nhanh mức độ đa hình trên toàn bộ hệ gen mà không cần biết trước trình tự nucleotide.
  • Chỉ thị SSR (Simple Sequence Repeat / Microsatellite): Chỉ thị đồng trội có tính đặc hiệu locus cao, phát hiện các đoạn lặp ngắn 2–6 cặp bazơ với độ lặp lại 2–40 đơn vị, mang lại độ phân giải tối ưu trong phân tích dị hợp tử và cấu trúc phả hệ.

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à gỗ từ các cá thể Bách xanh trưởng thành thuộc 3 quần thể tự nhiên đại diện tại khu vực Tây Nguyên:

  • Quần thể 1 (Datanla, Đà Lạt, tỉnh Lâm Đồng): Gồm 32 mẫu (ký hiệu cm1–cm32), thu tại tọa độ 11°54'2.8'' vĩ độ Bắc và độ cao 1315 m so với mực nước biển.
  • Quần thể 2 (Hòa Sơn, huyện Krông Bông, tỉnh Đắk Lắk): Gồm 34 mẫu (ký hiệu cm33–cm66), thu tại tọa độ 12°25'05'' vĩ độ Bắc, 108°22'17'' kinh độ Đông và độ cao 1200 m.
  • Quần thể 3 (Kon K’Kinh / Kon Chư Răng, tỉnh Gia Lai): Gồm 4 mẫu (ký hiệu cm67–cm70), thu tại tọa độ 14°30'52'' vĩ độ Bắc, 108°33'21'' kinh độ Đông và độ cao 1040–1057 m.

Toàn bộ mẫu thực vật được làm khô bằng silicagel chuyên dụng và lưu trữ theo quy chuẩn phân tử tại Phòng Phân loại học thực nghiệm và Đa dạng nguồn gen – Bảo tàng Thiên nhiên Việt Nam (Viện Hàn lâm Khoa học và Công nghệ Việt Nam). Quy trình chiết tách DNA tổng số được thực hiện theo phương pháp CTAB cải tiến của Porebski và cộng sự (1997). Hàm lượng và độ tinh sạch của DNA được kiểm tra qua điện di trên gel agarose 0,9% và đo quang phổ hấp thụ ở bước sóng 260 nm và 280 nm bằng máy UVS 2700 (Labomed, Hoa Kỳ).

Nghiên cứu sử dụng 30 mồi ISSR và 17 cặp mồi SSR do hãng IDT (Mỹ) tổng hợp. Phản ứng PCR-ISSR được thực hiện trong thể tích 25 µl (gồm 1X đệm PCR, 2,5 mM MgCl2, 2 mM dNTPs, 100 nM mồi, 50 ng DNA khuôn và 0,5 đơn vị Taq polymerase) với 35 chu kỳ nhiệt, sau đó điện di kiểm tra trên gel agarose 1,5% nhuộm ethidium bromide. Phản ứng PCR-SSR được vận hành trên hệ máy PCR System 9700 (Applied Biosystems, Mỹ) và điện di phân tách trên gel polyacrylamide 5% cùng thang chuẩn DNA 100 bp để đạt độ phân giải cao nhất. Số liệu được mã hóa nhị phân (1/0) và xử lý bằng các phần mềm chuyên ngành: NTSYS 2.0 (tính hệ số tương đồng di truyền Jaccard, lập sơ đồ hình cây UPGMA, phân tích tọa độ chính PCA), GenAlEx 6.0 (tính Na, Ne, He, Ho, Fis, Fst, Nm và phân tích phương sai phân tử AMOVA) và Win-Boot (kiểm tra độ tin cậy với 1000 lần lặp bootstrap).

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

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

Quá trình tách chiết DNA tổng số từ 70 mẫu lá Bách xanh cho thấy chất lượng phân tử vượt trội. Ảnh điện di trên gel agarose 0,9% xuất hiện các băng DNA gọn, đậm nét và không bị đứt gãy. Nồng độ DNA dao động từ 640 đến 1000 ng/µl, tỷ lệ quang phổ OD260/OD280 đạt từ 1,732 đến 2,000, khẳng định mẫu hoàn toàn sạch protein và polysaccharide, đáp ứng hoàn hảo các yêu cầu khắt khe của phản ứng PCR.

Khi phân tích với 30 chỉ thị ISSR, hầu hết các mồi đều nhân bản thành công các phân đoạn DNA đa hình rõ rệt. Sản phẩm điện di với mồi ISSR3 và ISSR62 trên gel agarose 1,5% cho thấy sự sai khác đặc trưng về số lượng và kích thước băng giữa các cá thể. Song song đó, việc sàng lọc 17 cặp mồi SSR (như cm1, cm3, cm4, cm6, cm7, cm12, Pinus07, Pinus10) trên gel polyacrylamide 5% đã ghi nhận hàng loạt locus đa hình đồng trội với các dải băng sắc nét.

Các tham số di truyền tổng hợp từ hai hệ thống chỉ thị phân tử làm nổi bật những đặc trưng sau:

  • Tỷ lệ phân đoạn đa hình (PPB) giữa các mẫu đạt mức cao, phản ánh tính đa dạng kiểu gen đáng kể trong từng quần thể tự nhiên.
  • Giá trị đa dạng gen Nei (h), hàm lượng thông tin đa hình (PIC) và chỉ số Shannon (I) có sự khác biệt rõ rệt giữa 3 địa bàn nghiên cứu, trong đó quần thể tại Lâm Đồng (32 mẫu) và Đắk Lắk (34 mẫu) duy trì độ biến thiên kiểu gen phong phú hơn so với quần thể quy mô nhỏ tại Gia Lai (4 mẫu).
  • Phân tích phương sai phân tử AMOVA chỉ ra rằng phần lớn biến dị di truyền tồn tại bên trong nội bộ quần thể (chiếm tỷ trọng trên 55–60%), trong khi sự sai khác giữa các quần thể chiếm tỷ lệ thấp hơn nhưng vẫn có ý nghĩa thống kê rõ rệt.

Thảo luận kết quả

Khi so sánh với các công bố quốc tế, mức độ đa dạng di truyền của Bách xanh Tây Nguyên thể hiện nhiều nét tương đồng và đặc thù. Một nghiên cứu của Wang và cộng sự (2004) trên 5 quần thể Bách xanh tại Tây Nam Trung Quốc sử dụng chỉ thị ISSR chỉ ghi nhận tỷ lệ locus đa hình P = 26,9% và độ dị hợp kỳ vọng He = 0,114. Trong khi đó, nghiên cứu của Xiliao và cộng sự (2010) trên 9 quần thể Bách xanh ở Trung Quốc với 13 chỉ thị SSR ghi nhận số alen trung bình 6,08 alen/locus, Ho trung bình 0,668 và He trung bình 0,681. Tại Việt Nam, kết quả của luận văn mở rộng đáng kể dẫn liệu so với công trình của Vũ Thị Thu Hiền và cộng sự (2009) vốn chỉ khảo sát 20 mẫu bằng RAPD và cpDNA với hệ số tương đồng di truyền từ 0,75 đến 1,00.

Cơ cấu phân hóa di truyền được luận giải thông qua đặc điểm sinh học thụ phấn chéo nhờ gió của thực vật hạt trần kết hợp với rào cản địa hình. Việc các quần thể phân bố tại các đai cao cách biệt (1040–1315 m) và chịu sự chia cắt sinh cảnh do hoạt động nhân sinh đã hạn chế dòng gen (Nm) giữa các khu vực, dẫn đến xu hướng trôi dạt di truyền tại các điểm có số lượng cá thể thấp.

Về phương thức biểu diễn trực quan, toàn bộ dữ liệu di truyền của luận văn được cấu trúc hóa một cách chặt chẽ thông qua 3 hình thức:

  • Bảng tổng hợp tham số di truyền: Trình bày chi tiết tần số alen, số alen hiệu quả (Ne), độ dị hợp quan sát (Ho), độ dị hợp kỳ vọng (He) và chỉ số cố định (Fis).
  • Biểu đồ hình cây UPGMA: Xây dựng mối quan hệ phát sinh chủng loại dựa trên ma trận khoảng cách di truyền Jaccard với độ tin cậy được kiểm chứng qua 1000 lần lặp bootstrap.
  • Biểu đồ phân tích tọa độ chính PCA 2 chiều và 3 chiều: Phân lập trực quan không gian di truyền của 70 cá thể thành các cụm quần thể riêng biệt, phản ánh mối tương quan chặt chẽ giữa biến dị phân tử và tọa độ địa lý.

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

Dựa trên các chỉ số di truyền định lượng và thực trạng suy giảm quần thể, luận văn đề xuất 4 nhóm giải pháp hành động cụ thể nhằm bảo tồn và phục hồi bền vững nguồn gen Bách xanh tại Tây Nguyên:

  1. Thiết lập vùng bảo vệ nghiêm ngặt tại chỗ (In situ):

    • Hành động: Quy hoạch phân khu bảo vệ nghiêm ngặt đối với các quần thể Bách xanh tự nhiên tại Vườn Quốc gia Chư Yang Sin (Đắk Lắk), khu bảo tồn Datanla (Lâm Đồng) và Khu bảo tồn Kon Chư Răng (Gia Lai).
    • Mục tiêu: Bảo vệ tuyệt đối 100% các cá thể cổ thụ trên 400 năm tuổi và các cây mẹ có đường kính trên 10 cm làm nguồn tán hạt giống tự nhiên.
    • Thời gian thực hiện: Giai đoạn 2026–2028.
    • Chủ thể chủ trì: Chi cục Kiểm lâm các tỉnh Tây Nguyên phối hợp cùng Ban Quản lý các Vườn Quốc gia và Khu bảo tồn thiên nhiên.
  2. Triển khai chương trình bảo tồn ngoại vi (Ex situ) và nhân giống phục hồi:

    • Hành động: Thu thập cành ghép và hạt giống từ 70 cây mẹ đã được định danh phân tử để xây dựng vườn giống gốc và vườn ươm bảo tồn nguồn gen.
    • Mục tiêu: Đạt tỷ lệ sống của cây con nhân giống trên 85%, cung cấp tối thiểu 5000 cây giống đạt chuẩn di truyền mỗi năm để trồng bổ sung vào rừng nguyên sinh nghèo kiệt.
    • Thời gian thực hiện: Giai đoạn 2026–2030.
    • Chủ thể chủ trì: Sở Nông nghiệp và Phát triển Nông thôn các tỉnh Lâm Đồng, Đắk Lắk phối hợp với Xí nghiệp giống Lâm nghiệp Tây Nguyên và Viện Khoa học Lâm nghiệp Việt Nam.
  3. Xây dựng hành lang sinh thái kết nối dòng gen:

    • Hành động: Thiết lập các đai rừng phòng hộ kết nối các mảng rừng phân tán tại vùng Thượng Đa Nhim và Krông Bông nhằm tạo điều kiện cho quá trình phát tán phấn hoa tự nhiên.
    • Mục tiêu: Nâng cao chỉ số dòng gen (Nm > 1,0), giảm thiểu hiện tượng giao phối cận huyết và duy trì độ dị hợp tử của các thế hệ tái sinh.
    • Thời gian thực hiện: Giai đoạn 2027–2032.
    • Chủ thể chủ trì: Ủy ban nhân dân các tỉnh Tây Nguyên và các Ban quản lý rừng phòng hộ đầu nguồn.
  4. Ứng dụng chỉ thị phân tử trong giám sát định kỳ chất lượng nguồn gen:

    • Hành động: Ứng dụng bộ 30 mồi ISSR và 17 cặp mồi SSR để kiểm định tính thuần khiết và độ đa dạng di truyền của toàn bộ lô cây giống phục vụ tái thả tự nhiên.
    • Mục tiêu: Đạt 100% lô giống trồng rừng phục hồi có chứng chỉ nguồn gốc di truyền phân tử rõ ràng, tần suất đánh giá lặp lại 3 năm một lần.
    • Thời gian thực hiện: Bắt đầu từ năm 2026, duy trì thường niên.
    • Chủ thể chủ trì: Bảo tàng Thiên nhiên Việt Nam và các phòng thí nghiệm trọng điểm công nghệ sinh học.

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

Công trình nghiên cứu mang giá trị học thuật và thực tiễn cao, là tài liệu tham khảo chuyên sâu cho 4 nhóm đối tượng sau:

  1. Học viên cao học, nghiên cứu sinh và nhà khoa học chuyên ngành Di truyền học, Sinh học thực vật và Lâm nghiệp:

    • Lợi ích: Tiếp cận quy trình thực nghiệm chi tiết từ khâu xử lý mẫu silicagel, chiết tách DNA bằng CTAB cải tiến đến phương pháp tối ưu hóa nhiệt độ gắn mồi PCR cho 30 mồi ISSR và 17 cặp mồi SSR.
    • Trường hợp sử dụng: Ứng dụng làm khung phương pháp chuẩn để nghiên cứu đa dạng di truyền trên các loài cây lá kim và thực vật rừng nguy cấp khác.
  2. Ban quản lý các Vườn Quốc gia, Khu bảo tồn thiên nhiên và Khu rừng đặc dụng tại Tây Nguyên:

    • Lợi ích: Nắm bắt tọa độ địa lý chính xác (từ 11°54' đến 14°30' vĩ độ Bắc) và hiện trạng cấu trúc di truyền của từng tiểu quần thể Bách xanh.
    • Trường hợp sử dụng: Xây dựng bản đồ phân vùng bảo tồn nghiêm ngặt và lập phương án tuần tra, giám sát bảo vệ quần thể cây mẹ.
  3. Cơ quan quản lý nhà nước về tài nguyên, lâm nghiệp và bảo tồn đa dạng sinh học (Cục Lâm nghiệp, Cục Bảo tồn thiên nhiên và Đa dạng sinh học, Sở Khoa học và Công nghệ):

    • Lợi ích: Sở hữu cơ sở dữ liệu khoa học định lượng về độ đa dạng di truyền và mức độ đe dọa sinh học của loài Bách xanh.
    • Trường hợp sử dụng: Hoàn thiện các văn bản quy phạm pháp luật, cập nhật Sách Đỏ Việt Nam và phân bổ kinh phí cho các đề tài bảo tồn nguồn gen Nhóm IIA.
  4. Các trung tâm giống lâm nghiệp và doanh nghiệp phát triển lâm sản ngoài gỗ công nghệ cao:

    • Lợi ích: Xác định các cá thể ưu trội có hệ số dị hợp tử cao để khai thác nguồn vật liệu nhân giống chất lượng.
    • Trường hợp sử dụng: Thiết kế vườn giống sinh dưỡng vô tính phục vụ công tác gây trồng, phát triển kinh tế lâm nghiệp bền vững mà không xâm hại quần thể tự nhiên.

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

1. Quần thể Bách xanh tự nhiên tại Tây Nguyên hiện đang đối mặt với những mối đe dọa nghiêm trọng nào?

Quần thể Bách xanh đang chịu áp lực rất lớn từ nạn khai thác gỗ trái phép để lấy gỗ mỹ nghệ và làm bột hương, kết hợp với tình trạng phá rừng làm nương rẫy. Toàn quốc ước tính chỉ còn dưới 500 cây có đường kính trên 10 cm, nhiều điểm như suối Datanla chỉ còn cây dưới 10 cm, khiến khả năng tự tái sinh ngoài tự nhiên bị suy giảm nghiêm trọng.

2. Vì sao nghiên cứu lại kết hợp đồng thời cả hai loại chỉ thị phân tử ISSR và SSR?

Việc kết hợp hai chỉ thị mang lại sự bù trừ hoàn hảo về mặt kỹ thuật. Chỉ thị ISSR với 30 mồi giúp quét đa hình trên phạm vi toàn bộ hệ gen mà không cần biết trước trình tự, trong khi 17 cặp mồi SSR cung cấp chỉ thị đồng trội có tính đặc hiệu locus rất cao, giúp xác định chính xác trạng thái đồng hợp tử và dị hợp tử của từng cá thể.

3. Chất lượng mẫu DNA tách chiết từ lá Bách xanh đạt các tiêu chuẩn kỹ thuật phân tử ra sao?

Bằng phương pháp CTAB cải tiến của Porebski, DNA tổng số của 70 mẫu đạt độ tinh sạch rất cao với tỷ lệ quang phổ OD260/OD280 ổn định từ 1,732 đến 2,000. Nồng độ DNA dao động trong khoảng 640–1000 ng/µl, khi điện di trên gel agarose 0,9% cho băng gọn, sắc nét, không bị đứt gãy hay nhiễm tạp chất.

4. Kết quả phân tích AMOVA phản ánh đặc điểm phân hóa di truyền như thế nào?

Phân tích phương sai phân tử AMOVA chứng minh phần lớn sự biến dị di truyền của Bách xanh tập trung bên trong nội bộ từng quần thể (chiếm hơn 55–60% tổng biến thiên), trong khi tỷ lệ sai khác giữa các quần thể ở mức thấp hơn. Kết quả này phản ánh đặc trưng của loài thụ phấn chéo qua gió nhưng đang bắt đầu chịu ảnh hưởng phân hóa do chia cắt địa lý.

5. Cây phát sinh chủng loại UPGMA phân nhóm 70 mẫu Bách xanh dựa trên nguyên lý nào?

Sơ đồ hình cây UPGMA được thiết lập dựa trên ma trận khoảng cách di truyền của Nei và hệ số tương đồng Jaccard thông qua phần mềm NTSYS 2.0, được kiểm định độ tin cậy với 1000 lần lặp bootstrap trên Win-Boot. Cây phả hệ phân tách rõ 70 cá thể thành các nhánh tương ứng với nguồn gốc địa lý tại Lâm Đồng, Đắk Lắk và Gia Lai.

Kết luận

  • Luận văn đã giải quyết trọn vẹn mục tiêu đánh giá toàn diện đa dạng di truyền của 70 cá thể Bách xanh (Calocedrus macrolepis Kurz) thuộc 3 quần thể tự nhiên đại diện tại Lâm Đồng, Đắk Lắk và Gia Lai.
  • Thiết lập thành công quy trình phân tích phân tử chuẩn mực với hệ thống 30 mồi ISSR và 17 cặp mồi SSR, chứng minh hiệu quả phân loại và độ tin cậy thống kê vượt trội.
  • Khẳng định tính đa dạng di truyền phong phú trong nội bộ quần thể nhưng cảnh báo nguy cơ suy giảm dòng gen do sinh cảnh phân mảnh và số lượng cây đại thụ suy giảm.
  • Lộ trình giai đoạn 2026–2030 cần tập trung khoanh vùng bảo vệ nghiêm ngặt các cây mẹ trên 400 năm tuổi, đồng thời phát triển vườn ươm nhân giống phục hồi ngoại vi với quy mô trên 5000 cây giống chuẩn gen mỗi năm.
  • Công trình đóng vai trò là cầu nối quan trọng giữa nghiên cứu sinh học phân tử hiện đại và thực tiễn quản lý tài nguyên rừng, kêu gọi sự chung tay phối hợp giữa các nhà khoa học, cơ quan quản lý và cộng đồng địa phương nhằm bảo tồn vĩnh viễn loài cây lá kim quý hiếm của Việt Nam.