Tổng quan luận án

Luận án tiến sĩ "Computational discovery of viruses and their hosts" (Khám phá bằng tính toán về virus và vật chủ của chúng) được thực hiện bởi nghiên cứu sinh Cormac Michael Kinsella (ISBN: 978-94-6483-273-0), bảo vệ tại Đại học Amsterdam (Universiteit van Amsterdam - AMC-UvA, Khoa Y học) vào ngày 11 tháng 9 năm 2023. Công trình được hướng dẫn bởi Promotor là Tiến sĩ Lia van der Hoek (AMC-UvA), cùng hai Copromotor là Giáo sư Ben Berkhout (AMC-UvA) và Tiến sĩ Aldert Bart (Bệnh viện Tergooi). Nghiên cứu nhận tài trợ từ chương trình Horizon 2020 của Liên minh Châu Âu theo thỏa thuận tài trợ Marie Skłodowska-Curie và hỗ trợ xuất bản từ Amsterdam UMC.

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

Trong lịch sử vi sinh học, việc phát hiện virus truyền thống bắt đầu từ việc nhận diện bệnh cảnh lâm sàng trên vật chủ, sau đó phân lập tác nhân trong môi trường nuôi cấy tế bào ("host-first"). Sự phát triển của công nghệ giải trình tự thế hệ mới hiệu năng cao (High-Throughput Sequencing - HTS) kết hợp phương pháp siêu hệ gen (Metagenomics) đã đảo ngược quy trình này: phần lớn các virus mới được phát hiện dưới dạng trình tự bộ gen thô từ các mẫu môi trường hoặc mẫu lâm sàng hỗn hợp mà không xác định được vật chủ ("stray viruses").

Thực trạng này tạo ra hai khoảng trống nghiên cứu chính:

  1. Khoảng trống phương pháp luận phát hiện virus phân kỳ: Các thuật toán dò tìm độ tương đồng trình tự dựa trên căn chỉnh dữ liệu (như BLAST, UBLAST) hoặc k-mer (như Kraken) gặp giới hạn khi xử lý các trình tự phân kỳ sâu ("vùng tối trình tự virus" - virus sequence dark matter), đồng thời các quy trình phân tích HTS chưa tận dụng được đặc tính sinh học phân tử của các thư viện làm giàu virus.
  2. Khoảng trống xác định vật chủ cho virus siêu hệ gen: Phương pháp phân lập nuôi cấy truyền thống có thông lượng thấp và không khả thi với các sinh vật chưa thể nuôi cấy; các mô hình máy học (machine learning) dự đoán vật chủ phụ thuộc nặng nề vào tập dữ liệu huấn luyện đã biết; trong khi đó, ngành virus DNA mạch đơn dạng vòng mã hóa protein sao chép Rep (Cressdnaviricota) chứa hàng loạt dòng virus liên quan đến bệnh tật ở đường tiêu hóa người và động vật nhưng chưa xác định được vật chủ sinh học cụ thể.

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

  1. Xây dựng và kiểm định quy trình tin sinh học phân tích dữ liệu giải trình tự VIDISCA-NGS nhằm tối ưu hóa độ nhạy phát hiện virus đã biết và khám phá các virus phân kỳ cao không phụ thuộc vào căn chỉnh tương đồng chuỗi.
  2. Ứng dụng quy trình tính toán để nhận diện các dòng CRESS virus chưa phân loại trong mẫu phân người và xác định chính xác vật chủ ký sinh trùng của chúng.
  3. Thiết lập khung dự đoán vật chủ không cần dữ liệu huấn luyện và không yêu cầu bộ gen vật chủ lắp ráp hoàn chỉnh cho các cressdnavirus đường tiêu hóa liên quan đến bệnh cảnh lâm sàng.
  4. Điều tra hiện tượng chuyển gen ngang (Horizontal Gene Transfer - HGT) giữa các dòng virus khác nhau để xác định tính hướng vật chủ (tropism) đối với động vật có xương sống.
  5. Khảo sát virome RNA trên các chủng phân lập lâm sàng của các loài vi nấm gây bệnh ở người.

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

  • Đối tượng nghiên cứu: Các virus DNA mạch đơn dạng vòng (CRESS DNA viruses thuộc ngành Cressdnaviricota), các virus RNA của vi nấm (mycoviruses), thực khuẩn thể gokushovirus, cùng các vật chủ tiềm năng bao gồm ký sinh trùng đơn bào (Entamoeba, Giardia, Blastocystis, Dientamoeba, Endolimax), vi nấm gây bệnh ở người và động vật có xương sống.
  • Phạm vi dữ liệu và không gian:
    • Mẫu phân người từ hai đoàn hệ độc lập tại Amsterdam (Tổng $N = 374$): Đoàn hệ 1 gồm 194 mẫu thu thập năm 1994–1995 từ người trưởng thành nhiễm HIV-1 bị tiêu chảy chưa rõ nguyên nhân; Đoàn hệ 2 gồm 180 mẫu từ nghiên cứu đoàn hệ ACS (85 mẫu HIV-1 dương tính và 95 mẫu HIV-1 âm tính).
    • 29 mẫu huyết thanh lợn (13 lợn nái và 16 lợn con) nhiễm virus gây dịch tả heo không điển hình (Atypical Porcine Pestivirus - APPV).
    • Dữ liệu siêu hệ gen công khai toàn cầu từ GenBank, ENA (mã PRJEB35571) và các bộ dữ liệu metagenome từ Bắc Mỹ, Nam Mỹ, Châu Âu, Châu Phi, Châu Á.

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

