Tổng quan nghiên cứu

Lúa gạo (Oryza sativa L.) là nguồn lương thực thiết yếu nuôi sống gần 50% dân số toàn cầu, giữ vị trí then chốt trong việc bảo đảm an ninh lương thực tại khu vực châu Á. Tuy nhiên, sản lượng lúa luôn đối mặt với nhiều hiểm họa dịch bệnh nguy hiểm, nổi bật là bệnh cháy bìa lá (bạc lá) do vi khuẩn Xanthomonas oryzae pv. oryzae (Xoo) gây ra. Khi dịch bệnh bùng phát mạnh, thiệt hại về năng suất vụ mùa có thể lên tới 75%, làm suy giảm nghiêm trọng thu nhập của người nông dân. Tại Việt Nam, thống kê nông nghiệp ghi nhận diện tích lúa nhiễm bệnh trên cả nước năm 2013 đạt khoảng 135,4 nghìn ha, và riêng vùng trọng điểm Đồng bằng sông Cửu Long năm 2019 đã có tới 50 đến 60 nghìn ha lúa bị xâm nhiễm nặng nề.

Trước thực trạng các biện pháp hóa học truyền thống như thuốc trừ sâu gốc đồng hay chất kháng sinh bộc lộ nhiều hạn chế về ô nhiễm môi trường, nguy cơ tồn dư độc chất và kích thích vi khuẩn phát triển tính kháng thuốc, việc ứng dụng tác nhân sinh học trở thành xu thế tất yếu. Nghiên cứu của học viên Phan Đình Hải thuộc chuyên ngành Công nghệ sinh học tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh, thực hiện từ tháng 02/2022 đến tháng 12/2022, đã tập trung đánh giá hoạt tính xâm nhiễm và giải mã bộ gene của 03 chủng thực khuẩn thể bản địa M223, M624, M625. Công trình mang ý nghĩa khoa học và thực tiễn cao khi bổ sung nguồn dữ liệu genome hoàn chỉnh của thực khuẩn thể tại Việt Nam, chứng minh khả năng ức chế vi khuẩn Xoo đạt 100% trong điều kiện in vitro suốt 48 giờ, mở ra hướng đi triển vọng cho sản xuất chế phẩm sinh học an toàn.

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 được xây dựng dựa trên nền tảng lý thuyết tương tác giữa thực khuẩn thể (bacteriophage) và tế bào vi khuẩn ký chủ, cụ thể là chu trình tiêu tan (lytic cycle) và chu trình tiềm tan (lysogenic cycle). Để đánh giá khả năng tiêu diệt mầm bệnh, mô hình động học tăng trưởng một bước (one-step growth kinetics) của Ellis và Delbrück được áp dụng nhằm xác định các thông số định lượng chuẩn xác như thời gian tiềm ẩn (latent period) và hệ số nhân phóng thích (burst size). Bên cạnh đó, nghiên cứu vận dụng lý thuyết giải trình tự thế hệ mới (Next Generation Sequencing - NGS) trên nền tảng công nghệ tổng hợp mạch gắn huỳnh quang thuận nghịch của Illumina, kết hợp các thuật toán tin sinh học hiện đại để lắp ráp de novo và chú giải chức năng hệ gene.

Các khái niệm trọng tâm xuyên suốt công trình bao gồm:

  • Thực khuẩn thể (Bacteriophage): Tác nhân virus chuyên biệt ký sinh và tiêu diệt vi khuẩn mà không gây hại cho sinh vật bậc cao và hệ vi sinh có lợi.
  • Họ Siphoviridae: Nhóm thực khuẩn thể có cấu trúc vỏ capsid đối xứng hình khối và phần đuôi dài, không co rút, chứa vật chất di truyền DNA chuỗi đôi.
  • Tỷ lệ lây nhiễm (Multiplicity of Infection - MOI): Tỷ số định lượng giữa số lượng hạt thực khuẩn thể và số lượng tế bào vi khuẩn đích trong môi trường phản ứng.
  • Khung đọc mở (Open Reading Frame - ORF): Đoạn trình tự nucleotide liên tục có tiềm năng mã hóa chuỗi polypeptide hoàn chỉnh của virus.

