Insertion sequence content reflects genome plasticity in strains of the root nodule actinobacterium Frankia

Genomic characterization of closely related species in the Rumoiensis clade infers ecogenomic signatures to non-marine environments

Members of the family Vibrionaceae are generally found in marine and brackish environments, playing important roles in nutrient cycling. The Rumoiensis clade is an unconventional group in the genus Vibrio, currently comprising six species from different origins including two species isolated from non-marine environments. In this study, we performed comparative genome analysis of all six species in the clade using their complete genome sequences. We found that two non-marine species, Vibrio casei and Vibrio gangliei, lacked the genes responsible for algal polysaccharide degradation, while a number of glycoside hydrolase genes were enriched in these two species. Expansion of insertion sequences was observed in V. casei and Vibrio rumoiensis, which suggests ongoing genomic changes associated with niche adaptations. The genes responsible for the metabolism of glucosylglycerate, a compound known to play a role as compatible solutes under nitrogen limitation, were conserved across the clade. These characteristics, along with genes encoding species-specific functions, may reflect the habit expansion which has led to the current distribution of Rumoiensis clade species. Genome analysis of all species in a single clade give us valuable insights into the genomic background of the Rumoiensis clade species and emphasize the genomic diversity and versatility of Vibrionaceae.

Large genomic differences between Moraxella bovoculi isolates acquired from the eyes of cattle with infectious bovine keratoconjunctivitis versus the deep nasopharynx of asymptomatic cattle

Moraxella bovoculi is a recently described bacterium that is associated with infectious bovine keratoconjunctivitis (IBK) or “pinkeye” in cattle. In this study, closed circularized genomes were generated for seven M. bovoculi isolates: three that originated from the eyes of clinical IBK bovine cases and four from the deep nasopharynx of asymptomatic cattle. Isolates that originated from the eyes of IBK cases profoundly differed from those that originated from the nasopharynx of asymptomatic cattle in genome structure, gene content and polymorphism diversity and consequently placed into two distinct phylogenetic groups. These results suggest that there are genetically distinct strains of M. bovoculi that may not associate with IBK.

Dynamics of bacterial insertion sequences: can transposition bursts help the elements persist?

Background

Currently there is no satisfactory explanation for why bacterial insertion sequences (ISs) widely occur across prokaryotes despite being mostly harmful to their host genomes. Rates of horizontal gene transfer are likely to be too low to maintain ISs within a population. IS-induced beneficial mutations may be important for both prevalence of ISs and microbial adaptation to changing environments but may be too rare to sustain IS elements in the long run. Environmental stress can induce elevated rates of IS transposition activities; such episodes are known as ‘transposition bursts’. By examining how selective forces and transposition events interact to influence IS dynamics, this study asks whether transposition bursts can lead to IS persistence.

Results

We show through a simulation model that ISs are gradually eliminated from a population even if IS transpositions occasionally cause advantageous mutations. With beneficial mutations, transposition bursts create variation in IS copy numbers and improve cell fitness on average. However, these benefits are not usually sufficient to overcome the negative selection against the elements, and transposition bursts amplify the mean fitness effect which, if negative, simply accelerates the extinction of ISs. If down regulation of transposition occurs, IS extinctions are reduced while ISs still generate variation amongst bacterial genomes.

Conclusions

Transposition bursts do not help ISs persist in a bacterial population in the long run because most burst-induced mutations are deleterious and therefore not favoured by natural selection. However, bursts do create more genetic variation through which occasional advantageous mutations can help organisms adapt. Regulation of IS transposition bursts and stronger positive selection of the elements interact to slow down the burst-induced extinction of ISs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0560-5) contains supplementary material, which is available to authorized users.

Candidatus Frankia Datiscae Dg1, the Actinobacterial Microsymbiont of Datisca glomerata, Expresses the Canonical nod Genes nodABC in Symbiosis with Its Host Plant

