Tổng quan luận án

Luận án tiến sĩ nông nghiệp chuyên ngành Nuôi trồng Thủy sản (mã số 9 62 03 01) với đề tài "Ứng dụng di truyền phân tử và di truyền số lượng phục vụ chọn giống cá tra kháng bệnh gan thận mủ" được thực hiện tại Trường Đại học Nông Lâm Thành phố Hồ Chí Minh, hoàn thành vào tháng 12 năm 2022.

Về tính cấp thiết, cá tra (Pangasianodon hypophthalmus) là đối tượng thủy sản nước ngọt chủ lực của Việt Nam với diện tích ương nuôi giống đạt gần 4.000 ha, sản lượng nuôi đạt 1,56 triệu tấn và giá trị kim ngạch xuất khẩu đạt 1,61 tỉ USD vào năm 2021 theo số liệu của VASEP. Cá tra cũng là sản phẩm quốc gia được phê duyệt theo Quyết định số 439/QĐ-TTg của Thủ tướng Chính phủ ngày 16 tháng 04 năm 2012. Tuy nhiên, sự phát triển thâm canh đối mặt với áp lực dịch bệnh thường xuyên. Trong đó, bệnh gan thận mủ do vi khuẩn Edwardsiella ictaluri (E. ictaluri) gây tỷ lệ chết lên đến 90% nếu không xử lý kịp thời, xuất hiện chủ yếu ở cá dưới 4 tháng tuổi, đặc biệt là giai đoạn cá hương 21–30 ngày tuổi (tỷ lệ nhiễm 46%) và cá giống 40–90 ngày tuổi (tỷ lệ nhiễm 30%). Các biện pháp sử dụng kháng sinh tiềm ẩn rủi ro kháng thuốc và tồn dư hóa chất không đáp ứng tiêu chuẩn xuất khẩu; việc dùng chất kích thích miễn dịch hay vắc-xin đòi hỏi chi phí cao, tác dụng ngắn hạn và khó áp dụng đại trà. Do đó, chọn giống di truyền kháng bệnh là hướng tiếp cận bền vững. Quần thể cá tra bố mẹ chọn giống kháng bệnh gan thận mủ thế hệ G0 đã được Viện Nghiên cứu Nuôi trồng Thủy sản II thiết lập từ đề tài KC.06 với hệ số di truyền ước tính 0,27–0,37 và hiệu quả chọn lọc 8,3%.

Về khoảng trống nghiên cứu, công tác chọn giống cần được tiếp tục duy trì qua nhiều thế hệ để kiểm chứng hệ số di truyền và gia tăng hiệu quả chọn lọc. Bệnh gan thận mủ phát sinh mạnh từ giai đoạn cá hương, nhưng các thông số di truyền và tương quan di truyền giữa giai đoạn cá hương và cá giống chưa từng được công bố. Bên cạnh đó, phương pháp đánh dấu vật lý bằng thẻ từ PIT (Passive Integrated Transponder) áp dụng ở thế hệ trước buộc phải ương riêng rẽ các gia đình đến khối lượng 15–20 g trong 3–4 tháng, gây ra ảnh hưởng môi trường chung ($c^2$) làm giảm độ chính xác khi ước tính thông số di truyền. Các bộ chỉ thị microsatellite trước đây chỉ đạt tỷ lệ truy xuất 81,3%, chưa đủ điều kiện thay thế thẻ PIT. Đồng thời, việc đánh giá tính trạng kháng bệnh thuần túy qua kiểu hình sống/chết (nhị phân) khi kết thúc thí nghiệm cảm nhiễm chưa phản ánh toàn diện tương tác miễn dịch của cá thể, đòi hỏi phải tìm kiếm các chỉ tiêu đáp ứng miễn dịch sinh học bổ trợ.

Mục tiêu nghiên cứu của luận án:

  1. 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, 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.
  2. Mục tiêu cụ thể:
    • 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.
    • 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: Cá tra (Pangasianodon hypophthalmus) thế hệ G1 được sản xuất từ đàn bố mẹ G0 tại Viện Nghiên cứu Nuôi trồng Thủy sản II, cùng chủng vi khuẩn gây bệnh Edwardsiella ictaluri Gly09M.
  • Phạm vi nội dung và không gian: Thí nghiệm cảm nhiễm cá hương thực hiện trên 33 gia đình; thí nghiệm cảm nhiễm cá giống thực hiện trên 130 gia đình; nuôi đánh giá tăng trưởng thương phẩm trên 130 gia đình tại ao đất 2.000 $\text{m}^2$; thử nghiệm bộ 10 chỉ thị microsatellite trên 50 gia đình (90 cá bố mẹ G0 và 500 cá con G1); đánh giá 11 chỉ tiêu miễn dịch trên 119 cá thể thuộc 2 nhóm gia đình kháng bệnh cao và thấp (3 gia đình/nhóm). Nghiên cứu tiến hành tại Trại Nghiên cứu Thực nghiệm Thủy sản Gò Vấp, Trung tâm Quốc gia Giống Thủy sản Nước ngọt Nam Bộ, Trung tâm Quan trắc Môi trường và Bệnh 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 và Phòng thí nghiệm Đa dạng Di truyền và Tiến hóa thuộc Đại học KU Leuven (Vương quốc Bỉ).
  • Thời gian nghiên cứu: Nghiên cứu số liệu thực nghiệm thuộc thế hệ G1, hoàn thành phân tích và nghiệm thu 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 thành tựu chọn giống thủy sản trên thế giới và tại Việt Nam:

Trên thế giới, Neira (2010) thống kê có hơn 101 chương trình chọn giống thủy sản từ năm 1983, tập trung nhiều nhất ở cá rô phi, cá hồi Đại Tây Dương, cá hồi vân và cá rohu. Các tính trạng chọn lọc phổ biến bao gồm tăng trưởng, kháng bệnh, hệ số chuyển đổi thức ăn, tuổi thành thục và tỷ lệ philê (Gjedrem & Baranski, 2009). Chương trình chọn giống tăng trưởng thành công đầu tiên được thực hiện trên cá hồi Đại Tây Dương từ năm 1971, đạt mức tăng khối lượng 113% sau 5 thế hệ (Gjedrem, 2005). Chương trình chọn giống kháng bệnh đầu tiên triển khai từ năm 1990 trên cá hồi Đại Tây Dương kháng bệnh hoại tử huyết (Galina, 2017), tiếp nối bởi các nghiên cứu trên cá hồi coho (Yáñez et al., 2016), cá rô phi vằn (Shoemaker et al., 2017; Wonmongkol et al., 2018), cá tráp (Xiong et al., 2017) và tôm thẻ chân trắng (Argue et al., 2002; Trinh et al., 2019). Về phương pháp gây bệnh, ba kỹ thuật chính được áp dụng gồm tiêm vi khuẩn trực tiếp, ngâm trong dung dịch mầm bệnh và phương pháp sống chung (cohabitant), trong đó phương pháp cohabitant được xem là giải pháp tái hiện tự nhiên cơ chế lây truyền mầm bệnh (Nordmo et al., 1998; Gomez et al., 2010). Về truy xuất phả hệ, chỉ thị microsatellite đã được áp dụng trên cá trắm cỏ đạt độ chính xác 99,6% (Fu et al., 2013), cá tuyết đạt 97% (Ashton et al., 2016), cá rô phi vằn đạt 100% (Liu et al., 2017). Về chỉ tiêu đáp ứng miễn dịch, Camp et al. (2000) đã phân tích các chỉ tiêu bạch cầu, kháng thể và trung tâm đại thực bào sắc tố (TTĐTB) trên cá nheo Mỹ (Ictalurus punctatus) cảm nhiễm E. ictaluri.

Tại Việt Nam, các chương trình chọn giống khởi đầu từ năm 1984 trên cá chép (Trần Mai Thiên, 1998), chọn giống cá rô phi vằn dòng GIFT nhập từ Philippines năm 1999 (Nguyễn Công Dân et al., 2003; Luan, 2008; Nguyễn Hữu Ninh et al., 2010). Trên cá tra, Viện Nghiên cứu Nuôi trồng Thủy sản II đã thực hiện chọn giống nâng cao tốc độ tăng trưởng qua 3 thế hệ từ năm 2001 (Nguyễn Văn Sáng et al., 2009, 2012, 2014). Về tính trạng kháng bệnh, Viện Nghiên cứu Nuôi trồng Thủy sản III chọn lọc tôm thẻ chân trắng kháng bệnh đốm trắng (Trinh et al., 2019). Trên cá tra, Phạm Đình Khôi et al. (2010), Pham et al. (2020a, 2020b) đã thăm dò các mô hình ước tính thông số di truyền; đề tài KC.06 (2012–2015) do Trịnh Quốc Trọng et al. chủ trì đã xây dựng thành công đàn bố mẹ G0 gồm 177 gia đình, xác định hệ số di truyền tính trạng kháng bệnh gan thận mủ ở cá giống là 0,37. Về công nghệ phân tử, Nguyễn Hữu Ninh và Lưu Thị Hà Giang (2012) dùng 8 microsatellite truy xuất được 94% cá tra con; Bùi Thị Liên Hà et al. (2017) dùng 7 microsatellite đạt tỷ lệ 81,3%; Thanh et al. (2019) dùng 4 microsatellite tối ưu 2 multiplex PCR đạt tỷ lệ 90,7% (ở ngưỡng cho phép 1 mismatch) nhưng quy mô chỉ dừng lại ở 5 gia đình full-sib.

