Tổng quan luận án

Luận án tiến sĩ "Xác định các gen - alen đặc thù liên quan đến sự phát triển bộ rễ của các giống lúa Việt Nam" được thực hiện bởi nghiên cứu sinh Phùng Thị Phương Nhung, thuộc chuyên ngành Di truyền và Chọn giống cây trồng (mã số: 9 62 01 11) tại Học viện Nông nghiệp Việt Nam năm 2019, dưới sự hướng dẫn khoa học của GS. Đỗ Năng Vịnh và GS. Pascal Gantet.

Tính cấp thiết và khoảng trống nghiên cứu

Cây lúa (Oryza sativa L.) là cây lương thực thiết yếu nuôi sống hơn một nửa dân số thế giới với diện tích canh tác toàn cầu dao động từ 156 đến 161,7 triệu hecta. Trước các tác động của biến đổi khí hậu như hạn hán, xâm nhập mặn và ngập lụt kéo dài, bộ rễ giữ vai trò quyết định đến khả năng bám giữ, hấp thu nước, vận chuyển dinh dưỡng và duy trì năng suất cây trồng. Một bộ rễ ăn sâu, rễ dày và phân nhánh tốt giúp cây lúa khai thác nước ở các tầng đất sâu trong điều kiện khô hạn, trong khi các cơ chế giải phẫu như mô khí giúp rễ thích ứng với điều kiện ngập nước.

Tuy nhiên, việc chọn tạo giống lúa có kiến trúc rễ tối ưu gặp nhiều trở ngại do rễ nằm dưới mặt đất, khó quan sát trực tiếp. Các nghiên cứu di truyền truyền thống xác định vị trí tính trạng số lượng (QTLs) dựa trên các quần thể lai lưỡng thân (như RILs, DHLs) bộc lộ hạn chế lớn: khoảng tin cậy của đoạn nhiễm sắc thể mang QTL rất rộng (tính bằng mega-base), số lượng alen đánh giá bị giới hạn giữa hai bố mẹ, và thời gian tạo quần thể kéo dài. Phương pháp nghiên cứu lập bản đồ liên kết toàn hệ gen (GWAS) kết hợp công nghệ giải trình tự thế hệ mới (NGS) và phân tích kiểu gen qua giải trình tự (GBS) cho phép khắc phục nhược điểm này nhờ mật độ chỉ thị SNP cao và sự phân rã nhanh của hiện tượng mất cân bằng liên kết (LD), giúp thu hẹp độ phân giải QTL xuống mức kilo-base. Dù nguồn gen lúa bản địa Việt Nam rất phong phú và phân bố trên nhiều vùng sinh thái đặc thù, nguồn tài nguyên này chưa được khai thác toàn diện ở mức độ hệ gen để phục vụ xác định các QTLs và gen ứng viên kiểm soát sự phát triển bộ rễ.

Mục tiêu nghiên cứu

Mục tiêu nghiên cứu được tác giả xác định gồm 4 nội dung:

  1. Lựa chọn các mẫu giống phù hợp để phát triển bộ dữ liệu kiểu gen và dữ liệu kiểu hình phục vụ cho các nghiên cứu GWAS thông qua khảo sát sự đa dạng về đặc điểm nông sinh học cơ bản và di truyền của một tập đoàn các mẫu giống lúa Việt Nam.
  2. Xây dựng bộ dữ liệu kiểu gen (haplotype) của tập đoàn mẫu giống được chọn, làm cơ sở dữ liệu để phát triển các nghiên cứu GWAS với các tính trạng quan tâm.
  3. Xây dựng bộ dữ liệu kiểu hình của một số tính trạng chính liên quan đến sự phát triển bộ rễ của các mẫu giống lúa được chọn, làm cơ sở cho nghiên cứu GWAS để xác định các QTLs/gen ứng viên liên quan đến sự phát triển bộ rễ của lúa Việt Nam.
  4. Sử dụng phương pháp GWAS để lập bản đồ liên kết toàn hệ gen, từ đó xác định các QTLs và gen ứng viên liên quan đến sự phát triển bộ rễ của các giống lúa Việt Nam.

Câu hỏi và giả thuyết nghiên cứu

