Tổng quan luận án

Công trình nghiên cứu "Xác định và phân tích hoàn chỉnh trình tự hệ gen ty thể của 6 giống lợn bản địa tại một số tỉnh Miền Bắc Việt Nam" của tác giả Bùi Anh Tuấn (chuyên ngành Công nghệ sinh học, mã số 94 20 20 1) được thực hiện tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam dưới sự hướng dẫn khoa học của GS. Nghiêm Ngọc Minh và PGS. Võ Thị Bích Thủy.

Tính cấp thiết và khoảng trống nghiên cứu

Ngành chăn nuôi lợn tại Việt Nam đang đối mặt với sự suy giảm số lượng nghiêm trọng của các giống lợn bản địa do xu hướng nhập nội các giống ngoại năng suất cao và áp lực khai thác thương mại. Nhiều giống lợn nội quý rơi vào tình trạng suy thoái đồng huyết hoặc đứng trước nguy cơ tuyệt chủng (như giống lợn Ỉ với số lượng cá thể sinh sản suy giảm từ 2 triệu con năm 1969 xuống còn dưới 100 con tại Thanh Hóa và Viện Chăn nuôi). Dù Việt Nam sở hữu khoảng 26 giống lợn bản địa có nhiều đặc tính thích nghi, kháng bệnh và chất lượng thịt đặc trưng, phần lớn các nghiên cứu trước đây chỉ tập trung vào chỉ tiêu sinh trưởng, sinh sản hoặc khảo sát một vài đoạn gen đơn lẻ. Đến thời điểm tác giả thực hiện nghiên cứu, chưa có công trình nào giải mã và phân tích toàn diện, hoàn chỉnh hệ gen ty thể của các giống lợn bản địa phía Bắc Việt Nam để làm sáng tỏ nguồn gốc, quan hệ phát sinh chủng loại và xây dựng cơ sở dữ liệu phân tử phục vụ bảo tồn nguồn gen.

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

Luận án xác định ba mục tiêu cụ thể:

  1. Thu được dữ liệu hoàn chỉnh hệ gen ty thể của sáu giống lợn bản địa Việt Nam (lợn Ỉ, lợn Móng Cái, lợn Mường Lay, lợn Hương, lợn Mường Khương và lợn Hạ Lang), đăng ký trên Ngân hàng gen quốc tế (GenBank - NCBI).
  2. Xác định thành phần, cấu trúc hệ gen ty thể, so sánh sự sai khác trình tự, xác định đặc điểm di truyền đặc trưng của sáu giống lợn bản địa, qua đó đóng góp vào cơ sở dữ liệu phục vụ công tác nhận dạng và bảo tồn.
  3. Xác định mối quan hệ về di truyền, nhận định nguồn gốc và phát sinh chủng loại của sáu giống lợn bản địa Việt Nam trong tương quan với các giống lợn khu vực và trên thế giới.

Nội dung nghiên cứu

Luận án triển khai 6 nội dung kỹ thuật:

  1. Điều tra, khảo sát về giống, nơi cư trú, thu thập mẫu máu của 6 giống lợn bản địa.
  2. Giải trình tự toàn bộ hệ gen ty thể của 6 giống lợn bản địa.
  3. Lắp ráp, xác định trình tự hoàn chỉnh toàn bộ hệ gen ty thể và chú giải cấu trúc gen.
  4. Phân tích thành phần nucleotide và cấu trúc hệ gen.
  5. Nghiên cứu đa hình trình tự, so sánh trình tự hệ gen ty thể của 6 giống lợn này với một số giống lợn ở Châu Á và Châu Âu.
  6. Xây dựng cây phát sinh chủng loại dựa trên trình tự vùng D-loop và trình tự hoàn chỉnh của hệ gen ty thể, phân tích mối quan hệ nguồn gốc phát sinh chủng loại.

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

  • Đối tượng: Các cá thể thuộc 6 giống lợn bản địa Việt Nam gồm: lợn Ỉ, lợn Móng Cái, lợn Mường Khương, lợn Mường Lay, lợn Hương và lợn Hạ Lang. Riêng đối với quần thể lợn Ỉ tại Thanh Hóa, do không còn tìm thấy cá thể mang đặc điểm của nhóm Ỉ Mỡ nên toàn bộ mẫu khảo sát và thu thập thuộc nhóm Ỉ Pha.
  • Phạm vi không gian: Khảo sát và thu mẫu tại các địa phương thuộc miền Bắc Việt Nam:
    • Lợn Ỉ: Tỉnh Thanh Hóa (xã Hoằng Long - TP. Thanh Hóa, xã Hoằng Đức - huyện Hoằng Hóa, huyện Thiệu Hóa, huyện Yên Định, thị xã Bỉm Sơn, huyện Đông Sơn).
    • Lợn Móng Cái: Thành phố Hải Phòng (Trung tâm Giống Tràng Duệ, huyện An Dương).
    • Lợn Mường Lay: Tỉnh Điện Biên (phường La Nay, phường Sông Đà, xã Lay Nưa thuộc thị xã Mường Lay).
    • Lợn Hương và Lợn Hạ Lang: Tỉnh Cao Bằng (Trung tâm giống cây trồng, vật nuôi, thủy sản Cao Bằng, xã Bạch Đằng, huyện Hòa An).
    • Lợn Mường Khương: Tỉnh Lào Cai (thôn Lao Chải, thôn Sảng Chải, thôn Sín Pao Chải thuộc xã Tả Gia Khâu, huyện Mường Khương).
  • Phạm vi đối sánh: Trình tự hệ gen ty thể của các giống lợn thuộc 6 khu vực địa lý tải từ GenBank (NCBI): Đông Bắc Á, Lưu vực sông Mê Kông, Nam Trung Quốc, Lưu vực sông Dương Tử, Lưu vực sông Hoàng Hà và các nước Châu Âu.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trên thế giới

