Tổng quan về luận án

Biển ven bờ Việt Nam sở hữu chiều dài đường bờ biển trên 3.260 km trải dài qua 13 vĩ độ cùng vùng đặc quyền kinh tế hơn 1 triệu km², tạo điều kiện thuận lợi phát triển kinh tế thủy sản. Tuy nhiên, theo số liệu từ Bộ Thủy sản (2005), hàng năm có khoảng 40% đến 50% các cơ sở và trại nuôi thủy sản bị thiệt hại nghiêm trọng do các bệnh ký sinh trùng gây ra. Trong các tác nhân gây bệnh thủy sinh vật, ngành giun tròn (Nematoda) giữ vai trò nguy hiểm hàng đầu khi không chỉ làm suy giảm sức sinh trưởng, gây tổn thương mô nội tạng và làm chết cá nuôi, mà còn tiềm ẩn nguy cơ lây truyền các hội chứng bệnh truyền lây từ động vật sang người (zoonoses), điển hình như bệnh giun Anisakis (Anisakiasis) hoặc bệnh giun Capillaria từng làm chết 77 người tại Philippines giai đoạn 1964-1967.

Trong quần xã cá biển, bộ cá Vược (Perciformes) là bậc phân loại lớn nhất trong phân lớp cá vây tia (Actinopterygii), chiếm tới 40% tổng số loài cá có xương đã biết trên thế giới với khoảng 7.000 loài. Tại Việt Nam, bộ cá Vược ghi nhận 1.078 loài thuộc 352 giống, 78 họ, trong đó có khoảng 90 loài cá kinh tế giữ tỷ trọng khai thác và nuôi biển chủ lực (Nguyễn Tấn Trịnh và cs., 1996). Mặc dù có vai trò kinh tế - sinh thái tối quan trọng, các công trình nghiên cứu hệ thống học và dịch tễ học ký sinh trùng đường ruột trên bộ cá Vược tại Việt Nam trong nhiều thập kỷ qua còn phân tán, thiếu tính đồng bộ và chưa được khảo sát đa vùng địa lý bằng các công cụ khoa học tiên tiến.

Luận án tiến sĩ sinh học của nghiên cứu sinh Hoàng Văn Hiền (2019) với đề tài "Thành phần loài giun tròn (Nematoda) ký sinh ở một số loài cá thuộc bộ cá Vược (Perciformes) ở biển ven bờ Việt Nam" (Mã số: 9 42 01 05, chuyên ngành Ký sinh trùng học) dưới sự hướng dẫn của PGS. Hà Duy Ngọ và TS. Nguyễn Văn Đức tại 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 đã giải quyết khoảng trống học thuật này.

Nghiên cứu tập trung giải quyết hai câu hỏi nghiên cứu cốt lõi:

  1. Câu hỏi 1 (Q1): Cấu trúc thành phần loài và vị trí phân loại học của khu hệ giun tròn ký sinh trên các đại diện bộ cá Vược tại các vùng sinh thái biển ven bờ Việt Nam được xác lập như thế nào dưới góc độ kết hợp hình thái học siêu cấu trúc và chỉ thị phân tử?
  2. Câu hỏi 2 (Q2): Quy luật phân bố dịch tễ học (tỷ lệ nhiễm, cường độ nhiễm, tính đặc hiệu vật chủ) của các loài giun tròn biến thiên ra sao theo bậc phân loại họ/loài cá và theo các phân vùng địa lý sinh học biển Bắc - Trung - Nam?

Hai giả thuyết nghiên cứu tương ứng được kiểm định:

  • Giả thuyết 1 (H1): Khu hệ giun tròn ký sinh ở bộ cá Vược biển ven bờ Việt Nam có độ đa dạng loài rất cao, chứa đựng các ghi nhận mới cho khu hệ quốc gia và các taxon tiềm năng mới cho khoa học do vị trí giao thoa sinh thái nhiệt đới Tây Thái Bình Dương.
  • Giả thuyết 2 (H2): Tỷ lệ và cường độ nhiễm giun tròn có sự phân hóa rõ rệt theo tập tính dinh dưỡng, tầng nước sinh sống của từng họ cá và chịu ảnh hưởng trực tiếp bởi điều kiện thủy văn, hải lưu cục bộ của từng vùng biển ven bờ.

Khung lý thuyết của luận án được định hình vững chắc dựa trên Hệ thống phân loại hiện đại ngành Giun tròn của De Ley & Blaxter (2004) kết hợp phân loại sinh học kinh điển của Maggenti (1991), Skrjabin (1928, 1949), Yamaguti (1961) và Lý thuyết vòng đời - tương tác ký sinh vật chủ của Cox (1993) và Malakhov (1986).

