Tổng quan nghiên cứu

Nhiễm mặn đất canh tác đang là một trong những hiểm họa phi sinh học nghiêm trọng nhất đối với an ninh lương thực toàn cầu. Theo các báo cáo nông nghiệp quốc tế, thế giới đối mặt với áp lực phải gia tăng khoảng 70% sản lượng lương thực để cung cấp cho thêm 2,3 tỷ người vào năm 2050. Tuy nhiên, biến đổi khí hậu và nước biển dâng có thể khiến 50% diện tích đất nông nghiệp hiện tại bị nhiễm mặn vào năm 2050. Việt Nam sở hữu đường bờ biển dài 3.260 km cùng hệ thống sông ngòi dày đặc, khiến hai vựa lúa trọng điểm là đồng bằng sông Hồng và đồng bằng sông Cửu Long thường xuyên chịu ảnh hưởng nặng nề bởi hiện tượng xâm nhập mặn.

Về mặt sinh lý, đất được xem là nhiễm mặn khi độ dẫn điện của dung dịch đất vượt quá 4 dS/m (tương đương 40 mM NaCl). Khi độ mặn trung bình vượt ngưỡng 1,9 dS/m, năng suất hạt lúa bắt đầu suy giảm rõ rệt; đặc biệt ở độ mặn trên 3,0 dS/m, mật độ cây lúa có thể sụt giảm từ 30% đến 40%. Stress mặn gây mất cân bằng ion do tích lũy quá mức ion Na+ và Cl-, ức chế nghiêm trọng quá trình quang hợp và sinh trưởng. Trong khi lúa gạo nuôi sống hơn 90% dân số châu Á, việc phát triển các giống lúa chịu mặn mang tính cấp thiết và tiết kiệm chi phí hơn nhiều so với các giải pháp công trình cải tạo đất.

Luận văn thạc sĩ chuyên ngành Di truyền học của tác giả Đặng Ngọc Hoa, thực hiện dưới sự hướng dẫn của Tiến sĩ Đỗ Thị Phúc tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, tập trung giải quyết bài toán cốt lõi: Đánh giá tính đa hình trình tự nucleotide vùng promoter, vùng mã hóa và phân tích mức độ biểu hiện của gen OsHKT1;4 ở các giống lúa dưới các mức độ stress mặn khác nhau. Nghiên cứu cung cấp cơ sở dữ liệu di truyền học phân tử chuẩn xác, mở đường cho việc ứng dụng chỉ thị phân tử trong chọn tạo giống lúa kháng mặn tại Việt Nam.

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 các khung lý thuyết sinh học phân tử và di truyền học thực vật hiện đại:

Cơ chế cân bằng nội môi ion và thải loại Na+: Khả năng chống chịu mặn ở cây lúa (Oryza sativa) phụ thuộc chặt chẽ vào hoạt động của các hệ thống vận chuyển ion qua màng tế bào. Hệ thống này bao gồm kênh đối chuyển Na+/H+ trên màng sinh chất (SOS1) giúp đẩy Na+ ra khỏi tế bào, kênh trao đổi Na+/H+ trên màng không bào (NHX1) giúp cô lập Na+ vào không bào, và họ protein màng vận chuyển cation HKT (High-Affinity K+ Transporter).

Mô hình họ gen HKT ở lúa: Họ gen HKT ở lúa gồm 9 thành viên, chia thành 2 phân họ. Phân họ 1 gồm 5 gen (OsHKT1;1 đến OsHKT1;5) có tính chọn lọc cao với Na+. Trong đó, gen OsHKT1;4 nằm trên nhiễm sắc thể số 4, có kích thước 5.252 bp với khung đọc mở (ORF) dài 1.503 bp mã hóa chuỗi polypeptide gồm 500 axit amin. Protein OsHKT1;4 định vị chủ yếu tại mô xylem thân bẹ và lá non, đảm nhận chức năng thu hồi Na+ từ dòng mạch gỗ, ngăn chặn sự di chuyển và tích lũy ion độc hại này lên các mô quang hợp nhạy cảm.