Tác giả đã tổng hợp các công trình quốc tế về nguồn gốc thuần hóa và ứng dụng hệ gen ty thể (mitochondrial DNA - mtDNA):

  • Nguồn gốc và thuần hóa: Larson và cs. (2005, 2007) xác định lợn nhà bắt nguồn từ lợn rừng nguyên thủy ở Đông Nam Á và chỉ ra sự di cư của lợn theo các đợt mở rộng quần thể ra Thái Bình Dương cũng như sự thuần hóa độc lập tại Cận Đông, Châu Âu và Trung Quốc. Giuffra và cs. (2000) sử dụng đồng hồ phân tử chứng minh lợn nhà có nguồn gốc từ lợn rừng Á - Âu (Sus scrofa) với các nhánh phân ly độc lập.
  • Giải trình tự hệ gen ty thể hoàn chỉnh: Ursing và Arnason (1998) lần đầu tiên giải trình tự hoàn chỉnh mtDNA của Sus scrofa (chiều dài khoảng 16,6 kb). Wang và cs. (2014) công bố dữ liệu hệ gen ty thể giống lợn Duroc. Gongora và cs. (2011) xác định trình tự mtDNA hoàn chỉnh của lợn Visayan và lợn Java (16.475 bp).
  • Phân tích phát sinh chủng loại và đa dạng di truyền: Kim và cs. (2002) phân tích tính đa hình trình tự vùng D-loop dài 1.036 bp để phân định các nhánh lợn Châu Âu và Châu Á; Watanobe và cs. (2003) nghiên cứu cấu trúc quần thể lợn rừng Nhật Bản; Wu và cs. (2007) chỉ ra vị trí phân loại của các giống lợn lưu vực sông Mê Kông và phân nhóm lợn rừng Malaysia; Yu và cs. (2013) nghiên cứu sự phân ly của giống lợn Lanyu và phát triển bộ mồi khuếch đại toàn bộ mtDNA; Fang và cs. (2006) khảo sát mức độ biến thiên nucleotide trên 17 giống lợn bản địa Trung Quốc và 3 giống lợn Châu Âu.
  • Ứng dụng chỉ thị phân tử và mã vạch DNA: Spychaj và cs. (2016), Dawnay và cs. (2007) sử dụng gen ty thể Cytochrome c oxidase subunit I (COI) để định danh loài; Fernández và cs. (2008) phát hiện mối liên hệ giữa đa hình mtDNA với chất lượng thịt ở lợn Iberia.

Các hướng nghiên cứu tại Việt Nam

Tác giả điểm lại các công trình nghiên cứu trong nước:

  • Nghiên cứu hình thái và di truyền quần thể: Hongo và cs. (2002) phân loại nhóm lợn rừng to và lợn nhỏ Việt Nam dựa trên giải phẫu xương sọ; Lan và Shi (1993) so sánh khoảng cách tiến hóa giữa lợn Vân Nam và lợn Việt Nam; Ishihara và cs. (2010) đánh giá 15 giống lợn địa phương bằng chỉ thị SNP; các công trình sử dụng chỉ thị microsatellite (2006, 2014) khẳng định các nhánh lợn bản địa Việt Nam có lịch sử thuần hóa độc lập và lâu dài.
  • Nghiên cứu sinh học phân tử và chỉ thị gen nhân: Nguyễn Thị Phương Mai (2017) xác định vị trí phân loại của heo rừng Tây Nguyên dựa trên 3 SNP vùng D-loop và 3 SNP gen Cytochrome b; Lê Thành Long và cs. (2014) phân tích quan hệ di truyền lợn Tây Nguyên qua gen Cytochrome b; Trần Thị Thúy Nhiên (2016) công bố trình tự hoàn chỉnh hệ gen ty thể lợn Móng Cái (16.632 bp) và phân tích quan hệ phát sinh bằng phương pháp UPGMA. Nhiều nhóm tác giả nghiên cứu các gen chức năng đơn lẻ: Lê Thị Thúy và cs. (2004 - gen Leptin), Đỗ Võ Anh Khoa và cs. (2011 - gen MASP1), Trần Xuân Hoàn và cs. (2013 - gen Mc4R), Đặng Vũ Bình và cs. (2012 - gen Halothane), Nguyễn Thị Diệu Thúy và cs. (2004 - gen hormone sinh trưởng), Nguyễn Thị Hoa (2013 - gen H-FABP, Lpin1), Lê Thị Thu Phương (2003 - gen thụ thể Estrogen), Nguyễn Văn Cường và cs. (2009 - các gen kháng bệnh).

