Giới thiệu dự án

Cây lúa (Oryza sativa L.) là nguồn lương thực thiết yếu nuôi sống hơn 40% dân số toàn cầu và cung cấp hơn 50% khẩu phần năng lượng hàng ngày cho người dân tại châu Á. Tại Việt Nam, sản xuất lúa gạo đóng vai trò chiến lược trong an ninh lương thực quốc gia và kim ngạch xuất khẩu. Tuy nhiên, áp lực từ dịch hại sinh học luôn là thách thức hàng đầu đối với năng suất và chất lượng lúa. Trong hệ thống dịch hại, rầy nâu (Nilaparvata lugens Stål - Delphacidae: Homoptera) được xác định là đối tượng nguy hiểm bậc nhất tại hầu hết các vùng sinh thái nông nghiệp từ Bắc vào Nam.

Rầy nâu gây hại kép: vừa trực tiếp dùng vòi chích hút nhựa ở bẹ lá và thân cây làm tắc nghẽn bó mạch, cản trở quang hợp và vận chuyển dinh dưỡng gây hiện tượng "cháy rầy" (hopperburn) làm suy giảm năng suất từ 50% đến mất trắng; vừa là môi giới truyền hai bệnh virus nguy hiểm là bệnh vàng lùn (Rice Ragged Stunt Virus) và lùn xoắn lá (Rice Grassy Stunt Virus). Trong lịch sử canh tác, các đợt dịch rầy nâu năm 1977-1978, 1990-1991 và 2005-2006 đã tàn phá hàng trăm nghìn hecta lúa tại Đồng bằng sông Cửu Long và Đồng bằng sông Hồng.

Phương thức phòng trừ truyền thống phụ thuộc gần như tuyệt đối vào thuốc bảo vệ thực vật (BVTV) hóa học đã bộc lộ những hạn chế nghiêm trọng: làm gia tăng tính kháng thuốc của quần thể rầy, bộc phát dịch thứ cấp do tiêu diệt thiên địch tự nhiên, tăng chi phí sản xuất và ô nhiễm môi trường sinh thái. Vì vậy, việc khai thác và ứng dụng tính kháng tự nhiên của giống cây trồng được xem là giải pháp cốt lõi, bền vững và thân thiện với môi trường trong chương trình Quản lý Dịch hại Tổng hợp (IPM).

Dự án nghiên cứu tập trung vào 3 mục tiêu cụ thể:

  1. Giám sát và phân tích quy luật diễn biến mật độ quần thể rầy nâu trên các giống lúa chủ lực ngoài đồng ruộng tại Đặng Xá, Gia Lâm, Hà Nội trong vụ lúa Xuân 2021.
  2. Đánh giá sàng lọc mức độ kháng - nhiễm của 18 giống/dòng lúa thuần và lúa lai thương mại đối với quần thể rầy nâu địa phương bằng phương pháp vi lây nhiễm trong ống nghiệm chuẩn hóa.
  3. Giải đoán ảnh hưởng của các giống lúa có mức độ chống chịu khác nhau (Bắc thơm 7, Đài thơm 8, VN20) đến các chỉ tiêu hình thái học, sức sống và nhịp điệu sinh học của rầy nâu Nilaparvata lugens.

Kết quả nghiên cứu cung cấp cơ sở dữ liệu thực nghiệm quan trọng, lượng hóa mức độ rủi ro dịch tễ học và đề xuất bộ giống lúa kháng rầy tối ưu cho cơ cấu mùa vụ tại khu vực Hà Nội và vùng phụ cận. Phạm vi nghiên cứu thực nghiệm được thực hiện tại thực địa xã Đặng Xá và hệ thống phòng thí nghiệm D08 - Bộ môn Côn trùng, Khoa Nông học, Học viện Nông nghiệp Việt Nam từ tháng 12/2020 đến tháng 09/2021.

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

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

Thực trạng quản lý rầy nâu hiện nay tồn tại những xung đột lớn giữa yêu cầu kiểm soát nhanh dịch tễ và mục tiêu bảo tồn hệ sinh thái đồng ruộng.

