Tổng quan luận án

Cá tra nuôi (Pangasianodon hypophthalmus) là đối tượng thủy sản nuôi chủ lực tại Việt Nam theo Quyết định số 50/2018/QĐ-TTg ngày 13/12/2018 của Thủ tướng Chính phủ, phục vụ cho xuất khẩu và tiêu thụ nội địa. Tuy nhiên, ngành sản xuất cá tra đối mặt với nhiều thách thức, đặc biệt là tình trạng suy giảm chất lượng con giống và rủi ro dịch bệnh. Nâng cao chất lượng nguồn giống thông qua chọn giống là nhiệm vụ cấp thiết nhằm phát triển bền vững ngành hàng này. Tăng trưởng là tính trạng số lượng quan trọng nhất trong các chương trình chọn giống thủy sản. Tại Việt Nam, chương trình chọn giống cá tra theo hướng tăng trưởng nhanh dựa trên phương pháp di truyền số lượng đã được triển khai từ năm 2001. Việc kết hợp di truyền số lượng với di truyền phân tử thông qua phương pháp chọn giống nhờ chỉ thị (Marker-assisted selection - MAS), đặc biệt là chỉ thị đa hình đơn nucleotide (Single Nucleotide Polymorphism - SNP), cho phép rút ngắn chu kỳ chọn lọc và nâng cao độ chính xác.

Hệ thống các nhân tố tăng trưởng tương tự insulin (Insulin-like Growth Factor - IGF) bao gồm các phối tử (IGF1, IGF2), các thụ thể màng (IGFR, tiêu biểu là IGF1R) và các protein liên kết IGF (IGFBP, gồm IGFBP-1 đến -7). Các phân tử này đóng vai trò trung tâm trong trục điều hòa nội tiết và cận tiết của quá trình phân bào, biệt hóa, và phì đại sợi cơ ở cá xương. Mặc dù các nghiên cứu tìm kiếm chỉ thị SNP trên các gen thuộc hệ thống IGF đã được thực hiện trên một số loài cá, dữ liệu đa hình phân tử trên các gen này ở cá tra nuôi (P. hypophthalmus) trước đây chưa từng được công bố. Luận án được thực hiện nhằm khai thác thông tin từ hệ gen tham chiếu đã giải mã để lấp đầy khoảng trống nghiên cứu này.

Mục tiêu nghiên cứu: Nghiên cứu phát hiện các đa hình SNP trên một số gen thuộc hệ thống IGF và kiểm nghiệm sự liên quan của các SNP này với tính trạng tăng trưởng ở cá tra nuôi, đề xuất được một số chỉ thị phân tử tiềm năng hướng tới ứng dụng trong chương trình chọn giống cá tra theo hướng tăng trưởng.

Nội dung nghiên cứu cụ thể:

  1. Phân tích cấu trúc các gen IGF1, IGF2, IGF1R, IGFBP-1, -2, -3, -5, -6, -7 thuộc hệ thống IGF của cá tra nuôi, xác định các vùng quan trọng để phát hiện đa hình SNP và giải trình tự các vùng đó trên 1 cá thể bằng phương pháp Sanger nhằm kiểm tra tính xác thực so với trình tự tham chiếu.
  2. Giải trình tự các vùng quan trọng thuộc 9 gen nêu trên trên bộ mẫu khởi tạo gồm 10 cá thể tăng trưởng nhanh và 10 cá thể tăng trưởng chậm bằng phương pháp Sanger, so sánh với trình tự tham chiếu để phát hiện và sàng lọc các SNP.
  3. Xác định kiểu gen (genotyping) các SNP được sàng lọc trên bộ mẫu kiểm nghiệm gồm 70 cá thể tăng trưởng nhanh và 70 cá thể tăng trưởng chậm bằng phương pháp kéo dài một nucleotide (Single Base Extension - SBE).
  4. Phân tích dữ liệu SNP thu được từ toàn bộ 160 cá thể (80 cá thể tăng trưởng nhanh và 80 cá thể tăng trưởng chậm), đánh giá mối liên quan giữa các SNP được sàng lọc, các kiểu gen đơn bội (haplotype) và tổ hợp kiểu gen với tính trạng tăng trưởng của cá tra nuôi.

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

  • Đối tượng nghiên cứu: Cá tra nuôi (Pangasianodon hypophthalmus (Sauvage, 1878)).
  • Vật liệu nghiên cứu: Mẫu mô vây của 160 cá thể cá tra thuộc thế hệ chọn giống thứ 3 (G3) cộng gộp (thu năm 2015) từ chương trình chọn giống của Viện Nghiên cứu Nuôi trồng Thủy sản 2. Mẫu được chia thành hai tập hợp: bộ mẫu khởi tạo (20 cá thể: 10 cá thể tăng trưởng nhanh và 10 cá thể tăng trưởng chậm) và bộ mẫu kiểm nghiệm (140 cá thể: 70 cá thể tăng trưởng nhanh và 70 cá thể tăng trưởng chậm). Mẫu mô được bảo quản trong cồn tuyệt đối ở nhiệt độ -20°C.
  • Phạm vi nội dung: 9 gen chức năng thuộc hệ thống IGF gồm IGF1, IGF2, IGF1R isoform X1, IGFBP-1, IGFBP-2A, IGFBP-3, IGFBP-5, IGFBP-6, IGFBP-7.
  • Thời gian và địa điểm thực hiện: Nghiên cứu được hoàn thành và bảo vệ năm 2023 tại Viện Nghiên cứu hệ gen, trực thuộc Học viện Khoa học và Công nghệ – Viện Hàn lâm Khoa học và Công nghệ Việt Nam.

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

