Complex virome in feces from Amerindian children in isolated Amazonian villages

2024 taxonomic update for the families Naryaviridae, Nenyaviridae, and Vilyaviridae

The families Naryaviridae (order Rivendellvirales), Nenyaviridae (order Rohanvirales), and Vilyaviridae (order Cirlivirales), all within the class Arfiviricetes of the phylum Cressdnaviricota, include single-stranded DNA viruses associated with protozoan parasites of the genera Entamoeba and Giardia as well as viruses found in various environmental samples, also likely infecting protozoans. Here, we provide a taxonomic update for these three families, which were recently expanded with multiple new members. In particular, we established seven new genera and nine new species in the family Naryaviridae, one new genus with one new species in the family Nenyaviridae, and three new genera and nine new species in the family Vilyaviridae. We also summarize the genomic properties and protein characteristics, including conserved motifs of the rolling-circle replication initiation proteins, of the viruses in the three families. Notably, the high GC content of vilyavirids (51-61%) and considerably lower GC content of naryavirids and nenyavirids (33-44%) appear to represent an adaptation to their hosts, Giardia and Entamoeba species, respectively.

Alterations of the gut microbiome in HIV infection highlight human anelloviruses as potential predictors of immune recovery

BACKGROUND

HIV-1 infection is characterized by a massive depletion of mucosal CD4 T cells that triggers a cascade of events ultimately linking gut microbial dysbiosis to HIV-1 disease progression and pathogenesis. The association between HIV infection and the enteric virome composition is less characterized, although viruses are an essential component of the gut ecosystem. Here, we performed a cross-sectional analysis of the fecal viral (eukaryotic viruses and bacteriophages) and bacterial microbiome in people with HIV (PWH) and in HIV-negative individuals. To gain further insight into the association between the gut microbiome composition, HIV-associated immunodeficiency, and immune recovery, we carried out a longitudinal study including 14 PWH who initiated antiretroviral therapy (ART) and were followed for 24 months with samplings performed at baseline (before ART) and at 2, 6, 12, and 24 months post-ART initiation.

RESULTS

Our data revealed a striking expansion in the abundance and prevalence of several human virus genomic sequences (Anelloviridae, Adenoviridae, and Papillomaviridae) in stool samples of PWH with severe immunodeficiency (CD4 < 200). We also noted a decreased abundance of sequences belonging to two plant viruses from the Tobamovirus genus, a reduction in bacterial alpha diversity, and a decrease in Inoviridae bacteriophage sequences. Short-term ART (24 months) was linked to a significant decrease in human Anelloviridae sequences. Remarkably, the detection of Anellovirus sequences at baseline independently predicted poor immune recovery, as did low CD4 T cell counts. The bacterial and bacteriophage populations were unique to each PWH with individualized trajectories; we found no discernable pattern of clustering after 24 months on ART.

CONCLUSION

Advanced HIV-1 infection was associated with marked alterations in the virome composition, in particular a remarkable expansion of human anelloviruses, with a gradual restoration after ART initiation. In addition to CD4 T cell counts, anellovirus sequence detection might be useful to predict and monitor immune recovery. This study confirms data on the bacteriome and expands our knowledge on the viral component of the gut microbiome in HIV-1 infection. Video Abstract.

Genetic diversity and cross-species transmissibility of bat-associated picornaviruses from Spain

BACKGROUND

Emerging zoonotic diseases arise from cross-species transmission events between wild or domesticated animals and humans, with bats being one of the major reservoirs of zoonotic viruses. Viral metagenomics has led to the discovery of many viruses, but efforts have mainly been focused on some areas of the world and on certain viral families.

METHODS

We set out to describe full-length genomes of new picorna-like viruses by collecting feces from hundreds of bats captured in different regions of Spain. Viral sequences were obtained by high-throughput Illumina sequencing and analyzed phylogenetically to classify them in the context of known viruses. Linear discriminant analysis (LDA) was performed to infer likely hosts based on genome composition.

RESULTS

We found five complete or nearly complete genomes belonging to the family Picornaviridae, including a new species of the subfamily Ensavirinae. LDA suggested that these were true vertebrate viruses, rather than viruses from the bat diet. Some of these viruses were related to picornaviruses previously found in other bat species from distant geographical regions. We also found a calhevirus genome that most likely belongs to a proposed new family within the order Picornavirales, and for which genome composition analysis suggested a plant host.

CONCLUSIONS

Our findings describe new picorna-like viral species and variants circulating in the Iberian Peninsula, illustrate the wide geographical distribution and interspecies transmissibility of picornaviruses, and suggest new hosts for calheviruses.

Meta-Transcriptomic Analysis Reveals Novel RNA Viruses in Polychaetes Perinereis

Perinereis species are essential benthonic animals in coastal ecosystems and have significant roles as live feed in aquaculture, owing to their high-protein and low-fat nutritional profile. Despite their ecological importance, the viral communities associated with these organisms need to be better understood. In this study, we generated 2.6 × 108 reads using meta-transcriptomic sequencing and de novo assembled 5.3 × 103 virus-associated contigs. We identified 12 novel RNA viruses from two species, Perinereis aibuhitensis and P. wilsoni, which were classified into four major viral groups: Picobirnaviridae, Marnaviridae, unclassified Picornavirales, and unclassified Bunyavirales. Our findings revealed the hidden diversity of viruses and genome structures in Perinereis, enriching the RNA virosphere and expanding the host range of Picobirnaviridae, Marnaviridae, and Bunyavirales. This study also highlighted the potential biosecurity risk of the novel viruses carried by Perinereis to aquaculture.

High Detection Frequency of Vaccine-Associated Polioviruses and Non-Polio Enteroviruses in the Stools of Asymptomatic Infants from the Free State Province, South Africa

Enterovirus (EV) infections are widespread and associated with a range of clinical conditions, from encephalitis to meningitis, gastroenteritis, and acute flaccid paralysis. Knowledge about the circulation of EVs in neonatal age and early infancy is scarce, especially in Africa. This study aimed to unveil the frequency and diversity of EVs circulating in apparently healthy newborns from the Free State Province, South Africa (SA). For this purpose, longitudinally collected faecal specimens (May 2021-February 2022) from a cohort of 17 asymptomatic infants were analysed using metagenomic next-generation sequencing. Overall, seven different non-polio EV (NPEV) subtypes belonging to EV-B and EV-C species were identified, while viruses classified under EV-A and EV-D species could not be characterised at the sub-species level. Additionally, under EV-C species, two vaccine-related poliovirus subtypes (PV1 and PV3) were identified. The most prevalent NPEV species was EV-B (16/17, 94.1%), followed by EV-A (3/17, 17.6%), and EV-D (4/17, 23.5%). Within EV-B, the commonly identified NPEV types included echoviruses 6, 13, 15, and 19 (E6, E13, E15, and E19), and coxsackievirus B2 (CVB2), whereas enterovirus C99 (EV-C99) and coxsackievirus A19 (CVA19) were the only two NPEVs identified under EV-C species. Sabin PV1 and PV3 strains were predominantly detected during the first week of birth and 6-8 week time points, respectively, corresponding with the OPV vaccination schedule in South Africa. A total of 11 complete/near-complete genomes were identified from seven NPEV subtypes, and phylogenetic analysis of the three EV-C99 identified revealed that our strains were closely related to other strains from Cameroon and Brazil, suggesting global distribution of these strains. This study provides an insight into the frequency and diversity of EVs circulating in asymptomatic infants from the Free State Province, with the predominance of subtypes from EV-B and EV-C species. This data will be helpful to researchers looking into strategies for the control and treatment of EV infection.

Metagenomic survey of viral diversity obtained from feces of piglets with diarrhea