Khoảng trống luận án lựa chọn giải quyết

Tổng quan cho thấy các dữ liệu di truyền về lợn bản địa Việt Nam chủ yếu là các phân đoạn ngắn (Cytochrome b, một phần D-loop) hoặc chỉ tập trung trên từng đối tượng riêng lẻ (như lợn Móng Cái, heo rừng Tây Nguyên). Luận án của Bùi Anh Tuấn lấp đầy khoảng trống này bằng cách tiến hành giải mã đồng thời, hoàn chỉnh toàn bộ hệ gen ty thể ở cả 6 giống lợn đại diện cho miền Bắc Việt Nam, phân tích cấu trúc 37 gen thành phần, dự đoán cấu trúc bậc hai của 22 tRNA, khảo sát toàn bộ các biến dị SNP trên toàn hệ gen và tái dựng cây phát sinh chủng loại bằng phương pháp suy luận Bayes.


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

  • Cấu trúc hệ gen ty thể động vật có vú: mtDNA là phân tử DNA mạch kép dạng vòng, kích thước ~16,6 kb, gồm chuỗi nặng (H) và chuỗi nhẹ (L). Hệ gen mã hóa 13 polypeptide (protein thuộc chuỗi chuyền điện tử hô hấp), 2 rRNA (12S rRNA và 16S rRNA) và 22 tRNA. Vùng D-loop là vùng không mã hóa lớn nhất, chứa promoter phiên mã, điểm khởi đầu sao chép chuỗi H và các khối trình tự bảo thủ (CSB-I, CSB-II, CSB-III) cùng hai đoạn siêu biến (HV1 và HV2). Mã di truyền ty thể động vật có xương sống có sai khác với mã di truyền nhân (TGA mã hóa Tryptophan, ATA mã hóa Methionine, AGA/AGG là bộ ba kết thúc).
  • Lý thuyết phát sinh chủng loại phân tử: Dựa trên sự tích lũy thay thế nucleotide theo thời gian tiến hóa. Sử dụng các mô hình tiến hóa nucleotide (Kimura 2-parameter, Jukes-Cantor, Tamura) và phương pháp suy luận xác suất Bayes kết hợp chuỗi Markov Monte Carlo (MCMC).

Phương pháp thực tế

1. Phương pháp khảo sát ngoại hình và thu mẫu

  • Khảo sát ngoại hình 3 cấp độ: (1) Màu sắc lông da (phân nhóm A+, A, B); (2) Tầm vóc, khối lượng (đo cao vai, dài thân, vòng ngực bằng thước dây, cân khối lượng); (3) Hình dáng cơ thể và số lượng núm vú. Lựa chọn các cá thể cái thuần, điển hình nhất của mỗi giống.
  • Thu mẫu máu vô trùng theo Quy chuẩn kỹ thuật quốc gia QCVN 01-83:2011/BNNPTNT: Lấy 3–5 ml máu từ tĩnh mạch tai, tĩnh mạch cổ hoặc động mạch đuôi, bảo quản trong ống chống đông EDTA ở 10°C khi vận chuyển và lưu giữ tại -20°C.

