Tổng quan về luận án

Nghiên cứu sinh thái học quần thể và sinh học bảo tồn các loài thú biển nguy cấp là một trong những trụ cột then chốt nhằm duy trì tính toàn vẹn của đa dạng sinh học biển toàn cầu. Trong số các loài thú biển phân bố tại vùng Ấn Độ Dương - Tây Thái Bình Dương, cá heo Ông Sư hay cá heo Irrawaddy (Orcaella brevirostris Owen in Gray, 1866) thuộc họ Delphinidae, bộ Cetartiodactyla, là đối tượng đặc biệt nhạy cảm trước các tác động nhân sinh và biến đổi môi trường. Loài này đã được Liên minh Bảo tồn Thiên nhiên Quốc tế (IUCN) xếp hạng Nguy cấp (Endangered - EN) trên phạm vi toàn cầu, trong khi toàn bộ 5 quần thể phụ phân bố biệt lập tại các hệ thống sông và đầm phá nước lợ (sông MêKông, sông Mahakam, sông Ayeyarwady, hồ Songkhla và vịnh Malampaya Sound) đều đối mặt với nguy cơ tuyệt chủng ở mức Cực kỳ nguy cấp (Critically Endangered - CR).

Khoảng trống nghiên cứu (Research Gap): Tại Việt Nam, các tư liệu khoa học về loài cá heo Ông Sư trong lịch sử ghi nhận chủ yếu mang tính chất rời rạc và mô tả định tính. Kể từ các quan sát ban đầu của Llozen (1970s) tại lưu vực sông MêKông và quần đảo Hải Tặc, các ghi nhận tiếp theo chủ yếu là các cá thể ngẫu nhiên mắc lưới ngư dân hoặc dạt bờ như ghi nhận của Đào Tấn Hổ và cộng sự (2002), Beasley và cộng sự (2009). Mặc dù các công trình khảo sát bước đầu của Nguyễn Thị Nga (2003, 2009) và Vũ Long cùng cộng sự (2015, 2017) đã xác nhận sự tồn tại của loài tại vùng biển Tây Nam Bộ, nhưng y văn trong nước hoàn toàn thiếu vắng một nghiên cứu liên ngành toàn diện, tích hợp dài hạn giữa sinh thái học không gian (spatial ecology), kích thước và động thái quần thể (population dynamics), di truyền học phân tử (molecular phylogenetics), cùng với các chỉ tiêu sinh lý - huyết học (hematology) và sinh hóa máu (biochemistry), dinh dưỡng và tập tính học (ethology) trong điều kiện ex-situ (nuôi nhốt có kiểm soát).

Câu hỏi nghiên cứu và Giả thuyết khoa học:

  • RQ1: Cấu trúc không gian, đặc điểm môi trường thủy lý - thủy hóa và vùng phân bố sinh cảnh thích hợp của quần thể O. brevirostris tại vùng biển Tây Nam Bộ được định hình bởi những yếu tố sinh thái then chốt nào?
  • RQ2: Kích thước quần thể, mật độ cá thể và cấu trúc đàn của cá heo Ông Sư biến động ra sao theo thời gian tại quần đảo Bà Lụa, huyện Kiên Lương, tỉnh Kiên Giang?
  • RQ3: Vị trí phân loại phân tử của quần thể cá heo Kiên Giang trên cây phát sinh chủng loại ti thể (gen cytochrome b) có quan hệ di truyền như thế nào so với các quần thể O. brevirostris và các loài thuộc họ Delphinidae trên thế giới?
  • RQ4: Các chỉ số chuẩn sinh lý, sinh hóa máu, động thái tiêu thụ năng lượng/thức ăn và tập tính hô hấp của loài trong điều kiện nuôi nhốt biểu hiện những quy luật sinh học nào để phục vụ công tác cứu hộ và bảo tồn chuyển vị?
  • H1: Độ sâu thủy vực và khoảng cách tới đường bờ/các đảo đá vôi là hai biến số môi trường chi phối chính đến xác suất hiện diện sinh cảnh của cá heo Ông Sư tại vùng biển Tây Nam Bộ.
  • H2: Quần thể cá heo Ông Sư tại vùng biển quần đảo Bà Lụa duy trì quy mô cục bộ ổn định nhưng chịu áp lực cao từ hoạt động khai thác thủy sản và giao thông thủy.
  • H3: Trình tự gen ti thể cytochrome b của quần thể Kiên Giang khẳng định quan hệ gần gũi với phân họ Delphininae và tương đồng cao với các quần thể ven biển Đông Nam Á.

Khung lý thuyết nền tảng (Theoretical Framework): Luận án vận dụng sự tích hợp của Lý thuyết Ổ sinh thái (Ecological Niche Theory - Hutchinson, 1957), Lý thuyết Quần thể biến động / Quần thể phân mảnh (Metapopulation Theory - Levins, 1969), Lý thuyết Tìm kiếm thức ăn tối ưu (Optimal Foraging Theory - MacArthur & Pianka, 1966), Lý thuyết Lịch sử vòng đời (Life History Theory - Stearns, 1992) và Hệ thống học phát sinh loài phân tử (Phylogenetic Systematics - Hennig, 1966).

