Bối cảnh và vấn đề nghiên cứu

Ngành chăn nuôi gia cầm, đặc biệt là thủy cầm, đóng vai trò kinh tế quan trọng trong nền nông nghiệp Việt Nam, tạo sinh kế và cải thiện thu nhập cho các nông hộ. Việt Nam hiện là quốc gia đứng thứ hai thế giới về sản lượng nuôi thủy cầm, tập trung hơn 60% đàn vịt tại hai vùng đồng bằng sông Hồng và sông Cửu Long. Các giống vịt bản địa của Việt Nam sở hữu nguồn gen quý với khả năng thích nghi cao, chống chịu dịch bệnh tốt và chất lượng thịt, trứng đặc thù. Tuy nhiên, trước áp lực du nhập các giống ngoại cao sản và tập quán chăn thả tự do không kiểm soát, các giống vịt bản địa đang đối mặt với nguy cơ lai tạp, suy giảm độ thuần chủng và xói mòn nguồn gen.

Đề tài khóa luận tốt nghiệp "Polymorphism of partial D-loop mitochondrial DNA region and phylogenetics in five duck breeds in Vietnam" (Đa hình đoạn vùng D-loop DNA ty thể và phân tích phát sinh chủng loại ở năm giống vịt tại Việt Nam) do sinh viên Bùi Bình An thực hiện dưới sự hướng dẫn của ThS. Nguyễn Quốc Trung tại Học viện Nông nghiệp Việt Nam nhằm đánh giá tính đa dạng di truyền và xác định nguồn gốc dòng mẹ của các giống vịt bản địa.

Nhiệm vụ nghiên cứu được xác định cụ thể gồm 4 nội dung:

  1. Thu thập mẫu máu và mẫu chân lông vũ (bloodfeathers) từ 5 giống vịt, tiến hành tách chiết DNA tổng số.
  2. Khuếch đại đoạn gen vùng D-loop ty thể bằng kỹ thuật PCR và giải trình tự nucleotide.
  3. Xác định các điểm đa hình nucleotide đơn (SNPs) và nhận diện các haplotype trong các chuỗi trình tự thu được.
  4. Phân tích mối quan hệ phát sinh chủng loại (phylogenetics) và phân bố haplotype của 5 giống vịt nghiên cứu.

Nghiên cứu được triển khai trên 5 giống vịt tại Việt Nam: vịt Cỏ, vịt Cổ Lũng, vịt Nà Tấu, vịt Trời và vịt Đốm (vịt Pất Lài). Thời gian thực hiện từ tháng 08/2022 đến tháng 01/2023 tại Bộ môn Sinh học Phân tử và Công nghệ Sinh học Ứng dụng, Khoa Công nghệ Sinh học, Học viện Nông nghiệp Việt Nam.

Cơ sở lý thuyết và phương pháp

Cơ sở lý thuyết và chỉ thị phân tử

Hệ gen ty thể (mtDNA) ở gia cầm có cấu trúc mạch vòng khép kín kích thước khoảng 16.785 bp. Vùng kiểm soát (D-loop hay displacement loop) là vùng không mã hóa protein, có tốc độ đột biến và tiến hóa nhanh hơn các vùng khác trong hệ gen ty thể. Do đặc tính di truyền theo dòng mẹ và không xảy ra hiện tượng tái tổ hợp di truyền như ADN nhân, trình tự vùng D-loop là chỉ thị phân tử chuẩn xác để truy nguyên nguồn gốc thuần hóa, xác định cấu trúc di truyền quần thể và dựng cây phát sinh chủng loại. Các giả thuyết về nguồn gốc vịt nhà tại Đông Á tập trung vào hai loài tổ tiên chính là vịt cổ xanh (Anas platyrhynchos) và vịt cánh trắng (Anas zonorhyncha), cùng với sự du nhập của vịt bơ / vịt xiêm (Cairina moschata).

Phương pháp nghiên cứu thực nghiệm

