A Phylogenetic Framework to Simulate Synthetic Interspecies RNA-Seq Data.

Interspecies RNA-Seq datasets are increasingly common, and have the potential to answer new questions about the evolution of gene expression. Single-species differential expression analysis is now a well-studied problem that benefits from sound statistical methods. Extensive reviews on biological or synthetic datasets have provided the community with a clear picture on the relative performances of the available methods in various settings. However, synthetic dataset simulation tools are still missing in the interspecies gene expression context. In this work, we develop and implement a new simulation framework. This tool builds on both the RNA-Seq and the phylogenetic comparative methods literatures to generate realistic count datasets, while taking into account the phylogenetic relationships between the samples. We illustrate the usefulness of this new framework through a targeted simulation study, that reproduces the features of a recently published dataset, containing gene expression data in adult eye tissue across blind and sighted freshwater crayfish species. Using our simulated datasets, we perform a fair comparison of several approaches used for differential expression analysis. This benchmark reveals some of the strengths and weaknesses of both the classical and phylogenetic approaches for interspecies differential expression analysis, and allows for a reanalysis of the crayfish dataset. The tool has been integrated in the R package compcodeR, freely available on Bioconductor.

Transgenic expression of Rubisco accumulation factor2 and Rubisco subunits increases photosynthesis and growth in maize.

Carbon assimilation by Rubisco is often a limitation to photosynthesis and therefore plant productivity. We have previously shown that transgenic co-expression of the Rubisco large (LS) and small (SS) subunits along with an essential Rubisco accumulation factor, Raf1, leads to faster growth, increased photosynthesis, and enhanced chilling tolerance in maize (Zea mays). Maize also requires Rubisco accumulation factor2 (Raf2) for full accumulation of Rubisco. Here we have analyzed transgenic maize lines with increased expression of Raf2 or Raf2 plus LS and SS. We show that increasing Raf2 expression alone had minor effects on photosynthesis, whereas expressing Raf2 with Rubisco subunits led to increased Rubisco content, more rapid carbon assimilation, and greater plant height, most notably in plants at least 6 weeks of age. The magnitude of the effects was similar to what was observed previously for expression of Raf1 together with Rubisco subunits. Taken together, this suggests that increasing the amount of either assembly factor with Rubisco subunits can independently enhance Rubisco abundance and some aspects of plant performance. These results could also imply either synergy or a degree of functional redundancy for Raf1 and Raf2, the latter of whose precise role in Rubisco assembly is currently unknown.

CFM1, a member of the CRM-domain protein family, functions in chloroplast group II intron splicing in Setaria viridis.

The chloroplast RNA splicing and ribosome maturation (CRM) domain is a RNA-binding domain found in a plant-specific protein family whose characterized members play essential roles in splicing group I and group II introns in mitochondria and chloroplasts. Together, these proteins are required for splicing of the majority of the approximately 20 chloroplast introns in land plants. Here, we provide evidence from Setaria viridis and maize that an uncharacterized member of this family, CRM Family Member1 (CFM1), promotes the splicing of most of the introns that had not previously been shown to require a CRM domain protein. A Setaria mutant expressing mutated CFM1 was strongly disrupted in the splicing of three chloroplast tRNAs: trnI, trnV and trnA. Analyses by RNA gel blot and polysome association suggest that the tRNA deficiencies lead to compromised chloroplast protein synthesis and the observed whole-plant chlorotic phenotypes. Co-immunoprecipitation data demonstrate that the maize CFM1 ortholog is bound to introns whose splicing is disrupted in the cfm1 mutant. With these results, CRM domain proteins have been shown to promote the splicing of all but two of the introns found in angiosperm chloroplast genomes.

Setaria viridis chlorotic and seedling-lethal mutants define critical functions for chloroplast gene expression.