Luận án điểm lại lịch sử phát triển của các kỹ thuật phát hiện virus:

  • Giai đoạn màng lọc và nuôi cấy: Khái niệm tác nhân gây bệnh siêu lọc được thiết lập qua công trình của Friedrich Loeffler và Paul Frosch (1898) trên bệnh lở mồm long móng, và Martinus Beijerinck (1898) trên bệnh khảm thuốc lá (Tobacco mosaic virus - TMV) với khái niệm "chất lỏng sống truyền nhiễm" (contagium vivum fluidum). Frederick Twort (1915) và Félix d'Hérelle (1917) phát hiện thực khuẩn thể (bacteriophage) và đặt nền móng cho kỹ thuật mảng tan (plaque assay).
  • Giai đoạn sinh học phân tử và giải trình tự: Frederick Sanger công bố bộ gen virus đầu tiên vào năm 1977 ($\phi$X174). Thập niên 1980 chứng kiến sự ra đời của kỹ thuật PCR bởi Kary Mullis và các cộng sự, tiếp nối bằng các phương pháp metagenomics dùng enzyme cắt giới hạn và mồi ngẫu nhiên không phụ thuộc chuỗi như RDA (Representational Difference Analysis), SISPA (Sequence-Independent Single Primer Amplification), và VIDISCA (Virus Discovery based on cDNA-amplified Fragment Length Polymorphism) do nhóm nghiên cứu tại AMC-UvA phát triển.
  • Hệ thống phân loại Cressdnaviricota: Trước thời điểm nghiên cứu (tháng 11/2017), chỉ có 5 họ CRESS virus lây nhiễm sinh vật nhân thực được Ủy ban Quốc tế về Phân loại Virus (ICTV) công nhận: Bacilladnaviridae (nhiễm tảo silic), Circoviridae (nhiễm động vật có xương sống), Geminiviridae (nhiễm thực vật), Genomoviridae (nhiễm nấm), và Nanoviridae (nhiễm thực vật). Sau đó, các họ Smacoviridae, Redondoviridae, Metaxyviridae và các cụm CRESSV1–CRESSV6 được mô tả.

Luận án định vị nghiên cứu vào việc giải quyết nút thắt cốt lõi: kết nối dữ liệu giải trình tự thô với vật chủ sinh học thực sự thông qua các dấu hiệu phân tử tích hợp trong bộ gen (yếu tố virus nội sinh - EVEs, đáp ứng sRNA) và mạng lưới tái tổ hợp mô-đun gen (Rep-Cap recombination networks).


Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết

  • Cấu trúc bộ gen CRESS virus: Bộ gen circular ssDNA nhỏ (1.5–3.5 kb), chứa gốc sao chép có cấu trúc kẹp tóc (stem-loop) với motif nonanucleotide bảo tồn, mã hóa hai protein chủ đạo: protein khởi đầu sao chép (Rep) có miền endonuclease và helicase, và protein vỏ capsid (Cap).
  • Tiến hóa mô-đun và tái tổ hợp: Gen Rep và gen Cap có thể tiến hóa độc lập; hiện tượng hoán đổi mô-đun gen (gene swapping) giữa các dòng virus chỉ có thể xảy ra khi hai virus cùng đồng nhiễm trong một tế bào vật chủ chung.
  • Yếu tố virus nội sinh (EVEs): Quá trình tích hợp tình cờ của bộ gen CRESS virus vào DNA nhiễm sắc thể của tế bào dòng mầm vật chủ tạo thành các hóa thạch phân tử ghi lại lịch sử nhiễm virus cổ xưa.
  • Độ lệch sử dụng codon và hàm lượng GC: Áp lực tiến hóa khiến hàm lượng Guanine-Cytosine (%GC) của virus có xu hướng thích nghi và tương quan thuận với hàm lượng GC của hệ gen vật chủ do giới hạn của bộ máy tRNA.

