Tổng quan luận án

Tu hài (Lutraria rhynchaena Jonas, 1844; tên đồng danh Lutraria philippinarum Deshayes), hay còn gọi là tu hài vòi trắng, là loài động vật thân mềm hai mảnh vỏ có giá trị dinh dưỡng và kinh tế cao thuộc họ Vọp (Mactridae). Tại Việt Nam, nghề nuôi tu hài thương phẩm từng phát triển mạnh tại các tỉnh Quảng Ninh, Hải Phòng, Khánh Hòa với tổng diện tích mặt nước ước tính khoảng 1.000 ha, sản lượng đạt 2.621,6 tấn và giá trị doanh thu trên 200 tỷ đồng mỗi năm. Nhu cầu con giống thương phẩm phục vụ sản xuất ước tính đạt khoảng 100 triệu con giống cấp I mỗi năm. Tuy nhiên, nghề nuôi tu hài đối mặt với sự suy thoái nghiêm trọng do hiện tượng cận huyết trong quá trình sản xuất giống nhân tạo truyền thống, dẫn đến tốc độ sinh trưởng suy giảm, sức chống chịu kém và bùng phát dịch bệnh sưng vòi giai đoạn 2012–2017.

Nghị quyết và chính sách phát triển kinh tế biển (theo Quyết định số 1664/QĐ-TTg ngày 04/10/2021 của Thủ tướng Chính phủ phê duyệt Đề án phát triển nuôi trồng thủy sản trên biển đến năm 2030, tầm nhìn đến năm 2045) đặt ra yêu cầu phải làm chủ công nghệ sản xuất và chọn giống thủy sản bản địa theo hướng tăng trưởng nhanh và thích ứng môi trường. Trong khi các nghiên cứu di truyền học phân tử trên đối tượng tu hài tại Việt Nam và trên thế giới còn phân tán, chủ yếu dừng lại ở việc đánh giá đa dạng di truyền qua hệ gene ty thể hoặc chỉ thị vi vệ tinh (microsatellite), khoảng trống nghiên cứu cốt lõi là sự thiếu vắng một bộ dữ liệu hệ gene tham chiếu hoàn chỉnh, hệ thống mã vạch DNA chuẩn hóa phục vụ giám định loài và các chỉ thị phân tử đơn nucleotide (SNP) liên kết trực tiếp với tính trạng tăng trưởng để ứng dụng vào công tác chọn giống dựa trên chỉ thị phân tử (Marker-Assisted Selection - MAS).

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

  1. Đánh giá được thành phần loài tu hài, tiềm năng phát triển nghề nuôi tu hài của huyện Vân Đồn, tỉnh Quảng Ninh.
  2. Xây dựng được mã vạch DNA phục vụ nhận dạng chính xác tu hài (Lutraria rhynchaena Jonas, 1844).
  3. Giải trình tự hệ gene tu hài và phát triển được bộ chỉ thị SNP để phục vụ xây dựng cơ sở dữ liệu phân tử của tu hài (Lutraria rhynchaena).
  4. Phát triển được một số chỉ thị SNP liên quan đến tăng trưởng ở tu hài góp phần phục vụ công tác chọn tạo giống.

