2.2.1.2. Đánh số nucleotide và trình bày các đột biến
Mỗi trình tự nucleotide được lưu trữ trong cơ sở dữ liệu sẽ được sắp gióng cột với trình tự chuẩn (mã số GenBank U96639.2) bằng phần mềm ClustalW [45]. Kết quả sắp gióng cột sẽ được hiệu chỉnh thủ cơng để xếp cho đúng vị trí của một vài nucleotide đột biến thêm vào hay mất đi (indel) theo đề xuất của Pereira và cộng sự [57]. Do trong GenBank, các trình tự được đệ trình có nhiều kích thước khác nhau, nên để thống nhất, các nucleotide trong mỗi trình tự sẽ được đánh số tương ứng với vị trí của nucleotide trong trình tự chuẩn [57], nghĩa là lấy vị trí tuyệt đối của nucleotide trên tồn bộ DNA ty thể, khơng phải vị trí tương đối của chúng trong trình tự. Các đột biến so với trình tự chuẩn được trình bày theo định dạng XaY, với X là nucleotide trên trình tự chuẩn ở vị trí a và Y là nucleotide tại vị trí tương ứng trên trình tự khảo sát. Vị trí nucleotide bị trống (trong trường hợp đột biến thêm hay mất nucleotide) được biểu thị bằng dấu “-“. Ví dụ, T15639G biểu thị nucleotide T ở vị trí 15639 đã bị biến đổi thành nucleotide G, -15535.1C -15535.2C biểu thị hai nucleotide C được thêm vào sau nucleotide ở vị trí 15535. Tất cả các đột biến xuất hiện trên vùng 582 cặp base của mỗi trình tự được gọi là bộ đột biến của trình tự, được dùng để xác định haplotype của trình tự đó.
2.2.1.3. Xác định haplotype và hiệu chỉnh dữ liệu
Các trình tự và các haplotype đã được cơng bố sẽ được thu thập để phân tích, xác định các bộ đột biến của mỗi haplotype. Các bộ đột biến này sẽ được dùng để làm cơ sở để xác định haplotype của các trình tự chưa biết haplotype. Ví dụ, haplotype A9 có bộ đột biến là “A15627G T15639A C15814T”, nghĩa là so với trình tự chuẩn thì trình tự có haplotype A9 mang 3 đột biến A thành G, T thành A và C thành T ở 3 vị trí tương ứng là 15627, 15639 và 15814.
Những thông tin về haplotype của các trình tự DNA bị sai lệch hoặc khơng rõ ràng sẽ được phát hiện và điều chỉnh. Quá trình này sẽ được tiến hành với các luật sau:
1. Trình tự đầu tiên được dùng để công bố một haplotype sẽ được ghi nhận là trình tự chuẩn của haplotype đó.
2. Trình tự giống trình tự chuẩn của một haplotype đã công bố nhưng được dùng để công bố haplotype khác thì haplotype sau sẽ bị loại ra khỏi danh sách haplotype đã được công bố.
3. Trình tự khác với trình tự chuẩn của một haplotype đã công bố nhưng cơng bố tên giống với haplotype đó thì được xem là cơng bố sai, thơng tin cần phải được hiệu chỉnh.
Tồn bộ các trình tự có haplotype đã biết sẽ được sắp gióng cột để xác định các đột biến đặc trưng cho từng haplogroup A, B, C, D, E, F. Các trình tự có bộ đột biến mới, chưa xác định được haplotype, nếu mang những đột biến đặc trưng của haplogroup nào sẽ được ghi nhận là có haplotype mới thuộc haplogroup đó.
2.2.1.4. Cơng cụ xác định haplotype vùng HV1 của DNA ty thể
Công cụ xác định nhanh haplotype vùng HV1 của DNA ty thể ở chó được viết bằng ngơn ngữ lập trình Perl, với giao diện là một trang web, được đặt tên là “Haplotype identifier”. Bằng cách so sánh trình tự cần xác định haplotype (trình tự truy vấn) với trình tự chuẩn [37], Haplotype identifier nhận diện đoạn trình tự 582 cặp base cần phân tích. Nếu trình tự truy vấn khơng bao gồm hoặc chỉ là một phần của đoạn trình tự 582 cặp base, chương trình sẽ báo trình tự khơng thích hợp để xác định haplotype. Nếu trình tự truy vấn bao gồm tồn bộ đoạn trình tự 582 cặp base, chương trình sẽ xác định các đột biến của trình tự truy vấn so với trình tự chuẩn. Các thơng tin về đột biến này sẽ được kiểm tra với các haplotype đã được công bố để định loại haplotype nếu trùng khớp với một haplotype xác định đã biết. Nếu không trùng khớp với một haplotype đã biết, trình tự này được xác định là một haplotype mới. Tập hợp các đột biến của haplotype mới sẽ tiếp tục được so sánh với tập hợp các đột biến đặc trưng của từng haplogroup để xác định haplogroup của haplotype mới này (Hình 2.4).
Hình 2.4: Sơ đồ hoạt động của công cụ xác định haplotype "Haplotype identifier”
2.2.1.5. Xây dựng cây phát sinh chủng loại
Cây phát sinh chủng loại được thiết lập bằng chương trình Paup* [82], sử dụng phương pháp Neighbor-Joining với 2000 lần lặp lại để xác định độ tin cậy của các nhánh cây. Ngồi các trình tự DNA cần khảo sát, việc thiết lập cây cịn bao gồm trình tự của 319 haplotype thuộc 6 haplogroup A, B, C, D, E, F đã được công bố trước đây. Các trình tự có mã số GenBank NC_008093, KF661096, DQ480509 – DQ480511, EU789789 của Canis latrans được sử dụng làm nhóm ngoại (outgroup) để xác định gốc cây. Hình ảnh cây được thể hiện bởi phần mềm TreeExplorer của phần mềm MEGA6 [84].
2.2.2. Xác định trình tự nucleotide đoạn 582 cặp base vùng HV1 DNA ty thể của chó nhà Việt Nam và chó lưng xốy Phú Quốc chó nhà Việt Nam và chó lưng xốy Phú Quốc
2.2.2.1. Dụng cụ và thiết bị sử dụng
Các dụng cụ (Phụ lục 2) và thiết bị (Phụ lục 3) cơ bản được trang bị ở Phịng thí nghiệm Sinh học phân tử (Đại học Sài gòn), Viện bảo tồn gene (Đại học Nguyễn Tất Thành) như máy ly tâm, máy PCR… đã đáp ứng được hầu hết các nội dung cần thực hiện của đề tài. Phân tích dữ liệu và xây dựng công cụ định loại haplotype được thực hiện trên máy tính cá nhân.
2.2.2.2. Thu thập và xử lý mẫu
Đề tài sử dụng hai nhóm chó để khảo sát: chó lưng xốy Phú Quốc và chó nhà Việt Nam. Cụm từ “Chó nhà Việt Nam” sử dụng trong đề tài để chỉ nhóm chó nhà (chó cỏ) Việt nam khơng có các đặc điểm hình thái giống các giống chó ngoại nhập thơng thường như chó đốm, chó Shepherd của Đức,… và khơng bao gồm chó lưng xốy Phú Quốc. Các tiêu chí để nhận diện chó lưng xốy Phú Quốc sử dụng trong đề tài bao gồm:
(1) Có xốy lưng. Đây là đặc điểm khơng độc quyền của chó lưng xốy Phú Quốc (chó lưng xốy Thái Lan cũng có đặc điểm này) cũng như khơng phải tất cả các chó lưng xốy Phú Quốc đều có lưng xốy. Việc sử dụng đặc điểm có xốy lưng nhằm giảm thiểu việc thu mẫu nhầm với các giống chó khơng có lưng xốy khác.
(2) Có nguồn gốc từ đảo Phú Quốc (chó hoặc chó cha mẹ đang sống trên đảo hoặc được mang đi nơi khác từ đảo)
(3) Chó trưởng thành (hơn 1 tuổi), trơng khỏe mạnh, có thể trọng khoảng 15 kg, chiều cao khoảng từ 45 – 55 cm. Chiều cao, cân nặng của chó khơng phải là đặc điểm tiên quyết để phân biệt với chó lưng xốy Thái Lan. Tuy nhiên, việc kết hợp đặc điểm này với các đặc điểm (1), (2) cũng giúp giảm thiểu việc thu mẫu nhầm với chó lưng xốy Thái Lan.
(4) Đầu chó nhỏ, sọ dài, tai nhỏ và đứng. Đi chó cong vút cần câu tạo thành ẵ n ắ vũng trũn, õy l c im dễ thấy ở chó lưng xốy Phú Quốc, đã được ghi nhận bởi Nguyễn Hữu Chiếm và cộng sự, Đại học Cần Thơ.
(5) Các đặc điểm về đốm lưỡi, màng bơi, móng đeo nếu có thì sẽ góp phần khẳng định chó lưng xốy Phú Quốc, nhưng khơng được xem là đặc điểm phải có.
Để dễ dàng trong việc đề cập và so sánh với các giống chó khác, đề tài xem hai nhóm chó khảo sát là hai giống chó riêng biệt, dù hai giống chó này chưa được FCI công nhận là “giống” về mặt phân loại học.
Các nghiên cứu tương tự trên thế giới đối với các giống chó khác nhau thường được thực hiện với số lượng mẫu là 30 con chó hoặc nhiều hơn để có ý nghĩa thống
kê. Căn cứ trên tình hình nghiên cứu chung trên thế giới và nguồn lực của đề tài, 200 mẫu lơng chó bao gồm 100 con chó lưng xốy Phú Quốc và 100 con chó nhà Việt Nam trong nghiên cứu này được thu thập từ đảo Phú Quốc, thành phố Rạch Giá (Kiên Giang), Thành phố Hồ Chí Minh và các vùng phụ cận (Lâm Đồng, Bình Dương) (Hình 2.5; chi tiết ở Phụ lục 4, Phụ lục 5). Những mẫu chó nhà Việt Nam được chọn ở vùng ngoại ô hoặc trong rẫy cách xa khu vực trung tâm để giảm thiểu khả năng lai tạp từ các giống chó ngoại quốc. Những mẫu chó có mối quan hệ gần gũi theo dịng mẹ sẽ bị loại trừ.