Bối cảnh và vấn đề nghiên cứu

Chi nấm Cordyceps sensu lato (s.l) thuộc ngành nấm túi Ascomycota là nhóm nấm ký sinh trên côn trùng (ấu trùng và cá thể trưởng thành), nổi tiếng trong y học cổ truyền và dược liệu học hiện đại nhờ chứa các hoạt chất có giá trị sinh học cao như cordycepin, polysaccharide tan trong kiềm, adenosine cùng nhiều hợp chất điều hòa miễn dịch và kháng khối u. Nổi bật trong chi này là các loài như Ophiocordyceps sinensis hay Cordyceps militaris. Tuy nhiên, việc định danh và phân loại các loài nấm thuộc chi Cordyceps s.l gặp nhiều trở ngại lớn khi chỉ dựa trên hình thái học truyền thống:

  1. Hình thái và cấu trúc quả thể nấm dễ biến đổi theo môi trường sống, điều kiện nhiệt độ, độ ẩm và loại vật chủ ký sinh.
  2. Vòng đời của nhiều loài nấm tồn tại hai trạng thái sinh sản khác biệt gồm giai đoạn vô tính (anamorph) và giai đoạn hữu tính (teleomorph), dẫn đến việc mô tả trùng lặp hoặc phân loại nhầm lẫn giữa các chi liên quan.
  3. Mẫu nấm thu hái ngoài tự nhiên thường không còn nguyên vẹn, bị gãy vỡ, biến màu hoặc tổn hại sau thời gian dài bảo quản, làm giảm độ chính xác của khóa phân loại hình thái.

Để giải quyết vấn đề phân loại hình thái không đầy đủ, phương pháp định danh phân tử và mã vạch DNA (DNA barcode) kết hợp phân tích đa gen giữ nhà (housekeeping genes) trở thành công cụ phân loại chính xác, hỗ trợ sắp xếp lại quan hệ phát sinh loài của chi Cordyceps s.l.

Mục tiêu và nhiệm vụ nghiên cứu của đề tài:

  1. Đánh giá và kiểm tra in silico các cặp mồi khuếch đại ba vùng gen chỉ thị gồm vùng đệm phiên mã nội (ITS), tiểu đơn vị lớn ribosome nhân (nrLSU) và gen ty thể mã hóa ATP synthase subunit 6 (Atp6).
  2. Tách chiết, thu nhận DNA tổng số và tối ưu hóa phản ứng PCR khuếch đại thành công các đoạn gen ITS, nrLSU, Atp6 từ các mẫu nấm ký sinh côn trùng thu thập tại thực địa.
  3. Đọc, xử lý và hiệu chỉnh tín hiệu giải trình tự nucleotide của các sản phẩm PCR thu được.
  4. Xây dựng bộ cơ sở dữ liệu trình tự DNA barcode đơn gen và đa gen kết hợp (nrLSU-ITS-Atp6, nrLSU-ITS).
  5. Xây dựng và đánh giá cây phát sinh chủng loại phân tử bằng các thuật toán Neighbor-Joining (NJ), Maximum Likelihood (ML) và Maximum Parsimony (MP) để định danh chính xác các mẫu nấm khảo sát.

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

  • Đối tượng nghiên cứu: Các mẫu nấm ký sinh côn trùng ký hiệu DL0036, DL0091, DL0099 (thuộc chi Cordyceps s.l) do Tiến sĩ Trương Bình Nguyên (Khoa Sinh học, Trường Đại học Đà Lạt) thu thập và cung cấp.
  • Phạm vi không gian: Mẫu sinh học có nguồn gốc thu nhận tại khu vực vùng núi Langbiang, tỉnh Lâm Đồng.
  • Phạm vi thời gian và địa điểm thực nghiệm: Đề tài được thực hiện từ tháng 12/2020 đến tháng 06/2021 tại Phòng thí nghiệm Sinh học Phân tử, Trường Đại học Mở Thành phố Hồ Chí Minh (Cơ sở Bình Dương: số 68 Lê Thị Trung, TP. Thủ Dầu Một, tỉnh Bình Dương).

Cơ sở lý thuyết và phương pháp

Khung lý thuyết và hệ thống phân loại