Frankia strains are nitrogen-fixing soil actinobacteria that can form root symbioses with actinorhizal plants. Phylogenetically, symbiotic frankiae can be divided into three clusters, and this division also corresponds to host specificity groups. The strains of cluster II which form symbioses with actinorhizal Rosales and Cucurbitales, thus displaying a broad host range, show suprisingly low genetic diversity and to date can not be cultured. The genome of the first representative of this cluster, Candidatus Frankia datiscae Dg1 (Dg1), a microsymbiont of Datisca glomerata, was recently sequenced. A phylogenetic analysis of 50 different housekeeping genes of Dg1 and three published Frankia genomes showed that cluster II is basal among the symbiotic Frankia clusters. Detailed analysis showed that nodules of D. glomerata, independent of the origin of the inoculum, contain several closely related cluster II Frankia operational taxonomic units. Actinorhizal plants and legumes both belong to the nitrogen-fixing plant clade, and bacterial signaling in both groups involves the common symbiotic pathway also used by arbuscular mycorrhizal fungi. However, so far, no molecules resembling rhizobial Nod factors could be isolated from Frankia cultures. Alone among Frankia genomes available to date, the genome of Dg1 contains the canonical nod genes nodA, nodB and nodC known from rhizobia, and these genes are arranged in two operons which are expressed in D. glomerata nodules. Furthermore, Frankia Dg1 nodC was able to partially complement a Rhizobium leguminosarum A34 nodC::Tn5 mutant. Phylogenetic analysis showed that Dg1 Nod proteins are positioned at the root of both α- and β-rhizobial NodABC proteins. NodA-like acyl transferases were found across the phylum Actinobacteria, but among Proteobacteria only in nodulators. Taken together, our evidence indicates an Actinobacterial origin of rhizobial Nod factors.

Local hopping mobile DNA implicated in pseudogene formation and reductive evolution in an obligate cyanobacteria-plant symbiosis

Background

Insertion sequences (ISs) are approximately 1 kbp long “jumping” genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (“cut-and-paste”) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium ‘Nostoc azollae’ 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of ‘Nostoc azollae’ 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria.

Results

Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (“local hopping”) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of ‘Nostoc azollae’ 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style.

Conclusions

The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wide this likely explains the proximity between IS elements and pseudogenes. These findings suggest that ISs facilitate the reductive evolution in for instance in the symbiotic cyanobacterium ‘Nostoc azollae’ 0708 and in other obligate prokaryotic symbionts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1386-7) contains supplementary material, which is available to authorized users.

Mobile genetic elements in the bacterial phylum Acidobacteria

Analysis of the genome of Candidatus Solibacter usitatus Ellin6076, a member of the phylum Acidobacteria, revealed a large number of genes associated with mobile genetic elements. These genes encoded transposases, insertion sequence elements and phage integrases. When the amino acid sequences of the mobile element-associated genes were compared, many of them had high (90–100%) amino acid sequence identities, suggesting that these genes may have recently duplicated and dispersed throughout the genome. Although phage integrase encoding genes were prevalent in the Can. S. usitatus Ellin6076 genome, no intact prophage regions were found. This suggests that the Can. S. usitatus Ellin6076 large genome arose by horizontal gene transfer via ancient bacteriophage and/or plasmid-mediated transduction, followed by widespread small-scale gene duplications, resulting in an increased number of paralogs encoding traits that could provide selective metabolic, defensive and regulatory advantages in the soil environment. Here we examine the mobile element repertoire of Can. S. usitatus Ellin6076 in comparison to other genomes from the Acidobacteria phylum, reviewing published studies and contributing some new analyses. We also discuss the presence and potential roles of mobile elements in members of this phylum that inhabit a variety of environments.

Characterization of pseudogenes in members of the order Frankineae