2. Phương pháp sinh học phân tử trong phòng thí nghiệm

  • Tách chiết DNA tổng số: Áp dụng phương pháp phenol - chloroform - isoamyl alcohol (tỷ lệ 25:24:1) kết hợp dung giải tế bào bằng Lysis Buffer và Proteinase K (2,5 mg/ml) ở 65°C trong 2 giờ theo Sambrook và cs. Kiểm tra độ tinh sạch và định lượng bằng máy đo quang phổ NanoDrop Spectrophotometer và điện di trên gel agarose 1%.
  • Khuếch đại phân đoạn mtDNA bằng PCR: Sử dụng bộ 30 cặp mồi đặc hiệu kế thừa và hiệu chỉnh từ Yu và cs. (2013), đối chiếu với trình tự chuẩn NC000845.1 trên GenBank. Phản ứng thực hiện bằng kit GoTaq Green Master Mix 2X (Promega), thể tích 25 µl trên máy luân nhiệt GeneAmp PCR System 9700. Tinh sạch sản phẩm PCR bằng bộ kit QIAquick PCR Purification Kit (Qiagen).
  • Giải trình tự DNA: Giải trình tự phân đoạn (shotgun) theo nguyên lý kết thúc chuỗi Sanger bằng kit BigDye Terminator v3.1 trên máy giải trình tự tự động ABI 3500 Genetic Analyzer (Thermo Fisher).

3. Phương pháp tin sinh học và phân tích dữ liệu

  • Xử lý và lắp ráp hệ gen: Đánh giá chất lượng tín hiệu sắc ký đồ (.ab1), cắt gọt đoạn chất lượng thấp bằng BioEdit v7.0; lắp ráp các đoạn contig gối đầu lên nhau bằng phần mềm DNA Dragon v1.0 và EditSeq để tạo trình tự mạch vòng hoàn chỉnh dạng FASTA.
  • Dóng hàng đa trình tự: Sử dụng thuật toán MUSCLE tích hợp trong phần mềm MEGA7. Loại bỏ các đoạn lặp ngắn 5'-CGTGCGTACA-3' trong vùng D-loop do có hiện tượng dị chất (heteroplasmy) không mang tín hiệu phát sinh loài.
  • Chú giải hệ gen và cấu trúc RNA: Sử dụng công cụ trực tuyến MITOS Web Server và DOGMA, đối chiếu BLAST trên NCBI RefSeq. Dự đoán cấu trúc bậc hai của 22 gen tRNA bằng thuật toán mô hình hiệp phương sai MITFI. Trực quan hóa bản đồ vòng hệ gen ty thể bằng GenomeVx.
  • Phân tích thành phần và đa hình: Tính tỷ lệ nucleotide, độ lệch AT-skew [(A-T)/(A+T)] và GC-skew [(G-C)/(G+C)] bằng DAMBE v6.0. Xác định vị trí đa hình đơn nucleotide (SNP) bằng DnaSP v5.10.
  • Xây dựng cây phát sinh chủng loại: Tính toán khoảng cách cặp trình tự ($p$-distance) theo mô hình Kimura 2 tham số. Tái dựng cây phát sinh chủng loại bằng phương pháp suy luận Bayes trên phần mềm BEAST v1.8, tìm cây tối ưu qua TreeAnnotator, hiển thị bằng FigTree v1.4.2 và kiểm tra giá trị bootstrap (1000 lần lặp) trên MEGA7. Sử dụng trình tự lợn hoang Malaysia (Sus barbatus, mã số EF545592.1) làm nhóm ngoại đối chứng (outgroup) để định gốc cho cây.

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

  • Mẫu sinh học: 6 mẫu máu từ 6 cá thể cái đại diện của 6 giống lợn bản địa miền Bắc.
  • Dữ liệu thứ cấp: Trình tự hệ gen ty thể chuẩn của lợn (Sus scrofa, NC000845.1) và 25 trình tự tham chiếu của các giống lợn Châu Á, Châu Âu từ cơ sở dữ liệu GenBank (NCBI).
Khu vực địa lý Tên giống lợn (Tiếng Việt / Tiếng Anh) Mã số truy cập GenBank
Đông Bắc Á Lợn rừng Hàn Quốc (WB-Korea) AY574047.1
Lợn rừng Đông Bắc Trung Quốc (WB-China northeast) EU333163.1
Khu vực Mê Kông Lợn rừng Việt Nam (WB-Vietnam) EF545584.1
Lợn Banna mini (Banna mini) GQ220328.1
Lợn Đa-hê (Dahe) GQ220329.1
Lợn rừng Malaysia (WB-Malaysia - Outgroup) EF545592.1
Nam Trung Quốc Lợn rừng Vân Nam (WB-Yunnan) EF545573.1
Lợn rừng Phúc Kiến (WB-Fujian) EF545569.1
Lợn rừng Hải Nam (WB-Hainan) EF545572.1
Lợn Lang-tan (Lantang) KC250274
Lợn La-ny-u (Lanyu) DQ518915.2
Lưu vực sông Dương Tử Lợn A-ba (Aba) EF545578.1
Lợn rừng Jiangxi (WB-Jiangxi) EF545579.1
Lợn Kim Hoa (Jinhua) KC469586
Lợn Xi-ang (Xiang) KC250273
Lưu vực sông Hoàng Hà Lợn Bi-hu (Bihu) EF545591.1
Lợn Wei (Wei) EF545577.1
Lợn Ba-mây (Bamei) EF545583.1
Lợn Hu-du (Huzu) EF545588.1
Các nước Châu Âu Lợn Béc-sai (Berkshire) AY574045.1
Lợn Đu-rốc / Heo bò (Duroc) AY337045.1
Lợn đen I-bê-ria (Iberian) FJ236994.1
Lợn Ham-p-shơ (Hampshire) AY574046.1
Lợn Đại Bạch (Large-white) KC250275
Lợn Lan-đờ-rát (Landrace) AF034253.1
Lợn Pi-ê-tờ-ranh (Pietrain) KC469587
Lợn rừng Châu Âu (WB-European) FJ237000.1

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