Đề tài dựa trên hệ thống phân loại nấm hiện đại theo Kepler et al. (2017) và nghiên cứu phân loại phát sinh loài của Sung et al. (2007). Nhóm nấm Cordyceps s.l ban đầu thuộc họ Clavicipitaceae đã được sắp xếp lại thành ba họ chính dựa trên dữ liệu đa gen:

  • Họ Clavicipitaceae (chứa các chi Metacordyceps, Hypocrella, Regiocrella, Torrubiella).
  • Họ Cordycipitaceae (chứa chi Cordyceps).
  • Họ Ophiocordycipitaceae (chứa các chi Ophiocordyceps, Elaphocordyceps).

Khung mã vạch DNA được xây dựng dựa trên ba vùng locus bảo tồn và phân kỳ:

  • Vùng ITS (Internal Transcribed Spacer): Gồm ITS1-5.8S-ITS2, nằm trong cụm gen rDNA nhân. Vùng này có tốc độ tiến hóa nhanh, mức độ biến dị cao giữa các loài gần gũi, là mã vạch chuẩn quốc tế cho giới Nấm.
  • Vùng nrLSU (nuclear ribosomal large subunit - 28S rRNA): Mã hóa tiểu phần lớn của ribosome nhân, mang tính bảo tồn cao, phù hợp cho phân loại ở cấp độ chi, họ và bộ.
  • Gen Atp6 (mitochondrial ATP synthase subunit 6): Gen nằm trên hệ gen ty thể mã hóa cho tiểu đơn vị số 6 của phức hợp F1-F0 ATPase. Do nằm trên DNA ty thể và có kích thước phân tử nhỏ, gen Atp6 phản ánh rõ nét các đột biến nucleotide tích lũy trong quá trình tiến hóa, hỗ trợ giải quyết quan hệ phân loại sâu giữa các nhánh trong họ Clavicipitaceae.

Phương pháp nghiên cứu thực nghiệm và công cụ xử lý

Quy trình nghiên cứu kết hợp giữa phân tích tin sinh học (in silico) và thực nghiệm phân tử (in vitro):

Mẫu nấm thực địa (DL0036, DL0091, DL0099)
Tách chiết DNA (Phenol/Chloroform/Isoamyl Alcohol + CTAB + β-mercaptoethanol)
Kiểm tra độ tinh sạch và nồng độ (Quang phổ kế OD260/OD280)
Khuếch đại PCR (Cặp mồi đặc hiệu: ITS1/ITS4, LR0R/LR5, ATP6-C1A/ATP6-C2A)
Điện di gel agarose 1.5% & Gửi giải trình tự Sanger (Nam Khoa Biotech)
Hiệu chỉnh trình tự (Chromas Lite, NCBI Dot plot, BLASTn)
Xây dựng bộ dữ liệu đa gen & Dò tìm mô hình tiến hóa tối ưu (MEGA X)
Hạng mục phương pháp Quy trình / Hóa chất / Thiết bị cụ thể Mục đích kỹ thuật
Phân tích In Silico IDT OligoAnalyzer, BLASTn (NCBI), Annhyb Đánh giá thông số lý hóa mồi (Tm, %GC, cấu trúc kẹp tóc, dimer) và mô phỏng kích thước PCR
Tách chiết DNA Dung dịch đệm Tris-HCl, EDTA, SDS, NaCl, CTAB, β-mercaptoethanol, Proteinase K, PCI (25:24:1), NaOAc, Isopropanol, Ethanol 70% Phá vỡ vách tế bào chitin/glucan bền vững của nấm túi, loại bỏ protein và thu nhận DNA sạch
Kiểm tra DNA Đo quang phổ hấp thụ tại bước sóng 260 nm và 280 nm (tỷ số $OD_{260}/OD_{280}$) Đánh giá độ tinh sạch và nồng độ DNA trước khi thực hiện PCR
Phản ứng PCR Master Mix (7.5 µl), mồi xuôi (0.5 µl), mồi ngược (0.5 µl), DNA khuôn (2.0 µl), $ddH_2O$ (4.5 µl); luân nhiệt 35 chu kỳ Khuếch đại đặc hiệu các đoạn gen ITS, nrLSU, Atp6
Điện di kiểm tra Gel agarose 1.5% trong đệm TBE 1X, nhuộm GelRed, chạy ở hiệu điện thế 100V trong 30 phút Kiểm tra sự hiện diện, kích thước và độ sắc nét của băng sản phẩm PCR
Giải trình tự Giải trình tự theo phương pháp Sanger hai chiều (Forward và Reverse) Xác định chính xác trật tự nucleotide của các đoạn gen
Xử lý tín hiệu Phần mềm Chromas Lite, NCBI Dot plot Đọc tín hiệu peak huỳnh quang, đảo mạch bù (reverse complement), loại bỏ lỗi chèn/mất đoạn (indels)
Phân tích tiến hóa Phần mềm MEGA X, kiểm định mô hình thay thế nucleotide theo tiêu chí BIC/AIC Dò tìm mô hình tiến hóa tốt nhất và xây dựng cây phát sinh loài (NJ, ML, MP, Bootstrap 1000)

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