Nghiên cứu tập trung giải quyết vấn đề: Liệu trong tập đoàn giống lúa bản địa Việt Nam có tồn tại các biến dị di truyền, alen đặc thù và các locus tính trạng số lượng (QTLs) liên kết chặt chẽ với các chỉ tiêu hình thái, khối lượng và kiến trúc phát triển của bộ rễ hay không; và phương pháp GWAS dựa trên hệ thống chỉ thị SNP mật độ cao từ GBS có thể xác định chính xác các gen ứng viên chức năng liên quan đến sự phát triển bộ rễ hay không.

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

  • Đối tượng nghiên cứu: Tập đoàn nguồn gen lúa bao gồm 270 mẫu giống ban đầu (trong đó có 214 mẫu giống lúa địa phương thu thập từ nhiều vùng sinh thái của Việt Nam, 33 mẫu giống đối chứng đại diện cho tính đa dạng của Oryza sativa trên thế giới do CIRAD cung cấp, và 23 mẫu giống do Viện Di truyền Nông nghiệp cung cấp).
  • Phạm vi kiểu gen: 270 mẫu giống được phân tích đa dạng di truyền bằng 241 chỉ thị DArT. Sau đó, tập đoàn gồm 200 mẫu giống lúa đại diện được lựa chọn để phân tích kiểu gen bằng phương pháp GBS với hơn 50.000 chỉ thị SNP thô, tạo ma trận 25.971 chỉ thị SNP đa hình chất lượng cao.
  • Phạm vi kiểu hình: Đánh giá 18 tính trạng nông sinh học và đặc điểm bộ rễ trên 194 mẫu giống ở giai đoạn cây 6 tuần tuổi trong hệ thống ống rễ cải tiến đặt tại nhà lưới.
  • Phạm vi phân tích GWAS: Thực hiện trên 3 ma trận dữ liệu: toàn bộ tập đoàn (185 mẫu giống × 21.623 marker), phân nhóm phụ indica (115 mẫu giống × 13.842 marker), và phân nhóm phụ japonica (64 mẫu giống × 8.821 marker).

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước trong y văn

Tác giả đã hệ thống hóa các tài liệu khoa học liên quan đến cấu trúc, chức năng sinh lý và cơ sở di truyền phân tử của bộ rễ lúa:

  • Đặc điểm hình thái và giải phẫu rễ: Rebouillat et al. (2009) và Coudert et al. (2010) mô tả bộ rễ lúa gồm 5 loại chính: rễ mầm (radicle - ra), rễ bất định phôi (embryonic crown root - ecr), rễ bất định sau nảy mầm hình thành từ các đốt thân (crown root - cr), rễ bên lớn (large lateral root - llr) và rễ bên nhỏ (small lateral root - slr). Cấu trúc xuyên tâm gồm trung trụ (xylem, phloem, trụ bì, nội bì), lớp vỏ mô khí (aerenchyma), lớp cứng mô, ngoại bì và biểu bì. Vùng trung tâm tĩnh (Quiescent Center - QC) ở chóp rễ duy trì hoạt động phân chia mô phân sinh (Coudert et al., 2010). Yoshida và Hasegawa (1982) chỉ ra quy luật phát sinh rễ: khi lá thứ $n$ xuất hiện, rễ sẽ đơm ra từ đốt thứ $(n-3)$, mỗi đốt phát sinh từ 5 đến 25 rễ.
  • Các nghiên cứu QTLs rễ qua quần thể lập bản đồ truyền thống: Champoux et al. (1995) công bố nghiên cứu đầu tiên về kiến trúc rễ lúa. Kamoshita et al. (2008) nghiên cứu trên quần thể CT9993/IR622266; Khowaja et al. (2009) trên quần thể Bala/Azucena. Courtois et al. (2009) tổng hợp 675 QTLs từ 24 nghiên cứu trên 12 quần thể lập bản đồ (chủ yếu là RILs hoặc DHLs từ phép lai indica × japonica), chỉ ra các cụm QTLs tập trung trên nhiễm sắc thể 1 (30-40 Mb), NST 2 (25-35 Mb), NST 3 (0-5 Mb), NST 4 (30-35 Mb), NST 9 (15-20 Mb). Topp et al. (2013) áp dụng công nghệ dựng ảnh 3D để tăng độ chính xác đánh giá kiểu hình rễ.
  • Dòng đẳng gen (NILs) và phân lập gen từ QTLs: Gen PUP1 và PSTOL1 (PHOSPHORUS-STARVATION TOLERANCE 1) được Gamuyao et al. (2012) nhân dòng từ QTL, mã hóa thụ thể kinase tế bào chất giúp tăng trưởng rễ sớm và tăng 60% năng suất trong đất thiếu lân. Gen DRO1 (DEEPER ROOTING 1) được Uga et al. (2013, 2015) clone từ QTL kiểm soát tính hướng trọng lực của rễ thông qua tín hiệu auxin (ARF), làm tăng góc rễ đâm sâu.
  • Các gen và mạng lưới điều hòa phát triển rễ:
    • Gen cảm ứng auxin: CRL1/ARL1 mã hóa yếu tố phiên mã LBD (OsLBD3-2) kiểm soát phát sinh rễ bất định (Inukai et al., 2005; Liu et al., 2005); CRL4/OsGNOM1 điều khiển vận chuyển auxin qua màng liên kết với protein G và OsPIN1 (Kitomi et al., 2008; Liu et al., 2009); cặp gen OsWOX3A/NAL2/NAL3 (Cho et al., 2013); thụ thể OsTIR1 và OsAFB2 bị ức chế bởi OsMir393 (Bian et al., 2012; Xia et al., 2012); chaperone OsCYP2/LRL2 tương tác với OsIAA11 và OsIAA13 (Kang et al., 2013; Kitomi et al., 2012; Zhu et al., 2012).
    • Gen cảm ứng phytohormone khác: CRL5 kích hoạt các yếu tố ARRs ức chế cytokinin (Kitomi et al., 2011); WOX11 điều hòa rễ bất định thông qua tương tác auxin - cytokinin (Zhao et al., 2009); HOX1 điều hòa bởi auxin, jasmonic acid và NO (Chen et al., 2012; Hsu et al., 2013).
    • Mô phân sinh và trung tâm QC: OsIAA23 (Jun et al., 2011); QHB đồng đẳng với WOX5 (Kamiya et al., 2003a); thụ thể OsRPK1 (Zou et al., 2010); các yếu tố phiên mã GRAS gồm OsSCR1 và OsSHR1 (Ni et al., 2014); DOCS1 điều hòa biệt hóa ngoại bì và phản ứng kháng nhôm (Huang et al., 2009).
    • Rễ hậu phôi và kéo dài vách tế bào: OsCAND1 (Wang et al., 2011); OsORC3 (Chen et al., 2013); OsDGL1 (Qin et al., 2009); OsGLU3 (Zhang et al., 2012); OsEXPA8 (Ma et al., 2013); OsSPR1 (Jia et al., 2011); RSS1 và RSS3 duy trì phân bào rễ trong điều kiện mặn (Ogawa et al., 2011; Toda et al., 2013); các yếu tố phiên mã NAC như OsNAC5, OsNAC9, OsNAC10 giúp tăng đường kính rễ và chịu hạn (Jeong et al., 2010, 2013; Redillas et al., 2012).
  • Nguyên lý và ứng dụng của GBS và GWAS: Kỹ thuật NGS và GBS cắt enzyme giới hạn (ApeKI, PstI, SbfI) gắn adapter mã vạch (Elshire et al., 2011; Davey et al., 2011; Sonah et al., 2013). Các nghiên cứu GWAS về rễ lúa trên tập đoàn quốc tế của Famoso et al. (2011), Clark et al. (2013), Courtois et al. (2013).

