Tích hợp phương pháp đánh giá kiểu hình bằng ảnh với bản đồ QTL để làm tăng tính kháng về di truyền và cải tiến giống lúa kháng bệnh thối hạt do vi khuẩn
Nguồn: Jae-Ryoung Park, Gileung Lee, Hyun-Sook Lee, Seung Young Lee, Su-Kyung Ha, Kyeongmin Kang, Mina Jin, Jung-Pil Suh, Jeonghwan Seo, Young Mi Choi, Bong Choon Lee & Hyun-Su Park. 2026. Integrating
image-based phenotyping and QTL mapping to enhance genetic resistance
and accelerate breeding for bacterial grain rot resistance in rice. TAG; March 26 2026; vol. 139; article 108

Hình: Bệnh “bacterial grain rot” do vi khuẩn Burkholderia glumae
Bệnh thối hạt vi khuẩn BGR (bacterial grain rot) bởi Burkholderia glumae, là bệnh chính làm giảm năng suất lúa (Oryza sativa L.), do vậy, nó đe dọa an ninh lương thực. Phân tích kiểu hình theo phương pháp truyền thống hạn chế khi đánh giá tính kháng bệnh và hạn chế hiểu biết cơ sở di truyền. Các tác giả tích hợp phân tích kiểu hình trên cơ sở chụp ảnh với bản đồ QTL mapping để sàng lọc những QTLs cũng như gen ứng cử viên gắn với tính kháng B. glumae. Vi khuẩn B. glumae được chủng vào 189 cá thể cây con, dòng cận giao tái tổ hợp (RILs) của cặp lai Kele (kháng) và IS592BB (nhiễm), theo sau bởi kết quả phân tích bằng mắt và phân tích di truyền số lượng với công cụ “DAB staining”. Thông số kiểu hình, bao gồm cho điểm tính kháng ngoài đồng, tỷ lệ hoa lúa bị bệnh (%), cường độ nhuộm phẩm màu DAB, và tỷ số diện tích nhiễm bệnh (%), tất cả được đo đếm và được đưa số liệu vào chạy QTL mapping. Trên nhiễm sắc thể 1, trong quãng phân tử Chr01_24592710-Chr01_37274755, có 4 QTLs—qFRS1 [LOD: 5.98, giải thích biến thiên kiểu hình (PVE): 15.41%], qRDS1 (LOD: 5.29, PVE: 18.56%), qQDS1 (LOD: 9.58, PVE: 22.02%), và qRDA1 (LOD: 8.44, PVE: 31.51%) – được xác định là chồng lấp lên nhau. Sau khi tiến hành “fine-mapping”, người ta thu hẹp quãng phân tử Chr01_33472174-Chr01_33838140 và có tổng số 16 gen ứng cử viên được sàng lọc trong vùng này. Tho đó, OsBGq1 mã hóa protein có tên là “nucleotide-binding LRR receptor” (NLR) domain. OsBGq1 biểu hiện làm tăng đáng kể sự xâm nhiễm của B. glumae. Hơn nữa, loại hình RILs Kele của quãng phân tử Chr01_33472174-Chr01_33838140 làm tăng hoạt tính men “ROS-scavenging” và tích tụ phytoalexin khi B. glumae xâm nhiễm, làm bệnh tăng thêm. Tích hợp phương pháp đánh giá kiểu hình bằng ảnh với “QTL mapping” được đề nghị là một chỉ định có tính khách quan hơn nhằm xác định được gen đích điều khiển tính kháng BGR.
Xem: https://link.springer.com/article/10.1007/s00122-026-05216-7
GHI CHÚ
Hướng tiếp cận tích hợp di truyền xuôi (Forward Genetics) cùng với kiểu hình độ chính xác cao (High-throughput Phenotyping) và sinh hóa sinh lý là một quy trình chuẩn mực, thường thấy trên các công bố uy tín thuộc cơ sở dữ liệu NCBI (như Phytopathology, Theoretical and Applied Genetics). Tối ưu hóa nghiên cứu đánh giá tính kháng bệnh thối hạt vi khuẩn (BGR) do Burkholderia glumae trên lúa thì hướng cần học
Hiểu rõ vai trò và mối liên kết giữa 5 mắt xích cốt lõi:
Quần thể RILs (Recombinant Inbred Lines): Đây là quần thể cận giao tái tổ hợp (thường từ thế hệ F6 trở đi), được tạo ra từ cặp bố mẹ có kiểu hình tương phản (Bố kháng × Mẹ nhiễm hoặc ngược lại). Tính trạng kháng BGR là tính trạng số lượng (quantitative trait). Quần thể RILs có kiểu gen đồng hợp tử ổn định, cho phép lặp lại thí nghiệm ở nhiều môi trường để đánh giá chính xác tương tác Kiểu gen × Môi trường (G×E).
Kiểu hình dựa trên hình ảnh (Image-based Phenotyping): Thay vì đánh giá cảm quan bằng mắt thường dễ sai số, việc chụp ảnh kỹ thuật số (RGB hoặc Multispectral) giúp định lượng chính xác diện tích hạt bị thối, mức độ đổi màu (discoloration index) thông qua các thuật toán phân tách pixel (color segmentation).
Nhuộm DAB (3,3’-Diaminobenzidine Staining): Đây là công cụ di truyền số lượng ở cấp độ tế bào/sinh hóa. Khi vi khuẩn tấn công, cây sẽ kích hoạt phản ứng phòng vệ tạo ra các gốc oxy hóa phản ứng (ROS), đặc biệt là hydro peroxide (H₂O₂). DAB sẽ phản ứng với H₂O₂ dưới sự xúc tác của peroxidase nội bào tạo thành kết tủa màu nâu. Độ đậm và diện tích vết nâu phản ánh cường độ phản ứng miễn dịch (Oxidative Burst) của từng dòng RIL.
Bản đồ QTL (QTL Mapping): Kết hợp dữ liệu kiểu hình (tỷ lệ bệnh từ ảnh, mức độ ROS từ DAB) với dữ liệu kiểu gen (Marker SNPs/SSRs) để xác định các locus trên nhiễm sắc thể kiểm soát tính kháng.
Sàng lọc Gen ứng cử viên (Candidate Genes): Dựa trên tọa độ vật lý của vùng QTL, tra cứu cơ sở dữ liệu NCBI (như NCBI RefSeq, BLAST, Gene) để tìm các gen tiềm năng nằm trong vùng này (ví dụ: các gen mã hóa thụ thể nhận diện mẫu phân tử - PRRs, gen kháng NLR, hoặc các enzyme trong con đường phenylpropanoid).
(Step-by-Step Protocol)
Bước 1: Chuẩn bị vật liệu và Lây nhiễm bệnh (Inoculation)
Vật liệu: Trồng quần thể RILs và 2 dòng bố mẹ trong điều kiện nhà lưới đồng đều.
Chuẩn bị nguồn bệnh: Nuôi cấy chủng Burkholderia glumae kháng sinh học trên môi trường King’s B hoặc NA ở 37°C trong 24–48 giờ. Pha huyền phù vi khuẩn bằng nước cất vô trùng hoặc đệm PBS đạt mật độ OD600≈0.2 (khoảng 108 CFU/mL).
Lây nhiễm: Tiến hành lây nhiễm vào giai đoạn lúa trỗ bông rộ (Anthesis/Heading). Dùng phương pháp phun sương huyền phù vi khuẩn trực tiếp lên bông lúa hoặc tiêm tiêm trực tiếp vào hoa lúa. Giữ độ ẩm cao (>90%) và nhiệt độ cao (30–35°C) trong 24 giờ đầu để tạo điều kiện tối ưu cho vi khuẩn phát triển.
Bước 2: Phân tích Sinh hóa bằng Nhuộm DAB (Thời điểm sớm: 12 - 72 giờ sau nhiễm)
Mục đích là bắt giữ tín hiệu H₂O₂ trước khi mô bị hoại tử hoàn toàn.
Thu mẫu: Thu thập vỏ trấu (lemma/palea) hoặc hạt non ở các thời điểm: 12h, 24h, 48h, và 72h sau khi lây nhiễm (hpi).
Pha dung dịch: Hòa tan DAB trong nước (đã chỉnh pH về 3.8 bằng HCl) để đạt nồng độ 1 mg/mL. Thêm chất hoạt động bề mặt (như Tween-20 tỷ lệ 0.05%) để tăng độ thấm.
Ủ mẫu: Ngâm ngập mẫu lúa vào dung dịch DAB. Sử dụng bơm hút chân không (vacuum infiltration) trong 10-15 phút để dung dịch ngấm sâu. Ủ trong tối ở nhiệt độ phòng từ 8–12 giờ.
Tẩy màu: Đun sôi mẫu trong ethanol 96% hoặc hỗn hợp (Ethanol : Axit Acetic : Glycerol = 3:1:1) cho đến khi chất diệp lục bị tẩy sạch hoàn toàn, chỉ còn lại màu nâu của phức hợp DAB-H₂O₂.
Chụp ảnh định lượng: Chụp ảnh mẫu dưới kính hiển vi soi nổi hoặc máy quét phẳng độ phân giải cao. Dùng phần mềm ImageJ (Color Thresholding) để đo cường độ màu nâu (Mean Gray Value) hoặc % diện tích mô bị nhuộm nâu.
Bước 3: Thu thập và Xử lý Kiểu hình bằng Hình ảnh (Thời điểm muộn: 7 - 10 ngày sau nhiễm)
Thu hoạch: Khi bệnh biểu hiện rõ (hạt thối, lép đen), thu hoạch bông lúa, tách hạt.
Chụp ảnh dẹt (Flat-bed scanning hoặc Photo-box): Rải đều các hạt của từng dòng RIL lên một nền tương phản màu trắng hoặc đen chuẩn. Chụp ảnh bằng camera RGB có cố định điều kiện ánh sáng.
Phân tích ảnh (Digital Image Analysis): * Sử dụng công cụ mã nguồn mở như ImageJ (IJ-Macro), PlantCV (Python).
Chuyển đổi không gian màu sang hệ màu L∗a∗b∗ hoặc HSV để dễ phân tách màu “hạt lành” (vàng/xanh) và “mô bệnh” (nâu/đen/xám).
Chỉ số đầu ra: Tính toán tỷ lệ diện tích vết bệnh trên tổng diện tích hạt (Disease Severity=AreatotalArealesion×100%).
Bước 4: Phân tích Di truyền số lượng & Bản đồ QTL (QTL Mapping)
Để đảm bảo tính chính xác, tính toán hệ số di truyền theo nghĩa rộng (H2) nhằm xác định biến động kiểu hình có thực sự do di truyền hay không:
H2=Vg+eVge+reVeVg
Trong đó: Vg là biến dị kiểu gen, Vge là tương tác kiểu gen × môi trường, Ve là sai số ngẫu nhiên, e là số môi trường, r là số lần lặp lại.
Tích hợp dữ liệu: Lập bảng dữ liệu kiểu hình gồm: Tỷ lệ thối hạt (từ Bước 3) và Hàm lượng ROS/DAB (từ Bước 2) cho tất cả các dòng RILs.
Chạy phần mềm QTL: Sử dụng các gói phần mềm chuyên dụng như R/qtl, ICIMapping, hoặc QTL Cartographer.
Phương pháp: Áp dụng phương pháp Bản đồ khoảng cách tổ hợp (Composite Interval Mapping - CIM) hoặc Bản đồ khoảng cách toàn vẹn cải tiến (ICIM) để tìm các đỉnh LOD (Logarithm of odds) vượt ngưỡng ý nghĩa (thường thiết lập bằng kiểm tra hoán vị - Permutation test, N=1000).
Bước 5: Tra cứu NCBI sàng lọc Gen ứng cử viên (Candidate Gene Identification)
[Đỉnh QTL trên Bản đồ] ➔ [Tọa độ vật lý (bp) trên NST] ➔ [Tra cứu NCBI-Genome/RefSeq] ➔ [Sàng lọc chức năng Gen]
Xác định vùng neo (Anchoring): Lấy trình tự mồi (primer sequence) của các marker phân tử bám sát đỉnh QTL (Flanking markers).
Sử dụng NCBI BLAST: Khảo sát trình tự các marker này trên hệ gen chuẩn của lúa (Oryza sativa thuộc phân loài indica hoặc japonica, tùy thuộc vào nền di truyền của quần thể RIL của bạn). Xác định khoảng khoảng cách vật lý (ví dụ: từ nucleotide thứ 12,500,000 đến 13,200,000 trên Nhiễm sắc thể số 1).
Khai thác NCBI Gene & RefSeq: Truy cập vùng tọa độ này trên trình duyệt đồ họa hệ gen của NCBI (NCBI Genome Data Viewer). Liệt kê tất cả các locus gen nằm trong vùng này.
Sàng lọc chức năng sinh học: Ưu tiên chọn lọc các gen có chú giải chức năng liên quan tới:
Sản xuất hoặc phân hủy ROS (như Peroxidase, NADPH oxidase - liên kết trực tiếp với dữ liệu nhuộm DAB).
Các yếu tố điều hòa phiên mã (Transcription Factors như WRKY, MYB).
Các kinase truyền tín hiệu (như Receptor-like Kinases - RLKs).
TOOLS THỰC HIỆN
Phần Mềm Phân Tích Hình Ảnh (Kiểu hình & Nhuộm DAB) giúp bạn số hóa dữ liệu từ các bức ảnh chụp hạt lúa bị bệnh hoặc mẫu trấu nhuộm DAB thành các con số định lượng (diện tích, mật độ màu).
Phần Mềm Di Truyền Số Lượng & Bản Đồ QTL Sau khi có dữ liệu kiểu hình (từ ImageJ) và dữ liệu kiểu gen (từ các chỉ thị phân tử Marker), bạn cần các phần mềm này để tìm vị trí gen (QTL).
QTL IciMapping http://www.isbreeding.net
Mục đích: Phần mềm chạy trên hệ điều hành Windows do Viện Khoa học Cây trồng Trung Quốc (CAAS) và CIMMYT phát triển. Nó tích hợp thuật toán ICIM (Inclusive Composite Interval Mapping) vô cùng tối ưu cho quần thể RILs, giao diện đồ họa rất trực quan, dễ dùng cho người mới bắt đầu.
R / RStudio (Gói gói lệnh qtl và qtl2) https://rqtl.org hoặc https://kbroman.org/qtl2/
Mục đích: “Tiêu chuẩn vàng” trong các công bố quốc tế. Gói lệnh R/qtl mạnh mẽ trong việc xử lý dữ liệu lớn, sàng lọc marker lỗi, hoán vị (permutation test) để tìm ngưỡng LOD chính xác và vẽ bản đồ QTL chất lượng cao để in ấn.
Windows QTL Cartographer https://statgen.ncsu.edu/qtlcart/WQTLCart.htm
3. Hệ Thống Website & Cơ Sở Dữ Liệu Tra Cứu (Genomics Databases)
NCBI (National Center for Biotechnology Information https://www.ncbi.nlm.nih.go * NCBI Genome Data Viewer (GDV): Trình duyệt đồ họa giúp bạn nhìn thấy trực quan tất cả các gen nằm lân cận vùng đỉnh QTL. * NCBI BLAST: Dùng để dán trình tự mồi (primer) của chỉ thị phân tử phân tích QTL nhằm xác định vị trí vật lý chính xác (tọa độ cơ sở - base pair) trên hệ gen lúa.
RAP-DB (Rice Annotation Project Database) https://rapdb.dna.affrc.go.jp
MSU Rice Genome Annotation Project http://rice.uga.edu
RiceData (China Rice Data Center) http://www.ricedata.cn
Mục đích: Đây là một kho báu cho nghiên cứu lúa. Website này tổng hợp toàn bộ các gen lúa đã được nhân bản và chứng minh chức năng trên thế giới. Khi bạn tìm thấy một danh sách gen ứng cử viên từ NCBI, hãy ném mã gen vào RiceData để xem gen đó đã có ai nghiên cứu về tính kháng bệnh hay liên quan đến con đường oxy hóa (ROS) chưa.
Bảng Tổng Hợp Công Cụ Theo Tiến Trình
Mẹo nhỏ từ chuyên gia: Khi viết bài báo quốc tế, ở phần phương pháp nghiên cứu (Materials and Methods), bạn hãy ghi rõ phiên bản phần mềm (ví dụ: QTL IciMapping version 4.2) và đính kèm link URL hoặc trích dẫn (citation) bài báo gốc của phần mềm đó để tăng tính thuyết phục về mặt học thuật.
Trong các nhóm công cụ trên, bạn đã có kinh nghiệm sử dụng phần mềm nào chưa (ví dụ như đã quen dùng môi trường R hay thích sử dụng giao diện nhấp chuột như IciMappin