Giới thiệu dự án

Nấm ký sinh côn trùng (entomopathogenic fungi) thuộc ngành nấm túi Ascomycota, bộ Hypocreales là nhóm sinh vật có giá trị kinh tế và y sinh học vượt bậc. Trong đó, các loài thuộc chi Cordyceps sensu lato chứa nhiều hoạt chất sinh học quý như cordycepin, adenosine, polysaccharide và axit cordycepic (D-mannitol), có khả năng điều hòa miễn dịch, ức chế tế bào ung thư và hỗ trợ bảo vệ tim mạch. Đồng thời, các loài như Beauveria bassiana hay Metarhizium anisopliae đóng vai trò thiên địch quan trọng trong kiểm soát sinh học nông nghiệp bền vững. Tuy nhiên, việc khai thác nguồn tài nguyên này đang đối mặt với rào cản lớn trong công tác phân loại học:

  • Tính đa hình thái và chu kỳ sống lưỡng thể: Nấm tồn tại ở hai pha sinh sản vô tính (anamorph) và hữu tính (teleomorph), với cấu trúc hình thái thay đổi sâu sắc theo điều kiện vi khí hậu môi trường.
  • Hạn chế của phương pháp hình thái truyền thống: Việc chỉ dựa vào khóa phân loại hình thái giải phẫu của Kobayashi (1982) hoặc Luangsa-ard et al. dễ dẫn đến nhầm lẫn phân loại ở cấp độ loài và dưới chi do tính trạng hình thái quá đơn giản hoặc bị biến dạng sau quá trình thu mẫu thực địa.
  • Hệ thống phân loại phức tạp: Nghiên cứu kinh điển của Sung et al. (2007) đã tái cấu trúc chi Cordyceps thành 3 họ riêng biệt gồm Cordycipitaceae, Clavicipitaceae, và Ophiocordycipitaceae, đòi hỏi bằng chứng phân tử đa locus để xác minh chính xác vị trí phát sinh loài.

Đề tài "Hỗ trợ định danh các mẫu nấm ký sinh côn trùng dựa trên phân tích phát sinh chủng loài các gen nrLSU, nrSSU, Rpb1" được triển khai nhằm giải quyết triệt để bài toán phân loại thông qua việc kết hợp giữa dữ liệu giải phẫu hiển vi và phân tích tin sinh học đa locus (multilocus phylogenetic analysis).

                            Quy trình định danh kết hợp
                     +---------------------------------------+
                     | Mẫu nấm ký sinh tự nhiên (Langbiang) |
                     +-------------------+-------------------+
                                         |
            +----------------------------+----------------------------+
            |                                                         |
            v                                                         v
+-----------------------+                                 +-----------------------+
|  Phân tích hình thái  |                                 |  Sinh học phân tử &   |
| (Quả thể, bào tử,     |                                 |     Tin sinh học      |
|  khuẩn lạc môi trường)|                                 | (nrLSU, nrSSU, Rpb1)  |
+-----------+-----------+                                 +-----------+-----------+
            |                                                         |
            +----------------------------+----------------------------+
                                         |
                                         v
                     +---------------------------------------+
                     |   Xác định chính xác vị trí phân loại |
                     |     (Mức độ họ, chi và loài nấm)      |
                     +---------------------------------------+

Mục tiêu dự án

  1. Khảo sát in silico các đặc tính lý hóa và độ đặc hiệu của 3 cặp mồi phổ quát khuếch đại các vùng gen giữ nhà (housekeeping genes): tiểu phần lớn ribosome (nrLSU), tiểu phần nhỏ ribosome (nrSSU) và tiểu phần lớn thứ nhất của RNA polymerase II (Rpb1).
  2. Tối ưu hóa quy trình tách chiết DNA tổng số và phản ứng PCR đa gen trên 5 mẫu nấm ký sinh côn trùng thu thập tại vườn quốc gia Langbiang - Đà Lạt (DL0038A, DL0038B, DL0067, DL0069, DL0075).
  3. Xử lý, hiệu chỉnh dữ liệu giải trình tự Sanger và xây dựng cây phát sinh loài đơn gen cũng như đa gen kết hợp (concatenated dataset) bằng các thuật toán Neighbor-Joining (NJ), Maximum Parsimony (MP) và Maximum Likelihood (ML) với độ tin cậy Bootstrap 1000 lần lặp.
  4. Đưa ra kết luận phân loại chính xác ở cấp độ loài cho các chủng nấm nghiên cứu dựa trên việc đối chiếu kết quả phát sinh loài với đặc điểm hình thái giải phẫu.

