Bối cảnh và vấn đề nghiên cứu

Lúa (Oryza sativa) là một trong những cây lương thực quan trọng nhất thế giới, cung cấp nguồn năng lượng chính cho hơn một nửa dân số toàn cầu và có lịch sử thuần hóa lâu đời tại châu Á. Bên cạnh các giống lúa tẻ và lúa nếp có vỏ hạt màu trắng thông thường, các giống lúa có sắc tố vỏ cám (đỏ, tím, nâu, đen) tích lũy hàm lượng cao các hợp chất phi diệp lục như anthocyanin, flavanol và proanthocyanidin. Tại Việt Nam, lúa nếp cẩm (nếp tím) là nguồn gen bản địa đặc sản quý giá, được canh tác lâu đời tại nhiều vùng sinh thái, đặc biệt là các tỉnh miền núi phía Bắc như Yên Bái, Sơn La, Hòa Bình, và phân bố rải rác tại Phú Thọ, Ninh Bình hay Đồng bằng sông Cửu Long (Long An, Cần Thơ). Lúa nếp cẩm có giá trị dinh dưỡng cao (hàm lượng protein cao hơn 6,8%, chất béo cao hơn 20% so với gạo thường, giàu acid amin, carotene và các nguyên tố vi lượng như Fe, Zn), đồng thời giữ vai trò quan trọng trong đời sống văn hóa ẩm thực truyền thống.

Tuy nhiên, các nghiên cứu ở cấp độ phân tử về cơ chế tổng hợp sắc tố anthocyanin và proanthocyanidin trên các giống lúa nếp đặc sản bản địa tại Việt Nam còn hạn chế, gây khó khăn cho công tác bảo tồn, khai thác và định hướng chọn tạo giống. Về mặt di truyền, gen Rc nằm trên nhiễm sắc thể số 7, mã hóa cho yếu tố điều hòa phiên mã thuộc họ basic helix-loop-helix (bHLH), đóng vai trò then chốt trong con đường tích lũy proanthocyanidin ở vỏ cám lúa. Hầu hết các giống lúa vỏ trắng (>97%) mang alen đột biến mất đoạn 14 bp làm dịch khung và tạo protein ngắn mất chức năng, hoặc mang alen Rc-s đột biến điểm chuyển dạng C sang A. Nhằm làm sáng tỏ đặc điểm di truyền phân tử của nguồn gen bản địa, khóa luận tốt nghiệp tập trung phân tích đa hình và cấu trúc gen Rc ở giống lúa Nếp Cẩm Yên Bái.

Nhiệm vụ nghiên cứu cụ thể gồm 4 nội dung:

  1. Tối ưu hóa điều kiện phản ứng PCR để nhân bản 18 phân đoạn bao phủ toàn bộ gen Rc.
  2. Giải trình tự 18 phân đoạn DNA và lắp ráp trình tự hoàn chỉnh của gen Rc.
  3. Phân tích cấu trúc exon - intron và đa hình trình tự nucleotide (SNP, indel) của gen Rc ở giống Nếp Cẩm Yên Bái.
  4. Xây dựng cây phát sinh chủng loại của giống Nếp Cẩm Yên Bái dựa trên trình tự đầy đủ của gen Rc so với các giống lúa đại diện trên thế giới.
  • Đối tượng nghiên cứu: Giống lúa nếp cẩm bản địa thu thập tại tỉnh Yên Bái (NepCam YB), do Bộ môn Sinh học phân tử và Công nghệ sinh học ứng dụng, Khoa Công nghệ sinh học, Học viện Nông nghiệp Việt Nam cung cấp.
  • Phạm vi nghiên cứu: Phân tích di truyền phân tử và tin sinh học đối với gen Rc của mẫu giống Nếp Cẩm Yên Bái, thực hiện tại Phòng thí nghiệm Bộ môn Sinh học phân tử, Khoa Công nghệ sinh học, Học viện Nông nghiệp Việt Nam, hoàn thành vào tháng 03/2022.

Cơ sở lý thuyết và phương pháp