Nội dung nghiên cứu

  1. Điều tra, đánh giá thành phần loài, hiện trạng và tiềm năng nghề nuôi tu hài của huyện Vân Đồn, tỉnh Quảng Ninh.
  2. Nghiên cứu xây dựng DNA barcoding cho tu hài (Lutraria rhynchaena).
  3. Nghiên cứu giải trình tự hệ gene tu hài và phát triển bộ chỉ thị SNP cung cấp cơ sở dữ liệu hệ gene tu hài (Lutraria rhynchaena).
  4. Nghiên cứu phát triển một số chỉ thị SNP liên quan đến tính trạng tăng trưởng ở tu hài.

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

  • Đối tượng nghiên cứu: Tu hài vòi trắng (Lutraria rhynchaena Jonas, 1844) thu thập tại các bãi tự nhiên và cơ sở nuôi tại huyện Vân Đồn, tỉnh Quảng Ninh.
  • Phạm vi không gian:
    • Khảo sát thực địa, điều tra sinh thái và thu thập mẫu sinh học tại vùng biển huyện Vân Đồn, tỉnh Quảng Ninh.
    • Phân tích hình thái, sinh hóa, tách chiết DNA tổng số và RNA tại Phòng Công nghệ sinh học, Trường Cao đẳng Kinh tế, Kỹ thuật và Thủy sản (Từ Sơn, Bắc Ninh).
    • Chuẩn bị thư viện giải trình tự và giải trình tự thế hệ mới (NGS) tại Trung tâm Công nghệ Genomics, Đại học Deakin (Geelong, Victoria, Australia).
    • Thu thập, đối chiếu dữ liệu trình tự tham chiếu của các loài thuộc giống Lutraria trên Ngân hàng Gene quốc tế (GenBank/NCBI) và Hệ thống Dữ liệu Mã vạch Sinh vật (BOLDSYSTEMS).
  • Phạm vi thời gian: Các nội dung khảo sát, thu mẫu và thực nghiệm phân tử được tiến hành từ tháng 3/2018 đến tháng 12/2021.

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 nghiên cứu về phân loại học, tình thái phân bố, kỹ thuật nuôi và chọn giống di truyền ở động vật hai mảnh vỏ. Về mặt phân loại, loài Lutraria rhynchaena Jonas, 1844 thuộc bộ Veneroida, họ Mactridae, phân bố tự nhiên ở các vùng biển ấm Tây Thái Bình Dương, Úc, Trung Quốc, Thái Lan, Philippines và Việt Nam. Tại Việt Nam, quần thể phân bố tự nhiên tập trung tại các vùng nước trong, sạch, nền đáy cát pha mảnh vỏ sinh vật ở vịnh Lan Hạ, đảo Cát Bà (Hải Phòng) đến vịnh Bái Tử Long, huyện Vân Đồn (Quảng Ninh) ở độ sâu 10–15 m, độ mặn ổn định 25–30‰.

Trong công tác chọn giống nhuyễn thể, luận án tổng kết hai khuynh hướng chính: chọn lọc dựa trên kiểu hình và chọn lọc dựa trên kiểu gene. Các nghiên cứu định lượng của Gjedrem (1997, 2005) chỉ ra rằng tính trạng tăng trưởng ở động vật thân mềm có hệ số di truyền từ trung bình đến cao (0,2–0,4), cho phép cải thiện năng suất 10–20% qua mỗi thế hệ chọn giống. Về mặt sinh lý, Koehn và Shumway (1982), Zouros (1987) đã chứng minh mối tương quan thuận giữa mức độ dị hợp tử của các gene mã hóa enzyme với tốc độ trao đổi chất và chuyển hóa protein. Ở mức độ phân tử, các nghiên cứu quốc tế đã định vị nhiều biến dị đơn nucleotide liên kết với tính trạng tăng trưởng ở các loài thân mềm hai mảnh vỏ:

  • Feng và cộng sự (2014) xác định gen PyIGFBP liên quan đến tính trạng sinh trưởng ở sò điệp, phát hiện đa hình SNP vị trí 1054 A > G có tương quan thuận rõ rệt với kích thước vỏ và khối lượng thân.
  • Phân tích gene điều hòa phát triển cơ Myostatin (CfMSTN) trên sò điệp Zhikong (Chlamys farreri), các tác giả đã định danh 4 vị trí SNP, trong đó alen mang kiểu gene TT tại vị trí g-1162 G > T cho giá trị chiều dài vỏ và khối lượng cơ vân cao hơn hẳn kiểu gene GT và GG.
  • Nghiên cứu của Meng và cộng sự trên gene PyE2F3 ở sò cho thấy các vị trí PyE2F3-1 và PyE2F3-2 tương quan trực tiếp đến chiều dài, chiều cao vỏ và khối lượng cơ thể.
  • Trên đối tượng hàu Thái Bình Dương (Crassostrea gigas), phân tích vùng điều hòa phiên mã của gene CgTPH xác định 4 SNP (-904, -522, -272, -262 bp) và 2 haplotype liên kết với tốc độ tăng trưởng, đặc biệt biến dị -272 C/T can thiệp vào quá trình liên kết với các yếu tố phiên mã FoxO và CUX.

Về ứng dụng chỉ thị phân tử và mã vạch DNA, khởi đầu từ các nghiên cứu nền tảng của Hebert và cộng sự (2003) cùng quy chuẩn của Consortium for the Barcode of Life (CBOL, 2004) và Taberlet và cộng sự (2007), đoạn trình tự gene ty thể Cytochrome c Oxidase subunit I (COI, độ dài ~650 bp) và 16S rRNA được sử dụng phổ biến làm hệ thống giám định loài. Liu và cộng sự (2018) đã sử dụng COI và 16S rRNA kết hợp gene nhân (18S, 28S, H3) để phân loại phát sinh chủng loại của 56 loài ngao; Setiamarga và cộng sự (2019) xây dựng mã vạch cho 5 loài sò tại Nhật Bản.

