Tổng quan luận án

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

Thông caribê (Pinus caribaea Morelet) là loài cây gỗ lớn bản địa của vùng Trung Mỹ và vùng biển Caribê, bao gồm 3 biến chủng chính: Pinus caribaea var. hondurensis (PCH), Pinus caribaea var. caribaea (PCC) và Pinus caribaea var. bahamensis (PCB). Đây là loài cây sinh trưởng nhanh, thân thẳng, cành nhánh nhỏ hơn so với các loài thông bản địa như Thông ba lá (Pinus kesiya), Thông đuôi ngựa (Pinus massoniana) và Thông nhựa (Pinus merkusii). Gỗ Thông caribê có thớ thẳng, dùng cho xẻ gỗ xây dựng, gỗ trụ mỏ, sản xuất ván dăm, ván bóc, ván ép LVL, bột giấy, đóng tàu thuyền, container và đồ nội thất. Do biên độ sinh thái rộng, sinh trưởng tốt trên đất thoái hóa nghèo dinh dưỡng, loài cây này được Bộ Nông nghiệp và Phát triển nông thôn đưa vào danh mục cây trồng rừng chính theo Quyết định số 4961/QĐ-BNN-TCLN năm 2014 và Thông tư số 22/2021/TT-BNNPTNT ngày 29/12/2021, khuyến khích gây trồng tại 7/9 vùng sinh thái lâm nghiệp của Việt Nam.

Mặc dù có nhiều ưu điểm, diện tích trồng Thông caribê tại Việt Nam đến năm 2020 mới đạt khoảng hơn 2.000 ha trong tổng số 400.000 ha rừng thông cả nước. Nguyên nhân chủ yếu là do khả năng kết hạt của loài kém, chỉ đậu hạt ổn định ở một số vùng như Ba Vì (Hà Nội), Đại Lải (Vĩnh Phúc), Huế và Đà Nẵng. Phần lớn các công trình nghiên cứu trước đây tại Việt Nam mới dừng lại ở khảo nghiệm loài, biến chủng và xuất xứ; việc tiếp cận nhập khẩu hạt giống từ các xuất xứ nguyên sản gặp nhiều trở ngại. Các khảo nghiệm hậu thế trước năm 2017 có số lượng gia đình hạn chế (như khảo nghiệm năm 2003 chỉ có 60–70 gia đình), chưa đại diện đầy đủ cho các nguồn giống hiện có. Đồng thời, các cây mẹ lấy giống trong nước phần lớn thu hái từ các rừng trồng chưa rõ lai lịch xuất xứ ban đầu. Do đó, việc nghiên cứu xác định mức độ biến dị, các tham số di truyền về sinh trưởng, chất lượng thân cây, khối lượng riêng của gỗ và đánh giá đa dạng di truyền bằng chỉ thị phân tử trong quần thể chọn giống là yêu cầu cấp thiết nhằm cung cấp luận cứ khoa học định hướng chiến lược cải thiện giống Thông caribê phục vụ trồng rừng gỗ lớn.

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

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

  1. Nghiên cứu mức độ biến dị về sinh trưởng, một số chỉ tiêu chất lượng thân cây và khối lượng riêng giữa các nguồn giống và các gia đình Thông caribê trong các khảo nghiệm hậu thế.
  2. Nghiên cứu khả năng di truyền, tương quan di truyền và tương tác kiểu gen – hoàn cảnh của các tính trạng sinh trưởng, chất lượng thân cây và khối lượng riêng ở Thông caribê.
  3. Đánh giá mức độ đa dạng di truyền một số nguồn giống Thông caribê ở nước ta bằng chỉ thị phân tử.
  4. Đề xuất một số định hướng nghiên cứu cải thiện giống Thông caribê tại Việt Nam.

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

  • Đối tượng nghiên cứu: Các gia đình hậu thế thụ phấn tự do của các cây trội chọn lọc từ 6 nguồn giống (rừng trồng và rừng giống) gồm: Rừng giống SPA Ba Vì 1 (trồng năm 1996), Rừng giống SPA Ba Vì 2 (trồng năm 1993), Rừng giống SPA Đại Lải (Vĩnh Phúc), nguồn giống Hương Thủy (Thừa Thiên Huế), Tứ Hạ (Thừa Thiên Huế) và Hải Vân (Đà Nẵng).
  • Phạm vi nội dung: Đánh giá biến dị và khả năng di truyền về sinh trưởng ($D_{00}, D_{1.3}, H_{vn}$), chất lượng thân cây (độ thẳng thân, độ nhỏ cành, phát triển ngọn, sức khỏe, khoảng cách vòng cành, tỷ lệ đuôi chồn) đến giai đoạn 5,5 tuổi tại 3 khảo nghiệm; đánh giá khối lượng riêng gỗ tại 2 khảo nghiệm (Ba Vì và Cam Lộ) ở tuổi 4,5 và 5 (40 gia đình/khảo nghiệm, 4 mẫu/gia đình); đánh giá đa dạng di truyền bằng chỉ thị phân tử ISSR đối với 100 gia đình tại khảo nghiệm Ba Vì (5 mẫu lá/gia đình, tổng số 500 mẫu lá).
  • Phạm vi không gian: Nghiên cứu thực hiện tại 3 địa điểm đại diện cho 3 vùng sinh thái lâm nghiệp:
    1. Ba Vì – Hà Nội (vùng Trung tâm Bắc Bộ).
    2. Ngọc Thanh, Phúc Yên – Vĩnh Phúc (vùng Đông Bắc Bộ).
    3. Cam Lộ – Quảng Trị (vùng Bắc Trung Bộ).
  • Phạm vi thời gian: Kế thừa số liệu theo dõi giai đoạn 2016–2020 từ đề tài cấp Bộ do ThS. Cấn Thị Lan làm chủ nhiệm và trực tiếp thu thập, phân tích dữ liệu bổ sung đến năm 2023 (cây đạt 5,5 tuổi).

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