Pseudogenes are defined as non-functional relatives of genes whose protein-coding abilities are lost and are no longer expressed within cells. They are an outcome of accumulation of mutations within a gene whose end product is not essential for survival. Proper investigation of the procedure of pseudogenization is relevant for estimating occurrence of duplications in genomes. Frankineae houses an interesting group of microorganisms, carving a niche in the microbial world. This study was undertaken with the objective of determining the abundance of pseudogenes, understanding strength of purifying selection, investigating evidence of pseudogene expression, and analysing their molecular nature, their origin, evolution and deterioration patterns amongst domain families. Investigation revealed the occurrence of 956 core pFAM families sharing common characteristics indicating co-evolution. WD40, Rve_3, DDE_Tnp_IS240 and phage integrase core domains are larger families, having more pseudogenes, signifying a probability of harmful foreign genes being disabled within transposable elements. High selective pressure depicted that gene families rapidly duplicating and evolving undoubtedly facilitated creation of a number of pseudogenes in Frankineae. Codon usage analysis between protein-coding genes and pseudogenes indicated a wide degree of variation with respect to different factors. Moreover, the majority of pseudogenes were under the effect of purifying selection. Frankineae pseudogenes were under stronger selective constraints, indicating that they were functional for a very long time and became pseudogenes abruptly. The origin and deterioration of pseudogenes has been attributed to selection and mutational pressure acting upon sequences for adapting to stressed soil environments.

Different dynamics of genome content shuffling among host-specificity groups of the symbiotic actinobacterium Frankia

BACKGROUND

Frankia is a genus of soil actinobacteria forming nitrogen-fixing root-nodule symbiotic relationships with non-leguminous woody plant species, collectively called actinorhizals, from eight dicotyledonous families. Frankia strains are classified into four host-specificity groups (HSGs), each of which exhibits a distinct host range. Genome sizes of representative strains of Alnus, Casuarina, and Elaeagnus HSGs are highly diverged and are positively correlated with the size of their host ranges.

RESULTS

The content and size of 12 Frankia genomes were investigated by in silico comparative genome hybridization and pulsed-field gel electrophoresis, respectively. Data were collected from four query strains of each HSG and compared with those of reference strains possessing completely sequenced genomes. The degree of difference in genome content between query and reference strains varied depending on HSG. Elaeagnus query strains were missing the greatest number (22-32%) of genes compared with the corresponding reference genome; Casuarina query strains lacked the fewest (0-4%), with Alnus query strains intermediate (14-18%). In spite of the remarkable gene loss, genome sizes of Alnus and Elaeagnus query strains were larger than would be expected based on total length of the absent genes. In contrast, Casuarina query strains had smaller genomes than expected.

CONCLUSIONS

The positive correlation between genome size and host range held true across all investigated strains, supporting the hypothesis that size and genome content differences are responsible for observed diversity in host plants and host plant biogeography among Frankia strains. In addition, our results suggest that different dynamics of shuffling of genome content have contributed to these symbiotic and biogeographic adaptations. Elaeagnus strains, and to a lesser extent Alnus strains, have gained and lost many genes to adapt to a wide range of environments and host plants. Conversely, rather than acquiring new genes, Casuarina strains have discarded genes to reduce genome size, suggesting an evolutionary orientation towards existence as specialist symbionts.

Elucidation of insertion elements carried on plasmids and in vitro construction of shuttle vectors from the toxic cyanobacterium Planktothrix