Về di truyền phân tử trên đối tượng tu hài, các nghiên cứu trước đây còn rất giới hạn:

  • Gan và cộng sự (2016) giải trình tự hoàn chỉnh hệ gene ty thể của L. rhynchaena dài 16.927 bp (tỷ lệ A+T chiếm 69%, G+C chiếm 31%, mã hóa 12 protein).
  • Thái Thanh Bình và cộng sự (2016) phân tích đoạn gene 16S rRNA trên 30 mẫu tu hài tại Quảng Ninh, Hải Phòng, Khánh Hòa, phát hiện 4 haplotype với độ đa dạng haplotype là 0,570 ± 0,09; tiếp đó công bố 12 locus microsatellite (2017).
  • Trần Thị Trang và cộng sự (2017) đánh giá đa dạng di truyền quần thể tu hài bằng microsatellite, ghi nhận số alen dao động từ 6 đến 12 (trung bình 9 alen/locus; locus mumLr10 có 12 alen, mumLr6 có 6 alen; dị hợp tử quan sát $H_o$ từ 0,371 đến 0,459; dị hợp tử kỳ vọng $H_e$ từ 0,402 đến 0,598; hệ số cận huyết $F_{is}$ từ 0,058 đến 0,248; mức độ phân hóa quần thể $F_{st}$ từ 0,015 đến 0,121).
  • Nghiên cứu metagenome tu hài tại Cát Bà (2020) công bố hệ gene có kích thước khoảng 542 Mbp lắp ráp thành 2.405 scaffolds với $N_{50}$ đạt 0,9 Mbp và 30.533 gene dự đoán.

Vị trí của luận án: Luận án là công trình tiếp nối có hệ thống, giải quyết khoảng trống về hệ gene tham chiếu của loài Lutraria rhynchaena, thiết lập bộ mã vạch DNA tích hợp hai vùng gene 16S rRNA và COI gắn mã QR, đồng thời lần đầu tiên áp dụng kỹ thuật giải trình tự liên kết vị trí cắt giới hạn (ezRAD-Seq) để sàng lọc và xác thực các chỉ thị SNP phân dị giữa hai nhóm tu hài phân ly về giá trị giống ước tính (EBV) đối với tính trạng tăng trưởng.

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

  • Lý thuyết di truyền số lượng và chọn giống phân tử (MAS): Tính trạng tăng trưởng là tính trạng số lượng chịu sự kiểm soát của nhiều gene (polygenic) kết hợp với ảnh hưởng của các yếu tố ngoại cảnh. Ứng dụng giá trị chọn giống ước tính (Estimated Breeding Value - EBV) thu được từ mô hình phả hệ giúp phân tách chính xác nhóm cá thể có tiềm năng di truyền tăng trưởng nhanh và tăng trưởng chậm.
  • Lý thuyết mã vạch DNA (DNA barcoding): Dựa trên nguyên lý khoảng cách di truyền giữa các loài (interspecific divergence) lớn hơn khoảng cách di truyền nội bộ loài (intraspecific divergence). Trình tự nucleotide chuẩn của vùng gene ty thể (COI và 16S rRNA) được sử dụng làm cơ sở nhận diện phân loại học và truy xuất nguồn gốc sinh học.
  • Lý thuyết phân tích hệ gene và kỹ thuật ezRAD (Enzyme Restriction-site Associated DNA): Kỹ thuật ezRAD sử dụng các enzyme cắt hạn chế phổ biến nhằm phân cắt hệ gene thành các phân đoạn ngắn, gắn adapter chuẩn bị thư viện giải trình tự Illumina mà không đòi hỏi thông tin trình tự mồi ban đầu, hỗ trợ phát hiện các đa hình nucleotide đơn (SNP) trên quy mô toàn hệ gene.

