Giới thiệu dự án

Nghiên cứu nguồn gốc tiến hóa và đa dạng di truyền của động vật thuần hóa đóng vai trò nền tảng trong bảo tồn nguồn gen bản địa và chọn giống sinh học. Loài chó (Canis lupus familiaris) được thuần hóa từ chó sói xám cổ đại cách đây khoảng 11.000 đến 100.000 năm, hiện có hơn 400 giống chó thuần chủng cùng hàng triệu biến thể lai trên toàn cầu. Tại Việt Nam, "Tứ đại danh khuyển" bao gồm chó Phú Quốc, chó cộc đuôi H'Mông, chó Bắc Hà và chó Lài (Dingo Đông Dương) đại diện cho nguồn tài nguyên di truyền đặc hữu quý giá. Tuy nhiên, sự gia tăng nhập khẩu các giống chó ngoại (Becgie, Husky, Ngao Tây Tạng, Chihuahua) cùng tập quán nuôi thả tự do đã dẫn đến nguy cơ lai tạp, xói mòn nguồn gen bản địa.

Phần lớn các nghiên cứu trong nước trước đây tập trung vào bệnh lý học thú y hoặc đặc điểm di truyền riêng biệt của chó xoáy Phú Quốc trên đảo. Quần thể chó nhà (Canis lupus familiaris) sinh sống tại đất liền khu vực Tây Nam Bộ – đặc biệt là thành phố Rạch Giá (tỉnh Kiên Giang), cửa ngõ giao thương trực tiếp với đảo Phú Quốc – vẫn chưa được phân tích di truyền phân tử hệ thống. Khoảng trống dữ liệu này làm hạn chế khả năng đối sánh di truyền giữa quần thể chó đất liền và chó đảo Phú Quốc, gây khó khăn cho việc truy nguyên dòng di cư và xác thực tính đặc hữu của nguồn gen chó bản địa.

       QUẦN THỂ CHÓ ĐẤT LIỀN (RẠCH GIÁ)                 QUẦN THỂ ĐẢO PHÚ QUỐC
  [Chó nhà bản địa: Đa hình kiểu hình]              [Chó Phú Quốc: Lưng xoáy đặc thù]
      Tách chiết DNA từ gốc lông                        Tách chiết DNA từ gốc lông
      Khuếch đại PCR vùng CR (1.232 bp)                 Khuếch đại PCR vùng CR (1.232 bp)
      Giải trình tự Sanger hai chiều                     Giải trình tự Sanger hai chiều
                           Phát hiện dòng Haplotype E4 hiếm
                           Chứng minh giao lưu nguồn gen qua lại

Đề tài tốt nghiệp Kỹ sư Công nghệ Sinh học "Đánh giá sự đa dạng di truyền của quần thể chó nhà tại Rạch Giá dựa vào vùng trình tự CR trên bộ gen ty thể" được thực hiện bởi tác giả Lê Quế Phương dưới sự hướng dẫn của ThS. Huỳnh Văn Hiếu tại Viện Bảo tồn Nguồn gen NTT (Đại học Nguyễn Tất Thành). Mục tiêu nghiên cứu bao gồm:

  1. Thu thập 23 mẫu lông chó nhà phân bố ngẫu nhiên, không cùng huyết thống tại TP. Rạch Giá, Kiên Giang.
  2. Chuẩn hóa quy trình tách chiết DNA tổng số chất lượng cao từ gốc chân lông (không xâm lấn).
  3. Tối ưu hóa phản ứng PCR khuếch đại phân đoạn 1.232 bp chứa vùng siêu biến 1 (HV1) thuộc vùng kiểm soát (Control Region - CR) trên hệ gen ty thể (mtDNA).
  4. Giải trình tự hai chiều Sanger phân đoạn HV1 và hiệu chỉnh chromatogram loại bỏ nhiễu kỹ thuật.
  5. Định vị các đa hình đơn nucleotide (SNP) trên đoạn lõi 582 bp và phân loại Haplotype/Haplogroup thông qua đối sánh dữ liệu NCBI GenBank.
  6. Tái dựng cây phát sinh chủng loài bằng giải thuật Neighbor-Joining (NJ) nhằm làm sáng tỏ mối quan hệ di truyền giữa chó nhà Rạch Giá, chó Phú Quốc và các quần thể chó trên thế giới.