Chương 1: Tổng quan tài liệu

Chương 1 hệ thống hóa toàn bộ lý thuyết nền tảng về chi nấm Cordyceps s.l, lịch sử phân loại học và các phương pháp định danh sinh học phân tử:

  • Đặc điểm sinh học và dược lý: Trình bày chi tiết chu kỳ ký sinh 3 giai đoạn của nấm trên côn trùng (xâm nhiễm mùa hè, sinh trưởng mùa thu giết chết vật chủ biến thành "Đông trùng", hình thành quả thể chứa stroma vào mùa xuân-hè năm sau thành "Hạ thảo"). Tổng hợp các nghiên cứu về hoạt tính sinh học từ Ophiocordyceps sinensis và Cordyceps militaris trong hỗ trợ điều trị ung thư, điều hòa hệ miễn dịch, bảo vệ chức năng gan, thận, tim mạch.
  • Hệ thống phân loại: Khái quát các tiêu chuẩn phân loại từ khóa hình thái của Kobayasi (1982) đến khóa phân loại đa gen phân tử của Sung et al. (2007) và Kepler et al. (2017), phân tách nhóm Cordyceps truyền thống thành ba họ độc lập.
  • Đặc điểm cấu trúc thành tế bào nấm: Nêu rõ thành phần vách nấm Ascomycota gồm hơn 80% polysaccharide ($\alpha, \beta$-glucan đan xen chitin liên kết chéo với protein), giải thích nguyên nhân gây khó khăn khi ly giải tế bào để tách chiết DNA và cơ chế bổ trợ của CTAB cùng $\beta$-mercaptoethanol nhằm phá vỡ liên kết disulfide của protein.
  • Đặc điểm các vùng gen chỉ thị: Trình bày chi tiết cấu trúc vùng rDNA (ITS1-5.8S-ITS2, nrLSU 28S rRNA) và gen ty thể Atp6, phân tích ưu thế của từng vùng gen trong nghiên cứu phả hệ.
  • Các bước xây dựng cây phát sinh loài: Hệ thống hóa 4 bước chuẩn gồm: (1) Sắp giống cột trình tự (Multiple sequence alignment) loại bỏ indels; (2) Dò tìm mô hình tiến hóa tối ưu (dựa trên tiêu chí BIC và AIC); (3) Dựng cây phả hệ (so sánh thuật toán khoảng cách NJ với thuật toán đặc tính MP, ML); (4) Đánh giá độ tin cậy của các nút phân nhánh qua chỉ số tái lấy mẫu Bootstrap (giá trị $\ge 50%$ được ủng hộ, $\ge 70%$ có ý nghĩa thống kê cao).
  • Tình hình nghiên cứu trong nước: Điểm qua các công trình tại Việt Nam của Đái Duy Ban et al. (2009), Lê Thị Anh Đào et al. (2010), Phạm Thị Hạnh et al. (2011), Đinh Minh Hiệp et al. (2014) và Vũ Tiến Luyện et al. (2014) về việc phân lập, nghiên cứu hoạt chất và định danh các loài nấm Cordyceps tại Langbiang - Lâm Đồng.

Chương 2: Vật liệu và phương pháp nghiên cứu

