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

Khoai tây (Solanum tuberosum L.) thuộc họ Cà (Solanaceae), là cây lương thực ngắn ngày lấy củ có giá trị dinh dưỡng cao, tiềm năng kinh tế và xuất khẩu lớn. Trên thế giới, khoai tây là cây trồng quan trọng thứ 5 sau lúa gạo, ngô, lúa mì và đậu tương. Tại Việt Nam, khoai tây được canh tác phổ biến trong vụ đông ở các tỉnh phía Bắc, đóng vai trò chủ lực trong cơ cấu cây màu lương thực. Để nâng cao năng suất và chất lượng khoai tây một cách bền vững, hướng tiếp cận sinh học dựa trên hệ vi sinh vật nội sinh (endophytic bacteria) đang được quan tâm nghiên cứu sâu rộng.

Vi khuẩn nội sinh là nhóm vi sinh vật cư trú bên trong các mô thực vật mà không gây ra bất kỳ triệu chứng bệnh lý hoặc tác động tiêu cực tức thời nào cho cây chủ. Nhóm vi sinh vật này mang lại nhiều lợi ích quan trọng cho cây trồng thông qua các cơ chế: cố định đạm sinh học, sinh tổng hợp hormone kích thích sinh trưởng thực vật (như Indole-3-acetic acid - IAA), hòa tan các hợp chất photpho khó tan, sản sinh siderophore giúp hấp thu sắt và đối kháng với các tác nhân gây bệnh thực vật (như nấm Phytophthora infestans gây bệnh mốc sương).

Trước đó, các mẫu vi khuẩn nội sinh đã được phân lập thành công từ rễ cây khoai tây và một số chủng đã được đánh giá bước đầu về tính đa dạng di truyền. Tuy nhiên, việc định danh chính xác ở mức độ loài cho các chủng phân lập này vẫn chưa được thực hiện hoàn chỉnh. Do đó, khóa luận tốt nghiệp tập trung vào đề tài: "Identification and sequence analysis of partial 16S rRNA of endophytic bacteria isolated from potato root" (Định danh và phân tích trình tự đoạn gen 16S rRNA của vi khuẩn nội sinh phân lập từ rễ khoai tây).

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

  1. Định danh và phân tích trình tự một phần đoạn gen 16S rRNA của các chủng vi khuẩn nội sinh phân lập từ rễ cây khoai tây.

Nhiệm vụ nghiên cứu cụ thể

  1. Định danh chính xác 12 chủng vi khuẩn nội sinh được phân lập từ mô rễ khoai tây.
  2. Phân tích đa dạng trình tự và dự đoán cấu trúc bậc hai của chuỗi 16S rRNA vi sinh vật nội sinh bằng các công cụ tin sinh học R2DT và UnaFold.

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

  • Đối tượng nghiên cứu: 12 mẫu vi khuẩn nội sinh phân lập từ rễ khoai tây (ký hiệu từ isolate_1 đến isolate_12).
  • Phạm vi không gian: Bộ môn Công nghệ sinh học Thực vật, Khoa Công nghệ sinh học, Học viện Nông nghiệp Việt Nam.
  • Thời gian hoàn thành: Năm 2022 (khóa luận bảo vệ tháng 09/2022).

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

Khung lý thuyết và khái niệm khoa học

  • Khái niệm vi khuẩn nội sinh: Được định nghĩa theo Bacon và White (2000), vi khuẩn nội sinh thực vật là những vi khuẩn cư trú bên trong mô sống của cây mà không gây ra các tác hại tức thời rõ rệt, thiết lập mối quan hệ cộng sinh hỗ trợ dinh dưỡng và bảo vệ cây chủ khỏi các tác nhân gây stress sinh học và phi sinh học.
  • Kỹ thuật giải trình tự gen 16S rRNA: Gen mã hóa 16S rRNA là thành phần có tính bảo thủ cao trong bộ máy phiên mã của mọi sinh vật nhân sơ (Prokaryote). Cấu trúc của gen 16S rRNA gồm các vùng bảo thủ xen kẽ với các vùng biến đổi cao (variable regions). Các cặp mồi phổ quát (universal primers) được thiết kế bám vào vùng bảo thủ để nhân bản gen qua PCR, trong khi việc phân tích trình tự ở vùng biến đổi cho phép phân loại và định danh chính xác đến mức chi và loài.
  • Các vùng chỉ thị phân tử đối chứng trong vi sinh: Bên cạnh 16S rRNA dùng cho vi khuẩn, khóa luận cũng hệ thống hóa các chỉ thị như đoạn ITS (Internal Transcribed Spacer), gen 18S rRNA và vùng D1/D2 của gen 26S rRNA thường dùng cho định danh nấm và sinh vật nhân thực.

