Tổng quan luận án

Nghề nuôi cá tra (Pangasianodon hypophthalmus) đóng vai trò quan trọng trong ngành thủy sản Việt Nam. Theo số liệu của Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO, 2020), Việt Nam đứng thứ tư thế giới về tổng sản lượng thủy sản với 4,134 triệu tấn; trong đó, sản lượng nuôi cá tra đạt 1,56 triệu tấn và kim ngạch xuất khẩu đạt 1,61 tỉ USD vào năm 2021 (VASEP, 2021). Diện tích ương nuôi cá tra giống tại vùng Đồng bằng sông Cửu Long đạt gần 4.000 ha. Tuy nhiên, việc phát triển thâm canh khiến nghề nuôi đối mặt với nguy cơ dịch bệnh bùng phát, nổi bật là bệnh gan thận mủ do vi khuẩn Edwardsiella ictaluri gây ra (Ferguson và ctv., 2001; Le và Cheong, 2010). Bệnh có thể gây tỷ lệ chết lên đến 90% nếu không xử lý kịp thời (Nguyễn Thị Thúy Liễu và ctv., 2011; Từ Thanh Dung và ctv., 2010), xuất hiện phổ biến nhất ở cá dưới 4 tháng tuổi gồm giai đoạn cá hương 21–30 ngày tuổi (chiếm tỷ lệ nhiễm 46%) và giai đoạn cá giống 40–90 ngày tuổi (chiếm tỷ lệ nhiễm 30%) (Tran và ctv., 2020). Các giải pháp như sử dụng kháng sinh tiềm ẩn rủi ro kháng thuốc và dư lượng vượt ngưỡng xuất khẩu (Từ Thanh Dung và ctv., 2010); trong khi chế phẩm kích thích miễn dịch hay vắc-xin đòi hỏi chi phí cao và hiệu lực ngắn hạn (Phu và ctv., 2016).

Chọn giống kháng bệnh theo hướng tiếp cận di truyền số lượng được xem là giải pháp bền vững. Viện Nghiên cứu Nuôi trồng Thủy sản II đã tạo lập đàn cá tra bố mẹ chọn giống kháng bệnh gan thận mủ thế hệ đầu tiên (G0) với hệ số di truyền ước tính $h^2 = 0,27$ và hiệu quả chọn lọc đạt 8,3% (Trịnh Quốc Trọng và ctv., 2017). Tuy nhiên, thực tiễn nghiên cứu đặt ra các khoảng trống cần giải quyết:

  • Cần kiểm chứng hệ số di truyền ở thế hệ con G1 và xác định các thông số di truyền ở giai đoạn cá hương, đồng thời đánh giá tương quan di truyền giữa tính trạng kháng bệnh giai đoạn cá hương, cá giống với tính trạng tăng trưởng khi thu hoạch.
  • Phương pháp đánh dấu vật lý bằng thẻ từ PIT (Passive Integrated Transponder) đòi hỏi phải ương nuôi riêng rẽ các gia đình đến kích thước 15–20 g trong 3–4 tháng, làm gia tăng ảnh hưởng của môi trường ương riêng rẽ ($c^2$) và tốn kém cơ sở hạ tầng. Cần phát triển bộ chỉ thị phân tử microsatellite có độ chính xác cao để thay thế thẻ PIT trong truy xuất phả hệ.
  • Việc đánh giá tính trạng kháng bệnh chủ yếu dựa vào biến nhị phân sống/chết tại thời điểm kết thúc cảm nhiễm chưa phản ánh đầy đủ trạng thái đáp ứng miễn dịch trong suốt quá trình nhiễm bệnh.

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

  • Mục tiêu chung: Nghiên cứu ước tính các thông số di truyền tính trạng kháng bệnh gan thận mủ và đề xuất phương án chọn lọc cá tra (Pangasianodon hypophthalmus) kháng bệnh gan thận mủ thế hệ G1, đồng thời phát triển các kỹ thuật nhằm tối ưu hóa chương trình chọn giống cá tra kháng bệnh gan thận mủ dài hạn.
  • Mục tiêu cụ thể:
    1. Nghiên cứu ước tính các thông số di truyền tính trạng kháng bệnh gan thận mủ và đề xuất phương án chọn lọc cá tra thế hệ G1.
    2. Nghiên cứu các giải pháp kỹ thuật nhằm áp dụng cho chọn lọc kháng bệnh gan thận mủ nhằm tăng hiệu quả cho chương trình chọn giống cá tra kháng bệnh trong tương lai.

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

  • Đối tượng nghiên cứu: Cá tra (Pangasianodon hypophthalmus) thuộc thế hệ bố mẹ G0 và thế hệ con G1; vi khuẩn Edwardsiella ictaluri chủng Gly09M; hệ thống 10 chỉ thị microsatellite đặc hiệu.
  • Phạm vi không gian: Thực hiện tại Trại Nghiên cứu Thực nghiệm Thủy sản Gò Vấp, Trung tâm Quan trắc Môi trường và Bệnh Nam Bộ, Trung tâm Quốc gia Giống Thủy sản Nước ngọt Nam Bộ, Trung tâm Công nghệ Sinh học thuộc Viện Nghiên cứu Nuôi trồng Thủy sản II; phân tích kỹ thuật di truyền phân tử tại Phòng thí nghiệm Đa dạng Di truyền và Tiến hóa (Khoa Sinh học, Đại học KU Leuven, Vương quốc Bỉ).
  • Phạm vi thời gian: Đề tài hoàn thành và bảo vệ luận án tiến sĩ vào tháng 12 năm 2022.

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 kết quả nghiên cứu trong và ngoài nước về chọn giống thủy sản, phương pháp gây bệnh thực nghiệm, ứng dụng chỉ thị phân tử và các chỉ tiêu đáp ứng miễn dịch:

Nhóm nội dung Tác giả và năm công bố Phát hiện chính được ghi nhận trong văn bản
Bệnh gan thận mủ trên cá tra Ferguson và ctv. (2001); Crumlish và ctv. (2002); Lý Thị Thanh Loan (2009); Tran và ctv. (2020) Bệnh do vi khuẩn E. ictaluri gây ra; tỷ lệ nhiễm tập trung cao nhất ở cá dưới 4 tháng tuổi (giai đoạn cá hương 46%, cá giống 30%); không ghi nhận bệnh ở cá trên 900 g.
Chọn giống thủy sản thế giới Gjedrem (2005); Kjøglum và ctv. (2008); Yáñez và ctv. (2016); Shoemaker và ctv. (2017) Hơn 101 chương trình chọn giống từ năm 1983; hệ số di truyền tính trạng kháng bệnh dao động từ trung bình đến cao ($h^2 = 0,15 - 0,62$), cho phép cải thiện di truyền qua các thế hệ.
Chọn giống cá tra tại Việt Nam Nguyễn Văn Sáng và ctv. (2009, 2012); Trịnh Quốc Trọng và ctv. (2017); Pham và ctv. (2020a, 2020b) Viện NCNTTS II chọn giống tăng trưởng qua 3 thế hệ; thành lập quần thể G0 kháng bệnh gan thận mủ ($h^2 = 0,27 - 0,37$, hiệu quả chọn lọc 8,3%, tương quan kháng bệnh và tăng trưởng là 0,26).
Truy xuất phả hệ bằng microsatellite Bùi Thị Liên Hà và ctv. (2017); Thanh và ctv. (2019); Nguyễn Văn Sáng và ctv. (2020) Sử dụng 7–11 microsatellite đạt tỷ lệ truy xuất đúng 81,3% đến 89,3%–90,7%; phát triển 10 locus mới dựa trên dữ liệu giải trình tự hệ gen cá tra (VN_pangasius, NCBI: GCA_003671635.1).
Chỉ tiêu miễn dịch trong chọn giống Camp và ctv. (2000); Galina (2017); Trần Thị Phương Dung và ctv. (2019) Số lượng bạch cầu, trung tâm đại thực bào sắc tố và hiệu giá kháng thể phản ánh mức độ đáp ứng miễn dịch và khả năng chịu đựng của cá khi nhiễm khuẩn E. ictaluri.

Khoảng trống học thuật mà luận án lựa chọn giải quyết là đánh giá di truyền đồng thời trên hai giai đoạn phát triển (cá hương và cá giống) của cùng một quần thể G1; kiểm định khả năng thay thế hoàn toàn thẻ từ PIT bằng bộ 9 chỉ thị microsatellite mới thiết kế; và xác định tổ hợp chỉ tiêu sinh lý - miễn dịch tại các thời điểm sau cảm nhiễm nhằm phục vụ đánh giá tính trạng kháng bệnh.

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

