Diversity and Long-Term Dynamics of Human Blood Anelloviruses.

Anelloviruses (AVs) are commensal members of the human blood virome. Even though it was estimated that over 90% of the human population carries AVs, the dynamics of the AV virome ("anellome") are unknown. We investigated the dynamics of blood anellomes in two healthy people followed up for more than 30 years. Both subjects were positive for AVs in the majority of samples. Alphatorquevirus (torque teno virus [TTV]) was the most common genus in both subjects, followed by Betatorquevirus (torque teno minivirus [TTMV]) and Gammatorquevirus (torque teno midivirus [TTMDV]). Almost five times more lineages were found in subject 1 than in subject 2, and the anellomes differed phylogenetically. Both anellomes remained compositionally stable, and 9 out of 64 AV lineages were detected in over half of the time points. We confirmed the long-term and short-term persistence of 13 lineages by specific quantitative PCR (qPCR). AV lineages were detected in blood for over 30 years. Noticeable differences in anellome richness were found between the tested subjects, but both anellomes remained compositionally stable over time. These findings demonstrate that the human blood anellome is personal and that AV infection is chronic and potentially commensal. IMPORTANCE Knowledge of the persistence of AVs in humans is crucial to our understanding of the nature of AV infection (chronic or acute) and the role of AV in the host. We therefore investigated the dynamics of anellovirus infection in two healthy people followed up for 30 years. Our findings suggest that the human blood anellovirus virome (anellome) remains stable and personal for decades.

Novel Orthobunyavirus Identified in the Cerebrospinal Fluid of a Ugandan Child With Severe Encephalopathy.

A Ugandan child with an unexplained encephalitis was investigated using viral metagenomics. Several sequences from all segments of a novel orthobunyavirus were found. The S-segment, used for typing, showed 41% amino acid diversity to its closest relative. The virus was named Ntwetwe virus, after the hometown of the patient.

A novel genus in the order Picornavirales detected in human stool.

Recently, four new viruses belonging to an unassigned family within the order Picornavirales were identified in excrements of healthy carp (fisavirus) and pigs (posavirus 1, 2 and 3). We report the detection and characterization of a fifth virus present in human faeces. The virus, named human stool-associated RNA virus (husavirus), contains a single ORF encoding a putative 2993 AA polyprotein, with a Hel-Pro-Pol replication block, typical for the Picornavirales. Phylogenetic analysis revealed that the closest relative to husavirus is posavirus 1, and together they cluster with fisavirus, posavirus 2 and 3 and a roundworm (Ascaris suum) derived virus. Husavirus was detected in eight human stool samples collected in 1984 (n52), 1985 (n54), 1995 (n51) and 2014 (n51). From three strains of husavirus from 1984 and 1985 the full genome sequence was determined, showing less than 5% intraspecies variation in the nucleotide composition. The host of this virus remains to be determined.

Isolation and characterization of current human coronavirus strains in primary human epithelial cell cultures reveal differences in target cell tropism.

The human airway epithelium (HAE) represents the entry port of many human respiratory viruses, including human coronaviruses (HCoVs). Nowadays, four HCoVs, HCoV-229E, HCoV-OC43, HCoV-HKU1, and HCoV-NL63, are known to be circulating worldwide, causing upper and lower respiratory tract infections in nonhospitalized and hospitalized children. Studies of the fundamental aspects of these HCoV infections at the primary entry port, such as cell tropism, are seriously hampered by the lack of a universal culture system or suitable animal models. To expand the knowledge on fundamental virus-host interactions for all four HCoVs at the site of primary infection, we used pseudostratified HAE cell cultures to isolate and characterize representative clinical HCoV strains directly from nasopharyngeal material. Ten contemporary isolates were obtained, representing HCoV-229E (n = 1), HCoV-NL63 (n = 1), HCoV-HKU1 (n = 4), and HCoV-OC43 (n = 4). For each strain, we analyzed the replication kinetics and progeny virus release on HAE cell cultures derived from different donors. Surprisingly, by visualizing HCoV infection by confocal microscopy, we observed that HCoV-229E employs a target cell tropism for nonciliated cells, whereas HCoV-OC43, HCoV-HKU1, and HCoV-NL63 all infect ciliated cells. Collectively, the data demonstrate that HAE cell cultures, which morphologically and functionally resemble human airways in vivo, represent a robust universal culture system for isolating and comparing all contemporary HCoV strains.