Chương 1: Tổng quan tài liệu

Chương 1 hệ thống hóa cơ sở lý thuyết và các nghiên cứu liên quan đến lợn nhà và hệ gen ty thể:

  • Lịch sử thuần hóa lợn trải qua 5 sự kiện chính: hình thành loài Sus tại Đông Nam Á (~4-2 triệu năm trước), phân ly thành dòng Châu Âu và Châu Á (~1,2 triệu năm trước), thuần hóa độc lập (~10.000 năm trước), tạp giao giữa lợn Châu Âu và Châu Á (~200 năm trước) và quá trình chọn tạo các giống lợn hiện đại.
  • Đặc điểm cấu trúc hệ gen ty thể ở động vật có vú: DNA vòng kép, kích thước ~16,6 kb, chứa 37 gen (13 gen mã hóa protein, 22 tRNA, 2 rRNA) và vùng kiểm soát D-loop. Sự khác biệt của mã di truyền ty thể và cơ chế phiên mã, sao chép độc lập với nhân.
  • Ứng dụng của mtDNA trong phân loại học, mã vạch DNA (sử dụng gen COI, 12S rRNA, 16S rRNA, Cytochrome b) và nghiên cứu phát sinh chủng loại học phân tử.
  • Các phương pháp dựng cây tiến hóa: nhóm dựa trên khoảng cách (UPGMA, Neighbor-Joining), Maximum Parsimony, Maximum Likelihood và suy luận Bayes (MCMC).
  • Điểm lại các công trình nghiên cứu mtDNA trên thế giới và tại Việt Nam.
  • Tổng quan hiện trạng phân bố, đặc điểm ngoại hình và khả năng sinh trưởng của 26 giống lợn bản địa Việt Nam, tập trung vào 6 giống nghiên cứu: lợn Ỉ, Móng Cái, Hạ Lang, Mường Lay, Mường Khương và Hương.

Chương 2: Đối tượng và phương pháp nghiên cứu

Chương 2 mô tả chi tiết vật liệu, hóa chất, thiết bị và các bước thực nghiệm:

  • Tiêu chí lựa chọn cá thể đại diện dựa trên khảo sát hình thái lông da, tầm vóc và hình dáng cơ thể tại các tỉnh Thanh Hóa, Hải Phòng, Điện Biên, Cao Bằng, Lào Cai.
  • Hóa chất, sinh phẩm và danh mục 30 cặp mồi PCR đặc hiệu khuếch đại phân đoạn phủ kín toàn bộ mtDNA lợn (nhiệt độ gắn mồi $T_a$ từ 53°C đến 55°C).
  • Quy trình ly trích DNA bằng Phenol:Chloroform:Isoamyl alcohol (25:24:1) và định lượng Nanodrop.
  • Quy trình phản ứng PCR (chu trình 25 chu kỳ) và tinh sạch bằng cột silica.
  • Quy trình giải trình tự Sanger phân đoạn trên máy ABI 3500 với kit BigDye Terminator v3.1.
  • Quy trình tin sinh học: BioEdit đánh giá sắc ký đồ; DNA Dragon và EditSeq lắp ráp contig; MITOS và DOGMA chú giải gen; MITFI dự đoán cấu trúc bậc hai tRNA; DAMBE tính tỷ lệ A, T, G, C và AT/GC-skew; DnaSP phân tích đa hình nucleotide; BEAST, TreeAnnotator, FigTree và MEGA7 dựng và đánh giá cây phát sinh chủng loại.

Chương 3: Kết quả và thảo luận