Lý thuyết điều hòa biểu hiện gen qua yếu tố cis-acting: Mức độ phiên mã của gen đáp ứng stress phi sinh học chịu sự kiểm soát của các trình tự điều hòa trên vùng promoter, nơi các yếu tố phiên mã (như GT-1, WRKY, MYB, GATA) gắn vào để kích hoạt hoặc ức chế biểu hiện gen dưới tác động của tín hiệu môi trường.

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

Nguồn dữ liệu và vật liệu chọn mẫu: Nghiên cứu sử dụng cỡ mẫu gồm 12 giống lúa được tuyển chọn từ Viện Di truyền Nông nghiệp và Học viện Nông nghiệp Việt Nam. Mẫu nghiên cứu bao gồm 2 giống chuẩn quốc tế (giống mẫn cảm Nipponbare và giống kháng mặn Pokkali) cùng 10 giống lúa bản địa Việt Nam (Nước mặn dạng 1, Chiêm rong, Cườm dạng 1, Nước mặn dạng 2, Cườm dạng 2, Chăm biển, Nếp cúc, Chăm, Tép lai, Chiêm đen). Phương pháp chọn mẫu có chủ đích giúp đối chiếu đa dạng di truyền giữa các nguồn gen bản địa thích nghi cục bộ với giống chuẩn quốc tế.

Quy trình phân tích đa hình gen: DNA tổng số được tách chiết từ mô lá non 50 mg bằng phương pháp CTAB cải tiến, kiểm tra độ tinh sạch qua đo tỷ lệ quang phổ A260/A280 và điện di trên gel agarose 1%. Phản ứng PCR khuếch đại vùng promoter (khoảng 1,8 kb) cùng vùng exon 2 (~700 bp) và exon 3 (~650 bp) sử dụng các cặp mồi đặc hiệu thiết kế qua công cụ Primer-BLAST. Sản phẩm PCR tinh sạch qua cột GeneJET được giải trình tự Sanger tự động trên hệ thống ABI PRISM 3730xl. Trình tự được gióng hàng bằng phần mềm MUTALIN và phân tích yếu tố cis qua cơ sở dữ liệu PLACE/Signal Scan.

Quy trình phân tích biểu hiện gen bán định lượng: Hai giống đại diện (Nipponbare và Pokkali) được trồng thủy canh trong dung dịch dinh dưỡng Yoshida (pH 5,5) suốt 14 ngày trước khi xử lý sốc mặn ở 2 mức nồng độ 50 mM và 100 mM NaCl. Mô lá được thu thập tại 3 thời điểm: 24 giờ, 48 giờ và 72 giờ. RNA tổng số được tách chiết bằng kit GeneJET Plant RNA, xử lý loại bỏ DNA bằng enzym DNase I, và tổng hợp cDNA chuỗi một qua enzym RevertAid M-MuLV. Phản ứng RT-PCR sử dụng mồi đặc hiệu RT-HKT1;4 (vắt qua ranh giới exon 1 và exon 2) được tiến hành đồng thời với gen nội chuẩn Actin (ACT). Mật độ quang của các băng điện di trên gel agarose 2% được số hóa và định lượng bán định lượng thông qua phần mềm ImageJ.

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

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

Phân tích đa hình vùng promoter: Quá trình giải trình tự đoạn promoter kích thước khoảng 1,8 kb trên 12 giống lúa và so sánh với cơ sở dữ liệu Rice Genome Annotation Project đã xác định được 15 vị trí đa hình. Trong đó có 14 vị trí đa hình đơn nucleotide (SNP) và 1 vị trí đột biến mất/thêm đoạn (indel) dài 12 bp (trình tự 5’-TCAAATGATTCA-3’ tại vị trí -1553). Dựa trên sự tương đồng trình tự, 12 giống lúa được phân chia thành 2 nhóm rõ rệt: Nhóm I gồm 3 giống (Nipponbare, Nếp cúc, Cườm dạng 2) có trình tự promoter tương đồng 100% với dữ liệu chuẩn; Nhóm II gồm 9 giống (Pokkali, Nước mặn dạng 1, Chiêm rong, Chiêm đen, Tép lai, Chăm, Cườm dạng 1, Chăm biển, Nước mặn dạng 2) mang trọn vẹn 15 sai khác di truyền.