Phương pháp nghiên cứu thực tế

  • Phương pháp điều tra sinh thái học và nghề nuôi:
    • Khảo sát thực địa, điều tra ngư dân và hộ nuôi tại huyện Vân Đồn thông qua phiếu điều tra kết hợp báo cáo thống kê của Phòng Nông nghiệp và Phát triển Nông thôn huyện Vân Đồn.
    • Thu thập số liệu khí tượng thủy văn, chất lượng môi trường nước biển, thành phần động - thực vật phù du theo tài liệu của Dự án Hợp phần SUMA và quan trắc địa phương.
  • Phương pháp sinh học phân tử xây dựng mã vạch DNA:
    • Tách chiết DNA tổng số từ mô cơ màng áo tu hài.
    • Khuếch đại PCR các đoạn gene chỉ thị: đoạn gene 16S rRNA (chiều dài phân tích 437 bp) và đoạn gene COI (chiều dài phân tích 615 bp).
    • Tinh sạch sản phẩm PCR, giải trình tự Sanger hai chiều.
    • Phân tích tương đồng chuỗi bằng công cụ BLAST/Megablast trên cơ sở dữ liệu NCBI. Căn chỉnh đa trình tự bằng phần mềm Clustal Omega.
    • Xây dựng cây phát sinh chủng loại bằng thuật toán Maximum Likelihood (ML) trên phần mềm MEGA X và thuật toán Suy luận Bayes (Bayesian Inference - BI) bằng phần mềm BEAST.
    • Mã hóa dữ liệu trình tự thành mã phản hồi nhanh (QR code) phục vụ giám định nhanh.
  • Phương pháp giải trình tự toàn bộ hệ gene (WGS) và hệ gene phiên mã (Transcriptome):
    • Mẫu WGS: Thu nhận 50 mg mô cơ màng áo của 01 cá thể tu hài vòi trắng trưởng thành; tách chiết DNA bằng kit thương mại, kiểm tra độ tinh sạch và nồng độ bằng điện di gel agarose và đo quang phổ.
    • Mẫu Transcriptome: Thu nhận 50 mg mô tiêu hóa của cá thể tu hài; tách chiết RNA tổng số, tinh sạch hạt từ tính poly-T để thu nhận mRNA.
    • Chuẩn bị thư viện giải trình tự NGS theo quy trình chuẩn của Illumina (TruSeq Library Prep).
    • Giải trình tự trên hệ thống giải trình tự thế hệ mới của Illumina tại Trung tâm Công nghệ Genomics, Đại học Deakin (Australia).
    • Xử lý tin sinh học: Kiểm tra chất lượng đọc bằng FastQC; loại bỏ adapter và đoạn đọc kém chất lượng bằng Trimmomatic; ước tính kích thước hệ gene thông qua phân tích phân bố tần số k-mer; lắp ráp de novo hệ gene và transcriptome; đánh giá độ hoàn chỉnh của bộ lắp ráp bằng công cụ BUSCO (Benchmarking Universal Single-Copy Orthologs).
  • Phương pháp sàng lọc chỉ thị SNP tăng trưởng bằng ezRAD-Seq:
    • Thiết kế mẫu: Thu nhận 60 mẫu mô cơ màng áo (50 mg/mẫu) chia thành 2 nhóm đối lập về tính trạng tăng trưởng:
      • Nhóm tăng trưởng nhanh: 30 cá thể chọn ngẫu nhiên từ 63 gia đình tu hài có giá trị EBV cao nhất.
      • Nhóm tăng trưởng chậm: 30 cá thể chọn ngẫu nhiên từ 43 gia đình tu hài có giá trị EBV thấp nhất.
      • Nguồn gốc: Cung cấp bởi đề tài "Ứng dụng di truyền số lượng chọn giống tu hài sinh trưởng nhanh" do Trường Cao đẳng Kinh tế, Kỹ thuật và Thủy sản nuôi giữ tại Trung tâm giống hải sản Vân Đồn.
    • Tách chiết DNA tổng số, xử lý cắt kép hệ gene bằng hai enzyme giới hạn MboI và Sau3AI.
    • Lựa chọn kích thước phân đoạn DNA mục tiêu từ 350 bp đến 550 bp bằng hạt tinh sạch SPB (Sample Purification Beads).
    • Gắn adapter của bộ kit TruSeq Nano DNA HT (Illumina), tiến hành khuếch đại thư viện bằng phản ứng PCR.
    • Giải trình tự Illumina thu nhận dữ liệu đọc đôi (paired-end reads).
    • Phân tích tin sinh học: Sử dụng hai đường ống phân tích (pipeline) độc lập gồm dDocent và Stacks để căn chỉnh các đoạn đọc lên hệ gene tham chiếu, lọc biến dị chất lượng cao, xác định các locus ngoại lai (outlier loci) chịu chọn lọc phân kỳ giữa hai nhóm.
    • Khảo sát các SNP nằm trong 12 gene chức năng liên quan đến quá trình tăng trưởng.
    • Thiết kế mồi PCR đặc hiệu qua phần mềm Primer3 để nhân đoạn trình tự chứa 11 SNP tiềm năng nhằm kiểm chứng tính đa hình.