Chương 3 trình bày các kết quả thực nghiệm và biện giải khoa học:

  • Kết quả khảo sát ngoại hình và chọn mẫu: Phân loại cá thể theo 3 cấp độ; lựa chọn 6 cá thể cái mang đầy đủ tính trạng đặc trưng nhất của giống (thuộc nhóm A/A+) để tách chiết DNA.
  • Tách chiết DNA tổng số và khuếch đại PCR: DNA tổng số của 6 giống lợn có độ tinh sạch cao (tỷ lệ $A_{260}/A_{280}$ đạt từ 1,8 đến 2,0), điện di gel agarose cho vệt DNA gọn, không đứt gãy. Cả 30 cặp mồi đều khuếch đại thành công các phân đoạn phân tử có kích thước từ 384 bp đến 1.200 bp với độ đặc hiệu cao.
  • Lắp ráp và xác định kích thước hệ gen ty thể hoàn chỉnh: Đã lắp ráp thành công hệ gen ty thể dạng vòng hoàn chỉnh cho cả 6 giống lợn bản địa:
    • Lợn Móng Cái: 16.632 bp
    • Lợn Mường Lay: 16.635 bp
    • Lợn Mường Khương: 16.679 bp
    • Lợn Hạ Lang: 16.637 bp
    • Lợn Hương: 16.637 bp
    • Lợn Ỉ: 16.639 bp
  • Thành phần nucleotide và độ lệch: Hệ gen ty thể của cả 6 giống lợn giàu A+T (chiếm khoảng 60,3% đến 60,5%), tỷ lệ các loại nucleotide: A (~34,6%), T (~25,8%), C (~26,2%), G (~13,4%). Giá trị AT-skew dương (dao động 0,145 đến 0,147) và GC-skew âm (dao động -0,323 đến -0,326), phản ánh sự bất đối xứng thành phần nucleotide điển hình giữa chuỗi H và chuỗi L ở động vật có vú.
  • Chú giải cấu trúc hệ gen: Hệ gen ty thể của cả 6 giống lợn đều bảo toàn cấu trúc chuẩn gồm 37 gen (13 gen mã hóa protein, 22 gen tRNA, 2 gen rRNA) và 1 vùng điều hòa D-loop.
    • 14 gen tRNA, 2 gen rRNA và 12 gen mã hóa protein nằm trên chuỗi nặng (H); riêng gen ND6 và 8 gen tRNA (tRNA-Gln, tRNA-Ala, tRNA-Asn, tRNA-Cys, tRNA-Tyr, tRNA-Ser(UCN), tRNA-Glu, tRNA-Pro) nằm trên chuỗi nhẹ (L).
    • Khảo sát bộ ba mở đầu và kết thúc: Đa số các gen bắt đầu bằng ATG, một số gen dùng ATA hoặc GTG; kết thúc bằng TAA, TAG hoặc bộ ba không hoàn chỉnh T-- được bổ sung đuôi polyA sau phiên mã.
    • Cấu trúc bậc hai của 22 tRNA: Hầu hết có cấu trúc nhánh 4 cánh điển hình hình lá chẽ ba (cloverleaf), riêng tRNA-Ser(AGY) thiếu nhánh dihydrouridine (DHU-arm).
  • Phân tích đa hình trình tự (SNP):
    • Vùng D-loop: Xác định được nhiều vị trí biến dị đơn nucleotide và các đoạn chèn/mất (indels), tập trung tại các vùng siêu biến HV1 và HV2.
    • Vùng mã hóa: Xác định các điểm đột biến nucleotide phân biệt giữa các giống lợn bản địa Việt Nam với các dòng lợn Châu Âu và Trung Quốc (như các vị trí biến dị đặc thù tại nucleotide 2250, 11318).
  • Phân tích cây phát sinh chủng loại:
    • Cây phát sinh dựa trên vùng D-loop và cây dựa trên trình tự hoàn chỉnh hệ gen ty thể đều phân chia rõ rệt thành hai nhánh chính: nhánh lợn Châu Âu và nhánh lợn Châu Á.
    • Cả 6 giống lợn bản địa Việt Nam đều nằm trọn vẹn trong nhánh lợn Châu Á, có quan hệ di truyền gần gũi với các giống lợn thuộc lưu vực sông Mê Kông và Nam Trung Quốc.
    • Trình tự hệ gen ty thể của giống lợn Hương và lợn Hạ Lang có mức độ tương đồng trình tự rất cao (gần như trùng khít về kích thước 16.637 bp, vị trí các biến dị và cùng nằm trên một nhánh phát sinh với hệ số bootstrap/xác suất hậu nghiệm đạt 100%), từ đó tác giả xây dựng giả thuyết hai giống này có thể có cùng một nguồn gốc di truyền dòng mẹ.

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