Pigs are natural host to various zoonotic pathogens including viruses. In this study, we analyzed the viral communities in the feces of 89 piglets with diarrhea under one month old which were collected from six farms in Jiangsu Province of the Eastern China, using the unbiased virus metagenomic method. A total of 89 libraries were constructed, and 46937894 unique sequence reads were generated by Illumina sequencing. Overall, the family Picornaviridae accounted for the majority of the total reads of putative mammalian viruses. Ten novel virus genomes from different family members were discovered, including Parvoviridae (n = 2), Picobirnaviridae (n = 4) and CRESS DNA viruses (n = 4). A large number of phages were identified, which mainly belonged to the order Caudovirales and the family Microviridae. Moreover, some identified viruses were closely related to viruses found in non-porcine hosts, highlighting the potential for cross-species virus dissemination. This study increased our understanding of the fecal virus communities of diarrhea piglets and provided valuable information for virus monitoring and preventing.

Identification of small circular DNA viruses in coyote fecal samples from Arizona (USA)

Coyotes (Canis latrans) have a broad geographic distribution across North and Central America. Despite their widespread presence in urban environments in the USA, there is limited information regarding viruses associated with coyotes in the USA and in particular the state of Arizona. To explore viruses associated with coyotes, particularly small DNA viruses, 44 scat samples were collected (April-June 2021 and November 2021-January 2022) along the Salt River near Phoenix, Arizona (USA), along 43 transects (500 m). From these samples, we identified 11 viral genomes: two novel circoviruses, six unclassified cressdnaviruses, and two anelloviruses. One of the circoviruses is most closely related to a circovirus sequence identified from an aerosolized dust sample in Arizona, USA. The second circovirus is most closely related to a rodent-associated circovirus and canine circovirus. Of the unclassified cressdnaviruses, three encode replication-associated proteins that are similar to those found in protists (Histomonas meleagridis and Monocercomonoides exilis), implying an evolutionary relationship with or a connection to similar unidentified protist hosts. The two anelloviruses are most closely related to those found in rodents, and this suggests a diet-related identification.

Characterization of Rotavirus Strains Responsible for Breakthrough Diarrheal Diseases among Zambian Children Using Whole Genome Sequencing

The occurrence of rotavirus (RV) infection among vaccinated children in high-burden settings poses a threat to further disease burden reduction. Genetically altered viruses have the potential to evade both natural infection and vaccine-induced immune responses, leading to diarrheal diseases among vaccinated children. Studies characterizing RV strains responsible for breakthrough infections in resource-limited countries where RV-associated diarrheal diseases are endemic are limited. We aimed to characterize RV strains detected in fully vaccinated children residing in Zambia using next-generation sequencing. We conducted whole genome sequencing on Illumina MiSeq. Whole genome assembly was performed using Geneious Prime 2023.1.2. A total of 76 diarrheal stool specimens were screened for RV, and 4/76 (5.2%) were RV-positive. Whole genome analysis revealed RVA/Human-wt/ZMB/CIDRZ-RV2088/2020/G1P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2 and RVA/Human-wt/ZMB/CIDRZ-RV2106/2020/G12P[4]-I1-R2-C2-M2-A2-N1-T2-E1-H2 strains were mono and multiple reassortant (exchanged genes in bold) respectively, whilst RVA/Human-wt/ZMB/CIDRZ-RV2150/2020/G12P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 was a typical Wa-like strain. Comparison of VP7 and VP4 antigenic epitope of breakthrough strains and Rotarix strain revealed several amino acid differences. Variations in amino acids in antigenic epitope suggested they played a role in immune evasion of neutralizing antibodies elicited by vaccination. Findings from this study have the potential to inform national RV vaccination strategies and the design of highly efficacious universal RV vaccines.

Isolation of a rhesus calicivirus that can replicate in human cells

Recoviruses (rhesus enteric caliciviruses) are members of the Caliciviridae family. They are a valuable model for studying human caliciviruses such as noroviruses. It has been suggested that some recoviruses may infect humans, which necessitates detailed studies on the cell type tropism of recoviruses. For the recoviruses that have been cultured to date, successful growth has only been reported in monkey kidney cell lines, precluding their use to study virus interactions with human cells. We isolated and characterized a new recovirus, Recovirus Mo/TG30/2012, from monkey stool which grew efficiently in the monkey kidney cell line LLC-MK2. Notably, the virus can infect and replicate in several human cell lines derived from different organs. The ability to infect a human cell culture system with a recovirus expands our understanding of the potential for spillover to humans as well as increases the value of recoviruses as a model of human caliciviruses.

Altered vaginal eukaryotic virome is associated with different cervical disease status

Viruses are important components of the human body. Growing evidence suggests that they are engaged in the physiology and disease status of the host. Even though the vaginal microbiome is involved in human papillomavirus (HPV) infection and cervical cancer (CC) progression, little is known about the role of the vaginal virome. In this pilot exploratory study, using unbiased viral metagenomics, we aim to investigate the vaginal eukaryotic virome in women with different levels of cervical lesions, and examine their associations with different cervical disease status. An altered eukaryotic virome was observed in women with different levels of lesions and Lactobacillus profiles. Anelloviruses and papillomaviruses are the most commonly detected eukaryotic viruses of the vaginal virome. Higher abundance and richness of anelloviruses and papillomaviruses were associated with low-grade squamous intraepithelial lesion (LSIL) and CC. Besides, higher anellovirus abundance was also associated with lactobacillus-depleted microbiome profiles and bacterial community state (CST) type IV. Furthermore, increased correlations between Anelloviridae and Papillomaviridae occurred in the women with increased cervical disease severity level from LSIL to CC. These data suggest underlying interactions between different microbes as well as the host physiology. Higher abundance and diversity of both anelloviruses and papillomaviruses shared by LSIL and CC suggest that anellovirus may be used as a potential adjunct biomarker to predict the risk of HPV persistent infection and/or CC. Future studies need to focus on the clinical relevance of anellovirus abundance with cervical disease status, and the evaluation of their potential as a new adjunct biomarker for the prediction and prognoses of CC.

Naryaviridae, Nenyaviridae, and Vilyaviridae: three new families of single-stranded DNA viruses in the phylum Cressdnaviricota

The phylum Cressdnaviricota includes viruses with circular single-stranded DNA (ssDNA) genomes and icosahedral capsids. These viruses display global environmental distribution and infect diverse eukaryotic hosts, including animals, plants, and fungi. Here, we report on the formal creation of two new orders, Rivendellvirales and Rohanvirales, and three new families, Naryaviridae, Nenyaviridae, and Vilyaviridae, of ssDNA viruses associated with protozoan parasites belonging to the genera Entamoeba and Giardia. We describe a sequence-based taxonomic framework, which was used to classify 60 ssDNA viruses into 12 genera (with 18 species) within the family Vilyaviridae; four genera (with five species) within the family Naryaviridae; and five genera (with six species) within the family Nenyaviridae. We also highlight the challenges associated with the classification of chimeric virus genomes, such as those in the families Naryaviridae and Nenyaviridae, where the replication initiation and capsid protein genes have undergone several independent non-orthologous replacements. The described taxonomic changes have been ratified by the International Committee on Taxonomy of Viruses (ICTV) and expand the phylum Cressdnaviricota to eight orders and 11 families.

Expanding known viral diversity in plants: virome of 161 species alongside an ancient canal

BACKGROUND

Since viral metagenomic approach was applied to discover plant viruses for the first time in 2006, many plant viruses had been identified from cultivated and non-cultivated plants. These previous researches exposed that the viral communities (virome) of plants have still largely uncharacterized. Here, we investigated the virome in 161 species belonging to 38 plant orders found in a riverside ecosystem.

RESULTS

We identified 245 distinct plant-associated virus genomes (88 DNA and 157 RNA viruses) belonging to 27 known viral families, orders, or unclassified virus groups. Some viral genomes were sufficiently divergent to comprise new species, genera, families, or even orders. Some groups of viruses were detected that currently are only known to infect organisms other than plants. It indicates a wider host range for members of these clades than previously recognized theoretically. We cannot rule out that some viruses could be from plant contaminating organisms, although some methods were taken to get rid of them as much as possible. The same viral species could be found in different plants and co-infections were common.

CONCLUSIONS

Our data describe a complex viral community within a single plant ecosystem and expand our understanding of plant-associated viral diversity and their possible host ranges.