Nguồn tư liệu và dữ liệu

  • Dữ liệu sơ cấp: Dữ liệu điều tra nghề nuôi tại Vân Đồn; dữ liệu trình tự thô WGS, Transcriptome và ezRAD-Seq từ 60 mẫu tu hài phân tích tại Đại học Deakin.
  • Dữ liệu thứ cấp: Trình tự nucleotide của các loài thuộc giống Lutraria trên GenBank (NCBI); báo cáo kinh tế - kỹ thuật của Phòng Nông nghiệp và Phát triển Nông thôn huyện Vân Đồn năm 2019; số liệu khảo sát thủy sinh vật của Dự án SUMA.

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

Chương 1: Tổng quan nghiên cứu

Chương 1 tổng thuật các cơ sở lý luận và thực tiễn về đối tượng tu hài Lutraria rhynchaena, bao gồm vị trí phân loại học, đặc điểm hình thái thích nghi với đời sống vùi cát (vỏ hình bầu dục dài 7–12 cm, hai mép vỏ hở, vòi xi phông phát triển kéo dài, khối lượng trung bình 100 g/con), đặc tính dinh dưỡng lọc mùn bã và thực vật phù du, tốc độ di chuyển trong nền đáy 0,8–1,2 m/giờ. Tác giả tổng hợp bức tranh sản xuất thủy sản hai mảnh vỏ toàn cầu theo thống kê của FAO (sản lượng nuôi đạt 41,8 triệu tấn năm 2020) và đánh giá chu kỳ thăng trầm của nghề nuôi tu hài tại Việt Nam từ giai đoạn phát triển nóng (2005–2012) đến khủng hoảng dịch bệnh sưng vòi làm giảm 2/3 sản lượng (2012–2017) và quá trình phục hồi hiện nay.

Phần trọng tâm của chương phân tích bản chất di truyền của tính trạng tăng trưởng ở nhuyễn thể biển, khẳng định tính trạng này chịu chi phối bởi hệ số di truyền mức độ trung bình đến cao và chịu ảnh hưởng trực tiếp từ cấu trúc đa hình kiểu gene (như các minh chứng cụ thể trên các gene IGFBP, MSTN, E2F3, TPH, Amylase ở sò điệp, hàu). Tác giả so sánh nguyên lý, ưu thế và hạn chế của các thế hệ chỉ thị phân tử: từ chỉ thị lai RFLP; chỉ thị dựa trên PCR ngẫu nhiên như RAPD, AFLP; chỉ thị vi vệ tinh SSR; cho đến bước chuyển sang chỉ thị đa hình đơn nucleotide (SNP) thông lượng cao bằng công nghệ NGS (Roche 454, Illumina SBS, SOLiD, PacBio, Ion Torrent). Đồng thời, tác giả luận giải sâu về quy trình xây dựng mã vạch DNA theo chuẩn quốc tế BOLDSYSTEMS và kỹ thuật giải trình tự giảm độ phức tạp hệ gene ezRAD, từ đó xác lập tính cấp thiết của việc áp dụng NGS để giải mã hệ gene và phát hiện SNP trên đối tượng L. rhynchaena.

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

Chương 2 mô tả chi tiết quy trình chuẩn bị vật liệu, thiết kế thí nghiệm và phân tích dữ liệu:

  • Bố trí thu mẫu khảo sát thành phần loài và thực trạng nuôi tại 8 xã, thị trấn ven biển của huyện Vân Đồn.
  • Quy chuẩn chọn mẫu sinh học: 01 cá thể tu hài đại diện thu tại Vân Đồn dùng tách chiết DNA mô cơ màng áo phục vụ giải trình tự toàn bộ hệ gene (WGS) và tách RNA mô tiêu hóa phục vụ giải trình tự hệ gene phiên mã (Transcriptome).
  • Quy chuẩn chọn mẫu nghiên cứu tăng trưởng: Thu nhận 30 mẫu nhóm tăng trưởng nhanh (chọn ngẫu nhiên từ 63 gia đình có giá trị EBV cao nhất) và 30 mẫu nhóm tăng trưởng chậm (chọn ngẫu nhiên từ 43 gia đình có giá trị EBV thấp nhất) từ đàn tu hài tuyển chọn của Trường Cao đẳng Kinh tế, Kỹ thuật và Thủy sản.
  • Chi tiết quy trình phòng thí nghiệm: Phản ứng cắt hạn chế kép bằng MboI và Sau3AI, tinh sạch bằng hạt từ SPB loại bỏ phân đoạn <350 bp và >550 bp, cấu trúc adapter Illumina, nồng độ hóa chất trong phản ứng PCR khuếch đại thư viện ezRAD.
  • Quy trình phân tích tin sinh học: Thiết lập thông số xử lý dữ liệu đọc trình tự trên các phần mềm BWA, SAMtools, dDocent, Stacks, đánh giá chất lượng qua FastQC, phân tích cây phát sinh chủng loại bằng MEGA X (thuật toán ML) và BEAST (thuật toán BI), ước tính kích thước hệ gene qua phân tích k-mer và đánh giá độ hoàn chỉnh lắp ráp qua bộ dữ liệu chuẩn BUSCO.