Phương pháp nghiên cứu

Luận án kết hợp chặt chẽ giữa kỹ thuật thực nghiệm phân tử và các thuật toán phân tích tin sinh học:

[Mẫu phân / Huyết thanh] 
[Chiết tách Boom + Phiên mã ngược với mồi non-rRNA hexamer]
[Khuếch đại thư viện VIDISCA-NGS + Giải trình tự Ion Torrent / Illumina]
[SortMeRNA v2.1: Lọc bỏ rRNA]   [Centrifuge v1.3: Metagenomics]  [CD-HIT v4.7: Gom cụm 100%]
[UBLAST: Khớp DB Protein Virus]  [Recentrifuge: Trực quan HTML] [Cluster-profiling Analysis]
[BLASTn nt: Lọc dương tính giả]                              [Nhận diện Virus copy cao]
  • Quy trình phòng thí nghiệm VIDISCA-NGS:

    • Xử lý mẫu: Ly tâm mẫu phân/huyết thanh loại bỏ cặn tế bào ($5000 \times g$, 10 phút), xử lý dịch nổi bằng $20\ \mu\text{l}$ TURBO DNase ($37^\circ\text{C}$, 30 phút) để phân giải DNA tế bào trần.
    • Tách chiết và phiên mã: Tách axit nucleic bằng phương pháp Boom; phiên mã ngược RNA sử dụng mồi hexamer phi ribosome (non-ribosomal hexamers) để giảm thiểu đọc rRNA của động vật có vú.
    • Xử lý enzyme Mse1: Cắt dsDNA bằng enzyme cắt giới hạn $Mse1$ (nhận biết vị trí TTAA), gắn adapter đặc hiệu vào hai đầu dính, chọn lọc kích thước đoạn chèn từ 200–600 bp, giải trình tự trên hệ máy Ion Torrent PGM / Ion Proton hoặc Illumina MiSeq ($2 \times 150$ bp).
    • Thử nghiệm xác thực vỏ capsid: Lọc dịch nổi qua màng lọc kích thước lỗ $1200\text{ nm}$ và $200\text{ nm}$, ủ với TURBO DNase, sau đó khuếch đại bằng PCR để chứng minh DNA virus được bảo vệ bên trong hạt virion nguyên vẹn.
  • Quy trình tính toán và công cụ tin sinh học:

    • Lọc và tiền xử lý: SortMeRNA v2.1 loại bỏ đọc rRNA; CD-HIT v4.7 gom cụm đọc không phải rRNA ở ngưỡng nhận diện 100% nucleotide.
    • Phân loại virus và lọc dương tính giả: Căn chỉnh UBLAST (USEARCH v10, thiết lập -mincodons 15, -accel 0.8, -evalue 1e-4) vào cơ sở dữ liệu protein virus rút gọn (NCBI Identical Protein Groups gom cụm ở 95% danh tính bằng CD-HIT). Các đọc khớp protein virus được đối chiếu lại cơ sở dữ liệu nucleotide toàn phần NCBI nt bằng BLASTn v2.0 để loại bỏ triệt để các kết quả dương tính giả.
    • Metagenomics toàn phần: Centrifuge v1.3 khớp chỉ mục FM-index nucleotide NCBI (tháng 2/2018), trực quan hóa bằng Recentrifuge và KronaTools v2.7.
    • Phân tích biên dạng cụm (Cluster-profiling analysis): Trích xuất các tệp .clstr từ CD-HIT, tính toán tỷ lệ sao chép của các đoạn cắt Mse1 lặp lại sinh học để phát hiện các bộ gen có số lượng bản sao cao (high copy number) mà không cần căn chỉnh tương đồng.
    • Lắp ráp và phát sinh loài: Lắp ráp de novo bằng SPAdes v3.12 sau khi cắt lọc chất lượng bằng Trimmomatic v0.38. Dựng cây phát sinh loài Maximum-Likelihood (ML) cho protein Rep và Cap với đánh giá bootstrap. Cấu trúc kẹp tóc gốc sao chép được dự đoán bằng MFOLD.

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