Cơ sở lý thuyết và khung phân tích

  • Lịch sử tiến hóa và thuần hóa lúa: Quá trình thuần hóa từ tổ tiên hoang dại Oryza rufipogon và Oryza nivara sang lúa trồng Oryza sativa (gồm hai phân loài chính indica và japonica) gắn liền với các đột biến làm giảm rụng hạt, thu ngắn râu hạt (awn) và thay đổi màu sắc vỏ cám.
  • Con đường sinh tổng hợp Anthocyanin và Proanthocyanidin (ABP): Sắc tố proanthocyanidin (condensed tannin) là một nhánh của con đường tổng hợp flavonoid/anthocyanin, khởi đầu từ phenylalanine hoặc malonyl-CoA tạo chalcone qua hàng loạt phản ứng enzyme.
  • Vai trò của gen Rc và Rd: Màu sắc vỏ cám được kiểm soát bởi tương tác bổ trợ giữa hai locus: kiểu gen RcRd cho màu đỏ, Rcrd cho màu nâu, còn rcRd và rcrd cho màu trắng. Locus Rd mã hóa enzyme Dihydroflavonol-4-reductase (DFR), trong khi locus Rc mã hóa yếu tố phiên mã bHLH. Sáu alen chính của Rc đã được ghi nhận trên thế giới gồm Rc (kiểu dại), rc (mất đoạn 14 bp), Rc-s (đột biến điểm C→A), Rc-g và Rcr (đột biến hồi biến phục hồi khung đọc), và Rc-g1 (đột biến A→T đặc thù ở lúa châu Phi).

Phương pháp nghiên cứu thực nghiệm

  • Tách chiết DNA tổng số: Tách chiết DNA từ lá non theo quy trình Dellaporta (1983). Mẫu lá khô chân không được nghiền bằng máy rung hạt sắt (Multi-beads shocker, Yasui Kikai) ở tốc độ 1800 vòng/phút trong 60 giây (lặp lại 2 lần), ủ trong dung dịch đệm chiết (1M Tris-HCl pH 8,0; 5M NaCl; 10% SDS) ở 65°C, kết tủa protein bằng 5M potassium acetate, tủa DNA bằng isopropanol, rửa sạch bằng ethanol 70% và hòa tan trong đệm TE (10mM Tris-HCl pH 8,0; 1mM EDTA pH 8,0). Kiểm tra độ sạch và nồng độ DNA bằng điện di gel agarose 1%.
  • Nhân bản phân đoạn gen (PCR): Sử dụng 18 cặp mồi xuôi và ngược (ký hiệu AF1 đến AF18, thiết kế gối đầu dọc theo gen Rc theo Li et al., 2014) để khuếch đại các đoạn DNA có kích thước dưới 1000 bp. Thành phần phản ứng 10 µL gồm: 5 µL 2X PCR Master Mix, 1 µL hỗn hợp mồi xuôi và ngược, 1 µL DNA khuôn và 3 µL nước khử nuclease. Chu trình nhiệt: biến tính ban đầu 95°C trong 5 phút; 35 chu kỳ (95°C trong 30 giây, bắt mồi ở 45°C – 62°C trong 30 giây tùy cặp mồi, kéo dài ở 72°C trong 1 phút); kéo dài lần cuối 72°C trong 7 phút; giữ ở 15°C. Sản phẩm PCR được kiểm tra bằng điện di gel agarose 2% với thang chuẩn KAPPA Universal Ladder 100 bp.
  • Giải trình tự DNA: Sản phẩm PCR đạt chất lượng được tinh sạch và giải trình tự theo phương pháp Sanger tại Công ty 1st BASE (Singapore).

