Host and body site-specific adaptation of Lactobacillus crispatus genomes

Vaginal microbiomes show ethnic evolutionary dynamics and positive selection of Lactobacillus adhesins driven by a long-term niche-specific process

The vaginal microbiome's composition varies among ethnicities. However, the evolutionary landscape of the vaginal microbiome in the multi-ethnic context remains understudied. We perform a systematic evolutionary analysis of 351 vaginal microbiome samples from 35 multi-ethnic pregnant women, in addition to two validation cohorts, totaling 462 samples from 90 women. Microbiome alpha diversity and community state dynamics show strong ethnic signatures. Lactobacillaceae have a higher ratio of non-synonymous to synonymous polymorphism and lower nucleotide diversity than non-Lactobacillaceae in all ethnicities, with a large repertoire of positively selected genes, including the mucin-binding and cell wall anchor genes. These evolutionary dynamics are driven by the long-term evolutionary process unique to the human vaginal niche. Finally, we propose an evolutionary model reflecting the environmental niches of microbes. Our study reveals the extensive ethnic signatures in vaginal microbial ecology and evolution, highlighting the importance of studying the host-microbiome ecosystem from an evolutionary perspective.

Time series of chicken stool metagenomics and egg metabolomics in changing production systems: preliminary insights from a proof-of-concept

BACKGROUND

Different production systems of livestock animals influence various factors, including the gut microbiota.

METHODS

We investigated whether changing the conditions from barns to free-range chicken farming impacts the microbiome over the course of three weeks. We compared the stool microbiota of chicken from industrial barns after introducing them either in community or separately to a free-range environment.

RESULTS

Over the six time points, 12 taxa-mostly lactobacilli-changed significantly. As expected, the former barn chicken cohort carries more resistances to common antibiotics. These, however, remained positive over the observed period. At the end of the study, we collected eggs and compared metabolomic profiles of the egg white and yolk to profiles of eggs from commercial suppliers. Here, we observed significant differences between commercial and fresh collected eggs as well as differences between the former barn chicken and free-range chicken.

CONCLUSION

Our data indicate that the gut microbiota can undergo alterations over time in response to changes in production systems. These changes subsequently exert an influence on the metabolites found in the eggs. The preliminary results of our proof-of-concept study motivate larger scale observations with more individual chicken and longer observation periods.

Sunscreens can preserve human skin microbiome upon erythemal UV exposure

OBJECTIVE

Ultraviolet radiation (UVR) is a known environmental key factor for premature skin ageing. Only few scientific evidence is available to support the effects of UVR on the skin microbiome. This in vivo pilot study aimed to evaluate the impact on the skin microbiome upon erythemal UV exposure and the protection of UV-exposed skin microbiome by UV filters.

METHODS

Ten female volunteers were treated with an sun protection factor (SPF) 20 sunscreen and placebo formulation (without UV filters) on their upper middle backs and irradiated with an erythemal dose (2 MED) by a solar simulator. Skin swabbing samples from four zones (i.e., unexposed, exposed, sunscreen- and placebo-treated on exposed skin) were collected for the microbiome analysis before and 2 h after UV exposure, respectively, and processed via shallow 16S rRNA Amplicon and Shotgun metagenomic sequencing. An in vitro UV method was developed to confirm the protection of isolated bacterial strains by single UV filters and combinations.

RESULTS

Alpha diversity was impacted by significant inter-individual differences and by treatment rather than by irradiation. Cutibacterium acnes was found to be the most abundant and a confounding factor for diversity. On a species level, Lactobacillus crispatus was negatively associated with UVR and placebo treatment, whereas there was a positive association with sunscreen treatment. The sunscreen treatment also favoured an interaction network with central Micrococcus genus. The in vitro results showed that both single UV filters and combinations had specific effects on the survival rates of L. crispatus, C. acnes, and Staphylococcus epidermidis.

CONCLUSION

We identified potential microorganisms and bacterial interactions that were associated with an SPF 20 sunscreen treatment. The specific protection of L. crispatus as a key player in the UV-exposed skin microbiome and reduction of C. acnes population by UV filters might lead to new cosmetic concepts for photoprotection.

Sequencing and culture-based characterization of the vaginal and uterine microbiota in beef cattle that became pregnant or remained open following artificial insemination

IMPORTANCE