Phương pháp tiếp cận dựa trên hệ gen ty thể (mtDNA) với chiều dài 16.727 bp mang ưu thế vượt trội: tốc độ đột biến cao gấp 10–25 lần so với DNA nhân, di truyền đơn bội theo dòng mẹ tuyệt đối (matrilineal inheritance) và không xảy ra hiện tượng tái tổ hợp di truyền. Nghiên cứu tập trung vào vùng kiểm soát CR (dài ~1.270 bp), đặc biệt là vùng siêu biến HV1 (Hypervariable Region 1, vị trí 15.458–16.130 bp) chứa đoạn đột biến tập trung (hotspot 15.460–15.520 bp).

Nghiên cứu ghi nhận 100% mẫu phân tích đạt độ tinh sạch cao ($OD_{260}/OD_{280} = 1{,}62 - 2{,}10$), khuếch đại thành công dải băng đơn 1.232 bp, xác định 9 haplotype (A11, A17, A18, A73, B1, B3, C2, C3, E4) thuộc 4 haplogroup chính (A: 60%, B: 10%, C: 20%, E: 10%). Kết quả phát hiện Haplotype E4 (tỷ lệ 10%) – dòng gen hiếm toàn cầu chỉ chiếm <2% (chủ yếu tại Đông Á) – khẳng định sự giao lưu di truyền chặt chẽ giữa chó đất liền Rạch Giá và chó xoáy Phú Quốc. Nghiên cứu giới hạn trong phạm vi 23 cá thể đại diện tại Rạch Giá và chỉ thị phân tử dòng mẹ mtDNA.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Để đánh giá mức độ đa dạng di truyền và tiến hóa dưới loài, nhiều hệ thống chỉ thị phân tử đã được áp dụng trong sinh học phân tử hiện đại:

Chỉ thị / Phương pháp Độ phân giải phân loại Cơ chế di truyền Ưu điểm kỹ thuật Nhược điểm & Hạn chế
Gen 16S rRNA (mtDNA) Bậc loài / chi Dòng mẹ Vùng bảo thủ cao, dễ khuếch đại PCR Độ biến dị thấp, không phân biệt được các thứ/nòi dưới loài
Microsatellite (STRs nhân) Bậc cá thể / quần thể Đồng trội (Mendel) Phân tích quan hệ huyết thống 2 dòng bố mẹ Bị tái tổ hợp, khó truy nguyên phát sinh chủng loài cổ xưa
Toàn bộ hệ gen (WGS/NGS) Siêu phân giải Toàn diện Phát hiện toàn bộ biến dị cấu trúc, SNP nhân Chi phí rất cao, đòi hỏi hạ tầng tính toán phức tạp
Vùng HV1 mtDNA (Giải pháp chọn) Bậc nòi / Haplotype Dòng mẹ không tái tổ hợp Đột biến nhanh (10-25x), chi phí tối ưu, giàu SNP chỉ thị Chỉ phản ánh dòng mẹ, không phát hiện dòng lai từ bố

Mô hình yêu cầu hệ thống theo phân loại MoSCoW:

  • Must have (Bắt buộc): Tách chiết DNA nguyên vẹn từ mẫu sinh thiết không xâm lấn (gốc lông); khuếch đại thành công phân đoạn CR 1.232 bp; đọc chuỗi Sanger 2 chiều; nhận diện chính xác các điểm SNP trên phân đoạn chuẩn 582 bp; định danh Haplotype theo chuẩn quốc tế.
  • Should have (Nên có): Thiết lập cây phát sinh chủng loài Neighbor-Joining với hệ số lặp lại Bootstrap $\ge 1.000$; so sánh tần số haplogroup với dữ liệu toàn cầu của Pang et al. (2009).
  • Could have (Có thể có): Mở rộng khảo sát đoạn lặp VNTR (Variable Number Tandem Repeat) giữa vùng HV1 và HV2 (~300 bp).
  • Won't have (Không làm đợt này): Giải mã toàn bộ hệ gen nhân hoặc phân tích nhiễm sắc thể giới tính Y.