Đóng góp đột phá của luận án được lượng hóa qua quy mô khảo sát mẫu: tiến hành mổ khám toàn diện 3.735 cá thể cá thuộc 129 loài, 38 họ cá Vược tại 9 tỉnh ven biển đại diện cho 3 miền Bắc - Trung - Nam trong giai đoạn 2014-2017. Luận án đã thu thập và giám định 7.951 mẫu giun tròn (gồm 2.717 cá thể trưởng thành và 5.234 ấu trùng) trên 825 cá thể cá bị nhiễm bệnh (tỷ lệ nhiễm chung đạt 22,09%); định danh thành công 37 loài thuộc 2 lớp (Enoplea, Chromadorea), 2 bộ (Trichinellida, Rhabditida), 9 họ và 21 giống. Công trình đã phát hiện 19 loài bổ sung mới cho khu hệ giun tròn Việt Nam, 1 loài tiềm năng mới cho khoa học thuộc giống Cucullanus (Cucullanus (Cucullanus) sp.), đồng thời xác lập 26 loài cá biển thuộc bộ cá Vược lần đầu tiên được ghi nhận là vật chủ mới của giun tròn ký sinh.


Literature Review và Positioning

Lịch sử nghiên cứu giun sán ký sinh ở động vật biển trên thế giới khởi đầu từ các công trình kinh điển của Linnaeus (1758), Rudolphi (1808, 1819), Dujardin (1845) và Schneider (1866). Đến giữa thế kỷ 20, hệ thống phân loại học giun tròn đạt bước tiến lớn với các công trình quy mô của Cobb (1919), Filipjev (1934), Skrjabin (1949-1970) và Yamaguti (1961). Tại khu vực châu Á - Thái Bình Dương, nghiên cứu ký sinh trùng cá biển phát triển mạnh mẽ qua các khảo sát tại Trung Quốc (Wu, 1927; Peng và cs., 2011 với 90 loài giun tròn biển), Philippines (Arthur và cs., 1997 với 20 loài giun tròn), Thái Lan (Purivirojkul, 2004; Thanapon Yooyen và cs., 2011), Pakistan (Akhtar và cs., 2006; Kazmi và cs., 2013 với 90 loài giun tròn), Malaysia (Arthur và cs., 2015 với 24 loài) và Nhật Bản (Layman, 1930; Yamaguti, 1941; Nagasawa và cs., 2007).

Tại Việt Nam, các dẫn liệu đầu tiên về ký sinh trùng cá biển được thực hiện chủ yếu qua các chuyến khảo sát hợp tác Xô - Việt giai đoạn 1961-1987 bởi các nhà khoa học Liên Xô như Oshmarin, Parukhin, Mamaev và Lebedev (1963-1987). Nhóm nghiên cứu đã mổ khám 1.118 cá thể thuộc 53 loài cá biển (trong đó có 888 cá thể của 44 loài bộ cá Vược), phát hiện 453 loài ký sinh trùng với 27 loài giun tròn. Mamaev (1970) khảo sát 541 cá thể của 28 loài cá vịnh Bắc Bộ, phát hiện 11 loài giun tròn; Lebedev (1970) nghiên cứu 582 cá thể cá di cư đại dương, thu được 7 loài giun tròn. Trong nước, các nghiên cứu chuyên sâu về giun tròn cá biển ghi nhận dấu ấn của Lê Văn Hòa và cs. (1972) với việc công bố 2 loài mới thuộc giống Bulbocephalus, Hà Duy Ngọ và cs. (2009), Nguyễn Văn Hà (2011) mô tả loài mới Ascarophis moraveci, Võ Thế Dũng (2010) và Đặng Nguyễn Anh Tuấn và cs. (2015).

                Bộ TRICHINELLIDA                                         Bộ RHABDITIDA
                   Giống Capillaria                         Phân bộ Spirurina       Phân bộ Cucullanina
                                                 Anisakidae  Raphidascarididae Camallanidae Cucullanidae

Trong bức tranh tổng quan y sinh học, y văn tồn tại ba luồng quan điểm tranh biện học thuật lớn:

  1. Tranh biện về bậc phân loại thượng đẳng: Quan điểm truyền thống xếp Nematoda là một lớp thuộc ngành Nemathelminthes (Schneider, 1886) đối lập hoàn toàn với trường phái hiện đại khẳng định Nematoda là một ngành độc lập (Phylum) được phân chia thành hai lớp EnopleaChromadorea dựa trên bằng chứng phát sinh chủng loại phân tử SSU rDNA (Maggenti, 1991; De Ley & Blaxter, 2004).
  2. Tranh biện về phân định địa lý sinh học và gộp taxon: Trong danh lục khu hệ giun tròn cá biển Trung Quốc, Peng và cs. (2011) đã gộp toàn bộ 19 loài giun tròn do các tác giả Liên Xô thu thập tại vịnh Bắc Bộ và biển Đông Việt Nam vào khu hệ biển Trung Quốc. Luận án đã phản biện và tái định vị ranh giới phân bố sinh học độc lập, khẳng định vị thế đặc trưng của khu hệ cá biển thềm lục địa Việt Nam.
  3. Tranh luận về định loại ấu trùng họ Anisakidae: Việc chỉ dựa vào đặc điểm hình thái học quang học đối với ấu trùng giai đoạn 3 (L3) của phức hợp loài Anisakis thường dẫn đến nhận định sai lầm khi quy chụp tất cả mẫu vật về loài Anisakis simplex (gây bệnh ở vùng ôn đới). Bằng việc giải trình tự vùng ITS2-rDNA, luận án đã chứng minh sự hiện diện ưu thế của Anisakis typica Diesing, 1860 tại vùng biển nhiệt đới ven bờ Việt Nam.