Phương pháp nghiên cứu thực nghiệm và tin sinh học

Nghiên cứu sử dụng kết hợp giữa kỹ thuật sinh học phân tử thực nghiệm và các công cụ phân tích tin sinh học:

  1. Định lượng và kiểm tra độ sạch DNA: Sử dụng máy quang phổ đo độ hấp thụ (BioPhotometer plus) với cuvette Hellma traycell ở bước sóng $\lambda = 260\text{ nm}$ và $\lambda = 280\text{ nm}$. Độ tinh sạch của mẫu DNA được xác định qua tỷ số $\text{OD}{260}/\text{OD}{280}$ (chuẩn tinh sạch đạt xấp xỉ 1,8 - 2,0).
  2. Kỹ thuật phản ứng chuỗi Polymerase (PCR): Khuếch đại đoạn gen 16S rRNA bằng cặp mồi phổ quát 27F và 1492R với bộ hóa chất 2X PCR Master Mix (i-Taq).
  3. Điện di gel agarose: Kiểm tra kích thước và độ đặc hiệu của sản phẩm PCR trên gel agarose 1% trong đệm 1X TAE, nhuộm bằng Ethidium bromide, chạy ở điện áp 80V trong 45 phút và chụp ảnh bản gel dưới hệ thống ghi hình UVP BioDoc-it.
  4. Giải trình tự và định danh phân tử: Giải trình tự đoạn gen 16S rRNA và đối soát độ tương đồng trình tự nucleotide trên ngân hàng dữ liệu NCBI bằng công cụ BLASTN.
  5. Phân tích phát sinh loài (Phylogenetic analysis): Xây dựng cây phát sinh loài bằng phần mềm MEGA7 theo thuật toán Maximum Likelihood (ML) với giá trị lặp lại (Bootstrap) là 1000 lần.
  6. Dự đoán cấu trúc bậc hai: Sử dụng hai công cụ trực tuyến UnaFold và R2DT để dự đoán cấu trúc không gian bậc hai của phân tử RNA từ chuỗi nucleotide, xác định các điểm đa hình cấu trúc giữa các chủng vi khuẩn.
  7. Phân tích mức độ bảo thủ trình tự: Dóng hàng trình tự chuỗi (multiple sequence alignment) trên phần mềm MEGA7 với thông số toggle 50%, gap open -400 và gap extend 0 để xác định các vị trí nucleotide đột biến và các vùng bảo thủ đặc trưng của chi.
Quy trình Hóa chất và vật liệu chính
Đo quang phổ DNA Mẫu DNA tổng số, đệm TE buffer (tỷ lệ 1:10), cuvette Hellma traycell 1mm (factor 10)
Phản ứng PCR 2X PCR master mix solution kit (i-Taq), DNA khuôn mẫu (5 ng/μl), mồi xuôi 27F, mồi ngược 1492R, nước khử ion
Điện di bản gel Agarose 1%, đệm 1X TAE, Ethidium Bromide, 6X Loading dye, DNA ladder chuẩn 1kb
Tin sinh học NCBI BLAST, phần mềm MEGA7, máy chủ phân tích R2DT, công cụ UnaFold

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

Part I: Introduction (Mở đầu)