Chương 2 mô tả chi tiết nguồn vật liệu, hóa chất, thiết bị và các quy trình thao tác kỹ thuật:

  • Mẫu sinh học: Sử dụng 3 mẫu nấm ký sinh côn trùng ký hiệu DL0036, DL0091, DL0099 thu thập tại vùng núi Langbiang (Lâm Đồng).
  • Thông số mồi PCR: Sử dụng ba cặp mồi phổ quát đã được công bố quốc tế để khuếch đại ba locus gen đích:
Vùng gen Tên mồi Trình tự nucleotide ($5' \rightarrow 3'$) Chiều dài Tài liệu tham khảo
ITS ITS1 (Xuôi)
ITS4 (Ngược)		 TCC GTA GGT GAA CCT GCG G
TCC TCC GCT TAT TGA TAT GC		 19 bp
20 bp		 White et al., 1990
nrLSU LR0R (Xuôi)
LR5 (Ngược)		 GTA CCC GCT GAA CTT AAG C
ATC CTG AGG GAA ACT TC		 19 bp
17 bp		 Vilgalys & Sun, 1994
Atp6 ATP6-C1A (Xuôi)
ATP6-C2A (Ngược)		 AGA WCA ATT YGA ART RAG AG
ACA AAY ACT TGW GCT TGK ATW AAI GC		 20 bp
26 bp		 Castlebury et al., 2004

(Trong đó ký hiệu degenerate: W = A/T; Y = C/T; R = A/G; K = G/T; I = inosine/bất kỳ base nào).

  • Quy trình tách chiết DNA: Nghiền mẫu mô nấm trong $100\ \mu\text{l}\ ddH_2O$, bổ sung $70\ \mu\text{l}$ Tris-HCl, $14\ \mu\text{l}$ EDTA, $140\ \mu\text{l}$ SDS, $140\ \mu\text{l}$ NaCl 2.5M, $70\ \mu\text{l}$ CTAB, $7\ \mu\text{l}\ \beta$-mercaptoethanol và $30\ \mu\text{l}$ Proteinase K; ủ qua đêm ở 65°C. Tiến hành chiết bằng Phenol/Chloroform/Isoamyl alcohol (25:24:1), kết tủa DNA bằng NaOAc 3M và isopropanol ở 4°C trong 2 giờ, rửa tủa bằng ethanol 70%, làm khô và hòa tan trong $50\ \mu\text{l}$ đệm TE.
  • Điều kiện chu nhiệt PCR: Thể tích phản ứng $15\ \mu\text{l}$. Chu trình nhiệt: Biến tính ban đầu 95°C trong 5 phút; 35 chu kỳ gồm [biến tính 95°C trong 30 giây, bắt cặp mồi ở nhiệt độ tối ưu $T_a$ trong 30 giây, kéo dài 72°C trong 1 phút]; kéo dài lần cuối 72°C trong 5 phút; bảo quản ở 4°C.

Chương 3: Kết quả và thảo luận