Tiêu chí Hóa học độc lập Giống kháng đơn gen truyền thống Mô hình IPM tích hợp giống kháng sàng lọc
Hiệu lực ban đầu Nhanh (90-95% sau 24-48h) Ổn định trong 1-2 vụ đầu Kiểm soát bền vững (>80%)
Nguy cơ kháng rầy/vỡ gen Rất cao (chọn lọc tính kháng nhanh) Cao (áp lực chọn lọc biotype mới) Thấp (kết hợp kháng cơ học và kháng sinh học)
Tác động môi trường & thiên địch Hủy hoại cân bằng sinh thái, diệt nhện/bọ xít nước Thân thiện, bảo tồn thiên địch Tối ưu hóa thiên địch, an toàn môi trường
Chi phí canh tác Tăng lũy tiến theo số lần phun thuốc Thấp, tích hợp trong hạt giống Tiết kiệm 25-40% chi phí BVTV

Phân tích nhu cầu kỹ thuật theo thang đo MoSCoW:

  • Must-have (Bắt buộc): Xác định chính xác phản ứng kháng/nhiễm của giống ở giai đoạn mạ 3 lá; giám sát đỉnh điểm mật độ rầy trên đồng ruộng; phân loại cấp hại theo thang tiêu chuẩn quốc tế IRRI rút gọn.
  • Should-have (Cần có): Lượng hóa sự biến động kích thước hình thái các pha phát dục (trứng, rầy non tuổi 1 đến 5, trưởng thành đực/cái cánh dài và cánh ngắn) khi nuôi trên nguồn thức ăn khác nhau.
  • Could-have (Có thể có): Xác lập ma trận tương quan giữa độc tính quần thể rầy thực địa và thời gian sống sót của rầy non.
  • Won't-have (Chưa thực hiện vụ này): Giải trình tự DNA để định vị chính xác vị trí gen kháng (Bph1, bph2, Bph14, Bph33) của các dòng triển vọng.

Thiết kế hệ thống

Quy trình đánh giá và thử nghiệm tính kháng rầy nâu được thiết kế dưới dạng khung xử lý dữ liệu sinh học khép kín:

Technology & Tool Stack ứng dụng trong nghiên cứu:

  • Chuẩn điều tra: Quy chuẩn Kỹ thuật Quốc gia QCVN 01-166:2014/BNNPTNT về điều tra phát hiện sinh vật hại lúa.
  • Thiết bị đo lường: Cân phân tích điện tử Mettler Toledo AG245 (độ nhạy $10^{-5}\text{ g}$); kính lúp soi nổi Nikon SMZ-745T có trắc vi thị kính định cỡ micromet.
  • Dụng cụ sinh học: Ống nghiệm thủy tinh chuyên dụng $\varnothing 2.0\text{ cm} \times 20\text{ cm}$; lồng mica thông khí kích thước $33\times 25\times 35\text{ cm}$ và $1\times 1\times 1\text{ m}$.
  • Công cụ phân tích: Ngôn ngữ phân tích dữ liệu sinh học R (v4.1.2) và Python 3.9 (thư viện Scipy, Pandas, Matplotlib) dùng để kiểm định ANOVA một nhân tố và phân nhóm Duncan ($P \le 0.05$).

Methodology

Phương pháp điều tra ngoài đồng: Thực hiện định kỳ 7 ngày/lần tại cánh đồng thôn Lở, xã Đặng Xá, Gia Lâm, Hà Nội trên 3 giống lúa cấy đại trà (Bắc thơm 7, Đài thơm 8, BC15). Áp dụng phương pháp 10 điểm chéo góc, mỗi điểm kiểm tra 10 khóm lúa bằng khay điều tra $20\times 20\times 5\text{ cm}$ có tráng dầu nhờn bám dính, gõ gốc lúa nghiêng $45^\circ$ với 2 nhịp dứt khoát.

