Transmission and dynamics of mother-infant gut viruses during pregnancy and early life

Early development of the gut ecosystem is crucial for lifelong health. While infant gut bacterial communities have been studied extensively, the infant gut virome remains under-explored. To study the development of the infant gut virome over time and the factors that shape it, we longitudinally assess the composition of gut viruses and their bacterial hosts in 30 women during and after pregnancy and in their 32 infants during their first year of life. Using shotgun metagenomic sequencing applied to dsDNA extracted from Virus-Like Particles (VLPs) and bacteria, we generate 205 VLP metaviromes and 322 total metagenomes. With this data, we show that while the maternal gut virome composition remains stable during late pregnancy and after birth, the infant gut virome is dynamic in the first year of life. Notably, infant gut viromes contain a higher abundance of active temperate phages compared to maternal gut viromes, which decreases over the first year of life. Moreover, we show that the feeding mode and place of delivery influence the gut virome composition of infants. Lastly, we provide evidence of co-transmission of viral and bacterial strains from mothers to infants, demonstrating that infants acquire some of their virome from their mother's gut.

Discovery, diversity, and functional associations of crAss-like phages in human gut metagenomes from four Dutch cohorts

The crAss-like phages are a diverse group of related viruses that includes some of the most abundant viruses of the human gut. To explore their diversity and functional role in human population and clinical cohorts, we analyze gut metagenomic data collected from 1,950 individuals from the Netherlands. We identify 1,556 crAss-like phage genomes, including 125 species-level and 32 genus-level clusters absent from the reference databases used. Analysis of their genomic features shows that closely related crAss-like phages can possess strikingly divergent regions responsible for transcription, presumably acquired through recombination. Prediction of crAss-like phage hosts points primarily to bacteria of the phylum Bacteroidetes, consistent with previous reports. Finally, we explore the temporal stability of crAss-like phages over a 4-year period and identify associations between the abundance of crAss-like phages and several human phenotypes, including depletion of crAss-like phages in inflammatory bowel disease patients.

Stability of the human gut virome and effect of gluten-free diet

The human gut microbiome consists of bacteria, archaea, eukaryotes, and viruses. The gut viruses are relatively underexplored. Here, we longitudinally analyzed the gut virome composition in 11 healthy adults: its stability, variation, and the effect of a gluten-free diet. Using viral enrichment and a de novo assembly-based approach, we demonstrate the quantitative dynamics of the gut virome, including dsDNA, ssDNA, dsRNA, and ssRNA viruses. We observe highly divergent individual viral communities, carrying on an average 2,143 viral genomes, 13.1% of which were present at all 3 time points. In contrast to previous reports, the Siphoviridae family dominates over Microviridae in studied individual viromes. We also show individual viromes to be stable at the family level but to vary substantially at the genera and species levels. Finally, we demonstrate that lower initial diversity of the human gut virome leads to a more pronounced effect of the dietary intervention on its composition.

Discovery of an essential nucleotidylating activity associated with a newly delineated conserved domain in the RNA polymerase-containing protein of all nidoviruses

RNA viruses encode an RNA-dependent RNA polymerase (RdRp) that catalyzes the synthesis of their RNA(s). In the case of positive-stranded RNA viruses belonging to the order Nidovirales, the RdRp resides in a replicase subunit that is unusually large. Bioinformatics analysis of this non-structural protein has now revealed a nidoviral signature domain (genetic marker) that is N-terminally adjacent to the RdRp and has no apparent homologs elsewhere. Based on its conservation profile, this domain is proposed to have nucleotidylation activity. We used recombinant non-structural protein 9 of the arterivirus equine arteritis virus (EAV) and different biochemical assays, including irreversible labeling with a GTP analog followed by a proteomics analysis, to demonstrate the manganese-dependent covalent binding of guanosine and uridine phosphates to a lysine/histidine residue. Most likely this was the invariant lysine of the newly identified domain, named nidovirus RdRp-associated nucleotidyltransferase (NiRAN), whose substitution with alanine severely diminished the described binding. Furthermore, this mutation crippled EAV and prevented the replication of severe acute respiratory syndrome coronavirus (SARS-CoV) in cell culture, indicating that NiRAN is essential for nidoviruses. Potential functions supported by NiRAN may include nucleic acid ligation, mRNA capping and protein-primed RNA synthesis, possibilities that remain to be explored in future studies.

Arterivirus nsp12 versus the coronavirus nsp16 2'-O-methyltransferase: comparison of the C-terminal cleavage products of two nidovirus pp1ab polyproteins

The 3'-terminal domain of the most conserved ORF1b in three of the four families of the order Nidovirales (except for the family Arteriviridae) encodes a (putative) 2'-O-methyltransferase (2'-O-MTase), known as non structural protein (nsp) 16 in the family Coronaviridae and implicated in methylation of the 5' cap structure of nidoviral mRNAs. As with coronavirus transcripts, arterivirus mRNAs are assumed to possess a 5' cap although no candidate MTases have been identified thus far. To address this knowledge gap, we analysed the uncharacterized nsp12 of arteriviruses, which occupies the ORF1b position equivalent to that of the nidovirus 2'-O-MTase (coronavirus nsp16). In our in-depth bioinformatics analysis of nsp12, the protein was confirmed to be family specific whilst having diverged much further than other nidovirus ORF1b-encoded proteins, including those of the family Coronaviridae. Only one invariant and several partially conserved, predominantly aromatic residues were identified in nsp12, which may adopt a structure with alternating α-helices and β-strands, an organization also found in known MTases. However, no statistically significant similarity was found between nsp12 and the twofold larger coronavirus nsp16, nor could we detect MTase activity in biochemical assays using recombinant equine arteritis virus (EAV) nsp12. Our further analysis established that this subunit is essential for replication of this prototypic arterivirus. Using reverse genetics, we assessed the impact of 25 substitutions at 14 positions, yielding virus phenotypes ranging from WT-like to non-viable. Notably, replacement of the invariant phenylalanine 109 with tyrosine was lethal. We concluded that nsp12 plays an essential role during EAV replication, possibly by acting as a co-factor for another enzyme.