Phương pháp nghiên cứu

Nghiên cứu sử dụng mẫu vi khuẩn chỉ thị Xanthomonas oryzae pv. oryzae LA1+ cung cấp bởi Khoa Nông nghiệp, Trường Đại học Cần Thơ, nuôi cấy đạt mật độ chuẩn 2x10^7 CFU/mL (tương đương giá trị đo OD600 là 0,1). Cỡ mẫu thực nghiệm bao gồm 03 chủng thực khuẩn thể M223, M624 và M625 được phân lập trực tiếp từ các mẫu lá lúa nhiễm bệnh thu thập tại các cánh đồng thuộc Đồng bằng sông Cửu Long.

Quy trình khảo sát hoạt tính xâm nhiễm được thực hiện qua phương pháp tạo vòng sinh tan (plaque assay) trên đĩa thạch hai lớp. Động học chu kỳ lây nhiễm được thiết lập ở tỷ lệ MOI là 0,01, tiến hành lấy mẫu định kỳ 5 phút một lần từ thời điểm 0 đến 120 phút với 3 lần lặp lại độc lập (n = 3). Khả năng ức chế vi khuẩn trong môi trường lỏng TSB được kiểm tra ở tỷ lệ MOI là 2,0, theo dõi mật độ quang OD600 liên tục mỗi giờ trong 48 giờ.

Về phương pháp phân tích hệ gene, DNA tổng số của thực khuẩn thể được giải trình tự bằng hệ thống Illumina Paired-End với chiều dài đoạn đọc từ 35 đến 151 bp, bảo đảm điểm chất lượng trình tự vượt ngưỡng Q30 (độ chính xác đạt 99,99%). Dữ liệu thô được kiểm tra chất lượng qua FastQC v0.11 và lọc bỏ adapter bằng Trimmomatic v0.39. Quá trình lắp ráp de novo thực hiện bằng Unicycler v0.4.8 tích hợp thuật toán đồ thị de Bruijn. Chú giải cấu trúc và chức năng gene được thực hiện trên máy chủ RAST v2.0 kết hợp tra cứu dữ liệu tương đồng qua BLASTP, cơ sở dữ liệu miền bảo tồn NCBI CDD. Đánh giá tính an toàn sinh học sử dụng hệ thống RGI phân tích cơ sở dữ liệu gen kháng kháng sinh CARD và các nhân tố độc lực PATRIC. Cây phát sinh loài được xây dựng bằng phần mềm IQ-TREE theo phương pháp Maximum Likelihood với giá trị kiểm định 1000 bootstrap dựa trên 02 protein bảo tồn là Terminase Large Subunit (TerL) và Major Capsid Protein (MCP).

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

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

Quan sát dưới kính hiển vi điện tử truyền qua (TEM) kết hợp đo đạc hình thái vòng sinh tan kích thước 1 đến 2 mm cho thấy cả 3 chủng thực khuẩn thể M223, M624, M625 đều thuộc họ Siphoviridae. Về động học xâm nhiễm, các chủng thể hiện chu kỳ lây nhiễm hiệu quả cao:

  • Chủng M223 có thời gian tiềm ẩn từ 25 đến 30 phút, hệ số nhân đạt 13 ± 2 hạt phage trên mỗi tế bào bị xâm nhiễm.
  • Chủng M624 có thời gian tiềm ẩn từ 30 đến 35 phút, hệ số nhân đạt 16 ± 3 hạt phage trên mỗi tế bào.
  • Chủng M625 thể hiện ưu thế vượt trội với thời gian tiềm ẩn từ 30 đến 35 phút, hệ số nhân bùng phát đạt 28 ± 5 hạt phage trên mỗi tế bào, chu kỳ xâm nhiễm hoàn tất trong 100 phút (năng suất nhân bản của M625 cao hơn M223 tới 115,4%).