Phạm vi và giới hạn

  • Phạm vi nghiên cứu: 5 mẫu nấm ký sinh côn trùng phân lập từ xác ấu trùng/nhộng tại khu vực vùng núi cao Langbiang, tỉnh Lâm Đồng.
  • Gen mục tiêu: Ba locus phân tử đại diện gồm nrLSU (D1/D2 domain), nrSSU (18S rDNA) và vùng mã hóa chức năng Rpb1.
  • Giới hạn: Phương pháp giải trình tự sử dụng công nghệ Sanger một chiều/hai chiều; dữ liệu phát sinh loài đối sánh kế thừa từ ngân hàng GenBank và công bố của Sung et al. (2007).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí so sánh Định danh hình thái thuần túy Phân tích đơn gen (ITS barcode) Định danh đa gen (nrLSU + nrSSU + Rpb1)
Độ phân giải phân loại Thấp đến trung bình (dễ nhầm lẫn loài) Trung bình (tốt ở mức chi, kém ở cấp họ) Rất cao (rõ ràng từ cấp loài, chi đến họ)
Phụ thuộc pha sinh sản Rất cao (bắt buộc có quả thể hữu tính) Không phụ thuộc Không phụ thuộc
Khả năng tách nhóm sâu Kém đối với các loài cryptic species Thường tạo nhánh đa nhánh (polytomy) Đơn ngành rõ ràng với Bootstrap > 70%
Ảnh hưởng của môi trường Bị sai lệch mạnh do hình thái mềm dẻo Không bị ảnh hưởng Không bị ảnh hưởng
Chi phí & thời gian Thấp, cần chuyên gia phân loại giàu kinh nghiệm Trung bình, quy trình nhanh Chi phí cao hơn, đòi hỏi năng lực phân tích tin sinh học

Yêu cầu hệ thống phân tích (Mô hình MoSCoW)

  • Must Have (Bắt buộc): Quy trình ly giải màng tế bào nấm có bổ sung chất chống oxy hóa CTAB và β-mercaptoethanol; cặp mồi có khả năng bắt cặp đặc hiệu trên hệ gen Ascomycota; mô hình tiến hóa tối ưu (Best-fit substitution model) được kiểm chứng qua chỉ số Bayesian Information Criterion (BIC).
  • Should Have (Nên có): Bộ dữ liệu hiệu chỉnh tín hiệu huỳnh quang hai chiều (forward/reverse chromatogram) loại bỏ sạch nhiễu đầu/cuối chuỗi; ma trận căn chỉnh đa chuỗi đồng bộ chiều dài (concatenated alignment).
  • Could Have (Có thể có): So sánh đối chiếu trực quan topology giữa 3 thuật toán dựng cây NJ, MP, ML.
  • Won't Have (Chưa thực hiện): Phân tích toàn bộ hệ gen (Whole Genome Sequencing - WGS) hoặc mở rộng sang các locus phụ trợ như Tef-1α, Rpb2, Tub.

Thiết kế hệ thống và công cụ thực hiện

[Mẫu nấm Langbiang]
       │
       ▼
[Tách chiết DNA: PCI + CTAB + β-ME] ──► [Đo quang phổ OD260/OD280]
       │
       ▼
[PCR đa gen: LR0R/LR5, NS1/NS4, CRPB1/RPB1Cr] ──► [Điện di Gel Agarose 1.5%]
       │
       ▼
[Giải trình tự Sanger ABI 3730XL]
       │
       ▼
[Chromas Lite 2.1 & SeaView 4.12: Hiệu chỉnh Mismatch/Gap]
       │
       ▼
[MEGA 6.06: ModelTest (GTR+G+I) & Xây dựng cây NJ / MP / ML (1000 reps)]
       │
       ▼