Several gene clusters that are responsible for toxin synthesis in bloom-forming cyanobacteria have been found to be associated with transposable elements (TEs). In particular, insertion sequence (IS) elements were shown to play a role in the inactivation or recombination of the genes responsible for cyanotoxin synthesis. Plasmids have been considered important vectors of IS element distribution to the host. In this study, we aimed to elucidate the IS elements propagated on the plasmids and the chromosome of the toxic cyanobacterium Planktothrix agardhii NIVA-CYA126/8 by means of high-throughput sequencing. In total, five plasmids (pPA5.5, pPA14, pPA50, pPA79, and pPA115, of 5, 6, 50, 79, and 120 kbp, respectively) were elucidated, and two plasmids (pPA5.5, pPA115) were found to propagate full IS element copies. Large stretches of shared DNA information between plasmids were constituted of TEs. Two plasmids (pPA5.5, pPA14) were used as candidates to engineer shuttle vectors (named pPA5.5SV and pPA14SV, respectively) in vitro by PCR amplification and the subsequent transposition of the Tn5 cat transposon containing the R6Kγ origin of replication of Escherichia coli. While pPA5.5SV was found to be fully segregated, pPA14SV consistently co-occurred with its wild-type plasmid even under the highest selective pressure. Interestingly, the Tn5 cat transposon became transferred by homologous recombination into another plasmid, pPA50. The availability of shuttle vectors is considered to be of relevance in investigating genome plasticity as a consequence of homologous recombination events. Combining the potential of high-throughput sequencing and in vitro production of shuttle vectors makes it simple to produce species-specific shuttle vectors for many cultivable prokaryotes.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

The Nocardia cyriacigeorgica GUH-2 genome shows ongoing adaptation of an environmental Actinobacteria to a pathogen's lifestyle

Background

Nocardia cyriacigeorgica is recognized as one of the most prevalent etiological agents of human nocardiosis. Human exposure to these Actinobacteria stems from direct contact with contaminated environmental matrices. The full genome sequence of N. cyriacigeorgica strain GUH-2 was studied to infer major trends in its evolution, including the acquisition of novel genetic elements that could explain its ability to thrive in multiple habitats.

Results

N. cyriacigeorgica strain GUH-2 genome size is 6.19 Mb-long, 82.7% of its CDS have homologs in at least another actinobacterial genome, and 74.5% of these are found in N. farcinica. Among N. cyriacigeorgica specific CDS, some are likely implicated in niche specialization such as those involved in denitrification and RuBisCO production, and are found in regions of genomic plasticity (RGP). Overall, 22 RGP were identified in this genome, representing 11.4% of its content. Some of these RGP encode a recombinase and IS elements which are indicative of genomic instability. CDS playing part in virulence were identified in this genome such as those involved in mammalian cell entry or encoding a superoxide dismutase. CDS encoding non ribosomal peptide synthetases (NRPS) and polyketide synthases (PKS) were identified, with some being likely involved in the synthesis of siderophores and toxins. COG analyses showed this genome to have an organization similar to environmental Actinobacteria.

Conclusion

N. cyriacigeorgica GUH-2 genome shows features suggesting a diversification from an ancestral saprophytic state. GUH-2 ability at acquiring foreign DNA was found significant and to have led to functional changes likely beneficial for its environmental cycle and opportunistic colonization of a human host.

A novel human-infection-derived bacterium provides insights into the evolutionary origins of mutualistic insect-bacterial symbioses

Despite extensive study, little is known about the origins of the mutualistic bacterial endosymbionts that inhabit approximately 10% of the world's insects. In this study, we characterized a novel opportunistic human pathogen, designated “strain HS,” and found that it is a close relative of the insect endosymbiont Sodalis glossinidius. Our results indicate that ancestral relatives of strain HS have served as progenitors for the independent descent of Sodalis-allied endosymbionts found in several insect hosts. Comparative analyses indicate that the gene inventories of the insect endosymbionts were independently derived from a common ancestral template through a combination of irreversible degenerative changes. Our results provide compelling support for the notion that mutualists evolve from pathogenic progenitors. They also elucidate the role of degenerative evolutionary processes in shaping the gene inventories of symbiotic bacteria at a very early stage in these mutualistic associations.

Many insects harbor symbiotic bacteria that perform diverse functions within their hosts. However, the origins of these associations have been difficult to define. In this study we isolate a novel bacterium from a human infection and show that this bacterium is a close relative of the Sodalis-allied clade of insect symbionts. Comparative genomic analyses reveal that this organism maintains many genes that have been inactivated and lost independently in derived insect symbionts as a result of rapid genome degeneration. Our work also shows that recently derived Sodalis-allied symbionts maintain a significant population of “cryptic” pseudogenes that are assumed to have no beneficial function in the symbiosis but have not yet accumulated mutations that disrupt their translation. Taken together, our results show that genome degeneration proceeds rapidly following the onset of symbiosis. They also highlight the potential for diverse insect taxa to acquire closely related insect symbionts as a consequence of vectoring bacterial pathogens to plants and animals.