Khoảng trống luận án giải quyết:

  • Luận án lấp đầy khoảng trống dữ liệu về thông số di truyền tính trạng kháng bệnh gan thận mủ ở giai đoạn cá hương và mối tương quan di truyền giữa cá hương với cá giống và tăng trưởng thương phẩm.
  • Khắc phục hạn chế về tỷ lệ truy xuất phả hệ thấp của các nghiên cứu trước đây bằng cách phát triển bộ chỉ thị microsatellite mới từ dữ liệu giải trình tự toàn bộ hệ gen cá tra (Kim et al., 2018).
  • Giải quyết bài toán thiếu hụt chỉ thị sinh học chức năng thông qua việc xác định tập hợp chỉ tiêu miễn dịch có khả năng phân biệt cá thể kháng bệnh cao và thấp.

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

Cơ sở lý thuyết

  • Lý thuyết di truyền số lượng: Phân tích các tính trạng số lượng dựa trên phương sai di truyền cộng gộp ($\sigma^2_A$), phương sai môi trường ương riêng rẽ ($\sigma^2_C$), phương sai sai số dư ($\sigma^2_E$) và phương sai kiểu hình ($\sigma^2_P$). Hệ số di truyền ($h^2$) được xác định bằng tỷ số giữa phương sai di truyền cộng gộp trên tổng phương sai kiểu hình ($h^2 = \sigma^2_A / \sigma^2_P$). Tương quan di truyền ($r_G$) biểu thị mức độ ảnh hưởng đồng thời của các nhóm gen đa hiệu lên hai tính trạng. Mô hình tuyến tính hỗn hợp cá thể (Animal Mixed Model) được sử dụng để ước tính giá trị giống (EBV) và hiệu quả chọn lọc ($R$).
  • Cơ sở di truyền phân tử: Ứng dụng chỉ thị chuỗi lặp lại đơn (Simple Sequence Repeat - SSR hay microsatellite) có tính đồng trội và tuân theo quy luật di truyền Mendel. Sử dụng phương pháp xác định quan hệ bố mẹ dựa trên khả năng (likelihood-based assignment) nhằm kiểm soát lỗi đột biến và lỗi đọc kiểu gen, kết hợp kiểm định cân bằng Hardy-Weinberg (HWE) và đánh giá tần số alen rỗng (null-allele).
  • Cơ sở miễn dịch học thủy sản: Phản ứng tự vệ của vật chủ trước mầm bệnh E. ictaluri thông qua hệ thống miễn dịch không đặc hiệu (tế bào máu, thực bào của đại thực bào, phản ứng melanin hóa tại trung tâm đại thực bào sắc tố) và miễn dịch đặc hiệu (hiệu giá kháng thể dịch thể).

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

Luận án bố trí thành 2 nội dung chính với các quy trình kỹ thuật cụ thể:

1. Phương pháp di truyền số lượng cho đàn cá G1:

  • Nuôi vỗ và sinh sản: Sử dụng đàn cá bố mẹ G0 nuôi vỗ tại Trung tâm Quốc gia Giống Thủy sản Nước ngọt Nam Bộ, kích dục tố sử dụng HCG. Thực hiện phối ghép cặp theo thiết kế bán giai thừa và thứ bậc, tạo ra 155 gia đình cá G1 (bao gồm cả full-sib và half-sib). Ương riêng rẽ từng gia đình trong các giai/bể đến kích cỡ đánh dấu, ghi nhận 130 gia đình sống sót.
  • Đánh dấu cá thể: Cấy thẻ từ PIT (Passive Integrated Transponder) vào xoang bụng của 7.664 cá thể cá giống thuộc 130 gia đình phục vụ cảm nhiễm và 5.838 cá thể nuôi ao đất.
  • Thử nghiệm cảm nhiễm cá giống: Phân bố 7.664 cá giống vào 2 bể composite (dung tích 16.000 lít/bể), mật độ 3.832 con/bể (tương đương 1 con/4,1 lít nước), nhiệt độ 26–28°C. Chủng vi khuẩn sử dụng là E. ictaluri Gly09M. Phương pháp cảm nhiễm: Cohabitant kết hợp (cá cảm nhiễm được tiêm vi khuẩn liều $2 \times 10^5\text{ CFU/cá}$ và thả ghép với tỷ lệ 35%; sau 2 ngày bổ sung vi khuẩn vào nước bể với nồng độ $10^5\text{ CFU/mL}$). Vớt cá chết mỗi 3 giờ/lần, kiểm tra mã PIT, ghi nhận giờ chết. Thí nghiệm dừng khi cá ngừng chết 5 ngày liên tục.
  • Thử nghiệm cảm nhiễm cá hương: Thực hiện trên 33 gia đình (50 cá thể/gia đình) bằng phương pháp ngâm vi khuẩn để xác định tỷ lệ sống và thời gian sống.
  • Đánh giá tăng trưởng: Thả nuôi chung 5.838 cá thể của 130 gia đình (bình quân 45 con/gia đình) trong ao đất 2.000 $\text{m}^2$ trong 156 ngày, thu hoạch đo đạc khối lượng (HW), chiều dài (HL) và tỷ lệ sống (SUR).
  • Xử lý số liệu di truyền: Dùng mô hình tuyến tính hỗn hợp cá thể trên phần mềm ASReml, phân tích tính trạng sống/chết (SUR) và thời gian sống (TIME) tại các điểm cắt ngang tỷ lệ sống khác nhau (trong đó có mốc 50% và 25%).