Các hướng nghiên cứu trên thế giới

  1. Phân loại và phân bố tự nhiên: Loài Pinus caribaea được Morelet mô tả lần đầu năm 1851 (dẫn theo Dieters và Nikles, 1997). Ba biến chủng có vùng phân bố tách biệt: P. caribaea var. bahamensis (PCB) phân bố ở vĩ độ 22°–27°N (quần đảo Bahamas, Caicos, bán đảo Yucatan); P. caribaea var. caribaea (PCC) phân bố ở vĩ độ 16°–20°N (Cuba, đảo Juventus); P. caribaea var. hondurensis (PCH) phân bố ở vĩ độ 12°–16°N (Belize, Honduras, Guatemala, Nicaragua) theo các công trình của Greaves và cộng sự (1977), Poynton (1977), Lamb (1973), Perry (1991), Camacho (2018).
  2. Lịch sử di thực và khảo nghiệm xuất xứ: Xuất khẩu hạt giống bắt đầu từ Belize sang Nam Phi và California năm 1927 (Luckhoff, 1964), tiếp tục du nhập vào Queensland – Úc năm 1930. Làn sóng du nhập thứ hai diễn ra từ năm 1962; đến năm 1971 có 35 xuất xứ được khảo nghiệm trên 50 quốc gia nhiệt đới (Lamb, 1973; Baylis và Barnes, 1989; Dvorak và cộng sự, 1993). Viện Lâm nghiệp Oxford xây dựng hệ thống hơn 180 khảo nghiệm trên 53 nước (Greaves, 1980, 1981), xác nhận PCH có tốc độ sinh trưởng nhanh nhất, tiếp đến là PCB và PCC. Eismann và cộng sự (1983) kết luận xếp hạng xuất xứ có sự thay đổi rõ rệt giữa tuổi 8 và 10, do đó cần theo dõi dài hạn trên nhiều lập địa.
  3. Khảo nghiệm hậu thế và tham số di truyền: Dean và cộng sự (1986), Woolaston và cộng sự (1990, 1991), Hodge và Dvorak (2001) tại Úc và Nam Mỹ ghi nhận hệ số di truyền nghĩa hẹp ($h^2$) của các tính trạng sinh trưởng dao động từ 0,10 đến 0,40 và tăng dần theo tuổi cây. Zakaria Ibrahim (1984) tại Malaysia nghiên cứu hiện tượng đuôi chồn (foxtailing), chỉ ra cây đuôi chồn có ưu thế chiều cao ở tuổi non nhưng khối lượng riêng gỗ không sai khác rõ rệt. Moura và cộng sự (2000), Moura và Dvorak (2001) đánh giá mạng lưới CAMCORE tại Brazil, Colombia, Venezuela ở tuổi 12 ghi nhận $h^2$ thể tích thân cây là 0,28–0,34, tính trạng đuôi chồn đạt 0,63. Sebbenn và cộng sự (2006) tại Brazil (tuổi 15), Kalanzi và cộng sự (2014) tại Uganda, Souza và cộng sự (2017) tại Brazil ($h^2 = 0,06–0,51$) và Belaber và cộng sự (2018) tại Argentina trên con lai PCH $\times$ Pinus elliottii ($h^2 = 0,24–0,53$) đều chỉ ra biến dị di truyền lớn giữa các gia đình.
  4. Khối lượng riêng gỗ: Khiari và cộng sự (1989) tại Zanzibar xác định khối lượng riêng Thông caribê 12 tuổi đạt 460 kg/m³, tăng dần từ ruột ra vỏ. Cown (1981) tại Fiji ghi nhận tuổi và độ cao lập địa chi phối 85% biến động khối lượng riêng. Wright (1989) và Oluwadare và cộng sự (2007) khẳng định tính chất gỗ bắt đầu ổn định sau 15 năm tuổi.
  5. Đa dạng di truyền phân tử: Shepherd và cộng sự (2002) phát triển 23/29 cặp mồi microsatellite đa hình. Zheng và Ennos (1999) sử dụng Isozyme chỉ ra dị hợp tử mong đợi $H_e = 0,26$, $F_{st} = 0,08$ ở PCB và 0,02 ở PCC. Tambarussi và cộng sự (2010) dùng 20 cặp mồi SSR chứng minh các gia đình hậu thế là dạng bán noãn phối (half-sibs). Sanchez và cộng sự (2014), Delgado và cộng sự (2011) ($H_o = 0,422$) và Camacho và cộng sự (2018) trên 15 quần thể tự nhiên ($H_e = 0,571$) đã phân tích cấu trúc di truyền và mức độ phân hóa của 3 biến chủng.