Khoảng trống nghiên cứu luận án giải quyết

Mặc dù các nghiên cứu quốc tế đã xây dựng một số bộ dữ liệu GWAS, các tập đoàn đó chỉ đại diện cho một phần nhỏ tính đa dạng di truyền toàn cầu. Nguồn gen lúa địa phương Việt Nam phân bố trên nhiều tiểu vùng sinh thái (tưới tiêu, nước trời, nương rẫy vùng cao, ngập úng, mặn ven biển) chứa đựng các alen thích nghi đặc thù nhưng chưa từng được giải trình tự quy mô lớn để lập bản đồ haplotype SNP mật độ cao và thực hiện GWAS đối với toàn bộ các chỉ tiêu hình thái, khối lượng của bộ rễ. Luận án giải quyết khoảng trống này bằng việc xây dựng bộ dữ liệu kiểu gen GBS và kiểu hình chuẩn hóa để lập bản đồ liên kết toàn hệ gen cho tập đoàn lúa Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và mô hình phân tích

  • Cơ sở di truyền số lượng và mất cân bằng liên kết (Linkage Disequilibrium - LD): Đánh giá mức độ liên kết không ngẫu nhiên giữa các alen tại các locus khác nhau trên 12 nhiễm sắc thể, xác định tốc độ suy giảm LD theo khoảng cách vật lý (kb/Mb) để xác định kích thước vùng tin cậy chứa gen ứng viên.
  • Mô hình lập bản đồ liên kết toàn hệ gen (GWAS): Sử dụng Mô hình tuyến tính hỗn hợp (Mixed Linear Model - MLM) nhằm kiểm soát đồng thời cấu trúc quần thể (ma trận Q thu được từ phân tích STRUCTURE hoặc PCA) và quan hệ họ hàng di truyền thân thuộc giữa các dòng (ma trận Kinship - K), hạn chế tối đa sai số dương tính giả.