graph TD
    A["Thu thập mẫu lông tại Rạch Giá (n=23)"] --> B["Tách chiết DNA (Proteinase K + PCI 25:24:1)"]
    B --> C["Kiểm tra OD260/OD280 & Điện di Gel Agarose 1%"]
    C --> D["Khuếch đại PCR Vùng CR (Cặp mồi 15412F/16625R ~1.232 bp)"]
    D --> E["Giải trình tự Sanger 2 chiều (Biotex Nam Khoa)"]
    E --> F["Hiệu chỉnh Chromatogram (FinchTV v1.4.0)"]
    F --> G["Lắp ráp Consensus & Gióng hàng (SeaView v4.4)"]
    G --> H["Truy vấn BLAST đối chiếu NCBI GenBank"]
    H --> I["Phân tích SNP & Định danh Haplotype"]
    I --> J["Tái dựng cây phát sinh loài (MEGA v7.0 - NJ / Kimura-2P)"]

Thiết kế hệ thống

Toàn bộ quy trình từ mẫu sinh học đến cây phân loại phát sinh loài được chuẩn hóa thành pipeline tích hợp phần cứng và công cụ tin sinh học:

Công cụ và thiết bị sử dụng:

  • Hóa chất & Enzyme: HS Prime Taq DNA Polymerase 2X MasterMix (GeNet Bio, Hàn Quốc - Mã Cat: G-7000); Proteinase K 20 mg/mL (Genet Bio); Phenol:Chloroform:Isoamyl alcohol (25:24:1, Merck KGaA, Đức); Agarose M (Merck KGaA); Thang chuẩn DNA Ladder 100 bp - 1.500 bp (Bio Basic Inc.).
  • Cặp mồi chuyên biệt khuếch đại vùng CR:
    • Mồi xuôi 15412F: 5'-CCACTATCAGCACCCAAAG-3' (19 nt, $T_m = 53{,}2^\circ\text{C}$).
    • Mồi ngược 16625R: 5'-AGACTACGAGACCAAATGCG-3' (20 nt, $T_m = 54{,}6^\circ\text{C}$).
  • Phần mềm Tin sinh học:
    • FinchTV v1.4.0 (Geospiza Inc.): Kiểm tra đồ thị huỳnh quang (electropherogram), phân tích đỉnh tín hiệu, cắt tỉa vùng nhiễu đầu/cuối mồi, hiệu chỉnh lỗi base-calling.
    • SeaView v4.4: Căn chỉnh đa chuỗi (Multiple Sequence Alignment - ClustalW/Muscle engine), trích xuất trình tự đồng thuận (Consensus sequence) từ 2 mạch xuôi-ngược.
    • MEGA v7.0 (Molecular Evolutionary Genetics Analysis): Đánh giá mô hình thay thế nucleotide tối ưu theo tiêu chuẩn AICc (Akaike Information Criterion corrected), ước tính khoảng cách di truyền Kimura 2-parameter (K2P), xây dựng cây quan hệ phát sinh loài bằng thuật toán Neighbor-Joining (NJ) với 1.000 bước lặp Bootstrap.

Cấu trúc dữ liệu phân tích chuẩn (FASTA Header Format):

>Sample_ID|Species|Location|Target_Region|Length
>R01_Canis_lupus_familiaris|Rach_Gia_VN|mtDNA_HV1|582bp
CCACTATCAGCACCCAAAGCTGAAATTCTACTTAAACTATTCCCTGACACCCCTACATTC
...

Methodology