Các hướng nghiên cứu tại Việt Nam

  • Gây trồng và khảo nghiệm xuất xứ: Trồng thử nghiệm từ năm 1963 tại Đà Lạt (Lê Đình Khả và Hồ Viết Sắc, 1990). Giai đoạn 1980–1990 nhập hạt phục vụ vùng nguyên liệu giấy (Phí Quang Điện, 1989; Stahl, 1988). Các khảo nghiệm tại Tuyên Quang, Vĩnh Phúc, Quảng Trị, Đồng Nai khẳng định biến chủng PCH sinh trưởng nhanh nhất (Lê Đình Khả, Phí Quang Điện, Đoàn Văn Nhưng, 1988; Lê Đình Khả và cộng sự, 1999, 2001; Phan Thanh Hương, 2000). Ba xuất xứ Poptun, Byfield, Cardwell được Bộ Nông nghiệp và PTNT công nhận là giống tiến bộ kỹ thuật theo Quyết định số 3614/QĐ-BNN-KHCN ngày 08/08/2001. Nguyễn Đình Hải và cộng sự (2002), Nguyễn Ngọc Đích và cộng sự (2004), Dieters và cộng sự (2006) tiếp tục đánh giá sự thích ứng của các xuất xứ tại nhiều vùng lập địa.
  • Khảo nghiệm hậu thế: Năm 2003, dự án Sida-SAREC xây dựng 6 ha khảo nghiệm hậu thế tại Ba Vì, Đồng Hới, Pleiku với 60–70 gia đình (Hà Huy Thịnh, 2006); khảo nghiệm Ba Vì đã được chuyển hóa thành rừng giống năm 2021. Giai đoạn 2016–2020, đề tài cấp Bộ do Cấn Thị Lan chủ nhiệm xây dựng 4 khảo nghiệm hậu thế kết hợp vườn giống thế hệ 1 tại Ba Vì, Đại Lải, Đông Hà, Hương Thủy, ghi nhận sai khác rõ rệt ở tuổi 3 giữa các gia đình và nguồn giống ($Fpr < 0,001$).
  • Chế biến và sử dụng gỗ: Nguyễn Xuân Quyền và cộng sự (2009) nghiên cứu độ bền dán dính keo PVAc sau khi luộc xử lý nhựa; Hà Tiến Mạnh và cộng sự (2021) nghiên cứu tạo ván LVL phối hợp Thông caribê và Bạch đàn Urô; Nguyễn Bảo Ngọc (2020) khảo sát thực trạng sử dụng gỗ thông tại các doanh nghiệp.
  • Đa dạng di truyền phân tử trên các loài thông: Ngô Thị Minh Duyên (2004) dùng chỉ thị SSR trên Pinus merkusii; Đinh Thị Phòng và cộng sự (2015) dùng ISSR trên Pinus dalatensis; Trần Thị Liễu và cộng sự (2015) dùng ISSR trên Pinus krempfii; Nguyễn Minh Tâm và cộng sự (2015) dùng ISSR trên Pinus armandii subsp. phanxiphagensis; Mai Thị Phương Thúy và cộng sự (2020) dùng ISSR và RAPD trên Pinus kwangtungensis.

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

Các công trình trong nước trước đây thiếu dữ liệu đánh giá có hệ thống về tham số di truyền ($h^2, CV_a, r_g$, tương tác $G \times E$, tăng thu di truyền) của Thông caribê ở cấp độ gia đình qua các mốc tuổi liên tục đến 5,5 tuổi trên các vùng sinh thái khác nhau. Đồng thời, chưa có công trình nào đánh giá mối liên hệ giữa sinh trưởng với khối lượng riêng gỗ ở cấp độ gia đình, cũng như chưa có nghiên cứu sử dụng chỉ thị sinh học phân tử (ISSR) để xác định mức độ đa dạng và cấu trúc di truyền của các nguồn giống Thông caribê đang lưu hành tại Việt Nam.

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 vận dụng lý thuyết di truyền số lượng trong lâm nghiệp (Quantitative Forest Genetics) và di truyền quần thể phân tử (Molecular Population Genetics). Khung phân tích gồm:

  • Mô hình tuyến tính phân tích phương sai (ANOVA) để tách các thành phần phương sai di truyền gia đình ($\sigma_f^2$), phương sai lặp lại, phương sai ô thí nghiệm và phương sai sai số.
  • Đánh giá hệ số di truyền nghĩa hẹp ($h^2$) cấp cá thể và cấp gia đình, hệ số biến động di truyền lũy tích ($CV_a$).
  • Đánh giá tương tác kiểu gen – hoàn cảnh ($G \times E$) thông qua phân tích phương sai gộp giữa các địa điểm khảo nghiệm.
  • Phân tích hệ số tương quan kiểu gen ($r_g$) và tương quan kiểu hình ($r_p$) giữa các tính trạng sinh trưởng, chất lượng thân cây và khối lượng riêng gỗ; tương quan giữa các độ tuổi (tuổi 3, tuổi 4,5 và tuổi 5,5).
  • Dự tính tăng thu di truyền lý thuyết ($G$) theo cường độ chọn lọc.
  • Khung phân tích phân tử: Kiểm định định luật Hardy-Weinberg, phân tích phương sai phân tử (AMOVA), tính toán số alen trung bình ($N_a$), số alen hiệu lực ($N_e$), chỉ số Shannon ($I$), độ đa dạng gen Nei ($h$), tỷ lệ đoạn DNA đa hình ($PPB$), dị hợp tử mong đợi ($H_e$), hệ số phân hóa di truyền ($F_{st}$), khoảng cách di truyền Nei và dựng cây phả hệ UPGMA, tọa độ chính PCoA.

