Tổng quan nghiên cứu

Lúa gạo là nguồn lương thực thiết yếu nuôi sống hơn 3,5 tỷ người trên toàn cầu, đóng góp khoảng 20% tổng lượng calo tiêu thụ hàng ngày của nhân loại. Tại Việt Nam, cây lúa giữ vị trí chiến lược hàng đầu khi chiếm tới 90% tổng sản lượng ngũ cốc và giải quyết việc làm cho gần 80% lao động nông thôn. Mặc dù sản lượng lúa cả nước năm 2017 đạt 42,76 triệu tấn với kim ngạch xuất khẩu 2,62 tỷ USD (tương đương 5,79 triệu tấn), ngành lúa gạo vẫn đối mặt với thách thức lớn khi năng suất thực tế bình quân chỉ đạt 5,55 tấn/ha, thấp hơn nhiều so với tiềm năng sinh học 10 tấn/ha. Một trong những nguyên nhân hàng đầu gây thất thoát năng suất là bệnh đạo ôn do nấm Magnaporthe oryzae gây ra, với mức độ thiệt hại có thể lên đến 80% tại các vùng dịch bệnh nghiêm trọng và ước tính gây tổn thất hơn 55 triệu USD mỗi năm tại khu vực Đông Nam Á.

Trước thực trạng biến đổi khí hậu và sự biến đổi liên tục của các chủng nấm bệnh, việc phát triển các giống lúa có tính kháng bền vững thông qua quy tụ nhiều gen kháng là giải pháp sinh học cấp thiết nhất. Đề tài tập trung giải quyết vấn đề nhận diện và khai thác nguồn gen kháng quý trong tập đoàn lúa bản địa Việt Nam. Mục tiêu cụ thể của luận văn là tầm soát hai gen kháng chủ lực gồm gen Pit và gen Pi21 trên 17 giống lúa bản địa đã được giải mã hệ gen, đồng thời thiết kế thành công các cặp chỉ thị ADN đặc hiệu phục vụ công tác chọn giống phân tử. Nghiên cứu được triển khai tại Bộ môn Kỹ thuật di truyền thuộc Viện Di truyền Nông nghiệp trong giai đoạn từ tháng 5/2018 đến tháng 4/2019. Kết quả đề tài mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp hệ thống chỉ thị phân tử có độ tin cậy đạt 100% trên tập mẫu thử nghiệm, giúp rút ngắn từ 3 đến 5 thế hệ chọn giống so với phương pháp lai truyền thống.

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 thuyết tương tác gen đối gen (Gene-for-gene hypothesis) của Flor và lý thuyết trung tâm (Central Dogma) của sinh học phân tử, khẳng định sự tương tác đặc hiệu giữa sản phẩm gen kháng của cây chủ và gen độc tố của nấm bệnh. Cấu trúc của hầu hết các gen kháng đạo ôn thuộc nhóm protein chứa vùng liên kết nucleotide (NBS) và trình tự lặp giàu leucine (LRR), đóng vai trò thụ thể nhận diện tác nhân gây bệnh. Nghiên cứu tập trung vào 4 khái niệm nền tảng:

Đa hình đơn nucleotide (SNP) và đột biến chèn/mất đoạn (InDel) đại diện cho các biến dị di truyền phong phú nhất trong hệ gen với tần suất phân bố dày đặc, đóng vai trò làm chỉ thị phân tử thế hệ mới. Chọn giống nhờ chỉ thị phân tử (MAS) và lai hồi giao kết hợp chỉ thị (MABC) là phương pháp chọn lọc cá thể mang gen mục tiêu ở giai đoạn cây con dựa trên kiểu gen thay vì đánh giá kiểu hình ngoài đồng. Hệ thống phân loại gen kháng gồm gen kháng chính (Major gene - Pit) và gen kháng phụ hay gen lặn (Minor/Recessive gene - Pi21), trong đó gen Pit nằm trên nhiễm sắc thể số 1 với vùng mã hóa dài 2970 cặp bazơ nitơ (bp), và gen Pi21 nằm trên nhiễm sắc thể số 12 giúp hạn chế sự phát triển xâm nhiễm của nấm bệnh. Công nghệ giải trình tự thế hệ mới (NGS) cung cấp nền tảng giải mã song song hàng triệu đoạn ADN với độ chính xác đạt trên 99,9% và điểm chất lượng dữ liệu vượt ngưỡng Q30.

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