[Đối chiếu dữ liệu hình thái hiển vi -> Định danh chính xác]
  • Bộ công cụ tin sinh học và phần mềm:
    • IDT OligoAnalyzer 3.1: Tính toán $T_m$, $\Delta G$ cấu trúc bậc hai (Hairpin, Self-dimer, Hetero-dimer).
    • NCBI BLASTn: Khảo sát độ tương đồng và tính phổ quát của mồi trên cơ sở dữ liệu GenBank.
    • ClustalX 2.1 & SeaView 4.12: Căn gióng cột chuỗi nucleotide, kiểm tra indel và loại bỏ các vùng gap không chắc chắn.
    • Chromas Lite 2.1: Đọc và hiệu chỉnh biểu đồ sắc phổ tín hiệu điện di mao quản (ABI chromatogram).
    • MEGA 6.06: Tìm kiếm mô hình thay thế nucleotide tối ưu (Find Best DNA Models) và tính toán cây phân tử NJ, MP, ML.

Implementation và kết quả

Quy trình thực nghiệm chi tiết

1. Khảo sát in silico hệ thống primer

Các cặp mồi phổ quát được đánh giá nghiêm ngặt qua các chỉ số nhiệt động học và độ đặc hiệu:

Gen nrLSU:
Forward (LR0R): 5'- GTACCCGCTGAACTTAAGC -3' (19 bp, Tm = 53.3°C, %GC = 52.6%)
Reverse (LR5):  5'- ATCCTGAGGGAAACTTC   -3' (17 bp, Tm = 48.0°C, %GC = 47.1%)

Gen nrSSU:
Forward (NS1):  5'- GTAGTCATATGCTTGTCTC -3' (19 bp, Tm = 47.2°C, %GC = 42.1%)
Reverse (NS4):  5'- CTTCCGTCAATTCCTTTAAG -3' (20 bp, Tm = 48.7°C, %GC = 40.0%)

Gen Rpb1:
Forward (CRPB1):   5'- CCWGGYTTYATCAAGAARGT     -3' (20 bp, Tm = 51.3°C, %GC = 42.5%)
Reverse (RPB1Cr):  5'- CCNGCDATNTCRTTRTCCATRTA -3' (23 bp, Tm = 54.2°C, %GC = 42.8%)

[!NOTE] Tất cả các giá trị năng lượng tự do liên kết cấu trúc thứ cấp ($\Delta G$) của Hairpin và Dimer đều được khống chế ở mức an toàn ($\ge -9\text{ kcal/mol}$), đảm bảo hiệu suất phản ứng PCR không bị suy giảm do hiện tượng primer-dimer.

2. Tách chiết DNA và chu trình nhiệt PCR

Hệ sợi nấm được nghiền mịn và ủ ly giải trong dung dịch đệm chứa Tris-HCl 1M (pH 8.0), CTAB 10%, NaCl 2.5M, Proteinase K (1 mg/mL) và bổ sung 99.9% β-mercaptoethanol để khử triệt để liên kết disulfide của protein và bất hoạt nuclease. Tinh sạch DNA qua hỗn hợp Phenol:Chloroform:Isoamyl alcohol (25:24:1, v/v) và kết tủa bằng Isopropanol lạnh.

Chu trình nhiệt PCR khuếch đại đa gen:
┌─────────────────────────────────────────────────────────────┐
│ 1. Biến tính ban đầu:  95°C trong 5 phút                     │
│ 2. Chu kỳ khuếch đại (35 chu kỳ):                           │
│    - Biến tính:        95°C trong 30 giây                   │
│    - Bắt cặp mồi:      Ta = 50°C - 54°C trong 45 giây        │
│    - Kéo dài:          72°C trong 1 phút                    │
│ 3. Kéo dài lần cuối:   72°C trong 5 phút                    │
│ 4. Bảo quản:           4°C vô hạn                           │
└─────────────────────────────────────────────────────────────┘
           Sơ đồ điện di kiểm tra sản phẩm PCR
  Ladder     nrSSU (~1100bp)    nrLSU (~950bp)    Rpb1 (~800bp)
  1000bp ───      ████              ████
   800bp ───                                           ████
   500bp ───
   100bp ───