Luận án vận dụng lý thuyết di truyền số lượng để phân tách các thành phần phương sai di truyền cộng gộp ($\sigma_A^2$), phương sai môi trường chung ($\sigma_C^2$), phương sai số dư ($\sigma_E^2$) và phương sai kiểu hình ($\sigma_P^2$) thông qua mô hình tuyến tính hỗn hợp cá thể (Linear Mixed Model). Lý thuyết di truyền phân tử được áp dụng trong phân tích tính đa hình di truyền, cân bằng Hardy-Weinberg (HWE), tần số dị hợp tử quan sát ($H_O$), dị hợp tử mong đợi ($H_E$), hệ số đa hình thông tin (PIC) và phương pháp xác định quan hệ bố mẹ dựa trên khả năng (likelihood-based parentage assignment). Lý thuyết miễn dịch học thủy sản được dùng để giải thích biến thiên các chỉ số huyết học, thực bào và mô học dưới áp lực xâm nhiễm của mầm bệnh.

Phương pháp thực nghiệm chi tiết

1. Phương pháp cho Nội dung 1 (Di truyền số lượng và chọn lọc G1)

  • Sản xuất gia đình: Nuôi vỗ đàn cá bố mẹ G0, tiến hành phối ghép thứ bậc và giai thừa một phần để tạo ra 155 gia đình full-sib và half-sib. Ương nuôi riêng rẽ từng gia đình đến kích cỡ gắn thẻ từ PIT, chọn lọc được 130 gia đình đạt yêu cầu.
  • Thử nghiệm cảm nhiễm cá giống: Đánh dấu PIT cho 7.664 con cá giống thuộc 130 gia đình, chia đều vào 2 bể composite thể tích 16.000 L (3.832 con/bể), mật độ nuôi 1 con/4,1 L nước, nhiệt độ duy trì 26–28°C. Sử dụng phương pháp cohabitant kết hợp: ghép 35% cá cohabitant (được tiêm vi khuẩn E. ictaluri Gly09M liều $2 \times 10^5$ CFU/cá) sống chung với cá thí nghiệm, đồng thời bổ sung vi khuẩn vào nước bể ở mật độ $10^5$ CFU/mL sau 2 ngày.
  • Thử nghiệm cảm nhiễm cá hương: Đánh giá trên 33 gia đình, quy mô 50 con/gia đình bằng phương pháp ngâm trong dung dịch vi khuẩn.
  • Nuôi đánh giá tăng trưởng: Sử dụng 5.838 cá thể thuộc 130 gia đình (trung bình 45 con/gia đình) thả nuôi chung trong ao đất diện tích 2.000 m2. Sau 156 ngày, thu hoạch và ghi nhận các chỉ tiêu khối lượng (HW), chiều dài tổng (HL) và tỷ lệ sống (SUR).
  • Phân tích thống kê di truyền: Sử dụng phần mềm ASReml để ước tính hệ số di truyền ($h^2$), tương quan di truyền ($r_G$), giá trị chọn giống ước tính (EBV) và hiệu quả chọn lọc ($R$) tại các mốc cắt ngang tỷ lệ sống 50%, 25% và khi kết thúc cảm nhiễm.
Mô hình tuyến tính hỗn hợp phân tích kiểu hình cá thể:
y = Xb + Za + Wc + e
Trong đó:
y: Vector giá trị kiểu hình quan sát (sống/chết, thời gian sống, HW, HL)
b: Vector các ảnh hưởng cố định (bể cảm nhiễm, đợt sinh sản, giới tính)
a: Vector ảnh hưởng di truyền cộng gộp ngẫu nhiên ~ N(0, A*sigma_A^2)
c: Vector ảnh hưởng môi trường riêng rẽ ngẫu nhiên ~ N(0, I*sigma_C^2)
e: Vector sai số ngẫu nhiên ~ N(0, I*sigma_E^2)
X, Z, W: Các ma trận thiết kế tương ứng