Nghiên cứu sử dụng tổng số 125 mẫu sinh học thuộc 5 giống vịt do Khoa Chăn nuôi, Học viện Nông nghiệp Việt Nam cung cấp. Quy trình thực nghiệm và phân tích tin sinh học bao gồm:

  • Tách chiết DNA tổng số:
    • Đối với mẫu máu (66 mẫu lấy từ tĩnh mạch cánh, bảo quản bằng ống chống đông EDTA): Tách chiết theo quy trình của Kulkarni và cộng sự (2015) sử dụng đệm lysis, SDS 10%, NaCl 5M và tủa bằng isopropanol (không dùng phenol hay enzyme phân giải protein).
    • Đối với mẫu chân lông vũ (59 mẫu): Tách chiết theo phương pháp cải tiến của Natalia và cộng sự (2001) với bước ủ Proteinase K (175 mg/ml) tại 60°C trong 4 giờ và tinh sạch bằng Phenol : Chloroform : Isoamylalcohol (25:24:1).
  • Điện di kiểm tra DNA: Điện di trên gel agarose 1% dưới ánh sáng tử ngoại (UV), định lượng sơ bộ bằng thang chuẩn Lamda DNA (100 ng/µL).
  • Phản ứng PCR: Khuếch đại đoạn vùng D-loop bằng cặp mồi L78 (5'-GTT ATT TGG TTA TGC ATA TCG TG-3') và H774 (5'-CCA TAT ACG CCA ACC GTC TC-3') (Sorenson & Fleischer, 1996). Thể tích phản ứng 40 µL sử dụng MyTaq Mix 2X trên máy luân nhiệt Agilent SureCycler 8800.
  • Giải trình tự: Sản phẩm PCR đạt chuẩn được gửi tinh sạch và giải trình tự theo phương pháp Sanger tại 1st BASE (Singapore).
  • Xử lý và phân tích dữ liệu:
    • Kiểm tra chất lượng sắc ký đồ bằng Chromas 2.4; căn chỉnh và cắt bỏ vùng nhiễu đầu/cuối bằng BioEdit 5.9 và ClustalX 2.1.
    • Tính toán các chỉ số đa dạng di truyền (số vị trí biến dị $S$, số haplotype $NHap$, độ đa dạng haplotype $Hd$, độ đa dạng nucleotide $Pi$, kiểm định trung tính Tajima's D) bằng phần mềm DnaSP 6.
    • Xây dựng cây phát sinh chủng loại bằng phương pháp Maximum Likelihood (ML) trên phần mềm MEGA11 với phương pháp kiểm tra Bootstrap 100 lần lặp.
    • Trình tự đối chứng từ GenBank (NCBI): Anas platyrhynchos (JN811041.1), Anas zonorhyncha (GU246018.1), Anas americana (HM063480.1) và Cairina moschata (GQ922096.1).

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

Chapter I: Introduction (Mở đầu)

Chương I nêu rõ vị trí của ngành chăn nuôi vịt tại Việt Nam và trên thế giới (châu Á chiếm 82,6% sản lượng thịt vịt toàn cầu). Tác giả làm rõ tính cấp thiết của việc bảo tồn nguồn gen 5 giống vịt bản địa (Cỏ, Cổ Lũng, Nà Tấu, Trời, Đốm) và xác định mục tiêu sử dụng vùng D-loop mtDNA làm công cụ giải mã đa hình di truyền và nguồn gốc tiến hóa.

Chapter II: Literature Review (Tổng quan tài liệu)

Chương này tổng hợp các công trình nghiên cứu trong nước và quốc tế về:

  • Hiện trạng phân bố, sản lượng vịt trên thế giới và vị thế của Việt Nam (tổng đàn 82,5 triệu con năm 2020).
  • Tổng quan đa dạng sinh học các giống vịt Việt Nam và các chính sách bảo tồn nguồn gen tại địa phương (vịt Bầu Bến, Bầu Quỳ, Cổ Lũng, Nà Tấu, Đốm...).
  • Đặc điểm sinh học và ngoại hình của 5 giống nghiên cứu: vịt Cỏ (thích nghi chăn thả chạy đồng, lông không đồng nhất), vịt Cổ Lũng (đặc sản Bá Thước, Thanh Hóa), vịt Nà Tấu (vịt bầu cổ ngắn Điện Biên), vịt Trời (Anas platyrhynchos thuần hóa), vịt Đốm (vịt Pất Lài Lạng Sơn).
  • Cơ sở lý thuyết về hệ gen ty thể và hiệu quả của vùng D-loop trong phân tích phả hệ dòng mẹ.

Chapter III: Materials and Methods (Vật liệu và phương pháp nghiên cứu)

Mô tả chi tiết 125 mẫu nghiên cứu, bao gồm mã mẫu, giống, giới tính và loại mẫu (máu hoặc chân lông vũ). Cung cấp đầy đủ công thức hóa chất, nồng độ thành phần phản ứng PCR và chu trình nhiệt (biến tính 95°C trong 5 phút; 35 chu kỳ: 95°C - 30 giây, 60°C - 30 giây, 72°C - 1 phút; kéo dài cuối 72°C - 5 phút). Trình bày chi tiết các thông số thiết lập trên các phần mềm DnaSP 6 và MEGA11.

Chapter IV: Results and Discussion (Kết quả và thảo luận)

Chương IV trình bày các kết quả thực nghiệm và thảo luận chuyên sâu:

  • Kết quả tách chiết DNA và PCR: DNA tách từ chân lông vũ theo phương pháp Natalia cho vạch điện di sắc nét, độ tinh sạch cao và ít lẫn tạp hơn phương pháp tách từ máu. Sản phẩm PCR đạt kích thước dự kiến nằm trong khoảng 600 - 700 bp. Mẫu VD09F không lên vạch PCR nên bị loại.
  • Sàng lọc trình tự: Trong số các mẫu giải trình tự, 35 mẫu bị loại do tín hiệu nhiễu hoặc độ dài dưới 200 bp; 6 mẫu bị loại do sai khác quá lớn khi căn gióng; 84 mẫu đạt chuẩn (đoạn 658 bp vùng 5') được đưa vào phân tích tin sinh học.
  • Đa dạng di truyền quần thể: Phát hiện 41 vị trí biến dị (SNPs), chia 84 cá thể thành 17 haplotype. Vịt Cỏ có độ đa dạng haplotype cao nhất ($Hd = 0,9130$), trong khi vịt Cổ Lũng có độ đa dạng haplotype thấp nhất ($Hd = 0,5263$) nhưng độ đa dạng nucleotide ($Pi$) cao nhất.
  • Phân tích phả hệ từng giống:
    • Vịt Cỏ (24 mẫu): Phân thành 15 haplotype đơn giống / 7 haplotype chung, thể hiện nguồn gốc lai tạp di truyền cao; có quan hệ gần với Anas platyrhynchosCairina moschata.
    • Vịt Cổ Lũng (19 mẫu): Chia thành 2 nhóm rõ rệt, nhóm 6 cá thể mái có quan hệ gần với Anas zonorhyncha, nhóm 13 cá thể trống có quan hệ gần với Cairina moschata.
    • Vịt Nà Tấu (12 mẫu) và vịt Đốm (16 mẫu): Cho thấy tính đồng nhất di truyền cao trong nội bộ giống khi phân tích toàn bộ quần thể, có chung nguồn gốc dòng mẹ với Anas zonorhynchaCairina moschata.
    • Vịt Trời (13 mẫu): Có 4 vị trí biến dị, liên hệ gần với Cairina moschataAnas americana.

Kết quả và đóng góp

Số liệu thực nghiệm chính

Quá trình phân tích 84 trình tự vùng D-loop ty thể chiều dài 658 bp ở 5 giống vịt bản địa ghi nhận các thông số di truyền cụ thể:

Giống vịt Số lượng mẫu ($N$) Số vị trí biến dị ($S$) Số haplotype ($NHap$) Độ đa dạng haplotype ($Hd$) Độ đa dạng nucleotide ($Pi$)
Vịt Cỏ 24 9 15 (đơn giống) / 7 (chung) 0,9130 Trung bình
Vịt Cổ Lũng 19 17 3 0,5263 Cao nhất
Vịt Nà Tấu 12 5 5 (đơn giống) / 1 (chung) 0,8017 Trung bình
Vịt Trời 13 4 4 (đơn giống) / 5 (chung) 0,8017 Thấp nhất
Vịt Đốm 16 9 3 (đơn giống) / 1 (chung) 0,8017 Trung bình
Toàn bộ quần thể 84 41 17 0,8975 0,0085

Đóng góp học thuật và thực tiễn

  • Minh chứng nguồn gốc tiến hóa dòng mẹ: Kết quả cây phát sinh chủng loại khẳng định 5 giống vịt bản địa Việt Nam có quan hệ di truyền nguồn gốc dòng mẹ gần gũi với hai loài chính là Anas zonorhyncha (vịt cánh trắng) và Cairina moschata (vịt bơ/vịt xiêm), đồng thời có sự phân nhóm sâu sắc giữa các cá thể đực và cái ở giống vịt Cổ Lũng.
  • So sánh với các nghiên cứu tiền nhiệm: Nghiên cứu chỉ ra chỉ số đa dạng haplotype ($Hd = 0,8975$) cao hơn đáng kể so với các công bố trước đây của N.V. Ba (2021) trên vịt Cỏ ($Hd = 0,00113$) và vịt Cổ Lũng ($Hd = 0,00371$), điều này được tác giả lý giải do dung lượng mẫu phân tích trên từng giống trong nghiên cứu này lớn hơn và phản ánh thực tế tập quán chăn thả tự do kéo dài.
  • Tối ưu hóa phương pháp tách mẫu: Khẳng định việc sử dụng mẫu chân lông vũ kết hợp ly giải bằng Proteinase K và tinh sạch Phenol:Chloroform cho hiệu quả thu nhận DNA chất lượng cao hơn mẫu máu bảo quản lạnh, giúp giảm thiểu xâm lấn trên cá thể gia cầm sống.

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

  • Hạn chế về chất lượng mẫu ban đầu: Tỷ lệ mẫu giải trình tự bị loại bỏ còn tương đối cao (35 mẫu bị nhiễu hoặc ngắn dưới 200 bp, 1 mẫu hỏng PCR) do điều kiện bảo quản mẫu máu dễ bị đông và suy thoái chất lượng trong quá trình lưu trữ trước khi xử lý.
  • Phạm vi chỉ thị: Đề tài mới chỉ khảo sát một đoạn ngắn 658 bp trên vùng D-loop của hệ gen ty thể, phản ánh quan hệ phả hệ theo dòng mẹ đơn độc mà chưa phản ánh được toàn diện biến dị di truyền từ phía dòng bố.
  • Hướng nghiên cứu tiếp: Cần mở rộng giải trình tự toàn bộ hệ gen ty thể (mitogenome) kết hợp với các chỉ thị gen nhân (microsatellite hoặc giải trình tự SNP toàn hệ gen) trên quy mô mẫu rộng hơn tại các vùng sinh thái khác nhau để đánh giá toàn diện mức độ lai tạp và phục vụ công tác chọn giống chuẩn xác.

Giá trị tham khảo

  • Đối tượng hưởng lợi: Khóa luận là tài liệu tham khảo trực tiếp cho sinh viên, học viên cao học và nghiên cứu viên các ngành Công nghệ Sinh học, Di truyền học Động vật, Chăn nuôi - Thú y và Bảo tồn Đa dạng Sinh học.
  • Nội dung giá trị chuyên môn cao:
    • Quy trình chuẩn tách chiết DNA từ chân lông vũ gia cầm không gây tổn thương cá thể.
    • Bộ thông số chu trình nhiệt PCR cho cặp mồi L78/H774 khuếch đại vùng D-loop gia cầm.
    • Phương pháp tích hợp các công cụ tin sinh học (BioEdit, DnaSP 6, MEGA11) trong phân tích đa hình SNPs, phân định haplotype và dựng cây phát sinh loài Maximum Likelihood.

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

1. Đề tài khảo sát những giống vịt nào tại Việt Nam và cỡ mẫu phân tích cuối cùng là bao nhiêu?

Nghiên cứu thu thập 125 mẫu từ 5 giống vịt: Cỏ, Cổ Lũng, Nà Tấu, Trời và Đốm (Pất Lài). Sau các bước kiểm tra chất lượng PCR và loại bỏ các trình tự ngắn/nhiễu, tổng cộng 84 trình tự D-loop đạt chuẩn (độ dài 658 bp) được sử dụng để phân tích di truyền (gồm 24 vịt Cỏ, 19 vịt Cổ Lũng, 12 vịt Nà Tấu, 13 vịt Trời và 16 vịt Đốm).

2. Cặp mồi nào được sử dụng để khuếch đại vùng D-loop ty thể trong nghiên cứu?

Nghiên cứu sử dụng cặp mồi L78 (trình tự xuôi: 5'-GTT ATT TGG TTA TGC ATA TCG TG-3') và H774 (trình tự ngược: 5'-CCA TAT ACG CCA ACC GTC TC-3') theo Sorenson & Fleischer (1996), cho sản phẩm PCR kích thước dự kiến từ 600 đến 700 bp.

3. Phương pháp tách chiết DNA nào cho kết quả tốt hơn giữa mẫu máu và mẫu chân lông vũ?

Phương pháp tách chiết từ chân lông vũ (sử dụng Proteinase K và Phenol:Chloroform:Isoamylalcohol theo Natalia et al., 2001) cho kết quả vượt trội hơn so với tách từ máu (theo Kulkarni et al., 2015). DNA từ lông vũ có độ tinh sạch cao, vạch điện di sắc nét, ít vết nhòe và cho kết quả khuếch đại PCR đậm rõ hơn.

4. Giống vịt nào có độ đa dạng di truyền cao nhất và thấp nhất trong nghiên cứu?

Vịt Cỏ có độ đa dạng haplotype ($Hd$) cao nhất đạt 0,9130 (91,3%) với 15 haplotype đơn giống được ghi nhận. Ngược lại, vịt Cổ Lũng có độ đa dạng haplotype thấp nhất ($Hd = 0,5263$) với 3 haplotype, nhưng lại ghi nhận số vị trí biến dị ($S = 17$) và độ đa dạng nucleotide ($Pi$) cao nhất trong 5 giống.

5. Cây phát sinh chủng loại cho thấy 5 giống vịt Việt Nam có nguồn gốc dòng mẹ gần với các loài nào trên NCBI?

Phân tích cây phát sinh chủng loại Maximum Likelihood chỉ ra 5 giống vịt Việt Nam có quan hệ phát sinh dòng mẹ gần gũi với loài Anas zonorhyncha (vịt cánh trắng, mã tiếp nhận GU246018.1) và Cairina moschata (vịt xiêm/vịt bơ, mã tiếp nhận GQ922096.1), riêng một phân nhóm vịt Trời có liên hệ gần với Anas americana (HM063480.1).

Kết luận

Khóa luận tốt nghiệp của tác giả Bùi Bình An đã phân tích thành công đoạn trình tự 658 bp vùng D-loop DNA ty thể trên 84 cá thể thuộc 5 giống vịt bản địa Việt Nam, phát hiện 41 vị trí đa hình nucleotide đơn (SNPs) phân bổ trên 17 haplotype với chỉ số đa dạng haplotype quần thể cao ($Hd = 0,8975$). Kết quả nghiên cứu xác định mối quan hệ di truyền dòng mẹ chặt chẽ giữa các giống vịt bản địa với Anas zonorhynchaCairina moschata, cung cấp dẫn liệu phân tử xác thực hỗ trợ công tác phân loại, chọn giống và xây dựng chiến lược bảo tồn nguồn gen thủy cầm tại Việt Nam.