Emerging evidence suggests that microbiome-targeted approaches may provide a novel opportunity to reduce the incidence of reproductive failures in cattle. To develop such microbiome-based strategies, one of the first logical steps is to identify reproductive microbiome features related to fertility and to isolate the fertility-associated microbial species for developing a future bacterial consortium that could be administered before breeding to enhance pregnancy outcomes. Here, we characterized the vaginal and uterine microbiota in beef cattle that became pregnant or remained open via artificial insemination and identified microbiota features associated with fertility. We compared similarities between vaginal and uterine microbiota and between heifers and cows. Using culturing, we provided new insights into the culturable fraction of the vaginal and uterine microbiota and their antimicrobial resistance. Overall, our findings will serve as an important basis for future research aimed at harnessing the vaginal and uterine microbiome for improved cattle fertility.

The mobilome of Lactobacillus crispatus M247 includes two novel genetic elements: Tn7088 coding for a putative bacteriocin and the siphovirus prophage ΦM247

Lactobacillus crispatus is a member of the vaginal and gastrointestinal human microbiota. Here we determined the complete genome sequence of the probiotic strain M247 combining Nanopore and Illumina technologies. The M247 genome is organized in one circular chromosome of 2 336 109 bp, with a GC content of 37.04 % and 2303 ORFs, of which 1962 could be annotated. Analysis of the M247 mobilome, which accounts for 14 % of the whole genome, revealed the presence of: (i) Tn7088, a novel 14 105 bp long integrative and mobilizable element (IME) containing 16 ORFs; (ii) ΦM247, a novel 42 510 bp long siphovirus prophage containing 52 ORFs; (iii) three clustered regularly interspaced short palindromic repeats (CRISPRs); and (iv) 226 insertion sequences (ISs) belonging to 14 different families. Tn7088 has a modular organization including a mobilization module encoding FtsK homologous proteins and a relaxase, an integration/excision module coding for an integrase and an excisionase, and an adaptation module coding for a class I bacteriocin and homologous to the listeriolysin S (lls) locus of Listeria monocytogenes. Genome-wide homology search analysis showed the presence of Tn7088-like elements in 12 out of 23 L. crispatus complete public genomes. Mobilization and integration/excision modules are essentially conserved, while the adaptation module is variable since it is the target site for the integration of different ISs. Prophage ΦM247 contains genes for phage structural proteins, DNA replication and packaging, lysogenic and lytic cycles. ΦM247-like prophages are present in seven L. crispatus complete genomes, with sequence variability mainly due to the integration of ISs. PCR and sequencing showed that the Tn7088 IME excises from the M247 chromosome producing a circular form at a concentration of 4.32×10-5 copies per chromosome, and reconstitution of the Tn7088 chromosomal target site occurred at 6.65×10-4 copies per chromosome. The ΦM247 prophage produces an excised form and a reconstituted target site at a level of 3.90×10-5 and 2.48×10-5 copies per chromosome, respectively. This study identified two novel genetic elements in L. crispatus. Tn7088 represents the first example of an IME carrying a biosynthetic gene cluster for a class I bacteriocin in L. crispatus.

Whole-body microbiota of newborn calves and their response to prenatal vitamin and mineral supplementation

Early life microbial colonization and factors affecting colonization patterns are gaining interest due to recent developments suggesting that early life microbiome may play a role in Developmental Origins of Health and Disease. In cattle, limited information exists on the early microbial colonization of anatomical sites involved in bovine health beyond the gastrointestinal tract. Here, we investigated 1) the initial microbial colonization of seven different anatomical locations in newborn calves and 2) whether these early life microbial communities and 3) serum cytokine profiles are influenced by prenatal vitamin and mineral (VTM) supplementation. Samples were collected from the hoof, liver, lung, nasal cavity, eye, rumen (tissue and fluid), and vagina of beef calves that were born from dams that either received or did not receive VTM supplementation throughout gestation (n = 7/group). Calves were separated from dams immediately after birth and fed commercial colostrum and milk replacer until euthanasia at 30 h post-initial colostrum feeding. The microbiota of all samples was assessed using 16S rRNA gene sequencing and qPCR. Calf serum was subjected to multiplex quantification of 15 bovine cytokines and chemokines. Our results indicated that the hoof, eye, liver, lung, nasal cavity, and vagina of newborn calves were colonized by site-specific microbiota, whose community structure differed from the ruminal-associated communities (0.64 ≥ R² ≥ 0.12, p ≤ 0.003). The ruminal fluid microbial community was the only one that differed by treatment (p < 0.01). However, differences (p < 0.05) by treatment were detected in microbial richness (vagina); diversity (ruminal tissue, fluid, and eye); composition at the phylum and genus level (ruminal tissue, fluid, and vagina); and in total bacterial abundance (eye and vagina). From serum cytokines evaluated, concentration of chemokine IP-10 was greater (p = 0.02) in VTM calves compared to control calves. Overall, our results suggest that upon birth, the whole-body of newborn calves are colonized by relatively rich, diverse, and site-specific bacterial communities. Noticeable differences were observed in ruminal, vaginal, and ocular microbiota of newborn calves in response to prenatal VTM supplementation. These findings can derive future hypotheses regarding the initial microbial colonization of different body sites, and on maternal micronutrient consumption as a factor that may influence early life microbial colonization.