Phương pháp nghiên cứu thực tế

Bố trí thí nghiệm và lập địa nghiên cứu

Các khảo nghiệm hậu thế được thiết kế theo khối ngẫu nhiên hàng – cột với 8 lần lặp lại, 3 cây/ô (24 cây/gia đình/địa điểm), cự ly trồng $3 \times 2$ m (mật độ 1.660 cây/ha). Biện pháp kỹ thuật: xử lý thực bì cày toàn diện, đào hố $40 \times 40 \times 40$ cm, bón lót 2 kg phân hữu cơ vi sinh + 300g NPK/hố, bón thúc 200g NPK/cây; chăm sóc 3 năm đầu (phát dọn, xới vun gốc, bảo vệ).

STT Địa điểm khảo nghiệm Số lượng gia đình Thời điểm trồng Đặc điểm khí hậu (Nhiệt độ TB / Lượng mưa TB) Đặc điểm đất đai
1 Ba Vì (Hà Nội) 100 08/2017 23,1°C / 1.700 mm/năm (mưa tập trung tháng 4–10) Feralit màu nâu vàng trên sa thạch, pH 3,5–3,6; nghèo dinh dưỡng
2 Đại Lải (Vĩnh Phúc) 80 08/2017 23,8°C / 1.040 mm/năm (mưa tập trung tháng 8–1) Feralit vàng nhạt trên sa phiến thạch, pH 4,02–4,23; tầng đất mỏng 20–50cm
3 Cam Lộ (Quảng Trị) 80 12/2017 25,1°C / 2.509 mm/năm (mưa nhiều cuối năm) Feralit đỏ vàng trên phiến thạch sét, pH 3,76–4,07; tầng đất dày

Nguồn gốc vật liệu nghiên cứu

Vật liệu gồm hạt giống thu từ các cây mẹ chọn lọc thuộc 6 nguồn giống:

  • SPA Ba Vì 1 (trồng 1996): 20 gia đình (B1–B126).
  • SPA Ba Vì 2 (trồng 1993): 19 gia đình (B130–B260).
  • SPA Đại Lải (Vĩnh Phúc): 23 gia đình (ĐL).
  • Hương Thủy (Huế): 19 gia đình (H).
  • Tứ Hạ (Huế): 6 gia đình (TH).
  • Hải Vân (Đà Nẵng): 13 gia đình (Đ).

Phương pháp thu thập và đo đếm số liệu

  1. Chỉ tiêu sinh trưởng: Đo đường kính gốc ($D_{00}$), đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$) toàn bộ cây theo quy trình điều tra rừng của Vũ Tiến Hinh và Phạm Ngọc Giao (1997) ở các thời điểm 3 tuổi, 4,5 tuổi và 5,5 tuổi.
  2. Chỉ tiêu chất lượng thân cây (mộc trắc 5 cấp):
    • Độ thẳng thân ($D_{tt}$): Cấp 1 (rất cong) đến Cấp 5 (rất thẳng) theo Lê Đình Khả và cộng sự (2003).
    • Độ nhỏ cành ($D_{nc}$): Cấp 1 ($d_{cành} \ge 1/2 d_{thân}$) đến Cấp 5 ($d_{cành} < 1/10 d_{thân}$) theo TCVN 8761:2017.
    • Phát triển ngọn ($P_{tn}$): Cấp 1 (mất ngọn, cụt ngọn) đến Cấp 5 (1 ngọn chính phát triển mạnh, tán tháp nhọn cân đối).
    • Sức khỏe ($S_k$): Cấp 1 (rất kém, lá úa khô) đến Cấp 5 (rất tốt, tán xanh cân đối, không sâu bệnh) theo Lê Đình Khả và Dương Mộng Hùng (2003).
    • Chỉ tiêu tổng hợp ($I_{cl}$): Tính theo công thức $I_{cl} = (D_{tt} + D_{nc} + P_{tn} + S_k)/4$.
    • Khoảng cách vòng cành ($K_{cvc}$): Đếm số vòng cành/cây theo DeYoung (2016), tính bằng $H_{vn} / \text{số vòng cành}$.
    • Tỷ lệ đuôi chồn: Xác định theo Zakaria Ibrahim (1984), ghi nhận giá trị 1 (đuôi chồn) hoặc 0 (bình thường).
  3. Xác định khối lượng riêng của gỗ (KLR): Lấy mẫu ngẫu nhiên tại 2 khảo nghiệm Ba Vì và Cam Lộ ở tuổi 4,5 và tuổi 5 (40 gia đình/khảo nghiệm, 4 mẫu/gia đình) bằng phương pháp thể tích – khối lượng khô tuyệt đối.
  4. Phân tích đa dạng di truyền bằng chỉ thị phân tử ISSR:
    • Thu 500 mẫu lá non (100 gia đình $\times$ 5 mẫu lá tại 5 lần lặp ở Ba Vì).
    • Tách chiết ADN tổng số, kiểm tra độ tinh sạch và nồng độ.
    • Sàng lọc các cặp mồi ISSR (mồi UBC855, UBC856,...), thực hiện phản ứng khuếch đại PCR, điện di trên gel agarose, đọc bản gel và lập ma trận nhị phân (1: có băng, 0: không có băng).

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

