Evolution of lactic acid bacteria in the order Lactobacillales as depicted by analysis of glycolysis and pentose phosphate pathways

Lactic acid bacteria (LAB) represent a functional group of bacteria that are fundamental in human nutrition because of their prominent role in fermented food production and their presence as commensals in the gut. LAB co-evolution and niche-adaptation have been analyzed in several phylogenomic studies due to the availability of complete genome sequences. The aim of this study was to provide novel insights into LAB evolution through the comparative analysis of the metabolic pathways related to carbohydrate metabolism. The analysis was based on 42 LAB genome sequences of representative strains belonging to Enterococcaceae, Lactobacillaceae, Leuconostocaceae and Streptococcaceae. A reference phylogenetic tree was inferred from concatenation of 42 ribosomal proteins; then 42 genes belonging to the Embden–Meyerhof–Parnas (or glycolysis; EMPP) and pentose phosphate (PPP) pathways were analyzed in terms of their distribution and organization in the genomes. Phylogenetic analyses confirmed the paraphyly of the Lactobacillaceae family, while the distribution and organization of the EMPP and PPP genes revealed the occurrence of lineage-specific trends of gene loss/gain within the two metabolic pathways examined. In addition, the investigation of the two pathways as structures resulting from different evolutionary processes provided new information concerning the genetic bases of heterofermentative/homofermentative metabolism.     

Enhanced Mutation Rate,Relaxed Selection,and the “Domino Effect” are associated with Gene Loss in Blattabacterium,A Cockroach Endosymbiont

Intracellular endosymbionts have reduced genomes that progressively lose genes at a timescale of tens of million years. We previously reported that gene loss rate is linked to mutation rate in Blattabacterium, however, the mechanisms causing gene loss are not yet fully understood. Here, we carried out comparative genomic analyses on the complete genome sequences of a representative set of 67 Blattabacterium strains, with sizes ranging between 511 and 645 kb. We found that 200 of the 566 analyzed protein-coding genes were lost in at least one lineage of Blattabacterium, with the most extreme case being one gene that was lost independently in 24 lineages. We found evidence for three mechanisms influencing gene loss in Blattabacterium. First, gene loss rates were found to increase exponentially with the accumulation of substitutions. Second, genes involved in vitamin and amino acid metabolism experienced relaxed selection in Cryptocercus and Mastotermes, possibly triggered by their vertically inherited gut symbionts. Third, we found evidence of epistatic interactions among genes leading to a “domino effect” of gene loss within pathways. Our results highlight the complexity of the process of genome erosion in an endosymbiont.     

A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes

Background

Sequencing the genomes of multiple, taxonomically diverse eukaryotes enables in-depth comparative-genomic analysis which is expected to help in reconstructing ancestral eukaryotic genomes and major events in eukaryotic evolution and in making functional predictions for currently uncharacterized conserved genes.

Results

We examined functional and evolutionary patterns in the recently constructed set of 5,873 clusters of predicted orthologs (eukaryotic orthologous groups or KOGs) from seven eukaryotic genomes: Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens, Arabidopsis thaliana, Saccharomyces cerevisiae, Schizosaccharomyces pombe and Encephalitozoon cuniculi. Conservation of KOGs through the phyletic range of eukaryotes strongly correlates with their functions and with the effect of gene knockout on the organism's viability. The approximately 40% of KOGs that are represented in six or seven species are enriched in proteins responsible for housekeeping functions, particularly translation and RNA processing. These conserved KOGs are often essential for survival and might approximate the minimal set of essential eukaryotic genes. The 131 single-member, pan-eukaryotic KOGs we identified were examined in detail. For around 20 that remained uncharacterized, functions were predicted by in-depth sequence analysis and examination of genomic context. Nearly all these proteins are subunits of known or predicted multiprotein complexes, in agreement with the balance hypothesis of evolution of gene copy number. Other KOGs show a variety of phyletic patterns, which points to major contributions of lineage-specific gene loss and the 'invention' of genes new to eukaryotic evolution. Examination of the sets of KOGs lost in individual lineages reveals co-elimination of functionally connected genes. Parsimonious scenarios of eukaryotic genome evolution and gene sets for ancestral eukaryotic forms were reconstructed. The gene set of the last common ancestor of the crown group consists of 3,413 KOGs and largely includes proteins involved in genome replication and expression, and central metabolism. Only 44% of the KOGs, mostly from the reconstructed gene set of the last common ancestor of the crown group, have detectable homologs in prokaryotes; the remainder apparently evolved via duplication with divergence and invention of new genes.