Tổng quan tài liệu của luận án hệ thống hóa các nghiên cứu trong nước và quốc tế về di truyền số lượng, di truyền phân tử và hệ thống tín hiệu IGF ở động vật thủy sản:

1. Các nghiên cứu chọn giống và di truyền trên cá tra tại Việt Nam:

  • Nghiên cứu chỉ thị phân tử bước đầu: Phạm Anh Tuấn và Nguyễn Hữu Ninh (2003) sử dụng 4 chỉ thị microsatellite để tìm kiếm mối tương quan với màu sắc thịt cá tra. Nguyễn Minh Thành và cộng sự (2015) nghiên cứu hệ gen chức năng ở mô thận nhằm xác định gen liên quan đến tính chịu mặn. Các nghiên cứu khác ứng dụng kỹ thuật RAPD, AFLP và microsatellite để đánh giá tính đa dạng di truyền quần đàn cá tra. Marnis và cộng sự (2018) kiểm nghiệm 5 chỉ thị microsatellite và phát hiện 3 chỉ thị liên quan đến chiều dài và khối lượng thân.
  • Chương trình chọn giống di truyền số lượng của Viện Nghiên cứu Nuôi trồng Thủy sản 2 (2001–2022): Trải qua các giai đoạn: chọn lọc cá thể theo tính trạng tăng trưởng (2001–2005); chọn lọc gia đình nâng cao tỷ lệ philê và tăng trưởng (2006–2008); đánh giá hiệu quả chọn giống (2010–2012); ứng dụng di truyền phân tử kết hợp di truyền số lượng (2013–2016); và chọn giống thế hệ G4 (2019–2022). Nghiên cứu xác định hệ số di truyền ước tính cho tính trạng tăng trưởng ở thế hệ G2 đạt mức 0,24–0,30; ở thế hệ G3 đạt 0,34–0,52 (hệ số di truyền thực tế 0,24–0,38). Đàn G3 chọn lọc gồm 1.230 cá thể thuộc 157 gia đình đạt hiệu quả chọn lọc thực tế cao hơn 20,4% so với đàn chưa qua chọn lọc.
  • Dữ liệu hệ gen tham chiếu cá tra: Kim Thị Phương Oanh và cộng sự (2018) giải mã và lắp ráp toàn bộ hệ gen cá tra nuôi bằng hệ thống Illumina với kích thước ước tính khoảng 700 Mbp, chú giải 28.600 gen mã hóa protein (mã số BioProject PRJNA448819, mã lắp ráp QUXB00000000). Nhóm tác giả đồng thời phân tích transcriptome mô cơ và mô não bằng RNA-seq giữa hai nhóm cá tăng trưởng nhanh và chậm.