Phương pháp thực nghiệm và thu thập dữ liệu

  • Đánh giá đặc điểm nông sinh học cơ bản: Khảo sát 270 mẫu giống theo tiêu chuẩn của IRRI (2002), ghi nhận các chỉ tiêu: thời gian sinh trưởng, chiều cao cây, số nhánh, số nhánh hữu hiệu, hình dạng hạt (tỷ lệ dài/rộng), tính chất nội nhũ (gạo tẻ, gạo nếp).
  • Chiết tách ADN và phân tích đa dạng di truyền bằng chỉ thị DArT: Chiết tách ADN tổng số từ lá non, lai với mảng vi mảng DArT gồm 6.144 clone, thu nhận 241 chỉ thị DArT đa hình rõ nét để phân tích cấu trúc di truyền lưỡng cực và xây dựng cây phân loại.
  • Phân tích kiểu gen qua giải trình tự (GBS): Chiết tách ADN của 200 mẫu giống lúa chọn lọc, xử lý bằng enzyme cắt giới hạn, gắn adapter mã vạch (barcode dài 4-8 nucleotide), giải trình tự trên hệ thống Illumina thế hệ mới. Dữ liệu thô được xử lý qua pipeline lọc chất lượng, căn chỉnh với hệ gen đối chứng, loại bỏ các marker có tần số alen nhỏ (MAF < 5%) và nội suy các điểm dữ liệu khuyết bằng thuật toán BEAGLE, tạo thành ma trận 25.971 chỉ thị SNP.
  • Đánh giá kiểu hình bộ rễ: Thí nghiệm tiến hành trong nhà lưới với 194 mẫu giống lúa, bố trí theo khối ngẫu nhiên hoàn chỉnh với 3 lần lặp lại. Cây lúa được trồng trong hệ thống ống rễ cải tiến chuyên dụng. Thu thập mẫu rễ ở giai đoạn cây 6 tuần tuổi, phân tích hình ảnh và đo đếm 18 tính trạng bằng phần mềm RASTA kết hợp cân sinh khối khô.
  • Phương pháp xác định gen ứng viên: Dựa trên các đỉnh SNP có ý nghĩa thống kê vượt ngưỡng ($P < 10^{-4}$) trên đồ thị Manhattan Plot, khoanh vùng khoảng cách vật lý theo bán kính suy giảm LD xung quanh QTL, tra cứu cơ sở dữ liệu hệ gen lúa để lọc các gen mã hóa và đối chiếu với 24 gen chức năng đã được chứng minh trong y văn.

Hệ thống các tính trạng kiểu hình theo dõi

  1. NCR (Number of Crown Root): Số lượng rễ bất định.
  2. THK (Root Thickness): Độ dày (đường kính) của rễ.
  3. MRL (Maximum Root Length): Chiều dài rễ tối đa (cm).
  4. DEPTH: Độ ăn sâu của rễ trong ống nghiệm.
  5. RDW (Root Dry Weight): Tổng khối lượng khô của bộ rễ (g).
  6. SDW (Shoot Dry Weight): Khối lượng khô của phần thân cây (g).
  7. PDW (Plant Dry Weight): Tổng khối lượng khô của cả cây (g).
  8. R_S (Root-Shoot Ratio): Tỷ lệ khối lượng khô giữa rễ và thân.
  9. DRW (Deep Root Weight): Khối lượng khô phần rễ ăn sâu.
  10. DRP (Deep Root Proportion): Phần trăm khối lượng khô của phần rễ ăn sâu so với tổng khối lượng rễ.
  11. SRP (Shallow Root Proportion): Phần trăm khối lượng khô của phần rễ ăn nông.
  12. DW0020: Khối lượng khô phần rễ ở tầng đất từ 0 đến 20 cm.
  13. DW2040: Khối lượng khô phần rễ ở tầng đất từ 20 đến 40 cm.
  14. DW4060: Khối lượng khô phần rễ ở tầng đất từ 40 đến 60 cm.
  15. DWB60: Khối lượng khô phần rễ ở độ sâu lớn hơn 60 cm.
  16. TIL (Number of Tillers): Số nhánh/cây.
  17. LLGTH (Longest Leaf Length): Chiều dài lá dài nhất (chiều dài thân).
  18. NR_T (Number of crown roots per Tiller): Số lượng rễ bất định trung bình trên một nhánh.

Nguồn dữ liệu và phần mềm phân tích

  • Phần mềm di truyền và tin sinh học: DArTsoft 7.4, DArWin 5, STRUCTURE, BEAGLE, TASSEL, RASTA, phần mềm R.
  • Phần mềm phân tích thống kê: SAS 9 (phân tích phương sai ANOVA, phân tích thành phần chính PCA, hệ số tương quan Pearson, hệ số di truyền nghĩa rộng $H^2$).
  • Cơ sở dữ liệu lưu trữ: Toàn bộ dữ liệu kiểu gen SNP được đăng tải và công khai trên hệ thống cơ sở dữ liệu TropGeneDB (CIRAD).

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

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

Chương này tổng kết toàn diện các cơ sở lý luận về hình thái học, giải phẫu học và sinh lý học của bộ rễ lúa. Tác giả phân tích chi tiết vai trò của hệ thống rễ trong việc thích ứng với các điều kiện bất lợi phi sinh học (hạn hán, ngập lụt, ngộ độc khoáng chất). Đồng thời, chương 1 điểm cứu chi tiết hơn 600 QTLs đã được công bố trên thế giới, các mạng lưới gen điều hòa phát triển rễ cảm ứng phytohormone (đặc biệt là auxin và cytokinin), các yếu tố phiên mã kiểm soát vùng trung tâm tĩnh (QC), mô phân sinh và quá trình sinh trưởng vách tế bào rễ. Tác giả cũng phân tích sâu nguyên lý kỹ thuật giải trình tự NGS, công nghệ GBS, phương pháp GWAS và vai trò của nguồn gen lúa bản địa Việt Nam.