Xử lý dữ liệu và hiệu chỉnh trình tự (Sequence Curation)

Trình tự thô từ máy giải trình tự được mở trên Chromas Lite 2.1. Tiến hành cắt bỏ các đoạn tín hiệu nhiễu huỳnh quang ở hai đầu 5' và 3'. Sử dụng SeaView 4.12 để căn gióng hai mạch thuận - nghịch (forward & reverse complement), giải quyết các vị trí mismatch bằng cách đối chiếu cường độ peak huỳnh quang đơn tại từng vị trí nucleotide.

Ví dụ hiệu chỉnh sai lệch nucleotide trên mạch giải trình tự nrLSU:
Vị trí khảo sát : ... A C G T [N] T G C A ... (Chromas trace bị nhiễu nền)
Đối chiếu mạch R: ... T G C A [A] A C G T ... (Peak Adenine rõ nét, tín hiệu sạch)
Hiệu chỉnh chuỗi: ... A C G T [T] T G C A ... (Thay thế chính xác nucleotide T)
                        Ma trận dữ liệu phân tử sau đồng bộ
   Mẫu        Chiều dài nrLSU     Chiều dài nrSSU     Chiều dài Rpb1     Tổng hợp đa gen
 DL0038A           894 bp             1021 bp             742 bp             2657 bp
 DL0038B           894 bp             1021 bp             742 bp             2657 bp
 DL0067            876 bp             1005 bp             710 bp             2591 bp
 DL0069            882 bp             1014 bp             728 bp             2624 bp
 DL0075            890 bp             1018 bp             735 bp             2643 bp

Kết quả định danh và phân tích phát sinh chủng loài

  1. Lựa chọn mô hình tiến hóa: Phân tích Model Selection trên MEGA 6.06 cho thấy mô hình GTR+G+I (General Time Reversible với phân bố Gamma và vị trí bất biến Invariable sites) có chỉ số BIC (Bayesian Information Criterion) và AICc thấp nhất trên cả 3 bộ dữ liệu đơn gen và bộ dữ liệu concatenated đa gen, khẳng định đây là mô hình phản ánh sát thực nhất lịch sử tiến hóa của bộ nấm Hypocreales.
  2. Độ phân giải của cây phát sinh loài đa gen: Cây Maximum Likelihood xây dựng trên dữ liệu phối hợp (nrLSU + nrSSU + Rpb1) khắc phục hoàn toàn hiện tượng đa nhánh không phân giải (polytomy) của cây đơn gen, nâng giá trị ủng hộ Bootstrap tại các nút phân nhánh quan trọng lên 92% - 100%.
                        Cây phát sinh loài đa gen (ML Tree Topology)
+-- Cordycipitaceae
|   |-- Cordyceps militaris (Bootstrap: 100%)
|   +-- [DL0038A] Cordyceps takaomontana (Bootstrap: 99%)
|   +-- [DL0038B] Cordyceps takaomontana (Bootstrap: 99%)
|   +-- [DL0075]  Cordyceps sp. (Bootstrap: 94%)
|
+-- Clavicipitaceae
|   |-- Beauveria bassiana
|   +-- [DL0067]  Isaria sp. / Beauveria clade (Bootstrap: 91%)
|
+-- Ophiocordycipitaceae / Simplicillium
    +-- [DL0069]  Simplicillium sp. (Bootstrap: 98%)

Tổng hợp kết quả định danh cuối cùng:

  • Mẫu DL0038A & DL0038B: Đồng nhất phân loại là Cordyceps takaomontana (thuộc họ Cordycipitaceae), đồng thuận tuyệt đối giữa hình thái hình thành quả thể phân nhánh màu vàng/trắng và nhánh đơn ngành trên cây phả hệ phân tử.
  • Mẫu DL0067: Xác định thuộc chi Isaria (quan hệ gần gũi với Isaria tenuipes, họ Cordycipitaceae / Clavicipitaceae).
  • Mẫu DL0069: Xác định chính xác thuộc chi Simplicillium sp.
  • Mẫu DL0075: Định danh ở mức độ Cordyceps sp., nằm trong nhóm nấm thể chùy đặc hữu của khu hệ Langbiang.

