Comparative Metagenome-Assembled Genome Analysis of “Candidatus Lachnocurva vaginae”, Formerly Known as Bacterial Vaginosis-Associated Bacterium−1 (BVAB1)

Introduction

Bacterial vaginosis (BV) is a common vaginal infection affecting approximately 30% of US reproductive-aged women, with both African- and Mexican-Americans disproportionately afflicted (Allsworth and Peipert, 2007; Koumans et al., 2007). Though treatable with antibiotics, BV has a high rate of recurrence (Bradshaw et al., 2006; Schwebke and Desmond, 2007). The microbiological diagnosis of BV is established using Nugent scoring of Gram-stained vaginal smears and is defined by a low abundance of Lactobacillus spp. morphotypes and a wide array of strict and facultative Gram-negative anaerobes (Nugent et al., 1991). The clinical diagnosis of BV is made when 3 of the 4 Amsel's criteria are met: vaginal pH > 4.5, homogenous vaginal discharge, a positive whiff test, and the presence of clue cells upon wet mount examination (Amsel et al., 1983). Aside from the burdensome symptoms of vaginal discharge and fishy odor, BV is also associated with increased risk to adverse health outcomes including preterm birth (Leitich et al., 2003), increased risk of sexually transmitted infections acquisition and transmission, including HIV (Taha et al., 1998; Cherpes et al., 2003; Ness et al., 2005; Atashili et al., 2008; Cohen et al., 2012) and pelvic inflammatory disease (Ness et al., 2004).

A critical step along the path to understanding the ecology and pathogenic potential of a bacterial species is the characterization of its genome. Yet many of the bacteria associated with BV have thus far been uncultivable, further complicating genome sequencing efforts. BV-associated bacterium 1 (BVAB1) is one such organism for which there is limited genetic information available. BVAB1 was first identified by Fredricks et al. (2005) using 16S rRNA gene amplicons Sanger sequencing of samples associated with BV, and has eluded cultivation efforts since. Based on the sequence of its 16S rRNA gene, BVAB1 belongs to the family Lachnospiraceae (Muzny et al., 2014) and has often been misidentified as belonging to the genus Shuttleworthia (Lamont et al., 2011; Petrova et al., 2013). Interestingly, Gram-negative curved rods designated Mobiluncus morphotypes on Gram stain in Nugent scoring have been shown to likely be BVAB1 (Srinivasan et al., 2013). Further, vaginal communities in which BVAB1 16S rRNA gene sequence is detected have been associated with vaginal inflammation and persistent BV in African women (Lennard et al., 2018). BVAB1 remains uncultured and aside from detection of this species via partial 16S rRNA gene amplicon sequencing, little is known about its metabolism, pathogenic potential, or ecology in the vaginal environment, especially during BV. Further understanding of the genetic and physiological properties of BVAB1 will help to dissect the complex etiology of BV. Previously, a 94 contig BVAB1 metagenome-assembled genome from short-read sequencing was produced (Fettweis et al., 2019).

In this study, we characterize the first circularized metagenome assembled genome (cMAG) of BVAB1 constructed from a metagenome sequenced using the PacBio Sequel II long read platform. In addition, we compare this cMAG to an additional six metagenome-assembled genomes (MAGs). All genomes originate from different women with symptomatic or asymptomatic BV. Based on phylogenetic analysis of full-length 16S rRNA gene sequences obtained from the genomic assemblies, we propose to rename this bacterium “Candidatus Lachnocurva vaginae”.

Materials and Methods

Sample Collection

Vaginal samples used in this study were identified as containing a high relative abundance (> 60%) of BVAB1 using 16S rRNA gene amplicon sequencing of the V3–V4 regions as previously reported (Holm et al., 2019). Cervicovaginal lavages from six participants that were collected as part of the NIH Longitudinal Study for Vaginal Flora (LSVF) (Klebanoff et al., 2004) by washing the vaginal walls with 3 mL sterile, deionized water, aspiration from the vaginal vault via pipette and stored at −80°C, were included in this study. Gram stain smears were prepared for Nugent scoring as previously described (Nugent et al., 1991). DNA was extracted from 200 μL of lavage fluid using the MagAttract Microbial DNA Kit (QIAGEN Inc., Germantown MD) automated on a Hamilton Star robotic platform according the manufacturer recommendations. DNA was eluted in a final volume of 110 μL nuclease-free water.

An additional swab sample collected as part of the UMB-HMP study was used in this study (Ravel et al., 2013). The swab was self-collected by a participant using a Copan Eswab re-suspended into 1 mL Amies transport medium (ESwab, Copan Diagnostics Inc.), frozen at −20°C for no more than a week, and then transferred to −80°C until analyzed. High molecular weight DNA was extracted from this sample using the MasterPure DNA purification kit (Lucigen) with two phenol/chloroform cleanups prior to DNA precipitation. DNA extraction was quantified on a TapeStation 2200 instrument run with a Genomic DNA tape (Agilent).

Metagenomic Library Construction and Sequencing on the PacBio Sequel II Platform

The extracted DNA from the swab collected as part of the UMB-HMP study (Ravel et al., 2013) was found to be of sufficient concentration (5.49 ng/μL in 200 μL) for long-read sequencing using the Pacific Biosciences Sequel II platform (Pacific Biosciences). The sequencing library was prepared with SMRTBell Template Prep Kit 1.0 and was size-selected on a BluePippen (Sage Science) with a cutoff of 5 kb. The library was barcoded and sequenced as part of a multiplexed run with four other unrelated samples. Sequencing was performed on a PacBio Sequel II instrument with an 8M cell loaded at 60 pM at the Genomic Resource Center of the University Maryland School of Medicine.