Conclusions

The KOG analysis reveals a conserved core of largely essential eukaryotic genes as well as major diversification and innovation associated with evolution of eukaryotic genomes. The results provide quantitative support for major trends of eukaryotic evolution noticed previously at the qualitative level and a basis for detailed reconstruction of evolution of eukaryotic genomes and biology of ancestral forms.     

Effect of topology of quorum sensing-related genes in Pectobacterium atrosepticum on their expression

In prokaryotic genomes, the neighboring genes are often located on the complementary DNA strands and adjoin each other by their 5′- or 3′-ends, or even overlap with their open reading frames. It was suggested that this gene topology has a functional purpose of regulating their expression. For the genes that overlap by their coding 3′-end encoding regions, this assumption has not been confirmed experimentally. In a broad group of bacteria that belong to proteobacteria, this convergent gene arrangement is typical for functionally connected quorum sensing-related genes “I” and “R,” which encode synthases of autoinducers, such as N-acyl homoserine lactones and their sensors, respectively. In the present study on the example of overlapping quorum sensing-related genes of plant pathogenic bacterium Pectobacterium atrosepticum SCRI1043, expI, and expR, it was shown that the topology of these genes determines the regulation of their expression.     

A Genome-Wide Profiling Strategy as an Aid for Searching Unique Identification Biomarkers for <Emphasis Type="Italic">Streptococcus</Emphasis>

The use of rrs (16S rRNA) gene is widely regarded as the “gold standard” for identifying bacteria and determining their phylogenetic relationships. Nevertheless, multiple copies of this gene in a genome is likely to give an overestimation of the bacterial diversity. In each of the 50 Streptococcus genomes (16 species, 50 strains), 4–7 copies of rrs are present. The nucleotide sequences of these rrs genes show high similarity within and among genomes, which did not allow unambiguous identification. A genome-wide search revealed the presence of 27 gene sequences common to all the Streptococcus species. Digestion of these 27 gene sequences with 10 type II restriction endonucleases (REs) showed that unique RE digestion in purH gene is sufficient for clear cut identification of 30 genomes belonging to 16 species. Additional gene-RE combinations allowed identification of another 15 strains belonging to S. pneumoniae, S. pyogenes, and S. suis. For the rest 5 strains, a combination of 2 genes was required for identifying them. The proposed strategy is likely to prove helpful in proper detection of pathogens like Streptococcus.

Electronic supplementary material

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

Random mutagenesis in Corynebacterium glutamicum ATCC 13032 using an IS6100-based transposon vector identified the last unknown gene in the histidine biosynthesis pathway

Background

Chlamydiae species are of much importance from a clinical viewpoint. Their diversity both in terms of their numbers as well as clinical involvement are presently believed to be significantly underestimated. The obligate intracellular nature of chlamydiae has also limited their genetic and biochemical studies. Thus, it is of importance to develop additional means for their identification and characterization.

Results

We have carried out analyses of available chlamydiae genomes to identify sets of unique proteins that are either specific for all Chlamydiales genomes, or different Chlamydiaceae family members, or members of the Chlamydia and Chlamydophila genera, or those unique to Protochlamydia amoebophila, but which are not found in any other bacteria. In total, 59 Chlamydiales-specific proteins, 79 Chlamydiaceae-specific proteins, 20 proteins each that are specific for both Chlamydia and Chlamydophila and 445 ORFs that are Protochlamydia-specific were identified. Additionally, 33 cases of possible gene loss or lateral gene transfer were also detected.

Conclusion

The identified chlamydiae-lineage specific proteins, many of which are highly conserved, provide novel biomarkers that should prove of much value in the diagnosis of these bacteria and in exploration of their prevalence and diversity. These conserved protein sequences (CPSs) also provide novel therapeutic targets for drugs that are specific for these bacteria. Lastly, functional studies on these chlamydiae or chlamydiae subgroup-specific proteins should lead to important insights into lineage-specific adaptations with regards to development, infectivity and pathogenicity.     