So với các công trình quốc tế như nghiên cứu của Peng và cs. (2011) tại Trung Quốc (90 loài giun tròn trên quần xã cá biển rộng) hay Kazmi và cs. (2013) tại Pakistan (90 loài), luận án của Hoàng Văn Hiền (2019) có định vị vượt trội khi đào sâu chuyên biệt vào một bộ cá duy nhất (Perciformes) với cỡ mẫu mổ khám thực địa kỷ lục (3.735 cá thể), tạo lập hệ thống dữ liệu đối sánh hình thái và dịch tễ học chuẩn xác dọc theo chiều dài bờ biển nhiệt đới Đông Nam Á.


Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án củng cố và mở rộng khung phân loại học tiến hóa của De Ley & Blaxter (2004) bằng việc chứng minh cấu trúc phân loại của khu hệ giun tròn ký sinh cá biển nhiệt đới bao gồm hai nhánh tiến hóa song song: lớp Enoplea (bộ Trichinellida, họ Capillariidae) và lớp Chromadorea (bộ Rhabditida với các họ Anisakidae, Raphidascarididae, Cystidicolidae, Spinitectidae, Philometridae, Camallanidae, Cucullanidae, Quimperiidae).

Về mặt lý thuyết sinh thái ký sinh trùng, luận án xác thực Lý thuyết chu trình truyền qua lưới thức ăn (Trophic Transmission Theory) của Malakhov (1986) và Cox (1993). Bằng chứng thực nghiệm từ dữ liệu mổ khám chứng minh các loài giun tròn ký sinh sinh học (bio-helminths - phát triển gián tiếp qua vật chủ trung gian là giáp xác, thân mềm hoặc cá mồi) chiếm ưu thế tuyệt đối (100% các loài phát hiện) so với giun tròn ký sinh địa học (geo-helminths).

Luận án xác lập 3 mệnh đề lý thuyết cốt lõi:

  • Mệnh đề 1 (Đặc hiệu vật chủ hẹp - Stenoxenous specificity): Các loài giun tròn trưởng thành có cấu tạo miệng và bao cơ quan sinh dục phức tạp (như Philometroides atropi, Spinitectus echenei, Buckleyella buckleyi) thể hiện tính gắn kết chặt chẽ với một hoặc vài loài vật chủ thuộc cùng một họ cá chuyên biệt.
  • Mệnh đề 2 (Đặc hiệu vật chủ rộng - Euryxenous specificity): Các loài giun tròn thuộc họ Raphidascarididae (điển hình là Hysterothylacium aduncum) có khả năng thích nghi chuyển dịch vật chủ cực cao, ký sinh trên 13 loài cá thuộc nhiều họ khác nhau nhờ khả năng xâm nhập qua các mắt xích dinh dưỡng phong phú.
  • Mệnh đề 3 (Cơ chế tích lũy ấu trùng - Paratenic accumulation): Tỷ lệ nhiễm ấu trùng giun tròn họ Anisakidae tăng vọt ở các loài cá vược ăn thịt bậc cao (như Trichiurus lepturus, Priacanthus hamrur, Siganus canaliculatus) do hiện tượng tích lũy sinh học qua chuỗi thức ăn đáy và cá tầng nổi.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp ba trường phái phương pháp luận:

  1. Hình thái học chẩn loại kinh điển (Morphotaxonomy): Dựa trên khóa phân loại của Yamaguti (1961), Anderson và cs. (2009), Moravec (2001, 2004), phân tích chi tiết cấu trúc môi, cánh đầu, xoang miệng, tỷ lệ thực quản cơ/tuyến, cấu trúc gai giao phối (spiculum), gai điều chỉnh (gubernaculum) và núm sinh dục đuôi đực.
  2. Kỹ thuật hiển vi điện tử quét (Ultrastructural SEM): Ứng dụng quy trình Fagerholm (1982) phân tích siêu cấu trúc bề mặt cutin, núm cảm giác đầu và gai cutin thân ở cấp độ nanomet.
  3. Phân loại học phân tử (Molecular Phylogenetics): Khuyếch đại và giải trình tự vùng gen ITS2-rDNA (Zhu và cs., 2000), phân tích khoảng cách di truyền Kimura-2-Parameter (K2P) và mô hình hóa cây phát sinh loài Maximum Likelihood (ML) trên phần mềm MEGA v.6 (Tamura và cs., 2013).