2. Phương pháp nghiên cứu giải pháp kỹ thuật hỗ trợ:

  • Sàng lọc bộ chỉ thị microsatellite: Sử dụng bộ dữ liệu hệ gen P. hypophthalmus (GenBank: GCA_003671635.1), khai thác 10 chỉ thị microsatellite được tối ưu trong một phản ứng multiplex PCR đơn (gồm Pahy-01, Pahy-02, Pahy-03, Pahy-04, Pahy-06, Pahy-10, Pahy-13, Pahy-15, Pahy-17, Pahy-18). Tách chiết ADN tổng số, khuếch đại PCR gắn huỳnh quang, điện di mao quản trên máy tự động. Kiểm tra độ ổn định, tính đa hình trên 50 cá thể bố mẹ G0 và 50 con G1, sau đó thử nghiệm truy xuất trên 50 gia đình (90 cá bố mẹ G0 và 500 cá con G1) bằng phần mềm chuyên dụng áp dụng thuật toán likelihood-based.
  • Xác định chỉ tiêu miễn dịch chức năng: Chọn 2 nhóm gia đình tương phản gồm 3 gia đình có giá trị giống ước tính kháng bệnh cao (KBC) và 3 gia đình kháng bệnh thấp (KBT) dựa trên EBV tại mốc cắt ngang tỷ lệ sống 50%. Tái cảm nhiễm 119 cá thể của 2 nhóm bằng vi khuẩn E. ictaluri. Thu mẫu máu và mô tạng 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). Định lượng 11 chỉ tiêu miễn dịch: 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), hiệu giá kháng thể (HGKT) và diện tích/số lượng trung tâm đại thực bào sắc tố (TTĐTB) ở gan, thận, lách. Xử lý thống kê qua mô hình đơn biến, đa biến, mô hình tiên lượng trung bình Bayes (Bayesian Model Average - BMA) trên phần mềm R và phân tích diện tích dưới đường cong ROC (AUC).

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

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

Chương 1 cung cấp cơ sở dữ liệu chuyên khảo về đặc điểm sinh học và kinh tế của cá tra tại Đồng bằng sông Cửu Long; lịch sử phát sinh, dịch tễ học và bệnh lý học của bệnh gan thận mủ do vi khuẩn Edwardsiella ictaluri. Điểm luận chi tiết các thành tựu chọn giống thủy sản quốc tế và trong nước, phân tích ưu nhược điểm của các phương pháp đánh dấu cá thể (cắt vây, nhuộm màu VIE, gắn chip PIT) và các phương pháp cảm nhiễm thực nghiệm. Chương này thiết lập cơ sở lý thuyết chứng minh rằng việc áp dụng công nghệ phân tử (microsatellite multiplex PCR) và sinh học miễn dịch là hướng đi tất yếu để khắc phục những nhược điểm cố hữu của phương pháp chọn giống truyền thống dựa trên dấu từ PIT và kiểu hình nhị phân.

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

Chương 2 trình bày toàn bộ thiết kế thực nghiệm, quy trình phối ghép bố mẹ G0, quy trình ương nuôi phân tách gia đình, kỹ thuật gắn thẻ từ PIT, phương pháp cảm nhiễm chuẩn hóa bằng cohabitant kết hợp và cảm nhiễm ngâm. Chương này đặc tả các mô hình toán học thống kê: mô hình động vật tuyến tính hỗn hợp ước tính các phương sai thành phần ($\sigma^2_A, \sigma^2_C, \sigma^2_E, \sigma^2_P$), hệ số di truyền ($h^2$), tương quan di truyền ($r_A$ hoặc $r_G$); mô tả quy trình phản ứng multiplex PCR cho 10 locus microsatellite; quy trình nhuộm Giemsa đếm tế bào máu, kỹ thuật mô học xác định TTĐTB, phương pháp ngưng kết kháng nguyên - kháng thể và các giải thuật tiên lượng Bayes (BMA), đường cong ROC.

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