A New Circular Single-Stranded DNA Virus Related with Howler Monkey Associated Porprismacovirus 1 Detected in Children with Acute Gastroenteritis

Putative replication-associated protein (REP) and capsid-like (CAP) proteins are encoded by circular single-stranded DNA viruses (CRESS DNA), which have been found in samples from most eukaryotic groups. However, the details of these viruses' life cycles and their significance in diseases have yet to be established. We presented and analyzed two full-length CRESS DNA genomes acquired from two children diagnosed with acute gastroenteritis (GI) in the northeast state of Tocantins, Brazil, using next-generation sequencing and a virus-like filtration approach. Both sequences (named SmaCV3BR08 and SmaCV3BR291) are closely similar to a prior CRESS DNA sequence discovered in the feces of a new world monkey (Alouatta caraya) from the United States in 2009 and termed Howler monkey-associated porprismacovirus 1 (Genbank ID: NC 026317). According to our comparative study, these porprismacovirus genomes deviate by 10% at the nucleotide level. For comparative reasons, the divergence between our sequences (SmaCV3BR08 and SmaCV3BR291) and a porprismacovirus recently identified in a human fecal sample from Peru is 37%. These data suggest that there is a great diversity of porprismacoviruses in South America, perhaps more than two species. In addition, the finding of closely related sequences of porprismacoviruses in humans and native monkeys highlights the zoonotic potential of these viruses.

Virome in the cloaca of wild and breeding birds revealed a diversity of significant viruses

BACKGROUND

Wild birds may harbor and transmit viruses that are potentially pathogenic to humans, domestic animals, and other wildlife.

RESULTS

Using the viral metagenomic approach, we investigated the virome of cloacal swab specimens collected from 3182 birds (the majority of them wild species) consisting of > 87 different species in 10 different orders within the Aves classes. The virus diversity in wild birds was higher than that in breeding birds. We acquired 707 viral genomes from 18 defined families and 4 unclassified virus groups, with 265 virus genomes sharing < 60% protein sequence identities with their best matches in GenBank comprising new virus families, genera, or species. RNA viruses containing the conserved RdRp domain with no phylogenetic affinity to currently defined virus families existed in different bird species. Genomes of the astrovirus, picornavirus, coronavirus, calicivirus, parvovirus, circovirus, retrovirus, and adenovirus families which include known avian pathogens were fully characterized. Putative cross-species transmissions were observed with viruses in wild birds showing > 95% amino acid sequence identity to previously reported viruses in domestic poultry. Genomic recombination was observed for some genomes showing discordant phylogenies based on structural and non-structural regions. Mapping the next-generation sequencing (NGS) data respectively against the 707 genomes revealed that these viruses showed distribution pattern differences among birds with different habitats (breeding or wild), orders, and sampling sites but no significant differences between birds with different behavioral features (migratory and resident).

CONCLUSIONS

The existence of a highly diverse virome highlights the challenges in elucidating the evolution, etiology, and ecology of viruses in wild birds. Video Abstract.

Benchmarking metagenomics classifiers on ancient viral DNA: a simulation study

Owing to technological advances in ancient DNA, it is now possible to sequence viruses from the past to track down their origin and evolution. However, ancient DNA data is considerably more degraded and contaminated than modern data making the identification of ancient viral genomes particularly challenging. Several methods to characterise the modern microbiome (and, within this, the virome) have been developed; in particular, tools that assign sequenced reads to specific taxa in order to characterise the organisms present in a sample of interest. While these existing tools are routinely used in modern data, their performance when applied to ancient microbiome data to screen for ancient viruses remains unknown. In this work, we conducted an extensive simulation study using public viral sequences to establish which tool is the most suitable to screen ancient samples for human DNA viruses. We compared the performance of four widely used classifiers, namely Centrifuge, Kraken2, DIAMOND and MetaPhlAn2, in correctly assigning sequencing reads to the corresponding viruses. To do so, we simulated reads by adding noise typical of ancient DNA to a set of publicly available human DNA viral sequences and to the human genome. We fragmented the DNA into different lengths, added sequencing error and C to T and G to A deamination substitutions at the read termini. Then we measured the resulting sensitivity and precision for all classifiers. Across most simulations, more than 228 out of the 233 simulated viruses were recovered by Centrifuge, Kraken2 and DIAMOND, in contrast to MetaPhlAn2 which recovered only around one third. Overall, Centrifuge and Kraken2 had the best performance with the highest values of sensitivity and precision. We found that deamination damage had little impact on the performance of the classifiers, less than the sequencing error and the length of the reads. Since Centrifuge can handle short reads (in contrast to DIAMOND and Kraken2 with default settings) and since it achieve the highest sensitivity and precision at the species level across all the simulations performed, it is our recommended tool. Regardless of the tool used, our simulations indicate that, for ancient human studies, users should use strict filters to remove all reads of potential human origin. Finally, we recommend that users verify which species are present in the database used, as it might happen that default databases lack sequences for viruses of interest.

Gut Virome: Role and Distribution in Health and Gastrointestinal Diseases

The study of the intestinal microbiome is an evolving field of research that includes comprehensive analysis of the vast array of microbes – bacterial, archaeal, fungal, and viral. Various gastrointestinal (GI) diseases, such as Crohn’s disease and ulcerative colitis, have been associated with instability of the gut microbiota. Many studies have focused on importance of bacterial communities with relation to health and disease in humans. The role of viruses, specifically bacteriophages, have recently begin to emerge and have profound impact on the host. Here, we comprehensively review the importance of viruses in GI diseases and summarize their influence in the complex intestinal environment, including their biochemical and genetic activities. We also discuss the distribution of the gut virome as it relates with treatment and immunological advantages. In conclusion, we suggest the need for further studies on this critical component of the intestinal microbiome to decipher the role of the gut virome in human health and disease.

Viruses associated with ill health in wild chimpanzees

Viral infection is a major cause of ill health in wild chimpanzees (Pan troglodytes), but most evidence to date has come from conspicuous disease outbreaks with high morbidity and mortality. To examine the relationship between viral infection and ill health during periods not associated with disease outbreaks, we conducted a longitudinal study of wild eastern chimpanzees (P. t. schweinfurthii) in the Kanyawara and Ngogo communities of Kibale National Park, Uganda. We collected standardized, observational health data for 4 years and then used metagenomics to characterize gastrointestinal viromes (i.e., all viruses recovered from fecal samples) in individual chimpanzees before and during episodes of clinical disease. We restricted our analyses to viruses thought to infect mammals or primarily associated with mammals, discarding viruses associated with nonmammalian hosts. We found 18 viruses (nine of which were previously identified in this population) from at least five viral families. Viral richness (number of viruses per sample) did not vary by health status. By contrast, total viral load (normalized proportion of sequences mapping to viruses) was significantly higher in ill individuals compared with healthy individuals. Furthermore, when ill, Kanyawara chimpanzees exhibited higher viral loads than Ngogo chimpanzees, and males, but not females, exhibited higher infection rates with certain viruses and higher total viral loads as they aged. Post-hoc analyses, including the use of a machine-learning classification method, indicated that one virus, salivirus (Picornaviridae), was the main contributor to health-related and community-level variation in viral loads. Another virus, chimpanzee stool-associated virus (chisavirus; unclassified Picornavirales), was associated with ill health at Ngogo but not at Kanyawara. Chisavirus, chimpanzee adenovirus (Adenoviridae), and bufavirus (Parvoviridae) were also associated with increased age in males. Associations with sex and age are consistent with the hypothesis that nonlethal viral infections cumulatively reflect or contribute to senescence in long-lived species such as chimpanzees.

Composition of Eukaryotic Viruses and Bacteriophages in Individuals with Acute Gastroenteritis