Chương 1: Tổng quan các vấn đề nghiên cứu

Chương 1 (23 trang) hệ thống hóa các tài liệu khoa học trong và ngoài nước về Thông caribê. Tác giả tổng thuật chi tiết về đặc điểm thực vật học, sinh thái phân bố tự nhiên của 3 biến chủng (PCH, PCC, PCB), quá trình di thực toàn cầu qua hai làn sóng chính và kết quả khảo nghiệm xuất xứ, hậu thế từ mạng lưới quốc tế (Viện Lâm nghiệp Oxford, CAMCORE) tại Úc, Nam Phi, Brazil, Colombia, Venezuela, Malaysia, Uganda, Argentina. Chương này cũng phân tích các nghiên cứu về tính chất gỗ, hiện tượng đuôi chồn và ứng dụng chỉ thị phân tử (Isozyme, SSR). Đối với Việt Nam, chương 1 điểm lại lịch sử trồng khảo nghiệm từ năm 1963 tại Đà Lạt, các đợt khảo nghiệm xuất xứ (công nhận giống Poptun, Byfield, Cardwell năm 2001), khảo nghiệm hậu thế năm 2003 (Sida-SAREC), đề tài cấp Bộ giai đoạn 2016–2020 và các nghiên cứu sử dụng gỗ thông làm ván ghép thanh, ván LVL. Qua đó, chương 1 làm nổi bật khoảng trống nghiên cứu về tham số di truyền cấp gia đình và di truyền phân tử tại Việt Nam.

Chương 2: Nội dung, vật liệu và phương pháp nghiên cứu

Chương 2 (16 trang) trình bày chi tiết về 4 nội dung nghiên cứu, mô tả 6 nguồn vật liệu với 100 gia đình thử nghiệm tại 3 vùng sinh thái (Ba Vì, Đại Lải, Cam Lộ). Chương này cung cấp số liệu khí tượng (nhiệt độ, lượng mưa, độ bốc hơi) và lý hóa tính đất (pH, mùn, đạm, lân, kali, cấp hạt cát/sét) tại từng địa điểm. Toàn bộ quy trình điều tra sinh trưởng, phương pháp cho điểm mộc trắc chất lượng thân cây theo tiêu chuẩn TCVN 8761:2017, phương pháp xác định khối lượng riêng gỗ và quy trình sinh học phân tử sử dụng chỉ thị ISSR (tách chiết ADN, phản ứng PCR, xử lý ma trận số liệu phân tử) được quy định chặt chẽ, tạo cơ sở cho việc phân tích ở chương sau.

Chương 3: Kết quả nghiên cứu và thảo luận

Chương 3 (78 trang) là nội dung trung tâm của luận án, trình bày và thảo luận các kết quả đạt được:

Biến dị về sinh trưởng, chất lượng thân cây và khối lượng riêng

  • Biến dị giữa các nguồn giống: Tại Ba Vì (5,5 tuổi), Đại Lải (5,5 tuổi) và Cam Lộ (4 và 5 tuổi), có sự sai khác rõ rệt giữa 6 nguồn giống về các chỉ tiêu $D_{1.3}, H_{vn}$, thể tích thân cây ($VOL$), độ thẳng thân ($D_{tt}$), độ nhỏ cành ($D_{nc}$), chỉ tiêu tổng hợp ($I_{cl}$) và tỷ lệ đuôi chồn. Các nguồn giống SPA Ba Vì 1, SPA Ba Vì 2 và Hải Vân đạt sinh trưởng tốt hơn so với các nguồn giống còn lại.
  • Biến dị giữa các gia đình: Sự sai khác về sinh trưởng giữa các gia đình trong cùng nguồn giống đạt mức tin cậy thống kê cao ($Fpr < 0,001$). Tồn tại sự phân hóa rõ nét giữa nhóm gia đình sinh trưởng tốt, trung bình và xấu.
  • Hiện tượng đuôi chồn: Tỷ lệ cây đuôi chồn có sự biến động tùy thuộc vào nguồn giống, gia đình và địa điểm khảo nghiệm. Cây có dạng hình đuôi chồn thể hiện ưu thế về tăng trưởng chiều cao ở tuổi non nhưng hình thái ngọn kém phân cành ngang.
  • Khối lượng riêng gỗ: Có sự sai khác có ý nghĩa thống kê về khối lượng riêng gỗ giữa các gia đình được kiểm tra tại Ba Vì và Cam Lộ ở giai đoạn 4,5–5 tuổi.