Phương pháp phân tích tin sinh học và nguồn dữ liệu

  • Kiểm tra chất lượng sắc phổ giải trình tự bằng phần mềm Chromas; loại bỏ nhiễu "N".
  • Ghép nối trình tự (assembly) và gióng hàng (alignment) bằng phần mềm BioEdit và MEGA X.
  • Xác định cấu trúc exon - intron và hiển thị mô hình gen bằng công cụ trực tuyến GSDS 2.0 (Gene Structure Display Server), đối chiếu với cơ sở dữ liệu Nipponbare GRASP5.0.
  • Phân tích đa hình nucleotide (SNP, indel) bằng DnaSP v6.0 và BioEdit.
  • Đánh giá thông số lý hóa protein (trọng lượng phân tử, điểm đẳng điện pI, chỉ số không bền vững Instability Index, chỉ số béo Aliphatic Index, chỉ số GRAVY) và vị trí định vị dưới tế bào bằng công cụ TargetP (mTP, cTP, SP).
  • Xây dựng cây phát sinh loài bằng phương pháp Maximum Likelihood (ML) trong phần mềm MEGA X.
  • Dữ liệu đối chiếu từ NCBI: Gồm trình tự tham chiếu chuẩn của lúa Oryza sativa (mã số KX549256) và 8 trình tự gen Rc đại diện cho 4 phân nhóm lúa (indica, aus, temperate japonica, tropical japonica) với cả hai kiểu hình hạt đỏ và hạt trắng:
STT Tên giống lúa Mã gen trên NCBI Màu sắc hạt Phân nhóm lúa Xuất xứ
1 Bengkongang DQ885810.1 Đỏ Indica Indonesia
2 Kasalath DQ885812.1 Đỏ Aus Ấn Độ
3 Gangdodo DQ885805.1 Đỏ Temperate japonica Hàn Quốc
4 Pae Daya Indolobye DQ885808.1 Đỏ Tropical japonica Indonesia
5 Dhala Shaitta DQ885821.1 Trắng Aus Bangladesh
6 Dee Geo Woo Gen DQ885818.1 Trắng Indica Đài Loan
7 Koshihikari DQ885803.1 Trắng Temperate japonica Nhật Bản
8 Jefferson DQ885802.1 Trắng Tropical japonica Hoa Kỳ

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

Chương 1: Tổng quan tài liệu (Review of Literature)

Chương này tổng kết có hệ thống các công trình nghiên cứu về:

  • Lịch sử tiến hóa và thuần hóa của chi lúa Oryza, sự phân tách giữa các nhóm indica và japonica, cùng những biến đổi hình thái như râu hạt, rụng hạt và màu hạt.
  • Đặc điểm nông sinh học và giá trị dinh dưỡng của lúa nếp cẩm, phân bố canh tác tại các nước châu Á (Trung Quốc chiếm 62% sản lượng, tiếp theo là Sri Lanka, Indonesia, Ấn Độ, Philippines, Thái Lan) và các vùng sinh thái ở Việt Nam.
  • Tổng quan cấu trúc hóa học của anthocyanin và proanthocyanidin, con đường biến dưỡng ABP, các gốc R khác nhau của anthocyanidin và ứng dụng trong công nghiệp thực phẩm, dược phẩm.
  • Cơ chế di truyền điều hòa màu sắc hạt bởi gen Rc và Rd, các kiểu đột biến mất đoạn 14 bp, đột biến điểm Rc-s, các đột biến hồi biến (Rc-g, Rcr), và sự phân bố của 6 alen Rc trên thế giới.
  • Đa dạng nguồn gen lúa bản địa tại Việt Nam: hiện Ngân hàng gen cây trồng quốc gia đang lưu giữ hơn 8.000 nguồn gen lúa (khu vực Tây Bắc chiếm 41% nguồn gen lúa nương thu thập). Điểm lại các nghiên cứu đánh giá đa dạng di truyền lúa nếp bằng chỉ thị phân tử SSR của Doãn Thanh Quỳnh (2016) và Ngô Thị Hồng Tươi (2014), cùng nghiên cứu cấu trúc quần thể 672 hệ gen lúa Việt Nam của Higgins et al. (2021).

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