Đặc điểm các yếu tố điều hòa cis: Phân tích tin sinh học phát hiện 17 vị trí yếu tố cis GT-1 (trình tự GRWAAW) trên vùng promoter. Yếu tố này chuyên biệt tham gia đáp ứng với stress phi sinh học và tín hiệu ánh sáng. Tại Nhóm II, đột biến thay thế C/G ở vị trí -1695 tạo thêm 1 vị trí GT-1, trong khi đột biến mất đoạn 12 bp ở vị trí -1553 làm mất đi 1 vị trí GT-1. Do đó, tổng số lượng vị trí nhận biết GT-1 vẫn được duy trì bảo tồn ở mức 17 vị trí trên cả hai nhóm giống.

Bảo tồn tuyệt đối vùng mã hóa exon 2 và exon 3: Kết quả khuếch đại và giải trình tự hai phân đoạn exon 2 (~700 bp) và exon 3 (~650 bp) cho thấy độ tương đồng đạt 100% giữa tất cả 12 giống lúa nghiên cứu và cơ sở dữ liệu tham chiếu quốc tế, không xuất hiện bất kỳ đột biến thay thế hay dịch khung nào.

Khác biệt về đáp ứng hình thái và mức độ biểu hiện phiên mã: Dưới nồng độ 100 mM NaCl, giống mẫn cảm Nipponbare xuất hiện hiện tượng xoăn và cháy đầu lá già rõ rệt chỉ sau 24 giờ, mức độ tổn thương tăng dần tại 48 giờ và 72 giờ. Ngược lại, giống kháng mặn Pokkali sinh trưởng bình thường, không ghi nhận tổn thương hình thái. Kết quả RT-PCR cho thấy mức độ tích lũy bản phiên mã mRNA của gen OsHKT1;4 ở mô lá giống Pokkali tăng mạnh và duy trì ổn định dưới xử lý mặn 50 mM và 100 mM NaCl tại các mốc 24h, 48h và 72h, trong khi giống Nipponbare biểu hiện yếu và giảm dần theo thời gian chịu stress.

Thảo luận kết quả

Sự bảo tồn 100% trình tự nucleotide tại exon 2 và exon 3 chỉ ra rằng cấu trúc protein vận chuyển OsHKT1;4 chịu áp lực chọn lọc tiến hóa cao nhằm duy trì tính toàn vẹn của lỗ kênh và các vòng xoắn xuyên màng vận chuyển Na+. Sự phân hóa khả năng chịu mặn giữa các giống lúa không xuất phát từ biến đổi cấu trúc protein mà chủ yếu được điều khiển bởi cơ chế điều hòa biểu hiện gen ở mức độ phiên mã thông qua vùng promoter.

Các dữ liệu nghiên cứu có thể được trực quan hóa hiệu quả thông qua hệ thống bảng ma trận đa hình nucleotide và biểu đồ cột biểu diễn cường độ biểu hiện gen tương đối từ ảnh điện di qua phân tích ImageJ. Sự tương đồng di truyền giữa 8 giống lúa bản địa Việt Nam (như Chiêm rong, Chăm biển, Nước mặn dạng 1) với giống chuẩn kháng mặn quốc tế Pokkali tại 15 vị trí đa hình promoter khẳng định nguồn gen bản địa của nước ta mang tiềm năng chống chịu mặn rất cao. Sự hiện diện dày đặc của 17 yếu tố cis GT-1 kết hợp cùng các motif ARFAT, WRKY71, GATABOX tạo nên mạng lưới điều hòa phiên mã linh hoạt, giúp kích hoạt nhanh chóng gen OsHKT1;4 ngay khi nồng độ ion Na+ ngoại bào tăng cao.

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