2. Phương pháp cho Nội dung 2 (Giải pháp kỹ thuật hỗ trợ)

  • Ứng dụng chỉ thị microsatellite: Sử dụng bộ 10 locus microsatellite (Pahy-01, Pahy-02, Pahy-03, Pahy-04, Pahy-06, Pahy-10, Pahy-13, Pahy-15, Pahy-17, Pahy-18) tối ưu trong 1 phản ứng Multiplex PCR. Tách chiết ADN tổng số, khuếch đại gen và chạy điện di mao quản trên máy giải trình tự tự động. Sàng lọc trên 50 mẫu bố mẹ G0 và 50 mẫu con G1; thử nghiệm truy xuất phả hệ trên 50 gia đình (gồm 90 cá bố mẹ G0 và 500 cá con G1) bằng phần mềm chuyên dụng theo thuật toán likelihood-based; tính toán tần số alen rỗng (null-allele), tỷ lệ lỗi kiểu gen và tỷ lệ mismatch.
  • Đánh giá chỉ tiêu đáp ứng miễn dịch: Dựa vào giá trị EBV kháng bệnh ở điểm cắt 50% tỷ lệ sống, chọn ra 2 nhóm gia đình: nhóm kháng bệnh cao (KBC) và nhóm kháng bệnh thấp (KBT), mỗi nhóm gồm 3 gia đình. Cảm nhiễm lại trên tổng số 119 cá thể. Thu mẫu máu và mô tại 5 thời điểm: trước cảm nhiễm (0 hpi), 24, 48, 264 và 312 giờ sau cảm nhiễm (hpi).
  • Các chỉ tiêu miễn dịch phân tích:
    • Miễn dịch không đặc hiệu: Tổng hồng cầu (THC), tổng bạch cầu (TBC), bạch cầu đơn nhân (MONO), bạch cầu trung tính (NEU), bạch cầu lympho (LYM), hoạt lực thực bào (PA), chỉ số thực bào (PI), số lượng trung tâm đại thực bào sắc tố (TTĐTB) trong mô gan, thận, lách.
    • Miễn dịch đặc hiệu: Hiệu giá kháng thể (HGKT).
    • Phân tích thống kê: Sử dụng mô hình đơn biến và mô hình trung bình Bayes (Bayesian Model Averaging - BMA) trên phần mềm R; đánh giá qua diện tích dưới đường cong ROC (AUC), độ nhạy (Sensitivity - Sen) và độ đặc hiệu (Specificity - Spe).

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 bối cảnh sản xuất và xuất khẩu cá tra tại Đồng bằng sông Cửu Long, đặc điểm dịch tễ học và bệnh lý của vi khuẩn Edwardsiella ictaluri. Tác giả phân tích cơ chế gây bệnh tích gan thận mủ, đánh giá hiệu quả và hạn chế của các biện pháp phòng trị truyền thống (hóa chất, kháng sinh, chất kích thích miễn dịch, vắc-xin ALPHAJECT Panga 1). Tiếp đó, chương điểm luận tiến trình chọn giống thủy sản trên thế giới từ năm 1983 đến nay trên các đối tượng như cá hồi Đại Tây Dương (Salmo salar), cá rô phi vằn (Oreochromis niloticus), cá chép (Cyprinus carpio), tôm thẻ chân trắng (Litopenaeus vannamei). Tác giả so sánh các thiết kế phối giống (cặp đơn, thứ bậc, giai thừa), các phương pháp đánh dấu cá thể (cắt vây, nhuộm huỳnh quang VIE, thẻ từ PIT) và các phương pháp cảm nhiễm thực nghiệm (tiêm, ngâm, cohabitant). Chương 1 cũng phân tích lịch sử các chương trình chọn giống thủy sản tại Việt Nam, sự phát triển của các công cụ di truyền phân tử (Allozyme, AFLP, Microsatellite, SNP) và tiềm năng ứng dụng các chỉ số miễn dịch học trong chọn giống kháng bệnh.

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

Chương 2 mô tả chi tiết vật liệu sinh học, địa điểm, trang thiết bị và các bước thực nghiệm cho hai nội dung nghiên cứu lớn:

  1. Quy trình nuôi vỗ 710 cá bố mẹ G0, kỹ thuật kích thích sinh sản bằng HCG, kỹ thuật thụ tinh nhân tạo và ấp nở để thu 155 gia đình G1.
  2. Quy trình ương nuôi riêng rẽ trong hệ thống bể và giai, kỹ thuật cấy thẻ từ PIT vào xoang bụng cá khi đạt khối lượng 15–20 g.
  3. Thiết lập hệ thống bể cảm nhiễm 16.000 L, quy trình nuôi cấy vi khuẩn E. ictaluri trong môi trường BHIB, chuẩn hóa phương pháp cohabitant kết hợp (35% cá cohabitant, liều bổ sung $10^5$ CFU/mL) và thí nghiệm ngâm trên cá hương.
  4. Bố trí ao nuôi tăng trưởng 2.000 m2 trong 156 ngày và phương pháp thu thập số liệu kiểu hình thu hoạch.
  5. Thiết lập phản ứng Multiplex PCR cho 10 locus microsatellite, quy trình điện di mao quản và các thuật toán tính toán phả hệ.
  6. Kỹ thuật thu mẫu máu từ tĩnh mạch cuống đuôi, đếm tế bào máu, nhuộm mô học xác định TTĐTB tại gan, thận, lách, kỹ thuật ngưng kết kháng thể và các mô hình toán học ASReml, BMA trong R.

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 chuyên sâu đối chiếu với y văn quốc tế:

1. Thông số di truyền và định hướng chọn lọc G1

  • Quá trình sản xuất tạo ra 155 gia đình G1, ương sống và đánh dấu PIT thành công 130 gia đình với tổng số 7.664 cá thể giống cho thí nghiệm cảm nhiễm và 5.838 cá thể cho thí nghiệm nuôi ao đất.
  • Thí nghiệm cảm nhiễm cá giống cho thấy tỷ lệ chết tích lũy đạt đỉnh trong khoảng 7–15 ngày sau cảm nhiễm.
  • Kết quả ước tính các thành phần phương sai và hệ số di truyền qua phần mềm ASReml được tổng hợp cụ thể như sau:
Giai đoạn / Tính trạng Điểm cắt ngang dữ liệu Phương sai di truyền cộng gộp ($\sigma_A^2$) Phương sai môi trường chung ($\sigma_C^2$) Hệ số di truyền ($h^2 \pm SE$) Mức độ di truyền
Cá giống (Kháng bệnh) Cắt ngang sống 50% 0,082 0,018 $0,38 \pm 0,06$ Cao
Cá giống (Kháng bệnh) Cắt ngang sống 25% 0,054 0,012 $0,22 \pm 0,05$ Trung bình
Cá giống (Kháng bệnh) Kết thúc cảm nhiễm 0,031 0,009 $0,14 \pm 0,04$ Trung bình thấp
Cá hương (Kháng bệnh) Thí nghiệm ngâm 0,046 - $0,18 \pm 0,05$ Trung bình
Thu hoạch (Tăng trưởng) Khối lượng thân (HW) 1245,6 112,4 $0,32 \pm 0,07$ Trung bình - Cao
Thu hoạch (Tăng trưởng) Chiều dài thân (HL) 4,82 0,35 $0,28 \pm 0,06$ Trung bình
Thu hoạch (Tỷ lệ sống) Sống sót nuôi ao (SUR) 0,021 0,004 $0,12 \pm 0,03$ Trung bình thấp
  • Tương quan di truyền ($r_G$): Tương quan di truyền giữa tính trạng kháng bệnh ở cá giống và tính trạng kháng bệnh ở cá hương là tương quan thuận ($r_G = 0,42 - 0,58$). Tương quan di truyền giữa tính trạng kháng bệnh giai đoạn cá giống và khối lượng thu hoạch (HW) là tương quan thuận dương ($r_G = 0,21 - 0,34$). Tương quan giữa kháng bệnh và tỷ lệ sống trong ao nuôi thương phẩm đạt $r_G = 0,45 - 0,61$.
  • Hiệu quả chọn lọc dự tính: Chọn lọc dựa trên giá trị EBV tại mốc cắt ngang 50% mang lại hiệu quả chọn lọc đạt 9,6–11,2% ở thế hệ tiếp theo.

2. Kết quả ứng dụng chỉ thị Microsatellite truy xuất phả hệ

  • Sàng lọc ban đầu trên 50 cá bố mẹ G0 và 50 cá con G1 cho thấy 10 locus đều có tính đa hình cao, số alen dao động từ 4 đến 14 alen/locus, PIC > 0,5.
  • Khi phân tích phả hệ trên 50 gia đình (90 bố mẹ, 500 cá con), locus Pahy-02 xuất hiện tần số alen rỗng cao (> 0,05), dẫn đến sai số nhận diện. Sau khi loại bỏ Pahy-02, bộ 9 locus còn lại cho kết quả:
Chỉ số truy xuất phả hệ Bộ 10 locus Microsatellite Bộ 9 locus (loại bỏ Pahy-02)
Tỷ lệ gán đúng cả bố và mẹ (Tổng thể) 88,2% 93,4%
Tỷ lệ gán đúng trên gia đình không có half-sib 89,0% 93,0%
Tỷ lệ gán đúng trên gia đình có half-sib 87,4% 94,0%
Tỷ lệ lỗi ghi nhận kiểu gen (Genotyping error) 2,8% 0,9%
Số lượng locus có tần số null-allele cao 1 (Pahy-02) 0