Chương 2: Vật liệu và phương pháp nghiên cứu

Chương 2 mô tả chi tiết nguồn gốc xuất xứ của 270 mẫu giống lúa nghiên cứu; quy trình chiết tách ADN; quy trình tạo thư viện GBS sử dụng enzyme giới hạn và hệ thống chỉ thị mã vạch; phương pháp đánh giá 18 tính trạng kiểu hình rễ trong hệ thống ống rễ nhà lưới 6 tuần tuổi; các thuật toán tính toán thống kê (ANOVA, PCA, hệ số tương quan, $F_{ST}$, LD decay) và các mô hình GLM/MLM trong GWAS trên phần mềm TASSEL và R.

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

1. Đặc điểm nông sinh học cơ bản và cấu trúc di truyền qua chỉ thị DArT

Tập đoàn 270 giống lúa thể hiện sự phân hóa nông sinh học: thời gian sinh trưởng, chiều cao cây, số nhánh hữu hiệu, hình dạng hạt và tính chất nội nhũ (gạo nếp và tẻ). Phân tích 241 chỉ thị DArT trên 270 giống lúa đã xác lập cây phân loại di truyền có cấu trúc lưỡng cực rõ rệt, chia tập đoàn thành hai nhóm loài phụ chính:

  • Nhóm loài phụ indica gồm 168 mẫu giống.
  • Nhóm loài phụ japonica gồm 88 mẫu giống.
  • Các giống còn lại thuộc nhóm Aus/Boro và Sadri/Basmati.

2. Kết quả phân tích kiểu gen qua giải trình tự (GBS) và cấu trúc phân nhóm bằng chỉ thị SNP

Từ tập đoàn chọn lọc, quá trình GBS đã xây dựng thành công ma trận haplotype gồm 25.971 chỉ thị SNP phân bố bao phủ toàn bộ 12 nhiễm sắc thể, có hàm lượng thông tin đa hình (PIC) trung bình đạt 32,0% (0,32).

  • Cấu trúc phân nhóm giống indica Việt Nam (114 mẫu giống): Phân chia thành 6 phân nhóm phụ độc lập, ký hiệu từ I1 đến I6, mang đặc trưng sinh thái và vùng thu thập riêng biệt.
  • Cấu trúc phân nhóm giống japonica Việt Nam (62 mẫu giống): Phân chia thành 4 phân nhóm phụ, ký hiệu từ J1 đến J4.
  • Phân tích mất cân bằng liên kết (Linkage Disequilibrium - LD): Sự suy giảm LD theo khoảng cách vật lý diễn ra nhanh chóng trên cả 12 nhiễm sắc thể, tạo cơ sở xác định các khoảng tin cậy hẹp xung quanh các locus liên kết.

3. Kết quả đánh giá kiểu hình các tính trạng bộ rễ

Phân tích phương sai (ANOVA) và thống kê mô tả 18 tính trạng rễ trên 194 mẫu giống lúa cho thấy sự biến dị lớn và có ý nghĩa thống kê ($P < 0,001$) giữa các giống:

  • Phân tích thành phần chính (PCA) và vòng tròn tương quan chỉ ra mối quan hệ mật thiết giữa các tính trạng khối lượng rễ ở các tầng đất với tổng khối lượng rễ (RDW) và độ ăn sâu (DEPTH).
  • So sánh giá trị trung bình giữa các phân nhóm cho thấy: Phân nhóm I3 và I6 (thuộc indica) cùng với phân nhóm J1 và J3 (thuộc japonica) sở hữu bộ rễ ăn sâu nhất và đường kính rễ dày nhất. Ngược lại, các phân nhóm I1, I4, J2, J4 có các chỉ số rễ thấp nhất.
  • Luận án đã xác định các giống lúa địa phương ưu tú thuộc phân nhóm I3 có bộ rễ dài, dày vượt trội gồm: Blề Blậu Chớ (G205), Tẻ nương (G153), Khẩu Nậm Rinh (G189) và Khẩu Pe Lạnh (G155).

4. Kết quả lập bản đồ liên kết toàn hệ gen (GWAS) và xác định gen ứng viên

Phân tích GWAS trên 3 ma trận dữ liệu (toàn tập đoàn, nhóm indica, nhóm japonica) với ngưỡng p-value $< 10^{-4}$ đã phát hiện:

  • 88 QTLs liên kết có ý nghĩa với 18 tính trạng phát triển bộ rễ.
  • 28 QTLs liên kết đồng thời với nhiều hơn một tính trạng kiểu hình rễ (QTLs đa hiệu).
  • 33 QTLs nằm trực tiếp bên trong vùng trình tự mã hóa của các gen chức năng.
  • 2 vùng QTLs liên kết mạnh đặc biệt:
    • Một vùng QTLs liên kết chặt chẽ với tính trạng số lượng rễ bất định (NCR) trên nhiễm sắc thể số 11 (vùng đỉnh q45).
    • Một vùng QTLs liên kết chặt chẽ với tính trạng độ dày rễ (THK) trên nhiễm sắc thể số 2.
  • Gen ứng viên: Từ vị trí của 88 QTLs và phạm vi phân rã LD, xác định được 889 gen ứng viên. Trong đó, 407 gen đã được xác định và phân loại nhóm chức năng giả định; 24 gen ứng viên đã có công bố quốc tế chứng minh chức năng sinh học và hóa sinh liên quan trực tiếp đến cơ chế hình thành và phát triển của bộ rễ lúa hoặc chỉnh thể trực giao (ortholog) ở Arabidopsis.