Strain-level genomic variation in natural populations of Lebetimonas from an erupting deep-sea volcano

Chemolithoautotrophic Epsilonproteobacteria are ubiquitous in sulfidic, oxygen-poor habitats, including hydrothermal vents, marine oxygen minimum zones, marine sediments and sulfidic caves and have a significant role in cycling carbon, hydrogen, nitrogen and sulfur in these environments. The isolation of diverse strains of Epsilonproteobacteria and the sequencing of their genomes have revealed that this group has the metabolic potential to occupy a wide range of niches, particularly at dynamic deep-sea hydrothermal vents. We expand on this body of work by examining the population genomics of six strains of Lebetimonas, a vent-endemic, thermophilic, hydrogen-oxidizing Epsilonproteobacterium, from a single seamount in the Mariana Arc. Using Lebetimonas as a model for anaerobic, moderately thermophilic organisms in the warm, anoxic subseafloor environment, we show that genomic content is highly conserved and that recombination is limited between closely related strains. The Lebetimonas genomes are shaped by mobile genetic elements and gene loss as well as the acquisition of novel functional genes by horizontal gene transfer, which provide the potential for adaptation and microbial speciation in the deep sea. In addition, these Lebetimonas genomes contain two operons of nitrogenase genes with different evolutionary origins. Lebetimonas expressed nifH during growth with nitrogen gas as the sole nitrogen source, thus providing the first evidence of nitrogen fixation in any Epsilonproteobacteria from deep-sea hydrothermal vents. In this study, we provide a comparative overview of the genomic potential within the Nautiliaceae as well as among more distantly related hydrothermal vent Epsilonproteobacteria to broaden our understanding of microbial adaptation and diversity in the deep sea.     

Knowledge-based discovery for designing CRISPR-CAS systems against invading mobilomes in thermophiles

Clustered regularly interspaced short palindromic repeats (CRISPRs) are direct features of the prokaryotic genomes involved in resistance to their bacterial viruses and phages. Herein, we have identified CRISPR loci together with CRISPR-associated sequences (CAS) genes to reveal their immunity against genome invaders in the thermophilic archaea and bacteria. Genomic survey of this study implied that genomic distribution of CRISPR-CAS systems was varied from strain to strain, which was determined by the degree of invading mobiloms. Direct repeats found to be equal in some extent in many thermopiles, but their spacers were differed in each strain. Phylogenetic analyses of CAS superfamily revealed that genes cmr, csh, csx11, HD domain, devR were belonged to the subtypes of cas gene family. The members in cas gene family of thermophiles were functionally diverged within closely related genomes and may contribute to develop several defense strategies. Nevertheless, genome dynamics, geological variation and host defense mechanism were contributed to share their molecular functions across the thermophiles. A thermophilic archaean, Thermococcus gammotolerans and thermophilic bacteria, Petrotoga mobilis and Thermotoga lettingae have shown superoperons-like appearance to cluster cas genes, which were typically evolved for their defense pathways. A cmr operon was identified with a specific promoter in a thermophilic archaean, Caldivirga maquilingensis. Overall, we concluded that knowledge-based genomic survey and phylogeny-based functional assignment have suggested for designing a reliable genetic regulatory circuit naturally from CRISPR-CAS systems, acquired defense pathways, to thermophiles in future synthetic biology.

Electronic supplementary material

The online version of this article (doi:10.1007/s11693-015-9176-8) contains supplementary material, which is available to authorized users.     

Functional Analysis of the Twin-Arginine Translocation Pathway in Corynebacterium glutamicum ATCC 13869

Compared to those of other gram-positive bacteria, the genetic structure of the Corynebacterium glutamicum Tat system is unique in that it contains the tatE gene in addition to tatA, tatB, and tatC. The tatE homologue has been detected only in the genomes of gram-negative enterobacteria. To assess the function of the C. glutamicum Tat pathway, we cloned the tatA, tatB, tatC, and tatE genes from C. glutamicum ATCC 13869 and constructed mutants carrying deletions of each tat gene or of both the tatA and tatE genes. Using green fluorescent protein (GFP) fused with the twin-arginine signal peptide of the Escherichia coli TorA protein, we demonstrated that the minimal functional Tat system required TatA and TatC. TatA and TatE provide overlapping function. Unlike the TatB proteins from gram-negative bacteria, C. glutamicum TatB was dispensable for Tat function, although it was required for maximal efficiency of secretion. The signal peptide sequence of the isomaltodextranase (IMD) of Arthrobacter globiformis contains a twin-arginine motif. We showed that both IMD and GFP fused with the signal peptide of IMD were secreted via the C. glutamicum Tat pathway. These observations indicate that IMD is a bona fide Tat substrate and imply great potential of the C. glutamicum Tat system for industrial production of heterologous folded proteins.     