Nghiên cứu kéo dài 7 tháng (02/2017 – 08/2017) tại Phòng thí nghiệm Viện Bảo tồn Nguồn gen NTT. Khung quản lý rủi ro và đảm bảo chất lượng:

  • Kiểm soát tạp nhiễm chéo: Sử dụng đầu tip có màng lọc (filter tips), khử trùng dụng cụ bằng cồn $96^\circ$, thao tác trong tủ cấy vô trùng HUY HOANG 475, chạy đối chứng âm ($H_2O$ thay DNA khuôn) trong toàn bộ mẻ PCR.
  • Xử lý mẫu suy thoái: Chân lông chứa hàm lượng DNA thấp; áp dụng phương pháp ủ ly giải 2 bước với Proteinase K ($55^\circ\text{C}$ qua đêm) kết hợp tăng nồng độ kết tủa NaCl 5M và Ethanol tuyệt đối lạnh gấp 3,5 lần thể tích dịch nổi.
  • Độ tin cậy giải trình tự: Giải trình tự độc lập hai chiều (Bidirectional Sanger Sequencing). Chỉ những vị trí nucleotide có đỉnh huỳnh quang rõ ràng, phân tách sắc nét và đồng nhất trên cả 2 mạch mới được đưa vào ma trận phân tích SNP.

Implementation và kết quả

Development process

Quá trình triển khai kỹ thuật được chia thành 4 pha thực nghiệm nghiêm ngặt:

Pha 1: Tách chiết DNA tổng số từ gốc lông (Cải tiến giao thức)

  1. Cắt 50–60 sợi lông (lấy đoạn 0,5–1,0 cm tính từ chân lông) đưa vào ống Eppendorf 2,0 mL.
  2. Bổ sung $200,\mu\text{L}$ đệm chiết (Tris-HCl 10 mM pH 8,0; Triton X-100 1%; SDS 1%; EDTA 10 mM), vortex 30 giây, ủ $50^\circ\text{C}$ trong 20 phút. Lặp lại với $100,\mu\text{L}$ đệm chiết.
  3. Bổ sung $5,\mu\text{L}$ Proteinase K ($20,\text{mg/mL}$), vortex 30 giây, ủ ấm qua đêm ($55^\circ\text{C}$).
  4. Chiết tách hữu cơ: Bổ sung 1 thể tích PCI (25:24:1), ly tâm $14.000,\text{rpm}$ trong 5 phút ở $4^\circ\text{C}$. Thu pha nước trên, chiết tiếp với 1 thể tích Chloroform:Isoamyl alcohol (24:1).
  5. Kết tủa: Bổ sung NaCl 5M và Ethanol 100% lạnh (tỷ lệ 1:3,5 thể tích dịch nổi), ủ lạnh sâu $-20^\circ\text{C}$ qua đêm.
  6. Thu cặn qua ly tâm $14.000,\text{rpm}$ trong 5 phút, rửa tủa bằng Ethanol 70% lạnh, làm khô tự nhiên và hoàn nguyên trong $30,\mu\text{L}$ nước cất vô trùng khử ion.

Pha 2: Thiết lập phản ứng PCR khuếch đại vùng CR (1.232 bp)

Thành phần phản ứng ($25,\mu\text{L}$): $2{,}5,\mu\text{L}$ 10X Reaction Buffer ($20,\text{mM MgCl}_2$), $0{,}625,\mu\text{L}$ dNTPs 10 mM, $0{,}25,\mu\text{L}$ HS Prime Taq DNA Polymerase ($2{,}5,\text{U}/\mu\text{L}$), $0{,}25,\mu\text{L}$ mồi 15412F ($10,\mu\text{M}$), $0{,}25,\mu\text{L}$ mồi 16625R ($10,\mu\text{M}$), $2{,}0,\mu\text{L}$ DNA khuôn ($~50-180,\text{ng}/\mu\text{L}$), $19{,}125,\mu\text{L}$ nước siêu sạch.

Do đặc tính sinh học của các mẫu lông, nhiệt độ bắt cặp ($T_a$) được phân nhóm và tối ưu hóa thành công:

  • Nhóm 1 (16 mẫu: R1, R2, R4–R9, R12, R13, R16–R21): $T_a = 49^\circ\text{C}$ (35 chu kỳ).
  • Nhóm 2 (4 mẫu: R10, R11, R14, R15): $T_a = 54^\circ\text{C}$ (35 chu kỳ).
  • Nhóm 3 (3 mẫu: R3, R22, R23): $T_a = 55^\circ\text{C}$ (35 chu kỳ).