Antibacterial potential of commercial and wild lactic acid bacteria strains isolated from ovine and caprine raw milk against Mycoplasma agalactiae

Introduction

The complexity of fighting contagious agalactia (CA) has raised the necessity of alternative antimicrobial therapies, such as probiotics. Lactic acid bacteria (LAB) are present in the mammary gland of small ruminants and their antimicrobial effect have been previously described against species like Mycoplasma bovis but never against Mycoplasma agalactiae (Ma). This in vitro study aims to evaluate the antimicrobial activity against Ma of ovine and caprine LAB strains and a human commercial probiotic (L2) of Lactobacillus spp.

Methods

A total of 63 possible LAB strains were isolated from nine ovine and caprine farms in Spain, three isolates (33B, 248D, and 120B) from the 63 strains were selected, based on their capacity to grow in a specific medium in vitro, for an in vitro experiment to assess their antimicrobial activity against Ma in Ultra High Temperature (UHT) processed goat milk (GM). A women commercial vaginal probiotic was also included in the study. The inoculum of L2 was prepared at a concentration of 3.24 × 10⁸  CFU/mL and the average concentration of the inoculum of the wild LAB varied from 7.9 × 10⁷ to 8.4 × 10⁸  CFU/mL.

Results

The commercial probiotic L2 significantly reduced the concentration of Ma to 0.000 log CFU/mL (p < 0.001), strain 33B reduced it from 7.185 to 1.279 log CFU/mL (p < 0.001), and 120B from 6.825 to 6.466 log CFU/mL (p < 0.05). Strain 248D presented a bacteriostatic effect in GM. Moreover, the three wild strains and the commercial probiotic produced a significative reduction of the pH (p < 0.001).

Discussion

This is the first in vivo report of the antimicrobial potential of LAB strains against Ma and its interaction. Our results support possible future alternative strategies to antibiotic therapy, previously not contemplated, to fight CA in small ruminants. Further studies are necessary to elucidate the action mechanisms through which these LAB are able to inhibit Ma and to assess the safety of using these strains in possible in vivo studies.

Female reproductive tract-organ axes

The female reproductive tract (FRT) and remote/versatile organs in the body share bidirectional communication. In this review, we discuss the framework of the "FRT-organ axes." Each axis, namely, the vagina-gut axis, uterus-gut axis, ovary-gut axis, vagina-bladder axis, vagina-oral axis, uterus-oral axis, vagina-brain axis, uterus-brain axis, and vagina-joint axis, is comprehensively discussed separately. Each axis could be involved in the pathogenesis of not only gynecological diseases but also diseases occurring apart from the FRT. Although the microbiota is clearly a key player in the FRT-organ axes, more quantitative insight into the homeostasis of the microbiota could be provided by host function measurements rather than current microbe-centric approaches. Therefore, investigation of the FRT-organ axes would provide us with a multicentric approach, including immune, neural, endocrine, and metabolic aspects, for understanding the homeostatic mechanism of women's bodies. The framework of the FRT-organ axes could also provide insights into finding new therapeutic approaches to maintain women's health.

The Differing Roles of Flavins and Quinones in Extracellular Electron Transfer in Lactiplantibacillus plantarum