Trình bày chi tiết về nguồn gốc mẫu lúa Nếp Cẩm Yên Bái, hóa chất, thiết bị và các bước thực nghiệm tách chiết DNA lá non. Chương này công bố chi tiết bảng 18 cặp mồi PCR với nhiệt độ bắt mồi ($T_a$) dao động từ 45°C đến 62°C và kích thước sản phẩm khuếch đại dự kiến từ 432 bp (AF8) đến 986 bp (AF15). Đồng thời, quy định chi tiết quy trình xử lý dữ liệu thô, gióng hàng chuỗi, phân tích cấu trúc gen và dựng cây phân loại bằng các phần mềm tin sinh học chuyên dụng.

Chương 3: Kết quả và thảo luận (Results and Discussion)

Chương này phản ánh toàn bộ dữ liệu thực nghiệm và phân tích tin sinh học:

  • Khuếch đại và giải trình tự 18 đoạn gen Rc: 18 cặp mồi đã nhân bản thành công các phân đoạn đơn băng, rõ nét trên gel agarose 2%. Ở lần giải trình tự đầu tiên, các đoạn AF4, AF7, AF17, AF18 chưa đạt chất lượng và AF12 gặp chuỗi lặp poly(AT); sau khi tối ưu và giải trình tự lại lần hai, toàn bộ 18 phân đoạn đã đạt độ nét cao trên Chromas và được ghép nối thành công.
  • Trình tự hoàn chỉnh và mô hình gen: Trình tự ghép nối gen Rc của Nếp Cẩm Yên Bái đạt tổng chiều dài 7500 bp (bao gồm cả vùng điều hòa 5' upstream và 3' downstream), có độ tương đồng 97% với trình tự chuẩn Oryza sativa KX549256. Mô hình gen hiển thị qua GSDS 2.0 xác định gen Rc ở Nếp Cẩm Yên Bái gồm 9 exon (hộp màu vàng) và 8 intron (đường kẻ đen).
  • Phân tích đột biến mất đoạn 14 bp: Khác với lúa trắng mang mất đoạn 14 bp ở vùng mã hóa gây mất chức năng protein, trên gen Rc của Nếp Cẩm Yên Bái phát hiện có một đoạn khuyết 14 bp nằm tại intron thứ 5 (từ vị trí nucleotide 5061 đến 5075). Phân tích đồng tuyến (collinearity) khẳng định giống Nếp Cẩm Yên Bái mang kiểu gen Rc hoạt động chức năng quy định hạt có sắc tố.
  • Đa hình nucleotide (SNP): So sánh với 8 giống lúa đối chứng quốc tế, phát hiện tổng cộng 54 vị trí SNP trên gen Rc của Nếp Cẩm Yên Bái, trong đó 33 SNP nằm ở các vùng exon và 21 SNP nằm ở các vùng intron.
  • Hàm lượng G-C và A-T: Tỷ lệ GC trung bình của 8 giống đối chứng dao động quanh mức 37%, trong khi tỷ lệ GC của gen Rc ở Nếp Cẩm Yên Bái đạt 35,55% (mức thấp nhất trong nhóm so sánh), phản ánh tính biến động và kém bền vững hơn của cấu trúc gen.
  • Đặc tính lý hóa và định vị dưới tế bào của protein Rc: Chuỗi polypeptide mã hóa dao động từ 445 đến 636 acid amin, trọng lượng phân tử từ 49615,34 đến 69796,71 Da. Chỉ số không bền vững (Instability Index) đạt giá trị > 40 (đạt tới 58,29), chứng minh protein Rc thuộc nhóm không bền (unstable). Chỉ số GRAVY mang giá trị âm (từ -0,553 đến -0,590), khẳng định đây là protein hình cầu ưa nước (hydrophilic globular protein). Công cụ TargetP dự đoán protein Rc của Nếp Cẩm Yên Bái định vị tại ty thể (mitochondria) với chỉ số peptide đích ty thể mTP = 0,535 (cao hơn hẳn cTP và SP).
  • Phát sinh chủng loại: Cây phát sinh loài theo phương pháp Maximum Likelihood cho thấy giống Nếp Cẩm Yên Bái có quan hệ di truyền gần gũi nhất với giống lúa Jefferson (thuộc phân nhóm Tropical japonica, xuất xứ Hoa Kỳ).