Nghiên cứu sử dụng cỡ mẫu gồm 17 giống lúa bản địa đặc sản của Việt Nam như Khấu điển lư, Thơm Lài, OM6377, Tám xoan Hải Hậu, Nàng quớt biển. Phương pháp chọn mẫu có chủ đích được áp dụng nhằm khai thác tối đa nguồn biến dị di truyền bản địa quý hiếm từ cơ sở dữ liệu của Trung tâm Tài nguyên Thực vật. Quy trình phân tích kết hợp giữa tin sinh học in silico và thực nghiệm in vitro:

Tầm soát dữ liệu hệ gen sử dụng phần mềm NextGENe phiên bản 2.3 và MEGA phiên bản 6.0 để dóng hàng đa chuỗi, căn chỉnh trình tự FastQ với chuỗi tham chiếu LOC_Os01g05620 (gen Pit). Lý do lựa chọn công cụ này là khả năng xử lý dữ liệu dung lượng lớn với tốc độ cao và phát hiện chính xác các đột biến SNP và InDel vi mô. Thiết kế mồi phân tử ứng dụng phần mềm Vector NTI Advance phiên bản 11.0 để thiết kế các cặp mồi có chiều dài 21-24 nucleotide, nhiệt độ nóng chảy tối ưu trong khoảng 55°C đến 60°C và hàm lượng GC từ 45% đến 55%.

Thực nghiệm kiểm chứng tiến hành tách chiết ADN tổng số từ lá lúa 3 tuần tuổi theo quy trình CTAB 1% cải tiến. Phản ứng PCR được thực hiện trong tổng thể tích 15 µl (chứa 10 ng ADN khuôn, 0,15 µM mỗi mồi, 0,2 mM dNTPs, 2,5 mM MgCl2 và 0,25 đơn vị Taq polymerase) trên máy luân nhiệt Veriti 96 giếng. Sản phẩm khuếch đại được phân tích qua điện di gel polyacrylamide 6,0% và gel agarose nhuộm ethidium bromide. Toàn bộ quy trình thực nghiệm được triển khai từ tháng 5/2018 đến tháng 4/2019.

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

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

Thứ nhất, phát hiện 100% (17/17) giống lúa bản địa mang trình tự tương đồng cao với locus gen Pit và gen Pi21 trên ngân hàng gen tham chiếu quốc tế NCBI.

Thứ hai, đối với gen Pit, nghiên cứu tách biệt rõ 2 nhóm kiểu gen chính trong tập đoàn lúa. Nhóm thứ nhất gồm 9 giống (chiếm tỷ lệ 52,94%) như Khấu điển lư, Thơm Lài, Tám xoan Hải Hậu chỉ thiếu 3 nucleotide ở vị trí 543-545 so với chuẩn 2970 bp, dịch mã thành protein hoàn chỉnh dài 987 amino acid. Nhóm thứ hai gồm 8 giống còn lại (chiếm tỷ lệ 47,06%) như Nếp mặn, Lúa gốc đỏ mang đột biến mất đoạn kép gồm 3 nucleotide đầu và 17 nucleotide ở vị trí 2647-2664 (tổng mất 20 bp), dịch mã thành chuỗi protein ngắn chỉ từ 980 đến 981 amino acid.