Những đóng góp mới của luận án

  1. Đã xác định và công bố trình tự hoàn chỉnh hệ gen ty thể của 6 giống lợn bản địa Việt Nam (lợn Ỉ, lợn Móng Cái, lợn Mường Khương, lợn Mường Lay, lợn Hương và lợn Hạ Lang) lên Ngân hàng Gen quốc tế (GenBank - NCBI).
  2. Đã phân tích thành phần nucleotide, chú giải đầy đủ trật tự vị trí 37 gen (13 gen mã hóa protein, 22 gen tRNA, 2 gen rRNA và vùng D-loop), dự đoán cấu trúc bậc hai của 22 phân tử tRNA cho cả 6 giống lợn bản địa.
  3. Cung cấp bộ dữ liệu so sánh đa hình nucleotide, xác định các vị trí SNP đặc trưng trên vùng D-loop và vùng mã hóa của hệ gen ty thể, làm sáng tỏ mức độ sai khác di truyền giữa lợn bản địa Việt Nam với các giống lợn Châu Á và Châu Âu.
  4. Đưa ra các bằng chứng phân tử xác định quan hệ phát sinh chủng loại dòng mẹ của 6 giống lợn bản địa miền Bắc, đồng thời xây dựng giả thuyết khoa học về nguồn gốc chung giữa giống lợn Hương và giống lợn Hạ Lang (tỉnh Cao Bằng).

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

  • Đưa bộ dữ liệu trình tự mtDNA hoàn chỉnh của 6 giống lợn vào hệ thống cơ sở dữ liệu quốc gia về nguồn gen vật nuôi để phục vụ công tác giám định, nhận dạng giống thuần và đánh giá suy thoái đồng huyết.
  • Sử dụng các chỉ thị SNP ty thể đã phát hiện để thiết lập các công cụ phân tử hỗ trợ chọn giống và loại bỏ các cá thể lai tạp trong các chương trình nhân thuần bảo tồn lợn bản địa (đặc biệt là quần thể lợn Ỉ và lợn Mường Lay).
  • Tiếp tục mở rộng các nghiên cứu giải mã hệ gen hoàn chỉnh trên các giống lợn bản địa còn lại của Việt Nam tại khu vực miền Trung, Tây Nguyên và Nam Bộ.

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

Hạn chế của nghiên cứu

  • Về cỡ mẫu giải trình tự hoàn chỉnh: Luận án mới chỉ thực hiện giải trình tự toàn bộ hệ gen ty thể trên 1 cá thể đại diện tiêu biểu (đã qua sàng lọc hình thái 3 cấp độ nghiêm ngặt) cho mỗi giống trong số 6 giống lợn nghiên cứu; chưa khảo sát biến dị mtDNA trên diện rộng quần thể của từng giống.
  • Về phạm vi đối tượng: Đề tài tập trung vào 6 giống lợn phân bố chủ yếu ở khu vực miền Bắc Việt Nam, chưa bao quát toàn bộ 26 giống lợn bản địa trên cả nước.
  • Về phạm vi vật chất di truyền: Nghiên cứu mới chỉ tiếp cận ở mức độ hệ gen ty thể (di truyền theo dòng mẹ), chưa kết hợp phân tích hệ gen nhân tế bào để đánh giá toàn diện mức độ lai tạp và hiện tượng du nhập dòng bố (paternal gene flow).

Hướng nghiên cứu tiếp theo

  • Mở rộng giải trình tự hệ gen ty thể hoàn chỉnh cho các cá thể khác trong cùng quần thể để đánh giá độ đa dạng haplotype nội giống.
  • Triển khai nghiên cứu giải mã hệ gen ty thể trên các giống lợn bản địa đặc thù khác ở miền Trung và miền Nam (như lợn Sóc, lợn Ba Xuyên, lợn Cỏ, lợn Vân Pa).
  • Kết hợp dữ liệu hệ gen ty thể với giải trình tự toàn bộ hệ gen nhân (whole-genome sequencing) để phân tích các gen chức năng liên quan đến chất lượng thịt thơm ngon, tính kháng bệnh và khả năng thích nghi sinh thái.

Giá trị tham khảo

  • Đối với nghiên cứu sinh và giảng viên chuyên ngành Di truyền học, Công nghệ sinh học và Chăn nuôi:
    • Tài liệu tham khảo về quy trình kỹ thuật 30 cặp mồi PCR khuếch đại bao phủ toàn bộ hệ gen ty thể động vật có vú.
    • Tham khảo quy trình tin sinh học tích hợp: từ xử lý sắc ký đồ Sanger, ghép nối contig, chú giải tự động (MITOS, DOGMA), dự đoán cấu trúc bậc hai RNA đến phân tích cây phát sinh chủng loại bằng suy luận Bayes trên BEAST.
  • Đối với các cơ quan quản lý và nhà làm chính sách bảo tồn nguồn gen vật nuôi:
    • Cung cấp cơ sở dữ liệu phân tử chuẩn mực để nhận diện nguồn gen thuần chủng, phục vụ công tác giám định giống, xây dựng đề án bảo tồn nguồn gen quốc gia cho các giống lợn nguy cấp (lợn Ỉ, lợn Mường Lay 14 vú).
    • Cung cấp dữ liệu phục vụ đăng ký sở hữu trí tuệ, chỉ dẫn địa lý và phát triển thương hiệu đặc sản bản địa cho các địa phương (Cao Bằng, Điện Biên, Lào Cai, Quảng Ninh).

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

