Tổng quan luận án
Tính cấp thiết và khoảng trống nghiên cứu
Tu hài (Lutraria rhynchaena Jonas, 1844), còn gọi là tu hài vòi trắng (tên tiếng Anh: Snout Otter Clam hay Otter Clam), là loài nhuyễn thể hai mảnh vỏ có giá trị kinh tế và dinh dưỡng cao, phân bố tự nhiên và được nuôi tại các vùng biển nhiệt đới ấm (18–30°C, độ mặn 25–30‰) như Quảng Ninh, Hải Phòng và Khánh Hòa. Tại Việt Nam, diện tích nuôi tu hài đạt khoảng 1.000 ha với sản lượng khoảng 2.621,6 tấn/năm và giá trị kinh tế trên 200 tỷ đồng/năm. Thị trường nuôi thương phẩm đòi hỏi 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 nhiều trở ngại:
- Công tác sản xuất con giống phần lớn dựa vào phương pháp truyền thống, chưa áp dụng các quy trình chọn giống phân tử, dẫn đến tình trạng cận huyết, thoái hóa giống, sinh trưởng không đồng đều và suy giảm sức đề kháng trước dịch bệnh (điển hình là dịch bệnh sưng vòi bùng phát giai đoạn 2012–2013).
- Nguồn giống phụ thuộc vào đánh bắt tự nhiên hoặc nhập tiểu ngạch không rõ nguồn gốc, tỉ lệ sống thấp.
- Về mặt di truyền học, dữ liệu hệ gene tham chiếu hoàn chỉnh (reference genome) của loài Lutraria rhynchaena chưa được giải mã đầy đủ; các nghiên cứu chỉ thị phân tử trước đây mới dừng lại ở phân tích đa dạng di truyền quần thể bằng chỉ thị ty thể hoặc microsatellite, chưa phát triển được hệ thống chỉ thị đa hình nucleotide đơn (SNP) gắn với các tính trạng kinh tế cụ thể như tốc độ tăng trưởng để phục vụ chọn giống hiện đại.
Đề tài được thực hiện nhằm đáp ứng định hướng tại Quyết định số 1664/QĐ-TTg ngày 04/10/2021 của Thủ tướng Chính phủ về Đề á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, chú trọng khai thác nguồn giống bản địa qua các chương trình chọn giống khoa học.
Mục tiêu nghiên cứu
Luận án xác định 4 mục tiêu nghiên cứu cụ thể:
- Đá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.
- 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).
- 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).
- 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.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Loài tu hài vòi trắng (Lutraria rhynchaena, Jonas 1844), thuộc chi Lutraria, họ Vọp (Mactridae), bộ Ngao (Veneroida), lớp Hai mảnh vỏ (Bivalvia), ngành Thân mềm (Mollusca).
- Phạm vi không gian: Vùng biển huyện đảo Vân Đồn, tỉnh Quảng Ninh (thu mẫu thực địa tại 8 xã, thị trấn); các thí nghiệm sinh học phân tử được thực hiện tại Phòng thí nghiệm Công nghệ sinh học - Trường Cao đẳng Kinh tế, Kỹ thuật và Thủy sản (Bắc Ninh) và Trung tâm Công nghệ Genomics - Đại học Deakin (Geelong, Victoria, Australia).
- Phạm vi thời gian: Nghiên cứu được triển khai từ tháng 3/2018 đến tháng 12/2021.
- Phạm vi nội dung: Điều tra sinh thái quần thể tự nhiên và nghề nuôi; xây dựng mã vạch DNA từ vùng gene 16S rRNA và COI; giải trình tự toàn bộ hệ gene (WGS) và hệ gene phiên mã (transcriptome); sàng lọc biến dị di truyền SNP liên kết tính trạng tăng trưởng trên hai nhóm cá thể phân tách theo giá trị chọn giống ước tính (EBV).
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu trước đây
Tác giả đã hệ thống hóa các công trình nghiên cứu trong nước và quốc tế theo các nhóm nội dung:
-
Nghiên cứu về phân loại học và mã vạch DNA (DNA barcoding):
- Paul Hebert (2003) khởi xướng phương pháp mã vạch DNA sử dụng vùng gene COI ty thể (~650 bp) để định danh gần 200 loài bộ Cánh vảy với độ chính xác cao; tổ chức CBOL (2004) chuẩn hóa hệ thống dữ liệu mã vạch toàn cầu.
- Taberlet và cs. (2007) xác lập 4 tiêu chí cho mã vạch DNA chuẩn.
- Liu và cs. (2018) kết hợp gene ty thể (COI, 16S rRNA) và gene nhân (18S, 28S, Histone H3) để phân loại 56 loài ngao, áp dụng ngưỡng sai khác di truyền 3% đối với động vật thân mềm.
- Setiamarga và cs. (2019) phân tích phát sinh loài cho 5 loài sò tại vùng biển Nada (Nhật Bản) bằng gene COI và Histone.
- Ứng dụng mã vạch trong giám định và truy xuất nguồn gốc các loài hai mảnh vỏ thuộc họ Mytilidae, Pectinidae, Ostreidae được phát triển nhằm chống gian lận thương mại thực phẩm.
-
Nghiên cứu chọn giống theo kiểu gene và chỉ thị SNP liên quan đến tăng trưởng:
- Feng và cs. (2014) giải trình tự đoạn gene PyIGFBP dài 10.193 bp ở sò điệp, xác định SNP 1054 A > G liên quan đến kích thước vỏ và khối lượng thân (kiểu gene GG có giá trị tăng trưởng cao nhất, AG thấp nhất).
- Nghiên cứu gene Myostatin (CfMSTN) ở sò điệp Zhikong (Chlamys farreri) phát hiện SNP g-1162 G > T vùng 5' điều hòa phát triển cơ vân, trong đó kiểu gene TT có các chỉ tiêu sinh trưởng vượt trội so với GT và GG.
- Nghiên cứu các SNP trên gene PyE2F3-1 và PyE2F3-2 ở sò; phân tích haplotype vùng điều hòa gene CgTPH ở hàu Thái Bình Dương (Crassostrea gigas) xác định SNP -272 C/T tác động đến vị trí gắn yếu tố phiên mã FoxO và CUX, ảnh hưởng đến tổng hợp 5-HT và tốc độ lớn.
- Chỉ thị SNP trong gene amylase liên quan đến sinh trưởng ở Crassostrea gigas, và SNP trong gene Hsp70 liên quan đến tính kháng bệnh ở bào ngư Haliotis rubra.
-
Nghiên cứu giải trình tự hệ gene ở động vật hai mảnh vỏ:
- Đến thời điểm nghiên cứu, hệ gene của 13 loài thân mềm hai mảnh vỏ đã được công bố: Scapharca broughtonii, Pinctada fucata, Chlamys farreri, Argopecten purpuratus, Ruditapes philippinarum, Mytilus galloprovincialis, Saccostrea glomerata, Venustaconcha ellipsiformis, Modiolus philippinarum, Bathymodiolus platifrons, Limnoperna fortunei, Patinopecten yessoensis, và Crassostrea gigas.
-
Nghiên cứu di truyền phân tử trên đối tượng tu hài:
- Gan và cs. (2016) công bố hệ gene ty thể của Lutraria rhynchaena (kích thước ghi nhận chứa 12 gene mã hóa protein, 22 tRNA, tỉ lệ A+T chiếm 69%, G+C chiếm 31%).
- Thái Thanh Bình và cs. (2015, 2016) sử dụng đoạn gene 16S rRNA phân tích 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 (đa dạng haplotype 0,570 ± 0,09) và phát triển 12 chỉ thị microsatellite.
- Trần Thị Trang và cs. (2017) đánh giá đa dạng di truyền tu hài bằng microsatellite: số alen từ 6–12 (trung bình 9 alen/locus), $H_o = 0,371–0,459$, $H_e = 0,402–0,598$, hệ số cận huyết $F_{is} = 0,058–0,248$, sai khác quần thể $F_{st} = 0,015–0,121$.
- Năm 2020, hệ gene metagenome của Lutraria philippinarum tại Cát Bà được lắp ráp với kích thước khoảng 542 Mbp, gồm 2.405 scaffolds ($N_{50} = 0,9\text{ Mbp}$), dự đoán 30.533 gene.
Khoảng trống luận án giải quyết
Mặc dù đã có các nghiên cứu sơ khởi về chỉ thị microsatellite và hệ gene ty thể, loài tu hài Lutraria rhynchaena chưa có hệ gene tham chiếu chuẩn phục vụ cho các phân tích di truyền học số lượng và chức năng. Đồng thời, chưa có công trình nào sàng lọc hệ thống chỉ thị SNP quy mô toàn hệ gene liên kết trực tiếp với giá trị chọn giống tăng trưởng (EBV) ở tu hài. Luận án giải quyết khoảng trống này bằng cách cung cấp dữ liệu hệ gene nhân, hệ gene phiên mã và phát triển tập hợp các chỉ thị SNP tiềm năng phục vụ chọn giống phân tử.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết và khái niệm áp dụng
- Lý thuyết di truyền số lượng: Tính trạng tăng trưởng ở thân mềm là tính trạng số lượng chịu sự kiểm soát của đa gene và tương tác với môi trường, có hệ số di truyền từ trung bình đến cao (hệ số di truyền thường cho phép cải thiện 10–20% qua mỗi thế hệ chọn lọc). Cá thể dị hợp tử thường có ưu thế trao đổi chất và tốc độ lớn nhanh hơn cá thể đồng hợp tử.
- Lý thuyết phân loại phân tử và mã vạch DNA: Dựa trên mức độ biến dị và bảo tồn của các đoạn DNA ty thể (16S rRNA và COI) để phân định loài độc lập với các biến đổi hình thái do môi trường.
- Công nghệ giải trình tự thế hệ mới (NGS) và kỹ thuật ezRAD: Sử dụng enzyme cắt giới hạn kết hợp giải trình tự thông lượng cao để phân tích đa hình toàn hệ gene mà không cần biết trước toàn bộ thông tin trình tự.
Phương pháp nghiên cứu thực nghiệm
QUY TRÌNH NGHIÊN CỨU TỔNG THỂ
[Điều tra sinh thái & nghề nuôi] [Xây dựng mã vạch DNA] [Giải trình tự WGS & RNA-Seq]
- 40 điểm tại 8 xã/thị trấn - Vùng gene 16S rRNA (437 bp) - Mô cơ màng áo (DNA: WGS)
- Khung định lượng 10 m² - Vùng gene COI (615 bp) - Mô tiêu hóa (RNA: RNA-Seq)
- Mật độ (A) & Trữ lượng (B) - Cây phát sinh loài (ML, BI) - Đánh giá BUSCO, k-mers
- Mã QR tích hợp - Lắp ráp & Chú thích gene
[Sàng lọc SNP Tăng trưởng]
- 30 mẫu nhanh (EBV cao)
- 30 mẫu chậm (EBV thấp)
- ezRAD-Seq (MboI + Sau3AI)
- Pipeline dDocent & Stacks
- Xác định 11 SNP tiềm năng
- Thiết kế mồi PCR kiểm chứng
-
Phương pháp điều tra thực địa:
- Thiết lập 40 điểm khảo sát tại 8 xã/thị trấn huyện Vân Đồn theo phương pháp vùng triều (English S., 1994).
- Dùng khung định lượng $10\text{ m}^2$ để thu mẫu, xác định mật độ $A$ ($\text{con/m}^2$) và trữ lượng tức thời $W = B \times S$ (với $B$ là trữ lượng trung bình $\text{g/m}^2$, $S$ là diện tích bãi triều).
- Đánh giá mức độ suy giảm nguồn lợi $D = P_t / P_n$.
-
Phương pháp xây dựng mã vạch DNA:
- Thu thập và đối chiếu trình tự 16S rRNA và COI của chi Lutraria trên GenBank bằng Megablast.
- Căn chỉnh trình tự bằng Clustal Omega.
- Xây dựng cây quan hệ phát sinh chủng loại bằng phương pháp Maximum Likelihood (ML trên MEGA X) và Bayesian Inference (BI trên BEAST).
- Tạo lập mã QR tích hợp thông tin mã vạch DNA nhận dạng loài.
-
Phương pháp giải trình tự WGS và Transcriptome:
- Chiết xuất DNA từ 50 mg mô cơ màng áo và RNA từ 50 mg mô tiêu hóa của 01 cá thể tu hài.
- Chuẩn bị thư viện TruSeq Nano DNA và RNA; giải trình tự trên hệ thống Illumina tại Đại học Deakin (Australia).
- Ước lượng kích thước hệ gene qua phân tích tần số $k\text{-mers}$, kiểm tra độ hoàn chỉnh bằng BUSCO, lắp ráp de novo và chú giải chức năng gene.
-
Phương pháp ezRAD-Seq phát hiện SNP liên quan tăng trưởng:
- Vật liệu mẫu: 30 cá thể nhóm tăng trưởng nhanh (chọn ngẫu nhiên từ 63 gia đình có EBV cao nhất) và 30 cá thể nhóm tăng trưởng chậm (chọn từ 43 gia đình có EBV thấp nhất) từ đàn chọn giống của Trường Cao đẳng Kinh tế, Kỹ thuật và Thủy sản tại Vân Đồn.
- Quy trình thư viện: DNA được xử lý methyl hóa Dam, phân cắt bằng hai enzyme giới hạn MboI và Sau3AI, tinh sạch và chọn kích thước đoạn chèn 350–550 bp bằng hạt từ SPB (Sample Purification Beads), gắn adapter Illumina và khuếch đại thư viện bằng PCR.
- Xử lý tin sinh: Kiểm tra chất lượng dữ liệu bằng FastQC; dóng hàng lên hệ gene tham chiếu bằng BWA; gọi biến dị và lọc locus ngoại lai (outlier loci) bằng pipeline
dDocent và Stacks; xác định các SNP đặc thù và thiết kế mồi PCR kiểm chứng.
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 kiến thức nền tảng về sinh học, phân loại học, hiện trạng phân bố và khai thác tu hài trên thế giới và tại Việt Nam. Về mặt phân loại, tu hài Lutraria rhynchaena (tên đồng danh L. philippinarum) thuộc họ Mactridae. Về sinh thái, loài phân bố ở vùng triều đáy cát, sỏi, mùn bã hữu cơ, lọc ăn sinh vật phù du và có khả năng vùi sâu 5–7 cm dưới đáy biển.
Chương này phân tích các yếu tố tự nhiên tại huyện Vân Đồn: diện tích mặt nước rộng, nguồn thức ăn thực vật phù du phong phú gồm 76 loài (thuộc hai ngành Bacillariophyta và Pyrrophyta với mật độ $2,5 \times 10^6$ đến $8 \times 10^6\text{ tế bào/m}^3$), biên độ triều đạt tối đa 4,8 m, tạo môi trường thuận lợi cho nuôi nhuyễn thể. Tuy nhiên, các yếu tố rủi ro gồm hiện tượng ngọt hóa đột ngột do mưa bão mùa hè và rét đậm mùa đông kèm bão lớn ảnh hưởng đến tỉ lệ sống.
Phần tổng quan lý thuyết di truyền chỉ ra rằng chọn lọc tính trạng số lượng qua kiểu hình bị giới hạn bởi tương tác môi trường, do đó cần tiếp cận qua kiểu gene bằng các công nghệ chỉ thị DNA (từ RFLP, RAPD, AFLP, SSR đến SNP và kỹ thuật RAD-Seq/ezRAD) nhằm nâng cao độ chính xác chọn giống.
CHƯƠNG 2. VẬT LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU
Chương 2 mô tả chi tiết nguồn vật liệu sinh học, địa điểm thu mẫu và quy trình kỹ thuật phòng thí nghiệm và tin sinh học:
- Địa bàn và mẫu điều tra: 40 điểm khảo sát tại 8 xã/thị trấn ven biển Vân Đồn.
- Mẫu giải trình tự hệ gene: Sử dụng mô cơ màng áo và mô đường tiêu hóa của tu hài vòi trắng thu tại Vân Đồn.
- Mẫu phân tích SNP tăng trưởng: 60 cá thể phân chia thành 2 nhóm hình thái - di truyền đối lập dựa trên giá trị EBV (30 mẫu sinh trưởng nhanh từ 63 gia đình EBV cao; 30 mẫu sinh trưởng chậm từ 43 gia đình EBV thấp).
- Quy trình kỹ thuật: Mô tả chi tiết các phản ứng hóa sinh gồm tách chiết DNA/RNA, phản ứng Dam methyl hóa, cắt hạn chế MboI/Sau3AI, gắn adapter TruSeq, quy trình PCR khuếch đại thư viện ezRAD và các tham số vận hành phần mềm chuyên dụng (Megablast, Clustal Omega, MEGA X, BEAST, dDocent, Stacks, BWA, SAMtools, BUSCO).
| Thông số nghiên cứu |
Quy mô / Đặc tính kỹ thuật |
| Số điểm khảo sát thực địa |
40 điểm tại 8 xã/thị trấn thuộc huyện Vân Đồn |
| Mẫu giải mã hệ gene (WGS & Transcriptome) |
01 cá thể (50 mg mô cơ màng áo; 50 mg mô tiêu hóa) |
| Mẫu sàng lọc SNP liên quan tăng trưởng |
60 cá thể (30 cá thể nhóm EBV cao; 30 cá thể nhóm EBV thấp) |
| Enzyme cắt hạn chế sử dụng trong ezRAD |
MboI và Sau3AI |
| Kích thước đoạn DNA thư viện ezRAD chọn lọc |
350 – 550 bp |
| Công cụ dóng hàng và xử lý SNP |
dDocent, Stacks, BWA, SAMtools |
CHƯƠNG 3. KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN
3.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
Kết quả điều tra ghi nhận loài phân bố tự nhiên và nuôi thương phẩm chủ yếu tại Vân Đồn là tu hài vòi trắng Lutraria rhynchaena. Tổng diện tích quy hoạch phát triển nuôi trồng thủy sản của huyện đạt 14.886 ha, trong đó nuôi nhuyễn thể chiếm 3.927 ha (diện tích nuôi tu hài chiếm khoảng 32%). Các hình thức nuôi chính gồm: nuôi rải đáy bãi triều, nuôi lồng (rổ) treo trên bè và nuôi lồng đặt đáy.
Nghề nuôi tu hài mang lại hiệu quả kinh tế cao nhưng đối mặt với thách thức lớn về thiếu hụt con giống tại chỗ (nguồn sản xuất tại địa phương chỉ đáp ứng 5–10% nhu cầu, phần còn lại phải nhập từ các tỉnh khác hoặc tiểu ngạch không kiểm soát được chất lượng) và nguy cơ dịch bệnh bùng phát.
3.2. Xây dựng mã vạch DNA cho tu hài Lutraria rhynchaena
Nghiên cứu đã khảo sát và phân tích trình tự hai vùng gene ty thể đặc trưng của chi Lutraria:
- Vùng gene 16S rRNA (chiều dài phân tích 437 bp).
- Vùng gene COI (chiều dài phân tích 615 bp).
So sánh mức độ tương đồng và sai khác nucleotide giữa L. rhynchaena với các loài khác trong cùng chi (L. australis, L. philippinarum, v.v.) trên cơ sở dữ liệu GenBank cho thấy độ đặc hiệu phân loại cao. Cây phát sinh chủng loại xây dựng theo hai thuật toán độc lập (Maximum Likelihood - ML và Bayesian Inference - BI) đều phân tách rõ ràng loài L. rhynchaena thành một nhánh đơn ngành độc lập với độ tin cậy thống kê cao. Dữ liệu trình tự mã vạch đã được số hóa và tích hợp thành công vào cấu trúc mã QR phục vụ nhận dạng và truy xuất nguồn gốc.
| Vùng gene mã vạch |
Kích thước phân tích (bp) |
Phương pháp phân tích phát sinh loài |
Ứng dụng nhận dạng |
| 16S rRNA |
437 bp |
Maximum Likelihood (ML) & Bayesian (BI) |
Phân định loài, tích hợp mã QR |
| COI |
615 bp |
Maximum Likelihood (ML) & Bayesian (BI) |
Giám định phân tử, phân loại học |
3.3. Giải trình tự hệ gene và hệ gene phiên mã của Lutraria rhynchaena
- Hệ gene nhân (WGS): Dữ liệu giải trình tự NGS toàn bộ hệ gene đạt chất lượng cao; phân tích phân bố $k\text{-mers}$ cho phép ước lượng kích thước hệ gene và mức độ dị hợp tử. Dữ liệu lắp ráp de novo được đánh giá độ hoàn chỉnh thông qua bộ chuẩn BUSCO. Kết quả dự đoán và chú giải chức năng gene cung cấp bộ dữ liệu hệ gene tham chiếu chuẩn cho loài L. rhynchaena trên ngân hàng gene quốc tế GenBank.
- Hệ gene phiên mã (Transcriptome): Giải trình tự RNA-Seq từ mô tiêu hóa thu được tập hợp contig biểu hiện cao, hỗ trợ việc chú thích chính xác các vùng mã hóa protein và các con đường chuyển hóa liên quan đến dinh dưỡng và tăng trưởng.
3.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
- Dữ liệu ezRAD-Seq: Quá trình giải trình tự ezRAD-Seq trên 60 cá thể thuộc hai nhóm tăng trưởng nhanh ($N = 30$) và tăng trưởng chậm ($N = 30$) tạo ra bộ dữ liệu đọc chất lượng cao ở cả hai chiều đọc xuôi và ngược.
- Xử lý biến dị: Dữ liệu được dóng hàng lên hệ gene tham chiếu thông qua hai công cụ
dDocent và Stacks. Kết quả ghi nhận hàng nghìn biến dị phân tử phân bố trên toàn hệ gene.
- Sàng lọc Locus ngoại lai và SNP chức năng: Phân tích tần số alen và sàng lọc outlier loci giữa hai nhóm kiểu hình cực đoan đã phát hiện các biến dị có sự sai khác mang ý nghĩa thống kê. Trong đó, nghiên cứu tập trung phân tích biến dị trong 12 gene chức năng liên quan trực tiếp đến các con đường điều hòa sinh trưởng, phát triển cơ và chuyển hóa năng lượng.
- Phát hiện SNP tiềm năng: Luận án đã xác định được 11 chỉ thị SNP tiềm năng có mối liên kết chặt chẽ với tính trạng tăng trưởng nhanh ở tu hài.
- Thiết kế mồi kiểm chứng: Đã thiết kế thành công các cặp mồi PCR đặc hiệu bao quanh các vị trí SNP tiềm năng, tiến hành khuếch đại và kiểm chứng sản phẩm PCR, tạo công cụ phân tử sẵn sàng ứng dụng cho chọn lọc cá thể mang kiểu gene mong muốn.
Kết quả và những đóng góp mới
Đóng góp mới về mặt khoa học và lý luận
- Dữ liệu phân loại và phát sinh chủng loại: Cung cấp dẫn liệu khoa học chuẩn xác về mã vạch DNA của loài tu hài Lutraria rhynchaena dựa trên hai đoạn gene ty thể 16S rRNA (437 bp) và COI (615 bp), làm rõ mối quan hệ phát sinh loài trong giống Lutraria.
- Dữ liệu hệ gene tham chiếu: Lần đầu tiên giải mã và công bố dữ liệu toàn bộ hệ gene (WGS) và hệ gene phiên mã (Transcriptome) của tu hài vòi trắng Lutraria rhynchaena trên cơ sở dữ liệu quốc tế GenBank, bổ sung nguồn tư liệu genomic cho lớp Hai mảnh vỏ (Bivalvia).
- Phát hiện chỉ thị phân tử tăng trưởng: Ứng dụng kỹ thuật ezRAD-Seq xác định được 11 chỉ thị SNP tiềm năng liên kết với tính trạng tốc độ tăng trưởng ở tu hài, đóng góp cơ sở lý luận cho di truyền học chức năng ở động vật thân mềm.
Đóng góp mới về mặt thực tiễn
- Công cụ định danh và truy xuất nguồn gốc: Bộ mã vạch DNA và mã QR tương ứng cung cấp công cụ kỹ thuật chính xác phục vụ công tác giám định loài, quản lý nguồn lợi và truy xuất nguồn gốc thương mại tu hài.
- Công cụ chọn giống năng suất cao: Tập hợp 11 chỉ thị SNP cùng hệ thống mồi PCR đặc hiệu được thiết kế là công cụ phân tử trực tiếp hỗ trợ các chương trình chọn giống tu hài tăng trưởng nhanh (Marker-Assisted Selection - MAS), giúp rút ngắn thời gian chọn giống so với phương pháp chọn lọc kiểu hình truyền thống.
- Cơ sở quy hoạch nghề nuôi: Cung cấp số liệu hiện trạng về phân bố, trữ lượng và tiềm năng mặt nước tại huyện Vân Đồn (Quảng Ninh), làm căn cứ cho việc xây dựng kế hoạch phát triển vùng nuôi an toàn sinh học và chủ động nguồn giống.
Kiến nghị và giải pháp đề xuất
- Ứng dụng ngay các cặp mồi và 11 chỉ thị SNP đã sàng lọc vào thực tế sản xuất giống tu hài để đánh giá, tuyển chọn đàn bố mẹ có giá trị di truyền cao về tốc độ tăng trưởng.
- Mở rộng vùng kiểm chứng các chỉ thị SNP trên các quần thể tu hài nuôi và tự nhiên tại các địa phương khác như Hải Phòng, Khánh Hòa, Kiên Giang.
- Kết hợp dữ liệu hệ gene với các nghiên cứu miễn dịch học nhằm tìm kiếm thêm các chỉ thị phân tử kháng bệnh (đặc biệt là bệnh sưng vòi) để hoàn thiện gói chỉ thị đa tính trạng (tăng trưởng nhanh kết hợp kháng bệnh).
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của đề tài
- Phạm vi địa lý thu mẫu hệ gene: Mẫu phân tích giải trình tự WGS tập trung tại khu vực huyện Vân Đồn, tỉnh Quảng Ninh, chưa bao quát toàn bộ các vùng phân bố sinh học khác của tu hài tại miền Trung và miền Nam.
- Kỹ thuật đọc ngắn: Công nghệ giải trình tự NGS sử dụng đoạn đọc ngắn (short reads) của Illumina tuy có độ chính xác cao nhưng gặp thách thức nhất định trong việc lắp ráp các vùng trình tự lặp lại phức tạp cao của hệ gene thân mềm so với các công nghệ đọc dài.
- Quy mô quần thể kiểm chứng: Số lượng cá thể phân tích ezRAD-Seq giới hạn ở 60 mẫu thuộc hai nhóm phân ly cực đoan về EBV trong điều kiện nuôi thực nghiệm tại một trung tâm giống.
Hướng nghiên cứu tiếp theo
- Bổ sung công nghệ giải trình tự đoạn đọc dài (như Oxford Nanopore Technologies - ONT hoặc PacBio) kết hợp kỹ thuật bắt cặp nhiễm sắc thể thông lượng cao (Hi-C) để nâng cấp hệ gene Lutraria rhynchaena lên mức độ nhiễm sắc thể (chromosome-level assembly).
- Mở rộng sàng lọc trên quy mô đàn thương phẩm lớn ngoài tự nhiên để xác định chính xác tần số alen và mức độ ảnh hưởng của từng SNP đến khối lượng và chiều dài vỏ ở các giai đoạn phát triển khác nhau.
- Xây dựng chip sinh học (SNP chip panel) chuyên dụng phục vụ công tác chọn giống tu hài quy mô công nghiệp.
Giá trị tham khảo
Luận án là tài liệu chuyên khảo có giá trị ứng dụng và tham khảo cho nhiều nhóm đối tượng:
- Nghiên cứu sinh, học viên cao học và nhà nghiên cứu: Thuộc các chuyên ngành Di truyền học, Sinh học phân tử, Công nghệ sinh học và Nuôi trồng thủy sản; có thể sử dụng quy trình xây dựng thư viện ezRAD, quy trình xử lý dữ liệu dDocent/Stacks và dữ liệu hệ gene tham chiếu trên GenBank cho các nghiên cứu tương tự trên động vật thân mềm.
- Giảng viên đại học: Sử dụng làm học liệu thực tế trong giảng dạy các học phần Di truyền số lượng, Công nghệ gen trong thủy sản, và Tin sinh học ứng dụng.
- Cơ sở sản xuất giống và cơ quan quản lý thủy sản: Các trung tâm giống thủy sản ven biển (tại Quảng Ninh, Hải Phòng, Khánh Hòa) có thể ứng dụng trực tiếp các chỉ thị SNP và quy trình PCR mồi đặc hiệu để kiểm định chất lượng đàn giống bố mẹ; các nhà hoạch định chính sách có thêm căn cứ dữ liệu để quy hoạch vùng nuôi biển bền vững theo Quyết định 1664/QĐ-TTg.
Câu hỏi thường gặp
1. Đối tượng nghiên cứu chính và địa bàn thực địa của luận án là gì?
Đối tượng nghiên cứu là loài tu hài vòi trắng (Lutraria rhynchaena Jonas, 1844), thuộc họ Vọp (Mactridae), bộ Ngao (Veneroida), lớp Hai mảnh vỏ (Bivalvia). Địa bàn điều tra thực địa, thu mẫu sinh thái và thử nghiệm nuôi được thực hiện tại 40 điểm thuộc 8 xã, thị trấn của huyện đảo Vân Đồn, tỉnh Quảng Ninh.
2. Hai vùng gene nào được sử dụng để xây dựng mã vạch DNA nhận dạng loài tu hài trong luận án?
Luận án sử dụng hai vùng gene ty thể đặc trưng:
- Đoạn gene 16S rRNA có kích thước phân tích là 437 bp.
- Đoạn gene COI (Cytochrome C Oxidase subunit I) có kích thước phân tích là 615 bp.
Cả hai vùng trình tự này đều cho phép phân tách chính xác loài Lutraria rhynchaena khỏi các loài họ hàng trong chi Lutraria trên cây phát sinh chủng loại xây dựng bằng thuật toán ML và BI.
3. Kỹ thuật di truyền phân tử nào được dùng để sàng lọc biến dị SNP giữa hai nhóm tu hài tăng trưởng nhanh và chậm?
Nghiên cứu sử dụng kỹ thuật giải trình tự DNA liên kết vị trí enzyme cắt hạn chế (ezRAD-Seq). DNA tổng số của 60 cá thể (30 cá thể nhóm tăng trưởng nhanh và 30 cá thể nhóm tăng trưởng chậm, phân loại dựa trên giá trị chọn giống ước tính EBV) được cắt đồng thời bằng hai enzyme MboI và Sau3AI, gắn adapter TruSeq và giải trình tự trên hệ thống máy Illumina.
4. Luận án đã phát hiện được bao nhiêu chỉ thị SNP tiềm năng liên quan đến tính trạng tăng trưởng ở tu hài?
Qua quy trình xử lý tin sinh học bằng pipeline dDocent và Stacks, đối chiếu với hệ gene tham chiếu và phân tích locus ngoại lai (outlier loci) trong 12 gene chức năng, luận án đã xác định và sàng lọc được 11 chỉ thị SNP tiềm năng có liên kết chặt chẽ với tính trạng tăng trưởng nhanh, đồng thời thiết kế thành công các cặp mồi PCR đặc hiệu để khuếch đại các phân đoạn chứa các SNP này.
5. Dữ liệu hệ gene của tu hài được giải mã bằng phương pháp nào và có ý nghĩa gì?
Toàn bộ hệ gene (WGS từ mô cơ màng áo) và hệ gene phiên mã (Transcriptome từ mô tiêu hóa) của Lutraria rhynchaena được giải mã bằng công nghệ giải trình tự thế hệ mới (NGS) trên nền tảng Illumina tại Đại học Deakin (Australia). Dữ liệu này lần đầu tiên cung cấp hệ gene tham chiếu (reference genome) hoàn chỉnh cho loài trên GenBank, phục vụ việc lắp ráp, chú thích gene và làm nền tảng dóng hàng dữ liệu ezRAD-Seq để phát hiện biến dị SNP.
Kết luận
Luận án của nghiên cứu sinh Triệu Anh Tuấn đã giải quyết trọn vẹn các mục tiêu nghiên cứu từ điều tra sinh thái thực địa đến ứng dụng công nghệ hệ gene học hiện đại trên đối tượng tu hài (Lutraria rhynchaena Jonas, 1844) tại Vân Đồn, Quảng Ninh. Công trình đã xây dựng thành công bộ mã vạch DNA chuẩn tích hợp mã QR nhận dạng loài từ hai đoạn gene 16S rRNA và COI, đồng thời lần đầu tiên công bố dữ liệu toàn bộ hệ gene và hệ gene phiên mã của loài trên cơ sở dữ liệu quốc tế GenBank. Quan trọng nhất, việc phát hiện và thiết kế mồi kiểm chứng cho 11 chỉ thị SNP tiềm năng liên kết với tính trạng tăng trưởng thông qua kỹ thuật ezRAD-Seq đã tạo ra công cụ phân tử phục vụ trực tiếp cho công tác chọn tạo giống tu hài năng suất cao, đóng góp thiết thực cho sự phát triển bền vững của ngành nuôi biển Việt Nam.