Bảng tổng hợp các thông số kỹ thuật và kết quả chính của luận án

Hạng mục / Chỉ số Giá trị định lượng trong luận án Ghi chú kỹ thuật / Chi tiết
Số lượng mẫu giống ban đầu 270 mẫu giống 214 giống VN, 33 giống đối chứng CIRAD, 23 giống Viện DTNN
Số lượng chỉ thị DArT phân tích 241 chỉ thị Phân chia 168 indica, 88 japonica, Aus/Boro và Basmati
Số lượng chỉ thị SNP từ GBS 25.971 SNPs Đã lọc MAF > 5% và nội suy dữ liệu khuyết (BEAGLE)
Hàm lượng thông tin đa hình (PIC) 32,0% (0,32) Giá trị trung bình trên toàn hệ gen
Phân nhóm cấu trúc giống indica 6 phân nhóm (I1 đến I6) Phân tích trên 114 mẫu giống lúa indica Việt Nam
Phân nhóm cấu trúc giống japonica 4 phân nhóm (J1 đến J4) Phân tích trên 62 mẫu giống lúa japonica Việt Nam
Số lượng tính trạng kiểu hình rễ 18 tính trạng NCR, THK, MRL, DEPTH, RDW, SDW, PDW, R_S, DRW, DRP,...
Số lượng giống đánh giá kiểu hình 194 mẫu giống Đánh giá ở giai đoạn cây 6 tuần tuổi trong hệ thống ống rễ
Tổng số QTLs xác định qua GWAS 88 QTLs Ngưỡng ý nghĩa thống kê $P < 10^{-4}$ trên mô hình MLM
Số QTLs liên kết đa tính trạng 28 QTLs Liên kết đồng thời với từ 2 tính trạng trở lên
Số QTLs nằm trong gen chức năng 33 QTLs Nằm trực tiếp trong vùng trình tự mã hóa
Vùng QTL liên kết số lượng rễ (NCR) Trên Nhiễm sắc thể số 11 Vùng đỉnh q45, liên kết chặt chẽ với số rễ bất định
Vùng QTL liên kết độ dày rễ (THK) Trên Nhiễm sắc thể số 2 Liên kết chặt chẽ với độ dày (đường kính) rễ
Tổng số gen ứng viên định vị 889 gen ứng viên Nằm trong vùng tin cậy phân rã LD của 88 QTLs
Số gen ứng viên có chức năng giả định 407 gen Đã được phân nhóm chức năng sinh học
Số gen ứng viên đã rõ chức năng rễ 24 gen Đã có công bố chứng minh vai trò hóa sinh/sinh học

Kết quả và những đóng góp mới

Đóng góp mới về mặt khoa học

  1. Đánh giá toàn diện, có hệ thống về đặc điểm nông sinh học cơ bản và mức độ đa dạng di truyền của tập đoàn 270 giống lúa (trong đó có 214 giống lúa địa phương Việt Nam) bằng công nghệ vi mảng DArT và giải trình tự GBS.
  2. Xây dựng và công bố bộ dữ liệu kiểu gen haplotype mật độ cao gồm 25.971 chỉ thị SNP đa hình cao (PIC trung bình đạt 32,0%), bao phủ toàn bộ 12 nhiễm sắc thể của cây lúa trên cơ sở dữ liệu quốc tế TropGeneDB.
  3. Xác lập cấu trúc di truyền chi tiết của nguồn gen lúa bản địa Việt Nam, phân định nhóm indica thành 6 phân nhóm (I1 - I6) và nhóm japonica thành 4 phân nhóm (J1 - J4) gắn liền với các vùng địa lý sinh thái.
  4. Cung cấp bộ dữ liệu kiểu hình chuẩn hóa gồm 18 tính trạng liên quan đến cấu trúc và khối lượng bộ rễ ở 194 giống lúa Việt Nam.
  5. Ứng dụng thành công phương pháp GWAS phát hiện 88 QTLs liên kết với 18 tính trạng rễ, xác định 889 gen ứng viên (trong đó có 407 gen có chức năng giả định và 24 gen đã được xác nhận vai trò sinh học với rễ), đồng thời phát hiện 2 vùng QTLs liên kết chặt chẽ: vùng liên kết với số lượng rễ bất định (NCR) trên NST 11 và vùng liên kết với độ dày rễ (THK) trên NST 2.