Chương 3 là phần trọng tâm của luận án, trình bày hai tuyến kết quả nghiên cứu độc lập nhưng bổ trợ chặt chẽ cho nhau:

1. Ứng dụng di truyền số lượng ước tính thông số di truyền thế hệ G1

  • Sản xuất đàn con G1: Đã ghép cặp sinh sản thành công tạo ra 155 gia đình cá G1. Giai đoạn ương từ bột lên hương và giống chịu hao hụt tự nhiên, còn lại 130 gia đình đạt kích thước cấy thẻ PIT với tổng số 7.664 cá thể cảm nhiễm và 5.838 cá thể nuôi ao.
  • Đánh giá cảm nhiễm cá hương và cá giống: Ở cá giống, triệu chứng gan thận mủ xuất hiện điển hình với các đốm trắng kích thước 1–3 mm trên gan, thận, tỳ tạng. Tỷ lệ chết tích lũy của 130 gia đình cá giống có sự biến thiên lớn giữa các gia đình, thể hiện tính đa dạng kiểu hình di truyền rõ rệt. Đường cong biểu diễn Kaplan-Meier phản ánh sự phân ly xác suất sống sót giữa các gia đình theo thời gian.
  • Nuôi tăng trưởng thương phẩm: Sau 156 ngày nuôi chung trong ao đất 2.000 $\text{m}^2$, tỷ lệ sống bình quân đạt mức cao, khối lượng và chiều dài thân phân bố chuẩn, chứng minh cá thể sau gắn PIT sinh trưởng ổn định.
Tính trạng và thời điểm cắt ngang Hệ số di truyền ($h^2$) Mức độ di truyền
Kháng bệnh cá hương (ngâm vi khuẩn) Thấp đến trung bình Có biến dị di truyền
Kháng bệnh cá giống tại mốc sống 50% 0,22 – 0,38 Trung bình đến cao
Kháng bệnh cá giống tại mốc sống 25% 0,22 – 0,38 Trung bình đến cao
Tăng trưởng khối lượng thu hoạch (HW) Thấp đến trung bình Ổn định kiểu hình
Tăng trưởng chiều dài thu hoạch (HL) Thấp đến trung bình Ổn định kiểu hình
Tỷ lệ sống thu hoạch ao đất (SUR) Thấp Chịu ảnh hưởng ngoại cảnh
  • Tương quan di truyền ($r_G$):
    • Tương quan giữa tính trạng kháng bệnh giai đoạn cá hương và cá giống mang giá trị dương (thuận).
    • Tương quan giữa tính trạng kháng bệnh giai đoạn cá giống và các chỉ tiêu tăng trưởng lúc thu hoạch (khối lượng HW, chiều dài HL) là tương quan thuận ($r_G > 0$). Phát hiện này có ý nghĩa then chốt: việc chọn giống nâng cao sức kháng bệnh gan thận mủ ở giai đoạn cá giống hoàn toàn không gây ảnh hưởng tiêu cực lên khả năng kháng bệnh ở giai đoạn cá hương, đồng thời duy trì hoặc cải thiện tốc độ sinh trưởng thương phẩm.
  • Hiệu quả chọn lọc dự kiến: Đạt mức cao khi áp dụng chọn lọc gia đình dựa trên giá trị EBV tại các điểm cắt ngang tỷ lệ sống 50% và 25%.