Identification and characterization of unrecognized viruses in stool samples of non-polio acute flaccid paralysis children by simplified VIDISCA.

BACKGROUND: The use of sequence independent methods combined with next generation sequencing for identification purposes in clinical samples appears promising and exciting results have been achieved to understand unexplained infections. One sequence independent method, Virus Discovery based on cDNA Amplified Fragment Length Polymorphism (VIDISCA) is capable of identifying viruses that would have remained unidentified in standard diagnostics or cell cultures. METHODS: VIDISCA is normally combined with next generation sequencing, however, we set up a simplified VIDISCA which can be used in case next generation sequencing is not possible. Stool samples of 10 patients with unexplained acute flaccid paralysis showing cytopathic effect in rhabdomyosarcoma cells and/or mouse cells were used to test the efficiency of this method. To further characterize the viruses, VIDISCA-positive samples were amplified and sequenced with gene specific primers. RESULTS: Simplified VIDISCA detected seven viruses (70%) and the proportion of eukaryotic viral sequences from each sample ranged from 8.3 to 45.8%. Human enterovirus EV-B97, EV-B100, echovirus-9 and echovirus-21, human parechovirus type-3, human astrovirus probably a type-3/5 recombinant, and tetnovirus-1 were identified. Phylogenetic analysis based on the VP1 region demonstrated that the human enteroviruses are more divergent isolates circulating in the community. CONCLUSION: Our data support that a simplified VIDISCA protocol can efficiently identify unrecognized viruses grown in cell culture with low cost, limited time without need of advanced technical expertise. Also complex data interpretation is avoided thus the method can be used as a powerful diagnostic tool in limited resources. Redesigning the routine diagnostics might lead to additional detection of previously undiagnosed viruses in clinical samples of patients.

Unexplained diarrhoea in HIV-1 infected individuals.

BACKGROUND: Gastrointestinal symptoms, in particular diarrhoea, are common in non-treated HIV-1 infected individuals. Although various enteric pathogens have been implicated, the aetiology of diarrhoea remains unexplained in a large proportion of HIV-1 infected patients. Our aim is to identify the cause of diarrhoea for patients that remain negative in routine diagnostics. METHODS: In this study stool samples of 196 HIV-1 infected persons, including 29 persons with diarrhoea, were examined for enteropathogens and HIV-1. A search for unknown and unexpected viruses was performed using virus discovery cDNA-AFLP combined with Roche-454 sequencing (VIDISCA-454). RESULTS: HIV-1 RNA was detected in stool of 19 patients with diarrhoea (66%) compared to 75 patients (45%) without diarrhoea. In 19 of the 29 diarrhoea cases a known enteropathogen could be identified (66%). Next to these known causative agents, a range of recently identified viruses was identified via VIDISCA-454: cosavirus, Aichi virus, human gyrovirus, and non-A non-B hepatitis virus. Moreover, a novel virus was detected which was named immunodeficiency-associated stool virus (IASvirus). However, PCR based screening for these viruses showed that none of these novel viruses was associated with diarrhoea. Notably, among the 34% enteropathogen-negative cases, HIV-1 RNA shedding in stool was more frequently observed (80%) compared to enteropathogen-positive cases (47%), indicating that HIV-1 itself is the most likely candidate to be involved in diarrhoea. CONCLUSION: Unexplained diarrhoea in HIV-1 infected patients is probably not caused by recently described or previously unknown pathogens, but it is more likely that HIV-1 itself plays a role in intestinal mucosal abnormalities which leads to diarrhoea.

Influenza-like illness symptoms due to endemic human coronavirus reinfections are not influenced by the length of the interval separating reinfections.