Ứng dụng 15 chỉ thị SNP/Indel trên promoter vào chọn giống nhờ chỉ thị phân tử (MAS): Viện Di truyền Nông nghiệp và các trung tâm giống cây trồng cần chuyển đổi 15 vị trí đa hình đã phát hiện thành các marker phân tử chuyên dụng (như CAPS hoặc dCAPS), giúp sàng lọc nhanh nguồn gen chịu mặn ở giai đoạn cây mạ, rút ngắn khoảng 30% đến 50% thời gian lai tạo giống lúa mới trong giai đoạn 2026-2030.

Nâng cấp phương pháp phân tích biểu hiện bằng Real-time RT-PCR: Các nhóm nghiên cứu sinh học phân tử cần mở rộng đánh giá định lượng tuyệt đối (Absolute Quantification Real-time PCR) mức độ biểu hiện gen OsHKT1;4 trên toàn bộ 10 giống lúa bản địa thuộc Nhóm II, khảo sát dải nồng độ mặn mở rộng từ 50 mM đến 150 mM NaCl trong vòng 12 đến 18 tháng tới.

Xác thực chức năng kênh vận chuyển qua công nghệ chỉnh sửa gen và chuyển gen: Các viện nghiên cứu phối hợp với trường đại học triển khai tạo dòng lúa mang đột biến định hướng (knock-out hoặc overexpression) đối với gen OsHKT1;4 trên nền giống lúa thuần chất lượng cao, nhằm đo lường trực tiếp tỷ lệ K+/Na+ nội bào và chỉ số tích lũy sinh khối trong điều kiện stress mặn thực tế trong chu kỳ 2 đến 3 năm.

Khảo nghiệm đồng ruộng và khai thác nguồn gen lúa bản địa tại các vùng ven biển: Sở Nông nghiệp và Phát triển nông thôn tại các tỉnh đồng bằng sông Cửu Long và ven biển miền Bắc cần phối hợp cùng các nhà khoa học đưa các giống lúa địa phương ưu tú thuộc Nhóm II (như giống Chăm biển, Nước mặn dạng 1, Chiêm rong) vào khảo nghiệm sinh thái tại các vùng đất có độ mặn từ 4 đến 6 dS/m trong 3 vụ canh tác liên tiếp để phục tráng và nhân rộng.

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

Học viên cao học và nghiên cứu sinh chuyên ngành Di truyền học, Sinh học phân tử: Tài liệu cung cấp quy trình thực nghiệm chuẩn mực từ thiết kế mồi RT-PCR gối ranh giới exon, xử lý enzym DNase I tinh sạch RNA, giải trình tự Sanger cho đến phân tích yếu tố cis bằng công cụ tin sinh học PLACE.

Các nhà chọn tạo giống cây trồng tại các viện nghiên cứu nông nghiệp: Cung cấp nguồn dữ liệu phân tử xác thực về 10 nguồn gen lúa bản địa Việt Nam, giúp định hướng chiến lược khai thác alen có lợi phục vụ các dự án tạo giống lúa chống chịu biến đổi khí hậu.

Giảng viên và cán bộ nghiên cứu tại các trường đại học khối Nông - Lâm - Sinh: Luận văn là tài liệu tham khảo giá trị cho các bài giảng chuyên đề về sinh lý chống chịu stress phi sinh học ở thực vật, di truyền học số lượng và cơ chế vận chuyển ion qua màng tế bào.

Cán bộ quản lý khuyến nông và hoạch định chính sách nông nghiệp địa phương: Giúp nắm bắt tiềm năng di truyền của các giống lúa địa phương để xây dựng kế hoạch bảo tồn quỹ gen thực vật quý hiếm và cơ cấu mùa vụ hợp lý trên hơn 1 triệu hecta đất nhiễm mặn ven biển.

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

Gen OsHKT1;4 có vai trò sinh học cụ thể như thế nào đối với khả năng chịu mặn ở lúa? Gen OsHKT1;4 mã hóa protein màng dài 500 axit amin, thuộc phân họ 1 của họ HKT. Protein này định vị tại thân bẹ và mạch gỗ, có chức năng thu hồi ion Na+ ra khỏi dòng xylem, ngăn chặn Na+ vận chuyển lên lá non và cơ quan quang hợp, từ đó duy trì tỷ lệ K+/Na+ tối ưu giúp cây sống sót ở nồng độ mặn trên 40 mM NaCl.