Đóng góp đột phá và Quy mô nghiên cứu: Luận án của nghiên cứu sinh Cù Nguyên Định (2022) đại diện cho công trình tiến sĩ chuyên khảo dài hạn đầu tiên tại Việt Nam giải quyết đồng bộ từ sinh thái học ngoài thực địa (in-situ) đến sinh thái - sinh lý học thực nghiệm (ex-situ). Dữ liệu nghiên cứu trải dài suốt 10 năm (2008 - 2018), trong đó dữ liệu sơ cấp thu thập liên tục từ 2011 đến 2018 tại vùng biển quần đảo Bà Lụa (diện tích trên 200 km² với khoảng 45 đảo đá vôi thuộc Khu dự trữ sinh quyển thế giới Kiên Giang) kết hợp hệ thống trạm thực nghiệm tại các Câu lạc bộ cá heo Âu Lạc (Quảng Ninh), Vinpearl (Nha Trang), Đại Nam (Bình Dương) và Trung tâm Nhiệt đới Việt - Nga.

Literature Review và Positioning

Khảo cứu y văn quốc tế và trong nước phản ánh các dòng nghiên cứu (research streams) đa chiều về chi Orcaella:

  1. Dòng nghiên cứu phân loại học và phát sinh chủng loại phân tử: Vị trí phân loại của Orcaella từng là tâm điểm tranh luận học thuật kéo dài hơn một thế kỷ. Kể từ mô tả đầu tiên của Owen (1866) từ mẫu vật bến cảng Vizagapatam (Ấn Độ), giống Orcaella từng được luân chuyển giữa các họ Orcellidae, Delphinapteridae và Monodontidae (do hình thái đầu tròn không mỏ tương đồng cá heo trắng Beluga Delphinapterus leucas). Đến cuối thế kỷ XX, Gretarsdottir & Arnason (1992), LeDuc và cộng sự (1999) thông qua giải trình tự gen cytochrome b đã chứng minh Orcaella thuộc họ Delphinidae và có quan hệ gần gũi với cá voi sát thủ (Orcinus orca). Bước ngoặt tiếp theo diễn ra khi Beasley, Robertson & Arnold (2005) tách quần thể Úc/Papua New Guinea thành loài mới - cá heo mũi hếch Úc (Orcaella heinsohni). Vilstrup và cộng sự (2011), Caballero và cộng sự (2019) tiếp tục khẳng định sự phân hóa sâu sắc giữa quần thể nước ngọt sông MêKông với các quần thể ven biển Vịnh Thái Lan và đầm Chilika (Ấn Độ).

  2. Dòng nghiên cứu mô hình hóa sinh cảnh và ổ sinh thái: Các nghiên cứu hiện đại chuyển dịch mạnh mẽ từ ghi nhận tọa độ đơn thuần sang mô hình hóa phân bố loài (Species Distribution Models - SDMs). Smith và cộng sự (2008) ứng dụng mô hình Generalized Additive Models (GAM) tại Bangladesh xác định độ sâu và độ mặn giải thích 36% biến động xuất hiện loài. Mahmud và cộng sự (2018) tại Vịnh Brunei, Peter và cộng sự (2016) tại Vịnh Kuching (Sarawak, Malaysia), Verutes và cộng sự (2020, 2021) và Jackson-Ricketts và cộng sự (2020) tại Vịnh Trat (Thái Lan) và Campuchia sử dụng mô hình Maxent và Hurdle chứng minh cá heo phân bố tập trung ở vùng nước nông cận bờ (độ sâu 2 - 10 m, nhiệt độ bề mặt biển SST 27 - 32°C, cự ly tới đảo < 7,5 km). Wang và cộng sự (2020) xác định sinh cảnh tối đa tại Vịnh Thái Lan là 28.930,2 km² với vùng lõi quan trọng chiếm 5.000 km², trong đó độ sâu đóng góp 91,6% vào mô hình dự đoán.

  3. Dòng nghiên cứu ước lượng quy mô quần thể: Ba phương pháp luận chính thống được phát triển gồm: Đếm trực tiếp (Direct counts) áp dụng tại sông MêKông (Baird & Beasley, 2005; 40 cá thể), sông Mahakam (Kreb, 2002; 34 cá thể), sông Ayeyarwady (Smith & Tun, 2007; 59 - 72 cá thể); Đánh dấu - Chụp ảnh nhận dạng vây lưng (Mark-Recapture / Photo-ID) áp dụng bởi Sutaria & Marsh (2011) tại hồ Chilika (107 - 109 cá thể), Kreb và cộng sự (2020) tại Vịnh Balikpapan (45 - 75 cá thể), Beasley (2007), Ryan và cộng sự (2011), Eam và cộng sự (2020) tại sông MêKông (71 - 136 cá thể); và Phương pháp Khoảng cách tuyến khảo sát (Distance Sampling) áp dụng bởi Smith và cộng sự (2004) tại Malampaya Sound (77 cá thể, CV = 27,4%), Smith và cộng sự (2008) tại Vịnh Bengal (5.383 cá thể, CV = 39,5%), Minton và cộng sự (2011) tại Vịnh Kuching (208 cá thể, CV = 29,1%), và Junchompoo và cộng sự (2017) tại Trat - Thái Lan (376 - 528 cá thể).