2. Các nghiên cứu đa hình gen liên quan đến tăng trưởng ở cá xương trên thế giới:

  • Gen GH và GHRH: Trên cá quế (Siniperca chuatsi), phát hiện 4 SNP trên gen GH (g.4948A>T, g.5045T>C...), trong đó 3 SNP tác động có ý nghĩa đến tính trạng tăng trưởng; nghiên cứu tiếp theo xác định 32 chỉ thị SNP/SSR hình thành 4 haplotype liên quan chặt chẽ đến sinh trưởng. Trên cá hồi Bắc cực (Salvelinus alpinus), phát hiện tương quan giữa SNP trên gen GHRH/PACAP2 với tốc độ tăng trưởng giai đoạn sớm.
  • Gen Myostatin (MTTN): Ở cá chép (Cyprinus carpio), phát hiện 2 SNP tại exon 3 (c.72C>T) trên gen MTTN liên quan đến khối lượng cơ thể. Ở cá ngão (Ancherythroculter nigrocauda), SNP g.1289G>A trên gen MTTN-1 liên quan đến tăng trưởng nhanh.
  • Các gen thuộc hệ thống IGF:
    • Cá hồi Đại Tây Dương (Salmo salar): Phát hiện 3 SNP trên gen IGF1 (g.5763G>T vùng promoter, g.7292C>T intron 1, g.4671A>C intron 3); trong đó 2 SNP liên quan chặt chẽ đến khối lượng cá thể trên quần thể 4.800 mẫu. Haplotype phát sinh từ 3 SNP liên quan đến khối lượng thu hoạch và tỷ lệ philê.
    • Cá chép (Cyprinus carpio): SNP g.7627T>A tại intron 2 của gen IGF1 liên quan đến khối lượng và chiều dài thân; cá thể mang kiểu gen AA có khối lượng trung bình cao hơn kiểu gen TT 5,9%.
    • Cá vược châu Âu (Dicentrarchus labrax): SNP g.46749C>T tại vùng 5'-UTR của gen IGF1 và SNP tại intron 3 của gen IGF2 liên quan có ý nghĩa đến chiều dài và khối lượng cơ thể.
    • Cá rô phi sông Nin (Oreochromis niloticus): Phát hiện 13 SNP trên gen IGF2, trong đó SNP G161A tại exon 3 liên quan đến tính trạng tăng trưởng.
    • Cá bản địa Trung Quốc (Odontobutis potamophila): Khảo sát 3 SNP trên intron của gen IGF1R (1777T>C, 1704G>A, 1208G>A), xác định kiểu gen dị hợp AG ở vị trí 1208G>A có ưu thế vượt trội về chiều dài, chiều rộng và khối lượng thân so với kiểu gen GG.

Khoảng trống nghiên cứu luận án giải quyết: Mặc dù hệ gen tham chiếu của cá tra đã được công bố, việc xác định các biến dị nucleotide đơn (SNP) trên các gen thuộc hệ thống IGF (IGF1, IGF2, IGF1R, IGFBP-1 đến -7) và phân tích mối liên kết giữa các SNP, haplotype, tổ hợp kiểu gen với tính trạng tăng trưởng (thông qua chỉ số giá trị chọn giống ước đoán - EBV và khối lượng thực tế) chưa từng được thực hiện trên đối tượng cá tra nuôi. Luận án định vị nghiên cứu vào việc khai thác trực tiếp nguồn dữ liệu hệ gen để cung cấp các chỉ thị phân tử phục vụ công tác chọn giống cá tra.

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

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

  • Cơ chế phát triển cơ ở cá xương: Tăng trưởng cơ sau giai đoạn phôi bắt đầu từ tế bào gốc chuyển hóa thành tế bào tiền thân sinh cơ (MPCs), tăng sinh, di cư và biệt hóa thành nguyên bào cơ (committed myoblasts). Quá trình diễn ra theo hai hướng: hợp nhất tạo sợi cơ mới (hyperplasia - tăng sản) hoặc sáp nhập vào sợi cơ hiện có (hypertrophy - phì đại, bồi tụ hạt nhân).
  • Hệ trục nội tiết điều hòa tăng trưởng: Tín hiệu môi trường kích hoạt vùng dưới đồi tiết GHRH (kích thích) hoặc GHIH (ức chế), điều hòa thùy trước tuyến yên tiết GH vào máu. GH liên kết thụ thể GHR ở mô gan kích thích sản sinh IGF1 và IGF2. Các phối tử IGF liên kết với protein vận chuyển IGFBP để tuần hoàn trong máu và tiếp cận thụ thể IGF1R trên bề mặt tế bào cơ.
  • Cấu trúc và con đường truyền tín hiệu hệ thống IGF:
    • Phối tử IGF1 và IGF2: Cấu trúc tiền thân gồm 5 domain (B-C-A-D-E). Sau khi cắt bỏ chuỗi peptit tín hiệu và domain E, phân tử trưởng thành gồm 4 domain B-C-A-D. IGF1 ở bộ cá nheo và cá chép bảo tồn thêm 2 axit amin (His, Asn) tại vị trí 39 và 40 của domain C.
    • Thụ thể IGF1R: Cấu trúc tetramer gồm hai chuỗi dị dimer $\alpha\beta$ liên kết bằng cầu nối disunfua. Tiểu phần $\beta$ chứa domain xuyên màng và domain tyrosine kinase nội bào mang motif bảo thủ bám ATP (GXGXXXG21XK), cụm 3 tyrosine tự phosphoryl hóa (YXXXYY), và motif NPEY liên kết cơ chất thụ thể insulin (IRS-1), hoạt hóa con đường PI3K/AKT/mTOR và MAPK.
    • Họ protein IGFBP: Gồm 2 domain chức năng bảo thủ: domain bám IGF ở đầu N (chứa 12 phân tử Cys bảo thủ và motif GCGCCXXC) và domain thyroglobulin type 1 ở đầu C (chứa 6 Cys và motif CWCV). Riêng IGFBP-7 (IGFBP-rP1) có ái lực yếu với IGF nhưng ái lực cao với insulin, tham gia điều hòa trục tín hiệu IIS.

