Chương 1. TỔNG QUAN NGHIÊN CỨU 1. Phương pháp phát sinh loài Phân tích phát sinh loài là thiết lập mối quan hệ tiến hóa giữa các gen hoặc đặc điểm của sinh vật. Và mặc dù sự sống có một nguồn gốc duy nhất và chỉ có một lịch sử tiến hóa thực sự, nhiệm vụ tái tạo lại lịch sử này có thể từ dễ dàng đến khó khăn tùy thuộc vào các đơn vị phân loại cụ thể hoặc các sinh vật được phân tích.
Những khó khăn phát sinh với việc tái tạo phát sinh loài do thực tế là hầu như tất cả các suy luận về sự tiến hóa và các mối quan hệ đều dựa trên nghiên cứu về các sinh vật hiện tại. Hiếm khi có thể lấp đầy những thiếu sót để xác minh thông tin về sự sống thực sự diễn ra như thế nào. Những vấn đề này được phóng đại ở vi khuẩn vì chúng có ít đặc điểm có thể nhìn thấy và hầu như không có hồ sơ hóa thạch. Hơn nữa, vi khuẩn thì cổ xưa và đã đa dạng hóa qua hàng tỷ năm và là đối tượng của các quá trình chuyển gen ngang (Horizontal gene transfer - HGT), có khả năng tạo ra bất kỳ đặc điểm nào trong bất kỳ dòng dõi nào và che khuất các mô hình tổ tiên truyền thống của Darwin, do đó làm cho lịch sử của chúng thậm chí còn khó xác định hơn [13].
Trong hàng trăm năm đầu tiên của ngành vi sinh học, vi khuẩn được phân loại, phân biệt và nhóm lại chủ yếu theo môi trường sống, đặc điểm sinh trưởng, thuộc tính sinh hóa và tiềm năng độc lực của chúng. Nhưng trong nửa sau của thế kỷ 20, đã có một sự thay đổi rõ rệt đối với việc sử dụng thông tin di truyền phân tử, thu được thông qua phân tích trình tự protein và axit nucleic. Điều này cung cấp ba lợi thế ngay lập tức. - Đầu tiên, các phân tích so sánh trở nên tập trung vào các ký tự được phân phối phổ biến, giúp loại bỏ các dòng được phân loại cùng nhau do thiếu các đặc điểm xác định cụ thể của chúng.
- Thứ hai, kiến thức chuyên sâu về các cơ chế vật lý mà theo đó các đại phân tử thông tin có thể thay đổi đã dẫn đến sự phát triển của các mô hình mạnh mẽ về tiến hóa trình tự. - Thứ ba, việc sử dụng các trình tự phân tử đã làm tăng đáng kể số lượng các ký tự rời rạc, theo đó các vi sinh vật có thể được so sánh với mỗi nucleotide hoặc axit amin có thể đóng vai trò là một đặc điểm thông tin [13].1 Phân tích phát sinh loài dựa vào trình tự phân tử Việc sử dụng dữ liệu trình tự phân tử để suy ra các mối quan hệ là trọng tâm về quá trình tiến hóa. Những thay đổi ở cấp độ phân tử đã được sử dụng để thiết lập phát sinh loài của nhiều loại sinh vật. Đối với vi khuẩn, đặc biệt, nơi các dấu hiệu kiểu hình cổ điển thường không hiệu quả hoặc không thực tế để nhóm các phân lập vi khuẩn, các phương pháp phát sinh loài phân tử đã mang tính cách mạng [13].
5 - Sắp xếp trình tự Điều kiện tiên quyết để tạo ra một kiểu phát sinh loài của sinh vật dựa trên trình tự phân tử là sự sắp xếp các trình tự trực giao (di truyền theo chiều dọc) gần đúng nhất với quá trình tiến hóa phân tử thực sự của chúng. Sự sắp xếp có thể được tạo bằng nhiều thuật toán, hầu hết trong số đó sử dụng một số dạng ma trận thay thế để định vị các nucleotide tương đồng hoặc axit amin để tối đa hóa số lượng vị trí giống hệt hoặc tương tự tại một vị trí nhất định. Clustal, hiện là phần mềm căn chỉnh phổ biến nhất, sử dụng ma trận IUB và ClustalW1.6 để sắp xếp DNA và ma trận PAM, BLOSUM và Gonnet để sắp xếp axit amin [13]. Khi một tập dữ liệu trình tự được đưa vào chương trình căn chỉnh, chẳng hạn như Clustal, thuật toán bắt đầu bằng cách tạo sự căn chỉnh theo cặp của tất cả các trình tự một cách độc lập.
Sau đó, thuật toán sẽ tính toán tất cả khoảng cách giữa các cặp trình tự, tạo ra một ma trận khoảng cách mà sau đó được chuyển thành biểu đồ cây (dendrogram). Như tên gọi của nó, cây này phục vụ như một hướng dẫn để bắt đầu căn chỉnh nhiều trình tự, với các trình tự giống nhau nhất được thêm vào trước, tiếp theo là bổ sung dần dần các trình tự khác nhau hơn. Loại phương pháp căn chỉnh trình tự lũy tiến theo cặp đơn giản và tương đối hiệu quả này đã được sử dụng rộng rãi trong nhiều nghiên cứu phát sinh gen, nhưng kém hiệu quả hơn đối với các bộ dữ liệu cụ thể. Điều này đã thúc đẩy sự phát triển của các thuật toán căn chỉnh cải tiến hơn, chẳng hạn như MUSCLE, T-Coffee, POA, DIALIGN, SAGA và MAFFT, hiệu quả và chính xác trong việc khôi phục các căn chỉnh tối ưu từ các bộ dữ liệu khác nhau, được thiết lập bằng cách so sánh với BAliBASE (cơ sở dữ liệu sắp xếp tham chiếu được tinh chỉnh thủ công).
Tuy nhiên, không phải tất cả các thuật toán này đều tương đương về thuật toán, chẳng hạn như MUSCLE, T-Coffee và MAFFT, thực hiện một loạt các phép toán lũy tiến, căn chỉnh và tối ưu hóa trong mỗi lần lặp lại để đảm bảo độ chính xác căn chỉnh cao hơn đáng kể. Sự liên kết tiến bộ, tinh tế được cung cấp bởi các phương pháp này là điều cần thiết để tái tạo chính xác các mối quan hệ phát sinh gen. Bất kể chương trình hoặc thuật toán được sử dụng để tạo liên kết nhiều trình tự, mỗi liên kết phải được đánh giá trước khi sử dụng để tái tạo phát sinh gen. Biến thể trình tự trên một căn chỉnh có thể không đồng nhất, dẫn đến các vùng được bảo tồn nằm xen kẽ giữa các vùng biến đổi cao.
Do sự liên kết của các vùng có thể siêu biến thường là vấn đề nên các phần này được che giấu hoặc loại bỏ tốt nhất. Điều này đặc biệt quan trọng bởi vì ngay cả những thay đổi nhỏ trong sự liên kết cũng sẽ bị nhiễu thông tin có khả năng che khuất sự phát sinh loài thực sự [13].2 Xây dựng cây phát sinh loài Khi một sự liên kết đáng tin cậy được tạo ra, nhiều phương pháp xây dựng cây có thể được sử dụng để chuyển đổi dữ liệu liên kết thành một cây phát sinh loài. Các 6 phương pháp xây dựng cây này được phân loại rộng rãi thành các phương pháp tiếp cận dựa trên khoảng cách, phân tích và xác suất [13].1 Cây khoảng cách Các cách tiếp cận dựa trên khoảng cách phân loại đại diện cụm dựa trên số lượng thay thế nucleotide hoặc axit amin giữa các trình tự. Một trong những phương pháp dựa trên khoảng cách đầu tiên và đơn giản nhất, được phát triển để tạo ra các bản sao kiểu hình, là UPGMA (Phương pháp nhóm cặp không trọng số với trung bình số học).
UPGMA sử dụng phương pháp phân cụm liên tiếp, theo đó một ma trận của tất cả các điểm tương đồng của trình tự theo cặp được tạo ra và hai đơn vị phân loại có khoảng cách nhỏ nhất được phân cụm trước. Ma trận khoảng cách được tính toán lại, coi các đơn vị phân loại đã được nhóm thành một đơn vị phân loại và đơn vị phân loại có khoảng cách nhỏ nhất tiếp theo được thêm vào cây. Điều này được lặp lại cho đến khi tất cả các đơn vị phân loại được thêm vào và một bộ phát sinh loài cuối cùng được hình thành. Do cách giải thích đơn giản này về ma trận khoảng cách và giả định về tốc độ tiến hóa không đổi (đồng hồ phân tử) cho tất cả các trình tự, UPGMA có xu hướng xây dựng các cây phát sinh loài không được hỗ trợ (unsupported phylograms).
Phương pháp Neighbor-Joining (NJ)[14], có cách tiếp cận tương tự để phân nhóm theo phân loại như UPGMA, trước tiên phân cụm các chuỗi (láng giềng) gần nhất và sau đó tính toán lại khoảng cách giữa các cặp lân cận. Do đó, NJ không gặp phải những hạn chế giống như UPGMA và đã trở thành phương pháp dựa trên khoảng cách được ưa thích để tái tạo phát sinh gen[13].2 Cây Likelihood Việc áp dụng các phương pháp likelihood để tái cấu trúc phát sinh loài đã trở nên ngày càng phổ biến, phần lớn là do độ chính xác và tính nhất quán cao hơn một chút của chúng trong việc khôi phục một phát sinh loài chính xác, và do sự gia tăng đáng kể về khả năng và tốc độ tính toán. Maximum Likelihood (ML) [15] và phương pháp tiếp cận Bayesian [16] đưa ra hai cách tiếp cận xác suất riêng biệt nhưng có liên quan để xác định kiểu phát sinh loài tốt nhất. Các phương pháp ML cố gắng xác định cấu trúc liên kết cây có xác suất cao nhất với dữ liệu chuỗi được cung cấp.
ML cần nhiều tính toán và không phù hợp với các tập dữ liệu rất lớn [13]. Phương pháp Bayesian để tái cấu trúc phát sinh loài được xem là một giải pháp thay thế lý tưởng cho phương pháp ML. Không giống như cách tiếp cận dựa trên khoảng cách, cách tiếp cận Bayes không tạo ra một cây duy nhất, nó lấy mẫu một loạt các cấu trúc liên kết cây có thể xảy ra dựa trên bộ dữ liệu và các phân bố xác suất trước (Prior probability). Bởi vì các xác suất sau đã được tính toán trong quá trình lấy mẫu, chúng được sử dụng để xác định giá trị độ tin cậy tại mỗi nút trên cây, không cần hỗ trợ thống kê thêm các giá trị như bootstrap hoặc jackknife [13].3 Cây Parsimony Các phương pháp phân tích để xây dựng cây phát sinh loài được đặt tiền đề để ủng hộ con đường tiến hóa với số lượng thay đổi ít nhất.
Không giống như các phương pháp dựa trên khoảng cách, tính phân tích dựa trên suy luận dựa trên ký tự của các kiểu phát sinh loài, sử dụng dữ liệu trình tự ở trạng thái ban đầu thay vì chuyển đổi nó thành khoảng cách. Một trong những phương pháp phân tích phổ biến hơn là Maximum Parsimony (MP)[17], một phương pháp tính toán các bước tiến hóa cho tất cả các cấu trúc liên kết có khả năng và trình bày cây có các bước tiến hóa ít nhất. Thông thường, có nhiều cây chi tiêu ngang nhau, một số có cấu trúc liên kết khác hẳn so với những cây khác.