Trong thử nghiệm ức chế vi khuẩn in vitro, cả 3 chủng phage ở tỷ lệ MOI là 2,0 đều cho thấy hiệu năng kiểm soát ấn tượng. Sau 3 đến 4 giờ bổ sung phage, tốc độ tăng sinh của vi khuẩn Xoo LA1+ bị chặn đứng hoàn toàn. Môi trường nuôi cấy chuyển từ đục sang trong suốt sau 7 đến 8 giờ và duy trì trạng thái ổn định, vi khuẩn hoàn toàn không có hiện tượng huyền phù hay tái sinh trở lại trong suốt 48 giờ theo dõi.

Phân tích bộ gene hoàn chỉnh xác định kích thước hệ gene của cả 3 chủng M223, M624 và M625 đều đạt 46.226 bp dạng DNA chuỗi đôi mạch vòng, tỷ lệ GC đạt khoảng 60% đến 61%, tỷ lệ nucleotide khuyết thiếu (N) là 0%. Số lượng khung đọc mở (ORF) dự đoán lần lượt là 74 ORF (M223), 76 ORF (M624) và 75 ORF (M625). Kết quả so sánh BLASTN cho thấy hệ gene tương đồng cao tới 99,13% với chủng Xanthomonas phage XPV2 đã công bố quốc tế. Đặc biệt, phân tích an toàn sinh học khẳng định 100% các chủng phage không mang bất kỳ gene kháng kháng sinh, nhân tố độc lực, gene sinh độc tố hay gene tích hợp integrase nào.

Thảo luận kết quả

Khả năng làm trong môi trường nuôi cấy vi khuẩn nhanh chóng sau 3 đến 4 giờ chứng minh tính chất tiêu tan mạnh mẽ của các dòng thực khuẩn thể bản địa. Khi biểu diễn diễn biến mật độ quang OD600 theo thời gian trên đồ thị đường thẳng, có thể thấy rõ sự phân kỳ triệt để: trong khi mẫu đối chứng không xử lý có mật độ vi khuẩn tăng vọt theo hàm mũ, các mẫu bổ sung phage M223, M624, M625 đều có đường biểu diễn giảm mạnh về sát trục hoành sau 8 giờ. Cơ chế này xuất phát từ hoạt tính của các module gene ly giải (lysin và holin) được mã hóa trong hệ gene, nhanh chóng phá vỡ màng tế bào vi khuẩn sau khi hoàn thành chu kỳ nhân bản.

Việc không phát hiện gene integrase hay protein ức chế trong bộ gene 46.226 bp khẳng định cả 3 thực khuẩn thể đều là phage độc lực tuyệt đối (obligate lytic phages), hoàn toàn không thể thực hiện chu trình tiềm tan để dung hợp vật chất di truyền vào hệ gene của vi khuẩn. Cây phát sinh loài xây dựng trên 1000 lần lặp bootstrap từ chuỗi polypeptide TerL và MCP xác định mối quan hệ di truyền chặt chẽ giữa các chủng phân lập tại Việt Nam với nhóm phage kiểm soát sinh học trên thế giới. Đây là bằng chứng khoa học vững chắc loại bỏ nguy cơ chuyển gene kháng thuốc hoặc độc lực sang các chủng vi sinh vật đất có lợi, bảo đảm an toàn sinh thái khi ứng dụng ngoài tự nhiên.

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