Điều kiện biên (Boundary conditions) của khung phân tích: Giới hạn trên các loài cá thuộc bộ cá Vược sinh sống tại vùng nước biển ven bờ có độ sâu từ 30 m trở vào (đối với vịnh Bắc Bộ và Nam Bộ) và từ 50 m trở vào (đối với vùng biển miền Trung), không áp dụng cho quần xã cá tầng khơi đại dương ngoài thềm lục địa.


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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng khoa học (Positivism) và chủ nghĩa hiện thực phân loại tích hợp (Integrative Taxonomic Realism). Thiết kế nghiên cứu cắt ngang đa tầng (Multi-level Cross-sectional Design) được triển khai đồng bộ qua 3 cấp độ:

  • Cấp độ 1 (Không gian địa lý): 9 tỉnh ven biển đại diện cho 3 vùng sinh thái hải văn đặc trưng: Miền Bắc (Quảng Ninh, Hải Phòng, Nam Định); Miền Trung (Nghệ An, Quảng Bình, Thừa Thiên Huế, Khánh Hòa); Miền Nam (Kiên Giang, Bạc Liêu).
  • Cấp độ 2 (Hệ sinh thái vật chủ): Khảo sát bao trùm các sinh cảnh rừng ngập mặn, vùng cửa sông, bãi triều đáy bùn cát và rạn san hô ngầm.
  • Cấp độ 3 (Bậc phân loại vật chủ): 3.735 cá thể cá thuộc 129 loài, 38 họ của bộ cá Vược. Cỡ mẫu cho các loài phổ biến đảm bảo tối thiểu 10 cá thể/loài/tỉnh.

Quy trình nghiên cứu rigorous

  • Quy trình mổ khám ký sinh trùng: Thực hiện theo phương pháp mổ khám ký sinh trùng toàn diện của Viện sĩ Skrjabin (1928). Cơ thể cá được giải phẫu tỉ mỉ từ xoang miệng, ống tiêu hóa (dạ dày, ruột trước, ruột sau, manh tràng), gan, mật, tim, thận, xoang cơ thể đến các lớp cơ và mắt dưới kính lúp lập thể; kết hợp phương pháp gạn lọc liên tục để thu thập triệt để các cá thể giun kích thước hiển vi.
  • Quy trình định hình và bảo quản mẫu: Mẫu giun sống được làm giãn cơ bằng nước nóng 60-70°C; sau đó chia làm hai phần: phần lớn cố định trong dung dịch Formalin 4% (đối với mẫu nhỏ) và Formalin 10% (đối với mẫu lớn) để nghiên cứu hình thái; phần còn lại cố định trong cồn tinh khiết 70-96% để phân tích sinh học phân tử.
  • Kỹ thuật làm tiêu bản hình thái:
    • Tiêu bản tạm thời: Làm trong mẫu bằng hỗn hợp dung dịch Glycerin + Axit Lactic + Nước cất theo tỷ lệ thể tích chuẩn 1:1:1.
    • Tiêu bản cố định lâu dài: Thực hiện theo quy trình bay hơi chậm của De Grisse (1969) qua ba bước dung dịch (Dung dịch I: Formalin 4% + Glycerin 99:1; Dung dịch II: Cồn 96% + Glycerin 95:5; Dung dịch III: Cồn 96% + Glycerin 50:50) trong tủ sấy 40°C, sau đó bao gắn tiêu bản trên lam kính bằng sáp ong tinh khiết.
  • Kỹ thuật hiển vi điện tử quét (SEM): Mẫu giun cố định được xử lý chuyển tiếp sang dung dịch Tetra-oxit Osimi (OsO4) 1% trong 1 giờ ở 40°C; khử nước qua thang nồng độ cồn tăng dần (80%, 96%, 100%); sấy khô tới hạn bằng khí CO2 lỏng; phủ màng vàng kim loại và chụp ảnh siêu cấu trúc bề mặt trên kính hiển vi điện tử quét FEIX L30 ESEMFEG (Fagerholm, 1982; Lee, 1992).
  • Quy trình tam giác giác đấu chuyên gia quốc tế (Investigator Triangulation): Kết quả giám định hình thái và mô tả loài được đối chiếu, thẩm định trực tiếp với các chuyên gia giun tròn hàng đầu thế giới: GS. Sergei Spiridonov (Viện Hàn lâm Khoa học Nga), GS. František Moravec (Viện Hàn lâm Khoa học CH Séc), TS. David González-Solís (Đại học Frontera, Mexico) và TS. Shokoofeh Shamsi (Đại học Charles Sturt, Australia).