Chương 1: General introduction and scope of this thesis

Chương này tổng kết lịch sử ngành virus học, các nguyên lý định nghĩa bản chất virus, sự phát triển của các công cụ phát hiện từ lọc, mảng tan, kính hiển vi điện tử, huyết thanh học, đến sinh học phân tử và HTS. Tác giả phân tích chi tiết "nghịch lý định danh vật chủ" (the host identity problem) khi công nghệ metagenomics tạo ra hàng loạt virus không rõ nguồn gốc. Khái quát đặc điểm ngành Cressdnaviricota và xác định phạm vi 7 chương của luận án.

Chương 2: Enhanced bioinformatic profiling of VIDISCA libraries for virus detection and discovery

Chương 2 giới thiệu quy trình tin sinh học chuyên biệt xử lý dữ liệu VIDISCA-NGS dựa trên cơ chế enzyme cắt giới hạn Mse1. Do bộ gen virus ngắn và có mật độ bản sao sinh học cao trong mẫu nhiễm, enzyme Mse1 tạo ra các đoạn cắt giống hệt nhau về điểm đầu và điểm cuối (biological replicates), dẫn đến hiện tượng đọc trình tự gom cụm cao sau giải trình tự.

Quy trình được kiểm định qua 3 bộ dữ liệu thực tế:

  1. 194 mẫu huyết thanh bệnh nhân HIV-1: Phát hiện HIV-1 trong 128/194 mẫu ($66%$) với số đọc trung bình 42,124 reads/mẫu. Trong số mẫu âm tính, 23/66 mẫu ($35%$) có tải lượng dưới ngưỡng phát hiện của NASBA; 9/128 mẫu dương tính VIDISCA cũng âm tính với NASBA. Tải lượng virus trung vị ở nhóm dương tính là 84 copies/$\mu\text{l}$ và âm tính là 14 copies/$\mu\text{l}$. Tương quan giữa tải lượng virus và tỷ lệ đọc HIV-1 đạt hệ số Spearman $\rho = 0.354$ ($p = 0.001$), mô hình hồi quy tuyến tính $F(1, 192) = 56.9$ giải thích $23%$ phương sai ($R^2 = 0.23$). Kiểm tra 11 mẫu ngoại lai (outliers, tải lượng $> 50\text{ copies}/\mu\text{l}$ nhưng VIDISCA âm tính): 3 mẫu bị áp đảo bởi Hepatitis B virus ($78\text{--}90%$ số đọc), 1 mẫu bị áp đảo bởi Torque Teno virus ($30%$), 6 mẫu chứa $80\text{--}95%$ DNA người/vi khuẩn do ly giải tế bào, 1 mẫu chứa $75%$ đọc không phân loại do nhiễm khuẩn mẫu lưu trữ.
  2. 194 mẫu phân người (phát hiện virus mới qua Cluster-profiling): Phân tích mẫu F115 cho thấy $12%$ trong số 16,160 đọc không phân loại gom cụm vào 100 chuỗi đại diện. Lắp ráp metagenome từ dữ liệu Illumina (ERR233419) tạo ra 9,157 contigs; 8/10 chuỗi đại diện VIDISCA khớp vào một contig duy nhất, phát hiện một gokushovirus mới dài 5 kb (mã truy cập NCBI MK263179) có $72%$ độ đồng nhất nucleotide với họ hàng gần nhất.
  3. 29 mẫu huyết thanh lợn nhiễm APPV: Phát hiện APPV trong 27/29 mẫu. Tương quan mạnh giữa tải lượng virus và tỷ lệ đọc APPV ($R^2 = 0.73$, $F(1,26) = 70.7$). Tải lượng virus tương quan chặt với mức độ gom cụm đọc trung bình (Spearman $\rho = 0.852$, $p < 0.001$). Phân tích kiểm soát sai số PCR cho thấy không có hiện tượng sai lệch khuếch đại template phong phú (Spearman $\rho = 0.169$, $p = 0.380$).

Chương 3: Entamoeba and Giardia parasites implicated as hosts of CRESS viruses

Chương 3 công bố phát hiện 20 trình tự mã hóa CRESS virus hoàn chỉnh (trong đó có 18 bộ gen vòng nguyên vẹn, mã truy cập MT293410–MT293429) từ mẫu phân người của 374 cá thể (thuộc 2 đoàn hệ). Kết hợp giải trình tự metagenomics, PCR nghịch và giải trình tự Sanger, nghiên cứu thiết lập bằng chứng xác định vật chủ ký sinh trùng:

Tình trạng nhiễm ký sinh trùng Tổng số mẫu ($N = 374$) Số mẫu dương tính CRESS virus
Chỉ dương tính với Entamoeba 130 18
Chỉ dương tính với Giardia 3 0
Đồng nhiễm Entamoeba và Giardia 8 3
Âm tính cả Entamoeba và Giardia 233 0
  • Phân tích dịch tễ phân tử: $100%$ ($21/21$) mẫu dương tính với CRESS virus đều nhiễm Entamoeba hoặc Giardia. Mối liên quan có ý nghĩa thống kê cao theo kiểm định Chi bình phương Pearson ($\chi^2 = 36.5$, $p = 1.5 \times 10^{-9}$). Các đơn bào đường ruột khác (Blastocystis, Dientamoeba, Cryptosporidium, Endolimax, Chilomastix, Pentatrichomonas, Retortamonas) đều âm tính ở phần lớn hoặc toàn bộ 21 mẫu này.
  • Yếu tố virus nội sinh (EVEs): Tìm thấy các tích hợp bộ gen CRESS virus hoàn chỉnh trong contig của Entamoeba dispar (contig AANV02000527.1 dài 11.6 kb chứa toàn bộ virus MT293413.1 với độ phủ $100%$, danh tính $84%$) và Giardia duodenalis (contig GL50803_contig0000007.1 dài 6.2 kb chứa virus MT293421.1 với độ phủ $99%$, danh tính $73%$). EVE bảo tồn $100%$ giữa các chủng E. histolytica phân lập tại Mexico năm 1972 (HM-3:IMSS) và Nhật Bản năm 2001 (KU27). Dữ liệu Nanopore (GCA_011634595.1) xác nhận một EVE lặp lại nối tiếp dài 10 kb trong hệ gen G. duodenalis.
  • Đáp ứng RNA can thiệp (RNAi): Dữ liệu giải trình tự sRNA miễn dịch kết tủa với protein AGO2-2 của E. histolytica cho thấy các đỉnh đọc sRNA khớp chính xác vào các vùng EVE và bộ gen CRESS virus ngoại sinh, chứng minh cơ chế ức chế phiên mã kháng virus của vật chủ.
  • Thiết lập 3 họ virus mới:
    • Naryaviridae (nhiễm Entamoeba): gồm 2 chi (Protoentvirus, Deuteroentvirus), khung đọc mở cùng chiều (sense ORFs), hàm lượng GC trung bình $37%$.
    • Nenyaviridae (nhiễm Entamoeba): gồm 2 chi (Tritoentvirus, Tetartoentvirus), lồng trong cụm CRESSV2, khung đọc mở cùng chiều, hàm lượng GC trung bình $42%$.
    • Vilyaviridae (nhiễm Giardia): gồm 3 chi (Protogiavirus, Deuterogiavirus, Tritogiavirus), khung đọc mở ngược chiều/hai chiều (ambisense/antisense), hàm lượng GC trung bình $59%$.
  • Mô hình hồi quy hàm lượng GC: Phân tích 79 trình tự bộ gen virus độc lập cho thấy mối tương quan tuyến tính thuận giữa %GC của CRESS virus và vật chủ ($r^2 = 0.69$, $p = 0.01$). Bộ gen virus ngoại sinh có %GC cao hơn vật chủ trung bình $12.6%$.
  • Mạng lưới tái tổ hợp mô-đun gen: Phân tích bộc lộ sự phân kỳ chiều dài bộ gen ở Naryaviridae và Nenyaviridae do mang hai nhóm protein vỏ Cap có nguồn gốc khác biệt: CCA1 (chiều dài trung bình 467 axit amin) và CCA2 (chiều dài trung bình 180 axit amin) nhưng cùng chia sẻ các gen Rep tương đồng. Vilyaviridae chỉ mang nhóm CCA3. Bằng chứng hoán đổi gen Rep-Cap chứng minh Naryaviridae và Nenyaviridae dùng chung vật chủ Entamoeba, hoàn toàn độc lập với Vilyaviridae.

Chương 4: Host prediction for disease-associated gastrointestinal cressdnaviruses