Lactiplantibacillus plantarum is a lactic acid bacterium that is commonly found in the human gut and fermented food products. Despite its overwhelmingly fermentative metabolism, this microbe can perform extracellular electron transfer (EET) when provided with an exogenous quinone, 1,4-dihydroxy-2-naphthoic acid (DHNA), and riboflavin. However, the separate roles of DHNA and riboflavin in EET in L. plantarum have remained unclear. Here, we seek to understand the role of quinones and flavins in EET by monitoring iron and anode reduction in the presence and absence of these small molecules. We found that addition of either DHNA or riboflavin can support robust iron reduction, indicating electron transfer to extracellular iron occurs through both flavin-dependent and DHNA-dependent routes. Using genetic mutants of L. plantarum, we found that flavin-dependent iron reduction requires Ndh2 and EetA, while DHNA-dependent iron reduction largely relies on Ndh2 and PplA. In contrast to iron reduction, DHNA-containing medium supported more robust anode reduction than riboflavin-containing medium, suggesting electron transfer to an anode proceeds most efficiently through the DHNA-dependent pathway. Furthermore, we found that flavin-dependent anode reduction requires EetA, Ndh2, and PplA, while DHNA-dependent anode reduction requires Ndh2 and PplA. Taken together, we identify multiple EET routes utilized by L. plantarum and show that the EET route depends on access to environmental biomolecules and on the electron acceptor. This work expands our molecular-level understanding of EET in Gram-positive microbes and provides additional opportunities to manipulate EET for biotechnology. IMPORTANCE Lactic acid bacteria are named because of their nearly exclusive fermentative metabolism. Thus, the recent observation of EET activity-typically associated with anaerobic respiration-in this class of organisms has forced researchers to rethink the rules governing microbial metabolic strategies. Our identification of multiple routes for EET in L. plantarum that depend on two different redox active small molecules expands our understanding of how microbes metabolically adapt to different environments to gain an energetic edge and how these processes can be manipulated for biotechnological uses. Understanding the role of EET in lactic acid bacteria is of great importance due to the significance of lactic acid bacteria in agriculture, bioremediation, food production, and gut health. Furthermore, the maintenance of multiple EET routes speaks to the importance of this process to function under a variety of environmental conditions.

Identification of shared bacterial strains in the vaginal microbiota of related and unrelated reproductive-age mothers and daughters using genome-resolved metagenomics

It has been suggested that the human microbiome might be vertically transmitted from mother to offspring and that early colonizers may play a critical role in development of the immune system. Studies have shown limited support for the vertical transmission of the intestinal microbiota but the derivation of the vaginal microbiota remains largely unknown. Although the vaginal microbiota of children and reproductive age women differ in composition, the vaginal microbiota could be vertically transmitted. To determine whether there was any support for this hypothesis, we examined the vaginal microbiota of daughter-mother pairs from the Baltimore metropolitan area (ages 14-27, 32-51; n = 39). We assessed whether the daughter's microbiota was similar in composition to their mother's using metataxonomics. Permutation tests revealed that while some pairs did have similar vaginal microbiota, the degree of similarity did not exceed that expected by chance. Genome-resolved metagenomics was used to identify shared bacterial strains in a subset of the families (n = 22). We found a small number of bacterial strains that were shared between mother-daughter pairs but identified more shared strains between individuals from different families, indicating that vaginal bacteria may display biogeographic patterns. Earlier-in-life studies are needed to demonstrate vertical transmission of the vaginal microbiota.

Viral biogeography of the mammalian gut and parenchymal organs

The mammalian virome has been linked to health and disease but our understanding of how it is structured along the longitudinal axis of the mammalian gastrointestinal tract (GIT) and other organs is limited. Here, we report a metagenomic analysis of the prokaryotic and eukaryotic virome occupying luminal and mucosa-associated habitats along the GIT, as well as parenchymal organs (liver, lung and spleen), in two representative mammalian species, the domestic pig and rhesus macaque (six animals per species). Luminal samples from the large intestine of both mammals harboured the highest loads and diversity of bacteriophages (class Caudoviricetes, family Microviridae and others). Mucosal samples contained much lower viral loads but a higher proportion of eukaryotic viruses (families Astroviridae, Caliciviridae, Parvoviridae). Parenchymal organs contained bacteriophages of gut origin, in addition to some eukaryotic viruses. Overall, GIT virome composition was specific to anatomical region and host species. Upper GIT and mucosa-specific viruses were greatly under-represented in distal colon samples (a proxy for faeces). Nonetheless, certain viral and phage species were ubiquitous in all samples from the oral cavity to the distal colon. The dataset and its accompanying methodology may provide an important resource for future work investigating the biogeography of the mammalian gut virome.

Genetic Elements Orchestrating Lactobacillus crispatus Glycogen Metabolism in the Vagina