Chu trình nhiệt chuẩn: Tiền biến tính $94^\circ\text{C}$ (10 phút); 35 chu kỳ gồm [Biến tính $94^\circ\text{C}$ (30 giây) $\rightarrow$ Bắt cặp $49-55^\circ\text{C}$ (30 giây) $\rightarrow$ Kéo dài $72^\circ\text{C}$ (1 phút)]; Kéo dài kết thúc $72^\circ\text{C}$ (5 phút); Giữ lạnh $4^\circ\text{C}$.

# BioPython script tự động xử lý chuỗi và phát hiện biến dị SNP từ file FASTA
from Bio import SeqIO
from Bio.Align import MultipleSeqAlignment

def analyze_snps(alignment_file, reference_id):
    alignment = SeqIO.parse(alignment_file, "fasta")
    records = {record.id: str(record.seq).upper() for record in alignment}
    
    ref_seq = records[reference_id]
    snps_found = {}
    
    for sample_id, seq in records.items():
        if sample_id == reference_id:
            continue
        diffs = []
        for pos, (ref_base, sample_base) in enumerate(zip(ref_seq, seq), start=1):
            if sample_base != ref_base and sample_base in ['A', 'C', 'G', 'T']:
                diffs.append((pos, ref_base, sample_base))
        snps_found[sample_id] = diffs
        
    return snps_found

# Ví dụ xuất kết quả phân tích SNP
# snp_results = analyze_snps("aligned_rachgia_dogs.fasta", "Canis_lupus_Ref")

Testing và validation

  • Chất lượng DNA tổng số: Điện di trên gel agarose 1% (đệm TBE 1X, 90V) cho thấy toàn bộ 23 mẫu xuất hiện băng gọn, đậm nét trên thang 10 kb, không có hiện tượng vệt mờ (smear) do đứt gãy DNA. Nồng độ quang phổ $OD_{260}/OD_{280}$ đạt ngưỡng tiêu chuẩn $1{,}62 - 2{,}10$ (mẫu R5 đạt cao nhất $181{,}3,\text{ng}/\mu\text{L}$, tỷ lệ 2,00).
  • Độ đặc hiệu sản phẩm PCR: Điện di sản phẩm PCR cho một băng đơn duy nhất có kích thước ~1.232 bp, đối chứng âm không xuất hiện băng lạ hay primer dimer.
  • Xác thực Sanger: 100% chromatogram được hiệu chỉnh trên FinchTV loại bỏ các đỉnh nhiễu, đối chiếu mạch thuận - nghịch tạo chuỗi đồng nhất chính xác tuyệt đối trên đoạn lõi 582 bp vùng siêu biến HV1.

Kết quả đạt được

Phân tích 582 bp vùng HV1 trên 23 cá thể chó nhà Rạch Giá xác định được 9 Haplotype phân bố trên 4 Haplogroup (A, B, C, E):

Nhóm Haplogroup Tần số xuất hiện (%) Các Haplotype xác định Số cá thể mang kiểu gen Tỷ lệ quần thể toàn cầu (Pang et al. 2009)
Nhóm A 60,0% A11, A17, A18, A73 14 cá thể 71,3% (Phổ biến toàn cầu)
Nhóm B 10,0% B1, B3 2 cá thể ~15,0% (Phổ biến toàn cầu)
Nhóm C 20,0% C2, C3 5 cá thể ~9,0% (Châu Á, Âu, Phi)
Nhóm E 10,0% E4 2 cá thể < 2,0% (Dòng cực hiếm Đông Á) 
Tần số phân bố Haplogroup tại Quần thể Chó nhà Rạch Giá (n=23):
██████████████████████████████ 60% Nhóm A (A11, A17, A18, A73)
█████ 10% Nhóm B (B1, B3)
██████████ 20% Nhóm C (C2, C3)
█████ 10% Nhóm E (E4 - Dòng gen cổ Đông Á quý hiếm)

Kết quả phân tích khoảng cách di truyền trên phần mềm MEGA 7.0 cho thấy các cá thể thuộc Haplotype E4 tại Rạch Giá có mức độ đồng nhất nucleotide tuyệt đối (100% identity) với dòng E4 phân lập trên chó Phú Quốc (Thái Kế Quân et al., 2016) và có quan hệ gần gũi với chó Pungsan (Triều Tiên) và các giống chó cổ Nhật Bản.