Chương 3 trình bày các số liệu thực nghiệm, kết quả tính toán tin sinh học và phân tích tiến hóa:

  • Khảo sát đặc tính lý hóa của mồi (In silico): Đánh giá qua IDT OligoAnalyzer xác nhận các cặp mồi đáp ứng đầy đủ tiêu chuẩn thiết kế:
    • Cặp ITS1/ITS4: Mồi ITS1 có chiều dài 19 bp, %GC = 63.2%, $T_m = 59.5^\circ\text{C}$; mồi ITS4 có chiều dài 20 bp, %GC = 45.0%, $T_m = 52.1^\circ\text{C}$; năng lượng liên kết cấu trúc thứ cấp $\Delta G \ge -9.82\ \text{kcal/mol}$. Mô phỏng trên Annhyb cho kích thước sản phẩm khuếch đại trên C. militaris đạt 567 bp (vị trí 2–568).
    • Cặp LR0R/LR5: Mồi LR0R có chiều dài 19 bp, %GC = 52.6%, $T_m = 53.3^\circ\text{C}$; mồi LR5 có chiều dài 17 bp, %GC = 47.1%, $T_m = 48.0^\circ\text{C}$; $\Delta G \ge -8.05\ \text{kcal/mol}$. Mô phỏng Annhyb cho kích thước khuếch đại 883 bp (vị trí 1–882 trên mã GQ249993.6).
    • Cặp ATP6-C1A/ATP6-C2A: Chiều dài lần lượt là 20 bp và 26 bp; $T_m$ biến thiên tương ứng 46.3–49.8°C và 55.9–58.3°C; $\Delta G \ge -10.01\ \text{kcal/mol}$. Mô phỏng Annhyb cho kích thước sản phẩm đạt 692 bp (vị trí 26252–26944 trên C. militaris).
  • Kết quả tách chiết và điện di sản phẩm PCR: DNA tổng số thu được từ các mẫu nấm đạt độ tinh sạch cao ($OD_{260}/OD_{280} \approx 1.8$). Điện di gel agarose 1.5% cho thấy các băng sản phẩm PCR của cả ba gen ITS, nrLSU và Atp6 đều đơn băng, sáng rõ, đúng kích thước dự đoán và không xuất hiện hiện tượng dimer hay khuếch đại không đặc hiệu.
  • Hiệu chỉnh trình tự và phân tích BLAST: Tín hiệu giải trình tự hai chiều được hiệu chỉnh bằng Chromas Lite và Dot plot. Kết quả tra cứu BLASTn trên NCBI cho thấy các đoạn trình tự sau hiệu chỉnh của mẫu DL0036, DL0091, DL0099 có độ tương đồng cao ($\ge 98-100%$) với các taxon thuộc chi Cordyceps và Ophiocordyceps.
  • Phân tích mô hình tiến hóa và xây dựng cây phát sinh loài:
    • Chức năng Find Best DNA Models trong MEGA X đã xác định các mô hình tiến hóa có điểm số BIC (Bayesian Information Criterion) thấp nhất cho từng tập dữ liệu đơn gen (nrLSU, ITS, Atp6) và tập dữ liệu đa gen kết hợp.
    • Cây phát sinh loài được dựng đồng thời bằng ba phương pháp NJ, ML và MP với 1000 lần lặp bootstrap. Cây đa gen kết hợp (nrLSU-ITS-Atp6 và nrLSU-ITS) thể hiện cấu trúc phân nhóm rõ ràng, giá trị bootstrap tại các nút phân nhánh chính đều đạt $\ge 50%$ đến $100%$, nâng cao độ tin cậy phân loại so với cây đơn gen. Kết quả định danh phân tử đồng nhất giữa các locus gen cho phép xếp các mẫu khảo sát vào các loài đại diện tương ứng trong chi Cordyceps s.l.

Kết quả và đóng góp

Kết quả chính

  1. Thiết lập quy trình chuẩn in silico: Đã kiểm tra và tối ưu hóa thành công các thông số lý hóa, kích thước sản phẩm và nhiệt độ bắt cặp cho 3 cặp mồi mã vạch (ITS1/ITS4, LR0R/LR5, ATP6-C1A/ATP6-C2A) với kích thước khuếch đại lý thuyết lần lượt là 567 bp, 883 bp và 692 bp.
  2. Thu nhận và khuếch đại thành công DNA chất lượng cao: Tách chiết thành công DNA tổng số từ mô nấm ký sinh tự nhiên bằng phương pháp Phenol/Chloroform có bổ sung CTAB và $\beta$-mercaptoethanol, khắc phục triệt để trở ngại từ vách polysaccharide/chitin của nấm. Phản ứng PCR cho sản phẩm khuếch đại đặc hiệu cho cả ba vùng gen.
  3. Xây dựng bộ dữ liệu mã vạch DNA: Giải trình tự, xử lý tín hiệu chromatogram và hoàn thiện tập dữ liệu DNA barcode cho 3 mẫu nấm thực địa (DL0036, DL0091, DL0099) trên cả ba locus ITS, nrLSU và Atp6.
  4. Định danh phân tử dựa trên cây phát sinh loài đa gen: Xác lập cây phả hệ phân tử bằng 3 thuật toán (NJ, ML, MP) với giá trị ủng hộ bootstrap cao, chứng minh hiệu quả của việc kết hợp đa gen (nrLSU-ITS-Atp6) trong việc tăng số lượng vị trí nucleotide phân tích, loại bỏ sai số phân loại đơn gen.