Deep insights into chloroplast biogenesis have been obtained by mutant analysis; however, in C4 plants a relevant mutant collection has only been developed and exploited for maize. Here, we report the initial characterization of an ethyl methyl sulfonate-induced mutant population for the C4 model Setaria viridis. Approximately 1000 M2 families were screened for the segregation of pale-green seedlings in the M3 generation, and a subset of these was identified to be deficient in post-transcriptional steps of chloroplast gene expression. Causative mutations were identified for three lines using deep sequencing-based bulked segregant analysis, and in one case confirmed by transgenic complementation. Using chloroplast RNA-sequencing and other molecular assays, we describe phenotypes of mutants deficient in PSRP7, a plastid-specific ribosomal protein, OTP86, an RNA editing factor, and cpPNP, the chloroplast isozyme of polynucleotide phosphorylase. The psrp mutant is globally defective in chloroplast translation, and has varying deficiencies in the accumulation of chloroplast-encoded proteins. The otp86 mutant, like its Arabidopsis counterpart, is specifically defective in editing of the rps14 mRNA; however, the conditional pale-green mutant phenotype contrasts with the normal growth of the Arabidopsis mutant. The pnp mutant exhibited multiple defects in 3' end maturation as well as other qualitative changes in the chloroplast RNA population. Overall, our collection opens the door to global analysis of photosynthesis and early seedling development in an emerging C4 model.

Rubisco production in maize mesophyll cells through ectopic expression of subunits and chaperones.

C4 plants, such as maize, strictly compartmentalize Rubisco to bundle sheath chloroplasts. The molecular basis for the restriction of Rubisco from the more abundant mesophyll chloroplasts is not fully understood. Mesophyll chloroplasts transcribe the Rubisco large subunit gene and, when normally quiescent transcription of the nuclear Rubisco small subunit gene family is overcome by ectopic expression, mesophyll chloroplasts still do not accumulate measurable Rubisco. Here we show that a combination of five ubiquitin promoter-driven nuclear transgenes expressed in maize leads to mesophyll accumulation of assembled Rubisco. These encode the Rubisco large and small subunits, Rubisco assembly factors 1 and 2, and the assembly factor Bundle sheath defective 2. In these plants, Rubisco large subunit accumulates in mesophyll cells, and appears to be assembled into a holoenzyme capable of binding the substrate analog CABP (carboxyarabinitol bisphosphate). Isotope discrimination assays suggest, however, that mesophyll Rubisco is not participating in carbon assimilation in these plants, most probably due to a lack of the substrate ribulose 1,5-bisphosphate and/or Rubisco activase. Overall, this work defines a minimal set of Rubisco assembly factors in planta and may help lead to methods of regulating the C4 pathway.

Systematic sequencing of chloroplast transcript termini from Arabidopsis thaliana reveals >200 transcription initiation sites and the extensive imprints of RNA-binding proteins and secondary structures.

Chloroplast transcription requires numerous quality control steps to generate the complex but selective mixture of accumulating RNAs. To gain insight into how this RNA diversity is achieved and regulated, we systematically mapped transcript ends by developing a protocol called Terminome-seq. Using Arabidopsis thaliana as a model, we catalogued >215 primary 5' ends corresponding to transcription start sites (TSS), as well as 1628 processed 5' ends and 1299 3' ends. While most termini were found in intergenic regions, numerous abundant termini were also found within coding regions and introns, including several major TSS at unexpected locations. A consistent feature was the clustering of both 5' and 3' ends, contrasting with the prevailing description of discrete 5' termini, suggesting an imprecision of the transcription and/or RNA processing machinery. Numerous termini correlated with the extremities of small RNA footprints or predicted stem-loop structures, in agreement with the model of passive RNA protection. Terminome-seq was also implemented for pnp1-1, a mutant lacking the processing enzyme polynucleotide phosphorylase. Nearly 2000 termini were altered in pnp1-1, revealing a dominant role in shaping the transcriptome. In summary, Terminome-seq permits precise delineation of the roles and regulation of the many factors involved in organellar transcriptome quality control.

Increased Rubisco content in maize mitigates chilling stress and speeds recovery.

Many C4 plants, including maize, perform poorly under chilling conditions. This phenomenon has been linked in part to decreased Rubisco abundance at lower temperatures. An exception to this is chilling-tolerant Miscanthus, which is able to maintain Rubisco protein content under such conditions. The goal of this study was to investigate whether increasing Rubisco content in maize could improve performance during or following chilling stress. Here, we demonstrate that transgenic lines overexpressing Rubisco large and small subunits and the Rubisco assembly factor RAF1 (RAF1-LSSS), which have increased Rubisco content and growth under control conditions, maintain increased Rubisco content and growth during chilling stress. RAF1-LSSS plants exhibited 12% higher CO2 assimilation relative to nontransgenic controls under control growth conditions, and a 17% differential after 2 weeks of chilling stress, although assimilation rates of all genotypes were ~50% lower in chilling conditions. Chlorophyll fluorescence measurements showed RAF1-LSSS and WT plants had similar rates of photochemical quenching during chilling, suggesting Rubisco may not be the primary limiting factor that leads to poor performance in maize under chilling conditions. In contrast, RAF1-LSSS had improved photochemical quenching before and after chilling stress, suggesting that increased Rubisco may help plants recover faster from chilling conditions. Relatively increased leaf area, dry weight and plant height observed before chilling in RAF1-LSSS were also maintained during chilling. Together, these results demonstrate that an increase in Rubisco content allows maize plants to better cope with chilling stress and also improves their subsequent recovery, yet additional modifications are required to engineer chilling tolerance in maize.