Data và phân tích

  • Quy trình sinh học phân tử vùng ITS2-rDNA: Tách chiết DNA tổng số bằng DNeasy Tissue Kit (QIAGEN, Đức). Phản ứng PCR khuếch đại vùng ITS2-rDNA sử dụng cặp mồi chuyên biệt NC5 (Mồi xuôi: 5'-GTAGGTGAACCTGCGGAAGGATCATT-3') và NC2 (Mồi ngược: 5'-TTAGTTTCTTCCTCCGCT-3') (Zhu và cs., 2000). Thành phần phản ứng 40 µl gồm: 20 µl Master Mix (2x), 1 µl Primer F (10 pmol), 1 µl Primer R (10 pmol), 2 µl DNA khuôn (10 ng) và 16 µl H2O cất khử ion. Chu trình nhiệt: Biến tính ban đầu 96°C trong 180 giây; 35 chu kỳ [95°C trong 30 giây, bắt cặp 48-51°C trong 40 giây, kéo dài 72°C trong 45 giây]; kéo dài kết thúc 72°C trong 10 phút. Điện di kiểm tra trên gel Agarose 1,0% nhuộm Ethidium Bromide, giải trình tự Sanger tự động tại Macrogen Inc. (Seoul, Hàn Quốc).
  • Phân tích số liệu phân tử và dịch tễ: Dãy trình tự nucleotide được căn chỉnh bằng BLASTn trên cơ sở dữ liệu GenBank (NCBI). Cây phát sinh chủng loại được thiết lập bằng thuật toán Maximum Likelihood (ML) kết hợp mô hình khoảng cách di truyền Kimura-2-Parameter (K2P) với 1.000 lần lặp bootstrap trên phần mềm MEGA v.6 (Tamura và cs., 2013).
  • Chỉ số dịch tễ học:
    • Tỷ lệ nhiễm (Prevalence - TLN %) = $(Số cá nhiễm giun tròn / Tổng số cá mổ khám) \times 100$.
    • Cường độ nhiễm (Intensity - CĐN) = Khoảng giá trị (Min - Max) số cá thể giun tròn ký sinh trên một cá thể vật chủ bị nhiễm.
    • Toàn bộ dữ liệu dịch tễ học được xử lý và phân tích trên phần mềm Microsoft Excel 2007; hình vẽ và sơ đồ hình thái được tái lập trên phần mềm đồ họa vector Adobe Illustrator.

Phát hiện đột phá và implications

Những phát hiện then chốt

Kết quả phân tích 7.951 mẫu giun tròn thu thập từ 3.735 cá thể cá thuộc 129 loài cá Vược đã làm sáng tỏ cấu trúc quần xã giun tròn ký sinh biển ven bờ Việt Nam với các phát hiện cụ thể:

STT Taxon Giun tròn Họ vật chủ tiêu biểu Tỷ lệ nhiễm (%) Cường độ nhiễm (con/cá) Vị trí ký sinh chính Ghi chú khoa học
1 Anisakis typica Carangidae, Priacanthidae, Trichiuridae 4,11 - 33,89% 1 - 18 Xoang cơ thể, màng treo ruột Khẳng định phân tử ITS2; nguy cơ zoonotic
2 Contracaecum osculatum Siganidae, Sparidae, Sciaenidae 5,56 - 100% 1 - 100 Dạ dày, ruột, xoang 100% ở Siganus canaliculatus (CĐN 23-100)
3 Hysterothylacium aduncum Carangidae, Sparidae, Serranidae 1,67 - 58,06% 1 - 30 Ruột non, dạ dày Phổ vật chủ rộng nhất (13 loài cá)
4 Hysterothylacium fabri Uranoscopidae (U. oligolepis) 58,82% 1 - 40 Dạ dày Ghi nhận mới cho khu hệ Việt Nam
5 Hysterothylacium longilabrum Siganidae (Siganus fuscescens) 41,53% 5 - 148 Ruột non Ghi nhận mới cho khu hệ Việt Nam
6 Philometra sciaenae Sciaenidae (Nibea albiflora) 42,86% 1 - 23 Tuyến sinh dục, xoang Gây teo tuyến sinh dục cá
7 Cucullanus (Cucullanus) sp. Mullidae, Priacanthidae 13,33 - 25,00% 1 - 6 Ruột Loài tiềm năng mới cho khoa học
  1. Xác lập đa dạng cấu trúc thành phần loài: Định danh 37 loài thuộc 21 giống, 9 họ, 2 bộ, 2 lớp. Trong đó có 24 loài đã định tên chính xác và 13 loài định tên đến giống. Đặc biệt, luận án đã bổ sung 19 loài giun tròn lần đầu tiên ghi nhận tại Việt Nam (như Contracaecum osculatum, Hysterothylacium fabri, H. longilabrum, Raphidascaris acus, Raphidascaroides nipponensis, Prospinitectus mollis, Cucullanus rastrelligeri, Cucullanus truttae...) và phát hiện 1 loài (Cucullanus (Cucullanus) sp.) mang đầy đủ đặc trưng hình thái của một loài mới cho khoa học.
  2. Xác lập 26 loài cá vật chủ mới: Bổ sung 26 loài cá biển thuộc bộ cá Vược lần đầu tiên được ghi nhận là vật chủ của các loài giun tròn trong khu hệ sinh học Việt Nam (tiêu biểu như Siganus fuscescens, Uranoscopus oligolepis, Nemipterus hexodon, Dentex tumifrons, Evynnis cardinalis...).
  3. Phát hiện dịch tễ nhiễm cực cao ở các họ cá chuyên biệt: Họ cá Dìa (Siganidae) và họ cá Trác (Priacanthidae) ghi nhận mức độ nhiễm ấu trùng và giun trưởng thành cao nhất. Điển hình loài cá Dìa hoa (Siganus canaliculatus) ghi nhận tỷ lệ nhiễm Contracaecum osculatum đạt 100% với cường độ nhiễm 23-100 giun/cá; cá Dìa (Siganus fuscescens) nhiễm tổng cộng 465 giun trưởng thành và 1.443 ấu trùng. Loài cá Trác đỏ (Priacanthus hamrur) nhiễm 31 giun trưởng thành và 400 ấu trùng giun tròn. Ngược lại, một số họ cá như họ cá Bướm (Chaetodontidae), họ cá Mó (Scaridae) có tỷ lệ nhiễm giun tròn bằng 0% hoặc rất thấp, phản ánh tính chọn lọc thức ăn san hô/rong biển không chứa vật chủ trung gian giáp xác.
  4. Giải mã căn nguyên phân tử loài giun gây bệnh Anisakis: Dữ liệu phân tích vùng ITS2-rDNA kết hợp cây phát sinh loài Maximum Likelihood đã chứng minh toàn bộ các mẫu giun thuộc giống Anisakis thu được trên cá Khế (Carangoides malabaricus), cá Nục (Decapterus macarellus), cá Trác (Priacanthus hamrur) và cá Hố (Trichiurus lepturus) tại biển Việt Nam đều là loài Anisakis typica Diesing, 1860 (với độ tương đồng di truyền >99% so với các chủng A. typica trên GenBank), bác bỏ các suy đoán trước đây về sự hiện diện phổ biến của Anisakis simplex ở cá ven bờ nhiệt đới.
  5. Quy luật phân hóa địa lý sinh học Bắc - Trung - Nam: Số lượng loài giun tròn phân bố giảm dần từ Bắc vào Nam: Vùng biển miền Bắc ghi nhận độ đa dạng loài phong phú nhất (29/37 loài), tiếp đến vùng biển miền Trung (23/37 loài) và thấp nhất ở vùng biển Tây Nam Bộ (12/37 loài). Sự suy giảm này tương thích chặt chẽ với sự khác biệt về độ sâu thềm lục địa, sự đa dạng của hệ sinh thái rừng ngập mặn vịnh Bắc Bộ và các dòng hải lưu ven bờ.