Vị thế khoa học của luận án (Positioning): Đặt trong tương quan với hai nghiên cứu quốc tế tương đồng nhất tại vùng Vịnh Thái Lan lân cận là công trình của Jackson-Ricketts và cộng sự (2020) tại tỉnh Trat (Thái Lan) và nghiên cứu của Peter và cộng sự (2016) tại Vịnh Kuching (Sarawak, Malaysia), luận án không chỉ tái xác lập mô hình sinh thái không gian cho vùng biển Tây Nam Bộ Việt Nam mà còn tạo nên bước đột phá khi lần đầu tiên cung cấp bộ dữ liệu sinh lý - sinh hóa máu chuẩn mực và mô hình tương quan dinh dưỡng ex-situ của loài, điều mà các công trình quốc tế nói trên chưa thể thực hiện được do giới hạn tiếp cận mẫu sinh học.

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

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

Luận án đóng góp trực tiếp vào việc mở rộng và kiểm chứng các lý thuyết sinh thái học và tiến hóa hiện đại:

  1. Mở rộng Lý thuyết Ổ sinh thái (Hutchinson, 1957): Xác lập các tham số không gian siêu chiều (multidimensional niche hypervolume) của O. brevirostris tại thủy vực nhiệt đới gió mùa Tây Nam Bộ, chỉ ra tính chất phụ thuộc chặt chẽ vào độ sâu nông (3 - 7 m) và cấu trúc đảo chắn sóng, xác định ranh giới sinh thái giới hạn sự phân bố loài.
  2. Đóng góp cho Lý thuyết Di truyền quần thể và Phát sinh loài (Kimura, 1968; Hennig, 1966): Cung cấp bằng chứng phân tử qua trình tự đoạn mã hóa gen cytochrome b (độ dài đọc chuẩn xác lập trên chủng dòng hóa E. coli DH5α-T1), củng cố vị trí phát sinh loài của Orcaella brevirostris trong họ Delphinidae, chứng minh mức độ sai khác di truyền nucleotide nội loài giữa các quần thể ven biển và nội địa.
  3. Bổ sung Lý thuyết Sinh lý học Thích nghi và Dinh dưỡng Năng lượng (Optimal Foraging Theory): Xây dựng phương trình tương quan định lượng giữa tỷ lệ thức ăn tiêu thụ hàng ngày trên khối lượng thân ($FI%$) với các nhóm tuổi và giới tính, bác bỏ giả định tuyến tính đơn giản và xác lập quy luật giảm dần phi tuyến tính của nhu cầu trao đổi chất tương đối theo sự tăng trưởng khối lượng cơ thể.

Mô hình định đề lý thuyết (Theoretical Propositions):

  • Proposition 1 (P1): Xác suất hiện diện của cá heo Ông Sư tỷ lệ nghịch với khoảng cách tới các đảo đá vôi và độ sâu thủy vực vượt quá 10 m.
  • Proposition 2 (P2): Quần thể cá heo vùng biển Bà Lụa biểu hiện cấu trúc phân hóa di truyền mtDNA đồng nhất với nhóm cá heo ven biển Đông Nam Á nhưng tách biệt rõ với phân quần thể nước ngọt sông MêKông.
  • Proposition 3 (P3): Hiệu suất đồng hóa và tỷ lệ tiêu thụ dinh dưỡng hàng ngày ở cá thể non ($IM$) cao hơn có ý nghĩa thống kê so với cá thể trưởng thành ($M$), tương ứng với sự biến thiên của các chỉ số sinh hóa máu (Urea, Creatinine, Men gan).

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

Khung phân tích của luận án kết hợp 3 nhánh tiếp cận liên hoàn:

  • Nhánh Phân tích Sinh thái Cảnh quan & Môi trường (Macro-level): Đo lường 25 thông số lý - hóa tầng nước mặt và trầm tích đáy bằng các thiết bị chuyên dụng (Van Dorn Sampler, Ekman Grab), chuẩn hóa theo TCVN, ISO và SMEWW, làm đầu vào cho mô hình hóa vùng phân bố thích hợp.
  • Nhánh Động thái Quần thể (Meso-level): Kết hợp điều tra thực địa trên ca-nô chuyên dụng 200CV và thuyền cao su 520AR, áp dụng song song phương pháp Khoảng cách (Distance Sampling) và Đánh dấu - Chụp ảnh nhận dạng vây lưng (Photo-ID Mark-Recapture).
  • Nhánh Sinh học Phân tử, Sinh lý & Dinh dưỡng Thực nghiệm (Micro-level): Tách chiết DNA bằng Kit Qiagen, khuếch đại PCR với cặp mồi chuyên biệt OBREF/OBRER trên máy PTC-100, dòng hóa plasmid vector pCR2.1, cắt enzyme giới hạn EcoRI, giải trình tự gen; kết hợp lấy mẫu máu tĩnh mạch vây đuôi định kỳ để xét nghiệm 15 chỉ số huyết học và sinh hóa lâm sàng.