2. Kết quả nghiên cứu các giải pháp kỹ thuật hỗ trợ

  • Nghiên cứu bộ chỉ thị microsatellite: Sàng lọc 10 locus (Pahy-01 đến Pahy-18) cho thấy các locus đều biểu hiện tính đa hình và ổn định cao trên 50 bố mẹ G0 và 50 cá con G1. Tuy nhiên, khi kiểm tra trên quy mô 50 gia đình (90 cá bố mẹ G0, 500 cá con G1), locus Pahy-02 phát hiện có tần số alen rỗng (null-allele) cao, có nguy cơ gây sai lệch kiểu gen. Sau khi loại bỏ Pahy-02, bộ chỉ thị còn lại gồm 9 locus microsatellite đạt tỷ lệ truy xuất phả hệ chính xác chung là 93,4%. Cụ thể: đối với các gia đình không có half-sib, tỷ lệ truy xuất chính xác đạt 93,0%; đối với các gia đình có cấu trúc half-sib, tỷ lệ truy xuất đạt 94,0%. Kết quả chứng minh bộ 9 microsatellite này đủ độ tin cậy để thay thế hoàn toàn thẻ PIT.
  • Nghiên cứu chỉ tiêu đáp ứng miễn dịch tiềm năng:
    • So sánh giữa 2 nhóm KBC và KBT (3 gia đình/nhóm, 119 cá thể) qua các mốc thời gian 0, 24, 48, 264 và 312 hpi cho thấy sự khác biệt có ý nghĩa thống kê ở nhiều chỉ tiêu miễn dịch.
    • Giai đoạn 24–48 hpi và 264–312 hpi là các mốc phản ánh rõ nhất sự phân hóa phản ứng miễn dịch giữa hai nhóm. Trong đó, giai đoạn 24–48 hpi có giá trị chẩn đoán nổi bật nhất với 3 chỉ tiêu then chốt: bạch cầu trung tính (NEU), hiệu giá kháng thể (HGKT) và số lượng/diện tích trung tâm đại thực bào sắc tố (TTĐTB) tại gan.
    • Mô hình BMA tối ưu trên phần mềm R cùng diện tích dưới đường cong ROC (AUC) khẳng định độ nhạy (Sensitivity) và độ đặc hiệu (Specificity) cao của mô hình tích hợp các chỉ tiêu miễn dịch này trong việc nhận diện cá thể có tiềm năng di truyền kháng bệnh.

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

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

  1. Lần đầu tiên trong công tác chọn giống cá tra kháng bệnh, nghiên cứu đã ước tính thành công các thông số di truyền (hệ số di truyền, tương quan di truyền) đối với tính trạng kháng bệnh gan thận mủ ở giai đoạn cá hương cho quần thể cá tra chọn giống tại Viện Nghiên cứu Nuôi trồng Thủy sản II.
  2. Chứng minh thực nghiệm mối tương quan di truyền thuận giữa tính trạng kháng bệnh gan thận mủ giai đoạn cá giống với tính trạng kháng bệnh giai đoạn cá hương và các tính trạng sinh trưởng (khối lượng, chiều dài) khi thu hoạch. Điều này cung cấp cơ sở khoa học khẳng định công tác chọn lọc cá giống kháng bệnh không làm suy thoái tiềm năng tăng trưởng thương phẩm.

Những đóng góp mới về mặt thực tiễn và kỹ thuật

  1. Phát triển và ứng dụng thành công bộ chỉ thị phân tử gồm 9 microsatellite (tối ưu hóa từ 10 locus ban đầu, loại bỏ Pahy-02 do lỗi null-allele) thực hiện trong 01 phản ứng multiplex PCR đơn, đạt độ chính xác truy xuất phả hệ bố mẹ 93,4% (đặc biệt đạt 94,0% trên gia đình half-sib). Bộ công cụ này cho phép loại bỏ phương pháp đánh dấu cơ học bằng thẻ PIT, giúp ương nuôi chung cá từ giai đoạn sớm, triệt tiêu ảnh hưởng của môi trường ương riêng rẽ ($c^2$) và giảm chi phí sản xuất giống.
  2. Lần đầu tiên xác lập tập hợp các chỉ tiêu đáp ứng miễn dịch (đặc biệt là tỷ lệ bạch cầu trung tính, hiệu giá kháng thể và diện tích trung tâm đại thực bào sắc tố ở gan tại thời điểm 24–48 giờ sau cảm nhiễm) có khả năng phân biệt chuẩn xác giữa nhóm gia đình kháng bệnh cao và kháng bệnh thấp, mở ra hướng ứng dụng biomarker sinh học trong chọn giống gián tiếp.

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

  • Đề xuất áp dụng tính trạng sống/chết tại các thời điểm cắt ngang tỷ lệ sống toàn bộ 50% và 25% trong thí nghiệm cảm nhiễm để tính toán giá trị giống ước tính (EBV) cho công tác chọn lọc đàn bố mẹ thế hệ kế tiếp.
  • Đề xuất đưa bộ 9 chỉ thị microsatellite multiplex PCR vào quy trình thường quy để quản lý phả hệ đàn cá tra chọn giống ở thế hệ G2 và các thế hệ tiếp theo.
  • Ứng dụng quy trình cảm nhiễm cohabitant kết hợp (35% cá tiêm ghép chung và bổ sung vi khuẩn nồng độ $10^5\text{ CFU/mL}$) làm tiêu chuẩn đánh giá tính trạng kháng bệnh gan thận mủ trên cá tra giống.

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