Implications đa chiều

  • Ý nghĩa học thuật và lý thuyết: Cung cấp bức tranh phân loại học hoàn chỉnh nhất từ trước đến nay về giun tròn cá biển Việt Nam; lấp đầy khoảng trống dữ liệu trong mạng lưới ngân hàng gen quốc tế GenBank/NCBI đối với khu hệ ký sinh trùng biển Đông Nam Á.
  • Đổi mới phương pháp luận: Thiết lập quy chuẩn kết hợp giải phẫu kinh điển Skrjabin với hiển vi điện tử quét SEM và mã vạch phân tử ITS2, tạo khung phương pháp chuẩn mực cho các nghiên cứu ký sinh trùng học thủy sản tiếp theo.
  • Ứng dụng thực tiễn trong nuôi biển: Cảnh báo sớm các nhóm giun tròn nguy hiểm gây tổn thương tuyến sinh dục (Philometra sciaenae), làm thủng vách ruột (Hysterothylacium longilabrum, Goezia sp.) và gây suy kiệt đàn cá nuôi thương phẩm có giá trị cao như cá Song (Epinephelus spp.), cá Vược (Lates calcarifer), cá Dìa (Siganus spp.).
  • Khuyến nghị chính sách và y tế công cộng: Cung cấp cơ sở khoa học để Cục Thú y, Tổng cục Thủy sản và Cục An toàn thực phẩm ban hành quy chuẩn kiểm soát ký sinh trùng trong chuỗi cung ứng hải sản xuất khẩu; đưa ra khuyến cáo y tế đối với thói quen ăn gỏi cá biển sống nhằm ngăn chặn nguy cơ nhiễm ấu trùng Anisakis typicaContracaecum osculatum.

Limitations và Future Research

Luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu:

  1. Giới hạn độ sâu không gian: Khảo sát mới tập trung ở dải ven bờ (<30 m ở vịnh Bắc Bộ/Nam Bộ và <50 m ở miền Trung), chưa bao quát được các ngư trường xa bờ, vùng biển sâu ngoài thềm lục địa và các quần đảo Hoàng Sa, Trường Sa.
  2. Định danh các ấu trùng chưa hoàn chỉnh: 13 loài mới dừng lại ở bậc định loại giống (sp.) do thu thập được ở giai đoạn ấu trùng hoặc thiếu mẫu vật đực trưởng thành để khảo sát đầy đủ cơ quan giao phối.
  3. Mức độ ứng dụng sinh học phân tử: Do giới hạn kinh phí, phân tích giải trình tự gen ITS2 mới chỉ áp dụng tập trung cho phức hợp loài Anisakis, chưa thực hiện barcoding toàn bộ 37 loài thu thập được.
  4. Tính thời điểm của mẫu dịch tễ: Dữ liệu thu thập theo các đợt khảo sát thực địa từ 2014 đến 2017 chưa phản ánh liên tục biến động mùa vụ theo chuỗi thời gian nhiều năm dưới tác động của biến đổi khí hậu biển.