Metagenomics based on the next-generation sequencing (NGS) technique is a target-independent assay that enables the simultaneous detection and genomic characterization of all viruses present in a sample. There is a limited amount of data about the virome of individuals with gastroenteritis (GI). In this study, the enteric virome of 250 individuals (92% were children under 5 years old) with GI living in the northeastern and northern regions of Brazil was characterized. Fecal samples were subjected to NGS, and the metagenomic analysis of virus-like particles (VLPs) identified 11 viral DNA families and 12 viral RNA families. As expected, the highest percentage of viral sequences detected were those commonly associated with GI, including rotavirus, adenovirus, norovirus (94.8%, 82% and 71.2%, respectively). The most common co-occurrences, in a single individual, were the combinations of rotavirus-adenovirus, rotavirus-norovirus, and norovirus-adenovirus (78%, 69%, and 62%, respectively). In the same way, common fecal-emerging human viruses were also detected, such as parechovirus, bocaporvirus, cosavirus, picobirnavirus, cardiovirus, salivirus, and Aichivirus. In addition, viruses that infect plants, nematodes, fungi, protists, animals, and arthropods could be identified. A large number of unclassified viral contigs were also identified. We show that the metagenomics approach is a powerful and promising tool for the detection and characterization of different viruses in clinical GI samples.

Virome of bat-infesting arthropods: highly divergent viruses in different vectors

Bats are reservoirs of important zoonotic viruses like Nipah and SARS viruses. However, whether the blood-sucking arthropods on the body surface of bats also carry these viruses, and the relationship between viruses carried by the blood-sucking arthropods and viruses carried by bats, have not been reported. This study collected 686 blood-sucking arthropods on the body surface of bats from Yunnan Province, China between 2012 and 2015, and they included wingless bat flies, bat flies, ticks, mites, and fleas. The viruses carried by these arthropods were analyzed using meta-transcriptomic approach, and 144 highly diverse positive-sense single-stranded RNA, negative-sense single-stranded RNA, and double-stranded RNA viruses were found, of which 138 were potentially new viruses. These viruses were classified into 14 different virus families or orders, including Bunyavirales, Mononegavirales, Reoviridae, and Picornavirales. Further analyses found that Bunyavirales were the most abundant virus group (84% of total virus RNA) in ticks, whereas narnaviruses were the most abundant (52-92%) in the bat flies and wingless bat flies libraries, followed by solemoviruses (1-29%) and reoviruses (0-43%). These viruses were highly structured based on the arthropod types. It is worth noting that no bat-borne zoonotic viruses were found in the virome of bat-infesting arthropod, seemly not supporting that bat surface arthropods are vectors of zoonotic viruses carried by bats. IMPORTANCE Bats are reservoir of many important viral pathogens. To evaluate whether bat-parasitic blood-sucking arthropods participate in the circulation of these important viruses, it is necessary to conduct unbiased virome studies on these arthropods. We evaluated five types of blood-sucking parasitic arthropods on the surface of bats in Yunnan, China and identified a variety of viruses, some of which had high prevalence and abundance level, although there is limited overlap in virome between distant arthropods. While most of the virome discovered here are potentially arthropod-specific viruses, we identified three possible arboviruses, including one orthobunyavirus and two vesiculoviruses (family Rhabdoviridae), suggesting bat-parasitic arthropods carry viruses with risk of spillage, which warrants further study.

Integrated diversity and shared species analyses of human viromes

Diversity analysis has been performed routinely on microbiomes, including human viromes. Shared species analysis has been conducted only rarely, but it can be a powerful supplement to diversity analysis. In the present study, we conducted integrated diversity and shared species analyses of human viromes by reanalyzing three published datasets of human viromes with more than 250 samples from healthy vs. diseased individuals and/or rural vs. urban individuals. We found significant differences in the virome diversity measured in the Hill numbers between the healthy and diseased individuals, with diseased individuals exhibiting higher virome diversity than healthy individuals, and rural individual exhibiting higher virome diversity than urban individuals. We applied both "read randomization" and "sample randomization" algorithms to perform shared species analysis. With the more conservative sample randomization algorithm, the observed number of shared species was significantly smaller than the expected shared species in 50% (8 of 16) of the comparisons. These results suggest that integrated diversity and shared species analysis can offer more comprehensive insights in comparing human virome samples than standard diversity analysis alone with potentially powerful applications in differentiating the effects of diseases or other meta-factors.

Virome of respiratory secretion from children with unknown etiological acute respiratory disease revealed recombinant human parechovirus and other significant viruses

Using viral metagenomics, viral nucleic acid in 30 respiratory secretion samples collected from children with unknown etiological acute respiratory disease were investigated. Sequences showing similarity to human parainfluenza virus 1, anellovirus, bocavirus, coxsackievirus A4, human parechovirus (HPeV), and alphaflexivirus were recovered from these samples. Complete genomes of one anellovirus, one coxsackievirus A4, three parechoviruses were determined from these libraries. The anellovirus (MW267851) phylogenetically clustered with an unpublished anellovirus (MK212032) from respiratory sample of a Vietnamese patient, forming a separate branch neighboring to strains within the genus Betatorquevirus. The genome of coxsackievirus A4 (MW267852) shares the highest sequence identity of 96.4% to a coxsackievirus A4 (MN964079) which was identified in clinical samples from children with Hand, Foot, and Mouth Disease (HFMD). Two (MW267853 and MW267854) of the three parechoviruses belong to HPeV-1 and the other one (MW267855) belongs to HPeV-6. Recombination analysis indicated that an HPeV-1 (MW267854) identified in this study is a putative recombinant occurred between HPeV-1 and HPeV-3. Whether these viruses have association with specific respiratory disease calls for further investigation.

Spatial heterogeneity analysis of the human virome with Taylor's power law

Spatial heterogeneity is a fundamental characteristic of organisms from viruses to humans. Measuring heterogeneity is challenging, especially for naked-eye invisible viruses, but of obvious importance. For example, spatial heterogeneity of virus distribution may strongly influence infection spreading and outbreaks in the case of pathogenic viruses; the spatial distribution (i.e., the inter-subject heterogeneity) of commensal viruses within/on our bodies can influence the competition, coexistence, and dispersal of viruses within or between our bodies. Taylor's power law (TPL) was first discovered in the 1960s to describe the spatial distributions of plant and/or animal populations, and since then it has been verified by numerous experimental and theoretical studies. Recently, TPL has been extended from population to community level and applied to bacterial communities. Here we report the first comprehensive testing of the TPL fitted to human virome datasets. It was found that the human virome follows the TPL as bacterial communities do. Furthermore, the TPL heterogeneity scaling parameter of human virome is virtually the same as that of the human bacterial microbiome (1.916 vs. 1.926). We postulate that the extreme closeness of human viruses and bacteria in heterogeneity scaling coefficients could be attributed to the fact that most of the viruses that were annotated in this study actually belong to bacteriophages (86% viral OTUs) that "piggyback" on their bacterial hosts, and their distributions are likely host-dependent. The scaling parameter, which measures the inter-subject heterogeneity changes, should be an innate property of human microbiomes including both bacteria and viruses. It is similar to the acceleration coefficient of the gravity (g = 9.8) as specified by Newton's law, which is invariant on the earth. Nevertheless, we caution that our postulation is contingent on an implicit assumption that the proportion of bacteriophages to total virome may not change significantly when more virus species can be identified in future.

Multiple clades of Husavirus in South America revealed by next generation sequencing

Husavirus (HuV) is an unclassified virus of the order Picornavirales that has already been identified worldwide in various locations. The genetic, epidemiological, and pathogenic characteristics are, however, little understood. In children with acute gastroenteritis, this study used next-generation sequencing to recognize unknown sources of viruses. In particular, 251 fecal samples obtained from individuals were sequenced in southern, northeastern, and northern Brazil. all samples were also analyzed using culture methods and parasitological tests to classify other enteric pathogens such as bacteria, parasites, and viruses. 1.9% of the samples tested positive for HuV, for a total of 5 positive children, with a mean age of 2 year, with three males and two females. Detailed molecular characterization of full genomes showed that Brazilian HuVs’ nucleotide divergence is less than 11%. The genetic gap between Brazilian sequences and the closest HuV reported previously, on the other hand, is 18%. The study showed that Brazilian sequences are closely related to the HuV defined in Viet Nam in 2013, further characterization based on phylogenetics. At least two divergent clades of HuV in South America were also seen in the phylogenetic study.