Điều kiện biên (Boundary Conditions): Kết quả khảo sát thực địa áp dụng trực tiếp cho vùng biển ven bờ có chế độ bán nhật triều thuộc Quần đảo Bà Lụa - Vịnh Rạch Giá - Hà Tiên; các thông số sinh lý - dinh dưỡng áp dụng chuẩn xác cho các cá thể trong điều kiện nuôi giữ nhân tạo tại vùng nhiệt đới với nguồn thức ăn tiêu chuẩn từ cá xương và thân mềm chân đầu.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết học Thực chứng (Positivism) kết hợp Thực nghiệm Sinh thái học (Empirical Ecological Paradigm), đảm bảo tính khách quan và khả năng tái lập cao. Thiết kế nghiên cứu hỗn hợp (Mixed-methods design) đan kết chặt chẽ giữa quan sát sinh thái học ngoài thiên nhiên và đo đạc sinh lý học trong phòng thí nghiệm.

Thiết kế đa cấp độ (Multi-level Design):

  • Cấp vĩ mô: Toàn bộ vùng biển Tây Nam Bộ và Vịnh Thái Lan (khảo sát không gian phân bố).
  • Cấp trung mô: Quần thể cá heo Ông Sư cư trú tại Quần đảo Bà Lụa, xã Sơn Hải và Bình An, huyện Kiên Lương, tỉnh Kiên Giang (tọa độ các trạm quan trắc điển hình từ $10.33230^\circ\text{N}$).
  • Cấp vi mô: Các cá thể cá heo được lưu giữ an toàn phục vụ nghiên cứu khoa học tại hệ thống lồng bè biển Kiên Lương, Nha Trang ($10\text{ m} \times 10\text{ m} \times 5\text{ m}$) và bể nuôi chuyên dụng trên bờ tại Bình Dương, Quảng Ninh ($25\text{ m} \times 15\text{ m} \times 6\text{ m}$).

Quy trình nghiên cứu rigorous

Chiến lược lấy mẫu và thu thập dữ liệu:

  • Dữ liệu môi trường nước và trầm tích: Mẫu nước mặt thu bằng thiết bị Van Dorn Model 3-1130-G42 (dung tích 3,2 lít); mẫu trầm tích đáy thu bằng gầu Ekman Model 3-196-B12 tại các điểm xuất hiện cá heo trong cả mùa mưa (MM) và mùa khô (MK) năm 2018. Phân tích 25 chỉ tiêu bao gồm: pH (TCVN 6492:2011), Độ mặn (SMEWW 2520B:2012), DO (TCVN 7325:2004), BOD5 (SMEWW 5210B:2012), COD (SMEWW 5220C:2012), các hợp chất nitơ ($N\text{-}NO_3^-, N\text{-}NO_2^-, N\text{-}NH_4^+$, Tổng N), phốt pho ($P\text{-}PO_4^{3-}$, Tổng P), kim loại nặng (Fe, Mn, Cu, Pb, Zn, Cd, As qua US EPA Method 6020A và SMEWW 3125B) và Tổng Coliforms (TCVN 6187-2:1996).

  • Điều tra khảo sát trên biển: Sử dụng ca-nô Yamaha 200CV (vận tốc hành trình 60 - 80 km/h) phối hợp thuyền cao su 520AR, định vị liên tục bằng GPS Garmin OREGON 650, dò tìm luồng cá và đo độ sâu bằng Navman Fish 4100, quan sát trực quan bằng ống nhòm Bosma 1050 ZCY, ghi nhận hình ảnh chất lượng cao phục vụ Photo-ID bằng máy ảnh Canon 600D và máy quay Sony HDR-XR100E.

  • Quy trình sinh học phân tử: Tách chiết DNA mô bằng Kit Qiagen Genomic DNA Extraction; phản ứng chuỗi Polymerase (PCR) khuếch đại gen cytochrome b bằng Kit AccuPower® ProFi Taq PCR Premix (Bioneer, Hàn Quốc) trên máy luân nhiệt PTC-100 (MJ Research, Mỹ); tinh sạch sản phẩm PCR bằng Kit QIAquick PCR Purification (Qiagen); điện di kiểm tra trên thạch agarose 1% với thang chuẩn GeneRuler 1kb DNA Ladder (Sigma, Mỹ); biến nạp sản phẩm vào tế bào khả biến Escherichia coli chủng DH5α-T1 thông qua vector dòng hóa pCR™2.1; cắt kiểm tra plasmid tái tổ hợp bằng enzyme cắt hạn chế EcoRI; giải trình tự và phân tích căn chỉnh chuỗi nucleotide và amino acid suy diễn bằng phần mềm chuyên dụng, so sánh đối chiếu với cơ sở dữ liệu NCBI GenBank.

  • Quy trình sinh lý, sinh hóa và dinh dưỡng: Lấy mẫu máu vô trùng từ tĩnh mạch vây đuôi của các cá thể cá heo thuần dưỡng bằng kim bướm chuyên dụng. Phân tích các chỉ số hồng cầu (RBC), bạch cầu (WBC), huyết sắc tố (Hb), hematocrit (HCT), tiểu cầu (PLT) và các chỉ số sinh hóa máu (Glucose, Cholesterol, Triglyceride, AST, ALT, Urea, Creatinine, Albumin, Protein toàn phần). Theo dõi khẩu phần dinh dưỡng hàng ngày, phân tích thành phần sinh hóa của các loại cá mồi (hàm lượng đạm, béo, khoáng, năng lượng trao đổi), theo dõi nhịp thở và thời gian lặn sâu.

  • Tam giác đạc (Triangulation): Kết hợp tam giác đạc dữ liệu (in-situ vs ex-situ), tam giác đạc phương pháp (Distance sampling vs Photo-ID Mark-recapture vs Ecological niche modeling) và tam giác đạc lý thuyết (sinh thái cảnh quan kết hợp tiến hóa phân tử).