Chương trình nghiên cứu tiếp nối trong giai đoạn tới:

  • Ứng dụng giải trình tự gen thế hệ mới (NGS) và đa chỉ thị phân tử (COI, 18S, 28S rDNA) để hoàn thiện mã vạch DNA toàn bộ các loài giun tròn biển Việt Nam.
  • Giải mã chi tiết vòng đời sinh học trong phòng thí nghiệm, xác định chính xác các loài giáp xác (Copepoda, Amphipoda) đóng vai trò vật chủ trung gian truyền bệnh.
  • Mở rộng đối tượng nghiên cứu sang các bộ cá kinh tế khác như bộ Cá Thu ngừ (Scombriformes), bộ Cá Trích (Clupeiformes) và bộ Cá Bơn (Pleuronectiformes).
  • Đánh giá cơ chế bệnh học tế bào (histopathology) và đáp ứng miễn dịch của cá nuôi lồng bè khi bị nhiễm giun tròn mật độ cao.

Tác động và ảnh hưởng

  • Tác động học thuật: Công trình là tài liệu tham khảo nền tảng phục vụ biên soạn bộ sách "Động vật chí Việt Nam - Chuyên khảo Giun tròn ký sinh ở cá biển". Kết quả nghiên cứu đã công bố trên các tạp chí chuyên ngành uy tín trong nước và quốc tế, thu hút trích dẫn từ cộng đồng ký sinh trùng học khu vực Tây Thái Bình Dương.
  • Tác động ngành thủy sản: Giúp các trung tâm giống và doanh nghiệp nuôi biển nhận diện chính xác tác nhân gây bệnh đường ruột, từ đó áp dụng các giải pháp quản lý chất lượng nước, khử trùng thức ăn tươi sống nhằm giảm thiểu tỷ lệ hao hụt cá giống và cá thương phẩm.
  • Tác động chính sách: Cung cấp bằng chứng thực chứng giúp Bộ Nông nghiệp và Phát triển Nông thôn xây dựng hàng rào kỹ thuật kiểm dịch ký sinh trùng trong xuất nhập khẩu thủy hải sản, đáp ứng các tiêu chuẩn khắt khe của thị trường EU, Hoa Kỳ và Nhật Bản.
  • Lợi ích xã hội và sức khỏe cộng đồng: Nâng cao nhận thức của cộng đồng ngư dân và người tiêu dùng về vệ sinh an toàn thực phẩm, giảm thiểu nguy cơ bùng phát các ca bệnh viêm loét dạ dày - ruột cấp tính do ấu trùng giun tròn biển gây ra.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học: Tiếp cận hệ thống dữ liệu mô tả hình thái học chi tiết, thông số đo đạc chuẩn tắc (bảng số đo gai giao phối, tỷ lệ thực quản, kích thước trứng) và phả hệ phát sinh loài để phát triển các đề tài chuyên sâu.
  • Các nhà khoa học phân loại học và sinh thái học: Kế thừa danh lục 37 loài giun tròn, 129 loài cá vật chủ và các ghi nhận phân bố địa lý sinh học để nghiên cứu tiến hóa và đa dạng sinh học biển Đông.
  • Kỹ sư và doanh nghiệp nuôi trồng thủy sản: Nắm bắt quy luật nhiễm bệnh theo loài cá (đặc biệt là cá Dìa, cá Đù, cá Lượng, cá Trác) để xây dựng phác đồ phòng trị bệnh giun tròn hiệu quả trong các hệ thống nuôi biển quy mô công nghiệp.
  • Cơ quan quản lý an toàn thực phẩm và dịch tễ thú y: Sử dụng dữ liệu tỷ lệ và cường độ nhiễm để khoanh vùng nguy cơ, thiết lập quy trình kiểm tra chỉ tiêu ký sinh trùng trên hải sản tươi sống lưu thông trên thị trường.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của công trình là gì và đã mở rộng lý thuyết nào? Trả lời: Luận án đã tích hợp thành công cấu trúc phân loại của 37 loài giun tròn ký sinh cá biển Việt Nam vào Hệ thống phát sinh chủng loại hiện đại của De Ley & Blaxter (2004). Công trình mở rộng Lý thuyết lây truyền qua lưới thức ăn (Trophic Transmission) của Malakhov (1986) bằng việc chứng minh cơ chế phân hóa đặc hiệu vật chủ: các loài giun họ Raphidascarididae (Hysterothylacium) thích nghi phổ vật chủ rộng thông qua đa dạng con mồi trung gian, trong khi các họ Philometridae, Spinitectidae chuyên hóa hẹp theo bậc phân loại của cá vật chủ.

  2. Tính đổi mới trong phương pháp nghiên cứu khi so sánh với các nghiên cứu trước đây? Trả lời: Khác biệt với các đợt khảo sát đơn thuần hình thái quang học trước năm 1987 của các tác giả Liên Xô (Parukhin, Mamaev) hay các công trình cục bộ của Lê Văn Hòa (1972) và Võ Thế Dũng (2010), luận án là công trình đầu tiên tại Việt Nam kết hợp đồng bộ "tam giác phương pháp": Mổ khám toàn diện Skrjabin (1928) trên cỡ mẫu lớn (3.735 cá thể) + Hiển vi điện tử quét SEM ESEMFEG chụp cấu trúc nano + Giải trình tự phân tử vùng ITS2-rDNA phân tích cây phát sinh loài ML.

  3. Phát hiện nào gây bất ngờ nhất và có dữ liệu chứng minh ra sao? Trả lời: Phát hiện bất ngờ nhất là loài cá Dìa hoa (Siganus canaliculatus) và cá Dìa (Siganus fuscescens) - loài cá ăn thực vật/rong biển - lại có cường độ và tỷ lệ nhiễm giun tròn cao nhất toàn bộ nghiên cứu. Cụ thể, 100% cá Dìa hoa nhiễm Contracaecum osculatum với cường độ 23-100 con/cá; cá Dìa nhiễm tới 1.443 ấu trùng giun tròn. Điều này chứng minh cá Dìa trong tự nhiên thường xuyên ăn kèm các giáp xác nhỏ bám trên rong biển, biến chúng thành vật chủ tích lũy ấu trùng khổng lồ. Đồng thời, nghiên cứu khẳng định loài giun gây bệnh Anisakis ở cá ven bờ biển Việt Nam là Anisakis typica, không phải Anisakis simplex.

  4. Luận án có cung cấp quy trình nhân bản nghiên cứu (Replication Protocol) hoàn chỉnh không? Trả lời: Luận án mô tả chi tiết từng bước quy trình kỹ thuật: tỷ lệ dung dịch làm trong Glycerin - Axit Lactic (1:1:1), quy trình rút nước 3 bước của De Grisse (1969), quy trình cố định mạ vàng SEM (Fagerholm, 1982), trình tự cặp mồi NC5/NC2, nồng độ hóa chất PCR, chu trình nhiệt luân nhiệt và thuật toán xử lý MEGA v.6, cho phép các phòng thí nghiệm quốc tế tái lập chính xác toàn bộ thí nghiệm.

  5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào? Trả lời: Định hướng 10 năm tới bao gồm: (1) Thiết lập ngân hàng mã vạch DNA toàn diện cho giun ký sinh thủy sản Việt Nam; (2) Giải mã trọn vẹn chu trình truyền lây tự nhiên của các loài giun tròn zoonotic; (3) Nghiên cứu sản xuất chế phẩm sinh học trị giun an toàn cho cá nuôi biển xuất khẩu; (4) Xây dựng hệ thống giám sát dịch tễ học ký sinh trùng biển quốc gia thích ứng với biến đổi khí hậu.