Đóng góp mới về mặt thực tiễn và kiến nghị

  • Cung cấp nguồn vật liệu chọn giống: Xác định được các giống lúa địa phương mang đặc tính rễ sâu và rễ dày vượt trội thuộc phân nhóm I3 gồm: Blề Blậu Chớ (G205), Tẻ nương (G153), Khẩu Nậm Rinh (G189), Khẩu Pe Lạnh (G155). Đây là nguồn vật liệu bố mẹ phục vụ các chương trình lai tạo cải tiến tính trạng bộ rễ nhằm nâng cao khả năng chống chịu điều kiện bất lợi.
  • Khai thác dữ liệu mở: Bộ dữ liệu 25.971 SNP trên TropGeneDB là nguồn tài nguyên mở để các nhà khoa học tiếp tục khai thác nghiên cứu GWAS cho các tính trạng nông học quan trọng khác như năng suất, chất lượng hạt, cấu trúc bông hoặc tính chống chịu sâu bệnh và stress phi sinh học.
  • Kiến nghị: Tiếp tục nghiên cứu sâu chức năng sinh học phân tử của các gen ứng viên tại vùng QTL trên nhiễm sắc thể số 11 và số 2 bằng các kỹ thuật chuyển gen, bất hoạt gen hoặc siêu biểu hiện gen; ứng dụng chỉ thị phân tử liên kết để chọn lọc hỗ trợ chỉ thị (MAS) trong công tác chọn giống lúa thích ứng với biến đổi khí hậu.

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

Hạn chế của nghiên cứu

  • Thí nghiệm đánh giá kiểu hình bộ rễ được thực hiện trên hệ thống ống rễ cải tiến trong điều kiện nhà lưới ở giai đoạn cây 6 tuần tuổi (sinh trưởng sinh dưỡng), chưa phản ánh toàn bộ động thái phát triển của bộ rễ trong điều kiện đồng ruộng tự nhiên qua các giai đoạn trỗ và chín.
  • Các gen ứng viên được xác định chủ yếu dựa trên vị trí liên kết thống kê của GWAS và phân tích tin sinh học trong vùng phân rã LD; luận án chưa thực hiện các thí nghiệm kiểm chứng chức năng sinh học trực tiếp (in vivo) như gây đột biến định hướng, knock-out gen hoặc tạo cây chuyển gen siêu biểu hiện cho từng gen ứng viên cụ thể.

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

  • Triển khai các nghiên cứu kiểm chứng chức năng hóa sinh và sinh học đối với các gen ứng viên triển vọng trong vùng QTL điều khiển số lượng rễ bất định (NCR) trên nhiễm sắc thể số 11 và độ dày rễ (THK) trên nhiễm sắc thể số 2.
  • Sử dụng các dòng lúa ưu tú thuộc phân nhóm I3 (Blề Blậu Chớ, Tẻ nương, Khẩu Nậm Rinh, Khẩu Pe Lạnh) làm vật liệu khởi đầu trong các chương trình lai tạo chọn giống cải tiến bộ rễ chống chịu hạn.
  • Tận dụng ma trận kiểu gen 25.971 SNP để mở rộng phân tích GWAS cho các tính trạng thích ứng khác như tính chịu mặn, chịu ngập, hiệu quả sử dụng dinh dưỡng khoáng và các yếu tố cấu thành năng suất.

Giá trị tham khảo

Luận án là tài liệu tham khảo trong nghiên cứu di truyền học và công nghệ sinh học nông nghiệp:

  • Đối với nghiên cứu sinh, học viên cao học và giảng viên: Cung cấp tài liệu tổng quan về các con đường tín hiệu di truyền điều hòa rễ lúa (mạng lưới auxin, cytokinin, các yếu tố phiên mã LBD, WOX, NAC, GRAS); đồng thời cung cấp quy trình phương pháp luận thực hiện nghiên cứu GWAS từ bước xử lý kiểu gen GBS, đánh giá kiểu hình đến mô hình thống kê MLM.
  • Đối với các nhà chọn tạo giống cây trồng: Cung cấp thông tin kiểu hình chi tiết của 194 nguồn gen lúa Việt Nam, danh mục các giống mang alen rễ sâu - dày, cùng hệ thống chỉ thị SNP và QTLs phục vụ công tác chọn giống nhờ chỉ thị phân tử (Marker-Assisted Selection - MAS).
  • Đối với các nhà quản lý và hoạch định chính sách bảo tồn nguồn gen: Cung cấp cơ sở dữ liệu về cấu trúc di truyền phân tử của nguồn gen lúa bản địa Việt Nam, hỗ trợ công tác bảo tồn và khai thác bền vững tài nguyên di truyền thực vật quốc gia.

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

1. Tập đoàn mẫu giống lúa trong luận án gồm bao nhiêu mẫu và được phân loại cấu trúc di truyền như thế nào?