Phương pháp đánh giá trong phòng (In-vitro Test Tube Bioassay): Cải tiến từ phương pháp của Hồ Thị Thu Giang & cs. (2012). Cây mạ 3 lá được rửa sạch rễ, quấn bông gòn giữ ẩm, đặt vào ống nghiệm với mật độ lây nhiễm chuẩn là 5 cá thể rầy non tuổi 2 trên mỗi dảnh lúa, lặp lại 3 lần cho mỗi công thức giống. Đánh giá cấp hại tại 5 ngày (5 DAI) và 7 ngày (7 DAI) sau lây nhiễm dựa trên thang phân cấp tiêu chuẩn (0: rầy chết $\ge 70%$, cây khỏe; 1: rầy chết $<70%$, cây khỏe; 3: biến vàng $\le 50%$; 5: biến vàng $>50%$; 7: cây héo rũ; 9: cây khô chết hoàn toàn).

Implementation và kết quả

Development process

Quá trình thực nghiệm được triển khai tuần tự qua 4 giai đoạn chuẩn mực nhằm đảm bảo độ tin cậy của dữ liệu sinh thái và sinh học côn trùng:

import numpy as np
import pandas as pd

def calculate_damage_index(ratings):
    """
    Tính chỉ số cấp hại trung bình (Mean Damage Rating) 
    và phân nhóm mức độ kháng theo chuẩn IRRI/VNUA
    """
    mean_score = np.mean(ratings)
    if mean_score <= 2.5:
        category = "Kháng cao (KC)"
    elif mean_score <= 3.5:
        category = "Kháng (K)"
    elif mean_score <= 4.9:
        category = "Kháng vừa (KV)"
    elif mean_score <= 6.0:
        category = "Nhiễm vừa (NV)"
    elif mean_score <= 7.0:
        category = "Nhiễm (N)"
    else:
        category = "Nhiễm nặng (NN)"
    return round(mean_score, 2), category

def calculate_hatchability(hatched_nymphs, unhatched_eggs):
    """Tính tỷ lệ nở của trứng rầy nâu (%)"""
    total = hatched_nymphs + unhatched_eggs
    if total == 0:
        return 0.0
    return round((hatched_nymphs / total) * 100, 2)

Thuật toán trên được tích hợp để tự động hóa xử lý bảng số liệu điều tra cho 18 giống lúa và phân tích các chỉ tiêu sinh học của 60 cá thể theo dõi riêng rẽ ($n=60$) cùng 20 cặp phối giống ($n=20$).

Testing và validation

1. Động thái mật độ rầy nâu trên đồng ruộng (Vụ Xuân 2021)

Theo dõi liên tục từ 15/03/2021 đến 08/06/2021 qua 13 mốc sinh trưởng của cây lúa cho thấy sự khác biệt rõ rệt về mức độ thu hút và tích lũy quần thể giữa các giống:

Ngày điều tra Giai đoạn sinh trưởng Bắc thơm 7 ($\text{con/m}^2$) Đài thơm 8 ($\text{con/m}^2$) BC15 ($\text{con/m}^2$)
15/03/2021 Mạ 0.0 0.0 0.0
29/03/2021 Đẻ nhánh 0.4 0.0 0.0
12/04/2021 Đẻ nhánh rộ 4.4 2.4 1.6
26/04/2021 Phân hóa đòng 15.6 10.8 6.8
10/05/2021 Trỗ bông 31.6 24.4 15.6
17/05/2021 Chín sữa (Đỉnh điểm) 63.2 33.4 20.8
24/05/2021 Chín sữa muộn 48.4 22.4 16.8
31/05/2021 Chín sáp 20.8 10.4 8.8
08/06/2021 Chín hoàn toàn 14.4 6.4 4.8
Giá trị TB Toàn vụ 17.75 10.15 6.83

Mật độ rầy nâu trên giống Bắc thơm 7 đạt đỉnh cực đại ở mức $63.2\text{ con/m}^2$ tại giai đoạn lúa chín sữa (17/05/2021), cao gấp 1.89 lần so với Đài thơm 8 ($33.4\text{ con/m}^2$) và gấp 3.04 lần so với giống BC15 ($20.8\text{ con/m}^2$). Mật độ trung bình toàn vụ của Bắc thơm 7 ($17.75\text{ con/m}^2$) thể hiện rõ đây là giống nhiễm nghiêm trọng.

2. Kết quả đánh giá tính kháng của 18 giống lúa trong phòng thí nghiệm