Kết quả và đóng góp

Kết quả chính kèm số liệu

  • Nhân bản và tối ưu hóa thành công điều kiện PCR cho 18 cặp mồi phân đoạn gen Rc.
  • Lắp ráp thành công trình tự hoàn chỉnh gen Rc ở giống Nếp Cẩm Yên Bái với chiều dài 7500 bp, đạt độ tương đồng 97% với trình tự chuẩn KX549256 trên NCBI.
  • Xác định cấu trúc gen Rc của Nếp Cẩm Yên Bái gồm 9 exon và 8 intron.
  • Phát hiện 54 điểm đa hình đơn nucleotide (SNP) (33 SNP tại exon, 21 SNP tại intron) và một đột biến mất đoạn 14 bp tại intron thứ 5 (vị trí 5061–5075).
  • Xác định hàm lượng GC của gen là 35,55%; protein Rc có tính chất ưa nước (GRAVY âm), không bền vững (Instability index > 40), định vị tại ty thể (mTP = 0,535).
  • Xác lập vị trí phát sinh loài của Nếp Cẩm Yên Bái có quan hệ gần gũi nhất với nhóm lúa Tropical japonica (cụ thể là giống Jefferson).

Đóng góp của nghiên cứu

  • Công bố bộ dữ liệu trình tự đầy đủ và cấu trúc phân tử hoàn chỉnh đầu tiên của gen Rc trên giống lúa nếp cẩm đặc sản Yên Bái.
  • Cung cấp dẫn liệu di truyền phân tử quan trọng phục vụ cho việc nghiên cứu nguồn gốc, phân loại học và quá trình thuần hóa các giống lúa bản địa Việt Nam.
  • Thiết lập hệ thống 18 cặp mồi và quy trình PCR chuẩn hóa hỗ trợ hiệu quả cho các nghiên cứu giải trình tự và sàng lọc gen Rc trên các nguồn gen lúa khác.

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

  • Hạn chế: Nghiên cứu mới dừng lại ở mức độ giải trình tự gen bộ gen (genomic DNA), phân tích tin sinh học về cấu trúc và dự đoán lý hóa protein; chưa tiến hành đánh giá mức độ biểu hiện gen ở mức phiên mã (mRNA/RT-PCR) và dịch mã trong từng giai đoạn phát triển hạt, cũng như chưa phân tích thực nghiệm hàm lượng anthocyanin cụ thể chịu ảnh hưởng của điều kiện sinh thái môi trường.
  • Hướng nghiên cứu tiếp:
    • Mở rộng giải trình tự và phân tích cấu trúc gen Rc trên nhiều giống lúa nếp cẩm bản địa khác thu thập từ các vùng sinh thái khác nhau tại Việt Nam.
    • Phân tích chức năng sinh học của gen Rc trong cơ chế sinh tổng hợp anthocyanin và đánh giá tương tác giữa yếu tố di truyền với điều kiện môi trường canh tác đối với hàm lượng sắc tố tích lũy trong hạt.

Giá trị tham khảo

  • Đối tượng hưởng lợi: Sinh viên, học viên cao học, nghiên cứu sinh và cán bộ nghiên cứu chuyên ngành Công nghệ sinh học, Di truyền học thực vật, Bảo tồn nguồn gen và Chọn giống cây trồng.
  • Nội dung tham khảo giá trị:
    • Bảng trình tự 18 cặp mồi PCR gối đầu và dải nhiệt độ bắt mồi tối ưu để khuếch đại toàn bộ gen Rc kích thước lớn (~6,4 – 7,5 kb).
    • Quy trình kỹ thuật kết hợp giữa sinh học phân tử thực nghiệm (tách chiết DNA, PCR, Sanger) và các phần mềm tin sinh học (BioEdit, MEGA X, GSDS 2.0, TargetP).
    • Bộ dữ liệu trình tự gen Rc 7500 bp của giống Nếp Cẩm Yên Bái làm trình tự đối chiếu cho các nghiên cứu đa hình di truyền lúa màu.

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