Mở rộng nghiên cứu phát sinh loài trên các dòng cressdnavirus chưa phân loại, phân giải các cụm từ CRESSV7 đến CRESSV39. Ứng dụng mô hình dự đoán tính toán không cần dữ liệu huấn luyện để xác định vật chủ cho 4 dòng virus liên quan đến bệnh đường tiêu hóa ở người và lợn:

  • Họ Redondoviridae được dự đoán lây nhiễm đơn bào hội sinh/gây bệnh vùng miệng-đường ruột Entamoeba gingivalis.
  • Họ Kirkoviridae lây nhiễm các vi sinh vật thuộc ngành Trùng roi cổ Parabasalids (bao gồm Dientamoeba).
  • Dòng CRESSV1 lây nhiễm Blastocystis.
  • Dòng CRESSV19 lây nhiễm đơn bào Endolimax.

Chương 5: Vertebrate-tropism of a cressdnavirus lineage implicated by poxvirus gene capture

Chương 5 mở rộng phân tích chuyển gen ngang từ virus sang virus (virus-to-virus HGT). Nghiên cứu chỉ ra dòng CRESSV3 đã chuyển gen Rep sang bộ gen của các avipoxvirus (nhóm virus DNA chuỗi đôi kích thước lớn gây bệnh trên chim và bò sát). Hiện tượng thu nhận gen này chứng minh dòng CRESSV3 có tính hướng vật chủ trên động vật có xương sống thuộc nhánh bò sát/chim (saurians), trở thành dòng CRESS virus thứ hai sau Circoviridae được xác định lây nhiễm động vật có xương sống. Nhóm tác giả đề xuất định danh dòng này thành họ mới Draupnirviridae, với bằng chứng tiến hóa cho thấy chúng đã lây nhiễm vật chủ saurian từ hơn 100 triệu năm trước.

Chương 6: Human clinical isolates of pathogenic fungi are host to diverse mycoviruses

Chương 6 tiến hành giải trình tự virome trên các chủng phân lập nấm gây bệnh lâm sàng ở người. Mặc dù không phát hiện thêm CRESS virus mới trên nấm, nghiên cứu đã phát hiện một sự đa dạng lớn các virus RNA (mycovirus) mới cư trú bên trong các mẫu nấm lâm sàng, trong đó có một dòng virus thuộc nhánh chưa từng được xác nhận lây nhiễm vi nấm trước đây.

Chương 7: General discussion