Đổi mới và đóng góp

  1. Phát hiện dòng Haplotype E4 quý hiếm trên đất liền: Haplogroup E là nhóm haplotype có tần số cực thấp trên thế giới (<2%), phân bố giới hạn tại Đông Á (Nhật Bản, Hàn Quốc, miền Nam Trung Quốc). Trước đây, tại Việt Nam, dòng E (E1, E4) chỉ được ghi nhận với tần số cao trên quần thể chó xoáy đảo Phú Quốc (16,7%). Việc ghi nhận Haplotype E4 chiếm tới 10% tại Rạch Giá là minh chứng di truyền học phân tử đầu tiên khẳng định dòng gen cổ này không chỉ tồn tại cô lập trên đảo mà đã có sự phân lưu, tiếp xúc lâu đời giữa đất liền Kiên Giang và đảo Phú Quốc.
  2. Cải tiến kỹ thuật tách chiết DNA từ gốc lông động vật: Tối ưu hóa thành công giao thức ly giải SDS-Proteinase K hai bước kết hợp tủa muối nồng độ cao (NaCl 5M), nâng hiệu suất thu hồi DNA đạt $57{,}0 - 181{,}3,\text{ng}/\mu\text{L}$, cho phép sử dụng mẫu gốc lông không xâm lấn thay thế hoàn toàn việc lấy máu hay mô xâm lấn gây đau đớn cho vật nuôi.
  3. Tối ưu chu trình nhiệt phân dải cao cho đoạn khuếch đại 1,2 kb: Khắc phục triệt để khó khăn khi thiết kế mồi do vùng lặp VNTR liền kề bằng việc khuếch đại trọn vẹn đoạn lớn 1.232 bp với 3 dải nhiệt độ bắt cặp ($49^\circ\text{C}$, $54^\circ\text{C}$, $55^\circ\text{C}$), đạt tỷ lệ khuếch đại thành công 100%.
Tiêu chí đối sánh Đinh Trần Mỹ Đức (2013) Trương N. T. Như Mai (2013) Nguyễn Văn Tú (2015) Nghiên cứu này (Lê Quế Phương, 2017)
Đối tượng nghiên cứu Chó Phú Quốc & Chó thường Chó Phú Quốc & Chó thường Chó Phú Quốc (tại đảo) Chó nhà Rạch Giá (đất liền)
Chỉ thị di truyền Vùng gen 16S rRNA 582 bp vùng HV1 mtDNA 582 bp vùng HV1 mtDNA 582 bp HV1 (thuộc đoạn 1.232 bp CR)
Kích thước mẫu ($n$) 19 mẫu 25 mẫu 30 mẫu 23 mẫu
Số SNP phát hiện 15 SNPs 16 SNPs 28 SNPs Phát hiện đầy đủ các SNP đặc trưng
Số Haplotype ghi nhận Không phân loại được 5 haplotypes 11 haplotypes 9 haplotypes (A, B, C, E)
Tỷ lệ dòng hiếm nhóm E Không xác định Có xuất hiện (E1, E4) 16,67% (E1, E4) 10,00% (Haplotype E4)

Ứng dụng thực tế và triển khai