After 3 years of its introduction to humans, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been declared as endemic. Little is known about the severity of the disease manifestation that future infections may cause, especially when reinfections occur after humoral immunity from a previous infection or vaccination has waned. Such knowledge could inform policymakers regarding the frequency of vaccination. Reinfections by endemic human coronaviruses (HCoVs) can serve as a model system for SARS-CoV-2 endemicity. We monitored 44 immunocompetent male adults with blood sampling every 6 months (for 17 years), for the frequency of HCoV (re-)infections, using rises in N-antibodies of HCoV-NL63, HCoV-29E, HCoV-OC43, and HCoV-HKU1 as markers of infection. Disease associations during (re-)infections were examined by comparison of self-reporting records of influenza-like illness (ILI) symptoms, every 6 months, by all participants. During 8,549 follow-up months, we found 364 infections by any HCoV with a median of eight infections per person. Symptoms more frequently reported during HCoV infection were cough, sore throat, and myalgia. Two hundred fifty-one of the 364 infections were species-specific HCoV-reinfections, with a median interval of 3.58 (interquartile range 1.92-5.67) years. The length of the interval between reinfections-being either short or long-had no influence on the frequency of reporting ILI symptoms. All HCoV-NL63, HCoV-229E, HCoV-OC43, and HCoV-HKU1 (re-)infections are associated with the reporting of ILIs. Importantly, in immunocompetent males, these symptoms are not influenced by the length of the interval between reinfections. IMPORTANCE: Little is known about the disease following human coronavirus (HCoV) reinfection occurring years after the previous infection, once humoral immunity has waned. We monitored endemic HCoV reinfection in immunocompetent male adults for up to 17 years. We found no influence of reinfection interval length in the disease manifestation, suggesting that immunocompetent male adults are adequately protected against future HCoV infections.

Identification of a new genotype of Torque Teno Mini virus.

BACKGROUND: Although human torque teno viruses (TTVs) were first discovered in 1997, still many associated aspects are not clarified yet. The viruses reveal a remarkable heterogeneity and it is possible that some genotypes are more pathogenic than others. The identification of all genotypes is essential to confirm previous pathogenicity data, and an unbiased search for novel viruses is needed to identify TTVs that might be related to disease. METHOD: The virus discovery technique VIDISCA-454 was used to screen serum of 55 HIV-1 positive injecting drug users, from the Amsterdam Cohort Studies, in search for novel blood-blood transmittable viruses which are undetectable via normal diagnostics or panvirus-primer PCRs. RESULTS: A novel torque teno mini virus (TTMV) was identified in two patients and the sequence of the full genomes were determined. The virus is significantly different from the known TTMVs (< 40% amino acid identity in ORF1), yet it contains conserved characteristics that are also present in other TTMVs. The virus is chronically present in both patients, and these patients both suffered from a pneumococcal pneumonia during follow up and had extremely low B-cells counts. CONCLUSION: We describe a novel TTMV which we tentatively named TTMV-13. Further research is needed to address the epidemiology and pathogenicity of this novel virus.

A Novel Simian Adenovirus Associating with Human Adeno-virus Species G Isolated from Long-Tailed Macaque Feces.

Metagenomics has demonstrated its capability in outbreak investigations and pathogen surveillance and discovery. With high-throughput and effective bioinformatics, many disease-causing agents, as well as novel viruses of humans and animals, have been identified using metagenomic analysis. In this study, a VIDISCA metagenomics workflow was used to identify potential unknown viruses in 33 fecal samples from asymptomatic long-tailed macaques (Macaca fascicularis) in Ratchaburi Province, Thailand. Putatively novel astroviruses, enteroviruses, and adenoviruses were detected and confirmed by PCR analysis of long-tailed macaque fecal samples collected from areas in four provinces, Ratchaburi, Kanchanaburi, Lopburi, and Prachuap Khiri Khan, where humans and monkeys live in proximity (total n = 187). Astroviruses, enteroviruses, and adenoviruses were present in 3.2%, 7.5%, and 4.8% of macaque fecal samples, respectively. One adenovirus, named AdV-RBR-6-3, was successfully isolated in human cell culture. Whole-genome analysis suggested that it is a new member of the species Human adenovirus G, closely related to Rhesus adenovirus 53, with evidence of genetic recombination and variation in the hexon, fiber, and CR1 genes. Sero-surveillance showed neutralizing antibodies against AdV-RBR-6-3 in 2.9% and 11.2% of monkeys and humans, respectively, suggesting cross-species infection of monkeys and humans. Overall, we reported the use of metagenomics to screen for possible new viruses, as well as the isolation and molecular and serological characterization of the new adenovirus with cross-species transmission potential. The findings emphasize that zoonotic surveillance is important and should be continued, especially in areas where humans and animals interact, to predict and prevent the threat of emerging zoonotic pathogens.

Correction: Kosoltanapiwat et al. A Novel Simian Adenovirus Associating with Human Adenovirus Species G Isolated from Long-Tailed Macaque Feces. Viruses 2023, 15, 1371.

After publication of the article, the authors received comments from a member of the Viruses editorial board who is an expert in the field of adenovirus concerning figures and references that should be included in the paper [...].