Nhằm chuyển hóa các kết quả nghiên cứu trong phòng thí nghiệm thành giải pháp thực tế bảo vệ cây lúa, các khuyến nghị can thiệp cụ thể được đề xuất như sau:

  • Phát triển công thức tổ hợp đa chủng (Phage Cocktail): Phối trộn 03 chủng M223, M624 và M625 theo tỷ lệ tối ưu 1:1:1 trong vòng 6 đến 12 tháng tới do Bộ môn Công nghệ sinh học chủ trì, nhằm mở rộng phổ lây nhiễm trên nhiều chủng Xoo đột biến và giảm thiểu nguy cơ hình thành tính kháng của vi khuẩn xuống dưới 10%.
  • Mở rộng khảo nghiệm trên đồng ruộng: Tiến hành thử nghiệm nhà lưới và mô hình cánh đồng mẫu lớn tại các vùng có áp lực bệnh cao như An Giang, Đồng Tháp trong lộ trình 2 đến 3 vụ mùa (2024 - 2026), hướng tới mục tiêu giảm chiều dài vết bệnh cháy bìa lá từ 65% đến 75% so với nhóm đối chứng không xử lý.
  • Nghiên cứu chất mang và phụ gia bảo vệ sinh học: Phối hợp cùng các doanh nghiệp nông nghiệp trong 18 tháng để nghiên cứu vỏ bọc phân tử hoặc phụ gia chống tia cực tím (UV-protectant), giúp duy trì hoạt tính xâm nhiễm của thực khuẩn thể trên bề mặt lá lúa kéo dài từ 4 giờ lên trên 72 giờ dưới ánh nắng nhiệt đới.
  • Xây dựng tiêu chuẩn kỹ thuật sản xuất sinh khối: Thiết lập quy trình chuẩn nuôi cấy lên men quy mô công nghiệp đạt mật độ tối thiểu 1x10^9 PFU/mL trong thời gian 24 tháng, hoàn thiện hồ sơ đăng ký chế phẩm sinh học bảo vệ thực vật theo quy chuẩn của Cục Trồng trọt và Bảo vệ thực vật.

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

Công trình nghiên cứu cung cấp khối lượng dữ liệu thực nghiệm và phương pháp luận chuyên sâu, mang lại giá trị thiết thực cho 4 nhóm đối tượng chính:

  • Nhà nghiên cứu và học viên cao học ngành Công nghệ sinh học, Bảo vệ thực vật: Tiếp cận quy trình chuẩn từ phân lập, đánh giá động học sinh tan đến các pipeline phân tích tin sinh học xử lý dữ liệu NGS (FastQC, Unicycler, RASTtk, IQ-TREE) trên đối tượng thực khuẩn thể.
  • Kỹ sư R&D tại các doanh nghiệp sản xuất chế phẩm sinh học: Sử dụng các thông số về hệ số nhân 28 ± 5 PFU/tế bào và dữ liệu an toàn bộ gene 46.226 bp làm cơ sở thiết kế quy trình sản xuất thuốc trừ bệnh sinh học nguồn gốc virus.
  • Cán bộ Chi cục Trồng trọt và Bảo vệ thực vật tại các địa phương: Tham khảo luận cứ khoa học để xây dựng mô hình quản lý dịch hại tổng hợp (IPM), từng bước thay thế thuốc bảo vệ thực vật hóa học và giảm thiểu 30% đến 50% chi phí xử lý bệnh bạc lá lúa cho nông dân.
  • Cơ quan quản lý an toàn sinh học và môi trường: Khai thác quy trình kiểm tra các nhân tố độc lực và gene kháng kháng sinh từ dữ liệu CARD làm tài liệu tham chiếu xây dựng khung đánh giá rủi ro sinh học đối với các chế phẩm vi sinh nông nghiệp.

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

Thực khuẩn thể M625 có ưu điểm gì vượt trội hơn so với các chủng còn lại?
Chủng M625 sở hữu hệ số nhân phóng thích cao nhất đạt 28 ± 5 hạt phage trên mỗi tế bào bị xâm nhiễm, cao hơn đáng kể so với M223 (13 ± 2) và M624 (16 ± 3). Hệ số nhân lớn kết hợp cùng chu kỳ lây nhiễm ngắn 100 phút giúp M625 nhanh chóng gia tăng số lượng trong môi trường, tối ưu hóa tốc độ tiêu diệt vi khuẩn Xoo.