[TRẠM BẢO TỒN / GIỐNG]                [PHÒNG LAB DI TRUYỀN]                [CƠ SỞ DỮ LIỆU GEN]
  • Bảo tồn và nhân giống có định hướng: Cung cấp công cụ chỉ thị phân tử giúp các hiệp hội bảo tồn (như VKA - Vietnam Kennel Association) và trại giống kiểm định nguồn gốc dòng mẹ, hạn chế giao phối cận huyết, chọn lọc cá thể mang haplotype quý (nhóm E) để nhân thuần.
  • Giám định pháp y thú y: Trình tự vùng HV1 mtDNA là tiêu chuẩn vàng trong pháp y sinh học nhằm xác định vết tích lông động vật tại hiện trường vụ việc, giải quyết tranh chấp quyền sở hữu chó cảnh bản địa giá trị cao.
  • Lộ trình mở rộng cơ sở dữ liệu:
    • Giai đoạn 1 (0–6 tháng): Số hóa 23 trình tự lên GenBank, tích hợp ma trận SNP vào cơ sở dữ liệu di truyền chó bản địa Việt Nam.
    • Giai đoạn 2 (6–18 tháng): Mở rộng phạm vi lấy mẫu sang các tỉnh lân cận thuộc Đồng bằng sông Cửu Long (Cà Mau, Bạc Liêu, Cần Thơ) với quy mô $n \ge 200$.
    • Giai đoạn 3 (18–36 tháng): Kết hợp giải trình tự đoạn gen nhân và microsatellite để xây dựng chip mã vạch DNA (DNA Barcoding) nhận diện nhanh giống chó bản địa.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật:
    • mtDNA chỉ phản ánh lịch sử tiến hóa dòng mẹ, không ghi nhận các biến cố lai tạo chéo từ dòng bố.
    • Kích thước mẫu ($n=23$) dù đủ ý nghĩa thống kê cho khảo sát ban đầu tại Rạch Giá nhưng cần mở rộng quy mô để đánh giá toàn diện tần số alen hiếm.
    • Phương pháp giải trình tự Sanger có chi phí tính trên mẫu đơn lẻ cao hơn so với Next-Generation Sequencing (NGS) khi triển khai quy mô hàng nghìn cá thể.
  • Hướng phát triển đề xuất:
    • Bổ sung chỉ thị trên nhiễm sắc thể giới tính Y (vùng MSY - Male Specific Region of Y chromosome) để làm sáng tỏ dòng tiến hóa theo dòng cha (patrilineal lineage).
    • Ứng dụng công nghệ Multiplex PCR hoặc giải trình tự gen thế hệ mới (Targeted NGS Panel) để phân tích đồng thời vùng HV1, HV2 và các locus microsatellite trong cùng một phản ứng.

Đối tượng hưởng lợi

  • Sinh viên & Học viên chuyên ngành Công nghệ Sinh học: Nắm vững giao thức thực nghiệm tách chiết DNA từ mẫu khó (gốc lông), kỹ thuật PCR tối ưu hóa gradient nhiệt và thành thạo bộ công cụ tin sinh học (FinchTV, SeaView, MEGA).
  • Các nhà nghiên cứu Di truyền & Động vật học: Bổ sung dữ liệu đối sánh phân tử quan trọng về nguồn gen chó đất liền Tây Nam Bộ, đóng góp vào bản đồ đa dạng sinh học của chi Canis tại Đông Nam Á.
  • Các nhà bảo tồn & Hiệp hội nhân giống vật nuôi: Sở hữu cơ sở khoa học chính xác để xây dựng tiêu chuẩn đăng bạ giống, truy nguyên huyết thống và bảo tồn nguồn gen "Tứ đại danh khuyển" Việt Nam.
  • Cơ quan kiểm định thú y & Pháp y sinh học: Ứng dụng quy trình phân tích ADN không xâm lấn vào công tác giám định, kiểm dịch động vật và bảo vệ nguồn gen bản địa.

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

1. Tại sao nghiên cứu lại chọn mẫu lông thay vì mẫu máu để tách chiết DNA?

Mẫu lông (đặc biệt là phần gốc chân lông chứa các tế bào nang lông) là phương pháp thu thập mẫu sinh học hoàn toàn không xâm lấn, không gây đau đớn hay tổn thương cho vật nuôi, không cần chuyên viên thú y thao tác lấy ven máu. Mẫu lông dễ dàng bảo quản khô trong túi zipper ở điều kiện thường khi di chuyển đường dài trước khi cấp đông $-20^\circ\text{C}$, giảm thiểu nguy cơ lây nhiễm mầm bệnh qua đường máu.

2. Ý nghĩa sinh học của việc phát hiện Haplotype E4 tại Rạch Giá là gì?