Chương tổng luận tổng hợp các kết quả đạt được, phân tích ý nghĩa của việc mở rộng phân loại học Cressdnaviricota từ 5 lên 8 họ được công nhận, đánh giá vai trò của mycovirus và virus ký sinh trùng đơn bào trong việc làm thay đổi độc lực của mầm bệnh (gây giảm độc lực - hypovirulence hoặc tăng độc lực - hypervirulence), và thảo luận các định hướng ứng dụng tin sinh học trong giải mã virome môi trường và lâm sàng.


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

  1. Quy trình tin sinh học và phân tích Cluster-profiling: Thiết lập quy trình xử lý dữ liệu VIDISCA-NGS tích hợp công cụ Centrifuge, UBLAST và cơ chế lọc dương tính giả hai bước. Sáng kiến phân tích "Cluster-profiling" khai thác sự lặp lại sinh học của enzyme Mse1, cho phép nhận diện các bộ gen virus có mật độ bản sao cao trong mẫu mà không phụ thuộc vào mức độ tương đồng trình tự sẵn có trong cơ sở dữ liệu.
  2. Khám phá và xác lập 3 họ CRESS virus mới: Phân lập và giải mã 20 bộ gen CRESS virus mới, đề xuất thành lập 3 họ virus mới được ICTV chính thức công nhận: Naryaviridae, Nenyaviridae (lây nhiễm Entamoeba) và Vilyaviridae (lây nhiễm Giardia duodenalis). Đây là những virus đầu tiên được chứng minh lây nhiễm trên hai chi ký sinh trùng đơn bào gây ra hơn 300 triệu ca bệnh mỗi năm ở người.
  3. Phương pháp luận mạng lưới tái tổ hợp mô-đun: Đưa ra khung phân tích mạng lưới tái tổ hợp giữa gen Rep và gen Cap (CRESS virus Cap assemblages - CCA1, CCA2, CCA3) để gom cụm các virus có chung phạm vi vật chủ mà không cần tri thức tiên nghiệm về vật chủ.
  4. Xác định vật chủ cho các dòng CRESS virus đường ruột: Dự đoán thành công vật chủ cho các họ và dòng Redondoviridae (Entamoeba gingivalis), Kirkoviridae (Parabasalids/Dientamoeba), CRESSV1 (Blastocystis), và CRESSV19 (Endolimax).
  5. Xác lập họ Draupnirviridae lây nhiễm động vật có xương sống: Chứng minh hiện tượng chuyển gen Rep từ cressdnavirus sang avipoxvirus, xác lập họ Draupnirviridae lây nhiễm động vật nhánh saurian với lịch sử tiến hóa trên 100 triệu năm.
  6. Bổ sung dữ liệu virome trên vi nấm lâm sàng: Khám phá các dòng mycovirus RNA mới trên các chủng nấm gây bệnh phân lập từ bệnh nhân.

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

  • Hạn chế:
    • Kỹ thuật phân tích Cluster-profiling chỉ nhận diện hiệu quả các tác nhân chưa biết khi chúng hiện diện ở tải lượng sinh học cao (high viral copy number/load); các trường hợp nhiễm virus tải lượng thấp nếu không có tương đồng chuỗi vẫn khó phát hiện.
    • Các kết luận về vật chủ của một số dòng (CRESSV1, CRESSV19, Kirkoviridae) và Draupnirviridae được suy luận dựa trên dữ liệu tính toán phát sinh loài, phân tích HGT và EVEs, chưa thực hiện được nuôi cấy đồng nhiễm in vitro do sự phức tạp trong việc duy trì hệ thống tế bào nuôi cấy ký sinh trùng kỵ khí.
  • Hướng nghiên cứu tiếp:
    • Thử nghiệm nuôi cấy in vitro Entamoeba và Giardia nhiễm virus để đánh giá tác động trực tiếp của Naryaviridae, Nenyaviridae và Vilyaviridae lên độc lực, khả năng gây bệnh và biểu hiện lâm sàng ở người.
    • Nghiên cứu cơ chế phân tử của hệ thống RNA can thiệp (AGO2-2) ở ký sinh trùng trong việc kiểm soát sự nhân lên của CRESS virus.
    • Ứng dụng quy trình tính toán trên các tập dữ liệu metagenome đại dương và đất quy mô lớn để tiếp tục giải mã "vùng tối trình tự" của giới virus.

Giá trị tham khảo

  • Đối với nghiên cứu sinh và nhà nghiên cứu Tin sinh học / Vi sinh học: Cung cấp tài liệu tham khảo về việc thiết kế pipeline xử lý dữ liệu metagenomics HTS, kỹ thuật gom cụm trình tự, phân tích phát sinh loài phân tử và phương pháp phát hiện chuyển gen ngang (HGT/EVEs).
  • Đối với chuyên gia Ký sinh trùng học và Bệnh truyền nhiễm: Cung cấp cơ sở dữ liệu về virome của các đơn bào đường ruột (Entamoeba, Giardia), mở ra hướng nghiên cứu về tương tác ba bên (Người – Ký sinh trùng – Virus) trong sinh bệnh học tiêu hóa.
  • Đối với các nhà phân loại học virus (ICTV): Bổ sung các bằng chứng phát sinh loài và cấu trúc gen vững chắc phục vụ việc thành lập và chuẩn hóa danh pháp các họ mới trong ngành Cressdnaviricota.

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

1. Phân tích Cluster-profiling trong quy trình VIDISCA hoạt động dựa trên nguyên lý nào?

Cluster-profiling khai thác cơ chế sinh học phân tử của enzyme cắt giới hạn $Mse1$ trong quy trình chuẩn bị thư viện VIDISCA. Vì enzyme $Mse1$ cắt DNA tại các vị trí nhận diện cố định (TTAA), các bộ gen virus ngắn có mật độ bản sao cao sẽ tạo ra hàng loạt đoạn cắt giống hệt nhau về điểm đầu và điểm cuối (biological replicates). Khi giải trình tự, các đoạn này tạo thành các cụm đọc (clusters) có mức độ lặp lại rất cao sau khi gom cụm bằng CD-HIT ở ngưỡng $100%$. Công cụ biểu diễn trực quan các cụm này dưới dạng biểu đồ cột phân loại mà không cần dựa vào mức độ tương đồng trình tự, cho phép nhận diện các virus mới phân kỳ cao bị thuật toán căn chỉnh bỏ sót.

2. Bằng chứng nào chứng minh các CRESS virus phát hiện trong phân người thực sự lây nhiễm ký sinh trùng Entamoeba và Giardia mà không phải lây nhiễm người?

