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

Cây đu đủ (Carica papaya L., thuộc họ Caricaceae) là cây ăn quả nhiệt đới có giá trị dinh dưỡng và hiệu quả kinh tế cao. Trong 100 g thịt quả chín chứa 86,6% nước, 12,1% tinh bột, 0,6% protein, 0,3% lipid, 0,7% xơ, cùng hàm lượng khoáng và vitamin đa dạng như kali (204 mg), canxi (34 mg), vitamin A (450 mg), vitamin C (74 mg). Tại Việt Nam, cây đu đủ được canh tác rộng rãi ở nhiều vùng sinh thái, đặc biệt là khu vực trung du và các tỉnh Đông Nam Bộ. Tuy nhiên, hoạt động sản xuất đu đủ đang đối mặt với thiệt hại nghiêm trọng do bệnh đốm vòng gây ra bởi Papaya ringspot virus (PRSV).

PRSV là loại virus lây lan nhanh qua côn trùng môi giới (rệp muội) và qua vết thương cơ học, không thể phòng trị bằng thuốc hóa học mà chủ yếu dựa vào việc sử dụng giống kháng và loại bỏ nguồn bệnh sớm. Việc chẩn đoán bệnh dựa trên triệu chứng hình thái thường gặp khó khăn ở giai đoạn đầu hoặc dễ nhầm lẫn do hiện tượng ẩn triệu chứng. Do đó, yêu cầu đặt ra là phải xây dựng được phương pháp xét nghiệm phân tử vừa có độ nhạy, độ đặc hiệu cao, vừa tiết kiệm chi phí và thời gian khi xử lý số lượng mẫu lớn.

Mục đích nghiên cứu: Xây dựng quy trình phát hiện và chẩn đoán virus PRSV gây bệnh đốm vòng trên cây đu đủ bằng phương pháp gộp mẫu (sample pooling) kết hợp kỹ thuật multiplex RT-PCR, nhằm nâng cao độ chính xác và hiệu quả kinh tế trong công tác quản lý dịch bệnh.