1. Kích thước hệ gen ty thể hoàn chỉnh của 6 giống lợn bản địa được xác định trong luận án là bao nhiêu?

Kích thước hệ gen ty thể hoàn chỉnh của 6 giống lợn bản địa Việt Nam dao động từ 16.632 bp đến 16.679 bp. Cụ thể: lợn Móng Cái (16.632 bp), lợn Mường Lay (16.635 bp), lợn Hương (16.637 bp), lợn Hạ Lang (16.637 bp), lợn Ỉ (16.639 bp) và lợn Mường Khương (16.679 bp).

2. Luận án đã sử dụng phương pháp và bộ mồi nào để giải mã toàn bộ hệ gen ty thể của 6 giống lợn?

Tác giả sử dụng phương pháp giải trình tự phân đoạn (shotgun) Sanger trên hệ thống máy giải trình tự tự động ABI 3500 với kit BigDye Terminator v3.1. Toàn bộ hệ gen được khuếch đại thông qua 30 phản ứng PCR sử dụng bộ 30 cặp mồi đặc hiệu gối đầu nhau, được hiệu chỉnh từ nghiên cứu của Yu và cs. (2013) dựa trên trình tự chuẩn NC000845.1.

3. Những phần mềm và công cụ tin sinh học chính nào được tác giả sử dụng trong phân tích hệ gen?

Luận án sử dụng hệ thống phần mềm gồm:

  • Xử lý và lắp ráp: BioEdit v7.0, DNA Dragon v1.0, EditSeq.
  • Chú giải gen và cấu trúc: MITOS Web Server, DOGMA, MITFI, GenomeVx.
  • Dóng hàng và tính toán thống kê: MEGA7 (thuật toán MUSCLE), DAMBE v6.0, DnaSP v5.10.
  • Phát sinh chủng loại: BEAST v1.8, TreeAnnotator, FigTree v1.4.2 và MEGA7.

4. Kết quả phân tích phát sinh chủng loại cho thấy quan hệ nguồn gốc của 6 giống lợn bản địa Việt Nam như thế nào?

Cây phát sinh chủng loại suy luận Bayes khẳng định cả 6 giống lợn bản địa Việt Nam đều thuộc nhánh lợn Châu Á, có mối quan hệ gần gũi với các nhóm lợn thuộc lưu vực sông Mê Kông và Nam Trung Quốc, tách biệt hoàn toàn với nhóm lợn Châu Âu. Đặc biệt, giống lợn Hương và lợn Hạ Lang có trình tự mtDNA gần như tương đồng hoàn toàn, chỉ ra giả thuyết hai giống này có chung nguồn gốc di truyền dòng mẹ.

5. Tại sao tác giả lựa chọn trình tự lợn hoang Malaysia (Sus barbatus) làm nhóm ngoại đối chứng (outgroup) khi xây dựng cây phát sinh chủng loại?

Trình tự mtDNA của lợn hoang Malaysia (Sus barbatus, mã số GenBank EF545592.1) được chọn làm nhóm ngoại vì loài này đã phân tách tiến hóa từ sớm so với các quần thể lợn hoang và lợn nhà Châu Á, Châu Âu. Đây là tiêu chuẩn đã được áp dụng trong các nghiên cứu phát sinh chủng loại quốc tế trước đó (như Wu và cs., 2007; Yu và cs., 2013) nhằm giúp định gốc cây tiến hóa chính xác.


Kết luận

Luận án của Bùi Anh Tuấn đã xác định và công bố hoàn chỉnh trình tự hệ gen ty thể của 6 giống lợn bản địa miền Bắc Việt Nam (lợn Ỉ, Móng Cái, Mường Khương, Mường Lay, Hương và Hạ Lang) với kích thước từ 16.632 bp đến 16.679 bp. Công trình đã chú giải chi tiết 37 gen thành phần, dự đoán cấu trúc bậc hai của 22 gen tRNA và phân tích toàn diện các vị trí đa hình đơn nucleotide trên vùng D-loop và vùng mã hóa. Kết quả phân tích phát sinh chủng loại bằng suy luận Bayes khẳng định các giống lợn bản địa Việt Nam thuộc dòng tiến hóa lợn Châu Á, đồng thời cung cấp bằng chứng di truyền học phân tử quan trọng phục vụ công tác nhận dạng, đánh giá phả hệ và bảo tồn bền vững nguồn gen vật nuôi bản địa.