Chương mở đầu trình bày tầm quan trọng của cây khoai tây trong nền nông nghiệp toàn cầu và tại Việt Nam. Tác giả nhấn mạnh xu hướng ứng dụng vi sinh vật nội sinh nhằm giảm thiểu phân bón hóa học và thuốc bảo vệ thực vật, hướng đến nông nghiệp bền vững. Nội dung nêu rõ lý do chọn đề tài, xác định mục tiêu định danh 12 chủng vi khuẩn nội sinh phân lập từ rễ khoai tây và ứng dụng tin sinh học để phân tích đa dạng cấu trúc bậc hai của đoạn gen 16S rRNA.

Part II: Literature Review (Tổng quan tài liệu)

Chương tổng quan tập trung hệ thống hóa 5 mảng kiến thức cốt lõi:

  • Đặc điểm thực vật học của khoai tây: Tên khoa học Solanum tuberosum L., hình thái thân thảo cao khoảng 60 cm, chu kỳ sống qua các giai đoạn: ngủ nghỉ, nảy mầm, phát triển thân lá, hình thành thân ngầm (tia củ xuất hiện 30-40 ngày sau trồng), phình củ và chín. Các yếu tố sinh thái ảnh hưởng gồm nhiệt độ thích hợp cho phát triển thân lá (20-25°C), nhiệt độ tạo củ (17-20°C) và tích lũy tinh bột (16-18°C); cường độ ánh sáng 40.000 - 60.000 lux; nhu cầu nước 500-700 mm; độ pH đất thích hợp từ 5,0 - 7,0 (tối ưu 6,0 - 6,5).
  • Tổng quan về vi khuẩn nội sinh: Lịch sử nghiên cứu từ năm 1926, sự phân bố trong mô thực vật (rễ, thân, lá, củ, hạt). Nguồn gốc xâm nhiễm chủ yếu từ vùng rễ (rhizosphere) qua các vết thương cơ học, chóp rễ, lông hút, khí khổng và gian bào.
  • Vai trò sinh học của vi khuẩn nội sinh: Tăng cường tổng hợp phytohormone auxin (IAA), cố định nitơ sinh học, hòa tan photphat khó tan, tiết siderophore thu hút ion sắt, cảm ứng tính kháng bệnh cho cây chủ đối với các loài vi nấm gây hại và giảm thiểu tác hại của biến đổi khí hậu.
  • Các phương pháp phân tử định danh vi sinh vật: Phân tích ưu nhược điểm của các vùng chỉ thị ITS, 16S rRNA, 18S rRNA và 26S rRNA (miền D1/D2).
  • Các nghiên cứu liên quan: Tổng hợp các công trình quốc tế sử dụng chỉ thị 16S rRNA nghiên cứu quần xã vi khuẩn nội sinh trên khoai tây, bao gồm kỹ thuật PCR-DGGE, T-RFLP và nuôi cấy truyền thống, ghi nhận sự ưu thế của các nhóm Proteobacteria, Firmicutes, Actinomycetes và Bacteroidetes.

Part III: Materials and Methods (Vật liệu và phương pháp nghiên cứu)

Trình bày chi tiết danh mục thiết bị phòng thí nghiệm, hóa chất và quy trình thực nghiệm chuẩn hóa:

  • Thành phần phản ứng PCR: Tổng thể tích phản ứng là 9 μl, bao gồm: 2X PCR master mix (i-Taq): 5 μl; mẫu DNA (5 ng/μl): 2 μl; mồi 27F (10 μM): 1 μl; mồi 1492R (10 μM): 1 μl; nước cất vô trùng ($H_2O$): 1 μl.
  • Chu trình nhiệt PCR:
    1. Giai đoạn biến tính ban đầu: 95°C trong 5 phút (1 chu kỳ).
    2. Giai đoạn chu kỳ (35 chu kỳ): Biến tính ở 95°C trong 30 giây; bắt mồi ở 50°C trong 30 giây; kéo dài ở 72°C trong 2 phút.
    3. Giai đoạn kéo dài cuối cùng: 72°C trong 7 phút (1 chu kỳ).
    4. Bảo quản mẫu: 20°C vô hạn.
  • Quy trình điện di và phân tích dữ liệu: Đổ gel agarose 1% pha trong đệm 1X TAE bổ sung 5 μl Ethidium Bromide, nạp 2 μl loading dye cùng mẫu DNA, chạy điện di ở 80V trong 45 phút. Quy trình phân tích dữ liệu BLAST, xuất cây phân loại học MEGA7 và chạy mô hình nếp gấp RNA trên máy chủ UnaFold và R2DT.