Về mặt hạn chế:

  • Thí nghiệm đánh giá trên giai đoạn cá hương mới triển khai được trên 33 gia đình bằng phương pháp cảm nhiễm ngâm, quy mô nhỏ hơn so với 130 gia đình ở giai đoạn cá giống.
  • Nghiên cứu chỉ tiêu miễn dịch chức năng mới thực hiện ở mức quy mô thực nghiệm đối chứng trên 2 nhóm gia đình tương phản (3 gia đình KBC và 3 gia đình KBT với tổng số 119 cá thể), chưa phân tích trên toàn bộ 130 gia đình do rào cản chi phí và độ phức tạp của xét nghiệm mô bệnh học, huyết học.
  • Chỉ thị Pahy-02 bộc lộ tần số null-allele cao trong quần thể nghiên cứu nên phải loại trừ, chỉ còn lại 9 locus trong multiplex PCR.

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

  • Tiếp tục thực hiện chọn lọc thế hệ G2 và các thế hệ tiếp nối từ đàn cá G1 đã được tuyển chọn nhằm kiểm chứng đáp ứng chọn lọc thực tế ($R$) so với giá trị lý thuyết.
  • Ứng dụng bộ 9 microsatellite để ương chung các gia đình ngay từ giai đoạn cá bột mới nở nhằm triệt tiêu hoàn toàn phương sai môi trường chung ($c^2$).
  • Mở rộng số lượng gia đình và cỡ mẫu cá thể trong đánh giá các chỉ tiêu miễn dịch (NEU, HGKT, TTĐTB gan ở 24–48 hpi), hướng tới tích hợp các chỉ tiêu này như một tính trạng thứ cấp trong mô hình chỉ số chọn giống (Selection Index) hoặc BLUP đa tính trạng.

Giá trị tham khảo

Luận án mang lại giá trị tham khảo thực tiễn và học thuật cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và cán bộ nghiên cứu chuyên ngành Nuôi trồng Thủy sản, Di truyền - Chọn giống động vật: Tham khảo phương pháp luận thiết kế thí nghiệm phối ghép bán giai thừa/thứ bậc, quy trình xử lý dữ liệu di truyền số lượng bằng mô hình tuyến tính hỗn hợp trên phần mềm ASReml, và kỹ thuật phân tích phả hệ dựa trên xác suất (likelihood-based parentage assignment).
  • Các nhà chọn giống và cơ sở sản xuất giống thủy sản tại Đồng bằng sông Cửu Long: Có thể tiếp nhận quy trình kỹ thuật bộ 9 microsatellite multiplex PCR để quản lý đàn bố mẹ, kiểm soát cận huyết và truy xuất nguồn gốc con giống mà không phụ thuộc vào thẻ từ PIT đắt tiền.
  • Các nhà nghiên cứu bệnh học thủy sản và miễn dịch học: Tham khảo bộ quy chuẩn cảm nhiễm vi khuẩn E. ictaluri chủng Gly09M bằng phương pháp cohabitant kết hợp, cùng khung đánh giá động học 11 chỉ tiêu miễn dịch qua các mốc thời gian sau nhiễm bệnh.
  • Cơ quan quản lý nhà nước (Bộ Nông nghiệp và Phát triển Nông thôn, Tổng cục Thủy sản): Dữ liệu của luận án cung cấp căn cứ khoa học phục vụ triển khai Đề án phát triển công nghệ sinh học ngành thủy sản và chiến lược tái cơ cấu sản xuất giống cá tra quốc gia theo hướng kháng bệnh, an toàn sinh học.

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

1. Hệ số di truyền ($h^2$) của tính trạng kháng bệnh gan thận mủ ở cá tra giống thế hệ G1 đạt mức bao nhiêu và có ý nghĩa như thế nào đối với công tác chọn giống? Trong thí nghiệm cảm nhiễm bằng phương pháp cohabitant kết hợp trên 130 gia đình cá giống, hệ số di truyền ước tính tại các thời điểm cắt ngang tỷ lệ sống toàn bộ 50% và 25% dao động trong khoảng từ 0,22 đến 0,38 (thuộc mức trung bình đến cao). Mức độ di truyền này chứng minh tính trạng kháng bệnh gan thận mủ ở cá tra chịu sự kiểm soát đáng kể của phương sai di truyền cộng gộp, cho thấy việc chọn giống cải thiện sức đề kháng của cá tra thông qua chọn lọc gia đình hoàn toàn khả thi và có thể đạt hiệu quả chọn lọc cao ở các thế hệ tiếp theo.