Vì sao phân tích bộ gene hoàn chỉnh là yêu cầu bắt buộc trước khi ứng dụng thực khuẩn thể?
Phân tích bộ gene nhằm loại trừ tuyệt đối sự hiện diện của gene tích hợp (integrase), gene kháng kháng sinh và các độc tố vi khuẩn. Với kích thước 46.226 bp được chú giải sạch hoàn toàn các yếu tố nguy hại, 03 chủng phage này được chứng minh chỉ thực hiện chu trình tiêu tan, bảo đảm an toàn sinh học cao nhất cho hệ sinh thái nông nghiệp.

Vi khuẩn gây bệnh cháy bìa lá có thể hình thành tính kháng lại thực khuẩn thể không?
Vi khuẩn có khả năng đột biến thụ thể bề mặt để kháng lại một chủng phage đơn lẻ. Tuy nhiên, việc áp dụng công thức phối hợp nhiều chủng như M223, M624 và M625 sẽ ngăn chặn nguy cơ này, do mỗi chủng có vị trí tiếp cận đa dạng, làm giảm xác suất vi khuẩn kháng thuốc xuống dưới mức 5% trong điều kiện đồng ruộng.

Sử dụng thực khuẩn thể có ảnh hưởng đến các vi khuẩn có ích khác trong đất và nước không?
Thực khuẩn thể có tính đặc hiệu vật chủ cực kỳ cao, chỉ nhận diện và lây nhiễm các thụ thể màng chuyên biệt của vi khuẩn Xanthomonas oryzae pv. oryzae. Các vi sinh vật có lợi trong môi trường đất như Bacillus hay vi khuẩn cố định đạm hoàn toàn không bị ảnh hưởng, giúp bảo tồn cân bằng sinh thái đất lúa.

Chế phẩm thực khuẩn thể có thời gian tác động kiểm soát vi khuẩn bao lâu trong phòng thí nghiệm?
Trong điều kiện in vitro ở nồng độ MOI là 2,0, thực khuẩn thể bắt đầu ức chế rõ rệt vi khuẩn sau 3 đến 4 giờ và làm dung dịch trong suốt hoàn toàn sau 7 đến 8 giờ. Mật độ tế bào vi khuẩn bị duy trì ở mức cực thấp và không ghi nhận hiện tượng bùng phát trở lại trong suốt 48 giờ theo dõi.

Kết luận

  • Luận văn đã phân lập và xác định thành công đặc tính hình thái của 03 chủng thực khuẩn thể bản địa M223, M624, M625 thuộc họ Siphoviridae có khả năng lây nhiễm đặc hiệu vi khuẩn Xoo LA1+.
  • Động học xâm nhiễm ghi nhận thời gian tiềm ẩn ngắn (25 đến 35 phút), trong đó chủng M625 đạt hệ số nhân bùng phát cao nhất là 28 ± 5 PFU/tế bào.
  • Thử nghiệm in vitro khẳng định hiệu quả kiểm soát tuyệt đối khi ức chế hoàn toàn vi khuẩn sau 3 đến 4 giờ và duy trì độ trong suốt của môi trường liên tục trong 48 giờ.
  • Hoàn thành giải mã và lắp ráp bộ gene với độ chính xác cao, kích thước 46.226 bp, chứa từ 74 đến 76 khung đọc mở và đạt độ tương đồng 99,13% với Xanthomonas phage XPV2.
  • Kiểm định an toàn sinh học xác nhận hệ gene không chứa bất kỳ gene kháng kháng sinh, nhân tố độc lực hay gene tích hợp integrase nào.

Đóng góp lớn nhất của công trình là cung cấp nguồn vật liệu sinh học giá trị cùng bộ dữ liệu bộ gene hoàn chỉnh của thực khuẩn thể kháng bệnh cháy bìa lá lúa tại Đồng bằng sông Cửu Long. Trong lộ trình 12 đến 24 tháng tới, các thử nghiệm thực địa và công nghệ tạo chế phẩm thương mại cần được đẩy mạnh triển khai. Các đơn vị nghiên cứu, trường đại học và doanh nghiệp nông nghiệp công nghệ cao hãy tăng cường hợp tác chuyển giao công nghệ để sớm đưa giải pháp sinh học tiên tiến này vào phục vụ canh tác lúa bền vững tại Việt Nam.