Kết luận

  1. Xác lập hệ thống dữ liệu hoàn chỉnh: Hoàn thành giám định 37 loài giun tròn thuộc 21 giống, 9 họ, 2 bộ, 2 lớp ký sinh trên 129 loài cá thuộc bộ cá Vược tại vùng biển ven bờ Việt Nam từ 7.951 mẫu giun thu thập.
  2. Bổ sung 19 ghi nhận loài mới: Phát hiện và mô tả 19 loài giun tròn lần đầu tiên ghi nhận cho khu hệ sinh học Việt Nam, trong đó có 1 loài tiềm năng mới cho khoa học (Cucullanus (Cucullanus) sp.).
  3. Mở rộng danh lục 26 loài vật chủ: Xác lập 26 loài cá biển thuộc bộ cá Vược lần đầu tiên được ghi nhận là vật chủ của giun tròn ký sinh tại Việt Nam.
  4. Giải mã căn nguyên phân tử loài Anisakis typica: Chứng minh bằng giải trình tự ITS2-rDNA loài Anisakis typica chiếm ưu thế tuyệt đối ở cá biển ven bờ, cung cấp luận cứ khoa học then chốt cho công tác dịch tễ y tế và an toàn thực phẩm.
  5. Làm sáng tỏ quy luật dịch tễ học và địa lý sinh học: Chứng minh tỷ lệ nhiễm chung đạt 22,09%, chỉ ra mức độ nhiễm cao đột biến ở họ cá Dìa (Siganidae), cá Trác (Priacanthidae), cá Lượng (Nemipteridae) và xác lập quy luật suy giảm độ đa dạng loài theo hướng Bắc - Nam.
  6. Mở ra các hướng nghiên cứu liên ngành: Đặt nền móng cho 3 hướng nghiên cứu tương lai gồm Ký sinh trùng học phân tử đại dương (Marine Molecular Parasitology), Dịch tễ học bệnh truyền lây từ hải sản (Seafood-borne Zoonoses Epidemiology) và Công nghệ kiểm soát dịch bệnh nuôi biển bền vững (Sustainable Mariculture Health Management).