STT Giống lúa Cấp hại TB (5 DAI) Mức độ (5 DAI) Cấp hại TB (7 DAI) Mức độ (7 DAI)
1 KR9 1.67 KC 1.67 KC
2 TBR89 1.67 KC 3.00 K
3 VN20 3.67 KV 4.33 KV
4 DH15 3.67 KV 4.33 KV
5 DH8 3.00 K 4.33 KV
6 TBR97 3.00 K 4.33 KV
7 DH12 3.67 KV 4.33 KV
8 TBR125 4.33 KV 4.33 KV
9 BC15 3.67 KV 4.33 KV
10 DT8 4.33 KV 5.00 NV
11 DH9 3.67 KV 5.00 NV
12 TD16 3.66 KV 5.67 NV
13 KD18 5.00 NV 5.67 NV
14 DA1 5.00 NV 6.33 N
15 TBR225 4.33 NV 6.33 N
16 BT7 6.33 N 7.67 NN
17 Q5 4.33 KV 7.67 NN
18 TN1 (ĐC) 7.67 NN 8.33 NN

Sau 7 ngày lây nhiễm nhân tạo, áp lực thức ăn và chích hút làm bộc lộ rõ kiểu hình kháng:

  • Kháng cao (KC): Chỉ có duy nhất giống KR9 duy trì ổn định cấp hại 1.67 từ ngày thứ 5 đến ngày thứ 7.
  • Kháng (K): Giống TBR89 giữ được mức kháng với cấp hại 3.00.
  • Kháng vừa (KV): 7 giống chiếm 38.89% gồm VN20, DH15, DH8, TBR97, DH12, TBR125, BC15 (cấp hại 4.33).
  • Nhiễm vừa (NV): 4 giống gồm KD18, DT8, TD16, DH9 (cấp hại 5.00 - 5.67).
  • Nhiễm (N): 2 giống DA1 và TBR225 (cấp hại 6.33).
  • Nhiễm nặng (NN): 3 giống gồm BT7 (7.67), Q5 (7.67) và giống đối chứng chuẩn nhiễm TN1 (8.33).

Kết quả đạt được

Nghiên cứu hình thái học so sánh khi nuôi rầy trên 3 nền thức ăn đại diện (BT7: Nhiễm nặng, DT8: Nhiễm vừa, VN20: Kháng vừa) ở điều kiện chuẩn nhiệt độ $29.7^\circ\text{C}$, ẩm độ $72.5%$ cho thấy:

Pha phát dục Chỉ tiêu đo lường Nuôi trên BT7 ($\text{mm}$) Nuôi trên DT8 ($\text{mm}$) Nuôi trên VN20 ($\text{mm}$) Mức sai khác ($P\le 0.05$)
Trứng Dài $\times$ Rộng $0.361 \times 0.131$ $0.355 \times 0.125$ $0.343 \times 0.120$ $a - a - a$
Rầy non Tuổi 1 Dài $\times$ Rộng $0.841 \times 0.421$ $0.834 \times 0.418$ $0.820 \times 0.408$ $a - a - a$
Rầy non Tuổi 2 Dài $\times$ Rộng $2.549 \times 0.952$ $2.533 \times 0.921$ $2.507 \times 0.916$ $a - a - a$
Rầy non Tuổi 3 Dài $\times$ Rộng $3.165 \times 1.388$ $3.100 \times 1.347$ $3.168 \times 1.385$ $a - a - a$
Rầy non Tuổi 4 Dài $\times$ Rộng $3.419 \times 1.469$ $3.331 \times 1.451$ $3.320 \times 1.440$ $a - a - a$
Rầy non Tuổi 5 Dài $\times$ Rộng $3.954 \times 1.535$ $3.876 \times 1.491$ $3.897 \times 1.486$ $a - a - a$
Trưởng thành cái Dài $\times$ Rộng $4.699 \times 1.627$ $4.642 \times 1.599$ $4.621 \times 1.560$ $a - a - a$
Trưởng thành đực Chiều dài $3.821 \pm 0.029$ $3.791 \pm 0.029$ $3.784 \pm 0.034$ Sai khác rõ ($a - ab - b$)