3. Kết quả đánh giá chỉ tiêu đáp ứng miễn dịch

  • Phân tích động thái miễn dịch giữa nhóm gia đình kháng bệnh cao (KBC) và kháng bệnh thấp (KBT) trên 119 cá thể qua 5 thời điểm:
Diễn biến đáp ứng miễn dịch theo thời gian sau cảm nhiễm (hpi):
- Giai đoạn 0 hpi: Chưa có sự khác biệt rõ rệt giữa hai nhóm KBC và KBT.
- Giai đoạn 24 - 48 hpi (Phản ứng sớm): 
  + Bạch cầu trung tính (NEU) nhóm KBC tăng nhanh gấp 1,8 - 2,4 lần so với nhóm KBT.
  + Số lượng TTĐTB ở mô gan nhóm KBC đạt mật độ cao vượt trội (p < 0,05).
  + Hiệu giá kháng thể (HGKT) nhóm KBC tăng sớm hơn rõ rệt.
- Giai đoạn 264 - 312 hpi (Phục hồi):
  + Nhóm KBC duy trì hoạt lực thực bào (PA) và phục hồi chỉ số hồng cầu (THC) nhanh hơn.
  • Mô hình Bayesian Model Averaging (BMA) đơn biến và đa biến trên phần mềm R xác định giai đoạn 24–48 hpi là thời điểm tối ưu nhất để phân biệt cá thể thuộc nhóm KBC hay KBT. Tổ hợp ba chỉ tiêu gồm Bạch cầu trung tính (NEU), Hiệu giá kháng thể (HGKT) và Trung tâm đại thực bào sắc tố (TTĐTB) ở gan cho diện tích dưới đường cong ROC đạt $\text{AUC} = 0,86 - 0,91$, độ nhạy $\text{Sen} = 84,2%$ và độ đặc hiệu $\text{Spe} = 88,5%$.

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

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

  1. Lần đầu tiên ước tính các thông số di truyền ở giai đoạn cá hương: Luận án cung cấp dẫn liệu khoa học đầu tiên về hệ số di truyền tính trạng kháng bệnh gan thận mủ ở giai đoạn cá hương ($h^2 = 0,18$) và chứng minh tương quan di truyền thuận giữa tính trạng kháng bệnh ở cá hương và cá giống ($r_G > 0$). Điều này chứng minh chọn lọc kháng bệnh ở giai đoạn cá giống không gây ảnh hưởng bất lợi đến khả năng kháng bệnh của cá ở giai đoạn nhỏ hơn.
  2. Chứng minh tương quan thuận giữa kháng bệnh và tăng trưởng: Xác lập luận cứ di truyền khẳng định tính trạng kháng bệnh gan thận mủ và tốc độ tăng trưởng khối lượng khi thu hoạch có tương quan di truyền dương ($r_G = 0,21 - 0,34$), loại bỏ lo ngại về sự đánh đổi di truyền tiêu cực giữa hai tính trạng kinh tế quan trọng này.
  3. Xác lập cơ sở khoa học của chỉ thị miễn dịch trong chọn giống: Lần đầu tiên chứng minh định lượng khả năng sử dụng các chỉ tiêu đáp ứng miễn dịch sớm (bạch cầu trung tính, hiệu giá kháng thể, trung tâm đại thực bào gan tại 24–48 hpi) làm tiêu chí hỗ trợ đánh giá giá trị chọn giống kháng bệnh trên cá tra.

Đóng góp mới về mặt thực tiễn

  1. Bộ chỉ thị phân tử thay thế thẻ từ PIT: Thiết lập thành công bộ 9 chỉ thị microsatellite (Pahy-01, Pahy-03, Pahy-04, Pahy-06, Pahy-10, Pahy-13, Pahy-15, Pahy-17, Pahy-18) với tỷ lệ truy xuất bố mẹ đạt độ chính xác 93,4%. Giải pháp này cho phép thả nuôi chung các gia đình ngay từ giai đoạn cá bột, loại bỏ chi phí thẻ PIT, giảm áp lực cơ sở hạ tầng bể ương và triệt tiêu ảnh hưởng của môi trường ương riêng rẽ ($c^2$).
  2. Chuẩn hóa phương pháp đánh giá kiểu hình: Đề xuất sử dụng dữ liệu sống chết tại điểm cắt ngang tỷ lệ sống 50% và 25% trong mô hình cohabitant kết hợp để ước tính EBV, giúp rút ngắn thời gian thí nghiệm cảm nhiễm nhưng vẫn đảm bảo độ tin cậy chọn lọc cao ($h^2 = 0,22 - 0,38$).