2. Việc chọn lọc nâng cao khả năng kháng bệnh có làm ảnh hưởng xấu đến tốc độ sinh trưởng của cá tra hay không? Kết quả phân tích di truyền đa tính trạng của luận án cho thấy mối tương quan di truyền giữa tính trạng kháng bệnh gan thận mủ ở giai đoạn cá giống với tính trạng tăng trưởng (khối lượng và chiều dài thân lúc thu hoạch sau 156 ngày nuôi ao) là tương quan thuận ($r_G > 0$). Đồng thời, tương quan giữa sức kháng bệnh ở cá giống và cá hương cũng mang giá trị dương. Do đó, việc chọn giống nâng cao khả năng kháng bệnh ở cá giống không làm suy giảm sức đề kháng của cá hương và không làm ảnh hưởng tiêu cực đến khả năng tăng trọng của cá thương phẩm.

3. Tại sao bộ chỉ thị microsatellite ban đầu gồm 10 locus lại bị rút gọn còn 9 locus, và tỷ lệ truy xuất phả hệ đạt được là bao nhiêu? Trong thử nghiệm truy xuất phả hệ trên 50 gia đình (90 cá bố mẹ G0 và 500 cá con G1), locus Pahy-02 bộc lộ tần số alen rỗng (null-allele) ở mức cao, gây sai lệch kiểu gen khi đối chiếu giữa bố mẹ và con. Nhằm đảm bảo độ chính xác tối đa, tác giả đã loại bỏ Pahy-02 khỏi mô hình phân tích. Bộ chỉ thị 9 microsatellite còn lại đạt tỷ lệ truy xuất chính xác cả bố và mẹ chung cho các gia đình là 93,4%, trong đó nhóm gia đình không có half-sib đạt 93,0% và nhóm gia đình có half-sib đạt 94,0%.

4. Quy trình cảm nhiễm thực nghiệm bệnh gan thận mủ trên cá giống được chuẩn hóa theo phương thức nào? Nghiên cứu sử dụng phương pháp cohabitant kết hợp: cá cảm nhiễm (cá cohabitant) được tiêm vi khuẩn E. ictaluri Gly09M với liều $2 \times 10^5\text{ CFU/cá}$ vào xoang bụng, sau 2 ngày được thả vào bể sống chung với cá thí nghiệm theo tỷ lệ ghép 35%. Sau đó 2 ngày, tiếp tục bổ sung vi khuẩn trực tiếp vào nước bể nuôi để đạt nồng độ $10^5\text{ CFU/mL}$. Mật độ nuôi đạt 1 con/4,1 lít nước ở nhiệt độ ổn định 26–28°C. Quy trình này phản ánh tự nhiên quá trình xâm nhập của mầm bệnh qua các hàng rào miễn dịch của cá thể.

5. Những chỉ tiêu miễn dịch nào và ở thời điểm nào cho phép phân biệt chuẩn xác nhất giữa nhóm cá thể kháng bệnh cao và kháng bệnh thấp? Trong số 11 chỉ tiêu đáp ứng miễn dịch khảo sát trên 119 cá thể thuộc 2 nhóm gia đình tương phản (KBC và KBT), giai đoạn 24–48 giờ sau cảm nhiễm (hpi) thể hiện khả năng phân biệt rõ nét nhất. Tại thời điểm này, 3 chỉ tiêu sinh học có giá trị dự đoán tối ưu nhất qua phân tích mô hình BMA và đường cong ROC là: tỷ lệ tế bào bạch cầu trung tính (NEU), hiệu giá kháng thể (HGKT) và số lượng/diện tích các trung tâm đại thực bào sắc tố (TTĐTB) tại mô gan.


Kết luận

Luận án tiến sĩ của tác giả đã giải quyết trọn vẹn mục tiêu ước tính các thông số di truyền tính trạng kháng bệnh gan thận mủ trên quần thể cá tra thế hệ G1, xác định hệ số di truyền ở mức trung bình đến cao ($h^2 = 0,22 - 0,38$) và chứng minh mối tương quan di truyền thuận giữa tính trạng kháng bệnh với sinh trưởng thương phẩm. Nghiên cứu đã tối ưu hóa thành công bộ 9 chỉ thị microsatellite multiplex PCR với độ chính xác truy xuất phả hệ đạt 93,4%, đồng thời xác định được tập hợp các chỉ tiêu đáp ứng miễn dịch (bạch cầu trung tính, hiệu giá kháng thể, trung tâm đại thực bào sắc tố ở gan tại 24–48 hpi) có tiềm năng chẩn đoán phân loại khả năng kháng bệnh. Các kết quả này cung cấp luận cứ khoa học và giải pháp công nghệ trực tiếp phục vụ chương trình chọn giống cá tra kháng bệnh bền vững tại Việt Nam.