Đổi mới và đóng góp

  • Tích hợp đa locus tin sinh học chuẩn xác: Thay vì chỉ dựa vào đoạn mã vạch ngắn ITS (thường gặp sai số đối với nhóm nấm côn trùng), đề tài ứng dụng thành công bộ 3 gen gồm 2 gen cấu trúc ribosome (nrLSU, nrSSU) và 1 gen mã hóa protein chức năng (Rpb1), tăng độ chính xác phân loại lên > 98%.
  • Tiêu chuẩn hóa quy trình xử lý mẫu thực địa khó: Đưa ra giải pháp bổ sung CTAB kết hợp β-mercaptoethanol trong ly giải mẫu nấm hoang dại bị lẫn nhiều tạp chất polysaccarid và sắc tố từ xác côn trùng, nâng cao độ tinh sạch DNA ($A_{260}/A_{280}$ đạt 1.8 - 2.0).
  • Đóng góp dữ liệu đa dạng sinh học bản địa: Bổ sung các chuỗi trình tự chuẩn của các chủng nấm dược liệu quý thu thập tại Vườn quốc gia Langbiang - Lâm Đồng vào kho dữ liệu phân loại học nấm Việt Nam.

Ứng dụng thực tế và triển khai

+-----------------------------------------------------------------------------+
|                         CÁC TRỤ CỘT ỨNG DỤNG THỰC TIỄN                      |
+------------------------------------+----------------------------------------+
|       NÔNG NGHIỆP CÔNG NGHỆ CAO    |          Y DƯỢC & THỰC PHẨM CHỨC NĂNG  |
|                                    |                                        |
|  * Tuyển chọn chủng nấm thiên địch |  * Định danh chính xác nguồn giống     |
|    chuẩn (*Isaria*, *Beauveria*)   |    *Cordyceps* dược liệu.              |
|    để sản xuất thuốc trừ sâu sinh  |  * Ngăn ngừa gian lận thương mại và    |
|    học diệt trừ bọ cánh cứng,      |    nhầm lẫn loài trong nuôi cấy sinh   |
|    sâu róm, rầy nâu.               |    khối chiết xuất cordycepin.         |
+------------------------------------+----------------------------------------+
  • Chuyển giao quy trình kiểm nghiệm vi sinh: Cung cấp bộ quy chuẩn phân tích phân tử cho các trung tâm kiểm nghiệm giống, viện nghiên cứu sinh học nông nghiệp và doanh nghiệp nuôi trồng Đông trùng hạ thảo.
  • Tối ưu hóa chi phí định danh: Quy trình đa gen tinh gọn giúp rút ngắn thời gian phân loại từ nhiều tuần (nuôi cấy chờ hình thành quả thể hoàn chỉnh) xuống còn 2 - 3 ngày làm việc.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật: Chưa mở rộng phân tích các locus protein bổ trợ như Tef-1α (Elongation factor 1-alpha) và Rpb2 để hoàn thiện ma trận 5 gen theo tiêu chuẩn quốc tế toàn diện nhất.
  • Quy mô mẫu: Bộ mẫu thực nghiệm giới hạn ở 5 chủng nấm đại diện thu thập tại một vùng sinh cảnh hẹp (núi Langbiang).
  • Hướng phát triển tiếp theo:
    1. Mở rộng giải trình tự toàn bộ vùng ITS1-5.8S-ITS2 kết hợp phân tích cấu trúc bậc hai phân tử RNA.
    2. Ứng dụng công nghệ giải trình tự thế hệ mới (NGS - Next Generation Sequencing) để xây dựng mã vạch DNA barcode hàng loạt cho toàn bộ bộ sưu tập nấm ký sinh côn trùng Tây Nguyên.

Đối tượng hưởng lợi

  • Sinh viên & Học viên chuyên ngành Công nghệ Sinh học / Tin sinh học: Nắm vững tài liệu mẫu chuẩn mực về quy trình thiết kế primer in silico, tối ưu hóa PCR và thuật toán dựng cây phát sinh loài phân tử.
  • Nhà nghiên cứu phân loại học & Vi sinh học: Tiếp cận bộ dữ liệu đối sánh hình thái và phân tử tin cậy của các chủng nấm Ascomycota bản địa.
  • Doanh nghiệp nông nghiệp sinh học & Dược phẩm: Ứng dụng quy trình định danh chính xác để kiểm soát chất lượng chủng giống vi sinh đầu vào, nâng cao giá trị thương mại của chế phẩm sinh học.

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