Thứ ba, đối với gen Pi21, ghi nhận 9/17 giống có chiều dài chuỗi tương đồng 100% với gen tham chiếu, trong khi 8/17 giống (tỷ lệ 47,06%) xuất hiện đột biến mất đoạn 24 nucleotide trong vùng mã hóa.

Thứ tư, thiết kế và kiểm chứng thành công hai cặp chỉ thị ADN đặc hiệu gồm cặp mồi Pit (PitF: TCTTAACGACTGCCCAAAACT / PitR: AGATTACAGAGATTGGAAATTCTC) và cặp mồi Pi21 (Pi21F: GACCCGGAGAAGCTGTGCAAGA / Pi21R: CTTTGGGCAGCACCACTGCC), cho kết quả nhân bản đạt độ đặc hiệu tuyệt đối 100% trên bản điện di phân tử.

Thảo luận kết quả

Sự sai khác về số lượng nucleotide trực tiếp dẫn đến biến đổi cấu trúc bậc một của protein kháng bệnh. Việc mất 3 nucleotide ở vị trí 543-545 đã làm tiêu biến 2 gốc amino acid quan trọng trong chuỗi polypeptide. Phân tích chi tiết thành phần 20 loại amino acid cho thấy Leucine chiếm tỷ lệ cao nhất với số lượng từ 140 đến 144 gốc amino acid (khoảng 14,3%), trong khi Tryptophan chiếm tỷ lệ thấp nhất với chỉ 15 đến 16 gốc (khoảng 1,6%). Cấu trúc giàu Leucine này hoàn toàn phù hợp với mô hình thụ thể LRR chuyên trách nhận diện nấm bệnh.

Dữ liệu nghiên cứu có thể được trực quan hóa thông qua biểu đồ phân bố tần suất 20 loại amino acid nhằm làm rõ sự khác biệt giữa 17 giống lúa bản địa và bảng ma trận đối sánh các đột biến InDel 17 bp và 24 bp. Đồng thời, ảnh chụp bản điện di gel polyacrylamide 6,0% thể hiện rõ các vạch phân đoạn ADN sắc nét, phân tách chính xác các alen kháng và nhiễm. Kết quả tầm soát gen Pit đạt tỷ lệ tương đồng 52,94% cao vượt trội so với nghiên cứu tầm soát gen SRWD2 trước đó trên 36 giống lúa Việt Nam (chỉ đạt 11,11%), chứng minh tập đoàn lúa bản địa là kho dự trữ gen kháng Pit vô cùng phong phú cho công tác lai tạo giống.

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

Để chuyển hóa các kết quả nghiên cứu thành giá trị thực tiễn phục vụ nền nông nghiệp, 4 giải pháp trọng tâm được đề xuất:

Ứng dụng ngay hệ thống chỉ thị phân tử Pit và Pi21 vào chương trình chọn giống lúa kháng đạo ôn bằng phương pháp MAS, hướng tới mục tiêu nâng cao độ chính xác chọn lọc kiểu gen lên 98% và giảm 50% chi phí khảo nghiệm đồng ruộng trong giai đoạn 2024 - 2026, do Bộ môn Kỹ thuật di truyền (Viện Di truyền Nông nghiệp) chủ trì phối hợp cùng các viện nghiên cứu cây lương thực.

Mở rộng quy mô tầm soát và giải mã hệ gen trên ít nhất 150 giống lúa bản địa khác tại các vùng sinh thái Tây Bắc và Tây Nguyên, nhằm xây dựng ngân hàng dữ liệu chứa 500 chỉ thị SNP và InDel liên kết gen kháng trong giai đoạn 2024 - 2027, do Viện Khoa học Nông nghiệp Việt Nam triển khai thực hiện.

Thực hiện chuyển nạp đồng thời cả 2 gen kháng Pit và Pi21 vào các giống lúa thuần chất lượng cao nhưng mẫn cảm với bệnh thông qua quy trình lai hồi giao MABC, đặt mục tiêu phục hồi trên 96% nền di truyền bố mẹ sau 4 thế hệ lai trong lộ trình 3 đến 5 năm, thực hiện bởi Trung tâm Khảo kiểm nghiệm Giống cây trồng Quốc gia.