Nghiên cứu sử dụng tập đoàn ban đầu gồm 270 mẫu giống lúa (214 giống địa phương Việt Nam từ nhiều vùng sinh thái, 33 giống đối chứng từ CIRAD và 23 giống từ Viện Di truyền Nông nghiệp). Phân tích 241 chỉ thị DArT chia tập đoàn thành cấu trúc lưỡng cực gồm 168 giống thuộc nhóm indica, 88 giống thuộc nhóm japonica, cùng một số giống thuộc nhóm Aus/Boro và Sadri/Basmati. Khi phân tích sâu bằng chỉ thị SNP từ GBS, nhóm indica Việt Nam phân thành 6 phân nhóm phụ (I1 - I6) và nhóm japonica phân thành 4 phân nhóm phụ (J1 - J4).

2. Bộ dữ liệu kiểu gen GBS được xây dựng có đặc điểm kỹ thuật gì?

Bộ dữ liệu kiểu gen GBS được xây dựng từ 200 mẫu giống lúa đại diện, sau quá trình lọc các alen có tần số nhỏ (MAF < 5%) và nội suy dữ liệu khuyết bằng phần mềm BEAGLE, đã thu được ma trận haplotype gồm 25.971 chỉ thị SNP đa hình chất lượng cao bao phủ 12 nhiễm sắc thể, có hàm lượng thông tin đa hình (PIC) trung bình đạt 32,0% và được công bố công khai trên cơ sở dữ liệu TropGeneDB.

3. Những phân nhóm và giống lúa nào thể hiện đặc tính bộ rễ vượt trội nhất?

Kết quả đánh giá 18 tính trạng kiểu hình rễ trên 194 giống lúa 6 tuần tuổi cho thấy phân nhóm I3 và I6 (thuộc nhóm indica) cùng phân nhóm J1 và J3 (thuộc nhóm japonica) có bộ rễ ăn sâu nhất và dày nhất. Trong đó, 4 giống lúa thuộc phân nhóm I3 thể hiện ưu thế vượt trội về chiều dài và độ dày rễ là: Blề Blậu Chớ (G205), Tẻ nương (G153), Khẩu Nậm Rinh (G189) và Khẩu Pe Lạnh (G155).

4. Kết quả GWAS đã phát hiện bao nhiêu QTLs và có những vùng liên kết nổi bật nào?

Phân tích GWAS trên mô hình tuyến tính hỗn hợp (MLM) với ngưỡng ý nghĩa $P < 10^{-4}$ đã xác định được 88 QTLs liên kết với 18 tính trạng rễ. Trong đó có 28 QTLs liên kết đa tính trạng và 33 QTLs nằm trong vùng mã hóa của gen chức năng. Hai vùng QTLs liên kết mạnh đặc biệt được ghi nhận gồm: vùng liên kết với số lượng rễ bất định (NCR) trên nhiễm sắc thể số 11 (vùng q45) và vùng liên kết với độ dày rễ (THK) trên nhiễm sắc thể số 2.

5. Luận án đã định vị được bao nhiêu gen ứng viên và bao nhiêu gen đã được xác nhận chức năng sinh học?

Dựa trên vị trí của 88 QTLs và phạm vi phân rã mất cân bằng liên kết (LD), nghiên cứu đã xác định được 889 gen ứng viên. Trong số này, 407 gen đã được phân nhóm chức năng giả định và 24 gen đã có các công bố khoa học quốc tế chứng minh chức năng hóa sinh, sinh học liên quan trực tiếp đến sự phát sinh và phát triển bộ rễ ở lúa hoặc cơ chế trực giao tương đương ở Arabidopsis.

Kết luận

Luận án của tác giả Phùng Thị Phương Nhung đã xây dựng thành công bộ dữ liệu kiểu gen gồm 25.971 chỉ thị SNP qua kỹ thuật GBS và bộ dữ liệu kiểu hình 18 tính trạng rễ trên tập đoàn giống lúa Việt Nam. Nghiên cứu đã xác lập cấu trúc di truyền gồm 6 phân nhóm indica và 4 phân nhóm japonica, đồng thời nhận diện được các giống lúa địa phương ưu tú có bộ rễ ăn sâu và rễ dày thuộc phân nhóm I3 như Blề Blậu Chớ, Tẻ nương, Khẩu Nậm Rinh và Khẩu Pe Lạnh. Thông qua phương pháp lập bản đồ GWAS, luận án đã xác định được 88 QTLs cùng 889 gen ứng viên liên quan đến sự phát triển bộ rễ, nổi bật với 2 vùng QTLs liên kết chặt trên nhiễm sắc thể số 11 (tính trạng số lượng rễ bất định) và nhiễm sắc thể số 2 (tính trạng độ dày rễ). Các kết quả này cung cấp cơ sở dữ liệu di truyền và nguồn vật liệu phục vụ công tác nghiên cứu chức năng gen và chọn tạo giống lúa thích ứng với biến đổi khí hậu.