The biology of Frankia sp. strains in the post-genome era

Progress in understanding symbiotic determinants involved in the N(2)-fixing actinorhizal plant symbioses has been slow but steady. Problems persist with studying the bacterial contributions to the symbiosis using traditional microbiological techniques. However, recent years have seen the emergence of several genomes from Frankia sp. strains and the development of techniques for manipulating plant gene expression. Approaches to understanding the bacterial side of the symbiosis have employed a range of techniques that reveal the proteomes and transcriptomes from both cultured and symbiotic frankiae. The picture beginning to emerge provides some perspective on the heterogeneity of frankial populations in both conditions. In general, frankial populations in root nodules seem to maintain a rather robust metabolism that includes nitrogen fixation and substantial biosynthesis and energy-generating pathways, along with a modified ammonium assimilation program. To date, particular bacterial genes have not been implicated in root nodule formation but some hypotheses are emerging with regard to how the plant and microorganism manage to coexist. In particular, frankiae seem to present a nonpathogenic presence to the plant that may have the effect of minimizing some plant defense responses. Future studies using high-throughput approaches will likely clarify the range of bacterial responses to symbiosis that will need to be understood in light of the more rapidly advancing work on the plant host.

The distribution of a phage-related insertion sequence element in the cyanobacterium, Microcystis aeruginosa

The cyanophage Ma-LMM01, specifically-infecting Microcystis aeruginosa, has an insertion sequence (IS) element that we named IS607-cp showing high nucleotide similarity to a counterpart in the genome of the cyanobacterium Cyanothece sp. We tested 21 strains of M. aeruginosa for the presence of IS607-cp using PCR and detected the element in strains NIES90, NIES112, NIES604, and RM6. Thermal asymmetric interlaced PCR (TAIL-PCR) revealed each of these strains has multiple copies of IS607-cp. Some of the ISs were classified into three types based on their inserted positions; IS607-cp-1 is common in strains NIES90, NIES112 and NIES604, whereas IS607-cp-2 and IS607-cp-3 are specific to strains NIES90 and RM6, respectively. This multiplicity may reflect the replicative transposition of IS607-cp. The sequence of IS607-cp in Ma-LMM01 showed robust affinity to those found in M. aeruginosa and Cyanothece spp. in a phylogenetic tree inferred from counterparts of various bacteria. This suggests the transfer of IS607-cp between the cyanobacterium and its cyanophage. We discuss the potential role of Ma-LMM01-related phages as donors of IS elements that may mediate the transfer of IS607-cp; and thereby partially contribute to the genome plasticity of M. aeruginosa.

Transcriptomes of Frankia sp. strain CcI3 in growth transitions

Background

Frankia sp. strains are actinobacteria that form N₂-fixing root nodules on angiosperms. Several reference genome sequences are available enabling transcriptome studies in Frankia sp. Genomes from Frankia sp. strains differ markedly in size, a consequence proposed to be associated with a high number of indigenous transposases, more than 200 of which are found in Frankia sp. strain CcI3 used in this study. Because Frankia exhibits a high degree of cell heterogeneity as a consequence of its mycelial growth pattern, its transcriptome is likely to be quite sensitive to culture age. This study focuses on the behavior of the Frankia sp. strain CcI3 transcriptome as a function of nitrogen source and culture age.

Results