Nguồn vật liệu và thiết kế thí nghiệm

Nghiên cứu sử dụng 160 mẫu cá tra thuộc thế hệ G3 chọn giống năm 2015 của Viện Nghiên cứu Nuôi trồng Thủy sản 2. Mẫu vây được phân tích theo 2 nhóm hình thái tăng trưởng dựa trên giá trị chọn giống ước đoán (EBV) cá thể và EBV gia đình:

Bộ mẫu Nhóm mẫu Ký hiệu Số lượng Khối lượng trung bình (g) EBV cá thể trung bình EBV gia đình trung bình Nguồn gia đình
Khởi tạo Tăng trưởng nhanh N1 – N10 10 $1.746,9 \pm 272,0$ $359,9 \pm 53,1$ $210,1 \pm 33,9$ 9 gia đình có EBV cao nhất
Tăng trưởng chậm C1 – C10 10 $512,2 \pm 148,8$ $-163,6 \pm 40,6$ $-78,0 \pm 49,5$ 9 gia đình có EBV thấp nhất
Kiểm nghiệm Tăng trưởng nhanh — 70 Phụ lục 1 Xếp hạng EBV cao nhất Xếp hạng EBV cao nhất 24 gia đình có EBV cao nhất
Tăng trưởng chậm — 70 Phụ lục 1 Xếp hạng EBV thấp nhất Xếp hạng EBV thấp nhất 31 gia đình có EBV thấp nhất
Tổng cộng — — 160 — — — —