Khả năng di truyền, tương quan di truyền và tăng thu di truyền

  • Hệ số di truyền ($h^2$) và hệ số biến động di truyền lũy tích ($CV_a$): Hệ số di truyền nghĩa hẹp của các tính trạng sinh trưởng ($D_{1.3}, H_{vn}, VOL$) ở mức trung bình đến khá ($h^2 = 0,15–0,45$) và có xu hướng tăng dần theo tuổi cây (từ tuổi 3 đến 5,5 tuổi). Hệ số di truyền của các tính trạng chất lượng thân cây ($D_{tt}, D_{nc}, I_{cl}$) đạt mức trung bình, trong khi tính trạng đuôi chồn có hệ số di truyền tương đối cao.
  • Tương tác kiểu gen – hoàn cảnh ($G \times E$): Phân tích gộp giữa 3 khảo nghiệm chỉ ra có sự tương tác kiểu gen – lập địa đối với sinh trưởng thể tích và độ thẳng thân, biểu hiện qua sự thay đổi thứ bậc xếp hạng của một số gia đình giữa vùng Bắc Trung Bộ (Cam Lộ) và vùng Đông Bắc Bộ/Trung tâm Bắc Bộ (Đại Lải, Ba Vì).
  • Tương quan di truyền ($r_g$) và kiểu hình ($r_p$): Tương quan kiểu gen giữa đường kính và chiều cao rất chặt chẽ ($r_g > 0,80$). Tương quan giữa sinh trưởng ở tuổi 3 với tuổi 5,5 đạt mức cao, hỗ trợ việc đánh giá sớm. Mối tương quan di truyền giữa sinh trưởng và khối lượng riêng gỗ ở giai đoạn 4,5–5 tuổi ở mức thấp đến trung bình.
  • Tăng thu di truyền lý thuyết ($G$): Khi áp dụng tỷ lệ chọn lọc 10–20% các gia đình tốt nhất, mức tăng thu di truyền lý thuyết dự tính về thể tích thân cây đạt từ 12% đến trên 20% tùy theo từng địa điểm khảo nghiệm.

Đa dạng di truyền và quan hệ di truyền bằng chỉ thị ISSR

  • Sàng lọc mồi: Các mồi ISSR đa hình (như UBC855, UBC856) nhân bản được nhiều phân đoạn DNA rõ nét, tỷ lệ phân đoạn đa hình ($PPB$) đạt mức cao.
  • Chỉ số đa dạng di truyền: Quần thể chọn giống duy trì mức đa dạng di truyền khá cao biểu hiện qua các thông số: số alen trung bình ($N_a$), số alen hiệu lực ($N_e$), độ đa dạng gen Nei ($h$) và chỉ số Shannon ($I$).
  • Phân tích AMOVA: Biến đổi di truyền phân tử chủ yếu xuất hiện bên trong các nguồn giống (chiếm tỷ lệ lớn), sự sai khác di truyền giữa 6 nguồn giống chiếm tỷ lệ nhỏ hơn nhưng có ý nghĩa thống kê.
  • Quan hệ di truyền: Khoảng cách di truyền Nei và sơ đồ UPGMA/PCoA phân nhóm rõ các gia đình và xác lập mối quan hệ phát sinh giữa 6 nguồn giống nghiên cứu (các nguồn giống Ba Vì, Đại Lải, miền Trung).
  • Định hướng cải thiện giống: Đề xuất lộ trình tỉa thưa di truyền để chuyển hóa khảo nghiệm Ba Vì và Đại Lải thành vườn giống hữu tính thế hệ 1, đồng thời lựa chọn các cây trội ưu việt phục vụ nhân giống sinh dưỡng và lai giống cho chu kỳ tiếp theo.

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

Đóng góp mới về mặt khoa học

  1. Xác định có hệ thống các tham số di truyền số lượng gồm hệ số di truyền nghĩa hẹp ($h^2$), hệ số biến động di truyền lũy tích ($CV_a$), tương quan kiểu gen ($r_g$) và tương quan kiểu hình ($r_p$) về sinh trưởng, chất lượng thân cây và khối lượng riêng gỗ của các gia đình Thông caribê đến 5,5 tuổi trên 3 vùng sinh thái đại diện tại Việt Nam.
  2. Cung cấp dẫn liệu khoa học về mức độ tương tác kiểu gen – hoàn cảnh ($G \times E$) và tương quan sớm – muộn (tuổi 3 so với tuổi 5,5), làm cơ sở cho công tác chọn lọc sớm trong cải thiện giống Thông caribê.
  3. Làm sáng tỏ mức độ đa dạng di truyền, biến đổi phân tử (AMOVA) và khoảng cách di truyền của 6 nguồn giống Thông caribê tại Việt Nam thông qua hệ thống chỉ thị phân tử ISSR.