Chuẩn hóa và tối ưu hóa quy trình phản ứng PCR Multiplex kết hợp đồng thời 2 cặp mồi Pit và Pi21 trong cùng một mẻ phản ứng 15 µl, giúp tiết kiệm 40% chi phí hóa chất tiêu hao và rút ngắn thời gian phân tích mẫu xuống dưới 4 giờ làm việc, áp dụng rộng rãi tại các phòng thí nghiệm công nghệ sinh học nông nghiệp từ năm 2024.

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

Tài liệu luận văn là nguồn tư liệu khoa học giá trị cao cho 4 nhóm đối tượng trọng tâm:

Các nhà chọn tạo giống cây trồng tại các viện và doanh nghiệp giống có thể nắm bắt chi tiết trình tự cặp mồi PitF/PitR và Pi21F/Pi21R cùng thông số nhiệt độ gắn mồi 55-58°C để tích hợp trực tiếp vào quy trình chọn giống phân tử MAS, rút ngắn 50% chu kỳ tạo giống kháng bệnh.

Giảng viên, học viên cao học và sinh viên ngành Công nghệ sinh học có thể sử dụng toàn bộ quy trình 8 bước tầm soát tin sinh học trên NextGENe phiên bản 2.3 và MEGA phiên bản 6.0 như một tài liệu chuẩn mực trong giảng dạy và thực hành phân tích dữ liệu giải trình tự NGS.

Kỹ thuật viên phòng thí nghiệm phân tử có thể tham khảo phác đồ tách chiết ADN bằng CTAB 1% và công thức phản ứng PCR 15 µl tối ưu (chứa 0,25 đơn vị Taq polymerase, 2,5 mM MgCl2) nhằm nâng cao tỷ lệ thành công của các phản ứng điện di lên 100%.

Cán bộ quản lý nguồn gen và bảo tồn tài nguyên thực vật có thể khai thác dữ liệu phân loại di truyền của 17 giống lúa bản địa (như Tám xoan Hải Hậu, Tép Thái Bình) để xây dựng chiến lược bảo tồn tại chỗ và khai thác bền vững nguồn gen quý phục vụ an ninh lương thực quốc gia.

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

Câu hỏi 1: Tại sao việc kết hợp cả hai gen Pit và Pi21 lại tạo ra tính kháng đạo ôn bền vững hơn? Trả lời: Gen Pit là gen kháng chính đơn diện có khả năng kích hoạt phản ứng miễn dịch siêu mẫn cực mạnh ngăn chặn nấm xâm nhập, trong khi Pi21 là gen lặn kháng phụ giúp hạn chế sự phát triển của sợi nấm trong mô bào. Khi tích hợp cả hai gen, giống lúa vừa có hàng rào phòng thủ ban đầu vừa có khả năng duy trì tính kháng phổ rộng, ngăn ngừa nguy cơ bị bẻ gãy tính kháng khi nấm nòi mới đột biến, với hiệu quả bảo vệ cây trồng đạt trên 95% trong thực tế canh tác.

Câu hỏi 2: Đột biến InDel 17 bp ở gen Pit có ý nghĩa gì đối với việc thiết kế chỉ thị phân tử? Trả lời: Đột biến mất đoạn 17 nucleotide ở vị trí từ 2647 đến 2664 phân chia rõ rệt 17 giống lúa thành hai nhóm kiểu gen (9 giống giữ nguyên và 8 giống mất đoạn). Độ chênh lệch 17 bp này tạo điều kiện thuận lợi để thiết kế cặp mồi nhân bản đoạn ADN mang đoạn chèn, giúp nhận diện chính xác alen kháng thông qua phương pháp điện di gel polyacrylamide 6,0% mà không cần phải giải trình tự tốn kém.