Evolutionary aspects of a unique internal mitochondrial targeting signal in nuclear-migrated rps19 of sugar beet (Beta vulgaris L.)

The endosymbiotic theory postulates that many genes migrated from endosymbionts to the nuclear genomes of their hosts. Some migrated genes lack presequences directing proteins to mitochondria, and their mitochondrial targeting signals appear to be inscribed in the core coding regions as internal targeting signals (ITSs). ITSs may have evolved after sequence transfer to nuclei or ITSs may have pre-existed before sequence transfer. Here, we report the molecular cloning of a sugar beet gene for ribosomal protein S19 (Rps19; the first letter is capitalized when the gene is a nuclear gene). We show that sugar beet Rps19 (BvRps19) is an ITS-type gene. Based on amino-acid sequence comparison, dicotyledonous rps19s (the first letter is lower-cased when the gene is a mitochondrial gene), such as tobacco rps19 (Ntrps19), resemble an ancestral form of BvRps19. We investigated whether differences in amino-acid sequences between BvRps19 and Ntrps19 were involved in ITS evolution. Analyses of the intracellular localization of chimaeric GFP-fusion proteins that were transiently expressed in Welsh onion cells showed that Ntrps19-gfp was not localized in mitochondria. When several BvRps19-type amino acid substitutions, none of which was seen in any other angiosperm rps19, were introduced into Ntrps19-gfp, the modified Ntrps19-gfp became localized in mitochondria, supporting the notion that an ITS in BvRps19 evolved following sequence transfer to nuclei. Not all of these substitutions were seen in other ITS-type Rps19s, suggesting that the ITSs of Rps19 are diverse.     

Gene decay in Shigella as an incipient stage of host-adaptation

Background

Many facultative bacterial pathogens have undergone extensive gene decay processes, possibly due to lack of selection pressure during evolutionary conversion from free-living to intracellular lifestyle. Shigella, the causative agents of human shigellosis, have arisen from different E. coli-like ancestors independently by convergent paths. As these bacteria all have lost large numbers of genes by mutation or deletion, they can be used as ideal models for systematically studying the process of gene function loss in different bacteria living under similar selection pressures.

Methodologies/Principal Findings

We compared the sequenced Shigella genomes and re-defined decayed genes (pseudogenes plus deleted genes) in these bacteria. Altogether, 85 genes are commonly decayed in the five analyzed Shigella strains and 1456 genes are decayed in at least one Shigella strain. Genes coding for carbon utilization, cell motility, transporter or membrane proteins are prone to be inactivated. Decayed genes tend to concentrate in certain operons rather than distribute averagely across the whole genome. Genes in the decayed operon accumulated more non-synonymous mutations than the rest genes and meanwhile have lower expression levels.

Conclusions

Different Shigella lineages underwent convergent gene decay processes, and inactivation of one gene would lead to a lesser selection pressure for the other genes in the same operon. The pool of superfluous genes for Shigella may contain at least two thousand genes and the gene decay processes may still continue in Shigella until a minimum genome harboring only essential genes is reached.     

Chromosomal mapping,differential origin and evolution of the S100 gene family