To study global transcription in Frankia sp. CcI3 grown under different conditions, complete transcriptomes were determined using high throughput RNA deep sequencing. Samples varied by time (five days vs. three days) and by culture conditions (NH₄⁺ added vs. N₂ fixing). Assembly of millions of reads revealed more diversity of gene expression between five-day and three-day old cultures than between three day old cultures differing in nitrogen sources. Heat map analysis organized genes into groups that were expressed or repressed under the various conditions compared to median expression values. Twenty-one SNPs common to all three transcriptome samples were detected indicating culture heterogeneity in this slow-growing organism. Significantly higher expression of transposase ORFs was found in the five-day and N₂-fixing cultures, suggesting that N starvation and culture aging provide conditions for on-going genome modification. Transposases have previously been proposed to participate in the creating the large number of gene duplication or deletion in host strains. Subsequent RT-qPCR experiments confirmed predicted elevated transposase expression levels indicated by the mRNA-seq data.

Conclusions

The overall pattern of gene expression in aging cultures of CcI3 suggests significant cell heterogeneity even during normal growth on ammonia. The detection of abundant transcription of nif (nitrogen fixation) genes likely reflects the presence of anaerobic, N-depleted microsites in the growing mycelium of the culture, and the presence of significantly elevated transposase transcription during starvation indicates the continuing evolution of the Frankia sp. strain CcI3 genome, even in culture, especially under stressed conditions. These studies also sound a cautionary note when comparing the transcriptomes of Frankia grown in root nodules, where cell heterogeneity would be expected to be quite high.

Identification of TTA codon containing genes in Frankia and exploration of the role of tRNA in regulating these genes

The TTA codon, one of the six available codons for the amino acid leucine, is the rarest codon among the high GC genomes of Actinobacteria including Frankia. This codon has been implicated in various regulatory mechanisms involving secondary metabolism and morphological development. TTA-mediated gene regulation is well documented in Streptomyces coelicolor, but that role has not been investigated in other Actinobacteria including Frankia. Among the various Actinomycetes with a GC content of more than 70%, Frankia genomes had the highest percentages of TTA-containing genes ranging from 5.2 to 10.68% of the genome. In contrast, TTA-bearing genes comprised 1.7, 3.4 and 4.1% of the Streptomyces coelicolor, S. avermitilis and Nocardia farcinia genomes, respectively. We analyzed their functional role, evolutionary significance, horizontal acquisition and the codon-anticodon interaction. The TTA-bearing genes were found to be well represented in metabolic genes involved in amino acid transport and secondary metabolism. A reciprocal Blast search reveal that many of the TTA-bearing genes have orthologs in the other Frankia genomes, and some of these orthologous genes also have a TTA codon in them. The gene expression level of TTA-containing genes was estimated by the use of the codon adaption index (CAI), and the CAI values were found to have a positive correlation with the GC3 (GC content at the 3rd codon position). A full-atomic 3D model of the leucine tRNA recognizing the TTA (UUA) codon was generated and utilized for in silico docking to determine binding affinity in codon-anticodon interaction. We found a proficient codon-anticodon interaction for this codon which is perhaps why so many genes hold on to this rare codon without compromising their translational efficiency.

Genomic insights into bifidobacteria

Since the discovery in 1899 of bifidobacteria as numerically dominant microbes in the feces of breast-fed infants, there have been numerous studies addressing their role in modulating gut microflora as well as their other potential health benefits. Because of this, they are frequently incorporated into foods as probiotic cultures. An understanding of their full interactions with intestinal microbes and the host is needed to scientifically validate any health benefits they may afford. Recently, the genome sequences of nine strains representing four species of Bifidobacterium became available. A comparative genome analysis of these genomes reveals a likely efficient capacity to adapt to their habitats, with B. longum subsp. infantis exhibiting more genomic potential to utilize human milk oligosaccharides, consistent with its habitat in the infant gut. Conversely, B. longum subsp. longum exhibits a higher genomic potential for utilization of plant-derived complex carbohydrates and polyols, consistent with its habitat in an adult gut. An intriguing observation is the loss of much of this genome potential when strains are adapted to pure culture environments, as highlighted by the genomes of B. animalis subsp. lactis strains, which exhibit the least potential for a gut habitat and are believed to have evolved from the B. animalis species during adaptation to dairy fermentation environments.