Chương 3: Kết quả nghiên cứu và thảo luận

Chương 3 trình bày các kết quả thực nghiệm tương ứng với 4 nội dung nghiên cứu chính:

1. Thành phần loài, hiện trạng và tiềm năng phát triển nghề nuôi tu hài tại huyện Vân Đồn, tỉnh Quảng Ninh

  • Kết quả khảo sát thực địa khẳng định loài tu hài phân bố tự nhiên và nuôi thương phẩm chủ yếu tại Vân Đồn là loài tu hài vòi trắng Lutraria rhynchaena.
  • Điều kiện sinh thái huyện Vân Đồn rất thích hợp cho tu hài: diện tích mặt nước tự nhiên rộng lớn, địa hình đáy biển tương đối bằng phẳng gồm cát và sỏi sạn mảnh vỏ, biên độ nhật triều cao (đạt tới 4,8 m), dòng chảy ven bờ 0,15–0,4 m/s. Khu hệ thực vật phù du phong phú với 76 loài thuộc hai ngành tảo khuê (Bacillariophyta) và tảo giáp (Pyrrophyta), mật độ đạt từ $2,5 \times 10^6$ đến $8 \times 10^6$ tế bào/$m^3$, cung cấp nguồn thức ăn tự nhiên dồi dào.
  • Hiện trạng nghề nuôi: Tính đến năm 2019, diện tích nuôi trồng thủy sản toàn huyện đạt 4.258 ha (riêng nhuyễn thể đạt 3.100 ha, sản lượng khai thác nhuyễn thể đạt 11.000 tấn). Quy hoạch tổng thể mặt nước biển đạt 14.886 ha, trong đó quy hoạch nuôi nhuyễn thể là 3.927 ha (diện tích nuôi tu hài chiếm 32%).
  • Thách thức được ghi nhận: Tỷ lệ tự túc con giống tại chỗ chỉ đạt 5–10%, còn lại phụ thuộc nguồn giống trôi nổi hoặc nhập tiểu ngạch; hiện tượng ngọt hóa đột ngột vào mùa mưa bão (chiếm >75% lượng mưa cả năm) kết hợp triều kiệt gây chết hàng loạt ở các bãi nuôi đáy nông; rét đậm kéo dài vào mùa đông gây ảnh hưởng đến sinh trưởng.

2. Xây dựng mã vạch DNA cho tu hài Lutraria rhynchaena

  • Trình tự đoạn gene 16S rRNA (kích thước phân tích 437 bp) và COI (kích thước phân tích 615 bp) của các cá thể tu hài Vân Đồn được khuếch đại và giải trình tự thành công.
  • Đối chiếu Megablast với các trình tự của giống Lutraria trên GenBank (Lutraria lutraria, Lutraria rhynchaena, Lutraria sieboldii, Lutraria bruuni, Lutraria australis):
    • Trình tự gene 16S rRNA của mẫu tu hài Vân Đồn tương đồng cao nhất với các mẫu L. rhynchaena đã công bố trên GenBank, đồng thời xuất hiện các vị trí sai khác nucleotide đặc trưng riêng biệt khi so sánh với 4 loài còn lại trong giống.
    • Trình tự gene COI phân lập thể hiện các vị trí nucleotide sai khác mang tính đặc trưng loài, cho phép phân định rõ ràng giữa L. rhynchaena và các loài đồng chi.
  • Cây phát sinh chủng loại xây dựng bằng cả hai phương pháp Maximum Likelihood (MEGA X) và Bayesian Inference (BEAST) trên cả hai vùng gene 16S rRNA và COI đều chỉ ra mẫu tu hài Vân Đồn tập trung cùng nhánh đơn dòng (monophyletic) với các mẫu L. rhynchaena, tách biệt hoàn toàn có ý nghĩa thống kê với các loài khác trong giống Lutraria.
  • Thiết lập thành công chuỗi mã vạch DNA chuẩn và chuyển hóa thành mã QR chứa thông tin nhận diện sinh học phục vụ truy xuất nguồn gốc loài.