Human Protoparvovirus DNA and IgG in Children and Adults with and without Respiratory or Gastrointestinal Infections

Three human protoparvoviruses, bufavirus (BuV), tusavirus (TuV) and cutavirus (CuV), have recently been discovered in diarrheal stool. BuV has been associated with diarrhea and CuV with cutaneous T-cell lymphoma, but there are hardly any data for TuV or CuV in stool or respiratory samples. Hence, using qPCR and IgG enzyme immunoassays, we analyzed 1072 stool, 316 respiratory and 445 serum or plasma samples from 1098 patients with and without gastroenteritis (GE) or respiratory-tract infections (RTI) from Finland, Latvia and Malawi. The overall CuV-DNA prevalences in stool samples ranged between 0-6.1% among our six patient cohorts. In Finland, CuV DNA was significantly more prevalent in GE patients above rather than below 60 years of age (5.1% vs 0.2%). CuV DNA was more prevalent in stools among Latvian and Malawian children compared with Finnish children. In 10/11 CuV DNA-positive adults and 4/6 CuV DNA-positive children with GE, no known causal pathogens were detected. Interestingly, for the first time, CuV DNA was observed in two nasopharyngeal aspirates from children with RTI and the rare TuV in diarrheal stools of two adults. Our results provide new insights on the occurrence of human protoparvoviruses in GE and RTI in different countries.

Persistent Detection of Cosavirus and Saffold Cardiovirus in Riachuelo River, Argentina

Cosaviruses (CoSV) and Saffold cardiovirus (SAFV) are novel members of the Picornaviridae family. The Matanza-Riachuelo river basin covers a total area of 2200 km² with approximately 60 km long. Its last section is called Riachuelo River. The aim of this study was to describe the circulation of both picornaviruses and their relationship with the environmental situation of the Riachuelo River using 274 samples collected from 2005 to 2015. CoSV and SAFV were investigated in samples available by two periods: 2005-2006 and 2014-2015 (103 and 101, respectively). Physicochemical and bacteriological parameters confirmed very high levels of human fecal contamination during the 11 years evaluated. CoSV was detected in 85.7% (66/77) and 65.4% (17/26) of the samples collected in 2005-2006 and 2014-2015 periods, respectively. Species A and D were identified, the first one being widely predominant: 74.1% (20/27) and 75.0% (3/4) in both periods. SAFV virus was detected in 47.1% (32/68) and 52.6% (10/19) in periods 2005-2006 and 2014-2015, respectively. SAFV-6 was the most identified genotype in the entire study, while SAFV-3 was predominant in 2005-2006. The contribution of genotypes 1, 2, 4 and 8 was minor. The high prevalence of CoSV and SAFV suggests that both viruses have been circulating in Argentina at least since 2005. Our results show that a watercourse with high rates of human fecal contamination can become a persistent source of new viruses which capacity to produce human diseases is unknown.

Preliminary assessment of viral metagenome from cancer tissue and blood from patients with lung adenocarcinoma

In this study, using a viral metagenomic method, we investigated the composition of virome in blood and cancer tissue samples that were collected from 25 patients with lung adenocarcinoma. Results indicated that virus sequences showing similarity to human pegivirus (HPgV), anellovirus, human endogenous retrovirus (HERV), and polyomavirus were recovered from this cohort. Three different complete genomes of HPgV were acquired from the blood samples and one complete genome of polyomavirus was determined from the cancer tissue sample. Phylogenetic analysis indicated that the three HPgV strains belonged to genotype 3 and the polyomavirus showed the highest sequence identity (99.73%) to trichodysplasia spinulosa-associated polyomavirus. PCR screening results indicated that the three HPgVs were present in 5 out of the 25 blood samples and the polyomavirus only existed in a cancer tissue sample pool. Whether infections with viruses have an association with lung cancer needs further study with a larger size of sampling.

Virome of Bat Guano from Nine Northern California Roosts

Bats are hosts to a large variety of viruses, including many capable of cross-species transmissions to other mammals, including humans. We characterized the virome in guano from five common bat species in 9 Northern California roosts and from a pool of 5 individual bats. Genomes belonging to 14 viral families known to infect mammals and 17 viral families infecting insects or of unknown tropism were detected. Nearly complete or complete genomes of a novel parvovirus, astrovirus, nodavirus, circular Rep-encoding single-stranded DNA (CRESS-DNA) viruses, and densoviruses, and more partial genomes of a novel alphacoronavirus and a bunyavirus were characterized. Lower numbers of reads with >90% amino acid identity to previously described calicivirus, circovirus, adenoviruses, hepatovirus, bocaparvoviruses, and polyomavirus in other bat species were also found, likely reflecting their wide distribution among different bats. Unexpectedly, a few sequence reads of canine parvovirus 2 and the recently described mouse kidney parvovirus were also detected and their presence confirmed by PCR; these possibly originated from guano contamination by carnivores and rodents. The majority of eukaryotic viral reads were highly divergent, indicating that numerous viruses still remain to be characterized, even from such a heavily investigated order as Chiroptera.IMPORTANCE Characterizing the bat virome is important for understanding viral diversity and detecting viral spillover between animal species. Using an unbiased metagenomics method, we characterize the virome in guano collected from multiple roosts of common Northern California bat species. We describe several novel viral genomes and report the detection of viruses with close relatives reported in other bat species, likely reflecting cross-species transmissions. Viral sequences from well-known carnivore and rodent parvoviruses were also detected, whose presence are likely the result of contamination from defecation and urination atop guano and which reflect the close interaction of these mammals in the wild.

CrAssphage and its bacterial host in cat feces

CrAssphages are a diverse group of related phages detected in human feces where they are the most prevalent and abundant prokaryotic virus. CrAssphages' cellular host has been identified as the anaerobic Bacteroides intestinalis. CrAssphage has also been reported in non-human primates and environmental samples and has been proposed as a marker of human fecal contamination. Here we describe crAssphage DNA in a feline fecal sample. 95% of the ~ 100 Kb genome could be assembled and classified in genus 1 of the recently proposed Alphacrassvirinae subfamily. The cat origin of the fecal sample was confirmed by partial mitochondrial DNA sequencing. High levels of Bacteroides intestinalis DNA could also be detected in this cat's feces. Fecal samples longitudinally collected over a 4-week period showed the continuous shedding of crAssphage DNA. We therefore report the first genome sequence-confirmed detection of crAssphage in fecal samples of a non-primate mammal.

Metagenomic sequencing of stool samples in Bangladeshi infants: virome association with poliovirus shedding after oral poliovirus vaccination

The potential role of enteric viral infections and the developing infant virome in affecting immune responses to the oral poliovirus vaccine (OPV) is unknown. Here we performed viral metagenomic sequencing on 3 serially collected stool samples from 30 Bangladeshi infants following OPV vaccination and compared findings to stool samples from 16 age-matched infants in the United States (US). In 14 Bangladeshi infants, available post-vaccination serum samples were tested for polio-neutralizing antibodies. The abundance (p = 0.006) and richness (p = 0.013) of the eukaryotic virome increased with age and were higher than seen in age-matched US infants (p < 0.001). In contrast, phage diversity metrics remained stable and were similar to those in US infants. Non-poliovirus eukaryotic virus abundance (3.68 log10 vs. 2.25 log10, p = 0.002), particularly from potential viral pathogens (2.78log10 vs. 0.83log10, p = 0.002), and richness (p = 0.016) were inversely associated with poliovirus shedding. Following vaccination, 28.6% of 14 infants tested developed neutralizing antibodies to all three Sabin types and also exhibited higher rates of poliovirus shedding (p = 0.020). No vaccine-derived poliovirus variants were detected. These results reveal an inverse association between eukaryotic virome abundance and poliovirus shedding. Overall gut virome ecology and concurrent viral infections may impact oral vaccine responsiveness in Bangladeshi infants.

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.