Phương pháp phân tích thực nghiệm và công cụ kỹ thuật

  1. Phân tích tin sinh học in silico: Sử dụng ProteinBLAST (blastp) đối chiếu cơ sở dữ liệu Protein Data Bank (PDB) để xác định các domain bảo thủ của 9 gen đích; phần mềm SignalP 5.0 dự đoán chuỗi peptit tín hiệu (Signal peptide).
  2. Tách chiết DNA tổng số: Mô vây được nghiền trong nitơ lỏng, ủ trong đệm ly giải ($0,01\text{ M EDTA}$, $0,01\text{ M Tris-HCl pH } 7,5$, $0,1\text{ M NaCl}$, $2,1%\text{ SDS}$, $100\ \mu\text{g/ml Proteinase K}$) ở $56^\circ\text{C}$ trong 3 giờ. Tách chiết bằng phenol/chloroform, kiểm tra độ tinh sạch và nồng độ bằng NanoDrop One Spectrophotometer và điện di gel agarose 1%.
  3. Thiết kế mồi và khuếch đại PCR: Thiết kế 55 cặp mồi bằng phần mềm Primer 3 (v.0.4.0) dựa trên hệ gen tham chiếu cá tra. Phản ứng PCR thể tích $25\ \mu\text{l}$ gồm $1\ \mu\text{l}$ DNA khuôn, $1\ \mu\text{l}$ mỗi mồi ($10\text{ pmol/}\mu\text{l}$), $12,5\ \mu\text{l}$ Taq 2X Master Mix (NEB). Chu trình nhiệt: biến tính ban đầu $95^\circ\text{C}$ (3 phút); 30 chu kỳ ($95^\circ\text{C}$ trong 30 giây, nhiệt độ gắn mồi tối ưu $T_a$, $68^\circ\text{C}$ trong 1 phút); kéo dài cuối cùng $68^\circ\text{C}$ (7 phút). Sản phẩm được tinh sạch bằng GeneJET PCR Purification Kit (Thermo Fisher Scientific).
  4. Giải trình tự Sanger: Giải mã trình tự các đoạn gen khuếch đại trên máy giải trình tự tự động ABI để kiểm tra tính xác thực so với trình tự tham chiếu và phát hiện các vị trí SNP trên 20 mẫu bộ khởi tạo.
  5. Kỹ thuật kéo dài một nucleotide (Single Base Extension - SBE): Sử dụng bộ kit thương mại SNaPshot Multiplex Kit (Applied Biosystems) với 4 loại ddNTP gắn huỳnh quang đặc trưng: dATP (dR6G - xanh lá), dTTP (dROX - đỏ), dCTP (dTAMRA - đen), dGTP (dR110 - xanh dương). Thiết kế các mồi SBE có phần đuôi 5' chênh lệch chiều dài từ 5 đến 8 nucleotide nhằm phân tách các đỉnh tín hiệu khi điện di mao quản. Dữ liệu kiểu gen được phân tích bằng phần mềm GeneMapper (ABI).
  6. Phương pháp phân tích di truyền thống kê: Tính toán tần số alen, tần số kiểu gen, kiểm định cân bằng Hardy-Weinberg ($p > 0,05$), phân tích mất cân bằng liên kết (Linkage Disequilibrium - LD qua các chỉ số $D'$ và $R^2$), phân tích haplotype và tổ hợp kiểu gen liên kết với khối lượng cơ thể và giá trị EBV trên toàn bộ 160 cá thể.

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

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

Chương 1 tổng thuật chi tiết về phân loại học, sinh học cá tra nuôi (Pangasianodon hypophthalmus), tình hình phát triển kinh tế theo Quyết định 50/2018/QĐ-TTg và kết quả xuất khẩu thủy sản đạt trên 1,61 tỷ USD năm 2021. Tác giả phân tích tiến trình chọn giống di truyền số lượng qua 4 thế hệ (G0 đến G4) tại Viện Nghiên cứu Nuôi trồng Thủy sản 2, làm rõ các thông số di truyền của tính trạng tăng trưởng (hệ số di truyền thế hệ G3 đạt 0,51 ở đàn cộng gộp; hiệu quả chọn lọc thực tế đạt 20,4%).

Chương này trình bày cơ sở sinh học phân tử của quá trình phát triển cơ ở cá xương, mô hình điều hòa trục GHRH/GHIH - GH - IGF, cấu trúc miền chức năng của các protein thuộc hệ thống IGF (IGF1, IGF2, thụ thể IGF1R, và các IGFBP-1 đến -7). Tác giả so sánh nguyên lý, ưu điểm và nhược điểm của 5 phương pháp định kiểu gen SNP (SNP genotyping): giải mã lại bằng Sanger, lai phân tử (SNP array), PCR đặc hiệu alen dựa trên FRET (ASQ), phương pháp kéo dài một nucleotide (SBE/SNaPshot), và phương pháp cắt bằng enzyme giới hạn (Invader assay). Đồng thời, chương 1 thống kê danh mục 15 gen thuộc hệ thống IGF đã được nhận diện trên hệ gen cá tra nuôi:

STT Tên gen Gene ID Vị trí Scaffold trên hệ gen Chiều dài Nucleotide Chiều dài Protein (aa)
1 IGF1 113532358 sc0000011 (phy_g435.t1) 736 245
2 IGF2 113540971 sc0000028 (phy_g25602.t1) 639 213
3 IGF1R isoform X1 113534424 sc0000015 (phy_g24431.t1) 4237 1412
4 IGFBP-1 113545355 sc0000042 (phy_g24144.t1) 796 265
5 IGFBP-2A 113545824 sc0000046 (phy_g8897.t1) 786 262
6 IGFBP-3 113545296 sc0000042 (phy_g24145.t1) 906 302
7 IGFBP-5 113525861 sc0000003 (phy_g8120.t1) 919 306
8 IGFBP-6 113530376 sc0000009 (phy_g1604.t1) 612 204
9 IGFBP-7 113546920 sc0000054 (phy_g24954.t1) 792 264

Chương 2: Nguyên vật liệu và phương pháp

Chương 2 mô tả chi tiết quy trình thu mẫu 160 cá thể cá tra thế hệ G3 từ Viện Nghiên cứu Nuôi trồng Thủy sản 2. Tác giả trình bày phương pháp luận phân chia mẫu thành 2 giai đoạn:

  • Bộ mẫu khởi tạo (20 cá thể: 10 tăng trưởng nhanh ký hiệu N1–N10, khối lượng trung bình $1.746,9 \pm 272,0\text{ g}$, EBV trung bình $359,9 \pm 53,1$; 10 tăng trưởng chậm ký hiệu C1–C10, khối lượng trung bình $512,2 \pm 148,8\text{ g}$, EBV trung bình $-163,6 \pm 40,6$).
  • Bộ mẫu kiểm nghiệm (140 cá thể: 70 nhanh thuộc 24 gia đình có EBV cao nhất và 70 chậm thuộc 31 gia đình có EBV thấp nhất).

Chương 2 chuẩn hóa toàn bộ quy trình thực nghiệm: ly giải mô vây và tách chiết DNA bằng phenol/chloroform; quy cách thiết kế 55 cặp mồi PCR đặc hiệu cho 9 gen đích bằng phần mềm Primer 3; chu trình nhiệt PCR dùng Taq 2X Master Mix (NEB); quy trình giải trình tự trực tiếp bằng phương pháp Sanger; quy trình phản ứng SBE với SNaPshot Multiplex Kit; thiết kế các mồi SBE đa phân đoạn; và quy trình xử lý tín hiệu huỳnh quang mao quản trên phần mềm GeneMapper kết hợp kiểm định di truyền số lượng.

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

Chương 3 trình bày các kết quả thực nghiệm và thảo luận học thuật:

  1. Phân tích cấu trúc 9 gen thuộc hệ thống IGF của cá tra:
    • Gen IGF1 và IGF2 bảo tồn cấu trúc 5 domain B-C-A-D-E ở dạng tiền thân, trong đó IGF1 cá tra nuôi chứa vị trí chèn đặc thù của bộ cá nheo.
    • Gen IGF1R mã hóa chuỗi polypeptide 1.412 axit amin, bảo tồn hoàn toàn domain tyrosine kinase nội bào, bao gồm motif liên kết ATP, cụm 3 tyrosine tự phosphoryl hóa và motif NPEY liên kết IRS-1.
    • Các gen IGFBP-1, -2A, -3, -5, -6 bảo tồn domain đầu N giàu cysteine và domain thyroglobulin type 1 ở đầu C. IGFBP-7 chứa domain IGFBP biến đổi có tính tương đồng cao với các loài cá xương khác.
  2. Kiểm tra tính xác thực của kết quả giải mã Sanger: Trình tự khuếch đại trên 1 cá thể kiểm tra cho độ tương đồng hoàn toàn với trình tự hệ gen tham chiếu công bố trên NCBI (BioProject PRJNA448819).
  3. Phát hiện và sàng lọc SNP trên bộ mẫu khởi tạo (20 cá thể): Giải trình tự Sanger 9 gen đích phát hiện hàng loạt các vị trí đa hình đơn nucleotide trên cả vùng exon và intron. Dựa trên tiêu chí sàng lọc (tần số alen thiểu số MAF, sự khác biệt phân bố kiểu gen giữa nhóm nhanh và nhóm chậm), luận án đã sàng lọc được các SNP mục tiêu trên từng gen.
  4. Định kiểu gen SBE trên bộ mẫu kiểm nghiệm (140 cá thể) và phân tích liên kết:
    • Kỹ thuật SNaPshot multiplex và phần mềm GeneMapper xác định chính xác kiểu gen của các mẫu kiểm nghiệm.
    • Tổng hợp phân tích trên 160 cá thể khẳng định 6 chỉ thị SNP tiềm năng có mối liên quan chặt chẽ với tính trạng tăng trưởng (khối lượng cơ thể và chỉ số EBV) của cá tra nuôi:
      • Gen IGF1: SNP $13680\text{ A>T}$ (IGF1.13680 A>T).
      • Gen IGF1R: Cụm 3 SNP (gồm IGF1R.83894 A>G và 2 SNP liên kết) hình thành các haplotype có tác động rõ rệt đến giá trị chọn giống.
      • Gen IGFBP-3: SNP $704\text{ C>G}$ (gây đột biến thay thế axit amin p.Leu8Val).
      • Gen IGFBP-5: SNP $525\text{ T>A}$ (gây đột biến thay thế axit amin p.Val16Glu).
      • Gen IGFBP-6: SNP $2278\text{ C>A}$.
      • Gen IGFBP-7: Cụm 3 SNP (tiêu biểu là IGFBP-7.4559 C>A) tạo thành haplotype liên quan đến tính trạng sinh trưởng.
    • Phân tích mất cân bằng liên kết (LD) qua hệ số $D'$ và $R^2$ xác lập được các khối liên kết di truyền chặt chẽ trên gen IGF1R và IGFBP-7. Đánh giá tác động cộng gộp của các haplotype và tổ hợp kiểu gen phát sinh từ 6 chỉ thị SNP cho thấy các cá thể mang tổ hợp alen ưu thế có khối lượng và chỉ số EBV vượt trội có ý nghĩa thống kê so với các tổ hợp kiểu gen còn lại.

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

Những đóng góp mới về mặt khoa học và lý luận:

  1. Cung cấp dữ liệu đầy đủ về đặc điểm cấu trúc gen và chuỗi protein suy diễn của 9 gen thuộc hệ thống IGF (IGF1, IGF2, IGF1R, IGFBP-1, IGFBP-2A, IGFBP-3, IGFBP-5, IGFBP-6, IGFBP-7) ở cá tra nuôi (Pangasianodon hypophthalmus), chứng minh mức độ bảo tồn cao của các domain chức năng so với các loài cá xương khác.
  2. Xác định bức tranh đa dạng di truyền nucleotide trên các gen thuộc hệ thống IGF ở cá tra, làm rõ mức độ mất cân bằng liên kết (LD) và cấu trúc haplotype của các locus gen chức năng này.

Những đóng góp mới về mặt thực tiễn:

  1. Phát hiện, phân tích và sàng lọc thành công 6 chỉ thị SNP tiềm năng trên các gen thuộc hệ thống IGF có liên quan chặt chẽ đến tính trạng tăng trưởng (khối lượng và giá trị chọn giống ước đoán EBV) ở cá tra nuôi, gồm:
    • SNP $13680\text{ A>T}$ trên gen IGF1.
    • Nhóm SNP trên gen IGF1R (tiêu biểu là SNP $83894\text{ A>G}$).
    • SNP $704\text{ C>G}$ (p.Leu8Val) trên gen IGFBP-3.
    • SNP $525\text{ T>A}$ (p.Val16Glu) trên gen IGFBP-5.
    • SNP $2278\text{ C>A}$ trên gen IGFBP-6.
    • Nhóm SNP trên gen IGFBP-7 (tiêu biểu là SNP $4559\text{ C>A}$).
  2. Thiết lập bộ mồi PCR (55 cặp mồi) và hệ thống mồi SBE đa kênh (multiplex SNaPshot) được tối ưu hóa để định kiểu gen nhanh các SNP đích trên đàn cá tra.
  3. Xác định các haplotype và tổ hợp kiểu gen ưu thế có liên quan đến khả năng tăng trưởng vượt trội, làm cơ sở khoa học để phát triển chỉ thị phân tử hỗ trợ chọn giống (MAS) nhằm nâng cao tốc độ sinh trưởng ở cá tra nuôi.

Kiến nghị của luận án:

  • Ứng dụng 6 chỉ thị phân tử SNP đã phát hiện vào công tác đánh giá và chọn lọc đàn cá bố mẹ trong các chương trình chọn giống cá tra tăng trưởng nhanh.
  • Tiếp tục mở rộng nghiên cứu kiểm chứng các chỉ thị SNP này trên các quần đàn cá tra thương phẩm và các đối tượng cá tra sinh sản tự nhiên ở các vùng sinh thái khác nhau.

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

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

  • Nghiên cứu tập trung đánh giá trên 160 cá thể thuộc thế hệ chọn giống thứ ba (G3) của một quần thể chọn giống tại Viện Nghiên cứu Nuôi trồng Thủy sản 2, chưa bao quát toàn bộ các dòng/quần đàn cá tra phân bố ở các vùng sinh thái khác.
  • Phạm vi khảo sát tập trung vào 9 gen đại diện trong hệ thống IGF; chưa khảo sát mở rộng sang toàn bộ các gen khác thuộc trục hormone sinh trưởng (như GHR, GHRH) hay các yếu tố điều hòa cơ sinh dưỡng (MyoD, MyoG, myf-5).

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

  • Kiểm nghiệm quy mô lớn các chỉ thị SNP, haplotype và tổ hợp kiểu gen đã xác định trên các quần thể cá tra tự nhiên và các thế hệ cá tra chọn giống tiếp theo (G4, G5).
  • Tích hợp các chỉ thị SNP tiềm năng này vào các tấm thẻ định kiểu gen (SNP array) hoặc panel đa gen để ứng dụng chọn giống dựa trên toàn hệ gen (Genomic Selection - GS).
  • Nghiên cứu cơ chế phân tử chức năng của các đột biến thay thế axit amin (như p.Leu8Val trên IGFBP-3 và p.Val16Glu trên IGFBP-5) đối với hoạt tính sinh học của protein tương ứng.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các đối tượng học thuật và thực tiễn:

  • Nghiên cứu sinh, học viên cao học và nhà khoa học chuyên ngành Công nghệ sinh học, Di truyền học, Thủy sản: Cung cấp phương pháp luận chuẩn xác về quy trình sàng lọc SNP từ hệ gen tham chiếu, kỹ thuật thiết kế thí nghiệm từ bộ mẫu khởi tạo (Sanger) đến bộ mẫu kiểm nghiệm mở rộng (SBE/SNaPshot), và các thuật toán phân tích di truyền liên kết (LD, Haplotype, EBV).
  • Cán bộ kỹ thuật và chuyên gia chọn giống thủy sản: Tham khảo trực tiếp danh mục 6 chỉ thị SNP tiềm năng và các tổ hợp kiểu gen ưu thế để tích hợp vào các chương trình chọn giống nhờ chỉ thị phân tử (MAS), giúp rút ngắn thời gian và chi phí chọn lọc đàn cá tra tăng trưởng nhanh.
  • Giảng viên đại học: Sử dụng làm tài liệu giảng dạy cho các học phần Di truyền số lượng, Công nghệ sinh học thủy sản, Hệ gen học chức năng và Đánh giá đa dạng di truyền phân tử.

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

1. Luận án đã lựa chọn những gen nào thuộc hệ thống IGF của cá tra để phân tích đa hình SNP?
Luận án lựa chọn 9 gen thuộc hệ thống IGF của cá tra nuôi (Pangasianodon hypophthalmus) dựa trên hệ gen tham chiếu đã giải mã, bao gồm: 2 gen phối tử (IGF1, IGF2), 1 gen thụ thể (IGF1R isoform X1), và 6 gen mã hóa protein bám (IGFBP-1, IGFBP-2A, IGFBP-3, IGFBP-5, IGFBP-6, IGFBP-7).

2. Thiết kế bộ mẫu nghiên cứu của luận án gồm bao nhiêu cá thể và dựa trên tiêu chí nào?
Nghiên cứu sử dụng tổng cộng 160 cá thể cá tra thế hệ G3 (năm 2015) từ Viện Nghiên cứu Nuôi trồng Thủy sản 2, chia làm 2 bộ mẫu dựa trên giá trị chọn giống ước đoán (EBV):

  • Bộ mẫu khởi tạo (20 cá thể): 10 cá thể tăng trưởng nhanh (khối lượng trung bình $1.746,9\text{ g}$, EBV trung bình $359,9$) và 10 cá thể tăng trưởng chậm (khối lượng trung bình $512,2\text{ g}$, EBV trung bình $-163,6$) thuộc các gia đình có EBV tương ứng cao nhất và thấp nhất.
  • Bộ mẫu kiểm nghiệm (140 cá thể): 70 cá thể tăng trưởng nhanh (từ 24 gia đình có EBV cao nhất) và 70 cá thể tăng trưởng chậm (từ 31 gia đình có EBV thấp nhất).

3. Phương pháp kỹ thuật nào được sử dụng để phát hiện và định kiểu gen SNP trong luận án?
Luận án kết hợp hai phương pháp chính:

  • Giải trình tự Sanger trực tiếp: Áp dụng trên bộ mẫu khởi tạo (20 cá thể) để kiểm tra tính xác thực trình tự so với hệ gen tham chiếu và phát hiện các vị trí SNP tiềm năng trên 9 gen đích.
  • Kéo dài một nucleotide (SBE) bằng SNaPshot Multiplex Kit: Kết hợp mồi SBE đa phân đoạn và ddNTP huỳnh quang, phân tích sản phẩm qua điện di mao quản trên hệ máy tự động ABI bằng phần mềm GeneMapper để định kiểu gen trên bộ mẫu kiểm nghiệm (140 cá thể).

4. Luận án đã xác định được bao nhiêu chỉ thị SNP có liên quan đến tính trạng tăng trưởng ở cá tra?
Luận án đã phát hiện, sàng lọc và kiểm nghiệm được 6 chỉ thị SNP tiềm năng liên quan có ý nghĩa đến tính trạng tăng trưởng (khối lượng và EBV) ở cá tra nuôi: SNP $13680\text{ A>T}$ (IGF1), cụm SNP trên gen IGF1R (tiêu biểu là SNP $83894\text{ A>G}$), SNP $704\text{ C>G}$ (IGFBP-3), SNP $525\text{ T>A}$ (IGFBP-5), SNP $2278\text{ C>A}$ (IGFBP-6), và cụm SNP trên gen IGFBP-7 (tiêu biểu là SNP $4559\text{ C>A}$).

5. Các phân tích haplotype và tổ hợp kiểu gen trong luận án mang lại kết quả gì?
Phân tích mất cân bằng liên kết (LD) xác định được các khối liên kết di truyền chặt chẽ trên các gen IGF1R và IGFBP-7. Đánh giá tương quan cho thấy các haplotype và tổ hợp kiểu gen phát sinh từ 6 chỉ thị SNP được lựa chọn có tác động cộng gộp rõ rệt, trong đó các cá thể mang tổ hợp kiểu gen ưu thế đạt khối lượng và chỉ số EBV vượt trội so với các kiểu gen đối chứng.

Kết luận

Luận án của nghiên cứu sinh Trần Thị Huyền Trang đã hoàn thành việc phân tích cấu trúc, phát hiện và kiểm nghiệm đa hình đơn nucleotide (SNP) trên 9 gen thuộc hệ thống IGF ở cá tra nuôi (Pangasianodon hypophthalmus). Nghiên cứu đã sàng lọc được 6 chỉ thị SNP tiềm năng trên các gen IGF1, IGF1R, IGFBP-3, IGFBP-5, IGFBP-6, IGFBP-7 cùng các haplotype và tổ hợp kiểu gen có mối liên quan chặt chẽ với tính trạng tăng trưởng và giá trị chọn giống ước đoán EBV. Kết quả đạt được cung cấp cơ sở dữ liệu di truyền phân tử có giá trị thực tiễn cao, hỗ trợ trực tiếp cho các chương trình chọn giống cá tra theo hướng tăng trưởng nhanh thông qua chọn giống nhờ chỉ thị phân tử (MAS).