Replication-dependent downregulation of cellular angiotensin-converting enzyme 2 protein expression by human coronavirus NL63.

Like severe acute respiratory syndrome coronavirus (SARS-CoV), human coronavirus (HCoV)-NL63 employs angiotensin-converting enzyme 2 (ACE2) as a receptor for cellular entry. SARS-CoV infection causes robust downregulation of cellular ACE2 expression levels and it has been suggested that the SARS-CoV effect on ACE2 is involved in the severity of disease. We investigated whether cellular ACE2 downregulation occurs at optimal replication conditions of HCoV-NL63 infection. The expression of the homologue of ACE2, the ACE protein not used as a receptor by HCoV-NL63, was measured as a control. A specific decrease for ACE2 protein level was observed when HCoV-NL63 was cultured at 34 °C. Culturing the virus at the suboptimal temperature of 37 °C resulted in low replication of the virus and the effect on ACE2 expression was lost. We conclude that the decline of ACE2 expression is dependent on the efficiency of HCoV-NL63 replication, and that HCoV-NL63 and SARS-CoV both affect cellular ACE2 expression during infection.

The first complete genome sequences of clinical isolates of human coronavirus 229E.

Human coronavirus 229E has been identified in the mid-1960s, yet still only one full-genome sequence is available. This full-length sequence has been determined from the cDNA-clone Inf-1 that is based on the lab-adapted strain VR-740. Lab-adaptation might have resulted in genomic changes, due to insufficient pressure to maintain gene integrity of non-essential genes. We present here the first full-length genome sequence of two clinical isolates. Each encoded gene was compared to Inf-1. In general, little sequence changes were noted, most could be attributed to genetic drift, since the clinical isolates originate from 2009 to 2010 and VR740 from 1962. Hot spots of substitutions were situated in the S1 region of the Spike, the nucleocapsid gene, and the non-structural protein 3 gene, whereas several deletions were detected in the 3'UTR. Most notable was the difference in genome organization: instead of an ORF4A and ORF4B, an intact ORF4 was present in clinical isolates.

Divergent Rhabdovirus Discovered in a Patient with New-Onset Nodding Syndrome.

A divergent rhabdovirus was discovered in the bloodstream of a 15-year-old girl with Nodding syndrome from Mundri West County in South Sudan. Nodding syndrome is a progressive degenerative neuropathy of unknown cause affecting thousands of individuals in Sub-Saharan Africa. The index case was previously healthy until she developed head-nodding seizures four months prior to presentation. Virus discovery by VIDISCA-NGS on the patient's plasma detected multiple sequence reads belonging to a divergent rhabdovirus. The viral load was 3.85 × 103 copies/mL in the patient's plasma and undetectable in her cerebrospinal fluid. Further genome walking allowed for the characterization of full coding sequences of all the viral proteins (N, P, M, U1, U2, G, U3, and L). We tentatively named the virus "Mundri virus" (MUNV) and classified it as a novel virus species based on the high divergence from other known viruses (all proteins had less than 43% amino acid identity). Phylogenetic analysis revealed that MUNV forms a monophyletic clade with several human-infecting tibroviruses prevalent in Central Africa. A bioinformatic machine-learning algorithm predicted MUNV to be an arbovirus (bagged prediction strength (BPS) of 0.9) transmitted by midges (BPS 0.4) with an artiodactyl host reservoir (BPS 0.9). An association between MUNV infection and Nodding syndrome was evaluated in a case-control study of 72 patients with Nodding syndrome (including the index case) matched to 65 healthy households and 48 community controls. No subject, besides the index case, was positive for MUNV RNA in their plasma. A serological assay detecting MUNV anti-nucleocapsid found, respectively, in 28%, 22%, and 16% of cases, household controls and community controls to be seropositive with no significant differences between cases and either control group. This suggests that MUNV commonly infects children in South Sudan yet may not be causally associated with Nodding syndrome.

Transmission of anelloviruses to HIV-1 infected children.

Anelloviruses (AVs) are widespread in the population and infect humans at the early stage of life. The mode of transmission of AVs is still unknown, however, mother-to-child transmission, e.g., via breastfeeding, is one of the likely infection routes. To determine whether the mother-to-child transmission of AVs may still occur despite the absence of natural birth and breastfeeding, 29 serum samples from five HIV-1-positive mother and child pairs were Illumina-sequenced. The Illumina reads were mapped to an AV lineage database "Anellometrix" containing 502 distinct ORF1 sequences. Although the majority of lineages from the mother were not shared with the child, the mother and child anellomes did display a significant similarity. These findings suggest that AVs may be transmitted from mothers to their children via different routes than delivery or breastfeeding.