The Husavirus Posa-Like Viruses in China, and a New Group of Picornavirales

Novel posa-like viral genomes were first identified in swine fecal samples using metagenomics and were designated as unclassified viruses in the order Picornavirales. In the present study, nine husavirus strains were identified in China. Their genomes share 94.1–99.9% similarity, and alignment of these nine husavirus strains identified 697 nucleotide polymorphism sites across their full-length genomes. These nine strains were directly clustered with the Husavirus 1 lineage, and their genomic arrangement showed similar characteristics. These posa-like viruses have undergone a complex evolutionary process, and have a wide geographic distribution, complex host spectrum, deep phylogenetic divergence, and diverse genomic organizations. The clade of posa-like viruses forms a single group, which is evolutionarily distinct from other known families and could represent a distinct family within the Picornavirales. The genomic arrangement of Picornavirales and the new posa-like viruses are different, whereas the posa-like viruses have genomic modules similar to the families Dicistroviridae and Marnaviridae. The present study provides valuable genetic evidence of husaviruses in China, and clarifies the phylogenetic dynamics and the evolutionary characteristics of Picornavirales.

Recent advances in understanding the neonatal microbiome

The neonatal developmental window represents a key time for establishment of the gut microbiota. First contact with these microbes within the infant gastrointestinal tract signifies the start of a critical mutualistic relationship, which is central for short- and longer-term health. Recent research has provided insights into the origin of these microbial pioneers, how they are maintained within the gut environment, and how factors such as antibiotics or preterm birth may disrupt the succession of beneficial microbes.  The acquisition, colonisation, and maintenance of the early life microbiota, and subsequent interactions with the host is a rapidly developing research area. In this review we explore some of these key topics which have been illuminated by recent research, and we highlight some of the important unresolved questions which currently limit our overall understanding of the neonatal gut microbiome.

Semen virome of men with HIV on or off antiretroviral treatment

OBJECTIVES

Improving immune status of people living with HIV through antiretroviral therapy (ART) may also reduce shedding of other viruses in semen. We characterized the seminal fluid virome of men with HIV and tested potential associations between viruses present and CD4 T-cell count, HIV viremia, and antiretroviral therapy (ART) status.

DESIGN AND METHODS

Metagenomics was used to enrich and sequence viral nucleic acids from the seminal fluid of 55 semen samples from 42 men living with HIV from San Francisco with a median age of 33 (IQR, 28.7-45) and median CD4 T-cell counts of 837 cells/μl (IQR, 258-1571 cells/μl). All samples were collected between 2005 and 2015, and ART status was ascertained from medical records.

RESULTS

Anelloviruses, cytomegalovirus (CMV), and multiple genotypes of human papillomaviruses were detected. Participants shed from 0 to 4 distinct human viruses. Longitudinally collected seminal fluid samples showed changes in the viruses shed. Viruses were more frequently shed by individuals with detectable HIV viremia (43.7 vs. 15.4%, P = 0.042). A trend was seen for increased shedding by individuals who were not on ART (42.8 vs. 17.8%, P = 0.082) or with CD4 T-cell count less than 350 cells/μl (35.3 vs. 20%, P = 0.27).

CONCLUSION

Seminal fluid from men with HIV from San Francisco contains nucleic acids from three different DNA viral families. A greater number of viruses, particularly CMV, were shed by participants with detectable HIV viremia (18.9 vs. 0%, P = 0.022). Control of viremia through ART may lower shedding of other viruses in semen in addition to HIV.

Identification and molecular characterization of the first complete genome sequence of Human Parechovirus type 15

Using a metagenomics approach, we have determined the first full-length genome sequence of a human parechovirus type 15 (HPeV15) strain, isolated from a child with acute flaccid paralysis and co-infected with EV-A71. HPeV15 is a rarely reported type. To date, no full-length genome sequence of HPeV15 is available in the GenBank database, where only limited VP1 sequences of this virus are available. Pairwise comparisons of the complete VP1 nucleotide and deduced amino acid sequences revealed that the study strain belongs to type 15 as it displayed 79.6% nucleotide and 93.4% amino acid identity with the HPeV15 prototype strain. Comparative analysis of available genomic regions and phylogenetic analysis using the P2 and P3 coding regions revealed low nucleotide identity to HPeV reference genomes. Phylogenetic and similarity plot analyses showed that genomic recombination events might have occurred in the UTRs and nonstructural region during HPeV15 evolution. The study strain has high similarity features with different variants of HPeV3 suggesting intertypic recombination. Our data contributes to the scarce data available on HPeVs in Africa and provides valuable information for future studies that aim to understand the evolutionary history, molecular epidemiology or biological and pathogenic properties of HPeV15.

Discovery and Prevalence of Divergent RNA Viruses in European Field Voles and Rabbits

The advent of unbiased metagenomic virus discovery has revolutionized studies of virus biodiversity and evolution. Despite this, our knowledge of the virosphere, including in mammalian species, remains limited. We used unbiased metagenomic sequencing to identify RNA viruses in European field voles and rabbits. Accordingly, we identified a number of novel RNA viruses including astrovirus, rotavirus A, picorna-like virus and a morbilli-like paramyxovirus. In addition, we identified a sobemovirus and a novel luteovirus that likely originated from the rabbit diet. These newly discovered viruses were often divergent from those previously described. The novel astrovirus was most closely related to a virus sampled from the rodent-eating European roller bird (Coracias garrulous). PCR screening revealed that the novel morbilli-like paramyxovirus in the UK field vole had a prevalence of approximately 4%, and shared common ancestry with other rodent morbilli-like viruses sampled globally. Two novel rotavirus A sequences were detected in a UK field vole and a French rabbit, the latter with a prevalence of 5%. Finally, a highly divergent picorna-like virus found in the gut of the French rabbit virus was only ~35% similar to an arilivirus at the amino acid level, suggesting the presence of a novel viral genus within the Picornaviridae.

Molecular-based investigation and genetic characterization of porcine stool-associated RNA virus (posavirus) lineages 1 to 3 in pigs in South Korea from 2017 to 2019

Recent results on the detection and genetic characterization of stool-associated RNA viruses from different species have increased the knowledge about the extreme genetic diversity of picornaviruses. This study aimed to investigate the presence of unclassified porcine stool-associated RNA viruses (posaviruses) in South Korea and to elucidate the molecular evolution of the viruses. By RT-PCR, posaviruses 1 and 3 were exclusively found in fecal samples and consistently detected in three consecutive years in six of eight provinces, with 148/697 (21.2%) and 33/84 (39.3%) positive samples and farms, respectively. Every age group but the older age groups (finisher, sow) had significantly higher positive rates of posavirus 1 than posavirus 3. An analysis of the RNA-dependent RNA polymerase sequences by likelihood mapping and maximum-likelihood-based phylogenetic analysis revealed that stool-associated RNA viruses formed four supergroups that were well separated from all recognized families of the order Picornavirales. Five genomes of Korean posaviruses generated in this study were phylogenetically grouped with posavirus 1 and posavirus 3 and were predicted to have the typical genome organization of picornaviruses.

The Coxsackievirus and Adenovirus Receptor, a Required Host Factor for Recovirus Infection, Is a Putative Enteric Calicivirus Receptor