The Evolution of Gene Expression Underlying Vision Loss in Cave Animals.

Dissecting the evolutionary genetic processes underlying eye reduction and vision loss in obligate cave-dwelling organisms has been a long-standing challenge in evolutionary biology. Independent vision loss events in related subterranean organisms can provide critical insight into these processes as well as into the nature of convergent loss of complex traits. Advances in evolutionary developmental biology have illuminated the significant role of heritable gene expression variation in the evolution of new forms. Here, we analyze gene expression variation in adult eye tissue across the freshwater crayfish, representing four independent vision-loss events in caves. Species and individual expression patterns cluster by eye function rather than phylogeny, suggesting convergence in transcriptome evolution in independently blind animals. However, this clustering is not greater than what is observed in surface species with conserved eye function after accounting for phylogenetic expectations. Modeling expression evolution suggests that there is a common increase in evolutionary rates in the blind lineages, consistent with a relaxation of selective constraint maintaining optimal expression levels. This is evidence for a repeated loss of expression constraint in the transcriptomes of blind animals and that convergence occurs via a similar trajectory through genetic drift.

Arabidopsis chloroplast mini-ribonuclease III participates in rRNA maturation and intron recycling.

RNase III proteins recognize double-stranded RNA structures and catalyze endoribonucleolytic cleavages that often regulate gene expression. Here, we characterize the functions of RNC3 and RNC4, two Arabidopsis thaliana chloroplast Mini-RNase III-like enzymes sharing 75% amino acid sequence identity. Whereas rnc3 and rnc4 null mutants have no visible phenotype, rnc3/rnc4 (rnc3/4) double mutants are slightly smaller and chlorotic compared with the wild type. In Bacillus subtilis, the RNase Mini-III is integral to 23S rRNA maturation. In Arabidopsis, we observed imprecise maturation of 23S rRNA in the rnc3/4 double mutant, suggesting that exoribonucleases generated staggered ends in the absence of specific Mini-III-catalyzed cleavages. A similar phenotype was found at the 3' end of the 16S rRNA, and the primary 4.5S rRNA transcript contained 3' extensions, suggesting that Mini-III catalyzes several processing events of the polycistronic rRNA precursor. The rnc3/4 mutant showed overaccumulation of a noncoding RNA complementary to the 4.5S-5S rRNA intergenic region, and its presence correlated with that of the extended 4.5S rRNA precursor. Finally, we found rnc3/4-specific intron degradation intermediates that are probable substrates for Mini-III and show that B. subtilis Mini-III is also involved in intron regulation. Overall, this study extends our knowledge of the key role of Mini-III in intron and noncoding RNA regulation and provides important insight into plastid rRNA maturation.

Activation of an Endoribonuclease by Non-intein Protein Splicing.

The Chlamydomonas reinhardtii chloroplast-localized poly(A)-binding protein RB47 is predicted to contain a non-conserved linker (NCL) sequence flanked by highly conserved N- and C-terminal sequences, based on the corresponding cDNA. RB47 was purified from chloroplasts in association with an endoribonuclease activity; however, protein sequencing failed to detect the NCL. Furthermore, while recombinant RB47 including the NCL did not display endoribonuclease activity in vitro, versions lacking the NCL displayed strong activity. Both full-length and shorter forms of RB47 could be detected in chloroplasts, with conversion to the shorter form occurring in chloroplasts isolated from cells grown in the light. This conversion could be replicated in vitro in chloroplast extracts in a light-dependent manner, where epitope tags and protein sequencing showed that the NCL was excised from a full-length recombinant substrate, together with splicing of the flanking sequences. The requirement for endogenous factors and light differentiates this protein splicing from autocatalytic inteins, and may allow the chloroplast to regulate the activation of RB47 endoribonuclease activity. We speculate that this protein splicing activity arose to post-translationally repair proteins that had been inactivated by deleterious insertions or extensions.