Data và phân tích

Phân tích thống kê và mô hình hóa nâng cao:

  • Sử dụng phần mềm Distance (phiên bản chuyên dụng cho Cetacean survey) để xử lý hàm phát hiện khoảng cách $g(y)$, tính toán hệ số biến thiên ($CV%$), mật độ cá thể ($D$, cá thể/km²), kích thước quần thể kỳ vọng ($\hat{N}$) kèm khoảng tin cậy 95% ($95% CI$).
  • Xây dựng mô hình phân bố loài (SDMs) sử dụng thuật toán Maxent kết hợp kiểm định diện tích dưới đường cong ROC ($AUC$) và chỉ số thống kê True Skill Statistic ($TSS$) để đánh giá độ chính xác dự báo.
  • Dựng cây phát sinh chủng loại di truyền theo phương pháp Maximum Likelihood và Neighbor-Joining với giá trị lặp lại Bootstrap 1.000 lần trên MEGA/BioEdit.
  • Kiểm định sai khác thống kê bằng ANOVA, Kruskal-Wallis và phân tích hồi quy phi tuyến tính trên SPSS/R.

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

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

Luận án đã đạt được 5 nhóm phát hiện khoa học đột phá mang tính hệ thống:

Thứ nhất: Xác lập đặc trưng sinh cảnh và mô hình phân bố không gian tại vùng biển Tây Nam Bộ. Kết quả phân tích thủy hóa tại Quần đảo Bà Lụa cho thấy các thông số môi trường nước mặt nằm trong ngưỡng đặc trưng của sinh cảnh cá heo ven biển nhiệt đới: pH dao động 7,9 - 8,2; độ mặn ổn định 28,5 - 32,4‰; oxy hòa tan DO đạt 5,8 - 6,8 mg/L; các chỉ tiêu kim loại nặng và độc chất đều dưới ngưỡng giới hạn cho phép. Mô hình hóa Maxent xác định độ sâu thủy vực và khoảng cách tới bờ/đảo là hai biến số đóng góp tỷ trọng cao nhất vào sự hiện diện của loài ($AUC > 0,85$; $TSS > 0,7$). Cá heo phân bố tập trung cao nhất tại vùng nước nông quanh các đảo Hòn Heo, Hòn Rễ Lớn, Hòn Rễ Nhỏ với độ sâu lý tưởng từ 3 - 7 m.

Thứ hai: Định lượng quy mô, mật độ và động thái cấu trúc đàn theo thời gian. Dữ liệu điều tra tuyến khoảng cách và định danh hình ảnh xác lập mật độ cá thể $D$ và kích thước quần thể $\hat{N}$ của cá heo Ông Sư tại Quần đảo Bà Lụa. Kích thước đàn phổ biến dao động từ 2 - 10 cá thể, có những đàn lớn tập hợp từ 10 - 20 cá thể vào các thời điểm kiếm ăn tập trung. Tần số bắt gặp đạt đỉnh vào mùa khô (tháng 12 đến tháng 3 năm sau), đặc biệt trong khung giờ từ 11:00 đến 15:00 hàng ngày, trùng hợp với pha triều và quy luật tập trung của các đàn cá mồi ven bờ.

Thứ ba: Giải mã và khẳng định vị trí phân loại phân tử bằng gen ti thể cytochrome b. Luận án đã thành công trong việc tách dòng và giải trình tự hoàn chỉnh đoạn gen cytochrome b từ hai mẫu đại diện ký hiệu OBRE1-VN và OBRE2-VN. Kết quả điện di sản phẩm PCR và cắt hạn chế bằng EcoRI trên vector pCR™2.1 cho băng DNA chuẩn xác. Phân tích gióng hàng chuỗi nucleotide và trình tự amino acid suy diễn khẳng định hai mẫu cá heo Kiên Giang có độ tương đồng 99,2 - 99,8% với các chuỗi Orcaella brevirostris công bố trên NCBI GenBank, hình thành một nhánh đơn dòng (monophyletic clade) vững chắc với nhóm cá heo ven biển Thái Lan và phân tách rõ rệt với phân loài cá heo mũi hếch Úc (Orcaella heinsohni) cũng như các đại diện thuộc họ Delphinidae khác.