Part IV: Results and Discussion (Kết quả và thảo luận)

Chương kết quả và thảo luận trình bày dữ liệu định danh 12 chủng vi khuẩn, phân tích cây phát sinh loài, dự đoán cấu trúc bậc hai và phân tích các vùng bảo thủ/biến đổi trên chuỗi nucleotide:

  • Sản phẩm PCR đoạn gen 16S rRNA với cặp mồi 27F/1492R cho vạch sáng rõ, kích thước đặc hiệu theo lý thuyết (~1400 - 1500 bp).
  • Kết quả tra cứu BLAST trên NCBI xác nhận cả 12 chủng vi khuẩn nội sinh phân lập từ rễ khoai tây đều thuộc chi Bacillus.
Ký hiệu mẫu Chủng tương đồng nhất (NCBI) Tên loài khoa học Max Score Total Score Query Cover E-value Per. Ident
isolate_1 Bacillus thuringiensis strain ERP1 Bacillus thuringiensis 756 756 99% 0.0 94.1%
isolate_2 Bacillus subtilis strain CUJ-41 Bacillus subtilis 887 887 100% 0.0 100.0%
isolate_3 Bacillus paralicheniformis strain HR-1 Bacillus paralicheniformis 889 889 100% 0.0 100.0%
isolate_4 Bacillus subtilis strain CUJ-41 Bacillus subtilis 887 887 100% 0.0 100.0%
isolate_5 Bacillus subtilis strain CUJ-41 Bacillus subtilis 887 887 100% 0.0 100.0%
isolate_6 Bacillus paralicheniformis strain AJVR1 Bacillus paralicheniformis 881 881 99% 0.0 99.1%
isolate_7 Bacillus subtilis strain CUJ-41 Bacillus subtilis 887 887 100% 0.0 100.0%
isolate_8 Bacillus sp. strain RWP_D7 Bacillus sp. - - - 0.0 -
isolate_9 Bacillus paralicheniformis strain HR-1 Bacillus paralicheniformis 883 883 100% 0.0 99.1%
isolate_10 Bacillus paralicheniformis strain HR-1 Bacillus paralicheniformis 889 889 100% 0.0 100.0%
isolate_11 Bacillus subtilis strain CUJ-41 Bacillus subtilis 887 887 100% 0.0 100.0%
isolate_12 Bacillus subtilis strain CUJ-41 Bacillus subtilis 887 887 100% 0.0 100.0%
  • Phân tích cây quan hệ phát sinh loài (Phylogenetic Tree): Cây Maximum Likelihood chia 12 mẫu thành 3 nhánh chính rõ rệt:
    • Nhánh 1 (Branch 1): Gồm 7 chủng (isolate_2, isolate_4, isolate_5, isolate_7, isolate_8, isolate_11, isolate_12), có mức độ tương đồng cao với Bacillus subtilis chủng CUJ-41.
    • Nhánh 2 (Branch 2): Gồm 4 chủng (isolate_3, isolate_6, isolate_9, isolate_10), thuộc nhóm Bacillus paralicheniformis. Chiều dài nhánh khác nhau giữa các chủng (dù cùng tương đồng với chủng HR-1) phản ánh sự phân kỳ và tích lũy các đột biến điểm trong chuỗi nucleotide.
    • Nhánh 3 (Branch 3): Gồm 1 chủng duy nhất là isolate_1 thuộc loài Bacillus thuringiensis.
  • Phân tích cấu trúc bậc hai (Secondary Structure):
    • Mô hình dự đoán qua R2DT cho thấy độ tương đồng hình thái cao giữa các chủng trong cùng một phân nhóm.
    • Mô hình dự đoán qua UnaFold có độ nhạy cao hơn, phát hiện rõ các điểm sai khác cấu trúc nếp gấp RNA tại các vị trí nucleotide 160, 200 và 240 ở Nhóm II (B. paralicheniformis), chứng minh sự đa dạng đa hình hình thái phân tử ngay cả khi trình tự sơ cấp tương đồng cao.
  • Tính bảo thủ và biến đổi của trình tự:
    • Vùng nucleotide 1 - 60: Duy nhất Bacillus thuringiensis (isolate_1) có sự sai khác nucleotide so với 11 chủng còn lại.
    • Vùng nucleotide 61 - 120: Xuất hiện các vị trí nucleotide phân tách đặc trưng giữa nhóm Bacillus paralicheniformis và nhóm Bacillus subtilis.