Glycogen in the female lower reproductive tract is a major carbon source for colonization and acidification by common vaginal Lactobacillus species, such as Lactobacillus crispatus. Previously, we identified the amylopullulanase encoding gene pulA of Lactobacillus crispatus to correlate with the ability to autonomously utilize glycogen for growth. Here, we further characterize genetic variation and differential regulation of pulA affecting the presence of its gene product on the outer surface layer. We show that alpha-glucan degrading activity dissipates when Lactobacillus crispatus is grown on glucose, maltose and maltotriose, in agreement with carbon catabolite repression elements flanking the pulA gene. Proteome analysis of the S-layer confirmed that the amylopullulanase protein is highly abundant in an S-layer enriched fraction, but not in a strain with a defective amylopullulanase variant or in an amylopullulanase-sufficient strain grown on glucose. In addition, we provide evidence that Lactobacillus crispatus pulA mutants are relevant in vivo, as they are commonly observed in metagenome datasets of human vaginal microbial communities. Analysis of the largest publicly available dataset of 1507 human vaginal metagenomes indicates that among the 270 samples that contain a Lactobacillus crispatuspulA gene, 62 samples (23%) had a defective variant of this gene. Taken together, these results demonstrate that both environmental, as well as genetic factors explain the variation of Lactobacillus crispatus alpha-glucosidases in the vaginal environment.

Lactobacillus, glycans and drivers of health in the vaginal microbiome

A microbiome consists of microbes and their genomes, encompassing bacteria, viruses, fungi, protozoa, archaea, and eukaryotes. These elements interact dynamically in the specific environment in which they reside and evolve. In the past decade, studies of various microbiomes have been prevalent in the scientific literature, accounting for the shift from culture-dependent to culture-independent identification of microbes using new high-throughput sequencing technologies that decipher their composition and sometimes provide insights into their functions. Despite tremendous advances in understanding the gut microbiome, relatively little attention has been devoted to the vaginal environment, notably regarding the ubiquity and diversity of glycans which denote the significant role they play in the maintenance of homeostasis. Hopefully, emerging technologies will aid in the determination of what is a healthy vaginal microbiome, and provide insights into the roles of Lactobacillus, glycans and microbiome-related drivers of health and disease.

Comparative Genomic Analyses and CRISPR-Cas Characterization of Cutibacterium acnes Provide Insights Into Genetic Diversity and Typing Applications

Cutibacterium acnes is an important member of the human skin microbiome and plays a critical role in skin health and disease. C. acnes encompasses different phylotypes that have been found to be associated with different skin phenotypes, suggesting a genetic basis for their impact on skin health. Here, we present a comprehensive comparative analysis of 255 C. acnes genomes to provide insights into the species genetic diversity and identify unique features that define various phylotypes. Results revealed a relatively small and open pan genome (6,240 genes) with a large core genome (1,194 genes), and three distinct phylogenetic clades, with multiple robust sub-clades. Furthermore, we identified several unique gene families driving differences between distinct C. acnes clades. Carbohydrate transporters, stress response mechanisms and potential virulence factors, potentially involved in competitive growth and host colonization, were detected in type I strains, which are presumably responsible for acne. Diverse type I-E CRISPR-Cas systems and prophage sequences were detected in select clades, providing insights into strain divergence and adaptive differentiation. Collectively, these results enable to elucidate the fundamental differences among C. acnes phylotypes, characterize genetic elements that potentially contribute to type I-associated dominance and disease, and other key factors that drive the differentiation among clades and sub-clades. These results enable the use of comparative genomics analyses as a robust method to differentiate among the C. acnes genotypes present in the skin microbiome, opening new avenues for the development of biotherapeutics to manipulate the skin microbiota.

Insight into phenotypic and genotypic differences between vaginal Lactobacillus crispatus BC5 and Lactobacillus gasseri BC12 to unravel nutritional and stress factors influencing their metabolic activity

The vaginal microbiota, normally characterized by lactobacilli presence, is crucial for vaginal health. Members belonging to L. crispatus and L. gasseri species exert crucial protective functions against pathogens, although a total comprehension of factors that influence their dominance in healthy women is still lacking. Here we investigated the complete genome sequence and comprehensive phenotypic profile of L. crispatus strain BC5 and L. gasseri strain BC12, two vaginal strains featured by anti-bacterial and anti-viral activities. Phenotype microarray (PM) results revealed an improved capacity of BC5 to utilize different carbon sources as compared to BC12, although some specific carbon sources that can be associated to the human diet were only metabolized by BC12, i.e. uridine, amygdalin, tagatose. Additionally, the two strains were mostly distinct in the capacity to utilize the nitrogen sources under analysis. On the other hand, BC12 showed tolerance/resistance towards twice the number of stressors (i.e. antibiotics, toxic metals etc.) with respect to BC5. The divergent phenotypes observed in PM were supported by the identification in either BC5 or BC12 of specific genetic determinants that were found to be part of the core genome of each species. The PM results in combination with comparative genome data provide insights into the possible environmental factors and genetic traits supporting the predominance of either L. crispatus BC5 or L. gasseri BC12 in the vaginal niche, giving also indications for metabolic predictions at the species level.