Câu hỏi 3: Vì sao luận văn sử dụng kết hợp phần mềm NextGENe phiên bản 2.3 và MEGA phiên bản 6.0 trong tầm soát gen? Trả lời: NextGENe phiên bản 2.3 chuyên trách xử lý các file dữ liệu đọc trình tự dạng FastQ từ hệ thống Illumina với độ chính xác trên 99,9%, lọc bỏ nhiễu và gọi đúng các điểm SNP và InDel. Tiếp đó, MEGA phiên bản 6.0 cho phép dóng hàng đa chuỗi, phân tích phả hệ tiến hóa và mô phỏng dịch mã sang 987 amino acid, giúp kiểm chứng toàn diện chức năng của locus gen trên 17 giống lúa.

Câu hỏi 4: Nồng độ hóa chất trong phản ứng PCR 15 µl được tối ưu hóa như thế nào? Trả lời: Phản ứng PCR được tối ưu hóa với 5 µl ADN khuôn nồng độ 10 ng, 0,15 µM mỗi mồi đặc hiệu, 0,2 mM dNTPs, 2,5 mM MgCl2 và 0,25 đơn vị Taq polymerase trong đệm 1X Buffer. Thành phần này đảm bảo độ khuếch đại tối đa của các đoạn gen mục tiêu trên máy luân nhiệt Veriti 96 giếng và không tạo ra bất kỳ băng phụ không đặc hiệu nào trên bản gel.

Câu hỏi 5: Các giống lúa bản địa mang gen Pit trong nghiên cứu có thể sử dụng trực tiếp làm vật liệu bố mẹ không? Trả lời: 9 giống lúa bản địa mang gen Pit nguyên vẹn như Thơm Lài, OM6377 và Tám xoan Hải Hậu là những nguồn vật liệu di truyền rất quý giá. Các nhà chọn giống có thể sử dụng ngay các giống này làm dòng cho gen trong quy trình lai hồi giao MABC để chuyển gen kháng vào các giống lúa cao sản thuần nông nghiệp trong thời gian từ 3 đến 4 vụ lai.

Kết luận

  • Tầm soát thành công 100% trình tự tương đồng của hai gen kháng đạo ôn chủ lực Pit và Pi21 trên toàn bộ 17 giống lúa bản địa Việt Nam.
  • Xác định chính xác 9/17 giống lúa bản địa (đạt tỷ lệ 52,94%) mang cấu trúc gen Pit chức năng mã hóa chuỗi protein gồm 987 amino acid.
  • Phát hiện các đột biến InDel đặc trưng gồm đoạn mất 17 bp trên gen Pit và 24 bp trên gen Pi21 làm cơ sở phân loại kiểu gen chính xác.
  • Thiết kế và chuẩn hóa thành công 2 cặp chỉ thị ADN đặc hiệu (PitF/PitR và Pi21F/Pi21R) cho độ khuếch đại chuẩn xác 100% trên điện di polyacrylamide 6,0%.
  • Cung cấp cơ sở dữ liệu hệ gen và công cụ phân tử hoàn chỉnh giúp rút ngắn 50% thời gian chọn tạo giống lúa kháng đạo ôn.

Đóng góp lớn nhất của luận văn là đã chuyển hóa dữ liệu giải mã hệ gen phức tạp thành các công cụ chỉ thị ADN thực nghiệm đơn giản, chi phí thấp nhưng có độ tin cậy tuyệt đối. Lộ trình triển khai tiếp theo trong giai đoạn 2024 - 2026 sẽ tập trung vào việc lai chuyển các gen này vào các giống lúa xuất khẩu chủ lực. Các nhà khoa học, chuyên gia chọn giống và đơn vị sản xuất nông nghiệp hãy khai thác toàn văn tài liệu và ứng dụng ngay hệ thống chỉ thị phân tử này vào các chương trình lai tạo giống lúa kháng bệnh công nghệ cao!