Kiến nghị của tác giả

  • Tiếp tục mở rộng thu thập thêm các mẫu nấm ký sinh côn trùng tại các khu vực rừng tự nhiên khác nhau của vùng núi Langbiang và cao nguyên Tây Nguyên nhằm làm phong phú bộ sưu tập nguồn gen.
  • Bổ sung thêm các locus gen mã hóa protein khác như Tef1 (elongation factor 1$\alpha$), Rpb1, Rpb2 (RNA polymerase II subunits) và $\beta$-tubulin (Tub) để hoàn thiện cây phát sinh loài ở mức độ siêu cấu trúc theo hệ thống phân loại toàn cầu.

Đóng góp khoa học

  • Cung cấp dữ liệu trình tự DNA barcode chuẩn xác cho các mẫu nấm Cordyceps s.l thu nhận tại vùng núi Langbiang - Lâm Đồng, bổ sung vào cơ sở dữ liệu nguồn gen nấm dược liệu của Việt Nam.
  • Khẳng định tính ưu việt và độ tin cậy của phương pháp phân tích đa gen giữ nhà kết hợp giữa gen nhân (nrLSU, ITS) và gen ty thể (Atp6) trong định danh nấm ký sinh côn trùng, khắc phục các hạn chế của phân loại hình thái truyền thống.

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

Hạn chế của nghiên cứu

  • Số lượng mẫu khảo sát thực nghiệm còn hạn chế (3 mẫu: DL0036, DL0091, DL0099), chưa đại diện cho toàn bộ tính đa dạng sinh học của chi Cordyceps s.l tại vùng núi Langbiang.
  • Nghiên cứu tập trung vào 3 vùng gen (nrLSU, ITS, Atp6); chưa kết hợp đầy đủ bộ 5–7 gen giữ nhà (nrSSU, Tef1, Rpb1, Rpb2, Tub) như các tiêu chuẩn phân loại quốc tế mở rộng của Sung et al. (2007).
  • Đề tài chỉ dừng lại ở bước định danh phân tử và phân tích phát sinh loài, chưa đi sâu phân tích thành phần hóa học hay đánh giá hoạt tính dược lý của các mẫu nấm thu nhận.

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

  • Mở rộng phạm vi thu thập mẫu nấm thực địa theo mùa và theo các đai cao khác nhau tại Lâm Đồng.
  • Phối hợp phân tích vi học giải phẫu chi tiết (đặc điểm nang, bào tử nang, thể quả) kết hợp song song với mã vạch đa gen để xây dựng khóa định danh tích hợp hình thái học - sinh học phân tử hoàn chỉnh.
  • Tiến hành phân lập thuần khiết, nhân nuôi sinh khối hệ sợi trong điều kiện nhân tạo và phân tích định lượng các hoạt chất sinh học quý (cordycepin, adenosine, polysaccharide) phục vụ ứng dụng y dược.

Giá trị tham khảo

  • Đối tượng tham khảo phù hợp:
    • Sinh viên, học viên cao học và nghiên cứu sinh chuyên ngành Công nghệ Sinh học, Sinh học Phân tử, Nấm học (Mycology) và Dược liệu học.
    • Các cán bộ nghiên cứu tại các viện, trung tâm công nghệ sinh học đang thực hiện các đề tài về phân loại vi sinh vật, bảo tồn nguồn gen quý và xây dựng dữ liệu mã vạch DNA.
  • Các phần nội dung có giá trị tham khảo cao:
    • Quy trình tách chiết DNA nấm khó ly giải: Phương pháp phối hợp Phenol/Chloroform với CTAB và $\beta$-mercaptoethanol xử lý hiệu quả mô nấm giàu polysaccharide và proteinase.
    • Phương pháp luận In Silico: Quy trình kiểm tra, thẩm định các thông số mồi, tính toán cấu trúc thứ cấp ($\Delta G$, dimer, hairpin) và mô phỏng phản ứng PCR trên phần mềm chuyên dụng.
    • Quy trình phân tích phát sinh loài: Các bước xử lý file giải trình tự bằng Chromas Lite, liên kết mạch, căn chỉnh đa trình tự và lựa chọn mô hình tiến hóa nucleotide tối ưu trên phần mềm MEGA X để dựng cây NJ, ML, MP.

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