1. Gen Rc ở giống Nếp Cẩm Yên Bái có chiều dài bao nhiêu và có cấu trúc exon - intron như thế nào?

Trình tự hoàn chỉnh của gen Rc (bao gồm cả vùng 5' upstream và 3' downstream) ở giống Nếp Cẩm Yên Bái sau khi lắp ráp và gióng hàng có chiều dài 7500 bp. Cấu trúc mô hình gen được xác định gồm 9 exon và 8 intron, tương đồng với mô hình cấu trúc gen Rc được công bố bởi Sweeney et al. (2007).

2. Đột biến mất đoạn 14 bp trên gen Rc của Nếp Cẩm Yên Bái nằm ở vị trí nào và có làm mất màu hạt không?

Trên gen Rc của giống Nếp Cẩm Yên Bái, đột biến mất đoạn 14 bp được phát hiện nằm tại intron thứ 5 (vị trí nucleotide từ 5061 đến 5075). Đột biến này không nằm ở vùng mã hóa exon nên không làm dịch khung đọc hay cắt cụt yếu tố phiên mã bHLH. Do đó, giống lúa này vẫn duy trì kiểu gen Rc hoạt động chức năng, tạo vỏ cám có màu sắc tím sẫm đặc trưng.

3. Phân tích đa hình chuỗi đã ghi nhận bao nhiêu điểm SNP trên gen Rc của Nếp Cẩm Yên Bái so với các giống đối chứng?

Khi so sánh với 8 trình tự đại diện thuộc các nhóm lúa indica, aus, temperate japonica và tropical japonica từ NCBI, nghiên cứu phát hiện tổng cộng 54 điểm đa hình đơn nucleotide (SNP) trên gen Rc của Nếp Cẩm Yên Bái, bao gồm 33 SNP nằm tại các vùng exon và 21 SNP nằm tại các vùng intron.

4. Các thông số lý hóa và phân tích tin sinh học dự đoán protein Rc của Nếp Cẩm Yên Bái có đặc tính gì?

Protein Rc của Nếp Cẩm Yên Bái là một protein hình cầu ưa nước (chỉ số GRAVY âm, từ -0,553 đến -0,590) và có tính chất không bền vững (chỉ số Instability Index > 40, đạt 58,29). Theo phân tích từ công cụ TargetP, protein Rc của giống lúa này được định vị trong ty thể (mitochondria) với chỉ số peptide đích ty thể mTP đạt 0,535.

5. Cây phát sinh chủng loại chỉ ra Nếp Cẩm Yên Bái có mối quan hệ di truyền gần nhất với giống lúa nào?

Phân tích cây phát sinh chủng loại bằng phương pháp Maximum Likelihood (ML) dựa trên trình tự đầy đủ của gen Rc cho thấy giống Nếp Cẩm Yên Bái có mối quan hệ di truyền gần gũi nhất với giống lúa Jefferson (một giống lúa hạt trắng thuộc phân nhóm Tropical japonica có xuất xứ từ Hoa Kỳ).


Kết luận

Khóa luận tốt nghiệp đã tối ưu hóa thành công quy trình PCR với 18 cặp mồi và xác định được trình tự hoàn chỉnh dài 7500 bp của gen Rc ở giống lúa Nếp Cẩm Yên Bái, với cấu trúc gồm 9 exon và 8 intron. Phân tích tin sinh học đã làm sáng tỏ 54 điểm SNP, phát hiện đột biến mất đoạn 14 bp nằm tại intron thứ 5, đồng thời xác định protein Rc có tính chất ưa nước, không bền vững và định vị tại ty thể. Cây phát sinh loài cho thấy mối quan hệ di truyền gần gũi của Nếp Cẩm Yên Bái với nhóm Tropical japonica, cung cấp nguồn dữ liệu phân tử quan trọng phục vụ công tác bảo tồn nguồn gen và nghiên cứu tiến hóa lúa bản địa tại Việt Nam.