Long Read Quality Filtering, Host-Read Removal, and cMAG Construction

Raw reads were demultiplexed with lima (version 1.9.0) using default parameters except for minimum barcode score set at 26 and a minimum read length of 50 bp after clipping of the barcode was enforced. Both tools are part of the SMRTLink 6.0.1 software package with updated CCS version 3.4.1. Human reads were detected using pbalign v0.4.1 (Pacific Biosciences) and the human genome build 38 (GRCh38.p12). Remaining reads were corrected and assembled via Canu v1.7 and the “-pacbio-raw” protocol (Koren et al., 2017). The largest resulting contig was ca. 1.6 Mb in length, much larger than the second longest contig (ca. 600 kb). Four copies of the 16S rRNA gene were detected on this contig and were identical to the existing BVAB1 16S rRNA gene reference AY724739 Fredricks et al., 2005 (NEJM). The contig had 6.5 kb of overlapping ends and was determined to be circular by Canu. This contig was then manually circularized using Geneious version 2019.2.1 (Galens et al., 2011), searched for the dnaA gene, and rotated so that dnaA was the first gene (Kearse et al., 2012). The circularized metagenome-assembled genome (cMAG) was annotated using the IGS Prokaryotic Annotation Pipeline (Galens et al., 2011). Translated gene sequences were assessed against KEGG using the BlastKOALA algorithm for insights into “Candidatus Lachnocurva vaginae” metabolism (Kanehisa et al., 2016). Bacteriophages were detected using PHASTER (Zhou et al., 2011; Arndt et al., 2016). Completeness and contamination of the cMAG was analyzed using CheckM version 1.0.18 (Parks et al., 2015) and the taxonomy_wf flag specifying Order Clostridiales and rerun specifying Family Lachnospiraceae. Metabolic reconstruction was examined with the cMAG using the metabolic modeling function in KBase (Arkin et al., 2018). Complete media for gapfilling and a Gram-positive template were used.

Metagenomic Library Construction and Sequencing on the Illumina HiSeq 4000 Platform

Metagenomic libraries for the six samples from the LSVF study were prepared using the KAPA HyperPlus Kit (Kapa Biosystems) with KAPA Single-Indexed Adapter Kit Set B. A fixed volume (35 μL) of genomic DNA was used as input, and libraries were prepared following the manufacturer's protocol with modifications based on their amount of input DNA as in Supplemental Data Sheet 1. For samples with 0.5 or 0.2 ng input DNA, the fragmentation enzyme was diluted 1:2 or 1:5 with water. All samples were fragmented at 37°C for 5 min. Adapter concentrations varied according to the input DNA as listed in Supplemental Data Sheet 1, and the adapter ligation was carried out overnight at 4°C for all samples. The post-ligation cleanup was performed with 0.8x Ampure XP beads (Beckman Coulter, Indianapolis IN) and 20 μL of sample was used in library amplification. Amplification library cycles varied by input DNA as listed in Supplemental Data Sheet 1. Post-amplification cleanup was performed with 1x Ampure XP beads; libraries with remaining adapter dimer peaks were cleaned a second time. The final elution was in 25 μL of nuclease-free water. Libraries were run on a TapeStation instrument with a D1000 tape (Agilent) to assess quality and concentration. Libraries were sequenced (8 libraries/lane) on an Illumina HiSeq 4000 instrument using the 150 bp paired-end protocol.

Short Read Quality Filtering, Host-Read Removal, and Metagenomic-Assembled Genome Reconstruction

Metagenomic reads were quality filtered using Trimmomatic v0.36 (Bolger et al., 2014) to remove sequencing adapters allowing for 2 mismatches, a palindromic clip threshold of 30, and a simple clip threshold of 10. Bases with quality scores < 3 were removed from the beginning and end of reads combined with a 4 bp sliding window which trimmed a read if the average quality score within that window fell below 15. Reads < 75 bp in length were removed. Human reads were detected by mapping to the human genome build 38 (GRCh38.p12) with Bowtie 2 v2.3.4.1 and default settings (Langmead and Salzberg, 2012), and removed using samtools v1.9 (Li et al., 2009) and bedtools v2.27.1 (Quinlan and Hall, 2010; Quinlan, 2014) (see Supplemental Data Sheet 2 for specific scripts). Metagenomic assemblies were produced using SPAdes genome assembler v3.13.0 (Bankevich et al., 2012) with the careful setting. Resulting contigs were aligned to the cMAG using NUCmer with the –mum setting and allowing for 5,000 bp breaklength, filtered to remove aligned contigs < 100 bp, and minimum contig coverage of 25%. MAGs were annotated with the Live Annotate & Predict tool from Geneious version 2019.2.1 (Kearse et al., 2012) using the “Ca. Lachnocurva vaginae” cMAG as annotation source (sequence similarity of ≥ 95% required). Genes with no similarity to the cMAG were annotated with Prokka v1.13 (Seemann, 2014). Phages were detected using PHASTER (Zhou et al., 2011). A circle plot was constructed with the BLAST Ring Image Generator v0.95 (Alikhan et al., 2011), and the annotated cMAG.

Metabolomic Analysis

Metabolomics analysis was conducted as previously described (Nelson et al., 2018) by Metabolon Inc. using 200 μL of each LSVF lavage sample used in this study or 200 μL of a frozen dry swab eluted in 1 ml of PBS for the sample from the UMB-HMP study (Ravel et al., 2013). The abundances of 561 compounds were quantified using Ultrahigh Performance Liquid Chromatography-Tandem Mass Spectroscopy (UPLC-MS/MS). Quantities were corrected for instrument block variability and reported as normalized area-under-the-curve estimates. Figures were generated using ggplot2 v3.2.0 (Wickham, 2009).