Human norovirus (HuNoV) is a leading cause of acute gastroenteritis in both developed and developing countries. Studies of HuNoV host cell interactions are limited by the lack of a simple, robust cell culture system. Due to their diverse HuNoV-like biological features, including histo-blood group antigen (HBGA) binding, rhesus enteric caliciviruses (ReCVs) are viable surrogate models for HuNoVs. In addition, several ReCV strains can be propagated to high titers in standard nonhuman primate cell lines while causing lytic infection and cell death. To identify the ReCV entry receptor, we performed CRISPR/Cas9 library screening in Vero cells, which identified the coxsackievirus and adenovirus receptor (CAR) as a candidate ReCV entry receptor. We showed that short interfering RNA, anti-human CAR (hCAR) monoclonal antibody RmcB treatment, and recombinant hCAR ectodomain blocked ReCV replication in LLC-MK2 cells. CRISPR/Cas9-targeted knockout of CAR in LLC-MK2 and Vero cells made these cell lines resistant to ReCV infection, and susceptibility to infection could be restored by transient expression of CAR. CHO cells do not express CAR or HBGAs and are resistant to ReCV infection. Recombinant CHO cells stably expressing hCAR or the type B HBGA alone did not support ReCV infection. However, CHO cells expressing both hCAR and the type B HBGA were susceptible to ReCV infection. In summary, we have demonstrated that CAR is required for ReCV infection and most likely is a functional ReCV receptor, but HBGAs are also necessary for infection.IMPORTANCE Because of the lack of a simple and robust human norovirus (HuNoV) cell culture system surrogate, caliciviruses still represent valuable research tools for norovirus research. Due to their remarkable biological similarities to HuNoVs, including the utilization of HBGAs as putative attachment receptors, we used rhesus enteric caliciviruses (ReCVs) to study enteric calicivirus host cell interactions. Using CRISPR/Cas9 library screening and functional assays, we identified and validated the coxsackievirus and adenovirus receptor (CAR) as a functional proteinaceous receptor for ReCVs. Our work demonstrated that CAR and HBGAs both are necessary to convert a nonsusceptible cell line to being susceptible to ReCV infection. Follow-up studies to evaluate the involvement of CAR in HuNoV infections are ongoing.

High Diversity and Novel Enteric Viruses in Fecal Viromes of Healthy Wild and Captive Thai Cynomolgus Macaques (Macaca fascicularis)

Cynomolgus macaques are common across South East Asian countries including Thailand. The National Primate Research Center of Thailand, Chulalongkorn University (NPRCT-CU) captures wild-borne cynomolgus macaque for research use. Limited information is available on the enteric viruses and possible zoonotic infections into or from cynomolgus macaques. We characterized and compare the fecal virome of two populations; healthy wild-originated captive cynomolgus macaques (n = 43) reared in NPRCT-CU and healthy wild cynomolgus macaques (n = 35). Over 90% of recognized viral sequence reads amplified from feces were from bacterial viruses. Viruses from seven families of mammalian viruses were also detected (Parvoviridae, Anelloviridae, Picornaviridae, Adenoviridae, Papillomaviridae, Herpesviridae, and Caliciviridae). The genomes of a member of a new picornavirus genus we named Mafapivirus, a primate chapparvovirus, and a circular Rep-encoding single-strand (CRESS) DNA virus were also characterized. Higher abundance of CRESS DNA viruses of unknown tropism and invertebrate-tropic ambidensovirus were detected in wild versus captive macaques likely reflecting dietary differences. Short term rearing in captivity did not have a pronounced effect on the diversity of mammalian viruses of wild cynomolgus macaques. This study is the first report of the fecal virome of cynomolgus macaques, non-human primates frequently used in biomedical research and vaccination studies.

Molecular detection of enteric viruses in the stool samples of children without diarrhea in Bangladesh

A number of molecular epidemiological studies reported the detection of enteric viruses in asymptomatic children. The role of these viruses in an asymptomatic infection remains unclear. This study investigated the enteric viruses in the stool samples collected from children without diarrhea. Stool samples were collected during June to October 2016, from 227 children who lived in Matlab, Bangladesh. Seventeen enteric viruses, including rotavirus A, B, and C (RVA, RVB, and RVC), norovirus GI (NoV GI), norovirus GII (NoV GII), sapovirus (SaV), adenovirus (AdV), human astrovirus (HAstV), Aichivirus (AiV), human parechovirus (HPeV), enterovirus (EV), human bocavirus (HBoV), Saffold virus (SAFV), human cosavirus (HCoSV), bufavirus (BufV), salivirus (SalV), and rosavirus (RoV), were investigated by RT-PCR method. One hundred and eighty-two (80.2%; 182/227) samples were positive for some of these viruses, and 19.8% (45/227) were negative. Among the positive samples, 46.7% (85/182) were a single infection, and 53.3% (97/182) were coinfection with multiple viruses. The HCoSV was the most prevalent virus (41.4%), followed by EV (32.2%), NoV GII (25.6%), HPeV (8.8%), RVA (6.2%), AdV (5.7%), AiV (5.3%), SAFV (4.4%), and SaV (2.6%). Each of NoV GI, HAstV, HBoV, and BufV was detected at 0.4%. However, RVB, RVC, SalV, and RoV were not detected in this study. Phylogenetic analysis showed that diverse HCoSV species and genotypes were circulating in Bangladesh, and four strains of species A are proposed to be new genotypes. The data indicated that non-diarrheal Bangladeshi children were asymptomatically infected with wide varieties of enteric viruses.

Effect of Geographic Isolation on the Nasal Virome of Indigenous Children

The influence of living in small remote villages on the diversity of viruses in the nasal mucosa was investigated in three Colombian villages with very different levels of geographic isolation. Metagenomic analysis was used to characterize viral nucleic acids in nasal swabs from 63 apparently healthy young children. Sequences from human virus members of the families Anelloviridae, Papillomaviridae, Picornaviridae, Herpesviridae, Polyomaviridae, Adenoviridae, and Paramyxoviridae were detected in decreasing proportions of children. The number of papillomavirus infections detected was greater among Hispanic children most exposed to outside contacts, while anellovirus infections were more common in the isolated indigenous villages. The diversity of the other human viruses detected did not differ among the villages. Closely related variants of rhinovirus A or B were identified in 2 to 4 children from each village, reflecting ongoing transmission clusters. Genomes of viruses not currently known to infect humans, including members of the families Parvoviridae, Partitiviridae, Dicistroviridae, and Iflaviridae and circular Rep-encoding single-stranded DNA (CRESS-DNA) virus, were also detected in nasal swabs, possibly reflecting environmental contamination from insect, fungal, or unknown sources. Despite the high levels of geographic and cultural isolation, the overall diversity of human viruses in the nasal passages of children was not reduced in highly isolated indigenous villages, indicating ongoing exposure to globally circulating viruses.IMPORTANCE Extreme geographic and cultural isolation can still be found in some indigenous South American villages. Such isolation may be expected to limit the introduction of otherwise common and widely distributed viruses. Very small population sizes may also result in rapid local viral extinction due to a lack of seronegative subjects to maintain transmission chains for rapidly cleared viruses. We compared the viruses in the nasal passages of young children in three villages with increasing levels of geographic isolation. We found that isolation did not reduce the overall diversity of viral infections. Multiple infections with nearly identical rhinoviruses could be detected within each village, likely reflecting recent viral introductions and transmission clusters among epidemiologically linked members of these very small communities. We conclude that, despite their geographic isolation, remote indigenous villages show evidence of ongoing exposure to globally circulating viruses.

Altered respiratory virome and serum cytokine profile associated with recurrent respiratory tract infections in children

Recurrent acute respiratory tract infections (ARTIs) affect a large population, yet the specific decisive factors are largely unknown. Here we study a population of 4407 children diagnosed with ARTI, comparing respiratory virome and serum cytokine profiles associated with multiple ARTIs and single ARTI during a six-year period. The relative abundance of Propionibacterium phages is significantly elevated in multiple ARTIs compared to single ARTI group. Serum levels of TIMP-1 and PDGF-BB are markedly increased in multiple ARTIs compared to single-ARTI and non-ARTI controls, making these two cytokines potential predictors for multiple ARTIs. The presence of Propionibacterium phages is associated with higher levels of TIMP-1 and PDGF-BB. Receiver operating characteristic (ROC) curve analyses show that the combination of TIMP-1, PDGF-BB and Propionibacterium phages could be a strong predictor for multiple ARTIs. These findings indicate that respiratory microbe homeostasis and specific cytokines are associated with the onset of multiple ARTIs over time.

Phylogenetic analysis of novel posaviruses detected in feces of Japanese pigs with posaviruses and posa-like viruses of vertebrates and invertebrates