Part V: Conclusions and Recommendations (Kết luận và kiến nghị)

Tổng kết toàn bộ các phát hiện chính của khóa luận về việc xác định thành phần loài vi khuẩn nội sinh và đưa ra kiến nghị học thuật phục vụ các nghiên cứu tiếp nối.


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

Kết quả chính

  • Khuếch đại thành công đoạn gen 16S rRNA từ 12 mẫu vi khuẩn nội sinh rễ khoai tây bằng phản ứng PCR với cặp mồi 27F/1492R.
  • Xác định 100% các mẫu phân lập (12/12 chủng) thuộc về chi Bacillus, phân bố trong 3 nhóm loài chính:
    • 7 chủng thuộc loài Bacillus subtilis (chiếm 58,33%).
    • 4 chủng thuộc loài Bacillus paralicheniformis (chiếm 33,33%).
    • 1 chủng thuộc loài Bacillus thuringiensis (chiếm 8,34%).
  • Thiết lập cây phát sinh loài Maximum Likelihood (Bootstrap 1000) phân tách 12 chủng thành 3 nhánh phân loại độc lập.
  • Chứng minh công cụ UnaFold có khả năng nhận diện các biến dị cấu trúc bậc hai nhạy hơn R2DT, chỉ ra các biến đổi không gian tại các vị trí nucleotide 160, 200 và 240 của chuỗi RNA.
  • Định vị được vùng bảo thủ chung của chi Bacillus và các vùng biến đổi đặc trưng nhận diện loài ở các đoạn nucleotide 1 - 60 và 61 - 120.

Đóng góp của đề tài

  • Bổ sung cơ sở dữ liệu trình tự nucleotide và phân loại chính xác các chủng vi khuẩn nội sinh có nguồn gốc từ rễ khoai tây tại Việt Nam, làm tiền đề cho việc sàng lọc các chủng vi khuẩn có lợi ứng dụng sản xuất chế phẩm sinh học đối kháng bệnh cây và kích thích sinh trưởng.
  • Kết hợp phương pháp dự đoán cấu trúc bậc hai (UnaFold/R2DT) vào đánh giá đa dạng di truyền vi sinh vật bên cạnh việc so sánh chuỗi sơ cấp 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

  • Phạm vi mẫu nghiên cứu dừng lại ở 12 chủng vi khuẩn nội sinh đã được phân lập trước đó, chưa bao quát toàn bộ hệ vi sinh vật nội sinh phức tạp trên các mô khác (như thân, lá, củ) hoặc các giống khoai tây khác nhau.
  • Chỉ giải trình tự một phần đoạn gen 16S rRNA (partial 16S rRNA sequence) nên khả năng phân giải dưới mức loài (phân loại dưới loài/chủng chuyên biệt) còn giới hạn.

Hướng nghiên cứu tiếp theo (Kiến nghị từ tác giả)

  • Ứng dụng bổ sung các hệ thống chỉ thị phân tử đa hình khác (molecular markers) để đánh giá sâu hơn sự phân kỳ di truyền giữa các chủng vi khuẩn nội sinh cùng loài.
  • Mở rộng nghiên cứu đánh giá hoạt tính sinh học in vitro và in vivo của các chủng Bacillus subtilis, Bacillus paralicheniformis và Bacillus thuringiensis đã định danh (khả năng sinh IAA, hòa tan lân, sinh siderophore và đối kháng nấm Phytophthora infestans).