A protein with an inactive pterin-4a-carbinolamine dehydratase domain is required for Rubisco biogenesis in plants.

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) plays a critical role in sustaining life by catalysis of carbon fixation in the Calvin-Benson pathway. Incomplete knowledge of the assembly pathway of chloroplast Rubisco has hampered efforts to fully delineate the enzyme's properties, or seek improved catalytic characteristics via directed evolution. Here we report that a Mu transposon insertion in the Zea mays (maize) gene encoding a chloroplast dimerization co-factor of hepatocyte nuclear factor 1 (DCoH)/pterin-4α-carbinolamine dehydratases (PCD)-like protein is the causative mutation in a seedling-lethal, Rubisco-deficient mutant named Rubisco accumulation factor 2 (raf2-1). In raf2 mutants newly synthesized Rubisco large subunit accumulates in a high-molecular weight complex, the formation of which requires a specific chaperonin 60-kDa isoform. Analogous observations had been made previously with maize mutants lacking the Rubisco biogenesis proteins RAF1 and BSD2. Chemical cross-linking of maize leaves followed by immunoprecipitation with antibodies to RAF2, RAF1 or BSD2 demonstrated co-immunoprecipitation of each with Rubisco small subunit, and to a lesser extent, co-immunoprecipitation with Rubisco large subunit. We propose that RAF2, RAF1 and BSD2 form transient complexes with the Rubisco small subunit, which in turn assembles with the large subunit as it is released from chaperonins.

Ribulose-1,5-bis-phosphate carboxylase/oxygenase accumulation factor1 is required for holoenzyme assembly in maize.

Most life is ultimately sustained by photosynthesis and its rate-limiting carbon fixing enzyme, ribulose-1,5-bis-phosphate carboxylase/oxygenase (Rubisco). Although the structurally comparable cyanobacterial Rubisco is amenable to in vitro assembly, the higher plant enzyme has been refractory to such manipulation due to poor understanding of its assembly pathway. Here, we report the identification of a chloroplast protein required for Rubisco accumulation in maize (Zea mays), RUBISCO ACCUMULATION FACTOR1 (RAF1), which lacks any characterized functional domains. Maize lines lacking RAF1 due to Mutator transposon insertions are Rubisco deficient and seedling lethal. Analysis of transcripts and proteins showed that Rubisco large subunit synthesis in raf1 plants is not compromised; however, newly synthesized Rubisco large subunit appears in a high molecular weight form whose accumulation requires a specific chaperonin 60 isoform. Gel filtration analysis and blue native gels showed that endogenous and recombinant RAF1 are trimeric; however, following in vivo cross-linking, RAF1 copurifies with Rubisco large subunit, suggesting that they interact weakly or transiently. RAF1 is predominantly expressed in bundle sheath chloroplasts, consistent with a Rubisco accumulation function. Our results support the hypothesis that RAF1 acts during Rubisco assembly by releasing and/or sequestering the large subunit from chaperonins early in the assembly process.

RNase J participates in a pentatricopeptide repeat protein-mediated 5' end maturation of chloroplast mRNAs.

Nucleus-encoded ribonucleases and RNA-binding proteins influence chloroplast gene expression through their roles in RNA maturation and stability. One mechanism for mRNA 5' end maturation posits that sequence-specific pentatricopeptide repeat (PPR) proteins define termini by blocking the 5'→3' exonucleolytic activity of ribonuclease J (RNase J). To test this hypothesis in vivo, virus-induced gene silencing was used to reduce the expression of three PPR proteins and RNase J, both individually and jointly, in Nicotiana benthamiana. In accordance with the stability-conferring function of the PPR proteins PPR10, HCF152 and MRL1, accumulation of the cognate RNA species atpH, petB and rbcL was reduced when the PPR-encoding genes were silenced. In contrast, RNase J reduction alone or combined with PPR deficiency resulted in reduced abundance of polycistronic precursor transcripts and mature counterparts, which were replaced by intermediately sized species with heterogeneous 5' ends. We conclude that RNase J deficiency can partially mask the absence of PPR proteins, and that RNase J is capable of processing chloroplast mRNAs up to PPR protein-binding sites. These findings support the hypothesis that RNase J is the major ribonuclease responsible for maturing chloroplast mRNA 5' termini, with RNA-binding proteins acting as barriers to its activity.