Tại sao nghiên cứu ghi nhận sự đa hình ở promoter nhưng exon 2 và exon 3 lại bảo tồn 100%? Vùng exon 2 và exon 3 mã hóa các domain chức năng xuyên màng cốt lõi của protein OsHKT1;4 nên được chọn lọc tự nhiên bảo tồn tuyệt đối 100% để duy trì hoạt tính vận chuyển ion. Ngược lại, vùng promoter xuất hiện 15 vị trí đa hình nhằm tạo ra các mức độ đáp ứng phiên mã linh hoạt khác nhau giữa các giống lúa trước áp lực sinh thái.

Nhóm I và Nhóm II trong nghiên cứu khác nhau cơ bản ở những điểm nào? Kết quả giải trình tự đoạn promoter dài khoảng 1,8 kb đã chia 12 giống thành 2 nhóm: Nhóm I gồm 3 giống (Nipponbare, Nếp cúc, Cườm dạng 2) giống hệt trình tự chuẩn; Nhóm II gồm 9 giống (chứa Pokkali và 8 giống địa phương) mang 14 điểm đột biến SNP và 1 đột biến mất đoạn 12 bp liên quan đến tính kháng mặn.

Yếu tố cis GT-1 trên vùng promoter có vai trò gì trong phản ứng chống chịu stress mặn? Yếu tố cis GT-1 là vị trí gắn đặc hiệu của các yếu tố phiên mã họ trihelix. Trên vùng promoter gen OsHKT1;4 có tới 17 vị trí GT-1, đóng vai trò công tắc sinh học nhận diện tín hiệu stress phi sinh học và cảm ứng gia tăng phiên mã khi đất bị nhiễm mặn vượt quá 4 dS/m.

Tại sao cần thiết kế mồi RT-PCR vắt qua ranh giới giữa exon 1 và exon 2? Việc thiết kế mồi ngược gối qua điểm nối exon 1 và exon 2 đảm bảo phản ứng RT-PCR chỉ khuếch đại phân tử cDNA đã hoàn tất quá trình cắt nối intron (splicing), loại bỏ hoàn toàn khả năng nhân bản nhầm lẫn từ DNA hệ gen (gDNA) tạp nhiễm, đảm bảo độ chính xác bán định lượng tuyệt đối khi so sánh với gen chuẩn Actin.

Kết luận

Xác lập thành công bản đồ đa hình di truyền của gen OsHKT1;4 trên 12 giống lúa, phát hiện 15 vị trí đa hình trên đoạn promoter dài khoảng 1,8 kb và chứng minh tính bảo tồn cấu trúc 100% tại vùng mã hóa exon 2 và exon 3.

Nhận diện 17 yếu tố cis GT-1 điều hòa đáp ứng stress phi sinh học trên vùng promoter, làm sáng tỏ cơ sở sinh học phân tử chi phối khả năng cảm ứng mặn ở lúa.

Chứng minh sự vượt trội về mức độ biểu hiện mRNA của gen OsHKT1;4 ở mô lá giống kháng mặn Pokkali so với giống mẫn cảm Nipponbare dưới các điều kiện stress 50 mM và 100 mM NaCl tại các thời điểm 24h, 48h và 72h.

Khẳng định 8 giống lúa bản địa Việt Nam thuộc Nhóm II mang cấu trúc promoter tương đồng với giống kháng Pokkali, mở ra nguồn vật liệu di truyền vô giá cho các chương trình chọn giống giai đoạn 2026-2030.

Cung cấp hệ thống cơ sở lý luận và chỉ thị phân tử tin cậy, đóng góp thiết thực cho chiến lược phát triển nền nông nghiệp lúa nước bền vững và thích ứng hiệu quả với biến đổi khí hậu toàn cầu.