Comparative Genome and Phylogenetic Analyses

Average nucleotide identities (ANI) were calculated using FastANI v1.1 with fragment lengths of 1,000 bp, which is about the mean length of coding sequences (Jain et al., 2018). Maximum-likelihood phylogenetic trees of full-length 16S rRNA gene sequences alignments were generated with MUSCLE v3.8.425 using default parameters (Edgar, 2004) and FastTree v2.1.11 using the Generalized Time Reversible model and default parameters (Price et al., 2009, 2010) in Geneious v11.0.3 (Kearse et al., 2012). The tree was rooted with Fusobacterium nucleatum (AJ133496) and Propionigenium modestum (X54275) (Domingo et al., 2009), and members of the genus Lachnoclostridium were included as neighbors (Yutin and Galperin, 2013) as well as 16S rRNA genes from Shuttleworthia satelles (NR_028827), and Lachnospiraceae bacterium 2_1_46FAA (NR_025127). In addition, the cMAG was submitted for taxonomic placement using the TrueBacID tool from EZBioCloud (Yoon et al., 2017). cMAG and MAG pseudomolecules were aligned and visualized using AliTV v1.0.6 (Ankenbrand et al., 2017). Genes were ordered by gene synteny of the cMAG. To examine homology to other reference sequences, the cMAG was aligned to the NCBI protein reference database (downloaded 1/28/2019), as well as the genomes of Shuttleworthia satelles (NZ_ACIP00000000), Lachnobacterium bovis (GCF_900107245.1), and Lachnospiraceae bacterium 2_1_46FAA (ADLB02000001.1) using BLAST v2.8.1+. The top 3 hits were chosen by lowest e-value, highest percent identity, and longest alignment length, in that order. Genomic and pathogenicity islands were explored using IslandViewer4 (Bertelli et al., 2017).

Results and Discussion

Participant Information

All six participants in the LSVF study were of reproductive age and had high Nugent scores (6–10) (Table 1). Four were diagnosed with asymptomatic Amsel-BV, one with symptomatic Amsel-BV, and one was negative for Amsel-BV (Mckinnon et al., 2019). The woman who participated in the UMB-HMP study (Ravel et al., 2013) had a Nugent score of 8 and was diagnosed with asymptomatic Nugent-BV. The vaginal microbiota of these 7 samples had >60% BVAB1 as defined by 16S rRNA gene V3-V4 amplicon sequencing (Ravel et al., 2013).

TABLE 1

Table 1. Participants demographics and cervicovaginal lavage microbial compositions for samples used in metagenomic reconstruction of the “Ca. Lachnocurva vaginae” circularized metagenome-assembled genome using long-read (*) and shotgun sequencing (*).

Phylogenetic Analysis

A phylogenetic analysis using the full-length 16S rRNA genes extracted from all MAGs revealed BVAB1 belongs in the Clostridales family Lachnospiraceae, and that S. satelles is the closest known relative (Figure 1), though nucleotide identity is only 89.2%. Similar results were obtained from EZBioCloud's TrueBacID (Supplemental Data Sheet 3). We therefore propose a new candidate species for BVAB1: “Candidatus Lachnocurva vaginae”, which represents the phylogenetic placement, the curved morphology of the cells, and the source of this cMAG. Considering the entire cMAG, the ANI between “Ca. Lachnocurva vaginae” and the S. satelles draft genome (NZ_ACIP00000000) was 81.85%. Relative to each other, “Ca. Lachnocurva vaginae” 16S rRNA genes were >99.7% identical and the cMAG and MAG average nucleotide identities were also high (98.6–99.2%, Table 2).

FIGURE 1

Figure 1. “Ca. Lachnocurva vaginae” is related to Shuttleworthia satelles; the full-length 16S rRNA gene sequences are 89–90% identical.

TABLE 2

Table 2. Average Nucleotide Identity (ANI) between the “Ca. Lachnocurva vaginae” circularized metagenome-assembled genome, each metagenome-assembled genome, and the Shuttleworthia satelles draft genome NZ_ACIP00000000.

Overall Genomic Features

The cMAG of “Ca. Lachnocurva vaginae” is 1,649,642 bp in size with 31.8% GC content, encodes 1,578 genes and was estimated to be 99.1% complete and 0% contaminated at the order level and 90.9% complete and 0% contaminated at the family level. The relatively low completion estimate at the family level was due to 22 missing marker genes out of 312 (Supplemental Table 2). Mean coverage of the cMAG assembly by long reads was 124X (range: [26–285]). The “Ca. Lachnocurva vaginae” MAGs constructed using the short-read Illumina HiSeq 4000 metagenomic sequence data were 1.48–1.62 Mb in size (mean: 1.57 Mb) with a number of contigs ranging from 26 to 152 (mean N50: 99,943 bp) and 31.6% GC (Table 3). Genome annotation indicated the presence of four complete rRNA operon copies (Figure 2) and 42 tRNA genes in the cMAG. As expected, partial rRNA operons were identified in the Illumina-based MAGs (Supplemental Image 1).

TABLE 3

Table 3. “Ca. Lachnocurva vaginae” MAG assembly characteristics.

FIGURE 2

Figure 2. Circle plot of the “Ca. Lachnocurva vaginae” circularized metagenome-assembled genome (cMAG). Circles are from inside to outside: Circle 1: cMAG position; Circle 2: GC content, 3: GC skew, 4: forward-direction coding sequences, 5: reverse-direction coding sequences, and 6: phage detected in cMAG (green bars) and rRNA operons (red bars).

Genomic Features of “Ca. Lachnocurva vaginae”