3. Giải trình tự hệ gene và hệ gene phiên mã của tu hài Lutraria rhynchaena

  • Dữ liệu WGS và Transcriptome được tạo ra bằng công nghệ giải trình tự thế hệ mới Illumina SBS.
  • Phân tích phân bố tần số k-mer cho phép ước lượng kích thước hệ gene, tỷ lệ dị hợp tử và nội dung các trình tự lặp lại trong hệ gene tu hài.
  • Dữ liệu giải trình tự transcriptome từ mô tiêu hóa được lắp ráp de novo, cung cấp bộ dữ liệu phiên mã phục vụ dự đoán và chú thích chức năng các gene mã hóa protein.
  • Kết quả đánh giá bằng BUSCO chứng minh tính hoàn chỉnh cao của bộ lắp ráp hệ gene khi so sánh với bộ dữ liệu tham chiếu chuẩn của động vật thân mềm. Dữ liệu hệ gene được so sánh với hệ gene của các loài hai mảnh vỏ đã công bố (Pinctada fucata, Chlamys farreri, Crassostrea gigas...), cung cấp hệ gene tham chiếu đại diện cho họ Mactridae trên ngân hàng dữ liệu GenBank.

4. Sàng lọc, tuyển chọn và phát triển chỉ thị SNP liên quan đến tính trạng tăng trưởng ở tu hài

  • Giải trình tự ezRAD-Seq trên 60 cá thể (30 cá thể tăng trưởng nhanh và 30 cá thể tăng trưởng chậm) tạo ra tập dữ liệu đọc đôi chất lượng cao sau khi kiểm định qua FastQC.
  • Quá trình dóng hàng trình tự lên hệ gene tham chiếu được thực hiện song song qua hai công cụ dDocent và Stacks, xác lập tỷ lệ định vị trình tự cao ở cả hai nhóm cá thể.
  • Phân tích phát hiện biến dị di truyền xác định hàng nghìn điểm đa hình SNP phân bố trên các contig/scaffolds của hệ gene.
  • Đánh giá sâu trên 12 gene chức năng được biết có liên quan đến tính trạng tăng trưởng ở động vật thân mềm, thống kê chính xác tần số xuất hiện của các SNP trong các vùng gene này.
  • Phân tích sàng lọc các locus ngoại lai (outlier loci) chịu áp lực chọn lọc phân kỳ giữa nhóm có giá trị EBV cao và nhóm có giá trị EBV thấp.
  • Luận án đã sàng lọc và xác định được 11 chỉ thị SNP tiềm năng có sự phân hóa tần số alen vượt trội giữa hai nhóm tăng trưởng nhanh và chậm.
  • Thiết kế thành công các cặp mồi PCR đặc hiệu cho các đoạn trình tự chứa 11 SNP tiềm năng này; kết quả chạy điện di kiểm chứng sản phẩm PCR cho thấy kích thước đoạn khuếch đại hoàn toàn trùng khớp với kích thước tính toán lý thuyết, xác nhận tính ứng dụng thực tế của bộ mồi trong phân tích kiểu gene.