Thứ tư: Thiết lập bộ chỉ số tham chiếu huyết học, sinh hóa máu chuẩn mực đầu tiên tại Việt Nam. Nghiên cứu ex-situ đã cung cấp bảng hằng số sinh lý - sinh hóa máu toàn diện phân theo nhóm tuổi (trưởng thành $M$ vs chưa trưởng thành $IM$) và giới tính (đực vs cái). Không có sự sai khác có ý nghĩa thống kê về số lượng hồng cầu (RBC: $3,8 - 4,6 \times 10^{12}/\text{L}$) và nồng độ huyết sắc tố (Hb: $135 - 162\text{ g/L}$) giữa hai giới, nhưng cá thể non biểu hiện hoạt độ enzyme trao đổi chất AST, ALT và chỉ số bạch cầu lympho cao hơn, phản ánh cường độ chuyển hóa và đáp ứng miễn dịch đang trong giai đoạn hoàn thiện.

Thứ năm: Mô hình hóa tiêu thụ dinh dưỡng, năng lượng sinh học và nhịp thở. Khảo sát thực nghiệm dinh dưỡng chỉ ra độ ưa thích thức ăn cao nhất đối với các loài cá nổi và cá sống đáy nhỏ cùng mực ống. Lượng thức ăn tiêu thụ hàng ngày tính theo phần trăm khối lượng thân ($FI% = \text{Lượng thức ăn} / \text{Khối lượng cơ thể} \times 100$) biểu hiện mối tương quan nghịch rõ rệt với tuổi và kích thước thân: cá thể non đạt $FI% = 6,5 - 8,2%/\text{ngày}$, trong khi cá thể trưởng thành giảm xuống còn $3,8 - 5,1%/\text{ngày}$. Về tập tính hô hấp, thời gian lặn thông thường duy trì từ 30 - 60 giây xen kẽ chu kỳ thở 3 lần nhanh liên tiếp; trong trạng thái kích động hoặc lặn sâu né tránh mối đe dọa, cá heo có khả năng nín thở kéo dài tới 12 phút với tốc độ bơi bứt phá đạt 25 km/h.

       BẢNG TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN

Implications đa chiều

  1. Ý nghĩa lý thuyết: Đóng góp cơ sở dữ liệu thực nghiệm quan trọng vào ngân hàng di truyền học thú biển quốc tế; khẳng định tính ứng dụng của lý thuyết ổ sinh thái cảnh quan biển đối với các loài thú biển có tập tính bơi cận bờ.
  2. Đổi mới phương pháp luận: Thiết lập quy trình chuẩn mực tích hợp giữa công nghệ định vị vệ tinh, trắc viễn thám, sinh học phân tử và phân tích sinh lý máu lâm sàng, có thể chuyển giao áp dụng cho nghiên cứu các loài thú biển khác tại Việt Nam như cá heo không vây (Neophocaena phocaenoides) và cá voi minke (Balaenoptera acutorostrata).
  3. Ứng dụng thực tiễn và Quản trị bảo tồn: Cung cấp hướng dẫn kỹ thuật cứu hộ, thuần dưỡng, chăm sóc thú y và dinh dưỡng cho cá heo Ông Sư bị mắc cạn hoặc lưu giữ tại các trung tâm bảo tồn biển; xác lập vùng ranh giới ưu tiên bảo tồn nghiêm ngặt tại Quần đảo Bà Lụa trong quy hoạch phân vùng chức năng Khu dự trữ sinh quyển thế giới Kiên Giang.
  4. Khuyến nghị chính sách: Đề xuất cấm hoàn toàn các hình thức khai thác thủy sản mang tính hủy diệt (lưới cào đáy mắt nhỏ, xung điện, bẫy rập dày đặc) trong hành lang di chuyển của cá heo; quy định tốc độ tối đa cho tàu cao tốc vận tải du lịch (< 10 hải lý/giờ khi qua vùng lõi Bà Lụa); thiết lập mạng lưới cứu nạn thú biển phối hợp giữa cộng đồng ngư dân và kiểm lâm/kiểm ngư.

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  1. Giới hạn mẫu di truyền: Dữ liệu giải trình tự gen cytochrome b mới thực hiện trên 2 mẫu đại diện thu nhận tại chỗ (OBRE1-VN và OBRE2-VN) do tính chất quý hiếm và rào cản đạo đức sinh học trong việc thu mẫu mô xâm lấn trên động vật hoang dã nguy cấp.
  2. Ranh giới không gian khảo sát thực địa: Phạm vi điều tra sinh thái học tập trung chính tại Quần đảo Bà Lụa và phụ cận huyện Kiên Lương, chưa thể bao phủ toàn diện 100% vùng biển mở tiếp giáp vùng biển Hà Tiên - Phú Quốc và biên giới hải phận Campuchia do hạn chế kinh phí và phương tiện hành trình biển xa.
  3. Độ phân giải thời gian của dữ liệu môi trường: Quan trắc môi trường nước và trầm tích mới phản ánh 2 mùa chính (mùa khô và mùa mưa) của năm 2018, chưa theo dõi được các biến động cực đoan của các chu kỳ khí hậu liên năm (ENSO/La Niña).

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Hướng 1: Mở rộng phân tích toàn bộ hệ gen ti thể (Mitogenome) và giải trình tự gen nhân độ phân giải cao (RAD-seq/Microsatellite) trên quy mô lớn hơn từ các nguồn mẫu rụng tự nhiên hoặc mẫu eDNA (Environmental DNA) thu thập trong nước biển.
  • Hướng 2: Triển khai công nghệ gắn máy phát tín hiệu vệ tinh (Satellite telemetry tagging) và thiết bị thu âm thụ động dưới nước (Passive Acoustic Monitoring - PAM) để theo dõi 24/7 đường di chuyển, tập tính phát sóng siêu âm định vị sinh học (echolocation clicks) và mức độ ô nhiễm tiếng ồn âm học đại dương tác động lên cá heo.
  • Hướng 3: Nghiên cứu mô hình tương tác tích hợp đa loài (Multi-species trophic modeling - Ecopath with Ecosim) nhằm lượng hóa sự cạnh tranh nguồn lợi thủy sản giữa cá heo Ông Sư và đội tàu khai thác của ngư dân tại Vịnh Thái Lan.
  • Hướng 4: Đánh giá độc học sinh thái (Ecotoxicology) chuyên sâu thông qua phân tích sự tích lũy sinh học của hạt vi nhựa (microplastics) và các chất ô nhiễm hữu cơ khó phân hủy (POPs, PCBs, DDTs) trong chuỗi thức ăn của cá heo.

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

  • Học thuật: Luận án tạo ra bộ dữ liệu tham chiếu tiêu chuẩn (benchmark dataset) về sinh học - sinh thái Orcaella brevirostris tại Việt Nam, mở ra tiềm năng trích dẫn cao trong các ấn phẩm quốc tế thuộc hệ thống ISI/Scopus về sinh học biển và bảo tồn đa dạng sinh học.
  • Chính sách môi trường: Kết quả luận án đóng vai trò luận cứ khoa học trực tiếp cho Ủy ban Nhân dân tỉnh Kiên Giang, Sở Nông nghiệp & Phát triển Nông thôn và Ban Quản lý Khu dự trữ sinh quyển Kiên Giang trong việc điều chỉnh quy hoạch bảo tồn biển, thành lập phân khu bảo tồn loài và sinh cảnh Bà Lụa.
  • Xã hội và Sinh kế: Thúc đẩy mô hình đồng quản lý nguồn lợi thủy sản (co-management), phát triển du lịch sinh thái quan sát cá heo có trách nhiệm (responsible dolphin-watching ecotourism), tạo sinh kế thay thế bền vững cho cộng đồng dân cư xã đảo Sơn Hải và Bình An.
  • Hợp tác quốc tế: Cung cấp dữ liệu kết nối mạng lưới bảo tồn cá heo Irrawaddy xuyên biên giới vùng Vịnh Thái Lan (Việt Nam - Campuchia - Thái Lan) theo định hướng của Nhóm chuyên gia Thú biển IUCN (IUCN Cetacean Specialist Group).

Đố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 phương pháp luận hoàn chỉnh từ thiết kế khảo sát biển, kỹ thuật xử lý ảnh Photo-ID, đến quy trình dòng hóa sinh học phân tử chuẩn mực.
  • Giảng viên và Nhà nghiên cứu Thú biển: Khai thác dữ liệu di truyền và sinh thái học không gian để so sánh đối chiếu với các quần thể cá heo trong khu vực Ấn Độ Dương - Tây Thái Bình Dương.
  • Bác sĩ thú y, Chuyên gia cứu hộ và Các thủy cung/câu lạc bộ chăm sóc thú biển: Sử dụng trực tiếp bảng tham chiếu huyết học, sinh hóa máu và mô hình tỷ lệ tiêu thụ thức ăn ($FI%$) để chẩn đoán bệnh lý, phục hồi thể lực và thiết lập khẩu phần dinh dưỡng chuẩn xác.
  • Nhà hoạch định chính sách và Cơ quan quản lý nhà nước: Sở Khoa học & Công nghệ, Sở Nông nghiệp & PTNT Kiên Giang có căn cứ chuẩn xác để ban hành chỉ thị bảo vệ thú biển và phân bổ ngân sách bảo tồn hiệu quả.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Ổ sinh thái (Ecological Niche Theory) của Hutchinson bằng cách lượng hóa không gian sinh thái siêu chiều của quần thể cá heo ven biển nhiệt đới chịu ảnh hưởng của chế độ hoàn lưu gió mùa. Luận án đã tích hợp thành công giữa các biến số hải dương học vật lý (độ sâu 3 - 7 m, nhiệt độ 27 - 32°C) với cấu trúc địa mạo đảo đá vôi phân tán, chứng minh ổ sinh thái của loài bị thu hẹp đáng kể vào mùa khô và định hình rõ rệt bởi sự phân bố của con mồi đáy.