Ribonuclease II preserves chloroplast RNA homeostasis by increasing mRNA decay rates, and cooperates with polynucleotide phosphorylase in 3' end maturation.

Ribonuclease R (RNR1) and polynucleotide phosphorylase (cpPNPase) are the two known 3'→5' exoribonucleases in Arabidopsis chloroplasts, and are involved in several aspects of rRNA and mRNA metabolism. In this work, we show that mutants lacking both RNR1 and cpPNPase exhibit embryo lethality, akin to the non-viability of the analogous double mutant in Escherichia coli. We were successful, however, in combining an rnr1 null mutation with weak pnp mutant alleles, and show that the resulting chlorotic plants display a global reduction in RNA abundance. Such a counterintuitive outcome following the loss of RNA degradation activity suggests a major importance of RNA maturation as a determinant of RNA stability. Detailed analysis of the double mutant demonstrates that the enzymes catalyze a two-step maturation of mRNA 3' ends, with RNR1 polishing 3' termini created by cpPNPase. The bulky quaternary structure of cpPNPase compared with RNR1 could explain this activity split between the two enzymes. In contrast to the double mutants, the rnr1 single mutant overaccumulates most mRNA species when compared with the wild type. The excess mRNAs in rnr1 are often present in non-polysomal fractions, and half-life measurements demonstrate a substantial increase in the stability of most mRNA species tested. Together, our data reveal the cooperative activity of two 3'→5' exoribonucleases in chloroplast mRNA 3' end maturation, and support the hypothesis that RNR1 plays a significant role in the destabilization of mRNAs unprotected by ribosomes.

Chloroplast RNase J compensates for inefficient transcription termination by removal of antisense RNA.

Ribonuclease J is an essential enzyme, and the Bacillus subtilis ortholog possesses both endoribonuclease and 5' → 3' exoribonuclease activities. Chloroplasts also contain RNase J, which has been postulated to participate, as both an exo- and endonuclease, in the maturation of polycistronic mRNAs. Here we have examined recombinant Arabidopsis RNase J and found both 5' → 3' exoribonuclease and endonucleolytic activities. Virus-induced gene silencing was used to reduce RNase J expression in Arabidopsis and Nicotiana benthamiana, leading to chlorosis but surprisingly few disruptions in the cleavage of polycistronic rRNA and mRNA precursors. In contrast, antisense RNAs accumulated massively, suggesting that the failure of chloroplast RNA polymerase to terminate effectively leads to extensive symmetric transcription products that are normally eliminated by RNase J. Mung bean nuclease digestion and polysome analysis revealed that this antisense RNA forms duplexes with sense strand transcripts and prevents their translation. We conclude that a major role of chloroplast RNase J is RNA surveillance to prevent overaccumulation of antisense RNA, which would otherwise exert deleterious effects on chloroplast gene expression.

Overaccumulation of the chloroplast antisense RNA AS5 is correlated with decreased abundance of 5S rRNA in vivo and inefficient 5S rRNA maturation in vitro.

Post-transcriptional regulation in the chloroplast is exerted by nucleus-encoded ribonucleases and RNA-binding proteins. One of these ribonucleases is RNR1, a 3'-to-5' exoribonuclease of the RNase II family. We have previously shown that Arabidopsis rnr1-null mutants exhibit specific abnormalities in the expression of the rRNA operon, including the accumulation of precursor 23S, 16S, and 4.5S species and a concomitant decrease in the mature species. 5S rRNA transcripts, however, accumulate to a very low level in both precursor and mature forms, suggesting that they are unstable in the rnr1 background. Here we demonstrate that rnr1 plants overaccumulate an antisense RNA, AS5, that is complementary to the 5S rRNA, its intergenic spacer, and the downstream trnR gene, which encodes tRNA(Arg), raising the possibility that AS5 destabilizes 5S rRNA or its precursor and/or blocks rRNA maturation. To investigate this, we used an in vitro system that supports 5S rRNA and trnR processing. We show that AS5 inhibits 5S rRNA maturation from a 5S-trnR precursor, and shorter versions of AS5 demonstrate that inhibition requires intergenic sequences. To test whether the sense and antisense RNAs form double-stranded regions in vitro, treatment with the single-strand-specific mung bean nuclease was used. These results suggest that 5S-AS5 duplexes interfere with a sense-strand secondary structure near the endonucleolytic cleavage site downstream from the 5S rRNA coding region. We hypothesize that these duplexes are degraded by a dsRNA-specific ribonuclease in vivo, contributing to the 5S rRNA deficiency observed in rnr1.