Giá trị tham khảo

  • Đối tượng tham khảo: Sinh viên, học viên cao học và cán bộ nghiên cứu thuộc các chuyên ngành Công nghệ sinh học Thực vật, Vi sinh học, Nông nghiệp sinh thái và Bảo vệ thực vật.
  • Nội dung giá trị:
    • Quy trình thực nghiệm chuẩn từ định lượng DNA, điều kiện chu trình nhiệt PCR đoạn gen 16S rRNA với mồi 27F/1492R, đến kỹ thuật điện di kiểm tra sản phẩm.
    • Phương pháp xử lý dữ liệu trình tự tin sinh học: BLAST đối soát cơ sở dữ liệu NCBI, thiết lập cây phát sinh loài Maximum Likelihood trên MEGA7 và phân tích cấu trúc không gian bậc hai RNA bằng R2DT/UnaFold.

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

1. 12 chủng vi khuẩn nội sinh rễ khoai tây trong nghiên cứu được định danh thuộc những loài nào?

Cả 12 chủng đều thuộc chi Bacillus, bao gồm 3 loài:

  • Bacillus subtilis (7 chủng: isolate_2, isolate_4, isolate_5, isolate_7, isolate_8, isolate_11, isolate_12).
  • Bacillus paralicheniformis (4 chủng: isolate_3, isolate_6, isolate_9, isolate_10).
  • Bacillus thuringiensis (1 chủng: isolate_1).

2. Cặp mồi và điều kiện nhiệt độ của phản ứng PCR khuếch đại gen 16S rRNA được thiết lập như thế nào?

  • Cặp mồi: 27F và 1492R (nồng độ 10 μM, mỗi mồi 1 μl trong phản ứng 9 μl).
  • Chu trình nhiệt: Biến tính ban đầu 95°C trong 5 phút; tiếp theo là 35 chu kỳ gồm [95°C trong 30 giây, 50°C trong 30 giây, 72°C trong 2 phút]; kéo dài cuối cùng ở 72°C trong 7 phút và giữ mẫu ở 20°C.

3. Cây phát sinh loài của 12 chủng vi khuẩn được xây dựng bằng phần mềm và thuật toán nào?

Cây phát sinh loài được xây dựng bằng phần mềm MEGA7, sử dụng thuật toán Maximum Likelihood (ML) với chỉ số kiểm tra lặp lại Bootstrap là 1000 lần, phân tách 12 chủng thành 3 nhánh phát sinh tương ứng với 3 loài vi khuẩn.

4. Sự khác biệt giữa hai công cụ dự đoán cấu trúc bậc hai R2DT và UnaFold trong nghiên cứu là gì?

Công cụ R2DT cho ra mô hình cấu trúc bậc hai có mức độ tương đồng rất cao giữa các chủng trong cùng một nhóm loài. Trong khi đó, công cụ UnaFold có độ phân giải nhạy hơn, chỉ ra các biến đổi đa hình nếp gấp RNA rõ rệt giữa các chủng, đặc biệt tại các vị trí nucleotide 160, 200 và 240 của nhóm Bacillus paralicheniformis.


Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Thị Ngọc Anh đã định danh thành công 12 chủng vi khuẩn nội sinh phân lập từ rễ khoai tây bằng kỹ thuật khuếch đại và giải trình tự một phần gen 16S rRNA. Toàn bộ 12 chủng đều thuộc chi Bacillus, phân bố vào 3 loài chính gồm Bacillus subtilis, Bacillus paralicheniformis và Bacillus thuringiensis. Việc kết hợp xây dựng cây phát sinh loài trên MEGA7 và mô hình hóa cấu trúc bậc hai qua công cụ UnaFold/R2DT đã làm sáng tỏ tính đa dạng phân tử và các vùng biến đổi nucleotide đặc trưng của hệ vi khuẩn nội sinh rễ khoai tây, cung cấp dữ liệu học thuật cho các nghiên cứu ứng dụng vi sinh vật có ích trong nông nghiệp sinh học.