Đóng góp mới về mặt thực tiễn

  1. Đánh giá và phân hạng được các nguồn giống tốt (SPA Ba Vì 1, SPA Ba Vì 2, Hải Vân) và tuyển chọn được tập hợp các gia đình cây trội ưu việt có năng suất cao, thân thẳng, cành nhỏ, sức khỏe tốt tại Ba Vì (Hà Nội), Đại Lải (Vĩnh Phúc) và Cam Lộ (Quảng Trị).
  2. Xác định được mức tăng thu di truyền lý thuyết khi chọn lọc, cung cấp vật liệu khởi đầu và cơ sở khoa học để tỉa thưa di truyền chuyển hóa các khảo nghiệm hậu thế thành vườn giống hữu tính cung cấp hạt giống chất lượng cao phục vụ trồng rừng gỗ lớn.

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

  • Tiếp tục quản lý, bảo vệ và theo dõi các khảo nghiệm hậu thế ở các giai đoạn tuổi lớn hơn (8–10 tuổi và cuối luân kỳ) để chuẩn xác hóa các tham số di truyền về sinh trưởng và phẩm chất gỗ.
  • Tiến hành tỉa thưa di truyền tại các khảo nghiệm hậu thế ở Ba Vì và Đại Lải, loại bỏ các gia đình kém sinh trưởng và có dạng hình xấu để chuyển đổi thành vườn giống hữu tính thế hệ 1.
  • Nhân rộng các gia đình ưu trội đã được tuyển chọn vào các chương trình trồng rừng sản xuất gỗ lớn tại các vùng sinh thái tương ứng (Bắc Trung Bộ, Trung tâm Bắc Bộ, Đông Bắc Bộ).

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

Hạn chế của luận án

  • Thời gian theo dõi sinh trưởng và chất lượng thân cây của các khảo nghiệm hậu thế mới dừng lại ở giai đoạn 5,5 tuổi, chưa đánh giá hết toàn bộ chu kỳ kinh doanh cây gỗ lớn.
  • Dữ liệu về khối lượng riêng gỗ mới được thu thập trên 40 gia đình tại 2 địa điểm khảo nghiệm (Ba Vì và Cam Lộ) ở mốc 4,5 và 5 tuổi, chưa khảo sát toàn diện các chỉ tiêu cơ lý khác như độ bền uốn tĩnh, mô đun đàn hồi hay độ co rút.
  • Phân tích đa dạng di truyền bằng chỉ thị phân tử ISSR mới thực hiện trên 100 gia đình thuộc khảo nghiệm Ba Vì; lai lịch ban đầu của một số rừng trồng chọn cây mẹ chưa được ghi nhận đầy đủ qua hồ sơ lưu trữ lịch sử.

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

  • Tiếp tục đo đếm, đánh giá các khảo nghiệm hậu thế ở tuổi 8–10 nhằm xác định độ ổn định xếp hạng gia đình và đánh giá toàn diện các tính chất gỗ thành thục.
  • Mở rộng phân tích di truyền phân tử bằng các chỉ thị đồng trội có độ lặp lại cao (như SSR hoặc SNP) trên toàn bộ mạng lưới khảo nghiệm để đánh giá dòng chảy gen và mức độ cận huyết trong vườn giống sau tỉa thưa.
  • Xây dựng khảo nghiệm hậu thế thế hệ 2 từ các cây trội được chọn lọc trong thế hệ 1 kết hợp khảo nghiệm nhân giống sinh dưỡng (hom vô tính).

Giá trị tham khảo

Đối tượng sử dụng

  • Giảng viên, nghiên cứu sinh, học viên cao học: Chuyên ngành Di truyền và chọn giống cây lâm nghiệp, Lâm học, Công nghệ sinh học lâm nghiệp sử dụng làm tài liệu tham khảo về phương pháp nghiên cứu di truyền số lượng và di truyền phân tử cây rừng.
  • Cán bộ nghiên cứu tại các viện/trung tâm lâm nghiệp: Vận dụng khung thiết kế thí nghiệm hàng – cột, phương pháp đánh giá mộc trắc chất lượng thân cây theo TCVN 8761:2017 và quy trình phân tích tham số di truyền trên cây gỗ lá kim.
  • Nhà quản lý lâm nghiệp và doanh nghiệp trồng rừng: Cơ quan quản lý cấp Bộ, Sở Nông nghiệp và PTNT các tỉnh và các công ty lâm nghiệp khai thác dữ liệu về khả năng thích ứng của các nguồn giống để quy hoạch vùng trồng rừng Thông caribê gỗ lớn phù hợp tại Đông Bắc Bộ, Trung tâm Bắc Bộ và Bắc Trung Bộ.

Các phần đáng tham khảo trong luận án

  • Chương 2 (Mục 2.5): Quy chuẩn thang điểm mộc trắc đánh giá chất lượng thân cây (độ thẳng thân, độ nhỏ cành, phát triển ngọn, sức khỏe, khoảng cách vòng cành, tỷ lệ đuôi chồn) và phương pháp tính chỉ tiêu tổng hợp $I_{cl}$.
  • Chương 3 (Mục 3.2): Bảng tổng hợp hệ số di truyền ($h^2$), hệ số biến động di truyền lũy tích ($CV_a$), ma trận tương quan di truyền và ước tính tăng thu di truyền ($G$) cho các tính trạng kinh tế.
  • Chương 3 (Mục 3.4): Quy trình sàng lọc mồi ISSR và bảng chỉ số đa dạng di truyền, khoảng cách di truyền giữa 6 nguồn giống Thông caribê tại Việt Nam.

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