Haplogroup E là một nhánh tiến hóa rất cổ xưa và hiếm gặp trên thế giới (<2% trên quy mô toàn cầu, chỉ xuất hiện ở một số giống chó cổ tại Nhật Bản, Hàn Quốc). Tại Việt Nam, trước đây chỉ có chó Phú Quốc trên đảo được ghi nhận mang dòng E với tỷ lệ cao. Việc phát hiện E4 chiếm 10% tại Rạch Giá chứng minh rằng dòng gen cổ quý hiếm này đã hiện diện trong quần thể chó đất liền Tây Nam Bộ, phản ánh lịch sử giao lưu nguồn gen lâu đời giữa bán đảo Kiên Giang và đảo Phú Quốc.

3. Vùng siêu biến HV1 có ưu điểm gì hơn so với các vùng gen ty thể khác như CytB hay 16S rRNA?

Vùng 16S rRNA và Cytochrome b (CytB) là các vùng gen mã hóa chức năng có tốc độ tiến hóa tương đối chậm, chỉ thích hợp để định danh phân loại ở cấp độ loài hoặc chi. Trong khi đó, HV1 nằm trong vùng kiểm soát CR (Control Region) không mã hóa, chịu áp lực chọn lọc thấp hơn nên tích lũy đột biến điểm (SNP) nhanh hơn 10–25 lần so với DNA nhân, tạo ra độ phân giải rất cao đủ để phân biệt các giống, nòi và dòng haplotype dưới cấp độ loài.

4. Cần xử lý như thế nào khi đồ thị huỳnh quang Sanger xuất hiện tín hiệu nhiễu ở hai đầu đoạn đọc?

Tín hiệu nhiễu ở 20–40 base đầu tiên thường do lượng mồi dư thừa hoặc enzyme gắn kết không ổn định; tín hiệu suy giảm ở cuối đoạn đọc (>800 bp) do cạn kiệt ddNTP huỳnh quang. Quy trình chuẩn là sử dụng phần mềm FinchTV để tỉa bỏ (trim) các vùng có chỉ số chất lượng thấp (Quality Score < 20), sau đó tiến hành gióng hàng 2 chiều mạch thuận (15412F) và mạch nghịch (16625R) trên SeaView để tạo chuỗi đồng thuận (Consensus) chính xác ở đoạn giữa 582 bp.

5. Chi phí và thời gian triển khai một quy trình định danh Haplotype như trong nghiên cứu?

Chi phí hóa chất tách chiết và phản ứng PCR trong phòng thí nghiệm dao động khoảng 50.000 – 70.000 VNĐ/mẫu; chi phí gửi dịch vụ giải trình tự Sanger 2 chiều khoảng 150.000 – 200.000 VNĐ/mẫu. Tổng chi phí phân tích phân tử dao động dưới 300.000 VNĐ/mẫu với thời gian hoàn thành từ 48 đến 72 giờ làm việc, hoàn toàn khả thi để ứng dụng rộng rãi trong các chương trình kiểm định di truyền thực tế.

Kết luận

Nghiên cứu của tác giả Lê Quế Phương đã đánh giá thành công sự đa dạng di truyền của quần thể chó nhà tại TP. Rạch Giá (Kiên Giang) thông qua phân tích vùng siêu biến HV1 trên hệ gen ty thể. Công trình đã chuẩn hóa quy trình tách chiết DNA hiệu suất cao từ gốc lông, tối ưu hóa dải nhiệt độ khuếch đại PCR phân đoạn 1.232 bp và giải mã trọn vẹn đoạn lõi 582 bp trên 23 cá thể.

Việc phát hiện 9 Haplotype thuộc 4 nhóm (A: 60%, B: 10%, C: 20%, E: 10%), đặc biệt là sự hiện diện của Haplotype E4 quý hiếm, đã cung cấp luận cứ khoa học vững chắc khẳng định mức độ đa dạng di truyền cao của chó nhà Tây Nam Bộ và mối liên hệ phát sinh loài mật thiết với chó xoáy Phú Quốc. Kết quả nghiên cứu đóng góp cơ sở dữ liệu quan trọng phục vụ công tác bảo tồn nguồn gen vật nuôi bản địa, chọn giống thú cảnh và phát triển công nghệ giám định sinh học phân tử tại Việt Nam.