Characterization of the first vaginal Lactobacillus crispatus genomes isolated in Brazil

BACKGROUND

Lactobacillus crispatus is the dominant species in the vaginal microbiota associated with health and considered a homeostasis biomarker. Interestingly, some strains are even used as probiotics. However, the genetic mechanisms of L. crispatus involved in the control of the vaginal microbiome and protection against bacterial vaginosis (BV) are not entirely known. To further investigate these mechanisms, we sequenced and characterized the first four L. crispatus genomes from vaginal samples from Brazilian women and used genome-wide association study (GWAS) and comparative analyses to identify genetic mechanisms involved in healthy or BV conditions and selective pressures acting in the vaginal microbiome.

METHODS

The four genomes were sequenced, assembled using ten different strategies and automatically annotated. The functional characterization was performed by bioinformatics tools comparing with known probiotic strains. Moreover, it was selected one representative strain (L. crispatus CRI4) for in vitro detection of phages by electron microscopy. Evolutionary analysis, including phylogeny, GWAS and positive selection were performed using 46 public genomes strains representing health and BV conditions.

RESULTS

Genes involved in probiotic effects such as lactic acid production, hydrogen peroxide, bacteriocins, and adhesin were identified. Three hemolysins and putrescine production were predicted, although these features are also present in other probiotic strains. The four genomes presented no plasmids, but 14 known families insertion sequences and several prophages were detected. However, none of the mobile genetic elements contained antimicrobial resistance genes. The genomes harbor a CRISPR-Cas subtype II-A system that is probably inactivated due to fragmentation of the genes csn2 and cas9. No genomic feature was associated with a health condition, perhaps due to its multifactorial characteristic. Five genes were identified as under positive selection, but the selective pressure remains to be discovered. In conclusion, the Brazilian strains investigated in this study present potential protective properties, although in vitro and in vivo studies are required to confirm their efficacy and safety to be considered for human use.

Bactopia: a Flexible Pipeline for Complete Analysis of Bacterial Genomes

Sequencing of bacterial genomes using Illumina technology has become such a standard procedure that often data are generated faster than can be conveniently analyzed. We created a new series of pipelines called Bactopia, built using Nextflow workflow software, to provide efficient comparative genomic analyses for bacterial species or genera. Bactopia consists of a data set setup step (Bactopia Data Sets [BaDs]), which creates a series of customizable data sets for the species of interest, the Bactopia Analysis Pipeline (BaAP), which performs quality control, genome assembly, and several other functions based on the available data sets and outputs the processed data to a structured directory format, and a series of Bactopia Tools (BaTs) that perform specific postprocessing on some or all of the processed data. BaTs include pan-genome analysis, computing average nucleotide identity between samples, extracting and profiling the 16S genes, and taxonomic classification using highly conserved genes. It is expected that the number of BaTs will increase to fill specific applications in the future. As a demonstration, we performed an analysis of 1,664 public Lactobacillus genomes, focusing on Lactobacillus crispatus, a species that is a common part of the human vaginal microbiome. Bactopia is an open source system that can scale from projects as small as one bacterial genome to ones including thousands of genomes and that allows for great flexibility in choosing comparison data sets and options for downstream analysis. Bactopia code can be accessed at https://www.github.com/bactopia/bactopiaIMPORTANCE It is now relatively easy to obtain a high-quality draft genome sequence of a bacterium, but bioinformatic analysis requires organization and optimization of multiple open source software tools. We present Bactopia, a pipeline for bacterial genome analysis, as an option for processing bacterial genome data. Bactopia also automates downloading of data from multiple public sources and species-specific customization. Because the pipeline is written in the Nextflow language, analyses can be scaled from individual genomes on a local computer to thousands of genomes using cloud resources. As a usage example, we processed 1,664 Lactobacillus genomes from public sources and used comparative analysis workflows (Bactopia Tools) to identify and analyze members of the L. crispatus species.