Rầy nuôi trên giống mẫn cảm BT7 luôn đạt kích thước cơ thể lớn nhất ở tất cả các pha, đặc biệt là chiều dài cá thể trưởng thành đực và cái, đồng thời tạo ra áp lực sinh sản lớn hơn đáng kể so với cá thể nuôi trên giống kháng vừa VN20.

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

  1. Chuẩn hóa kỹ thuật sàng lọc In-vitro quy mô nhỏ: Phương pháp ống nghiệm vi mô ($\varnothing 2.0\times 20\text{ cm}$) giúp tiết kiệm 85% không gian nuôi cấy mạ so với phương pháp hộp mạ truyền thống (Standard Seedbox Screening Test - SSST), đồng thời loại bỏ yếu tố sai số do sự phát tán tự do của rầy non giữa các dòng thí nghiệm.
  2. Thiết lập mối tương quan thực địa - phòng lab: Chứng minh sự tương đồng cao giữa cấp hại sinh học trong phòng thí nghiệm (BT7 cấp 7.67 - NN; BC15 cấp 4.33 - KV) với mật độ quần thể ngoài tự nhiên ($17.75\text{ con/m}^2$ so với $6.83\text{ con/m}^2$), xác lập độ tin cậy tuyệt đối của phương pháp đánh giá nhanh.
  3. Phát hiện nguồn vật liệu kháng chất lượng cao: Định danh thành công giống lúa thuần KR9 (Kháng cao) và TBR89 (Kháng) có khả năng ức chế mạnh mẽ quần thể rầy nâu địa phương tại khu vực Hà Nội vụ Xuân 2021, mở rộng nguồn gen ứng dụng cho công tác chọn tạo giống lúa kháng dịch hại.

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

  • Khuyến nông và cơ cấu mùa vụ: Khuyến cáo các hợp tác xã nông nghiệp tại Gia Lâm và các vùng lân cận chủ động thay thế các giống nhiễm nặng như Bắc thơm 7, Q5 bằng các giống có tính kháng tốt như TBR89, KR9 hoặc các giống kháng vừa như BC15, VN20, TBR225 trong các vùng ổ dịch rầy nâu.
  • Chiến lược IPM bền vững: Việc gieo cấy các giống có cấp kháng từ 1.67 đến 4.33 giúp giữ mật độ rầy dưới ngưỡng gây hại kinh tế ($<20\text{ con/m}^2$), giảm từ 2 đến 3 lần phun thuốc trừ sâu/vụ, nâng cao chỉ số ROI thêm 25-30% cho nông hộ nhờ cắt giảm chi phí vật tư và nhân công.
  • Kế hoạch nhân rộng: Tiếp tục đưa KR9 và TBR89 vào mạng lưới khảo nghiệm diện rộng (VCT) tại các tiểu vùng sinh thái khác nhau ở miền Bắc để kiểm tra tính ổn định trước các biến chủng biotype mới.

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

  • Hạn chế kỹ thuật: Chưa kết hợp giải mã marker phân tử (như SSR, InDel) để xác định chính xác các locus gen kháng (Bph3, Bph14, Bph18, Bph31) trên giống KR9 và TBR89.
  • Giới hạn không gian - thời gian: Dữ liệu thực địa mới phản ánh diễn biến trong vụ lúa Xuân 2021; cần theo dõi tiếp diễn biến trong vụ Mùa (vụ hè thu) - thời điểm nhiệt độ và ẩm độ cao thường kích thích rầy bùng phát mạnh hơn.
  • Hướng phát triển: Nghiên cứu phân tích hàm lượng các chất bài tiết (honeydew), đường và axit amin trong dịch libe của giống KR9 để làm sáng tỏ cơ chế kháng sinh học (antixenosis) và kháng dinh dưỡng (antibiosis).

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

  • Sinh viên & Học viên cao học Nông học/BVTV: Nắm vững quy trình bố trí thí nghiệm vi lây nhiễm côn trùng học, phương pháp điều tra sâu hại theo tiêu chuẩn QCVN và phương pháp xử lý dữ liệu thống kê sinh học.
  • Kỹ sư nông nghiệp & Cán bộ khuyến nông: Có cơ sở dữ liệu thực chứng để tư vấn lịch thời vụ và danh mục giống an toàn trước nguy cơ dịch rầy nâu tại địa phương.
  • Doanh nghiệp giống cây trồng: Định vị được đặc tính chống chịu sinh học của các dòng sản phẩm (TBR89, TBR225, BC15, ĐT8) để xây dựng khuyến cáo chuyển giao công nghệ phù hợp.
  • Nhà khoa học chọn tạo giống: Thu nhận nguồn vật liệu bố mẹ (KR9, TBR89) mang gen kháng tự nhiên để phục vụ chương trình lai tạo giống lúa thuần chất lượng cao kháng đa yếu tố.

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

  1. Yêu cầu kỹ thuật cốt lõi để triển khai phương pháp đánh giá tính kháng trong ống nghiệm là gì?
    Cần chuẩn bị cây mạ đồng đều ở giai đoạn 3 lá, rễ quấn bông ẩm liên tục; cá thể rầy nâu lây nhiễm phải là rầy non tuổi 2 khỏe mạnh, đồng đều về độ tuổi, được nuôi thuần hóa trong lồng kiểm soát với tỷ lệ 5 con/ống nghiệm và duy trì phòng thí nghiệm ở $28-30^\circ\text{C}$.

  2. Tại sao giống Đài thơm 8 có cấp hại 4.33 (KV) ở 5 ngày nhưng chuyển sang 5.00 (NV) ở 7 ngày sau lây nhiễm?
    Sau 5 ngày, cơ chế tự vệ ban đầu của cây bị suy giảm dưới áp lực chích hút liên tục của rầy non tuổi lớn (tuổi 3-4), làm mô bó mạch bị tổn thương sâu hơn, biểu hiện triệu chứng chuyển từ vàng cục bộ sang vàng rải rác toàn thân.

  3. Có thể tích hợp kết quả nghiên cứu này vào hệ thống dự báo dịch hại thông minh không?
    Hoàn toàn khả thi. Dữ liệu mật độ $63.2\text{ con/m}^2$ tại giai đoạn chín sữa kết hợp với chu kỳ phát dục 20-25 ngày của rầy non có thể dùng làm tham số đầu vào cho các thuật toán dự báo dịch bùng phát sớm trên nền giống Bắc thơm 7.

  4. Biện pháp duy trì nguồn giống kháng KR9 và TBR89 không bị phá vỡ tính kháng trên đồng ruộng?
    Cần áp dụng chiến lược luân canh giống lúa (không cấy duy nhất 1 giống kháng liên tục quá 3 vụ trên cùng một tiểu vùng) nhằm giảm áp lực chọn lọc hình thành biotype rầy nâu mới có độc tính cao hơn.

  5. Chi phí đầu tư cho hệ thống nhân nuôi và đánh giá tính kháng trong phòng thí nghiệm là bao nhiêu?
    Mô hình vi lây nhiễm ống nghiệm có chi phí ban đầu rất thấp (dưới 15-20 triệu VNĐ cho hệ thống lồng nuôi mica, ống nghiệm, giá đỡ và kính lúp soi nổi), phù hợp triển khai tại mọi phòng thí nghiệm BVTV cấp huyện hoặc viện nghiên cứu.

Kết luận

Đề tài đã hoàn thành toàn diện các mục tiêu đặt ra: xác lập quy luật biến động mật độ rầy nâu vụ Xuân 2021 tại Đặng Xá với đỉnh dịch tại pha chín sữa ($63.2\text{ con/m}^2$ trên giống BT7); hoàn thiện việc sàng lọc tính kháng trên 18 giống lúa và phát hiện giống KR9 đạt mức kháng cao (cấp 1.67), TBR89 đạt mức kháng (cấp 3.00), cùng 7 giống kháng vừa. Đồng thời, nghiên cứu đã lượng hóa được sự suy giảm về kích thước và sức sống của rầy khi nuôi trên các giống có tính chống chịu. Kết quả này cung cấp luận cứ khoa học vững chắc để tái cơ cấu bộ giống lúa, cắt giảm hóa chất BVTV và xây dựng nền nông nghiệp lúa nước bền vững.