Hạng mục dữ liệu / Phân tích Phương pháp / Công cụ Đối tượng / Cỡ mẫu Kết quả chính ghi nhận trong luận án
Khảo sát sinh thái & nghề nuôi Điều tra thực địa, phiếu hỏi, số liệu Phòng NN&PTNT Huyện Vân Đồn, Quảng Ninh 76 loài thực vật phù du ($2,5\text{--}8 \times 10^6$ tế bào/$m^3$); quy hoạch 3.927 ha nhuyễn thể (tu hài chiếm 32%); giống tại chỗ chỉ đáp ứng 5–10%.
Mã vạch DNA Khuếch đại PCR, Sanger, BLAST, MEGA X (ML), BEAST (BI) Đoạn gene 16S rRNA (437 bp) và COI (615 bp) Xác lập các sai khác nucleotide đặc trưng loài; cây phát sinh chủng loại khẳng định đơn dòng L. rhynchaena; lập mã QR.
Hệ gene & Transcriptome Illumina NGS, k-mer, BUSCO, dDocent, Stacks 50 mg cơ màng áo (WGS); 50 mg mô tiêu hóa (RNA-Seq) Lắp ráp de novo hệ gene tham chiếu; chú giải dữ liệu phiên mã; công bố hệ gene L. rhynchaena trên GenBank.
Chỉ thị SNP tăng trưởng ezRAD-Seq (enzyme MboI + Sau3AI, hạt SPB 350–550 bp) 30 cá thể tăng trưởng nhanh vs 30 cá thể tăng trưởng chậm (chọn theo EBV) Khảo sát SNP trong 12 gene chức năng; sàng lọc outlier loci; phát hiện 11 SNP tiềm năng; thiết kế cặp mồi PCR kiểm chứng thành công.

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

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

  • Cung cấp dẫn liệu khoa học cập nhật và có hệ thống về thành phần loài tu hài tự nhiên, đặc điểm sinh thái vùng phân bố, hiện trạng và tiềm năng phát triển nghề nuôi tại vùng biển đảo huyện Vân Đồn, tỉnh Quảng Ninh.
  • Xây dựng thành công cơ sở dữ liệu mã vạch DNA dựa trên hai chỉ thị gene ty thể 16S rRNA (437 bp) và COI (615 bp), cung cấp công cụ phân tử chuẩn xác để định danh và phân loại loài tu hài Lutraria rhynchaena, giải quyết triệt để sự nhầm lẫn hình thái học giữa các loài trong giống Lutraria.
  • Giải trình tự thành công toàn bộ hệ gene (WGS) và hệ gene phiên mã (Transcriptome) của loài tu hài Lutraria rhynchaena bằng kỹ thuật NGS, công bố dữ liệu hệ gene tham chiếu đầu tiên của loài này trên ngân hàng gene quốc tế GenBank, bổ sung nguồn dữ liệu cho nghiên cứu hệ gene học (genomics) ở lớp Bivalvia.
  • Xác định được biến dị di truyền ở mức độ toàn hệ gene giữa hai nhóm tu hài phân ly về tốc độ sinh trưởng; sàng lọc và phát hiện được 11 chỉ thị SNP liên kết với tính trạng tăng trưởng ở tu hài, đóng góp cơ sở dữ liệu cho di truyền học số lượng và sinh học phân tử thủy sản tại Việt Nam.

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

  • Mã vạch DNA và mã QR được thiết lập cung cấp giải pháp kỹ thuật chính xác cho các cơ quan kiểm định, quản lý giống và doanh nghiệp thủy sản nhằm truy xuất nguồn gốc giống, kiểm soát thương phẩm và bảo vệ bản quyền giống bản địa.
  • Bộ 11 chỉ thị SNP tiềm năng cùng các cặp mồi khuếch đại đặc hiệu cung cấp công cụ chẩn đoán kiểu gene phục vụ trực tiếp cho các chương trình chọn giống tu hài theo hướng tăng trưởng nhanh bằng phương pháp chọn lọc hỗ trợ bởi chỉ thị phân tử (MAS), giúp rút ngắn thời gian chọn giống và hạn chế hiện tượng thoái hóa giống do giao phối cận huyết.

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

  • Ứng dụng quy trình mã vạch DNA đã xây dựng vào công tác kiểm tra, chứng nhận nguồn gốc xuất xứ con giống tu hài tại các cơ sở sản xuất giống nhân tạo ở Quảng Ninh và các tỉnh ven biển.
  • Đưa bộ chỉ thị 11 SNP vào các chương trình chọn giống quy mô lớn của ngành thủy sản, kết hợp đánh giá giữa giá trị kiểu hình, giá trị EBV và phân tích kiểu gene để tuyển chọn đàn tu hài bố mẹ hậu bị có tốc độ tăng trưởng cao.
  • Chính quyền địa phương huyện Vân Đồn và tỉnh Quảng Ninh cần quy hoạch chi tiết vùng mặt nước nuôi biển, kiểm soát chất lượng con giống nhập khẩu tiểu ngạch, tăng cường đầu tư công nghệ sản xuất giống tại chỗ nhằm nâng cao tỷ lệ tự chủ con giống đạt chuẩn chất lượng.

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

  • Hạn chế:
    • Phạm vi thu mẫu phân tích thực nghiệm giải trình tự sâu và đánh giá chỉ thị SNP tập trung chủ yếu trên các nhóm gia đình tu hài nuôi giữ tại Trung tâm giống hải sản thuộc Trường Cao đẳng Kinh tế, Kỹ thuật và Thủy sản tại Vân Đồn, chưa mở rộng trên diện rộng đối với toàn bộ các quần đàn tự nhiên tại các vùng địa lý khác (như Cát Bà hay Khánh Hòa).
    • Số lượng 11 SNP tiềm năng bước đầu được xác thực qua thiết kế mồi và khuếch đại phân