Kiến nghị và giải pháp đề xuất

  • Áp dụng phương pháp chọn lọc gia đình kết hợp nhiều tính trạng (kháng bệnh và tăng trưởng) dựa trên giá trị EBV tại điểm cắt 50% để chọn cá bố mẹ cho thế hệ G2.
  • Đưa bộ 9 chỉ thị microsatellite vào quy trình sản xuất giống chọn lọc đại trà tại các trung tâm giống thủy sản vùng Đồng bằng sông Cửu Long.
  • Tiếp tục nghiên cứu tích hợp các chỉ tiêu miễn dịch nổi bật (NEU, HGKT, TTĐTB gan) vào chỉ số chọn giống nhằm nâng cao độ chính xác chọn lọc cá thể mà không cần cảm nhiễm gây chết toàn bộ đàn 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 ở việc đánh giá thế hệ chọn giống thứ nhất (G1), cần tiếp tục theo dõi biến động di truyền qua các thế hệ tiếp theo (G2, G3) để kiểm định sự ổn định của hiệu quả chọn lọc.
    • Thí nghiệm phân tích chỉ tiêu miễn dịch được thực hiện trên quy mô mẫu hạn chế (3 gia đình KBC và 3 gia đình KBT, tổng số 119 cá thể), chưa bao phủ toàn bộ 130 gia đình trong quần thể.
    • Bộ chỉ thị phân tử sử dụng là microsatellite; dù đạt tỷ lệ truy xuất 93,4% nhưng vẫn tồn tại 6,6% chưa định danh chính xác hoặc có sai số.
  • Hướng nghiên cứu tiếp:
    • Tiếp tục nhân dòng chọn lọc thế hệ G2 và G3 theo định hướng đa tính trạng (kháng bệnh kết hợp tăng trưởng và tỷ lệ philê).
    • Nghiên cứu phát triển chip SNP mật độ cao để nâng cao độ chính xác truy xuất phả hệ lên trên 98% và hướng tới chọn giống bằng giá trị chọn giống hệ gen (Genomic Selection - GS).
    • Ứng dụng các chỉ tiêu miễn dịch không xâm lấn hoặc phân tích biểu hiện gen miễn dịch để đánh giá trực tiếp trên cá thể sống dùng làm hậu bị sinh sản.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị cao đối với:

  • Nghiên cứu sinh, học viên cao học và sinh viên chuyên ngành Nuôi trồng Thủy sản, Di truyền học: Tham khảo phương pháp luận thiết kế thí nghiệm chọn giống gia đình, kỹ thuật cảm nhiễm cohabitant chuẩn hóa, quy trình phân tích mô hình hỗn hợp tuyến tính bằng phần mềm ASReml và phân tích mô hình BMA trong R.
  • Giảng viên và nhà nghiên cứu tại các Viện, Trường đại học: Sử dụng làm tài liệu giảng dạy chuyên đề Di truyền số lượng thủy sản, Ứng dụng công nghệ sinh học phân tử trong chọn giống động vật thủy sản.
  • Các Trung tâm giống thủy sản quốc gia và doanh nghiệp sản xuất giống cá tra: Áp dụng trực tiếp quy trình multiplex PCR 9 microsatellite để quản lý đàn cá bố mẹ, kiểm soát cận huyết và truy xuất phả hệ đàn giống thương phẩm, giúp tiết kiệm chi phí gắn thẻ PIT và nâng cao chất lượng kháng bệnh của cá tra giống cung ứng cho thị trường Đồng bằng sông Cửu Long.

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

1. Tại sao nghiên cứu đề xuất dùng mốc cắt ngang tỷ lệ sống 50% và 25% thay vì chờ kết thúc thí nghiệm cảm nhiễm?

Theo kết quả ước tính từ phần mềm ASReml trong luận án, hệ số di truyền ($h^2$) tính trạng kháng bệnh tại mốc cắt 50% ($h^2 = 0,38 \pm 0,06$) và 25% ($h^2 = 0,22 \pm 0,05$) đạt mức trung bình đến cao, cao hơn so với thời điểm kết thúc thí nghiệm ($h^2 = 0,14 \pm 0,04$). Việc xử lý số liệu tại các điểm cắt này phản ánh chính xác nhất biến thiên di truyền của quần thể khi áp lực bệnh ở đỉnh điểm, đồng thời giúp rút ngắn thời gian theo dõi cảm nhiễm và giảm thiểu các yếu tố gây nhiễu ngẫu nhiên ở giai đoạn muộn.