Atypical Antibody Dynamics During Human Coronavirus HKU1 Infections.

Human coronavirus HKU1 (HCoV-HKU1) is one of the four endemic coronaviruses. It has been suggested that there is a difference in incidence, with PCR-confirmed HCoV-NL63 and HCoV-OC43 infections occurring more commonly, whereas HCoV-HKU1 is the least seen. Lower incidence of HCoV-HKU1 infection has also been observed in serological studies. The current study aimed to investigate antibody dynamics during PCR-confirmed HCoV-HKU1 infections using serum collected during infection and 1 month later. We expressed a new HCoV-HKU1 antigen consisting of both the linker and carboxy-terminal domain of the viral nucleocapsid protein and implemented it in ELISA. We also applied a spike-based Luminex assay on serum samples from PCR-confirmed infections by the four endemic HCoVs. At least half of HCoV-HKU1-infected subjects consistently showed no antibody rise via either assay, and some subjects even exhibited substantial antibody decline. Investigation of self-reported symptoms revealed that HCoV-HKU1-infected subjects rated their illness milder than subjects infected by other HCoVs. In conclusion, HCoV-HKU1 infections reported in this study displayed atypical antibody dynamics and milder symptoms when compared to the other endemic HCoVs.

Early-Life Colonization by Anelloviruses in Infants.

Anelloviruses (AVs) are found in the vast majority of the human population and are most probably part of a healthy virome. These viruses infect humans in the early stage of life, however, the characteristics of the first colonizing AVs are still unknown. We screened a collection of 107 blood samples from children between 0.4 and 64.8 months of age for the presence of three AV genera: the Alpha-, Beta- and Gammatorquevirus. The youngest child that was positive for AV was 1.2 months old, and a peak in prevalence (100% of samples positive) was reached between the twelfth and eighteenth months of life. Intriguingly, the beta- and gammatorqueviruses were detected most at the early stage of life (up to 12 months), whereas alphatorqueviruses, the most common AVs in adults, increased in prevalence in children older than 12 months. To determine whether that order of colonization may be related to oral transmission and unequal presence of AV genera in breast milk, we examined 63 breast milk samples. Thirty-two percent of the breast milk samples were positive in a qPCR detecting beta- and gammatorqueviruses, while alphatorqueviruses were detected in 10% of the samples, and this difference was significant (p = 0.00654). In conclusion, we show that beta- and gammatorqueviruses colonize humans in the first months of life and that breastfeeding could play a role in AV transmission.

Culturing the unculturable: human coronavirus HKU1 infects, replicates, and produces progeny virions in human ciliated airway epithelial cell cultures.

Culturing newly identified human lung pathogens from clinical sample isolates can represent a daunting task, with problems ranging from low levels of pathogens to the presence of growth suppressive factors in the specimens, compounded by the lack of a suitable tissue culture system. However, it is critical to develop suitable in vitro platforms to isolate and characterize the replication kinetics and pathogenesis of recently identified human pathogens. HCoV-HKU1, a human coronavirus identified in a clinical sample from a patient with severe pneumonia, has been a major challenge for successful propagation on all immortalized cells tested to date. To determine if HCoV-HKU1 could replicate in in vitro models of human ciliated airway epithelial cell cultures (HAE) that recapitulate the morphology, biochemistry, and physiology of the human airway epithelium, the apical surfaces of HAE were inoculated with a clinical sample of HCoV-HKU1 (Cean1 strain). High virus yields were found for several days postinoculation and electron micrograph, Northern blot, and immunofluorescence data confirmed that HCoV-HKU1 replicated efficiently within ciliated cells, demonstrating that this cell type is infected by all human coronaviruses identified to date. Antiserum directed against human leukocyte antigen C (HLA-C) failed to attenuate HCoV-HKU1 infection and replication in HAE, suggesting that HLA-C is not required for HCoV-HKU1 infection of the human ciliated airway epithelium. We propose that the HAE model provides a ready platform for molecular studies and characterization of HCoV-HKU1 and in general serves as a robust technology for the recovery, amplification, adaptation, and characterization of novel coronaviruses and other respiratory viruses from clinical material.

Identification of a new human coronavirus.