S100 proteins are calcium-binding proteins, which exist only in vertebrates and which constitute a large protein family. The origin and evolution of the S100 family in vertebrate lineages remain a challenge. Here, we examined the synteny conservation of mammalian S100A genes by analysing the sequence of available vertebrate S100 genes in databases. Five S100A gene members, unknown previously, were identified by chromosome mapping analysis. Mammalian S100A genes are duplicated and clustered on a single chromosome while two S100A gene clusters are found on separate chromosomes in teleost fish, suggesting that S100A genes existed in fish before the fish-specific genome duplication took place. During speciation, tandem gene duplication events within the cluster of S100A genes of a given chromosome have probably led to the multiple members of the S100A gene family. These duplicated genes have been retained in the genome either by neofunctionalisation and/or subfunctionalisation or have evolved into non-coding sequences. However in vertebrate genomes, other S100 genes are also present i.e. S100P, S100B, S100G and S100Z, which exist as single copy genes distributed on different chromosomes, suggesting that they could have evolved from an ancestor different to that of the S100A genes.     

Comparative Analysis of the Mosaic Genomes of Tailed Archaeal Viruses and Proviruses Suggests Common Themes for Virion Architecture and Assembly with Tailed Viruses of Bacteria

Tailed double-stranded DNA viruses (order Caudovirales) represent the dominant morphotype among viruses infecting bacteria. Analysis and comparison of complete genome sequences of tailed bacterial viruses provided insights into their origin and evolution. Structural and genomic studies have unexpectedly revealed that tailed bacterial viruses are evolutionarily related to eukaryotic herpesviruses. Organisms from the third domain of life, Archaea, are also infected by viruses that, in their overall morphology, resemble tailed viruses of bacteria. However, high-resolution structural information is currently unavailable for any of these viruses, and only a few complete genomes have been sequenced so far. Here we identified nine proviruses that are clearly related to tailed bacterial viruses and integrated into chromosomes of species belonging to four different taxonomic orders of the Archaea. This more than doubled the number of genome sequences available for comparative studies. Our analyses indicate that highly mosaic tailed archaeal virus genomes evolve by homologous and illegitimate recombination with genomes of other viruses, by diversification, and by acquisition of cellular genes. Comparative genomics of these viruses and related proviruses revealed a set of conserved genes encoding putative proteins similar to virion assembly and maturation, as well as genome packaging proteins of tailed bacterial viruses and herpesviruses. Furthermore, fold prediction and structural modeling experiments suggest that the major capsid proteins of tailed archaeal viruses adopt the same topology as the corresponding proteins of tailed bacterial viruses and eukaryotic herpesviruses. Data presented in this study strongly support the hypothesis that tailed viruses infecting archaea share a common ancestry with tailed bacterial viruses and herpesviruses.     

Extensive intra-phylotype diversity in lactobacilli and bifidobacteria from the honeybee gut

Background

In the honeybee Apis mellifera, the bacterial gut community is consistently colonized by eight distinct phylotypes of bacteria. Managed bee colonies are of considerable economic interest and it is therefore important to elucidate the diversity and role of this microbiota in the honeybee. In this study, we have sequenced the genomes of eleven strains of lactobacilli and bifidobacteria isolated from the honey crop of the honeybee A. mellifera.

Results

Single gene phylogenies confirmed that the isolated strains represent the diversity of lactobacilli and bifidobacteria in the gut, as previously identified by 16S rRNA gene sequencing. Core genome phylogenies of the lactobacilli and bifidobacteria further indicated extensive divergence between strains classified as the same phylotype. Phylotype-specific protein families included unique surface proteins. Within phylotypes, we found a remarkably high level of gene content diversity. Carbohydrate metabolism and transport functions contributed up to 45% of the accessory genes, with some genomes having a higher content of genes encoding phosphotransferase systems for the uptake of carbohydrates than any previously sequenced genome. These genes were often located in highly variable genomic segments that also contained genes for enzymes involved in the degradation and modification of sugar residues. Strain-specific gene clusters for the biosynthesis of exopolysaccharides were identified in two phylotypes. The dynamics of these segments contrasted with low recombination frequencies and conserved gene order structures for the core genes. Hits for CRISPR spacers were almost exclusively found within phylotypes, suggesting that the phylotypes are associated with distinct phage populations.

Conclusions

The honeybee gut microbiota has been described as consisting of a modest number of phylotypes; however, the genomes sequenced in the current study demonstrated a very high level of gene content diversity within all three described phylotypes of lactobacilli and bifidobacteria, particularly in terms of metabolic functions and surface structures, where many features were strain-specific. Together, these results indicate niche differentiation within phylotypes, suggesting that the honeybee gut microbiota is more complex than previously thought.

Electronic supplementary material

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