2. Nguyên nhân nào khiến locus Pahy-02 bị loại khỏi bộ chỉ thị truy xuất phả hệ?

Trong quá trình sàng lọc và kiểm tra trên 50 gia đình (90 cá bố mẹ và 500 cá con), locus Pahy-02 có tần số alen rỗng (null-allele frequency) ở mức cao vượt ngưỡng cho phép (> 0,05). Sự hiện diện của alen rỗng dẫn đến hiện tượng sai lệch kiểu gen dị hợp thành đồng hợp giả, gây ra tỷ lệ lỗi mismatch cao. Khi loại bỏ Pahy-02, tỷ lệ gán đúng cả bố và mẹ của bộ 9 locus còn lại tăng từ 88,2% lên 93,4% và tỷ lệ lỗi ghi nhận kiểu gen giảm xuống chỉ còn 0,9%.

3. Việc chọn giống tăng tính kháng bệnh gan thận mủ có làm giảm tốc độ tăng trưởng của cá tra không?

Không. Kết quả phân tích tương quan di truyền trong luận án chỉ ra rằng mối tương quan giữa tính trạng kháng bệnh gan thận mủ giai đoạn cá giống và khối lượng thu hoạch (HW) là tương quan thuận ($r_G = 0,21 - 0,34$). Do đó, việc chọn lọc nâng cao khả năng kháng bệnh không gây ảnh hưởng tiêu cực mà còn đồng thời cải thiện tốc độ sinh trưởng của cá tra khi nuôi thương phẩm trong ao đất.

4. Những chỉ tiêu miễn dịch nào có tiềm năng nhất để chẩn đoán khả năng kháng bệnh của cá tra?

Thông qua mô hình đa biến Bayesian Model Averaging (BMA) tại giai đoạn sớm 24–48 giờ sau cảm nhiễm (24–48 hpi), ba chỉ tiêu cho giá trị dự báo chính xác nhất gồm: tỷ lệ bạch cầu trung tính (NEU), hiệu giá kháng thể (HGKT) và số lượng trung tâm đại thực bào sắc tố (TTĐTB) trong mô gan. Tổ hợp này đạt độ nhạy 84,2%, độ đặc hiệu 88,5% và diện tích dưới đường cong ROC đạt 0,86–0,91.

5. Phương pháp gây bệnh thực nghiệm nào được tác giả sử dụng trên cá giống G1 và các thông số kỹ thuật là gì?

Luận án sử dụng phương pháp cohabitant kết hợp (cá bệnh sống chung với cá khỏe và bổ sung vi khuẩn). Cá cohabitant được tiêm vi khuẩn E. ictaluri Gly09M liều $2 \times 10^5$ CFU/cá, ghép vào bể với tỷ lệ 35% so với cá thí nghiệm. Bể composite có thể tích 16.000 L, mật độ 1 con/4,1 L nước, nhiệt độ duy trì 26–28°C. Sau 2 ngày ghép cá cohabitant, bổ sung vi khuẩn vào nước bể với nồng độ $10^5$ CFU/mL.

Kết luận

Luận án tiến sĩ của tác giả đã ứng dụng thành công phương pháp tiếp cận tích hợp giữa di truyền số lượng, di truyền phân tử và sinh lý miễn dịch phục vụ công tác chọn giống cá tra (Pangasianodon hypophthalmus) kháng bệnh gan thận mủ thế hệ G1. Công trình đã xác định hệ số di truyền tính trạng kháng bệnh đạt mức trung bình đến cao ($h^2 = 0,22 - 0,38$), xác lập mối tương quan di truyền thuận giữa khả năng kháng bệnh ở cá hương, cá giống với tốc độ tăng trưởng thu hoạch. Luận án đã tối ưu hóa bộ 9 chỉ thị microsatellite đạt độ chính xác truy xuất phả hệ 93,4% thay thế thẻ PIT và xác định tổ hợp chỉ tiêu miễn dịch sớm (NEU, HGKT, TTĐTB gan tại 24–48 hpi) có giá trị dự báo cao. Các kết quả này cung cấp cơ sở khoa học và giải pháp kỹ thuật trực tiếp để nâng cao hiệu quả chọn lọc cá tra kháng bệnh cho các thế hệ tiếp theo.