Three human coronaviruses are known to exist: human coronavirus 229E (HCoV-229E), HCoV-OC43 and severe acute respiratory syndrome (SARS)-associated coronavirus (SARS-CoV). Here we report the identification of a fourth human coronavirus, HCoV-NL63, using a new method of virus discovery. The virus was isolated from a 7-month-old child suffering from bronchiolitis and conjunctivitis. The complete genome sequence indicates that this virus is not a recombinant, but rather a new group 1 coronavirus. The in vitro host cell range of HCoV-NL63 is notable because it replicates on tertiary monkey kidney cells and the monkey kidney LLC-MK2 cell line. The viral genome contains distinctive features, including a unique N-terminal fragment within the spike protein. Screening of clinical specimens from individuals suffering from respiratory illness identified seven additional HCoV-NL63-infected individuals, indicating that the virus was widely spread within the human population.

Integrative Transkingdom Analysis of the Gut Microbiome in Antibiotic Perturbation and Critical Illness.

Bacterial microbiota play a critical role in mediating local and systemic immunity, and shifts in these microbial communities have been linked to impaired outcomes in critical illness. Emerging data indicate that other intestinal organisms, including bacteriophages, viruses of eukaryotes, fungi, and protozoa, are closely interlinked with the bacterial microbiota and their host, yet their collective role during antibiotic perturbation and critical illness remains to be elucidated. We employed multi-omics factor analysis (MOFA) to systematically integrate the bacterial (16S rRNA), fungal (intergenic transcribed spacer 1 rRNA), and viral (virus discovery next-generation sequencing) components of the intestinal microbiota of 33 critically ill patients with and without sepsis and 13 healthy volunteers. In addition, we quantified the absolute abundances of bacteria and fungi using 16S and 18S rRNA PCRs and characterized the short-chain fatty acids (SCFAs) butyrate, acetate, and propionate using nuclear magnetic resonance spectroscopy. We observe that a loss of the anaerobic intestinal environment is directly correlated with an overgrowth of aerobic pathobionts and their corresponding bacteriophages as well as an absolute enrichment of opportunistic yeasts capable of causing invasive disease. We also observed a strong depletion of SCFAs in both disease states, which was associated with an increased absolute abundance of fungi with respect to bacteria. Therefore, these findings illustrate the complexity of transkingdom changes following disruption of the intestinal bacterial microbiome.IMPORTANCE While numerous studies have characterized antibiotic-induced disruptions of the bacterial microbiome, few studies describe how these disruptions impact the composition of other kingdoms such as viruses, fungi, and protozoa. To address this knowledge gap, we employed MOFA to systematically integrate viral, fungal, and bacterial sequence data from critically ill patients (with and without sepsis) and healthy volunteers, both prior to and following exposure to broad-spectrum antibiotics. In doing so, we show that modulation of the bacterial component of the microbiome has implications extending beyond this kingdom alone, enabling the overgrowth of potentially invasive fungi and viruses. While numerous preclinical studies have described similar findings in vitro, we confirm these observations in humans using an integrative analytic approach. These findings underscore the potential value of multi-omics data integration tools in interrogating how different components of the microbiota contribute to disease states. In addition, our findings suggest that there is value in further studying potential adjunctive therapies using anaerobic bacteria or SCFAs to reduce fungal expansion after antibiotic exposure, which could ultimately lead to improved outcomes in the intensive care unit (ICU).

Entamoeba and Giardia parasites implicated as hosts of CRESS viruses.

Metagenomic techniques have enabled genome sequencing of unknown viruses without isolation in cell culture, but information on the virus host is often lacking, preventing viral characterisation. High-throughput methods capable of identifying virus hosts based on genomic data alone would aid evaluation of their medical or biological relevance. Here, we address this by linking metagenomic discovery of three virus families in human stool samples with determination of probable hosts. Recombination between viruses provides evidence of a shared host, in which genetic exchange occurs. We utilise networks of viral recombination to delimit virus-host clusters, which are then anchored to specific hosts using (1) statistical association to a host organism in clinical samples, (2) endogenous viral elements in host genomes, and (3) evidence of host small RNA responses to these elements. This analysis suggests two CRESS virus families (Naryaviridae and Nenyaviridae) infect Entamoeba parasites, while a third (Vilyaviridae) infects Giardia duodenalis. The trio supplements five CRESS virus families already known to infect eukaryotes, extending the CRESS virus host range to protozoa. Phylogenetic analysis implies CRESS viruses infecting multicellular life have evolved independently on at least three occasions.