crAssphage abundance and genomic selective pressure correlate with altered bacterial abundance in the fecal microbiota of South African mother-infant dyads
Authors: Brown Bryan P. , Wendoh Jerome , Chopera Denis , Havyarimana Enock , Jaumdally Shameem , Nyangahu Donald D. , Gray Clive M. , Martin Darren P. , Varsani Arvind , Jaspan Heather B. 
Published: March 2019

crAssphages are a class of bacteriophages that are highly abundant in the human gastrointestinal tract. Accordingly, crAssphage genomes have been identified in most human fecal viral metagenome studies. However, we currently have an incomplete understanding of factors impacting the transmission frequencies of these phages between mothers and infants, and the evolutionary pressures associated with such transmissions. Here, we use metagenome sequencing of stool-associated virus-like particles to identify the prevalence of crAssphage across ten South African mother-infant dyads that are discordant for HIV infection. We report the identification of a complete 97kb crAssphage genome, parts of which are detected at variable levels across each mother-infant dyad. We observed average nucleotide sequence identities of >99% for crAssphages from related mother-infant pairs but ∼97% identities between crAssphages from unrelated mothers and infants: a finding strongly suggestive of vertical mother to infant transmission. We further analyzed patterns of nucleotide diversity across the crAssphage sequences described here, identifying particularly elevated positive selection in RNA polymerase and phage tail protein encoding genes, which we validated against a crAssphage genome from previous studies. Using 16S rRNA gene sequencing, we found that the relative abundances ofBacteroides thetaiotaomicronandParabacteroides merdae(Order: Bacteroidales) were differentially correlated with crAssphage abundance. Together, our results reveal that crAssphages may be vertically transmitted from mothers to their infants and that hotspots of selection within crAssphage RNA polymerase and phage tail protein encoding genes are potentially mediated by interactions between crAssphages and their bacterial partners.

Importance

crAssphages are an ubiquitous member of the human gut microbiome and modulate interactions with key bacterial associates within the order Bacteroidales. However, the role of this interaction in the genomic evolution of crAssphage remains unclear. Across a longitudinally sampled cohort of ten South African mother-infant dyads, we use metagenome sequencing of the fecal virome and 16S rRNA gene sequencing of the fecal bacterial microbiota to elucidate the ecological and evolutionary dynamics of these interactions. Here, we demonstrate elevated levels of crAssphage average nucleotide identity between related mother-infant dyads as compared to unrelated individuals, suggesting vertical transmission. We report strong positive selection in crAssphage RNA polymerase and phage tail protein genes. Finally, we demonstrate that crAssphage abundance is linearly correlated (P < 0.014) with the abundance of two bacterial taxa,Bacteroides thetaiotaomicronandParabacteroides merdae.These results suggest that phage-bacterial interactions may help shape ecological and evolutionary dynamics in the gut.