Bằng chứng bao gồm:

  1. Dịch tễ học phân tử: $100%$ mẫu phân dương tính với virus ($21/21$) đều dương tính với Entamoeba hoặc Giardia ($\chi^2 = 36.5$, $p = 1.5 \times 10^{-9}$), trong khi các đơn bào khác vắng mặt.
  2. Yếu tố virus nội sinh (EVEs): Tìm thấy các đoạn tích hợp nguyên vẹn cả gen Rep và gen Cap của virus nằm lọt trong các contig nhiễm sắc thể của E. dispar ($100%$ độ phủ, $84%$ danh tính) và G. duodenalis ($99%$ độ phủ, $73%$ danh tính), bảo tồn qua nhiều chủng phân lập độc lập từ các quốc gia và thời điểm khác nhau.
  3. Miễn dịch học phân tử: Các sRNA gắn với phức hợp Argonaute (AGO2-2) của E. histolytica khớp chính xác vào các vùng gen EVE và virus ngoại sinh, thể hiện phản ứng RNA can thiệp chống virus của ký sinh trùng.
  4. Tương quan nucleotide: Hàm lượng %GC của virus tương quan thuận mạnh mẽ với %GC của hệ gen ký sinh trùng ($r^2 = 0.69$, $p = 0.01$).

3. Tại sao sự hoán đổi mô-đun gen Cap giữa Naryaviridae và Nenyaviridae lại chứng minh chúng có chung vật chủ?

Quá trình tái tổ hợp và hoán đổi mô-đun gen cấu trúc (Cap) và gen sao chép (Rep) giữa hai dòng virus chỉ có thể diễn ra khi cả hai virus cùng xâm nhập và nhân lên đồng thời trong cùng một tế bào vật chủ. Phân tích phát sinh loài cho thấy các thành viên của cả hai họ Naryaviridae và Nenyaviridae đều có thể mang vỏ capsid loại CCA1 (dài ~467 aa) hoặc CCA2 (dài ~180 aa) với đầy đủ 4 tổ hợp Rep-Cap có thể có, chứng minh chúng cùng trao đổi vật liệu di truyền trong tế bào Entamoeba. Ngược lại, họ Vilyaviridae chỉ mang vỏ CCA3 và không hề có dấu hiệu tái tổ hợp với hai họ trên, chứng minh chúng thuộc một phân nhóm vật chủ riêng biệt (Giardia).

4. Dòng cressdnavirus CRESSV3 liên quan thế nào đến họ Draupnirviridae và động vật có xương sống?

Nghiên cứu phát hiện dòng CRESSV3 đã chuyển giao ngang (HGT) gen Rep của mình sang bộ gen của các virus thuộc chi Avipoxvirus (họ Poxviridae – các virus DNA chuỗi đôi lớn chuyên gây bệnh trên chim và các loài bò sát thuộc nhánh saurian). Do avipoxvirus chỉ nhân lên trong tế bào động vật có xương sống saurian, sự kiện bắt giữ gen này chứng minh tổ tiên của dòng CRESSV3 đã cùng tồn tại và lây nhiễm trên động vật saurian. Dòng này được đổi tên thành họ Draupnirviridae, đại diện cho dòng cressdnavirus thứ hai (cùng với Circoviridae) được xác định lây nhiễm động vật có xương sống.


Kết luận

Luận án tiến sĩ của Cormac Michael Kinsella đã xây dựng thành công quy trình tin sinh học tích hợp phân tích Cluster-profiling để tối ưu hóa việc phát hiện virus phân kỳ từ dữ liệu giải trình tự siêu hệ gen VIDISCA-NGS. Bằng việc kết hợp dịch tễ học phân tử, phân tích yếu tố virus nội sinh (EVEs), đáp ứng sRNA và mạng lưới tái tổ hợp mô-đun gen, tác giả đã giải quyết thành công bài toán xác định vật chủ cho các dòng CRESS virus chưa phân loại, chính thức thiết lập ba họ virus mới lây nhiễm ký sinh trùng đường ruột (Naryaviridae, Nenyaviridae, Vilyaviridae) và họ lây nhiễm động vật có xương sống (Draupnirviridae). Công trình cung cấp hệ thống phương pháp luận tính toán quan trọng và mở rộng hiểu biết về sự đa dạng, lịch sử tiến hóa và tính hướng vật chủ của ngành Cressdnaviricota.