Nhiệm vụ nghiên cứu:

  1. Xây dựng quy trình phát hiện virus PRSV bằng kỹ thuật multiplex RT-PCR kết hợp đối chứng nội.
  2. Xây dựng quy trình sample pooling nhằm xác định ngưỡng phát hiện virus PRSV trong mẫu gộp.
  3. Ứng dụng quy trình sample pooling kết hợp multiplex RT-PCR để phát hiện PRSV trên các mẫu đu đủ thu thập từ thực địa.

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

  • Đối tượng nghiên cứu: Cây đu đủ (Carica papaya), Papaya ringspot virus (PRSV), bộ gen RNA tổng số của đu đủ và PRSV, các cặp mồi đặc hiệu khuếch đại gen P3 của PRSV và gen actin của đu đủ.
  • Phạm vi thời gian: Từ tháng 8/2022 đến tháng 12/2022.
  • Phạm vi không gian: Phòng 204 – Phòng Sinh học phân tử, Viện Nghiên cứu Công nghệ Sinh học và Môi trường, Trường Đại học Nông Lâm TP. Hồ Chí Minh; vùng thu mẫu thực địa tại xã Bình Hòa Nam, huyện Đức Huệ, tỉnh Long An (tọa độ 10°48'59.8"E).

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

Khung lý thuyết và đối tượng virus

Papaya ringspot virus (PRSV) thuộc chi Potyvirus, họ Potyviridae, có cấu trúc hạt xoắn mềm dẻo kích thước 760–800 nm × 12 nm, vật chất di truyền là phân tử RNA sợi đơn dương (ssRNA(+)) kích thước khoảng 12 kb. PRSV được chia thành hai chủng chính: PRSV-p (gây hại trên đu đủ và cây họ Bầu bí) và PRSV-w (chỉ gây hại trên cây họ Bầu bí). Virus lây truyền không bền vững qua các loài rệp muội như rệp bông (Aphis gossypii) và rệp đào (Myzus persicae), có sự hỗ trợ của protein thể vùi vô định hình (AI).

Phương pháp nghiên cứu thực nghiệm

  • Ly trích RNA tổng số: Sử dụng bộ kit EZ-10 Spin Column Plant RNA Mini-Preps Kit (BioBasic, Canada) dựa trên nguyên lý hấp phụ màng silica. Nồng độ và độ tinh sạch RNA được xác định bằng máy đo quang phổ BioDrop thông qua tỷ số $A_{260}/A_{280}$.
  • Phiên mã ngược (Reverse Transcription): Tổng hợp mạch cDNA từ khuôn RNA tổng số bằng bộ hóa chất SensiFAST cDNA Synthesis Kit (Bioline, Anh) với chu trình nhiệt: $25^\circ\text{C}$ (10 phút) $\rightarrow 42^\circ\text{C}$ (15 phút) $\rightarrow 48^\circ\text{C}$ (15 phút) $\rightarrow 85^\circ\text{C}$ (5 phút).
  • Thiết kế mồi và kiểm tra phân tử: Thiết kế mồi đặc hiệu bằng phần mềm Primer3 và kiểm tra độ đặc hiệu bằng công cụ BLAST trên NCBI. Đoạn mồi khuếch đại gen P3 của virus PRSV và gen actin của đu đủ (đối chứng nội nhằm kiểm soát chất lượng tách chiết và phản ứng).
  • Tối ưu hóa phản ứng Multiplex RT-PCR: Thực hiện khảo sát gradient nhiệt độ bắt cặp ($52–62^\circ\text{C}$) và tối ưu hóa 9 tổ hợp nồng độ mồi ($0{,}1,\mu\text{M}$; $0{,}2,\mu\text{M}$; $0{,}4,\mu\text{M}$).
  • Phương pháp gộp mẫu (Sample pooling): Khảo sát ngưỡng phát hiện của virus PRSV ở các tỷ lệ pha loãng RNA bệnh trong nước DEPC và trong RNA mô đu đủ khỏe (từ mức gộp 10 đến 100.000 lần) ở nồng độ chuẩn hóa $50 \pm 2\text{ ng/}\mu\text{L}$.
  • Nguồn dữ liệu và vật liệu thực nghiệm: Mẫu đối chứng dương PRSV từ Viện Nghiên cứu Công nghệ Sinh học và Môi trường; 30 mẫu lá đu đủ thực địa thu thập tại Đức Huệ, Long An; cơ sở dữ liệu trình tự đối chiếu mã số MT470189.1 trên NCBI.

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

Chương 1: Đặt vấn đề

Chương này tổng quan về giá trị kinh tế của cây đu đủ và tác hại nghiêm trọng của bệnh đốm vòng do PRSV gây ra đối với vùng trồng đu đủ ở Việt Nam và thế giới. Tác giả chỉ ra tính cấp thiết của việc phát triển một công cụ chẩn đoán phân tử có khả năng phát hiện nhanh, chính xác mầm bệnh ở giai đoạn sớm với chi phí thấp, từ đó xác lập mục tiêu, 3 nội dung nghiên cứu cụ thể và phạm vi thực hiện của khóa luận.

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

Chương 2 cung cấp hệ thống cơ sở lý thuyết bao gồm:

  • Đặc điểm thực vật học, sinh thái học, giá trị dinh dưỡng, dược liệu và tình hình sản xuất đu đủ tại Việt Nam và trên thế giới (Ấn Độ dẫn đầu sản lượng toàn cầu với 38,61% thị phần xuất khẩu giai đoạn 2008–2010; năng suất đu đủ tại Việt Nam đạt trung bình khoảng 20 tấn/ha).
  • Danh mục sâu bệnh hại chính trên cây đu đủ gồm côn trùng (rệp sáp, rệp muội, nhện đỏ, ruồi đục quả), nấm bệnh (Pythium aphanidermatum, Phyllosticta sulata, Helminthosporium rostratum, Oidium caricae), tuyến trùng (Meloidogyne incognita, Rotylenchulus reniformis) và virus (Papaya mosaic virus, Papaya ringspot virus).
  • Phân loại học, cấu trúc bộ gen, đặc tính lan truyền và triệu chứng của PRSV.
  • Tổng quan các phương pháp chẩn đoán PRSV: quan sát triệu chứng, cây chỉ thị, kính hiển vi điện tử, huyết thanh học (ELISA, DIBA) và sinh học phân tử (PCR, RT-PCR, Multiplex RT-PCR).
  • Nguyên lý và lịch sử ứng dụng phương pháp sample pooling trong tầm soát bệnh tốn ít mẫu đơn và tiết kiệm chi phí.
Tên cặp primer Trình tự mồi ($5' \rightarrow 3'$) Kích thước sản phẩm Gen mục tiêu / Nguồn
PRSV613-F
PRSV613-R
TTGTGTACGACTTCTCACCGAA
CGAATGTCATCCAAAGACTGATGATAAAC
613 bp Gen P3 của PRSV (Tuo và ctv, 2014)
PRSV331-F
PRSV331-R
TGTGTACGACTTCTCACCGA
CTTAACGCACAACTCAGCGTAT
331 bp Gen P3 của PRSV (Tự thiết kế)
ACT-F
ACT-R
TGGTCGTACAACTGGTATTG
ACAATGGATGGACCAGACTC
675 bp Gen actin đu đủ - Đối chứng nội (Tự thiết kế)

Chương 3: Vật liệu và phương pháp

Chương 3 trình bày chi tiết danh mục hóa chất, sinh phẩm, thiết bị phòng thí nghiệm cùng các quy trình kỹ thuật:

  • Quy trình 10 bước tách chiết RNA tổng số qua cột lọc silica.
  • Quy trình tổng hợp cDNA thể tích $10,\mu\text{L}$ và phản ứng PCR đơn mồi thể tích $25,\mu\text{L}$.
  • Ma trận 9 tổ hợp nồng độ mồi cho phản ứng multiplex RT-PCR thể tích $12{,}5,\mu\text{L}$ (gồm $6{,}25,\mu\text{L}$ MyTaq Mix 2X, mồi, $1,\mu\text{L}$ cDNA và nước không chứa nuclease).
  • Bố trí thí nghiệm sample pooling: thiết lập các mức gộp thể tích từ 10 đến $100.000,\mu\text{L}$ ($1,\mu\text{L}$ RNA bệnh kết hợp với RNA khỏe hoặc nước DEPC tương ứng).
  • Sơ đồ tầm soát 30 mẫu thực địa theo thuật toán phân cấp: kiểm tra mẫu gộp tổng 30 $\rightarrow$ phân 6 nhóm gộp 5 $\rightarrow$ kiểm tra từng cá thể trong nhóm dương tính.
Nồng độ mồi PRSV ($\mu\text{M}$) Nồng độ mồi Actin $0{,}1,\mu\text{M}$ Nồng độ mồi Actin $0{,}2,\mu\text{M}$ Nồng độ mồi Actin $0{,}4,\mu\text{M}$
$0{,}1$ $0{,}1,/,0{,}1$ $0{,}1,/,0{,}2$ $0{,}1,/,0{,}4$
$0{,}2$ $0{,}2,/,0{,}1$ $0{,}2,/,0{,}2$ $0{,}2,/,0{,}4$
$0{,}4$ $0{,}4,/,0{,}1$ $0{,}4,/,0{,}2$ $0{,}4,/,0{,}4$

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

Chương 4 phân tích toàn bộ số liệu thu được từ thực nghiệm:

  • Đánh giá chất lượng RNA ly trích: mẫu dương đạt nồng độ $110{,}2\text{ ng/}\mu\text{L}$ với tỷ số $A_{260}/A_{280} = 2{,}096$, đáp ứng tiêu chuẩn tinh sạch cho phản ứng khuếch đại.
  • Điện di sản phẩm RT-PCR mồi PRSV613 cho một băng duy nhất kích thước 613 bp. Kết quả giải trình tự và gióng hàng BLAST đạt độ tương đồng 95,4% với trình tự PRSV mã số MT470189.1 trên NCBI (công bố bởi Tran và Yeh, 2020).
  • Tối ưu hóa nhiệt độ bắt cặp: Khảo sát dải nhiệt độ $52–62^\circ\text{C}$ xác định được nhiệt độ ủ tối ưu là $56^\circ\text{C}$. Ở nhiệt độ này, cả hai phân đoạn gen 331 bp (PRSV) và 675 bp (actin) khuếch đại đồng đều và không xuất hiện băng phụ.
  • Tối ưu hóa nồng độ mồi: Nồng độ tối ưu cho phản ứng multiplex RT-PCR là $0{,}2,\mu\text{M}$ mồi PRSV331-F/R và $0{,}2,\mu\text{M}$ mồi ACT-F/R.
  • Đánh giá giới hạn phát hiện của phương pháp sample pooling: Khi pha loãng trong nước DEPC, sản phẩm khuếch đại xuất hiện ở mức pha loãng lên đến 10.000 lần. Khi pha loãng trong RNA đu đủ khỏe (nồng độ đo lại dao động từ $49{,}3$ đến $51{,}4\text{ ng/}\mu\text{L}$), ngưỡng phát hiện đạt mức gộp $1/1.000$ (1 mẫu bệnh trong 999 mẫu khỏe). Mức pha loãng $1/10.000$ và $1/100.000$ không cho tín hiệu khuếch đại do sự ức chế và cạnh tranh từ lượng lớn RNA nền khỏe mạnh.
  • Kết quả kiểm tra 30 mẫu thực địa tại Long An: Mẫu gộp 30 cho kết quả dương tính $\rightarrow$ kiểm tra 6 nhóm gộp 5 phát hiện nhóm 2 (mẫu 6–10) và nhóm 6 (mẫu 26–30) dương tính $\rightarrow$ kiểm tra từng mẫu riêng lẻ xác định chính xác mẫu số 9 và mẫu số 27 dương tính với PRSV.
Mức gộp mẫu ($\mu\text{L}$) Thể tích RNA bệnh ($\mu\text{L}$) Thể tích RNA khỏe ($\mu\text{L}$) Nồng độ RNA đo lại ($\text{ng/}\mu\text{L}$) Kết quả khuếch đại Multiplex RT-PCR
10 1 9 51,1 Dương tính (Băng sáng rõ)
20 1 19 50,3 Dương tính (Băng sáng rõ)
30 1 29 51,0 Dương tính (Băng sáng rõ)
40 1 39 50,3 Dương tính (Băng sáng rõ)
50 1 49 49,3 Dương tính (Băng sáng rõ)
100 1 99 51,2 Dương tính (Băng sáng rõ)
150 1 149 50,3 Dương tính (Băng sáng rõ)
200 1 199 50,1 Dương tính (Băng sáng rõ)
250 1 249 49,3 Dương tính (Băng sáng rõ)
300 1 299 51,4 Dương tính (Băng sáng rõ)
1.000 1 999 50,6 Dương tính (Băng mờ)
10.000 1 9.999 49,3 Âm tính (Không xuất hiện băng)
100.000 1 99.999 49,8 Âm tính (Không xuất hiện băng)

Chương 5: Kết luận và đề nghị

Chương cuối tóm tắt toàn bộ các thông số kỹ thuật đã tối ưu hóa, khẳng định tính khả thi của quy trình gộp mẫu $1/1.000$ kết hợp multiplex RT-PCR, đồng thời đề xuất mở rộng thử nghiệm trên quy mô mẫu thực địa lớn hơn.


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

Kết quả định lượng chính

  • Thiết lập thành công phản ứng Multiplex RT-PCR phát hiện đồng thời 2 đoạn gen: đoạn gen P3 của PRSV (kích thước 331 bp) và đoạn gen actin nội chứng của cây đu đủ (kích thước 675 bp).
  • Xác định chu trình nhiệt tối ưu: 1 chu kỳ $94^\circ\text{C}$ trong 4 phút; 35 chu kỳ gồm biến tính $94^\circ\text{C}$ (40 giây), bắt cặp $56^\circ\text{C}$ (30 giây), kéo dài $72^\circ\text{C}$ (1 phút); kết thúc bằng 1 chu kỳ kéo dài $72^\circ\text{C}$ trong 5 phút.
  • Xác định nồng độ mồi tối ưu: $0{,}2,\mu\text{M}$ đối với cặp mồi PRSV331-F/R và $0{,}2,\mu\text{M}$ đối với cặp mồi ACT-F/R.
  • Xác định ngưỡng phát hiện gộp mẫu: Phát hiện được 1 mẫu đu đủ nhiễm PRSV trong 999 mẫu đu đủ khỏe (tỷ lệ gộp $1/1.000$) ở nồng độ RNA chuẩn hóa $50 \pm 2\text{ ng/}\mu\text{L}$.
  • Ứng dụng thực địa: Tầm soát thành công 30 mẫu lá đu đủ tại Bình Hòa Nam, Đức Huệ, Long An, phát hiện chính xác 2 mẫu dương tính (mẫu 9 và mẫu 27). Quy trình gộp mẫu chỉ cần thực hiện 17 phản ứng PCR, tiết kiệm 43,3% số phản ứng so với phương pháp RT-PCR từng mẫu đơn lẻ (30 phản ứng) mà vẫn cho kết quả định tính hoàn toàn tương đương.

Đóng góp kỹ thuật của tác giả

  • Tích hợp thành công đối chứng nội gen actin đu đủ vào phản ứng multiplex RT-PCR, giúp kiểm soát chất lượng mẫu tách chiết và loại bỏ hoàn toàn các trường hợp âm tính giả do lỗi kỹ thuật.
  • Chuẩn hóa quy trình sample pooling 5 bước (Ly trích RNA $\rightarrow$ Chuẩn hóa nồng độ $50 \pm 2\text{ ng/}\mu\text{L} \rightarrow$ Gộp mẫu tỷ lệ $\le 1/1.000 \rightarrow$ Tổng hợp cDNA $\rightarrow$ Multiplex RT-PCR và điện di), cung cấp giải pháp sàng lọc diện rộng tiết kiệm thời gian, hóa chất và nhân lực phục vụ công tác kiểm dịch thực vật.

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

Hạn chế

  • Dung lượng mẫu khảo nghiệm thực địa còn khiêm tốn (30 mẫu thu thập tại một địa điểm duy nhất ở huyện Đức Huệ, tỉnh Long An).
  • Phản ứng multiplex RT-PCR ở mức gộp cao ($1/1.000$) cho băng điện di tương đối mờ do ảnh hưởng ức chế cạnh tranh từ RNA nền của cây khỏe.

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

  • Thử nghiệm và mở rộng quy trình sample pooling trên các tập mẫu thực địa lớn hơn tại nhiều vùng trồng đu đủ trọng điểm khác nhau.
  • Đưa quy trình vào ứng dụng thực tế trong công tác kiểm soát giống và kiểm dịch thực vật để đánh giá độ ổn định và độ nhạy lâm sàng.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo chuyên khảo hữu ích cho:

  • Sinh viên và học viên cao học ngành Công nghệ sinh học, Bảo vệ thực vật, Nông học: Tham khảo quy trình thiết kế mồi, kỹ thuật ly trích RNA thực vật bằng cột silica, phương pháp tối ưu hóa multiplex RT-PCR và phương pháp luận thiết lập mô hình sample pooling.
  • Cán bộ kỹ thuật phòng xét nghiệm phân tử và trạm bảo vệ thực vật: Ứng dụng trực tiếp bộ thông số chu trình nhiệt, nồng độ mồi và quy trình sàng lọc phân cấp mẫu gộp để phát hiện sớm PRSV, giúp giảm thiểu đáng kể chi phí xét nghiệm trong các chương trình giám sát dịch tễ trên diện rộng.

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

1. Phân đoạn gen nào của PRSV và cây đu đủ được lựa chọn để khuếch đại trong phản ứng multiplex RT-PCR?

Phản ứng multiplex RT-PCR khuếch đại đồng thời 2 phân đoạn: phân đoạn gen P3 của virus PRSV có kích thước 331 bp (sử dụng cặp mồi PRSV331-F/PRSV331-R) và phân đoạn gen actin của bộ gen đu đủ có kích thước 675 bp (sử dụng cặp mồi ACT-F/ACT-R làm đối chứng nội).

2. Điều kiện tối ưu cho phản ứng multiplex RT-PCR phát hiện PRSV được xác định như thế nào?

Điều kiện tối ưu gồm nhiệt độ bắt cặp là $56^\circ\text{C}$, nồng độ mỗi cặp mồi (PRSV331 và Actin) là $0{,}2,\mu\text{M}$. Chu trình luân nhiệt gồm 1 chu kỳ $94^\circ\text{C}$ (4 phút); 35 chu kỳ: $94^\circ\text{C}$ (40 giây) – $56^\circ\text{C}$ (30 giây) – $72^\circ\text{C}$ (1 phút); và 1 chu kỳ kéo dài $72^\circ\text{C}$ (5 phút).

3. Ngưỡng phát hiện tối đa của phương pháp sample pooling khi gộp RNA bệnh vào RNA khỏe là bao nhiêu?

Khi RNA mẫu bệnh được pha loãng trong RNA mẫu khỏe ở nồng độ chuẩn hóa $50 \pm 2\text{ ng/}\mu\text{L}$, ngưỡng phát hiện tối đa đạt mức gộp $1/1.000$ (tương ứng 1 mẫu đu đủ bệnh trong 999 mẫu đu đủ khỏe). Ở các mức pha loãng cao hơn ($1/10.000$ và $1/100.000$), phản ứng không phát hiện được virus.

4. Hiệu quả phát hiện virus PRSV trên 30 mẫu thực địa bằng phương pháp sample pooling thể hiện ra sao?

Quy trình gộp mẫu qua 3 bước phân cấp đã phát hiện chính xác 2 mẫu dương tính (mẫu 9 và mẫu 27) trong tổng số 30 mẫu thực địa thu thập tại Đức Huệ, Long An. Kết quả này trùng khớp 100% với phương pháp RT-PCR từng mẫu đơn lẻ nhưng chỉ tiêu tốn 17 phản ứng thay vì 30 phản ứng.

5. Phân đoạn gen PRSV khuếch đại ban đầu đạt độ tương đồng bao nhiêu phần trăm khi tra cứu trên NCBI?

Sản phẩm khuếch đại 613 bp từ cặp mồi PRSV613-F/R sau khi giải trình tự đã cho độ tương đồng 95,4% với trình tự gen PRSV mang mã số MT470189.1 đã được công bố trên Ngân hàng Gen NCBI bởi Tran và Yeh (2020).


Kết luận

Khóa luận tốt nghiệp đã xây dựng thành công quy trình kỹ thuật phát hiện Papaya ringspot virus (PRSV) gây bệnh đốm vòng trên cây đu đủ bằng phương pháp sample pooling kết hợp multiplex RT-PCR có tích hợp gen đối chứng nội actin. Nghiên cứu đã xác lập các thông số tối ưu gồm nhiệt độ ủ $56^\circ\text{C}$, nồng độ mồi $0{,}2,\mu\text{M}$ và đạt ngưỡng phát hiện mẫu gộp ở tỷ lệ $1/1.000$. Thử nghiệm trên 30 mẫu lá đu đủ thực địa tại Long An khẳng định quy trình có độ chính xác tương đương phương pháp RT-PCR đơn lẻ nhưng giúp giảm 43,3% số phản ứng xét nghiệm. Đây là cơ sở thực nghiệm tin cậy để triển khai các chương trình chẩn đoán, sàng lọc diện rộng phục vụ quản lý và phòng trừ bệnh đốm vòng hại đu đủ.