2. Đổi mới phương pháp luận của luận án thể hiện thế nào khi so sánh với 2 nghiên cứu quốc tế tiêu biểu? So với nghiên cứu của Jackson-Ricketts và cộng sự (2020) tại Thái Lan (chỉ dùng mô hình Hurdle trên dữ liệu tọa độ đơn thuần) và nghiên cứu của Beasley (2007) trên sông MêKông (chỉ dùng Photo-ID Mark-recapture), luận án đã thực hiện bước đột phá phương pháp luận bằng cách tích hợp tam giác đạc: Distance Sampling khảo sát mật độ ngoài biển + Khuếch đại dòng hóa gen Cyt-b xác định loài + Xét nghiệm sinh hóa/huyết học và thực nghiệm dinh dưỡng ex-situ. Đây là thiết kế phương pháp toàn diện nhất từng được thực hiện cho loài này tại Đông Nam Á.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng số liệu đi kèm là gì? Phát hiện bất ngờ nhất là sự thích nghi dinh dưỡng và hô hấp cực đoan của cá heo trong điều kiện tự nhiên và nuôi giữ. Mặc dù cá heo có nhịp lặn trung bình ngắn (30 - 60 giây), nhưng khi chịu kích thích sợ hãi hoặc truy đuổi, chúng có khả năng mở rộng thời gian lặn sâu lên tới 12 phút (gấp 12 - 24 lần bình thường) với tốc độ bơi vọt đạt 25 km/h. Đồng thời, nhu cầu năng lượng trao đổi chất của cá heo con cao vượt bậc, thể hiện qua tỷ lệ tiêu thụ thức ăn $FI%$ đạt tới 8,2% khối lượng thân/ngày so với mức chỉ 3,8% ở cá thể trưởng thành.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không? Có. Toàn bộ quy trình được chuẩn hóa chi tiết: từ mã hiệu thiết bị (Van Dorn 3-1130-G42, Ekman 3-196-B12, GPS Oregon 650, máy dò Navman 4100), thông số hóa chất và kit tách chiết (Qiagen, Bioneer AccuPower), cặp mồi chuyên biệt OBREF/OBRER, chu trình nhiệt PCR trên máy PTC-100, vector pCR™2.1, tế bào khả biến E. coli DH5α-T1, đến các tiêu chuẩn phương pháp thử nghiệm nước (TCVN, ISO, SMEWW). Mọi phòng thí nghiệm sinh học phân tử và sinh thái biển đều có thể tái lập chính xác.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo ra sao? Chương trình tập trung vào 4 trọng tâm: (1) Ứng dụng giải trình tự eDNA đại dương kết hợp Mitogenomics; (2) Gắn thẻ vệ tinh Telemetry và giám sát âm học thụ động (PAM); (3) Xây dựng mô hình chuỗi thức ăn tích hợp Ecopath with Ecosim; (4) Đánh giá độc học tích lũy sinh học vi nhựa và kim loại nặng trong cơ thể thú biển.

Kết luận

Luận án tiến sĩ của Cù Nguyên Định (2022) đại diện cho một công trình khoa học mẫu mực, giải quyết trọn vẹn và chuyên sâu các mục tiêu nghiên cứu đã đề ra:

  1. Xác lập bức tranh toàn cảnh về sinh thái không gian, môi trường thủy lý - thủy hóa và mô hình phân bố sinh cảnh thích hợp của cá heo Ông Sư (Orcaella brevirostris) tại vùng biển Quần đảo Bà Lụa, Kiên Giang, với độ sâu 3 - 7 m và vùng nước quanh đảo đá vôi là sinh cảnh then chốt.
  2. Cung cấp dữ liệu định lượng chính xác về quy mô, mật độ, cấu trúc đàn và quy luật biến động theo mùa của quần thể cá heo trong suốt giai đoạn điều tra 2008 - 2018.
  3. Giải mã thành công đoạn gen ti thể cytochrome b (OBRE1-VN, OBRE2-VN), khẳng định vững chắc vị trí phát sinh loài của cá heo Kiên Giang thuộc họ Delphinidae trên cây phân loại quốc tế.
  4. Thiết lập lần đầu tiên tại Việt Nam bộ chỉ số chuẩn sinh lý, huyết học, sinh hóa máu lâm sàng phân theo lứa tuổi và giới tính phục vụ trực tiếp cho y học thú biển và công tác cứu hộ.
  5. Định lượng hóa phương trình tiêu thụ dinh dưỡng, nhu cầu năng lượng theo khối lượng cơ thể và ghi nhận chi tiết các đặc trưng tập tính hô hấp của loài.
  6. Đề xuất hệ thống giải pháp bảo tồn đồng bộ, khả thi, gắn liền với phân vùng chức năng Khu dự trữ sinh quyển thế giới Kiên Giang, mở ra hướng đi bền vững cho sự tồn vong của loài động vật nguy cấp quý hiếm này trong hải phận Việt Nam.