Metabolic modeling and reconstruction of the “Ca. Lachnocurva vaginae” cMAG contained 122 genes, 480 reactions, 581 metabolites, and no mass imbalance was found. Of the 1,578 genes detected, 255 genes had best blast matches of 80% sequence identity covering > 80% of the query gene. Of these, 165 matches were within the order Clostridiales, and 90 were not (Supplemental Image 2). Genes encoding transporters for mannose, fructose, and L-ascorbate were identified in all MAGs (Figure 3 and Supplemental Table 1, LCVA_199-201, LCVA_1491-1494), as were complete pathways for glycolysis, pyruvate oxidation, and the non-oxidative phase of the pentose phosphate pathway. Mannose was noticeably absent, or below the limit of detection, in the metabolome of all the samples, but present in other BV-like samples that did not contain a high relative abundance of the candidate species (Figure 4). We found “Ca. Lachnocurva vaginae” to have the genetic capability to uptake and metabolize mannose (LCVA_1184), suggesting mannose could also be a carbohydrate source for the candidate species (Figure 3).

FIGURE 3

Figure 3. Genomic features of “Ca. Lachnocurva vaginae”. Genes coding for the following functions were observed: Methyl-accepting chemotaxis (MCP) and subsequent flagella assembly, Sec-SRP secretion systems, choline import, and metabolism to trimethylamine (TMA), bacteriocin production and export, mannose, fructose, and L-ascorbate transport and metabolism, and D-lactate dehydrogenase.

FIGURE 4

Figure 4. Metabolomic signatures for biochemicals of interest in vaginal community state type (CST) IV-A (high abundance of “Ca. Lachnocurva vaginae”), CST I (dominated by Lactobacillus crispatus), and CST IV-B (high abundance of Gardnerella vaginalis).

The D-lactate dehydrogenase gene (LCVA_41), but not an L-lactate dehydrogenase gene was also observed in all MAGs. This result was unexpected, as the production of D-lactate in the vaginal environment is thought to be a key and somewhat unique feature of certain Lactobacillus spp. (Witkin et al., 2013; France et al., 2016). Metabolomic data for each sample were explored for the presence of lactate to provide evidence supporting this genomic finding. While lactate was observed in most samples, its abundance was substantially lower than that from representative samples that were dominated by Lactobacillus crispatus, a known D-lactate producer (Figure 4). Thus, either the enzyme exhibits lower activity in “Ca. Lachnocurva vaginae” or is instead used in the reverse reaction to consume Lactobacillus spp.-produced D-lactate to produce pyruvate. This is a strategy described for another vaginal anaerobic Gram-negative coccus, Veillonella parvula, which is able to grow on lactate as the sole source of carbon (Gronow et al., 2010), and would possibly contribute to an increase in vaginal pH, thereby creating a more favorable growth environment. However, for lactate to convert to pyruvate in anaerobic bacteria, an electron-bifurcating complex with an electron transfer flavoprotein alpha and beta subunit (EtfAB) would be necessary to make the reaction favorable (Weghoff et al., 2015). An EtfAB homolog was not observed in the “Ca. Lachnocurva vaginae” cMAG, suggesting that another mechanism may be involved. Instead, succinate was the most abundant short chain fatty acid in the communities from which “Ca. Lachnocurva vaginae” MAGs originated (Figure 4). Not surprisingly, succinate has been associated with bacterial vaginosis, a condition in which “Ca. Lachnocurva vaginae” is often found (Srinivasan et al., 2015). Interestingly, all “Ca. Lachnocurva vaginae” MAGs lacked genes for all steps of the tricarboxylic acid cycle, indicating it is unable to produce succinate via this pathway. It is possible that “Ca. Lachnocurva vaginae” produces metabolite(s) other species in the community can convert into succinate or produces it via a yet to be determined pathway. While “Ca. Lachnocurva vaginae” does not seem to be able to produce succinate, it has the genetic machinery to produce acetate from pyruvate (i.e., ackA and pta LCVA_939 and LCVA_1006). We were, however, unable to detect acetate in the metabolome through the methods used, as it is too small of a molecule.

A full suite of genes required for flagella assembly was observed in all MAGs, as well as genes required for methyl-accepting chemotaxis (LCVA_205, LCVA_992, LCVA_1053, LCVA_1122) and for the downstream signal process that mediates flagellar response, cheY (LCVA_362, LCVA_688) (Bren and Eisenbach, 2000). Flagella have yet to be visually observed on “Ca. Lachnocurva vaginae” and were not observed in the source samples used for this study (Figure 5).

FIGURE 5

Figure 5. Gram stains of LSVF samples containing >70% “Ca. Lachnocurva vaginae” and used in this study. Morphologically “Ca. Lachnocurva vaginae” is a curved rod. Images are labeled with the sample ID.

We found the cMAG and some MAGs to also include genes that encode protection against antibiotics including the tetO gene (LCVA_1310, also observed in MAG Y3255) and drug efflux pumps for macrolides (macB LCVA_500, also observed in MAGs Y2337 and Y2266) and fluoroquinolones (LCVA_1202, present in all MAGs). A gene for hemolysin, tlyA, was also observed in all MAGs (LCVA_1031), suggesting direct interaction between “Ca. Lachnocurva vaginae” and host tissues.