Posaviruses and posa-like viruses are unclassified viruses with sequence similarity to viruses of the order Picornavirales. They have been reported in various vertebrates and invertebrates. We identified 11 posavirus-like sequences in porcine feces and performed phylogenic analysis. Previously reported Japanese posaviruses and those identified in this study clustered with posavirus 1, 4, and 7 and husavirus 1, while five viruses branched into three independent lineages, tentatively named posavirus 10, 11, and 12. Interestingly, posaviruses, except for posavirus 8 and 9, husaviruses, and rasaviruses, formed a cluster consisting of viruses only from pigs, humans, and rats, while posavirus 8 and 9, fisavirus, and basaviruses clustered with posa-like viruses from invertebrates.

Caliciviridae Other Than Noroviruses

Besides noroviruses, the Caliciviridae family comprises four other accepted genera: Sapovirus, Lagovirus, Vesivirus, and Nebovirus. There are six new genera proposed: Recovirus, Valovirus, Bavovirus, Nacovirus, Minovirus, and Salovirus. All Caliciviridae have closely related genome structures, but are genetically and antigenically highly diverse and infect a wide range of mammalian host species including humans. Recombination in nature is not infrequent for most of the Caliciviridae, contributing to their diversity. Sapovirus infections cause diarrhoea in pigs, humans and other mammalian hosts. Lagovirus infections cause systemic haemorrhagic disease in rabbits and hares, and vesivirus infections lead to lung disease in cats, vesicular disease in swine, and exanthema and diseases of the reproductive system in large sea mammals. Neboviruses are an enteric pathogen of cattle, differing from bovine norovirus. At present, only a few selected caliciviruses can be propagated in cell culture (permanent cell lines or enteroids), and for most of the cultivatable caliciviruses helper virus-free, plasmid only-based reverse genetics systems have been established. The replication cycles of the caliciviruses are similar as far as they have been explored: viruses interact with a multitude of cell surface attachment factors (glycans) and co-receptors (proteins) for adsorption and penetration, use cellular membranes for the formation of replication complexes and have developed mechanisms to circumvent innate immune responses. Vaccines have been developed against lagoviruses and vesiviruses, and are under development against human noroviruses.

Do you cov me? Effect of coverage reduction on species identification and genome reconstruction in complex biological matrices by metagenome shotgun high-throughput sequencing

Shotgun metagenomics sequencing is a powerful tool for the characterization of complex biological matrices, enabling analysis of prokaryotic and eukaryotic organisms and viruses in a single experiment, with the possibility of reconstructing de novo the whole metagenome or a set of genes of interest. One of the main factors limiting the use of shotgun metagenomics on wide scale projects is the high cost associated with the approach. However, we demonstrate that-for some applications-it is possible to use shallow shotgun metagenomics to characterize complex biological matrices while reducing costs. We measured the variation of several summary statistics simulating a decrease in sequencing depth by randomly subsampling a number of reads. The main statistics that were compared are alpha diversity estimates, species abundance, detection threshold, and ability of reconstructing the metagenome in terms of length and completeness. Our results show that a classification of prokaryotic, eukaryotic and viral communities can be accurately performed even using very low number of reads, both in mock communities and in real complex matrices. With samples of 100,000 reads, the alpha diversity estimates were in most cases comparable to those obtained with the full sample, and the estimation of the abundance of all the present species was in excellent agreement with those obtained with the full sample. On the contrary, any task involving the reconstruction of the metagenome performed poorly, even with the largest simulated subsample (1M reads). The length of the reconstructed assembly was smaller than the length obtained with the full dataset, and the proportion of conserved genes that were identified in the meta-genome was drastically reduced compared to the full sample. Shallow shotgun metagenomics can be a useful tool to describe the structure of complex matrices, but it is not adequate to reconstruct-even partially-the metagenome.

Do you cov me? Effect of coverage reduction on species identification and genome reconstruction in complex biological matrices by metagenome shotgun high-throughput sequencing

Shotgun metagenomics sequencing is a powerful tool for the characterization of complex biological matrices, enabling analysis of prokaryotic and eukaryotic organisms and viruses in a single experiment, with the possibility of reconstructing de novo the whole metagenome or a set of genes of interest. One of the main factors limiting the use of shotgun metagenomics on wide scale projects is the high cost associated with the approach. We set out to determine  if it is possible to use shallow shotgun metagenomics to characterize complex biological matrices while reducing costs. We measured the variation of several summary statistics simulating a decrease in sequencing depth by randomly subsampling a number of reads. The main statistics that were compared are alpha diversity estimates, species abundance, and ability of reconstructing de novo the metagenome in terms of length and completeness. Our results show that diversity indices of complex prokaryotic, eukaryotic and viral communities can be accurately estimated with 500,000 reads or less, although particularly complex samples may require 1,000,000 reads. On the contrary, any task involving the reconstruction of the metagenome performed poorly, even with the largest simulated subsample (1,000,000 reads). The length of the reconstructed assembly was smaller than the length obtained with the full dataset, and the proportion of conserved genes that were identified in the meta-genome was drastically reduced compared to the full sample. Shallow shotgun metagenomics can be a useful tool to describe the structure of complex matrices, but it is not adequate to reconstruct-even partially-the metagenome.

Do you cov me? Effect of coverage reduction on species identification and genome reconstruction in complex biological matrices by metagenome shotgun high-throughput sequencing

Shotgun metagenomics sequencing is a powerful tool for the characterization of complex biological matrices, enabling analysis of prokaryotic and eukaryotic organisms in a single experiment, with the possibility of de novo reconstruction of the whole metagenome or a set of genes of interest. One of the main factors limiting the use of shotgun metagenomics on wide scale projects is the high cost associated with the approach. However, we demonstrate that-for some applications-it is possible to use shallow shotgun metagenomics to characterize complex biological matrices while reducing costs. Here we compared the results obtained on full size, real datasets with results obtained by randomly extracting a fixed number of reads. The main statistics that were compared are alpha diversity estimates, species abundance, and ability of reconstructing the metagenome in terms of length and completeness. Our results show that a classification of the communities present in a complex matrix can be accurately performed even using very low number of reads. With samples of 100,000 reads, the alpha diversity estimates were in most cases comparable to those obtained with the full sample, and the estimation of the abundance of all the present species was in excellent agreement with those obtained with the full sample. On the contrary, any task involving the reconstruction of the metagenome performed poorly, even with the largest simulated subsample (1M reads). The length of the reconstructed assembly was sensibly smaller than the length obtained with the full dataset, and the proportion of conserved genes that were identified in the meta-genome was drastically reduced compared to the full sample. Shallow shotgun metagenomics can be a useful tool to describe the structure of complex matrices, but it is not adequate to reconstruct de novo-even partially-the metagenome.

Do you cov me? Effect of coverage reduction on metagenome shotgun sequencing studies

Shotgun metagenomics sequencing is a powerful tool for the characterization of complex biological matrices, enabling analysis of prokaryotic and eukaryotic organisms and viruses in a single experiment, with the possibility of reconstructing de novo the whole metagenome or a set of genes of interest. One of the main factors limiting the use of shotgun metagenomics on wide scale projects is the high cost associated with the approach. We set out to determine if it is possible to use shallow shotgun metagenomics to characterize complex biological matrices while reducing costs. We used a staggered mock community to estimate the optimal threshold for species detection. We measured the variation of several summary statistics simulating a decrease in sequencing depth by randomly subsampling a number of reads. The main statistics that were compared are diversity estimates, species abundance, and ability of reconstructing de novo the metagenome in terms of length and completeness. Our results show that diversity indices of complex prokaryotic, eukaryotic and viral communities can be accurately estimated with 500,000 reads or less, although particularly complex samples may require 1,000,000 reads. On the contrary, any task involving the reconstruction of the metagenome performed poorly, even with the largest simulated subsample (1,000,000 reads). The length of the reconstructed assembly was smaller than the length obtained with the full dataset, and the proportion of conserved genes that were identified in the meta-genome was drastically reduced compared to the full sample. Shallow shotgun metagenomics can be a useful tool to describe the structure of complex matrices, but it is not adequate to reconstruct-even partially-the metagenome.