1. Luận án nghiên cứu những nguồn giống Thông caribê nào và bố trí khảo nghiệm tại những địa phương nào?

Luận án nghiên cứu 6 nguồn giống cây mẹ gồm: Rừng giống SPA Ba Vì 1 (trồng 1996), Rừng giống SPA Ba Vì 2 (trồng 1993), Rừng giống SPA Đại Lải (Vĩnh Phúc), nguồn giống Hương Thủy (Huế), Tứ Hạ (Huế) và Hải Vân (Đà Nẵng). Các nguồn giống này được bố trí trong 3 khảo nghiệm hậu thế kết hợp vườn giống tại:

  • Cẩm Lĩnh, Ba Vì, Hà Nội (100 gia đình).
  • Ngọc Thanh, Phúc Yên, Vĩnh Phúc (80 gia đình).
  • Cam Lộ, Quảng Trị (80 gia đình).

2. Thiết kế thí nghiệm tại các hiện trường khảo nghiệm hậu thế được thực hiện như thế nào?

Thí nghiệm được thiết kế theo sơ đồ khối ngẫu nhiên hàng – cột với 8 lần lặp lại, 3 cây/ô (mỗi gia đình có 24 cây/khảo nghiệm), cự ly trồng $3 \times 2$ m (mật độ 1.660 cây/ha). Biện pháp kỹ thuật bao gồm cày toàn diện, hố đào kích thước $40 \times 40 \times 40$ cm, bón lót 2 kg phân hữu cơ + 300g NPK/hố, bón thúc 200g NPK/cây và chăm sóc trong 3 năm đầu.

3. Các chỉ tiêu chất lượng thân cây được đánh giá bằng những phương pháp và tiêu chuẩn nào?

Chất lượng thân cây được đánh giá qua phương pháp mộc trắc cho điểm theo thang 5 cấp (từ 1 đến 5 điểm) cho 4 tính trạng:

  • Độ thẳng thân ($D_{tt}$) theo phương pháp của Lê Đình Khả và cộng sự (2003).
  • Độ nhỏ cành ($D_{nc}$) theo tiêu chuẩn TCVN 8761:2017.
  • Phát triển ngọn ($P_{tn}$) và Sức khỏe cây ($S_k$) theo Lê Đình Khả và Dương Mộng Hùng (2003).
  • Tính chỉ tiêu tổng hợp $I_{cl} = (D_{tt} + D_{nc} + P_{tn} + S_k)/4$.
  • Khoảng cách vòng cành ($K_{cvc}$) tính bằng $H_{vn} / \text{số vòng cành}$ theo DeYoung (2016).
  • Tỷ lệ cây đuôi chồn ghi nhận theo Zakaria Ibrahim (1984).

4. Loại chỉ thị phân tử nào được tác giả sử dụng để đánh giá đa dạng di truyền và quy mô mẫu ra sao?

Tác giả sử dụng chỉ thị phân tử ISSR (Inter Simple Sequence Repeat) với các mồi đa hình (như UBC855, UBC856,...). Nghiên cứu thực hiện trên 500 mẫu lá non thu từ 100 gia đình (mỗi gia đình thu 5 mẫu lá tại 5 lần lặp lại) thuộc khảo nghiệm hậu thế Thông caribê tại Ba Vì, Hà Nội.

5. Kết quả đánh giá hệ số di truyền ($h^2$) của các tính trạng sinh trưởng ở Thông caribê thể hiện xu hướng gì?

Hệ số di truyền nghĩa hẹp ($h^2$) của các tính trạng sinh trưởng ($D_{1.3}, H_{vn}, VOL$) ở các khảo nghiệm đạt mức từ trung bình đến khá ($h^2 = 0,15–0,45$) và có xu hướng tăng dần theo độ tuổi của cây (từ tuổi 3 đến tuổi 5,5). Tính trạng thể tích thân cây và đường kính thể hiện biến dị di truyền gia đình rõ nét, cho phép đạt mức tăng thu di truyền lý thuyết từ 12% đến trên 20% khi tiến hành chọn lọc gia đình ưu trội.

Kết luận

Luận án Tiến sĩ Lâm nghiệp của nghiên cứu sinh Trần Đắc Vượng (2023) đã giải quyết trọn vẹn các mục tiêu nghiên cứu về đặc điểm biến dị, tham số di truyền số lượng và cấu trúc di truyền phân tử của loài Thông caribê (Pinus caribaea Morelet) tại 3 vùng sinh thái đại diện của Việt Nam. Bằng việc kết hợp phương pháp khảo nghiệm hậu thế truyền thống ở giai đoạn đến 5,5 tuổi với kỹ thuật chỉ thị phân tử ISSR, công trình đã xác định được các tham số di truyền ($h^2, CV_a, r_g, G \times E$), làm rõ mức độ đa dạng di truyền của 6 nguồn giống và chọn lọc được các gia đình cây trội ưu việt. Kết quả nghiên cứu cung cấp hệ thống luận cứ khoa học và nguồn vật liệu di truyền phục vụ trực tiếp cho chiến lược cải thiện giống, tỉa thưa chuyển hóa rừng giống và phát triển rừng trồng Thông caribê cung cấp gỗ lớn tại Việt Nam.