Additionally, choline transporters, as well as the cutC and cutD genes were detected in all MAGs (LCVA_1194 and LCVA_1193, respectively). The cut genes metabolize choline and produce trimethylamine (TMA) (Martinez-Del Campo et al., 2015). Aside from transport of exogenous choline, another potential source of choline would be via phospholipase D hydrolysis of phosphotidyl choline into choline as seen in gut bacteria (Chittim et al., 2019), however no genes encoding such activity were found in the cMAG of “Ca. Lachnocurva vaginae”. Comparing the translated cutC amino acid sequence from “Ca. Lachnocurva vaginae” to those reported by Martinez-Del Campo et al. (2015), we observed three of the five active sites conserved in “Ca. Lachnocurva vaginae” (Cys489, Glu491, Gly821, data not shown) indicating that it is likely functional. TMA is one of the substances believed to be responsible for the fishy odor associated with BV (Brand and Galask, 1986; Wolrath et al., 2002), however it is rarely detected in metabolomics analyses as it is highly volatile, unless it is performed on freshly collected samples (Wolrath et al., 2002).

Two intact bacteriophages were detected in the cMAG (Figure 2). Partial matches to the same phages were observed in all MAGs (Supplemental Image 1). Best BLAST hits indicated that both bacteriophages belong to the Siphoviridae family of double-stranded viruses (NC_009552, NC_011167) which can exhibit both lytic and lysogenic phases. Bacteriophages of this family have previously been reported in the vaginal species Lactobacillus jensenii (Martin et al., 2010). Genomic islands were detected at multiple sites (Supplemental Table 1 and Figure 2). Analysis of proteins encoded in the largest island (65.5 kb, coordinates 1,331,959–1,391,037, LCVA_1280-1330) showed similarity to proteins from Shuttleworthia satelles, and other taxa from the Clostridiales (Supplemental Table 3 and Supplemental Image 2) and several transposons and integrases, as well as tetracycline resistance proteins (TetO, LCVA_1310), and ABC transporters. We ruled out the presence of genomic islands due to mis-assembly by analyzing the long-read coverage spanning the 5′ and 3′ ends of genomic islands. More than 100 long PacBio reads of a means size of 9 kb span the junctions of all detected genomic islands. Further, mean coverage of the genomic islands was similar to that of the mean cMAG coverage of 124X. Portions of these islands were observed in the other MAGs assembled in this study (Supplemental Image 1 and Supplemental Table 1). The stringency of the mapping would not recover reads to regions of sequence diversity, thus these results may indicate this “Ca. Lachnocurva vaginae” likely contains multiple regions of genetic fluidity or diversity. However, the missing regions in MAGs may also be an artifact of metagenomic assembly.

Conclusion

We present here a circularized MAG of “Ca. Lachnocurva vaginae” and six MAGs of the candidate species, previously known as BVAB1, an important member of the human vaginal microbiota associated with bacterial vaginosis and other adverse outcomes. Short-read metagenomic assembles do not perform well and lead to sub-optimal assemblies with missing regions, whereas long read metagenomic assemblies are promising and can generate circularized metagenome assembled genomes, as shown in this study. Our inability to culture this bacterium has limited our understanding of its ecological role in the vaginal environment and its relation to women's health. We have shown that “Ca. Lachnocurva vaginae” has the genomic potential for motility and chemotaxis, and is likely capable of resisting several antibiotics via drug efflux systems. Our analysis indicates this candidate species may contribute to the fishy odor characteristic of bacterial vaginosis through the production of TMA from choline. This crucial genomic data could be used in metabolic modeling experiments to define a culture medium suitable for the cultivation of “Ca. Lachnocurva vaginae”, a critical step to further understand its role in the vaginal microbiome.

Data Availability Statement

Sequence data have been submitted to NCBI under BioProject PRJNA562728.

Ethics Statement

Samples used in this study were archived and de-identified cervicovaginal lavages and swabs. The samples were originally collected after obtaining informed consent by all participants, who also provided consent for storage of the samples and use in future research studies related to women's health. The original study was approved by the University of Maryland School of Medicine Institutional Review Board.

Author Contributions

EM and CR performed the extraction and library preparations from the raw samples. JH, AM, LT, MF, and BM performed bioinformatics analyses to generate the MAGs. JH, RB, and JR designed the study. JH, MF, BM, and JR executed the study, contributed to the analyses. JH, RB, MF, BM, and JR wrote and edited the manuscript.

Funding

The research reported in this publication was supported in part by the National Institute of Allergy and Infectious Diseases of the National Institutes of Health under award numbers F32AI136400 (to JH), U19AI084044, R01AI116799, R01NR015495, and the Bill & Melinda Gates Foundation award OPP1189217.

Conflict of Interest

JR is co-founder of LUCA Biologics, a biotechnology company focusing on translating microbiome research into live biotherapeutics drugs for women's health.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors are grateful to the University of Maryland School of Medicine Genome Resource Center and Elizabeth O'Hanlon for providing the Gram-stained images. We are also grateful to Dr. Benjamin Woolston for critical reading of the manuscript and suggestions on the mechanisms for TMA production and D-lactate metabolism.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.00117/full#supplementary-material

Supplemental Data Sheet 1. Metagenomic library preparation modifications based on input amounts of DNA.

Supplemental Data Sheet 2. Processing scripts for analysis.

Supplemental Data Sheet 3. Taxonomic and phylogenetic placement of “Ca. Lachnocurva vaginae” (UAB071) from EZBioCloud TrueBacID system.

Supplemental Table 1. List of genes annotated in the “Ca. Lachnocurva vaginae” cMAG. Excludes rRNA and tRNA.

Supplemental Table 2. List of marker genes used in Family-level CheckM analysis and presence (1) or absence (0) in the “Ca. Lachnocurva vaginae” cMAG.

Supplemental Table 3. Top 3 best BLAST matches to the genes of the “Ca. Lachnocurva vaginae” cMAG.

Supplemental Image 1. Visualization of “Ca. Lachnocurva vaginae” whole genome alignments of cMAG and other MAGs.