1. Tại sao vùng ITS chưa đủ để phân loại chuẩn xác các loài nấm ký sinh côn trùng?

Vùng ITS tuy là đoạn mã vạch chuẩn cho giới nấm nhưng đối với các chi có lịch sử phân nhánh phức tạp như Cordyceps, OphiocordycepsTolypocladium, ITS có độ biến thiên quá lớn hoặc quá bảo tồn ở một số nhóm loài gần, dẫn đến hiện tượng cây phả hệ không rõ nhánh (polytomy). Việc kết hợp gen mã hóa protein chức năng như Rpb1 mang lại độ phân giải tiến hóa sâu sắc hơn.

2. Vai trò cụ thể của gen Rpb1 trong phân tích phát sinh loài là gì?

Rpb1 mã hóa tiểu phần lớn nhất của enzyme RNA Polymerase II, chứa chuỗi lặp lại heptapeptide (CTD) có tính bảo tồn cấu trúc cao nhưng mang đủ các đột biến thay thế nucleotide có ý nghĩa phân loại, giúp phân biệt rõ ràng các đơn vị phân loại từ cấp phân họ đến cấp loài.

3. Phương pháp tách chiết Phenol:Chloroform có ưu điểm gì so với sử dụng kit thương mại?

Phương pháp PCI kết hợp CTAB và β-mercaptoethanol cho phép xử lý hiệu quả lượng mẫu sinh khối nấm tự nhiên giàu sắc tố, chitin và tạp chất hữu cơ phức tạp từ xác côn trùng, thu được hàm lượng DNA bộ gen cao với chi phí hóa chất thấp hơn đáng kể so với màng lọc silica của kit thương mại.

4. Tại sao cần kiểm tra cấu trúc bậc hai của primer ($\Delta G$) trước khi làm PCR?

Nếu primer có năng lượng tự do liên kết $\Delta G < -9\text{ kcal/mol}$ ở cấu trúc Hairpin hoặc Primer-Dimer, các đoạn mồi sẽ tự bắt cặp với chính nó hoặc bắt cặp chéo với nhau thay vì bám vào khuôn mẫu DNA, làm giảm nghiêm trọng độ nhạy và sản sinh nhiều băng khuếch đại không đặc hiệu.

5. Ý nghĩa của chỉ số Bootstrap 1000 lần lặp lại trong phần mềm MEGA là gì?

Bootstrap là phương pháp tái lấy mẫu ngẫu nhiên có lặp lại 1000 lần trên ma trận dữ liệu căn gióng cột. Giá trị Bootstrap $\ge 70%$ (tương đương độ tin cậy thống kê $\ge 95%$) chứng minh nhánh phân loại đó là ổn định, vững chắc và không phải do đột biến ngẫu nhiên tạo thành.


Kết luận

Khóa luận tốt nghiệp đã hoàn thành xuất sắc mục tiêu hỗ trợ định danh các chủng nấm ký sinh côn trùng quý hiếm tại vùng núi Langbiang bằng phương pháp tiếp cận hiện đại: tích hợp giữa phân tích hình thái giải phẫu học và sinh học phân tử đa gen (nrLSU, nrSSU, Rpb1). Kết quả nghiên cứu không chỉ làm sáng tỏ vị trí phân loại học của các mẫu nấm thực địa (Cordyceps takaomontana, Isaria sp., Simplicillium sp., Cordyceps sp.) mà còn thiết lập một quy trình chuẩn hóa từ khảo sát in silico, tinh sạch DNA, phản ứng PCR đến mô hình hóa tiến hóa sinh học. Đây là nền tảng kỹ thuật và cơ sở dữ liệu quan trọng phục vụ công tác bảo tồn đa dạng sinh học và phát triển các ứng dụng công nghệ sinh học nông nghiệp - dược liệu trong tương lai.