Phages infecting Faecalibacterium prausnitzii belong to novel viral genera that help to decipher intestinal viromes
© The Author(s). 2018
Received: 24 October 2017
Accepted: 21 March 2018
Published: 3 April 2018
Viral metagenomic studies have suggested a role for bacteriophages in intestinal dysbiosis associated with several human diseases. However, interpretation of viral metagenomic studies is limited by the lack of knowledge of phages infecting major human gut commensal bacteria, such as Faecalibacterium prausnitzii, a bacterial symbiont repeatedly found depleted in inflammatory bowel disease (IBD) patients. In particular, no complete genomes of phages infecting F. prausnitzii are present in viral databases.
We identified 18 prophages in 15 genomes of F. prausnitzii, used comparative genomics to define eight phage clades, and annotated the genome of the type phage of each clade. For two of the phages, we studied prophage induction in vitro and in vivo in mice. Finally, we aligned reads from already published viral metagenomic data onto the newly identified phages.
We show that each phage clade represents a novel viral genus and that a surprisingly large fraction of them (10 of the 18 phages) codes for a diversity-generating retroelement, which could contribute to their adaptation to the digestive tract environment. We obtained either experimental or in silico evidence of activity for at least one member of each genus. In addition, four of these phages are either significantly more prevalent or more abundant in stools of IBD patients than in those of healthy controls.
Since IBD patients generally have less F. prausnitzii in their microbiota than healthy controls, the higher prevalence or abundance of some of its phages may indicate that they are activated during disease. This in turn suggests that phages could trigger or aggravate F. prausnitzii depletion in patients. Our results show that prophage detection in sequenced strains of the microbiota can usefully complement viral metagenomic studies.
Faecalibacterium prausnitzii is one of the most abundant bacterial species of the human gastrointestinal tract (GIT) microbiota. This species is an important producer of butyrate, a short-chain fatty acid with a broad range of beneficial properties on human health [1, 2]. F. prausnitzii also exerts anti-inflammatory activity through the production of the MAM protein  and provides protection against inflammation in different animal models [3–6]. In humans, a low abundance of this species has been regularly reported in the gut microbiota of inflammatory bowel disease (IBD) patients [6, 7]. Depletion of F. prausnitzii in patients has been associated with a higher risk of recurrence of flares [4, 8], possibly in relation with the capacity of this species to modulate T cell responses . Bacteriophages, or phages, the viruses infecting bacteria, might play a role in this depletion . But, although viral metagenomics studies have provided indications that F. prausnitzii phages are present in human fecal samples , no complete F. prausnitzii phage has been characterized to date. Phages can be either virulent or temperate: virulent phages perform only lytic cycles, during which new phage particles are produced and the infected bacteria are killed. Temperate phages can similarly replicate through lytic cycles, but they are also able to enter a latent lysogenic state. During lysogeny, the repressed phage, known as a prophage, is generally integrated into the bacterial chromosome, and the expression of most of its genes is repressed. Nevertheless, a few genes continue to be expressed, including the repressor of the lytic cycle and a class of genes known as morons whose expression often modifies the physiology of the host bacterium. In this way, prophages can modify numerous phenotypes of their bacterial host, such as resistance against phage attack, stress resistance, pathogenicity, or ability to form biofilms [11, 12]. In response to a specific signal, active prophages can switch back to lytic multiplication, a process known as induction. Because of this, active prophages can be highly detrimental to their host in some inducing conditions .
Therefore, not only virulent phages but also temperate ones have been suspected to play a role in dysbiosis of intestinal microbiota [14, 15]. To examine this possibility, several viral metagenomic studies of the intestinal microbiota have been performed in recent years [10, 16–18]. However, a major limitation of these studies is the difficulty of viral sequence analysis. Indeed, annotation of viral sequences, and in particular the prediction of the bacterial host, is mainly based on alignment to reference viral sequences in databases , and limitations in the viral databases render annotation of most viral sequences difficult if not impossible. Limitations are particularly flagrant concerning phages infecting important groups such as Ruminococcaceae, previously known as Clostridium cluster IV. There is, to our knowledge, no phage infecting members of this family whose complete genome has been deposited in viral databases. Isolating and describing phages infecting Clostridiaceae species such as F. prausnitzii is thus a prerequisite for investigation of their potential role in dysbiosis. In this study, we characterize for the first time complete F. prausnitzii temperate phages and show that they belong to eight new proposed viral genera. In addition, we show that some of them are either more prevalent or more abundant in the fecal samples of IBD patients compared to healthy controls, suggesting these phages might play a role in the disease.
Identification of 23 prophages corresponding to novel proposed genera
Main characteristics of the eight proposed F. prausnitzii phage genera
Strains with complete prophagea
Evidence of activity
Genome size, kb
Myoviridae type 1, cluster 6
A2-165, B. hansenii 20583, CM62/1
Myoviridae type1, cluster 8 “Saltoviridae”
A2-165, C_4543, C_4574
Ruminococcaceae bacterium D16
R. bacterium D16
Siphoviridae type 1, cluster 1
Myoviridae type 1, cluster 6
Myoviridae type 1, cluster 6
Myoviridae type 1, cluster 6
Myoviridae type 1, cluster 7
Siphoviridae type 1, cluster 2
Given the ubiquity of F. prausnitzii in human gut microbiota, we investigated whether similar phage sequences had been found in other studies. To this end, we compared the prophage genomes against the NCBI nr database. Prophages Lagaffe, Mushu, and ToutatisSL3/3 have homology to four large contigs assembled from reads obtained by sequencing DNA from virus-like particles (VLPs) isolated from human fecal samples  (Fig. 1). In addition, ToutatisL2/6 is homologous to a 9-kb contig assembled from reads of the VLPs of a human fecal sample. This contig had been predicted to correspond to a F. prausnitzii prophage . This provides further support for the proposition that these prophage regions are of viral nature. We then looked for CRISPR spacers matching the prophages in the 15 F. prausnitzii genomes, and found spacers for eight prophages (Table 1), which is indicative of recent infection of these strains by the respective phages . Overall, indications of activity were obtained for at least one prophage of each proposed genus, as summarized in Table 1.
For six of the eight prophage groups, the host range seemed to be restricted to F. prausnitzii, since no similar prophage could be detected by BLASTn in other genomes of the NCBI RefSeq database. On the other hand, prophages with 65% overall sequence identity to Mushu prophage were found in the unclassified Ruminococcaceae bacterium D16 (Fig. 1 and Additional file 3) and in the unclassified bacterium SS3/4 (Genbank sequence FP929062.1). Even more striking was the detection of a prophage almost identical to Lagaffe (98% nucleotide identity) in the genome of Blautia hansenii DSM 20583 (Fig. 1 and Additional file 3), a bacterium only related to F. prausnitzii A2-165 at the order level (they share only 29% of their proteins).
To classify these proposed new phage genera, since they greatly differ from known viruses, we used the Virfam webserver (biodev.cea.fr/virfam). Virfam uses virion protein remote homology and gene synteny to classify phages . All the phages were classified among the Neck type 1, the largest clade of the Caudoviridae order (Additional file 4). Four genera grouped within the cluster 6 of this neck type, which branches deeply in the overall tree and comprises phages that could not be attributed to better-resolved clusters. Phage Mushu belongs to cluster 8, which corresponds to the “Saltoviridae” proposed family, composed of Mu-like transposable phages . Accordingly, Mushu and its relative from the Ruminococcaceae bacterium D16 genome have been assigned to this proposed family [29, 30]. Phages from the proposed “FPlughvirus” genus and phage Oengus group in clusters 1 and 2 respectively, with other unclassified Siphoviridae phages infecting Firmicutes. Surprisingly, Brigit phage groups in cluster 7, which typically comprises phages infecting Enterobacteriaceae.
Functional annotation of prophage genes reveals a high prevalence of diversity-generating retroelements
A remarkable feature of these genomes is the high prevalence of diversity-generating retroelements (DGRs). DGRs generate variability in target genes through a reverse transcriptase (RT)-mediated mechanism that introduces nucleotide substitutions at defined locations in specific genes (reviewed in ). In this process, adenine residues are converted to random nucleotides in a unidirectional process from a donor template repeat (TR) to a recipient variable repeat (VR) present in the target gene. TR and VR sequences were identified in each prophage type coding for an RT (Additional file 5). As with other DGRs, TR/VR sequences were between 100 and 135 bp long, they differed almost exclusively at sites corresponding to TR adenines, and the VRs were located at the 3′ end of protein-coding genes. We then searched for analogues of Avd, a low-molecular-weight protein required for mutagenic homing in the prototype DGR of the Bordetella pertussis phage BPP1 . Homologs of Avd were detected in Mushu and Lagaffe genomes only (Fig. 2). In other phages, DGRs might function with another small protein, as has previously been suggested for other DGRs . In Lagaffe and Epona, the putative target gene encodes a protein belonging to the formylglycine-generating enzyme (FGE) superfamily, which includes the BPP1 tail protein MTD, implicated in bacterial recognition . Mushu and Brigit DGR target genes code for a Hoc-like capsid decoration protein, which in T4 improves the binding to intestinal mucus . In Toutatis and Taranis, no putative function could be attributed to the DGR target.
Lagaffe and Mushu produce virions
For both phages, reads spanning the VR region of the DGR are more variable than the rest of the genome, an evidence of DGR activity (Additional file 5C). The level of variability observed in Mushu is however remarkable, since the phage genomes sequenced result in theory from a single cycle of multiplication (they are produced following induction in lysogenic cells), and suggests a very high efficacy of TR-to-VR transfer in this DGR.
Lagaffe and Mushu virion concentration in culture supernatants and in mouse feces
Virions/ml or virions/g with F. prausnitzii A2165
Virions/ml in B. hansenii cultures
Healthy mouse gut
Inflamed mouse gut
3 ± 0.3 × 105
< 1 × 106
< 1 × 106
6 ± 3 × 105
8 ± 2 × 106
9 ± 2 × 106
7 ± 2 × 106
Lagaffe and Mushu prophages are not inducible by DNA damage
To increase virion titers and search for inducing signals, A2-165 F. prausnitzii cultures were treated with different concentrations of three potent inducers of phage lytic cycle: mitomycin C, hydrogen peroxide, and bile salts (cholate/deoxycholate [1:1]) [13, 42]. Mitomycin C in particular causes DNA damage, which is a general signal for prophage derepression. However, none of the treatments significantly increased virion concentrations (Additional file 7). For each compound tested, the highest dose tested almost abolished bacterial growth, indicating that they are toxic for F. prausnitzii and that the absence of prophage induction is not due to too low concentrations of potent inducers. Lagaffe and Mushu are thus probably not inducible by DNA damage, similar to most Mu-like phages.
Mushu virions are produced in the mouse gastro-intestinal tract but are not increased by colitis
We next investigated the production of A2-165 prophages in the mouse gastro-intestinal tract (GIT). Indeed, recent studies have shown increased prophage induction in the murine GIT, as compared to in vitro conditions [43, 44]. In order to obtain a colonization level of F. prausnitzii sufficient to detect virions in feces, we used dixenic mice colonized with F. prausnitzii A2-165 and Escherichia coli. In these mouse feces, the Lagaffe virion concentration was too low to be detected, but Mushu virions were present in similar concentrations to those observed in vitro (Table 2 and Additional file 8). Of note, however, the potential binding of virions to bacterial cells and intestinal mucus prevents the direct deduction of prophage induction rates from the concentrations of free virions. We then monitored the effect of moderate inflammation on the prophage induction rate. To mimic moderate inflammatory conditions, dextran sodium sulfate (DSS) was added at 2% to drinking water for 7 days. This treatment resulted in mild inflammation but did not increase virion concentrations in mouse feces (Additional file 8).
Some F. prausnitzii phages are more abundant in the gut viromes of IBD subjects
F. prausnitzii genomes are particularly rich in prophages: in this study, we found that 85% of them harbor at least one complete prophage, as opposed to only 46% of all completely sequenced bacterial genomes . The 23 identified complete prophages in this study correspond to 18 different phages, which can be grouped into eight new viral genera. Thanks to experimental and in silico investigations, we found evidence of activity for at least one member of each proposed genus. Their novelty and diversity are striking and greatly enrich the viral databases, which are at present highly skewed towards phages infecting proteobacteria or other pathogens, which together account for more than 75% of phages in the NCBI viral genome database (as of December 2017). The novelty of the viral proteins coded by F. prausnitzii phages should not only improve the analysis of viral metagenomic sequences but might also lead to biotechnology applications, such as genetic engineering of F. prausnitzii.
A very high proportion F. prausnitzii phages possesses DGRs: these diversification systems are present in 10 out of the 18 phages. By comparison, only 0.05% of the 876 phages present in the Phagonaute database (comprising all fully sequenced phage genomes available at the NCBI as of May 2014) possess DGRs. This high prevalence in F. prausnitzii phages might explain the previously observed high frequency of DGRs in intestinal viromes . By generating variants of target genes, DGRs presumably confer selective advantages to the phages, such as the ability to infect bacteria with different variants of the surface component used as receptor by the phage, as has been shown in the case of B. pertussis phage BPP1. In this phage, the DGR-mediated diversification of the MTD protein generates at low-frequency phage variants able to infect bacteria with receptors different from those of the original host . In phages Lagaffe and Epona, the proteins diversified by the DGR are homologous to BPP1 MTD, so their diversification probably also facilitates phage adaptation to different bacterial receptors. This property would provide a selective advantage to the phage either in the case of a preexisting high bacterial surface receptor diversity or in conditions of continually evolving bacteria in the context of an ongoing arms race between phage and their hosts. In phages Brigit and Mushu, by targeting Hoc-like capsid decoration proteins, the DGR system could improve the binding of virions to mucus, a strategy that has been suggested to increase the probability for phages to encounter their bacterial host . Whatever the role of these DGRs, they contribute to the generation of microdiversity regions in the genomes of intestinal phages, which notoriously complicates the assembly of reads, and hence the analysis of GIT viral metagenomes .
We found that the F. prausnitzii temperate phages described here are present in 5 to 40% of the human fecal samples from the study of Norman and colleagues . Despite their relatively low abundance in most samples of this dataset, they can represent up to 22% of total viral reads (Additional file 10). Most importantly, some of them are either more prevalent or more abundant in samples from IBD patients compared to that of healthy controls. This result is striking since IBD patients have generally a lower abundance of F. prausnitzii in their microbiota [6, 7]. This suggests that this higher prevalence or abundance is not just the consequence of a higher abundance of their bacterial host but rather that these F. prausnitzii prophages are more active in the microbiota of IBD patients as compared to healthy controls. This phage activity could aggravate or even trigger F. prausnitzii depletion. Since this depletion is highly suspected to participate in intestinal inflammation, our results suggest for the first time a role of F. prausnitzii phages in IBD.
Here, we describe for the first time complete F. prausnitzii temperate phage genomes. These phages belong to the new viral genera, and most of their proteins have only remote homology with previously described viral proteins. The novelty of F. prausnitzii phages warrants great enrichment of databases, which should in turn significantly improve annotation of viral gut metagenomic sequences. In addition, we show that some of these F. prausnitzii phages are either more prevalent or more abundant in the fecal samples of IBD patients compared to that of healthy controls, suggesting that these phages might play a role in the disease. Future mechanistic exploration of the nature of phage-bacteria interactions is now necessary to elucidate the role of temperate phages in dysbiosis.
Prophage detection and comparison
F. prausnitzii genomes listed in Additional file 1 were analyzed with PHASTER (http://phaster.ca). “Intact” and “questionable” regions were then validated by manual inspection. We retained regions with genes in each of the following modules: replication, capsid, lysis, and lysogeny, and no genes specific of integrative plasmids or insertion sequences. Of note, several prophage regions were on borders of contigs, suggesting that some complete prophages could have been missed because of splitting between two contigs and that some of the retained prophages might lack some genes. Prophage borders were refined by searching for phage integrase and phage-specific genes by BLASTp against the nr NCBI database. In several cases, the ORF corresponding to the integrase had not been detected automatically and was added manually. For type prophages, few other ORFs were also added. In addition, repeats on both sides (attL and attR sites) were searched for. They could be identified only in prophages encoding a tyrosine integrase, which is not surprising since tyrosine integrases generally use longer identical repeat than serine integrases. A multi-fasta file containing all prophage nucleotide sequences is available in the Additional file 11. The type phage genomes were deposited in the NCBI viral database (accession numbers in Additional file 1) and are consultable on the Phagonaute webserver (genome.jouy.inra.fr/phagonaute).
Genomic regions homologous to the 18 prophages were searched for in the nr/nt nucleotide collection of the NCBI by BLASTn, with the discontinuous megaBLAST default parameters. Sequences producing alignments covering more than 5% of the query were retained. Whole genome dotplots were realized with Gepard , using a length of 11 nt for minimal exact matches. CRISPR spacers were searched for with CRISPRDetect , using default parameters.
Shared proteins between prophage pairs
For each prophage pair to be compared, all proteins of the first phage were compared to the proteins of the second one (BLASTp, E value < 10−3), and best hits were considered shared proteins. The proportion of shared proteins was then calculated by dividing the number of shared proteins by the total protein number (minimal value of the pair), and the average identity percent of the shared proteins was computed. To select an outgroup phage for each proposed genus, all proteins of the type phage were compared to the viral protein database subset of the NCBI. The most frequent subject phage among the top 10 results for each protein was selected as the outgroup phage.
Gene product prediction
Gene function prediction for type prophages was achieved by combining Virfam  and Phagonaute  web tools. Concerning Virfam, default homology thresholds were used for recombinase and connector proteins. In Phagonaute, the proteins having a HHsearch hit against Pfam (P > 95%) were annotated with the Pfam information. Otherwise, the protein was searched against all other phage proteins (probability threshold of 95% and two iterations of HHsearch) and annotated using the most prevalent annotation. If no homolog with experimentally validated function was retrieved by the second iteration, the word “putative” was indicated in the Genbank file.
Virome read treatments
Reads from 171 viromes extracted from feces of healthy subjects and patients with inflammatory bowel diseases  were downloaded from the EBI (fastq files, accession numbers ESR698653-ESR699023). Paired-end reads were treated with trimmomatic  to remove low-quality reads and unclip primer sequences present at the 3′ ends of reads, using the following command ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 SLIDINGWINDOW:4:15 MINLEN:36. Reads were then dereplicated with fastx uniques from the usearch9 suite. Since ~ 30% of reads were no longer paired after the trimmomatic treatment, only one read from each pair was used for mapping analyses, to remain quantitative and avoid over-representation of some phage sequences compared to others. Average read length was 250 bp. To complement this analysis, we also analyzed the viral contigs generated in the Reyes et al. study . In this study, fecal viruses from 20 Malawian families were extracted at different time points and sequenced by 454 shotgun pyrosequencing. A cross-assembly strategy yielded 17,696 contigs, among which 780 had a size above 10 kb. These contigs were downloaded from ENA (PRJEB 9818) and compared with the prophages by megaBLAST alignments. Four of these contigs (> 10 kb) shared 93–97% identity to Epona, and together, the similarity regions covered 73% of Epona, but each of these contigs also contained large segments dissimilar to Epona. All four contigs originated from reads of family 57. These reads were therefore downloaded from ENA (PRJEB 9818, samples ERS795750 to 66), quality treated with trimmomatic, and dereplicated with fastx_uniques. This generated 694,650 reads, among which 680 mapped on Epona, whereas only 71 mapped to bacterial F. prausnitzii DNA (see the mapping method below and Fig. 4).
Quantification of F. prausnitzii phages in viromes
The proportion of reads homologous to each prophage was determined with BLASTn (E value < 10−9, > 75% identity). To account for different genome lengths, read counts were normalized by dividing them by genome length and then multiplying the result by the average genome length (50 kb). This allows comparing phage read proportions as if all phages had a genome of length 50 kb. Since temperate phages can also be present in the bacterial DNA as prophages, we verified the level of bacterial DNA contamination of samples. An analysis of reads mapping to 16S rRNA genes with Bowtie2  indicated that the vast majority of the 171 viromes had a low proportion of bacterial DNA: 166 out of 171 had less than 10−4 reads matching to 16S rRNA sequences. Two additional specific controls of F. prausnitzii DNA presence were included for the mapping analysis: a 50-kb segment corresponding to the beginning of the F. prausnitzii A2-165 Genbank file and a 50-kb segment encoding mainly metabolic genes (coordinates 1,451,689–1,501,688). Samples were classified as phage-positive for a given phage when the proportion of reads homologous to that phage was significantly higher than the proportion of reads homologous to these two bacterial regions (one tailed two-proportion z test, p < 0.05). As a control of the method, we included the highly prevalent crAssphage in the analysis and found that it is more prevalent in healthy controls than in IBD patients (Fig. 5a, two-tailed two-proportion z test, p < 0.01), in agreement with previously published results . Inspection of prophage coverage within each virome (number of reads per 500 bp interval) revealed three situations: (i) the prophage was largely covered (> 50% length, see Fig. 4), (ii) coverage was partial (Taranis and Lugh in Fig. 4); and (iii) matches were concentrated into 1–2 contiguous prophage intervals. We excluded these last cases (less than three non-overlapping windows of 500 bp, corresponding to at most 5% of the genome) that correspond to spurious matches from various phages, mapping to a conserved region, such as an anti-repressor or an integrase gene. Additionally, to minimize the impact of these spurious matches, the median and average values of matches among all 500 bp intervals with a non-null amount of matches were computed for each phage/sample pair, and in cases were average was 10-fold higher to median, intervals displaying a number of matches above 10 times the median were considered spurious and replaced by an arbitrary value of 10 times the median. The final abundance matrix and coverage value for each phage/virome pair is reported in Additional file 10.
Viral reads assembly for phage Oengus
Given the elevated abundance of Oengus in virome ERR843984, all reads from this virome mapping on Oengus were retrieved (n = 193,414) and assembled with SPAdes (version 3.9) , using the following options: -k 21,33,55,77,99,127 --careful --only-assembler. This generated a 53.7-kb contig 99% identical to Oengus, covering 93% of the prophage, with a mean coverage of 322.
Sequencing of the DNA present in virions produced by the A2-165 strain
Five hundred milliliters of A2-165 culture grown overnight in LY-BHI (brain-heart infusion medium supplemented with 0.5% yeast extract [Difco], cellobiose (1 mg/ml [Sigma]), maltose (1 mg/ml [Sigma]), and cysteine (0.5 mg/ml [Sigma])), in an anaerobic chamber filled with 90% N2, 5% CO2, and 5% H2, was centrifuged at 5,200g for 15 min at 4 °C. Supernatant was recovered and centrifuged at 5,200g for 30 min at 4 °C. This step was repeated with 1 h of centrifugation. Polyethylene glycol (PEG 8000) and NaCl (Sigma) were then added to the resulting supernatant to a final concentration of 10% and 1 M respectively. The resulting mixture was kept at 4 °C overnight and then centrifuged at 5200g for 1 h at 4 °C. Supernatant was removed, and the phage pellet was resuspended in 2 ml of SM buffer (100 mM NaCl, 8 mM MgSO4, 50 mM Tris pH 7.5). Sample was then treated with 0.25 μg of RNAse A (Sigma) and Dnase I (Sigma) at 37 °C for 1 h. Phage DNA was extracted using Promega kit Wizard® DNA Clean-up system. Sequencing was performed with the Ion proton sequencing technology. Alignment of obtained reads was performed using bowtie2 (-N 1 –L 32) and then visualized with Tablet, using default parameters . We obtained 204,837 reads mapping to the A2-165 genome, of which 19,702 and 30,274 reads aligned to the Lagaffe and Mushu prophages respectively.
Prophage induction assay and quantitative PCR
Quantification was made as previously described . Briefly, phage particles from mouse feces (see below) or from 10 ml of overnight culture in LY-BHI were recovered by PEG precipitation, as described above. For induction assays, overnight cultures of F. prausnitzii A2-165 were diluted in LY-BHI 25-fold. When culture absorbance reached 0.2, potent inducers were added at final concentrations of 1 or 5 μg/ml for mitomycin C (Sigma, Ref: M4287), 0.1 and 0.05% for bile salts (Sigma 48305), and 100 μM and 10 μm for hydrogen peroxide (Sigma H1009). Five hours after induction, phages were recovered from culture supernatants by PEG precipitation as described above. Samples were then treated with 10 U of Turbo DNase (Ambion) for 1 h at 37 °C and incubated at 95 °C for 30 min in 0.2 ml PCR tubes. qPCR on 100-fold diluted samples in water was performed using the Takyon ROX SYBR Mastermix blue dTTP kit (Eurogentec) and the StepOnePlus real-time PCR system (Applied Biosystem). Phage and bacterial concentrations were determined using specific primer pairs (Additional file 12). Copy numbers were calculated with a standard curve built with bacterial DNA containing the prophages, whose concentration was estimated with Qubit (Thermo Fisher). The reaction mix was the following: 7.5 μl Takyon mix, 0.9 μl dH2O, 0.3 μl of each primer (200 nM final concentration), and 6 μl of DNA diluted in H2O. The PCR conditions were the following: 95 °C 15 s, 58 °C 45 s, and 72 °C 30s for 45 cycles, 72 °C 5 min, followed by melting curves. Results were analyzed using the StepOne Software 2.3.
Mouse colonization experiments
Germfree 7- to 8-week-old C3H/HeN mice (female) were obtained from the germfree rodent breeding facility of Anaxem-Micalis (INRA, Jouy-en-Josas, France). Animals were kept in flexible-film isolators (Getinge-La Calhène) in standard cages (two to five mice/cage) with sterile wood shavings as animal bedding. Mice were given free access to autoclaved tap water and to standard diet R03-40 (Scientific Animal food and Engineering, Augy, France), sterilized by gamma irradiation at 45 kGy (IBA Mediris). Isolators were maintained under controlled conditions of light (12 h), temperature (20 to 23 °C), and humidity (45 to 55%). To obtain dixenic E. coli/F. prausnitzii-diassociated mice, germfree mice were orally inoculated with a fresh culture of E. coli LF82 (108 to 109 CFU/ml). Forty-eight and 72 h later, 109 CFU of F. prausnitzii A2-165 in 0.2 ml were administered by intragastric gavage to mice pretreated with sodium bicarbonate (0.2 M, 0.1 ml by intragastric gavage, given 10 min before inoculation of bacteria). To induce colitis, 21 days after F. prausnitzii colonization, mice were given drinking water supplemented with 2% (w/v) dextran sulfate sodium (DSS; MP Biomedicals) for 7 days. Body weight, blood in stool, and stool consistency were monitored daily. Phage particles were recovered from one or two frozen fecal pellets unfrozen on ice for 10 min and diluted 40-fold in cold PBS. Re-suspended feces were kept on ice for 5 min with regular agitation prior to centrifugation for 10 min at 5251g at 4 °C. Supernatants were recovered and filtered through a 0.2-μm filter (PALL Corporation Acrodisc PF syringe filter). Samples were supplemented with 10% PEG 8000 (Sigma) and 0.5 M NaCl (Sigma). After incubation overnight at 4 °C, phage particles were harvested by centrifugation (5250g for 1 h at 4 °C with a swinging rotor). Supernatant was carefully removed, and pellets were re-suspended in 100 μl of SM buffer (100 mM NaCl, 8 mM MgSO4, 50 mM Tris pH 7.5) for 1 h at 4 °C.
Transmission electronic microscopy
Phages were concentrated by ultracentrifugation (20,000g, 1 h, 4 °C), starting from 500 ml of F. prausnitzii A2-165 culture. The phage pellet was resuspended in 100 μl of MgSO4 10 mM. Droplets of this preparation were directly placed on Formvar carbon-coated grids for 5 min. The grids were stained with 1% uranyl acetate and then viewed for TEM using a HITACHI HT 7700 (Elexience, France) at 80 kV. Microphotographies were acquired with a charge-coupled device camera AMT.
We are very grateful to Chantal Bridonneau for her expertise in anaerobic culture; to Anne Foussier, Fatima Joly, and Aurélie Balvay from the Anaxem facility for their help with the animal handlings; to Christine Longin from the MIMA2 facilities (UMR 1313 GABI, INRA) for the TEM observations; to the Migale platform (INRA) for the bio-informatics environment; and to Colin Tinsley for his careful reading of the manuscript.
This work was supported by a grant from Association François Aupetit.
Availability of data and materials
The genomes of bacteria used are freely available on the NCBI website; the accession numbers are given in Additional file 1.
JC and MDP performed the prophage detection. LB provided the F. prausnitzii genomes. JC and EM performed the experiments. MAP performed the viral metagenomic analyses. MM performed the statistical analyses. JC, MAP, PL, HS, and MDP analyzed the results. MDP wrote the manuscript with the help of MAP and MM, and all authors revised and approved the manuscript.
Ethics approval and consent to participate
All animal procedures were carried out according to the European Community Rules of Animal Care and with authorization 1234-2015101315238694 from French services.
Consent for publication
The manuscript does not contain any individual personal data in any form.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Rios-Covian D, Ruas-Madiedo P, Margolles A, Gueimonde M, de Los Reyes-Gavilan CG, Salazar N. Intestinal short chain fatty acids and their link with diet and human health. Front Microbiol. 2016;7:185.View ArticlePubMedPubMed CentralGoogle Scholar
- Canani RB, Costanzo MD, Leone L, Pedata M, Meli R, Calignano A. Potential beneficial effects of butyrate in intestinal and extraintestinal diseases. World J Gastroenterol. 2011;17(12):1519–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Quevrain E, Maubert MA, Michon C, Chain F, Marquant R, Tailhades J, Miquel S, Carlier L, Bermudez-Humaran LG, Pigneur B, et al. Identification of an anti-inflammatory protein from Faecalibacterium prausnitzii, a commensal bacterium deficient in Crohn's disease. Gut. 2016;65(3):415–25.View ArticlePubMedGoogle Scholar
- Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermudez-Humaran LG, Gratadoux JJ, Blugeon S, Bridonneau C, Furet JP, Corthier G, et al. Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci U S A. 2008;105(43):16731–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Miquel S, Martin R, Lashermes A, Gillet M, Meleine M, Gelot A, Eschalier A, Ardid D, Bermudez-Humaran LG, Sokol H, et al. Anti-nociceptive effect of Faecalibacterium prausnitzii in non-inflammatory IBS-like models. Sci Rep. 2016;6:19399.View ArticlePubMedPubMed CentralGoogle Scholar
- Rossi O, van Berkel LA, Chain F, Tanweer Khan M, Taverne N, Sokol H, Duncan SH, Flint HJ, Harmsen HJ, Langella P, et al. Faecalibacterium prausnitzii A2-165 has a high capacity to induce IL-10 in human and murine dendritic cells and modulates T cell responses. Sci Rep. 2016;6:18507.View ArticlePubMedPubMed CentralGoogle Scholar
- Cao Y, Shen J, Ran ZH. Association between Faecalibacterium prausnitzii reduction and inflammatory bowel disease: a meta-analysis and systematic review of the literature. Gastroenterol Res Pract. 2014;2014:872725.PubMedPubMed CentralGoogle Scholar
- Sokol H, Seksik P, Furet JP, Firmesse O, Nion-Larmurier I, Beaugerie L, Cosnes J, Corthier G, Marteau P, Dore J. Low counts of Faecalibacterium prausnitzii in colitis microbiota. Inflamm Bowel Dis. 2009;15(8):1183–9.View ArticlePubMedGoogle Scholar
- Mills S, Shanahan F, Stanton C, Hill C, Coffey A, Ross RP. Movers and shakers: influence of bacteriophages in shaping the mammalian gut microbiota. Gut Microbes. 2013;4(1):4–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Minot S, Sinha R, Chen J, Li H, Keilbaugh SA, Wu GD, Lewis JD, Bushman FD. The human gut virome: inter-individual variation and dynamic response to diet. Genome Res. 2011;21(10):1616–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Bondy-Denomy J, Davidson AR. When a virus is not a parasite: the beneficial effects of prophages on bacterial fitness. J Microbiol. 2014;52(3):235–42.View ArticlePubMedGoogle Scholar
- Brussow H, Canchaya C, Hardt WD. Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion. Microbiol Mol Biol Rev. 2004;68(3):560–602. table of contentsView ArticlePubMedPubMed CentralGoogle Scholar
- Nanda AM, Thormann K, Frunzke J. Impact of spontaneous prophage induction on the fitness of bacterial populations and host-microbe interactions. J Bacteriol. 2015;197(3):410–9.View ArticlePubMedPubMed CentralGoogle Scholar
- De Paepe M, Leclerc M, Tinsley CR, Petit MA. Bacteriophages: an underestimated role in human and animal health? Front Cell Infect Microbiol. 2014;4:39.View ArticlePubMedPubMed CentralGoogle Scholar
- Manrique P, Dills M, Young MJ. The human gut phage community and its implications for health and disease. Viruses. 2017;9(6)Google Scholar
- Kim MS, Park EJ, Roh SW, Bae JW. Diversity and abundance of single-stranded DNA viruses in human feces. Appl Environ Microbiol. 2011;77(22):8062–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Manrique P, Bolduc B, Walk ST, van der Oost J, de Vos WM, Young MJ. Healthy human gut phageome. Proc Natl Acad Sci U S A. 2016;113(37):10400–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Norman JM, Handley SA, Baldridge MT, Droit L, Liu CY, Keller BC, Kambal A, Monaco CL, Zhao G, Fleshner P, et al. Disease-specific alterations in the enteric virome in inflammatory bowel disease. Cell. 2015;160(3):447–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Ren J, Ahlgren NA, Lu YY, Fuhrman JA, Sun F. VirFinder: a novel k-mer based tool for identifying viral sequences from assembled metagenomic data. Microbiome. 2017;5(1):69.View ArticlePubMedPubMed CentralGoogle Scholar
- Benevides L, Burman S, Martin R, Robert V, Thomas M, Miquel S, Chain F, Sokol H, Bermudez-Humaran LG, Morrison M, et al. New insights into the diversity of the genus Faecalibacterium. Front Microbiol. 2017;8:1790.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: a fast phage search tool. Nucleic Acids Res. 2011;39(Web Server issue):W347–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, Wishart DS. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44(W1):W16–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Kiss J, Szabo M, Olasz F. Site-specific recombination by the DDE family member mobile element IS30 transposase. Proc Natl Acad Sci U S A. 2003;100(25):15000–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Lavigne R, Seto D, Mahadevan P, Ackermann HW, Kropinski AM. Unifying classical and molecular taxonomic classification: analysis of the Podoviridae using BLASTP-based tools. Res Microbiol. 2008;159(5):406–14.View ArticlePubMedGoogle Scholar
- Adriaenssens E, Brister JR. How to name and classify your phage: an informal guide. Viruses. 2017;9(4)Google Scholar
- Minot S, Grunberg S, Wu GD, Lewis JD, Bushman FD. Hypervariable loci in the human gut virome. Proc Natl Acad Sci U S A. 2012;109(10):3962–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Heler R, Marraffini LA, Bikard D. Adapting to new threats: the generation of memory by CRISPR-Cas immune systems. Mol Microbiol. 2014;93(1):1–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lopes A, Tavares P, Petit MA, Guerois R, Zinn-Justin S. Automated classification of tailed bacteriophages according to their neck organization. BMC Genomics. 2014;15:1027.View ArticlePubMedPubMed CentralGoogle Scholar
- Hulo C, Masson P, Le Mercier P, Toussaint A. A structured annotation frame for the transposable phages: a new proposed family “Saltoviridae” within the Caudovirales. Virology. 2015;477:155–63.View ArticlePubMedGoogle Scholar
- Toussaint A. Transposable Mu-like phages in Firmicutes: new instances of divergence generating retroelements. Res Microbiol. 2013;164(4):281–7.View ArticlePubMedGoogle Scholar
- Delattre H, Souiai O, Fagoonee K, Guerois R, Petit MA. Phagonaute: a web-based interface for phage synteny browsing and protein function prediction. Virology. 2016;496:42–50.View ArticlePubMedGoogle Scholar
- Guo H, Arambula D, Ghosh P, Miller JF. Diversity-generating retroelements in phage and bacterial genomes. Microbiol Spectr. 2014;2(6)Google Scholar
- Liu M, Deora R, Doulatov SR, Gingery M, Eiserling FA, Preston A, Maskell DJ, Simons RW, Cotter PA, Parkhill J, et al. Reverse transcriptase-mediated tropism switching in Bordetella bacteriophage. Science. 2002;295(5562):2091–4.View ArticlePubMedGoogle Scholar
- Paul BG, Burstein D, Castelle CJ, Handa S, Arambula D, Czornyj E, Thomas BC, Ghosh P, Miller JF, Banfield JF, et al. Retroelement-guided protein diversification abounds in vast lineages of bacteria and archaea. Nat Microbiol. 2017;2:17045.View ArticlePubMedPubMed CentralGoogle Scholar
- Schillinger T, Zingler N. The low incidence of diversity-generating retroelements in sequenced genomes. Mob Genet Elements. 2012;2(6):287–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Barr JJ, Auro R, Furlan M, Whiteson KL, Erb ML, Pogliano J, Stotland A, Wolkowicz R, Cutting AS, Doran KS, et al. Bacteriophage adhering to mucus provide a non-host-derived immunity. Proc Natl Acad Sci U S A. 2013;110(26):10771–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Allet B. Mu insertion duplicates a 5 base pair sequence at the host inserted site. Cell. 1979;16(1):123–9.View ArticlePubMedGoogle Scholar
- Garneau JR, Depardieu F, Fortier LC, Bikard D, Monot M. PhageTerm: a tool for fast and accurate determination of phage termini and packaging mechanism using next-generation sequencing data. Sci Rep. 2017;7(1):8292.View ArticlePubMedPubMed CentralGoogle Scholar
- Weitz JS, Poisot T, Meyer JR, Flores CO, Valverde S, Sullivan MB, Hochberg ME. Phage-bacteria infection networks. Trends Microbiol. 2013;21(2):82–91.View ArticlePubMedGoogle Scholar
- de Leeuw M, Baron M, Brenner A, Kushmaro A. Genome analysis of a novel broad host range Proteobacteria phage isolated from a bioreactor treating industrial wastewater. Genes (Basel). 2017;8(1)Google Scholar
- Malki K, Kula A, Bruder K, Sible E, Hatzopoulos T, Steidel S, Watkins SC, Putonti C. Bacteriophages isolated from Lake Michigan demonstrate broad host-range across several bacterial phyla. Virol J. 2015;12:164.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim S, Ryu K, Biswas D, Ahn J. Survival, prophage induction, and invasive properties of lysogenic Salmonella typhimurium exposed to simulated gastrointestinal conditions. Arch Microbiol. 2014;196(9):655–9.View ArticlePubMedGoogle Scholar
- De Paepe M, Tournier L, Moncaut E, Son O, Langella P, Petit MA. Carriage of lambda latent virus is costly for its bacterial host due to frequent reactivation in Monoxenic mouse intestine. PLoS Genet. 2016;12(2):e1005861.View ArticlePubMedPubMed CentralGoogle Scholar
- Diard M, Bakkeren E, Cornuault JK, Moor K, Hausmann A, Sellin ME, Loverdo C, Aertsen A, Ackermann M, De Paepe M, et al. Inflammation boosts bacteriophage transfer between Salmonella spp. Science. 2017;355(6330):1211–5.View ArticlePubMedGoogle Scholar
- Bobay LM, Touchon M, Rocha EP. Manipulating or superseding host recombination functions: a dilemma that shapes phage evolvability. PLoS Genet. 2013;9(9):e1003825.View ArticlePubMedPubMed CentralGoogle Scholar
- De Paepe M, Hutinet G, Son O, Amarir-Bouhram J, Schbath S, Petit MA. Temperate phages acquire DNA from defective prophages by relaxed homologous recombination: the role of Rad52-like recombinases. PLoS Genet. 2014;10(3):e1004181.View ArticlePubMedPubMed CentralGoogle Scholar
- Touchon M, Bernheim A, Rocha EP. Genetic and life-history traits associated with the distribution of prophages in bacteria. ISME J. 2016;10(11):2744–54.Google Scholar
- Martinez-Hernandez F, Fornas O, Lluesma Gomez M, Bolduc B, de la Cruz Pena MJ, Martinez JM, Anton J, Gasol JM, Rosselli R, Rodriguez-Valera F, et al. Single-virus genomics reveals hidden cosmopolitan and abundant viruses. Nat Commun. 2017;8:15892.View ArticlePubMedPubMed CentralGoogle Scholar
- Krumsiek J, Arnold R, Rattei T. Gepard: a rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics. 2007;23(8):1026–8.View ArticlePubMedGoogle Scholar
- Biswas A, Staals RH, Morales SE, Fineran PC, Brown CM. CRISPRDetect: a flexible algorithm to define CRISPR arrays. BMC Genomics. 2016;17:356.View ArticlePubMedPubMed CentralGoogle Scholar
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Reyes A, Blanton LV, Cao S, Zhao G, Manary M, Trehan I, Smith MI, Wang D, Virgin HW, Rohwer F, et al. Gut DNA viromes of Malawian twins discordant for severe acute malnutrition. Proc Natl Acad Sci U S A. 2015;112(38):11941–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Milne I, Stephen G, Bayer M, Cock PJ, Pritchard L, Cardle L, Shaw PD, Marshall D. Using Tablet for visual exploration of second-generation sequencing data. Brief Bioinform. 2013;14(2):193–202.View ArticlePubMedGoogle Scholar
- Chassaing B, Srinivasan G, Delgado MA, Young AN, Gewirtz AT, Vijay-Kumar M. Fecal lipocalin 2, a sensitive and broadly dynamic non-invasive biomarker for intestinal inflammation. PLoS One. 2012;7(9):e44328.View ArticlePubMedPubMed CentralGoogle Scholar