Supplemental Image 2. Circle plot of “Ca. Lachnocurva vaginae” cMAG with best BLAST hits indicated as either Clostridiales or not (Ring 5). If hit was within Clostridiales and >80% identical to cMAG, the specific taxon of close relatives was indicated in Ring 4. Included BLAST hits covered at least 80% of the query gene.

References

Alikhan, N. F., Petty, N. K., Ben Zakour, N. L., and Beatson, S. A. (2011). BLAST ring image generator (BRIG): simple prokaryote genome comparisons. BMC Genomics 12:402. doi: 10.1186/1471-2164-12-402

PubMed Abstract | CrossRef Full Text | Google Scholar

Allsworth, J. E., and Peipert, J. F. (2007). Prevalence of bacterial vaginosis: 2001-2004 national health and nutrition examination survey data. Obstet. Gynecol. 109, 114–120. doi: 10.1097/01.AOG.0000247627.84791.91

PubMed Abstract | CrossRef Full Text | Google Scholar

Amsel, R., Totten, P. A., Spiegel, C. A., Chen, K. C., Eschenbach, D., and Holmes, K. K. (1983). Nonspecific vaginitis. Diagnostic criteria and microbial and epidemiologic associations. Am. J. Med. 74, 14–22. doi: 10.1016/0002-9343(83)91112-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Ankenbrand, M. J., Hohlfeld, S., Hackl, T., and Förster, F. (2017). AliTV—interactive visualization of whole comparisons. PeerJ. Comput. Sci. 3:e116. doi: 10.7717/peerj-cs.116

CrossRef Full Text | Google Scholar

Arkin, A. P., Cottingham, R. W., Henry, C. S., Harris, N. L., Stevens, R. L., Maslov, S., et al. (2018). KBase: the United states department of energy systems biology knowledgebase. Nat. Biotechnol. 36, 566–569. doi: 10.1038/nbt.4163

PubMed Abstract | CrossRef Full Text | Google Scholar

Arndt, D., Grant, J. R., Marcu, A., Sajed, T., Pon, A., Liang, Y., et al. (2016). PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 44, W16–W21. doi: 10.1093/nar/gkw387

PubMed Abstract | CrossRef Full Text | Google Scholar

Atashili, J., Poole, C., Ndumbe, P. M., Adimora, A. A., and Smith, J. S. (2008). Bacterial vaginosis and HIV acquisition: a meta-analysis of published studies. AIDS 22, 1493–1501. doi: 10.1097/QAD.0b013e3283021a37

PubMed Abstract | CrossRef Full Text | Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertelli, C., Laird, M. R., Williams, K. P., Simon Fraser University Research Computing, Group, Lau, B. Y., Hoad, G., et al. (2017). IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res. 45, W30–W35. doi: 10.1093/nar/gkx343

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradshaw, C. S., Morton, A. N., Hocking, J., Garland, S. M., Morris, M. B., Moss, L. M., et al. (2006). High recurrence rates of bacterial vaginosis over the course of 12 months after oral metronidazole therapy and factors associated with recurrence. J. Infect. Dis. 193, 1478–1486. doi: 10.1086/503780

PubMed Abstract | CrossRef Full Text | Google Scholar

Brand, J. M., and Galask, R. P. (1986). Trimethylamine: the substance mainly responsible for the fishy odor often associated with bacterial vaginosis. Obstet. Gynecol. 68, 682–685.

PubMed Abstract | Google Scholar

Bren, A., and Eisenbach, M. (2000). How signals are heard during bacterial chemotaxis: protein-protein interactions in sensory signal propagation. J. Bacteriol. 182, 6865–6873. doi: 10.1128/JB.182.24.6865-6873.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Cherpes, T. L., Meyn, L. A., Krohn, M. A., Lurie, J. G., and Hillier, S. L. (2003). Association between acquisition of herpes simplex virus type 2 in women and bacterial vaginosis. Clin. Infect. Dis. 37, 319–325. doi: 10.1086/375819

PubMed Abstract | CrossRef Full Text | Google Scholar

Chittim, C. L., Martinez Del Campo, A., and Balskus, E. P. (2019). Gut bacterial phospholipase Ds support disease-associated metabolism by generating choline. Nat. Microbiol. 4, 155–163. doi: 10.1038/s41564-018-0294-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, C. R., Lingappa, J. R., Baeten, J. M., Ngayo, M. O., Spiegel, C. A., Hong, T., et al. (2012). Bacterial vaginosis associated with increased risk of female-to-male HIV-1 transmission: a prospective cohort analysis among African couples. PLoS Med. 9:e1001251. doi: 10.1371/journal.pmed.1001251

PubMed Abstract | CrossRef Full Text | Google Scholar

Domingo, M. C., Huletsky, A., Boissinot, M., Helie, M. C., Bernal, A., Bernard, K. A., et al. (2009). Clostridium lavalense sp. nov., a glycopeptide-resistant species isolated from human faeces. Int. J. Syst. Evol. Microbiol. 59, 498–503. doi: 10.1099/ijs.0.001958-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340

PubMed Abstract | CrossRef Full Text | Google Scholar

Fettweis, J. M., Serrano, M. G., Brooks, J. P., Edwards, D. J., Girerd, P. H., Parikh, H. I., et al. (2019). The vaginal microbiome and preterm birth. Nat. Med. 25, 1012–1021. doi: 10.1038/s41591-019-0450-2

PubMed Abstract | CrossRef Full Text | Google Scholar