1. Tại sao nghiên cứu lại lựa chọn kết hợp ba vùng gen nrLSU, ITS và Atp6 thay vì chỉ sử dụng một vùng gen duy nhất?
Trả lời: Vùng gen ITS có tốc độ tiến hóa nhanh, rất nhạy để phân biệt ở cấp độ loài và dưới loài nhưng kém tin cậy ở cấp độ chi/họ. Ngược lại, vùng gen nrLSU có tính bảo tồn cao, phù hợp để xác định quan hệ ở cấp độ họ và bộ. Gen ty thể Atp6 mã hóa cho protein thiết yếu, có tần số tích lũy đột biến riêng biệt giúp phân định các nhánh gần nhau trong họ Clavicipitaceae. Việc kết hợp cả ba locus giúp tăng tổng số vị trí nucleotide đối sánh, làm tăng giá trị bootstrap và độ tin cậy cấu trúc của cây phả hệ phân tử.

2. Thành phần nào trong quy trình tách chiết đóng vai trò quyết định giúp thu nhận DNA nấm đạt chất lượng cao?
Trả lời: Vách tế bào nấm Ascomycota cấu tạo từ mạng lưới chitin và glucan phức tạp liên kết với protein màng. Việc bổ sung chất hoạt động bề mặt CTAB giúp phá vỡ màng lipid và tạo phức kết tủa polysaccharide, trong khi $\beta$-mercaptoethanol có vai trò cắt đứt các cầu nối disulfide trong cấu trúc bậc 3 và bậc 4 của protein, hỗ trợ Proteinase K phân giải triệt để protein và ức chế nuclease, giúp DNA thu được sạch và không bị đứt gãy.

3. Tiêu chí nào được sử dụng để đánh giá một cặp mồi đạt yêu cầu trong phân tích in silico trước khi tiến hành PCR?
Trả lời: Các thông số vật lý của mồi phải thỏa mãn: chiều dài từ 17–28 bp; hàm lượng %GC nằm trong khoảng 40–60%; nhiệt độ nóng chảy $T_m$ trong khoảng 50–65°C và chênh lệch $T_m$ giữa hai mồi không quá 5°C; đầu $3'$ có mặt G hoặc C để tăng độ bám đặc hiệu; trị tuyệt đối năng lượng tự do hình thành cấu trúc kẹp tóc (hairpin), tự bắt cặp (self-dimer) và bắt cặp chéo (hetero-dimer) phải có $\Delta G \ge -9\ \text{kcal/mol}$.

4. Ý nghĩa của giá trị bootstrap trong đánh giá cây phát sinh loài được thể hiện như thế nào trong khóa luận?
Trả lời: Giá trị bootstrap (tính theo tỷ lệ %) thể hiện tần suất xuất hiện của một nhánh phân loại khi lặp lại ngẫu nhiên tập dữ liệu 1000 lần. Giá trị bootstrap $\ge 50%$ thể hiện sự phân nhóm được dữ liệu ủng hộ, và khi đạt $\ge 70%$ (tương đương xác suất thống kê 95%) thì mối quan hệ phân loại đơn ngành giữa các taxon trong nhánh đó đạt mức tin cậy khoa học rất cao.

Kết luận

Khóa luận tốt nghiệp của tác giả Võ Thanh Nhàn đã hoàn thành việc xây dựng bộ dữ liệu mã vạch DNA cho các mẫu nấm ký sinh côn trùng thuộc chi Cordyceps s.l thu nhận tại vùng núi Langbiang (Lâm Đồng) trên ba locus gen chỉ thị nrLSU, ITS và Atp6. Đề tài đã chuẩn hóa thành công quy trình thực nghiệm từ phân tích in silico, tách chiết ly giải vách tế bào nấm, tối ưu phản ứng PCR đến xử lý tín hiệu giải trình tự Sanger hai chiều. Kết quả xây dựng cây phát sinh loài đa gen bằng các phương pháp NJ, ML và MP đã xác định chính xác vị trí phân loại phân tử của các mẫu nấm khảo sát, đóng góp cơ sở dữ liệu quan trọng phục vụ công tác giám định nguồn gen và định hướng nghiên cứu nấm dược liệu tại Việt Nam.