MRL1, a conserved Pentatricopeptide repeat protein, is required for stabilization of rbcL mRNA in Chlamydomonas and Arabidopsis.

We identify and functionally characterize MRL1, a conserved nuclear-encoded regulator of the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase. The nonphotosynthetic mrl1 mutant of Chlamydomonas reinhardtii lacks ribulose-1,5-bisphosphate carboxylase/oxygenase, and the resulting block in electron transfer is partially compensated by redirecting electrons toward molecular oxygen via the Mehler reaction. This allows continued electron flow and constitutive nonphotochemical quenching, enhancing cell survival during illumination in spite of photosystem II and photosystem I photoinhibition. The mrl1 mutant transcribes rbcL normally, but the mRNA is unstable. The molecular target of MRL1 is the 5 ' untranslated region of rbcL. MRL1 is located in the chloroplast stroma, in a high molecular mass complex. Treatment with RNase or deletion of the rbcL gene induces a shift of the complex toward lower molecular mass fractions. MRL1 is well conserved throughout the green lineage, much more so than the 10 other pentatricopeptide repeat proteins found in Chlamydomonas. Depending upon the organism, MRL1 contains 11 to 14 pentatricopeptide repeats followed by a novel MRL1-C domain. In Arabidopsis thaliana, MRL1 also acts on rbcL and is necessary for the production/stabilization of the processed transcript, presumably because it acts as a barrier to 5 ' >3 ' degradation. The Arabidopsis mrl1 mutant retains normal levels of the primary transcript and full photosynthetic capacity.

Amphibian myelopoiesis.

Macrophage-lineage cells are indispensable to immunity and physiology of all vertebrates. Amongst these, amphibians represent a key stage in vertebrate evolution and are facing decimating population declines and extinctions, in large part due to emerging infectious agents. While recent studies indicate that macrophages and related innate immune cells are critically involved during these infections, much remains unknown regarding the ontogeny and functional differentiation of these cell types in amphibians. Accordingly, in this review we coalesce what has been established to date about amphibian blood cell development (hematopoiesis), the development of key amphibian innate immune cells (myelopoiesis) and the differentiation of amphibian macrophage subsets (monopoiesis). We explore the current understanding of designated sites of larval and adult hematopoiesis across distinct amphibian species and consider what mechanisms may lend to these species-specific adaptations. We discern the identified molecular mechanisms governing the functional differentiation of disparate amphibian (chiefly Xenopus laevis) macrophage subsets and describe what is known about the roles of these subsets during amphibian infections with intracellular pathogens. Macrophage lineage cells are at the heart of so many vertebrate physiological processes. Thus, garnering greater understanding of the mechanisms responsible for the ontogeny and functionality of these cells in amphibians will lend to a more comprehensive view of vertebrate evolution.

Delivery mode impacts gut bacteriophage colonization during infancy.

Background: Cesarean section delivery is associated with altered early-life bacterial colonization and later adverse inflammatory and immune health outcomes. Although gut bacteriophages can alter gut microbiome composition and impact host immune responses, little is known about how delivery mode impacts bacteriophage colonization over time. To begin to address this we examined how delivery mode affected bacteriophage colonization over the first two years of life. Results: Shotgun metagenomic sequencing was conducted on 272 serial stool samples from 55 infants, collected at 1-2 days of life and 2, 6, 12 and 24 months. 33/55 (60%) infants were born by vaginal delivery. DNA viruses were identified, and by host inference, 94% of the viral sequences were found to be bacteriophages. Alpha diversity of the virome was increased in vaginally delivered infants compared to cesarean section delivered infants at 2 months (Shannon index, p=0.022). Beta diversity significantly differed by delivery mode at 2, 6, and 12 months when stratified by peripartum antibiotic use (Bray-Curtis dissimilarity, all p<0.05). Significant differentially abundant predicted bacteriophage hosts by delivery mode were seen at all time points. Moreover, there were differences in predicted bacteriophage functional gene abundances up to 24 months by delivery mode. Many of the functions considered to play a role in host response were increased in vaginal delivery. Conclusions: Clear differences in bacteriophage composition and function were seen by delivery mode over the first two years of life. Given that phages are known to affect host immune response, our results suggest that future investigation into how delivery mode may lead to adverse inflammatory outcomes should not only include bacterial microbial colonization but also the potential role of bacteriophages and transkingdom interactions.