France, M. T., Mendes-Soares, H., and Forney, L. J. (2016). Genomic comparisons of Lactobacillus crispatus and Lactobacillus iners reveal potential ecological drivers of community composition in the vagina. Appl. Environ. Microbiol. 82, 7063–7073. doi: 10.1128/AEM.02385-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Fredricks, D. N., Fiedler, T. L., and Marrazzo, J. M. (2005). Molecular identification of bacteria associated with bacterial vaginosis. N. Engl. J. Med. 353, 1899–1911. doi: 10.1056/NEJMoa043802

PubMed Abstract | CrossRef Full Text | Google Scholar

Galens, K., Orvis, J., Daugherty, S., Creasy, H. H., Angiuoli, S., White, O., et al. (2011). The IGS standard operating procedure for automated prokaryotic annotation. Stand. Genomic Sci. 4, 244–251. doi: 10.4056/sigs.1223234

PubMed Abstract | CrossRef Full Text | Google Scholar

Gronow, S., Welnitz, S., Lapidus, A., Nolan, M., Ivanova, N., Glavina Del Rio, T., et al. (2010). Complete genome sequence of Veillonella parvula type strain (Te3). Stand. Genomic Sci. 2, 57–65. doi: 10.4056/sigs.521107

PubMed Abstract | CrossRef Full Text | Google Scholar

Holm, J. B., Humphrys, M., Robinson, C. K., Settles, M. L., Ott, S., Fu, L., et al. (2019). Ultra-high throughput multiplexing and sequencing of >500 bp amplicon regions on the Illumina HiSeq 2500 platform. mSystems. 4:e00029-19. doi: 10.1128/mSystems.00029-19

CrossRef Full Text | Google Scholar

Jain, C., Rodriguez, R. L., Phillippy, A. M., Konstantinidis, K. T., and Aluru, S. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 9:5114. doi: 10.1038/s41467-018-07641-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Sato, Y., and Morishima, K. (2016). BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J. Mol. Biol. 428, 726–731. doi: 10.1016/j.jmb.2015.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., et al. (2012). Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649. doi: 10.1093/bioinformatics/bts199

PubMed Abstract | CrossRef Full Text | Google Scholar

Klebanoff, M. A., Schwebke, J. R., Zhang, J., Nansel, T. R., Yu, K. F., and Andrews, W. W. (2004). Vulvovaginal symptoms in women with bacterial vaginosis. Obstet. Gynecol. 104, 267–272. doi: 10.1097/01.AOG.0000134783.98382.b0

PubMed Abstract | CrossRef Full Text | Google Scholar

Koren, S., Walenz, B. P., Berlin, K., Miller, J. R., Bergman, N. H., and Phillippy, A. M. (2017). Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736. doi: 10.1101/gr.215087.116

PubMed Abstract | CrossRef Full Text | Google Scholar

Koumans, E. H., Sternberg, M., Bruce, C., Mcquillan, G., Kendrick, J., Sutton, M., et al. (2007). The prevalence of bacterial vaginosis in the United States, 2001-2004; associations with symptoms, sexual behaviors, and reproductive health. Sex. Transm. Dis. 34, 864–869. doi: 10.1097/OLQ.0b013e318074e565

PubMed Abstract | CrossRef Full Text | Google Scholar

Lamont, R. F., Sobel, J. D., Akins, R. A., Hassan, S. S., Chaiworapongsa, T., Kusanovic, J. P., et al. (2011). The vaginal microbiome: new information about genital tract flora using molecular based techniques. BJOG 118, 533–549. doi: 10.1111/j.1471-0528.2010.02840.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

Leitich, H., Bodner-Adler, B., Brunbauer, M., Kaider, A., Egarter, C., and Husslein, P. (2003). Bacterial vaginosis as a risk factor for preterm delivery: a meta-analysis. Am. J. Obstet. Gynecol. 189, 139–147. doi: 10.1067/mob.2003.339

PubMed Abstract | CrossRef Full Text | Google Scholar

Lennard, K., Dabee, S., Barnabas, S. L., Havyarimana, E., Blakney, A., Jaumdally, S. Z., et al. (2018). Microbial composition predicts genital tract inflammation and persistent bacterial vaginosis in South African adolescent females. Infect. Immun. 86:e00410–17. doi: 10.1128/IAI.00410-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, R., Escobedo, S., and Suarez, J. E. (2010). Induction, structural characterization, and genome sequence of Lv1, a prophage from a human vaginal Lactobacillus jensenii strain. Int. Microbiol. 13, 113–121. doi: 10.2436/20.1501.01.116

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez-Del Campo, A., Bodea, S., Hamer, H. A., Marks, J. A., Haiser, H. J., Turnbaugh, P. J., et al. (2015). Characterization and detection of a widely distributed gene cluster that predicts anaerobic choline utilization by human gut bacteria. MBio 6:e00042–15. doi: 10.1128/mBio.00042-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Mckinnon, L. R., Achilles, S. L., Bradshaw, C. S., Burgener, A., Crucitti, T., Fredricks, D. N., et al. (2019). The evolving facets of bacterial vaginosis: implications for HIV transmission. AIDS Res. Hum. Retroviruses. 35, 219–228. doi: 10.1089/aid.2018.0304

PubMed Abstract | CrossRef Full Text | Google Scholar

Muzny, C. A., Sunesara, I. R., Griswold, M. E., Kumar, R., Lefkowitz, E. J., Mena, L. A., et al. (2014). Association between BVAB1 and high Nugent scores among women with bacterial vaginosis. Diagn. Microbiol. Infect. Dis. 80, 321–323. doi: 10.1016/j.diagmicrobio.2014.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelson, T. M., Borgogna, J. C., Michalek, R. D., Roberts, D. W., Rath, J. M., Glover, E. D., et al. (2018). Cigarette smoking is associated with an altered vaginal tract metabolomic profile. Sci. Rep. 8:852. doi: 10.1038/s41598-017-14943-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ness, R. B., Hillier, S. L., Kip, K. E., Soper, D. E., Stamm, C. A., Mcgregor, J. A., et al. (2004). Bacterial vaginosis and risk of pelvic inflammatory disease. Obstet. Gynecol. 104, 761–769. doi: 10.1097/01.AOG.0000139512.37582.17

PubMed Abstract | CrossRef Full Text | Google Scholar

Ness, R. B., Kip, K. E., Soper, D. E., Hillier, S., Stamm, C. A., Sweet, R. L., et al. (2005). Bacterial vaginosis (BV) and the risk of incident gonococcal or chlamydial genital infection in a predominantly black population. Sex. Transm. Dis. 32, 413–417. doi: 10.1097/01.olq.0000154493.87451.8d

PubMed Abstract | CrossRef Full Text | Google Scholar

Nugent, R. P., Krohn, M. A., and Hillier, S. L. (1991). Reliability of diagnosing bacterial vaginosis is improved by a standardized method of gram stain interpretation. J. Clin. Microbiol. 29, 297–301. doi: 10.1128/JCM.29.2.297-301.1991

PubMed Abstract | CrossRef Full Text | Google Scholar

Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114

PubMed Abstract | CrossRef Full Text | Google Scholar

Petrova, M. I., van den Broek, M., Balzarini, J., Vanderleyden, J., and Lebeer, S. (2013). Vaginal microbiota and its role in HIV transmission and infection. FEMS Microbiol. Rev. 37, 762–792. doi: 10.1111/1574-6976.12029

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2009). FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26, 1641–1650. doi: 10.1093/molbev/msp077

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS ONE. 5:e9490. doi: 10.1371/journal.pone.0009490

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinlan, A. R. (2014). BEDTools: the Swiss-Army tool for genome feature analysis. Curr. Protoc. Bioinformatics 47, 11–34. doi: 10.1002/0471250953.bi1112s47

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinlan, A. R., and Hall, I. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. doi: 10.1093/bioinformatics/btq033

PubMed Abstract | CrossRef Full Text | Google Scholar

Ravel, J., Brotman, R. M., Gajer, P., Ma, B., Nandy, M., Fadrosh, D. W., et al. (2013). Daily temporal dynamics of vaginal microbiota before, during and after episodes of bacterial vaginosis. Microbiome. 1:29. doi: 10.1186/2049-2618-1-29

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwebke, J. R., and Desmond, R. A. (2007). A randomized trial of the duration of therapy with metronidazole plus or minus azithromycin for treatment of symptomatic bacterial vaginosis. Clin. Infect. Dis. 44, 213–219. doi: 10.1086/509577

PubMed Abstract | CrossRef Full Text | Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | CrossRef Full Text | Google Scholar

Srinivasan, S., Morgan, M. T., Fiedler, T. L., Djukovic, D., Hoffman, N. G., Raftery, D., et al. (2015). Metabolic signatures of bacterial vaginosis. MBio 6:e00204–15. doi: 10.1128/mBio.00204-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Srinivasan, S., Morgan, M. T., Liu, C., Matsen, F. A., Hoffman, N. G., Fiedler, T. L., et al. (2013). More than meets the eye: associations of vaginal bacteria with gram stain morphotypes using molecular phylogenetic analysis. PLoS ONE 8:e78633. doi: 10.1371/journal.pone.0078633

PubMed Abstract | CrossRef Full Text | Google Scholar

Taha, T. E., Hoover, D. R., Dallabetta, G. A., Kumwenda, N. I., Mtimavalye, L. A., Yang, L. P., et al. (1998). Bacterial vaginosis and disturbances of vaginal flora: association with increased acquisition of HIV. AIDS 12, 1699–1706. doi: 10.1097/00002030-199813000-00019

PubMed Abstract | CrossRef Full Text | Google Scholar

Weghoff, M. C., Bertsch, J., and Muller, V. (2015). A novel mode of lactate metabolism in strictly anaerobic bacteria. Environ. Microbiol. 17, 670–677. doi: 10.1111/1462-2920.12493

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickham, H. (2009). Ggplot2 : Elegant Graphics for Data Analysis. New York, NY: Springer.

Google Scholar

Witkin, S. S., Mendes-Soares, H., Linhares, I. M., Jayaram, A., Ledger, W. J., and Forney, L. J. (2013). Influence of vaginal bacteria and D- and L-lactic acid isomers on vaginal extracellular matrix metalloproteinase inducer: implications for protection against upper genital tract infections. mBio 4:e00460–13. doi: 10.1128/mBio.00460-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolrath, H., Borén, H., Hallén, A., and Forsum, U. (2002). Trimethylamine content in vaginal secretion and its relation to bacterial vaginosis. APMIS 110, 819–824. doi: 10.1034/j.1600-0463.2002.1101108.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, S. H., Ha, S. M., Kwon, S., Lim, J., Kim, Y., Seo, H., et al. (2017). Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 67, 1613–1617. doi: 10.1099/ijsem.0.001755

PubMed Abstract | CrossRef Full Text | Google Scholar

Yutin, N., and Galperin, M. Y. (2013). A genomic update on clostridial phylogeny: Gram-negative spore formers and other misplaced clostridia. Environ. Microbiol. 15, 2631–2641. doi: 10.1111/1462-2920.12173

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Liang, Y., Lynch, K. H., Dennis, J. J., and Wishart, D. S. (2011). PHAST: a fast phage search tool. Nucleic Acids Res. 39, W347–W352. doi: 10.1093/nar/gkr485

PubMed Abstract | CrossRef Full Text | Google Scholar