Phylogenetics and the cohesion of bacterial genomes

Research progress in isolation and identification of rumen probiotics

With the increasing research on the exploitation of rumen microbial resources, rumen probiotics have attracted much attention for their positive contributions in promoting nutrient digestion, inhibiting pathogenic bacteria, and improving production performance. In the past two decades, macrogenomics has provided a rich source of new-generation probiotic candidates, but most of these "dark substances" have not been successfully cultured due to the restrictive growth conditions. However, fueled by high-throughput culture and sorting technologies, it is expected that the potential probiotics in the rumen can be exploited on a large scale, and their potential applications in medicine and agriculture can be explored. In this paper, we review and summarize the classical techniques for isolation and identification of rumen probiotics, introduce the development of droplet-based high-throughput cell culture and single-cell sequencing for microbial culture and identification, and finally introduce promising cultureomics techniques. The aim is to provide technical references for the development of related technologies and microbiological research to promote the further development of the field of rumen microbiology research.

Molecular Accounting and Profiling of Human Respiratory Microbial Communities: Toward Precision Medicine by Targeting the Respiratory Microbiome for Disease Diagnosis and Treatment

The wide diversity of microbiota at the genera and species levels across sites and individuals is related to various causes and the observed differences between individuals. Efforts are underway to further understand and characterize the human-associated microbiota and its microbiome. Using 16S rDNA as a genetic marker for bacterial identification improved the detection and profiling of qualitative and quantitative changes within a bacterial population. In this light, this review provides a comprehensive overview of the basic concepts and clinical applications of the respiratory microbiome, alongside an in-depth explanation of the molecular targets and the potential relationship between the respiratory microbiome and respiratory disease pathogenesis. The paucity of robust evidence supporting the correlation between the respiratory microbiome and disease pathogenesis is currently the main challenge for not considering the microbiome as a novel druggable target for therapeutic intervention. Therefore, further studies are needed, especially prospective studies, to identify other drivers of microbiome diversity and to better understand the changes in the lung microbiome along with the potential association with disease and medications. Thus, finding a therapeutic target and unfolding its clinical significance would be crucial.

Genomic diversity, pathogenicity and antimicrobial resistance of Escherichia coli isolated from poultry in the southern United States

Escherichia coli (E. coli) are typically present as commensal bacteria in the gastro-intestinal tract of most animals including poultry species, but some avian pathogenic E. coli (APEC) strains can cause localized and even systematic infections in domestic poultry. Emergence and re-emergence of antimicrobial resistant isolates (AMR) constrain antibiotics usage in poultry production, and development of an effective vaccination program remains one of the primary options in E. coli disease prevention and control for domestic poultry. Thus, understanding genetic and pathogenic diversity of the enzootic E. coli isolates, particularly APEC, in poultry farms is the key to designing an optimal vaccine candidate and to developing an effective vaccination program. This study explored the genomic and pathogenic diversity among E. coli isolates in southern United States poultry. A total of nine isolates were recovered from sick broilers from Mississippi, and one from Georgia, with epidemiological variations among clinical signs, type of housing, and bird age. The genomes of these isolates were sequenced by using both Illumina short-reads and Oxford Nanopore long-reads, and our comparative analyses suggested data from both platforms were highly consistent. The 16 s rRNA based phylogenetic analyses showed that the 10 bacteria strains are genetically closer to each other than those in the public database. However, whole genome analyses showed that these 10 isolates encoded a diverse set of reported virulence and AMR genes, belonging to at least nine O:H serotypes, and are genetically clustered with at least five different groups of E. coli isolates reported by other states in the United States. Despite the small sample size, this study suggested that there was a large extent of genomic and serological diversity among E. coli isolates in southern United States poultry. A large-scale comprehensive study is needed to understand the overall genomic diversity and the associated virulence, and such a study will be important to develop a broadly protective E. coli vaccine.

Oral microbiota may predict the presence of esophageal squamous cell carcinoma

PURPOSE

Microbial imbalances have been well elucidated in esophageal adenocarcinoma (EAC), but few studies address the oral microbiota in esophageal squamous cell carcinoma (ESCC). In view of the fact, we aimed to explore the associations of oral microbiota with these patients suffering from ESCC.

METHODS

In our study, a total of 109 individuals were enrolled (control = 53, ESCC = 56). We profiled the microbiota in oral swabs from individuals with control (ConT) and ESCC (ESCCT). 16S rRNA gene sequencing was applied to analyze the microbiome. The α and β diversity differences were tested by Tukey Test and Partial Least Squares Discriminant Analysis (PLS-DA) respectively. Linear discriminant analysis effect size (LEfSe) analysis was performed to assess taxonomic differences between the two groups.

RESULTS

Our results showed that the microbial richness and diversity was a slightly higher in ESCCT groups than that in ConT groups. Bacteroidota, Firmicutes, Proteobacteria, Fusobacteria, Actinobacteria and Patescibacteria were the six dominant bacteria of oral flora in the two groups. When compared with control group, increased Fusobacterioa at phylum level, Neisseriaceae at family level and Leptotrichia at genus level were detected. LEfSe analysis indicated a greater abundance of Leptotrichiaceae, Leptotrichia, Fusobacteriales, Fusobacteria and Fusobacteriota in ESCC groups.

CONCLUSION

Our study suggests a potential association between oral microbiome dysbiosis and ESCC and provides insights on a potential screening marker for esophageal cancer.

Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies

BACKGROUND

The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82].

RESULTS

At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%).

CONCLUSIONS

Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.

Classes of explicit phylogenetic networks and their biological and mathematical significance

The evolutionary relationships among organisms have traditionally been represented using rooted phylogenetic trees. However, due to reticulate processes such as hybridization or lateral gene transfer, evolution cannot always be adequately represented by a phylogenetic tree, and rooted phylogenetic networks that describe such complex processes have been introduced as a generalization of rooted phylogenetic trees. In fact, estimating rooted phylogenetic networks from genomic sequence data and analyzing their structural properties is one of the most important tasks in contemporary phylogenetics. Over the last two decades, several subclasses of rooted phylogenetic networks (characterized by certain structural constraints) have been introduced in the literature, either to model specific biological phenomena or to enable tractable mathematical and computational analyses. In the present manuscript, we provide a thorough review of these network classes, as well as provide a biological interpretation of the structural constraints underlying these networks where possible. In addition, we discuss how imposing structural constraints on the network topology can be used to address the scalability and identifiability challenges faced in the estimation of phylogenetic networks from empirical data.

Horizontal gene transfer and adaptive evolution in bacteria

Horizontal gene transfer (HGT) is arguably the most conspicuous feature of bacterial evolution. Evidence for HGT is found in most bacterial genomes. Although HGT can considerably alter bacterial genomes, not all transfer events may be biologically significant and may instead represent the outcome of an incessant evolutionary process that only occasionally has a beneficial purpose. When adaptive transfers occur, HGT and positive selection may result in specific, detectable signatures in genomes, such as gene-specific sweeps or increased transfer rates for genes that are ecologically relevant. In this Review, we first discuss the various mechanisms whereby HGT occurs, how the genetic signatures shape patterns of genomic variation and the distinct bioinformatic algorithms developed to detect these patterns. We then discuss the evolutionary theory behind HGT and positive selection in bacteria, and discuss the approaches developed over the past decade to detect transferred DNA that may be involved in adaptation to new environments.

Modular evolution of secretion systems and virulence plasmids in a bacterial species complex

BACKGROUND

Many named species as defined in current bacterial taxonomy correspond to species complexes. Uncertainties regarding the organization of their genetic diversity challenge research efforts. We utilized the Agrobacterium tumefaciens species complex (a.k.a. Agrobacterium biovar 1), a taxon known for its phytopathogenicity and applications in transformation, as a study system and devised strategies for investigating genome diversity and evolution of species complexes.

RESULTS

We utilized 35 genome assemblies, including 14 newly generated ones, to achieve a phylogenetically balanced sampling of A. tumefaciens. Our genomic analysis suggested that the 10 genomospecies described previously are distinct biological species and supported a quantitative guideline for species delineation. Furthermore, our inference of gene content and core-genome phylogeny allowed for investigations of genes critical in fitness and ecology. For the type VI secretion system (T6SS) involved in interbacterial competition and thought to be conserved, we detected multiple losses and one horizontal gene transfer. For the tumor-inducing plasmids (pTi) and pTi-encoded type IV secretion system (T4SS) that are essential for agrobacterial phytopathogenicity, we uncovered novel diversity and hypothesized their involvement in shaping this species complex. Intriguingly, for both T6SS and T4SS, genes encoding structural components are highly conserved, whereas extensive diversity exists for genes encoding effectors and other proteins.

CONCLUSIONS

We demonstrate that the combination of a phylogeny-guided sampling scheme and an emphasis on high-quality assemblies provides a cost-effective approach for robust analysis in evolutionary genomics. We show that the T6SS VgrG proteins involved in specific effector binding and delivery can be classified into distinct types based on domain organization. The co-occurrence patterns of VgrG-associated domains and the neighboring genes that encode different chaperones/effectors can be used to infer possible interacting partners. Similarly, the associations between plant host preference and the pTi type among these strains can be used to infer phenotype-genotype correspondence. Our strategies for multi-level investigations at scales that range from whole genomes to intragenic domains and phylogenetic depths from between- to within-species are applicable to other bacteria. Furthermore, modularity observed in the molecular evolution of genes and domains is useful for inferring functional constraints and informing experimental works.

ToRQuEMaDA: tool for retrieving queried Eubacteria, metadata and dereplicating assemblies

TQMD is a tool for high-performance computing clusters which downloads, stores and produces lists of dereplicated prokaryotic genomes. It has been developed to counter the ever-growing number of prokaryotic genomes and their uneven taxonomic distribution. It is based on word-based alignment-free methods (k-mers), an iterative single-linkage approach and a divide-and-conquer strategy to remain both efficient and scalable. We studied the performance of TQMD by verifying the influence of its parameters and heuristics on the clustering outcome. We further compared TQMD to two other dereplication tools (dRep and Assembly-Dereplicator). Our results showed that TQMD is primarily optimized to dereplicate at higher taxonomic levels (phylum/class), as opposed to the other dereplication tools, but also works at lower taxonomic levels (species/strain) like the other dereplication tools. TQMD is available from source and as a Singularity container at [https://bitbucket.org/phylogeno/tqmd ].

Exploration of space to achieve scientific breakthroughs

Living organisms adapt to changing environments using their amazing flexibility to remodel themselves by a process called evolution. Environmental stress causes selective pressure and is associated with genetic and phenotypic shifts for better modifications, maintenance, and functioning of organismal systems. The natural evolution process can be used in complement to rational strain engineering for the development of desired traits or phenotypes as well as for the production of novel biomaterials through the imposition of one or more selective pressures. Space provides a unique environment of stressors (e.g., weightlessness and high radiation) that organisms have never experienced on Earth. Cells in the outer space reorganize and develop or activate a range of molecular responses that lead to changes in cellular properties. Exposure of cells to the outer space will lead to the development of novel variants more efficiently than on Earth. For instance, natural crop varieties can be generated with higher nutrition value, yield, and improved features, such as resistance against high and low temperatures, salt stress, and microbial and pest attacks. The review summarizes the literature on the parameters of outer space that affect the growth and behavior of cells and organisms as well as complex colloidal systems. We illustrate an understanding of gravity-related basic biological mechanisms and enlighten the possibility to explore the outer space environment for application-oriented aspects. This will stimulate biological research in the pursuit of innovative approaches for the future of agriculture and health on Earth.

Genome Sequencing of Mesonia algae K4-1 Reveals Its Adaptation to the Arctic Ocean

The special ecological environment of the Arctic has brought about a large number of salt-tolerant and psychrotolerant microorganisms. We isolated two culturable bacterial strains of the genus Mesonia; one from the Arctic ocean, Mesonia algae K4-1, and one from the tropical sea, Mesonia sp. HuA40. Our genome analysis and phenotypic experiments indicated that Mesonia algae K4-1 is a moderately halophilic and psychrophilic bacterium. Mesonia algae K4-1 can tolerate 3–14% NaCl and grow at a wide range of temperatures from 4 to 50°C. Mesonia sp. HuA40 is a mesophilic bacterium that can only grow with 3–9% NaCl. In addition, the salt adaptation strategy of Mesonia algae K4-1 accumulates organic osmolytes in the cell. RNA helicases, glutathione and organic compatible solutes may play important roles in maintaining the metabolism and physiological function of Mesonia algae K4-1 under cold stress. Moreover, the ability of Mesonia algae K4-1 to adapt to an oligotrophic marine environment is likely due to the synthesis of a large number of extracellular polysaccharides and the secretion of various families of extracellular proteases. This study systematically analyzed the relationship between genomic differentiation and environmental factors of the Mesonia genus and revealed the possible adaptation mechanism of Mesonia algae K4-1 in the extreme Arctic marine environment at the genomic level.

A Literature Review of Metagenomics and Culturomics of the Peri-implant Microbiome: Current Evidence and Future Perspectives

Background and objectives: In recent years, many different culture-independent molecular techniques have been developed with the aim of investigating the not yet cultivated part of the resident flora of the oral cavity and of analyzing the peri-implant and periodontal flora both in healthy and diseased sites. The most used technologies are Roche 454 pyrosequencing, Illumina HiSeq/MiSeq, ABI SOLiD and Ion Torrent. Due to these methods, two different approaches are available: Metagenomics and the 16S gene analysis. A complementary strategy was also recently developed: Culturomics. Culturomics consists of different culture conditions that allow a very rapid bacterial identification. The focused question of this review was developed in PICO format in order to investigate the role of metagenomics, 16S gene analysis and culturomics (interventions) in the differential study (comparison) of the peri-implant and periodontal microbiome (outcome) in humans (participants). The secondary aim was the characterization of currents limits and future applications of the three techniques. Methods: The authors performed a literature search on three databases (Web of Science, Scopus and PubMed) from 01/01/2003 to 31/06/2019. Date of last search was: 25/08/19. Any type of article dealing with the analysis of periodontal and peri-implant flora with metagenomic, culturomic or 16S gene analysis was included. No language restrictions were applied. Risk of bias for RCT was assessed using the Cochrane collaboration's tool whereas case-control and cohort studies were evaluated through the Newcastle-Ottawa scale. Results: The initial search resulted in 330 titles in total. After careful evaluation of all results no studies were found to satisfy the primary outcome of the present review. Hence a narrative review dealing with the secondary aim was performed. Conclusions: Metagenomic and 16S gene analysis approaches contributed in clarifying some crucial aspects of the oral microbiome. Based on the reported evidence some bacteria could be found around teeth and implants even in the absence of signs of inflammation and other species are more frequently found in supragingival peri-implant biofilm. Teeth and implants (even if adjacent) seem not to share the same microbiome and healthy teeth have a more diversified one. The same analyses also highlighted that the oral biofilm of smokers is composed by more periodontopathogen bacteria compared to non-smokers and that geographical location and ethnicity seem to play a role in bacterial composition. Culturomics, which has not yet been applied to the study of oral microbiota, consists of the use of different culture conditions and of the identification by matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) with the aim of increasing the bacterial repertoire and avoiding the limits of molecular methods. In order to better evaluate perspectives and limits of the all presented approaches further studies comparing the different molecular techniques are encouraged. This review received no funding.

A library of human gut bacterial isolates paired with longitudinal multiomics data enables mechanistic microbiome research

Our understanding of how the gut microbiome interacts with its human host has been restrained by limited access to longitudinal datasets to examine stability and dynamics, and by having only a few isolates to test mechanistic hypotheses. Here, we present the Broad Institute-OpenBiome Microbiome Library (BIO-ML), a comprehensive collection of 7,758 gut bacterial isolates paired with 3,632 genome sequences and longitudinal multi-omics data. We show that microbial species maintain stable population sizes within and across humans and that commonly used 'omics' survey methods are more reliable when using averages over multiple days of sampling. Variation of gut metabolites within people over time is associated with amino acid levels, and differences across people are associated with differences in bile acids. Finally, we show that genomic diversification can be used to infer eco-evolutionary dynamics and in vivo selection pressures for strains within individuals. The BIO-ML is a unique resource designed to enable hypothesis-driven microbiome research.

Molecular identification and phylogenetic relationships of clinical Nocardia isolates

Species identification of Nocardia is difficult because of a complex and rapidly evolving taxonomy. In this study, gene sequencing (16S rRNA, gyrB, secA1, hsp65, rpoB), single 16S rRNA gene sequence phylogenetic analysis, and MALDI-TOF analysis were used to accurately identify 46 clinical Nocardia isolates to the species level. This provided a basis for establishing a routine method of multi-locus sequence analysis (MLSA) for molecular identification of Nocardia species. Genetic polymorphism analysis showed that MLSA was a powerfully discriminating method compared with the 16S rRNA gene to identify clinical Nocardia isolates. However, five-locus (gyrB-16S rRNA-secA1-hsp65-rpoB) MLSA led to misidentifications of all of the five Nocardia abscessus, which were confirmed by digital DNA-DNA hybridization (DDH) analysis. Interestingly, four strains identified as Nocardia beijingensis by a 16S rRNA gene phylogenetic tree may be novel species as suggested by DDH studies. For the purpose of achieving both accuracy and discrimination, the data of MLSA were reanalyzed. A three-locus MLSA with concatenated gyrB-16S rRNA-secA1 sequences was used to construct the phylogenetic tree with high accuracy and powerful discrimination. Therefore, a routine method of MLSA was developed to identify clinical Nocardia species.

The synergistic effect of concatenation in phylogenomics: the case in Pantoea

With the increased availability of genome sequences for bacteria, it has become routine practice to construct genome-based phylogenies. These phylogenies have formed the basis for various taxonomic decisions, especially for resolving problematic relationships between taxa. Despite the popularity of concatenating shared genes to obtain well-supported phylogenies, various issues regarding this combined-evidence approach have been raised. These include the introduction of phylogenetic error into datasets, as well as incongruence due to organism-level evolutionary processes, particularly horizontal gene transfer and incomplete lineage sorting. Because of the huge effect that this could have on phylogenies, we evaluated the impact of phylogenetic conflict caused by organism-level evolutionary processes on the established species phylogeny for Pantoea, a member of the Enterobacterales. We explored the presence and distribution of phylogenetic conflict at the gene partition and nucleotide levels, by identifying putative inter-lineage recombination events that might have contributed to such conflict. Furthermore, we determined whether smaller, randomly constructed datasets had sufficient signal to reconstruct the current species tree hypothesis or if they would be overshadowed by phylogenetic incongruence. We found that no individual gene tree was fully congruent with the species phylogeny of Pantoea, although many of the expected nodes were supported by various individual genes across the genome. Evidence of recombination was found across all lineages within Pantoea, and provides support for organism-level evolutionary processes as a potential source of phylogenetic conflict. The phylogenetic signal from at least 70 random genes recovered robust, well-supported phylogenies for the backbone and most species relationships of Pantoea, and was unaffected by phylogenetic conflict within the dataset. Furthermore, despite providing limited resolution among taxa at the level of single gene trees, concatenated analyses of genes that were identified as having no signal resulted in a phylogeny that resembled the species phylogeny of Pantoea. This distribution of signal and noise across the genome presents the ideal situation for phylogenetic inference, as the topology from a ≥70-gene concatenated species phylogeny is not driven by single genes, and our data suggests that this finding may also hold true for smaller datasets. We thus argue that, by using a concatenation-based approach in phylogenomics, one can obtain robust phylogenies due to the synergistic effect of the combined signal obtained from multiple genes.

Genome-Wide Comparative Analysis of Phylogenetic Trees: The Prokaryotic Forest of Life

Genome-wide comparison of phylogenetic trees is becoming an increasingly common approach in evolutionary genomics, and a variety of approaches for such comparison have been developed. In this article we present several methods for comparative analysis of large numbers of phylogenetic trees. To compare phylogenetic trees taking into account the bootstrap support for each internal branch, the boot-split distance (BSD) method is introduced as an extension of the previously developed split distance (SD) method for tree comparison. The BSD method implements the straightforward idea that comparison of phylogenetic trees can be made more robust by treating tree splits differentially depending on the bootstrap support. Approaches are also introduced for detecting treelike and netlike evolutionary trends in the phylogenetic Forest of Life (FOL), i.e., the entirety of the phylogenetic trees for conserved genes of prokaryotes. The principal method employed for this purpose includes mapping quartets of species onto trees to calculate the support of each quartet topology and so to quantify the tree and net contributions to the distances between species. We describe the applications methods used to analyze the FOL and the results obtained with these methods. These results support the concept of the Tree of Life (TOL) as a central evolutionary trend in the FOL as opposed to the traditional view of the TOL as a "species tree."

Failure to Recover Major Events of Gene Flux in Real Biological Data Due to Method Misapplication

In prokaryotes, known mechanisms of lateral gene transfer (transformation, transduction, conjugation, and gene transfer agents) generate new combinations of genes among chromosomes during evolution. In eukaryotes, whose host lineage is descended from archaea, lateral gene transfer from organelles to the nucleus occurs at endosymbiotic events. Recent genome analyses studying gene distributions have uncovered evidence for sporadic, discontinuous events of gene transfer from bacteria to archaea during evolution. Other studies have used traditional models designed to investigate gene family size evolution (Count) to support claims that gene transfer to archaea was continuous during evolution, rather than involving occasional periodic mass gene influx events. Here, we show that the methodology used in analyses favoring continuous gene transfers to archaea was misapplied in other studies and does not recover known events of single simultaneous origin for many genes followed by differential loss in real data: plastid genomes. Using the same software and the same settings, we reanalyzed presence/absence pattern data for proteins encoded in plastid genomes and for eukaryotic protein families acquired from plastids. Contrary to expectations under a plastid origin model, we found that the methodology employed inferred that gene acquisitions occurred uniformly across the plant tree. Sometimes as many as nine different acquisitions by plastid DNA were inferred for the same protein family. That is, the methodology that recovered gradual and continuous lateral gene transfer among lineages for archaea obtains the same result for plastids, even though it is known that massive gains followed by gradual differential loss is the true evolutionary process that generated plastid gene distribution data. Our findings caution against the use of models designed to study gene family size evolution for investigating gene transfer processes, especially when transfers involving more than one gene per event are possible.

Evolutionary History of Copper Membrane Monooxygenases

Copper membrane monooxygenases (CuMMOs) oxidize ammonia, methane and some short-chain alkanes and alkenes. They are encoded by three genes, usually in an operon of xmoCAB. We aligned xmo operons from 66 microbial genomes, including members of the Alpha-, Beta-, and Gamma-proteobacteria, Verrucomicrobia, Actinobacteria, Thaumarchaeota and the candidate phylum NC10. Phylogenetic and compositional analyses were used to reconstruct the evolutionary history of the enzyme and detect potential lateral gene transfer (LGT) events. The phylogenetic analyses showed at least 10 clusters corresponding to a combination of substrate specificity and bacterial taxonomy, but with no overriding structure based on either function or taxonomy alone. Adaptation of the enzyme to preferentially oxidize either ammonia or methane has occurred more than once. Individual phylogenies of all three genes, xmoA, xmoB and xmoC, closely matched, indicating that this operon evolved or was consistently transferred as a unit, with the possible exception of the methane monooxygenase operons in Verrucomicrobia, where the pmoB gene has a distinct phylogeny from pmoA and pmoC. Compositional analyses indicated that some clusters of xmoCAB operons (for example, the pmoCAB in gammaproteobacterial methanotrophs and the amoCAB in betaproteobacterial nitrifiers) were compositionally very different from their genomes, possibly indicating recent lateral transfer of these operons. The combined phylogenetic and compositional analyses support the hypothesis that an ancestor of the nitrifying bacterium Nitrosococcus was the donor of methane monooxygenase (pMMO) to both the alphaproteobacterial and gammaproteobacterial methanotrophs, but that before this event the gammaproteobacterial methanotrophs originally possessed another CuMMO (Pxm), which has since been lost in many species.

Assessing the Bacterial Community Structure in the Rhizoplane of Wetland Plants

Plant-microorganism interaction in the rhizosphere is important for nutrient cycling, carbon sequestration in natural ecosystems, contaminant elimination and ecosystem functioning. Abundance of microbial communities and variation in species composition can be an imperative determinant of phytoremediation capability. In the present study we have assessed the bacterial community structure in the rhizoplane of wetland plants, Acorus calamus, Typha latifolia, and Phragmites karka using Terminal restriction fragment length polymorphism technique. The most dominant phylum, in the plants under study, was phylum Firmicutes, followed by Proteobacteria and Actinobacteria. Bacterial groups belonging to phylum Chloroflexi, Acidobacteria, Deferribacteres and Thermotogae also showed their presence in P. karka and T. latifolia but were absent in A. calamus. Diversity indices of bacterial community were assessed. The results of this study show the presence of bacterial phyla which play an important role in bioremediation of contaminants. Thus these plants can be used as potential candidates of phytoremediation.

Chromosome architecture constrains horizontal gene transfer in bacteria

Despite significant frequencies of lateral gene transfer between species, higher taxonomic groups of bacteria show ecological and phenotypic cohesion. This suggests that barriers prevent panmictic dissemination of genes via lateral gene transfer. We have proposed that most bacterial genomes have a functional architecture imposed by Architecture IMparting Sequences (AIMS). AIMS are defined as 8 base pair sequences preferentially abundant on leading strands, whose abundance and strand-bias are positively correlated with proximity to the replication terminus. We determined that inversions whose endpoints lie within a single chromosome arm, which would reverse the polarity of AIMS in the inverted region, are both shorter and less frequent near the replication terminus. This distribution is consistent with the increased selection on AIMS function in this region, thus constraining DNA rearrangement. To test the hypothesis that AIMS also constrain DNA transfer between genomes, AIMS were identified in genomes while ignoring atypical, potentially laterally-transferred genes. The strand-bias of AIMS within recently acquired genes was negatively correlated with the distance of those genes from their genome’s replication terminus. This suggests that selection for AIMS function prevents the acquisition of genes whose AIMS are not found predominantly in the permissive orientation. This constraint has led to the loss of at least 18% of genes acquired by transfer in the terminus-proximal region. We used completely sequenced genomes to produce a predictive road map of paths of expected horizontal gene transfer between species based on AIMS compatibility between donor and recipient genomes. These results support a model whereby organisms retain introgressed genes only if the benefits conferred by their encoded functions outweigh the detriments incurred by the presence of foreign DNA lacking genome-wide architectural information.

The potential success of horizontal gene transfer events is historically equated to the benefits conferred by encoded products. Here we show that gene transfer events are observed less frequently if the introduced genes disrupt important patterns of genomic information, suggesting that this disruption would confer an unacceptable cost. As a result, gene transfer events are less likely to be successful if the potential donor genomes have incompatible genome architecture. Because more distantly-related genes are less compatible, chromosome architecture serves as a mechanism to bias gene transfer events to those involving closer relatives, thereby providing a mechanism for the genotypic and phenotypic cohesion of higher taxonomic groups.

16S rRNA Gene Primer Validation for Bacterial Diversity Analysis of Vegetable Products

High-throughput sequencing of the 16S rRNA gene enhances understanding of microbial diversity from complex environmental samples. The 16S rRNA gene is currently the most important target in bacterial evolution and ecology studies, particularly for determination of phylogenetic relationships among taxa, exploration of bacterial diversity in a given environment, and quantification of the relative abundance of taxa at various levels. However, some parts of the conserved region of the bacterial 16S rRNA gene are similar to the conserved regions of plant chloroplasts and eukaryotic mitochondria. Therefore, if DNA contains a large amount of nontarget DNA, this nontarget DNA can be coamplified and consequently produce useless sequence reads. We experimentally assessed the primer pair 335f/769r and the widely used bacterial primer pair SD (S-D-Bact-0341-b-S-17/S-D-Bact-0785-a-A-21). The primer pair 335f/769r was examined for its ability to amplify bacterial DNA in plant and animal feed samples by using the single-strand confirmation polymorphism method. In our present study, these primer pairs were validated for microbial community structure analysis with complex food matrices by using next-generation sequencing. The sequencing results revealed that the primer pair 335f/769r successfully resulted in fewer chloroplast and mitochondrial sequence reads than generated by the universal primer pair SD and therefore is comparatively suitable for metagenomic analyses of complex food matrices, particularly those that are rich in plant DNA. Additionally, some taxonomic groups were missed entirely when only the SD primer pair was used.

Isolation and Characteristics of Biotechnologically Important Antagonistic Thermophilic Bacteria from Rhizosphere of Haloxylon salicornicum

Rhizobacteria are an active part of microbial population in the rhizosphere of plants. In this study, twenty rhizobacteria were isolated from the rhizosphere of a perennial grass, Haloxylon salicornicum, found in Cholistan desert, an arid landmass near Bahawalpur Pakistan, in one set of experimental conditions. Colony characteristics, biochemical and molecular analyses of these isolates were performed. All isolates were bacilli, gram positive with off-white colonies and exhibited typical bacilli colony morphology. None of the isolates was gelatinase, urease, indole, H2S and catalase producer. Eleven isolates were amylase producers and 8 isolates were acid producers. All isolates fermented glucose, 3 fermented lactose and 19 fermented fructose. Molecular data revealed that out of twenty isolates, 14 isolates showed 91-99% identity with Brevibacillus borstelensis, 4 with Bacillus subtilis (97-98%) and 2 with Bacillus licheniformis (94-99%) through BLAST analysis. All identified bacterial isolates cladded with their respective groups in the phylogenetic tree. Many (11-15 out of 20) of the isolates were more effective in inhibiting growth of the tested bacterial strains as compared to the positive control (Ampicillin 50 μg/disc). We conclude that bacilli are the predominant form populating rhizosphere of this desert grass. Among the isolated bacteria Brevibacillus borstelensis, Bacillus subtilis and Bacillus licheniformis are the most predominant species.

Genetic Competence Drives Genome Diversity in Bacillus subtilis

Prokaryote genomes are the result of a dynamic flux of genes, with increases achieved via horizontal gene transfer and reductions occurring through gene loss. The ecological and selective forces that drive this genomic flexibility vary across species. Bacillus subtilis is a naturally competent bacterium that occupies various environments, including plant-associated, soil, and marine niches, and the gut of both invertebrates and vertebrates. Here, we quantify the genomic diversity of B. subtilis and infer the genome dynamics that explain the high genetic and phenotypic diversity observed. Phylogenomic and comparative genomic analyses of 42 B. subtilis genomes uncover a remarkable genome diversity that translates into a core genome of 1,659 genes and an asymptotic pangenome growth rate of 57 new genes per new genome added. This diversity is due to a large proportion of low-frequency genes that are acquired from closely related species. We find no gene-loss bias among wild isolates, which explains why the cloud genome, 43% of the species pangenome, represents only a small proportion of each genome. We show that B. subtilis can acquire xenologous copies of core genes that propagate laterally among strains within a niche. While not excluding the contributions of other mechanisms, our results strongly suggest a process of gene acquisition that is largely driven by competence, where the long-term maintenance of acquired genes depends on local and global fitness effects. This competence-driven genomic diversity provides B. subtilis with its generalist character, enabling it to occupy a wide range of ecological niches and cycle through them.

Highly divergent 16S rRNA sequences in ribosomal operons of Scytonema hyalinum (Cyanobacteria)

A highly divergent 16S rRNA gene was found in one of the five ribosomal operons present in a species complex currently circumscribed as Scytonema hyalinum (Nostocales, Cyanobacteria) using clone libraries. If 16S rRNA sequence macroheterogeneity among ribosomal operons due to insertions, deletions or truncation is excluded, the sequence heterogeneity observed in S. hyalinum was the highest observed in any prokaryotic species thus far (7.3-9.0%). The secondary structure of the 16S rRNA molecules encoded by the two divergent operons was nearly identical, indicating possible functionality. The 23S rRNA gene was examined for a few strains in this complex, and it was also found to be highly divergent from the gene in Type 2 operons (8.7%), and likewise had nearly identical secondary structure between the Type 1 and Type 2 operons. Furthermore, the 16S-23S ITS showed marked differences consistent between operons among numerous strains. Both operons have promoter sequences that satisfy consensus requirements for functional prokaryotic transcription initiation. Horizontal gene transfer from another unknown heterocytous cyanobacterium is considered the most likely explanation for the origin of this molecule, but does not explain the ultimate origin of this sequence, which is very divergent from all 16S rRNA sequences found thus far in cyanobacteria. The divergent sequence is highly conserved among numerous strains of S. hyalinum, suggesting adaptive advantage and selective constraint of the divergent sequence.

Horizontal Gene Transfer and the History of Life

Microbes acquire DNA from a variety of sources. The last decades, which have seen the development of genome sequencing, have revealed that horizontal gene transfer has been a major evolutionary force that has constantly reshaped genomes throughout evolution. However, because the history of life must ultimately be deduced from gene phylogenies, the lack of methods to account for horizontal gene transfer has thrown into confusion the very concept of the tree of life. As a result, many questions remain open, but emerging methodological developments promise to use information conveyed by horizontal gene transfer that remains unexploited today.

Early Microbial Evolution: The Age of Anaerobes

In this article, the term "early microbial evolution" refers to the phase of biological history from the emergence of life to the diversification of the first microbial lineages. In the modern era (since we knew about archaea), three debates have emerged on the subject that deserve discussion: (1) thermophilic origins versus mesophilic origins, (2) autotrophic origins versus heterotrophic origins, and (3) how do eukaryotes figure into early evolution. Here, we revisit those debates from the standpoint of newer data. We also consider the perhaps more pressing issue that molecular phylogenies need to recover anaerobic lineages at the base of prokaryotic trees, because O2 is a product of biological evolution; hence, the first microbes had to be anaerobes. If molecular phylogenies do not recover anaerobes basal, something is wrong. Among the anaerobes, hydrogen-dependent autotrophs--acetogens and methanogens--look like good candidates for the ancestral state of physiology in the bacteria and archaea, respectively. New trees tend to indicate that eukaryote cytosolic ribosomes branch within their archaeal homologs, not as sisters to them and, furthermore tend to root archaea within the methanogens. These are major changes in the tree of life, and open up new avenues of thought. Geochemical methane synthesis occurs as a spontaneous, abiotic exergonic reaction at hydrothermal vents. The overall similarity between that reaction and biological methanogenesis fits well with the concept of a methanogenic root for archaea and an autotrophic origin of microbial physiology.

Detecting Horizontal Gene Transfer between Closely Related Taxa

Horizontal gene transfer (HGT), the transfer of genetic material between organisms, is crucial for genetic innovation and the evolution of genome architecture. Existing HGT detection algorithms rely on a strong phylogenetic signal distinguishing the transferred sequence from ancestral (vertically derived) genes in its recipient genome. Detecting HGT between closely related species or strains is challenging, as the phylogenetic signal is usually weak and the nucleotide composition is normally nearly identical. Nevertheless, there is a great importance in detecting HGT between congeneric species or strains, especially in clinical microbiology, where understanding the emergence of new virulent and drug-resistant strains is crucial, and often time-sensitive.

We developed a novel, self-contained technique named Near HGT, based on the synteny index, to measure the divergence of a gene from its native genomic environment and used it to identify candidate HGT events between closely related strains. The method confirms candidate transferred genes based on the constant relative mutability (CRM). Using CRM, the algorithm assigns a confidence score based on “unusual” sequence divergence. A gene exhibiting exceptional deviations according to both synteny and mutability criteria, is considered a validated HGT product. We first employed the technique to a set of three E. coli strains and detected several highly probable horizontally acquired genes. We then compared the method to existing HGT detection tools using a larger strain data set.

When combined with additional approaches our new algorithm provides richer picture and brings us closer to the goal of detecting all newly acquired genes in a particular strain.

The transfer of genetic material between organisms, usually denoted as horizontal (or lateral) gene transfer (HGT or LGT), is a prime mechanism in microbial evolution and responsible for genetic innovation and the evolution of genome architecture. Detecting HGT between closely related species or strains is imperative as drug-resistant pathogenic strains most often acquire their virulence from closely related bacteria. The proposed method combines two evolutionary signals that were not employed in the past for this task. One is the synteny index (SI), measuring the loss of synteny in an organism, and the other is a novel concept—constant relative mutability (CRM), maintaining that genes preserve their relative evolution rate along linages (although the latter ones may each change).

We show both in simulation and real biological data that the method is sound and, in the cases examined, provides stronger sensitivity than existing methods. We therefore believe this novel approach represents a significant advance, for the first time enabling the detection of previously ignored HGT events that will bring us closer to the goal of detecting all newly acquired genes in a particular strain.

Availability: The method is publicly available at http://research.haifa.ac.il/~ssagi/software/nearHGT.zip

Inferring gene duplications, transfers and losses can be done in a discrete framework

In the field of phylogenetics, the evolutionary history of a set of organisms is commonly depicted by a species tree-whose internal nodes represent speciation events-while the evolutionary history of a gene family is depicted by a gene tree-whose internal nodes can also represent macro-evolutionary events such as gene duplications and transfers. As speciation events are only part of the events shaping a gene history, the topology of a gene tree can show incongruences with that of the corresponding species tree. These incongruences can be used to infer the macro-evolutionary events undergone by the gene family. This is done by embedding the gene tree inside the species tree and hence providing a reconciliation of those trees. In the past decade, several parsimony-based methods have been developed to infer such reconciliations, accounting for gene duplications ([Formula: see text]), transfers ([Formula: see text]) and losses ([Formula: see text]). The main contribution of this paper is to formally prove an important assumption implicitly made by previous works on these reconciliations, namely that solving the (maximum) parsimony [Formula: see text] reconciliation problem in the discrete framework is equivalent to finding a most parsimonious [Formula: see text] scenario in the continuous framework. In the process, we also prove several intermediate results that are useful on their own and constitute a theoretical toolbox that will likely facilitate future theoretical contributions in the field.

Microorganisms-A Journal and a Unifying Concept for the Science of Microbiology

The MDPI journal Microorganisms is still very young, having been launched in 2013, but the concept of the microorganism has been in use for at least a century as a unifying principle for the discipline of microbiology, which was cemented firmly by the intellectual work of Roger Stanier and colleagues in their Microbial World and other general microbiology textbooks and related articles from the 1950s to the 1970s [1,2]. Merging the idea of the microscopic and the very small with the older idea of an organism as a living entity or cell, the concept of a microorganism enabled a real appreciation of the microbial world as one that is amenable to study using similar tools and approaches even though representing distinctly different types of reproductive units and cell organizations. In the late 20th century following the work of Carl Woese and other molecular evolutionists, biologists came to appreciate the commonality among all organisms, all being comprised of cells that bear a remarkable similarity to one another and that share a common evolutionary ancestry, and consequently with major features of a largely shared genetic code and molecular biology. In this sense microbiology and biology as a whole became unified as they never had been before.[...].

Goods-thinking vs. tree-thinking: Finding a place for mobile genetic elements

While it has become increasingly clear that the Tree of Life hypothesis has limitations in its ability to describe the evolution of all evolving entities on the planet, there has been a marked reluctance to move away from the tree-based language. Ironically, while modifying the idea of the Tree of Life to the extent that it is only very distantly related to its original descriptions, there has been a very careful attempt to retain the language of tree-thinking. The recent movement away from a tree-thinking language toward a goods-thinking language and perspective is a significant improvement. In this commentary, we describe how goods-thinking can provide better descriptions of evolution, can integrate evolution with environment more closely and can offer an equal place for Mobile Genetic Elements and chromosomal elements in discussions of evolutionary history.

The extent of genome flux and its role in the differentiation of bacterial lineages

Horizontal gene transfer (HGT) and gene loss are key processes in bacterial evolution. However, the role of gene gain and loss in the emergence and maintenance of ecologically differentiated bacterial populations remains an open question. Here, we use whole-genome sequence data to quantify gene gain and loss for 27 lineages of the plant-associated bacterium Pseudomonas syringae. We apply an extensive error-control procedure that accounts for errors in draft genome data and greatly improves the accuracy of patterns of gene occurrence among these genomes. We demonstrate a history of extensive genome fluctuation for this species and show that individual lineages could have acquired thousands of genes in the same period in which a 1% amino acid divergence accrues in the core genome. Elucidating the dynamics of genome fluctuation reveals the rapid turnover of gained genes, such that the majority of recently gained genes are quickly lost. Despite high observed rates of fluctuation, a phylogeny inferred from patterns of gene occurrence is similar to a phylogeny based on amino acid replacements within the core genome. Furthermore, the core genome phylogeny suggests that P. syringae should be considered a number of distinct species, with levels of divergence at least equivalent to those between recognized bacterial species. Gained genes are transferred from a variety of sources, reflecting the depth and diversity of the potential gene pool available via HGT. Overall, our results provide further insights into the evolutionary dynamics of genome fluctuation and implicate HGT as a major factor contributing to the diversification of P. syringae lineages.

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.

PhyBin: binning trees by topology

A major goal of many evolutionary analyses is to determine the true evolutionary history of an organism. Molecular methods that rely on the phylogenetic signal generated by a few to a handful of loci can be used to approximate the evolution of the entire organism but fall short of providing a global, genome-wide, perspective on evolutionary processes. Indeed, individual genes in a genome may have different evolutionary histories. Therefore, it is informative to analyze the number and kind of phylogenetic topologies found within an orthologous set of genes across a genome. Here we present PhyBin: a flexible program for clustering gene trees based on topological structure. PhyBin can generate bins of topologies corresponding to exactly identical trees or can utilize Robinson-Fould’s distance matrices to generate clusters of similar trees, using a user-defined threshold. Additionally, PhyBin allows the user to adjust for potential noise in the dataset (as may be produced when comparing very closely related organisms) by pre-processing trees to collapse very short branches or those nodes not meeting a defined bootstrap threshold. As a test case, we generated individual trees based on an orthologous gene set from 10 Wolbachia species across four different supergroups (A–D) and utilized PhyBin to categorize the complete set of topologies produced from this dataset. Using this approach, we were able to show that although a single topology generally dominated the analysis, confirming the separation of the supergroups, many genes supported alternative evolutionary histories. Because PhyBin’s output provides the user with lists of gene trees in each topological cluster, it can be used to explore potential reasons for discrepancies between phylogenies including homoplasies, long-branch attraction, or horizontal gene transfer events.

Prokaryotic phylogenies inferred from whole-genome sequence and annotation data

Phylogenetic trees are used to represent the evolutionary relationship among various groups of species. In this paper, a novel method for inferring prokaryotic phylogenies using multiple genomic information is proposed. The method is called CGCPhy and based on the distance matrix of orthologous gene clusters between whole-genome pairs. CGCPhy comprises four main steps. First, orthologous genes are determined by sequence similarity, genomic function, and genomic structure information. Second, genes involving potential HGT events are eliminated, since such genes are considered to be the highly conserved genes across different species and the genes located on fragments with abnormal genome barcode. Third, we calculate the distance of the orthologous gene clusters between each genome pair in terms of the number of orthologous genes in conserved clusters. Finally, the neighbor-joining method is employed to construct phylogenetic trees across different species. CGCPhy has been examined on different datasets from 617 complete single-chromosome prokaryotic genomes and achieved applicative accuracies on different species sets in agreement with Bergey's taxonomy in quartet topologies. Simulation results show that CGCPhy achieves high average accuracy and has a low standard deviation on different datasets, so it has an applicative potential for phylogenetic analysis.

The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses

16S ribosomal RNA currently represents the most important target of study in bacterial ecology. Its use for the description of bacterial diversity is, however, limited by the presence of variable copy numbers in bacterial genomes and sequence variation within closely related taxa or within a genome. Here we use the information from sequenced bacterial genomes to explore the variability of 16S rRNA sequences and copy numbers at various taxonomic levels and apply it to estimate bacterial genome and DNA abundances. In total, 7,081 16S rRNA sequences were in silico extracted from 1,690 available bacterial genomes (1-15 per genome). While there are several phyla containing low 16S rRNA copy numbers, in certain taxa, e.g., the Firmicutes and Gammaproteobacteria, the variation is large. Genome sizes are more conserved at all tested taxonomic levels than 16S rRNA copy numbers. Only a minority of bacterial genomes harbors identical 16S rRNA gene copies, and sequence diversity increases with increasing copy numbers. While certain taxa harbor dissimilar 16S rRNA genes, others contain sequences common to multiple species. Sequence identity clusters (often termed operational taxonomic units) thus provide an imperfect representation of bacterial taxa of a certain phylogenetic rank. We have demonstrated that the information on 16S rRNA copy numbers and genome sizes of genome-sequenced bacteria may be used as an estimate for the closest related taxon in an environmental dataset to calculate alternative estimates of the relative abundance of individual bacterial taxa in environmental samples. Using an example from forest soil, this procedure would increase the abundance estimates of Acidobacteria and decrease these of Firmicutes. Using the currently available information, alternative estimates of bacterial community composition may be obtained in this way if the variation of 16S rRNA copy numbers among bacteria is considered.

Understanding phylogenetic incongruence: lessons from phyllostomid bats

All characters and trait systems in an organism share a common evolutionary history that can be estimated using phylogenetic methods. However, differential rates of change and the evolutionary mechanisms driving those rates result in pervasive phylogenetic conflict. These drivers need to be uncovered because mismatches between evolutionary processes and phylogenetic models can lead to high confidence in incorrect hypotheses. Incongruence between phylogenies derived from morphological versus molecular analyses, and between trees based on different subsets of molecular sequences has become pervasive as datasets have expanded rapidly in both characters and species. For more than a decade, evolutionary relationships among members of the New World bat family Phyllostomidae inferred from morphological and molecular data have been in conflict. Here, we develop and apply methods to minimize systematic biases, uncover the biological mechanisms underlying phylogenetic conflict, and outline data requirements for future phylogenomic and morphological data collection. We introduce new morphological data for phyllostomids and outgroups and expand previous molecular analyses to eliminate methodological sources of phylogenetic conflict such as taxonomic sampling, sparse character sampling, or use of different algorithms to estimate the phylogeny. We also evaluate the impact of biological sources of conflict: saturation in morphological changes and molecular substitutions, and other processes that result in incongruent trees, including convergent morphological and molecular evolution. Methodological sources of incongruence play some role in generating phylogenetic conflict, and are relatively easy to eliminate by matching taxa, collecting more characters, and applying the same algorithms to optimize phylogeny. The evolutionary patterns uncovered are consistent with multiple biological sources of conflict, including saturation in morphological and molecular changes, adaptive morphological convergence among nectar-feeding lineages, and incongruent gene trees. Applying methods to account for nucleotide sequence saturation reduces, but does not completely eliminate, phylogenetic conflict. We ruled out paralogy, lateral gene transfer, and poor taxon sampling and outgroup choices among the processes leading to incongruent gene trees in phyllostomid bats. Uncovering and countering the possible effects of introgression and lineage sorting of ancestral polymorphism on gene trees will require great leaps in genomic and allelic sequencing in this species-rich mammalian family. We also found evidence for adaptive molecular evolution leading to convergence in mitochondrial proteins among nectar-feeding lineages. In conclusion, the biological processes that generate phylogenetic conflict are ubiquitous, and overcoming incongruence requires better models and more data than have been collected even in well-studied organisms such as phyllostomid bats.

Protein based molecular markers provide reliable means to understand prokaryotic phylogeny and support Darwinian mode of evolution

The analyses of genome sequences have led to the proposal that lateral gene transfers (LGTs) among prokaryotes are so widespread that they disguise the interrelationships among these organisms. This has led to questioning of whether the Darwinian model of evolution is applicable to prokaryotic organisms. In this review, we discuss the usefulness of taxon-specific molecular markers such as conserved signature indels (CSIs) and conserved signature proteins (CSPs) for understanding the evolutionary relationships among prokaryotes and to assess the influence of LGTs on prokaryotic evolution. The analyses of genomic sequences have identified large numbers of CSIs and CSPs that are unique properties of different groups of prokaryotes ranging from phylum to genus levels. The species distribution patterns of these molecular signatures strongly support a tree-like vertical inheritance of the genes containing these molecular signatures that is consistent with phylogenetic trees. Recent detailed studies in this regard on the Thermotogae and Archaea, which are reviewed here, have identified large numbers of CSIs and CSPs that are specific for the species from these two taxa and a number of their major clades. The genetic changes responsible for these CSIs (and CSPs) initially likely occurred in the common ancestors of these taxa and then vertically transferred to various descendants. Although some CSIs and CSPs in unrelated groups of prokaryotes were identified, their small numbers and random occurrence has no apparent influence on the consistent tree-like branching pattern emerging from other markers. These results provide evidence that although LGT is an important evolutionary force, it does not mask the tree-like branching pattern of prokaryotes or understanding of their evolutionary relationships. The identified CSIs and CSPs also provide novel and highly specific means for identification of different groups of microbes and for taxonomical and biochemical studies.

Detection of horizontal transfer of individual genes by anomalous oligomer frequencies

Background

Understanding the history of life requires that we understand the transfer of genetic material across phylogenetic boundaries. Detecting genes that were acquired by means other than vertical descent is a basic step in that process. Detection by discordant phylogenies is computationally expensive and not always definitive. Many have used easily computed compositional features as an alternative procedure. However, different compositional methods produce different predictions, and the effectiveness of any method is not well established.

Results

The ability of octamer frequency comparisons to detect genes artificially seeded in cyanobacterial genomes was markedly increased by using as a training set those genes that are highly conserved over all bacteria. Using a subset of octamer frequencies in such tests also increased effectiveness, but this depended on the specific target genome and the source of the contaminating genes. The presence of high frequency octamers and the GC content of the contaminating genes were important considerations. A method comprising best practices from these tests was devised, the Core Gene Similarity (CGS) method, and it performed better than simple octamer frequency analysis, codon bias, or GC contrasts in detecting seeded genes or naturally occurring transposons. From a comparison of predictions with phylogenetic trees, it appears that the effectiveness of the method is confined to horizontal transfer events that have occurred recently in evolutionary time.

Conclusions

The CGS method may be an improvement over existing surrogate methods to detect genes of foreign origin.

Unique core genomes of the bacterial family vibrionaceae: insights into niche adaptation and speciation

Background

The criteria for defining bacterial species and even the concept of bacterial species itself are under debate, and the discussion is apparently intensifying as more genome sequence data is becoming available. However, it is still unclear how the new advances in genomics should be used most efficiently to address this question. In this study we identify genes that are common to any group of genomes in our dataset, to determine whether genes specific to a particular taxon exist and to investigate their potential role in adaptation of bacteria to their specific niche. These genes were named unique core genes. Additionally, we investigate the existence and importance of unique core genes that are found in isolates of phylogenetically non-coherent groups. These groups of isolates, that share a genetic feature without sharing a closest common ancestor, are termed genophyletic groups.

Results

The bacterial family Vibrionaceae was used as the model, and we compiled and compared genome sequences of 64 different isolates. Using the software orthoMCL we determined clusters of homologous genes among the investigated genome sequences. We used multilocus sequence analysis to build a host phylogeny and mapped the numbers of unique core genes of all distinct groups of isolates onto the tree. The results show that unique core genes are more likely to be found in monophyletic groups of isolates. Genophyletic groups of isolates, in contrast, are less common especially for large groups of isolate. The subsequent annotation of unique core genes that are present in genophyletic groups indicate a high degree of horizontally transferred genes. Finally, the annotation of the unique core genes of Vibrio cholerae revealed genes involved in aerotaxis and biosynthesis of the iron-chelator vibriobactin.

Conclusion

The presented work indicates that genes specific for any taxon inside the bacterial family Vibrionaceae exist. These unique core genes encode conserved metabolic functions that can shed light on the adaptation of a species to its ecological niche. Additionally, our study suggests that unique core genes can be used to aid classification of bacteria and contribute to a bacterial species definition on a genomic level. Furthermore, these genes may be of importance in clinical diagnostics and drug development.

Variable copy number, intra-genomic heterogeneities and lateral transfers of the 16S rRNA gene in Pseudomonas

Even though the 16S rRNA gene is the most commonly used taxonomic marker in microbial ecology, its poor resolution is still not fully understood at the intra-genus level. In this work, the number of rRNA gene operons, intra-genomic heterogeneities and lateral transfers were investigated at a fine-scale resolution, throughout the Pseudomonas genus. In addition to nineteen sequenced Pseudomonas strains, we determined the 16S rRNA copy number in four other Pseudomonas strains by Southern hybridization and Pulsed-Field Gel Electrophoresis, and studied the intra-genomic heterogeneities by Denaturing Gradient Gel Electrophoresis and sequencing. Although the variable copy number (from four to seven) seems to be correlated with the evolutionary distance, some close strains in the P. fluorescens lineage showed a different number of 16S rRNA genes, whereas all the strains in the P. aeruginosa lineage displayed the same number of genes (four copies). Further study of the intra-genomic heterogeneities revealed that most of the Pseudomonas strains (15 out of 19 strains) had at least two different 16S rRNA alleles. A great difference (5 or 19 nucleotides, essentially grouped near the V1 hypervariable region) was observed only in two sequenced strains. In one of our strains studied (MFY30 strain), we found a difference of 12 nucleotides (grouped in the V3 hypervariable region) between copies of the 16S rRNA gene. Finally, occurrence of partial lateral transfers of the 16S rRNA gene was further investigated in 1803 full-length sequences of Pseudomonas available in the databases. Remarkably, we found that the two most variable regions (the V1 and V3 hypervariable regions) had probably been laterally transferred from another evolutionary distant Pseudomonas strain for at least 48.3 and 41.6% of the 16S rRNA sequences, respectively. In conclusion, we strongly recommend removing these regions of the 16S rRNA gene during the intra-genus diversity studies.

High expression hampers horizontal gene transfer

Horizontal gene transfer (HGT), the movement of genetic material from one species to another, is a common phenomenon in prokaryotic evolution. Although the rate of HGT is known to vary among genes, our understanding of the cause of this variation, currently summarized by two rules, is far from complete. The first rule states that informational genes, which are involved in DNA replication, transcription, and translation, have lower transferabilities than operational genes. The second rule asserts that protein interactivity negatively impacts gene transferability. Here, we hypothesize that high expression hampers HGT, because the fitness cost of an HGT to the recipient, arising from the 1) energy expenditure in transcription and translation, 2) cytotoxic protein misfolding, 3) reduction in cellular translational efficiency, 4) detrimental protein misinteraction, and 5) disturbance of the optimal protein concentration or cell physiology, increases with the expression level of the transferred gene. To test this hypothesis, we examined laboratory and natural HGTs to Escherichia coli. We observed lower transferabilities of more highly expressed genes, even after controlling the confounding factors from the two established rules and the genic GC content. Furthermore, expression level predicts gene transferability better than all other factors examined. We also confirmed the significant negative impact of gene expression on the rate of HGTs to 127 of 133 genomes of eubacteria and archaebacteria. Together, these findings establish the gene expression level as a major determinant of horizontal gene transferability. They also suggest that most successful HGTs are initially slightly deleterious, fixed because of their negligibly low costs rather than high benefits to the recipient.

Lateral gene transfer as a support for the tree of life

Lateral gene transfer (LGT), the acquisition of genes from other species, is a major evolutionary force. However, its success as an adaptive process makes the reconstruction of the history of life an intricate puzzle: If no gene has remained unaffected during the course of life's evolution, how can one rely on molecular markers to reconstruct the relationships among species? Here, we take a completely different look at LGT and its impact for the reconstruction of the history of life. Rather than trying to remove the effect of LGT in phylogenies, and ignoring as a result most of the information of gene histories, we use an explicit phylogenetic model of gene transfer to reconcile gene histories with the tree of species. We studied 16 bacterial and archaeal phyla, representing a dataset of 12,000 gene families distributed in 336 genomes. Our results show that, in most phyla, LGT provides an abundant phylogenetic signal on the pattern of species diversification and that this signal is robust to the choice of gene families under study. We also find that LGT brings an abundant signal on the location of the root of species trees, which has been previously overlooked. Our results quantify the great variety of gene transfer rates among lineages of the tree of life and provide strong support for the "complexity hypothesis," which states that genes whose products participate to macromolecular protein complexes are relatively resistant to transfer.

Phylogenetic analysis supports horizontal gene transfer of L-amino acid oxidase gene in Streptococcus oligofermentans

Phylogenetic analysis of 10 amino acid sequences from 19 Streptococcus species showed that S. oligofermentans clustered within the mitis group. However, the l-amino acid oxidase (LAAO) of S. oligofermentans showed a different clustering pattern from the other proteins analyzed implicating horizontal gene transfer (HGT) in the origin of the S. oligofermentans LAAO gene. LAAO of S. oligofermentans is known to confer ability to compete with other oral cavity bacteria, most notably S. mutans; therefore, the HGT event may have been important in extending the ecological niche occupied by this species, consistent with those of other studies suggesting that HGT can play a key role in enabling bacterial species to occupy new ecological niches.

Analysis of the global ocean sampling (GOS) project for trends in iron uptake by surface ocean microbes

Microbial metagenomes are DNA samples of the most abundant, and therefore most successful organisms at the sampling time and location for a given cell size range. The study of microbial communities via their DNA content has revolutionized our understanding of microbial ecology and evolution. Iron availability is a critical resource that limits microbial communities' growth in many oceanic areas. Here, we built a database of 2319 sequences, corresponding to 140 gene families of iron metabolism with a large phylogenetic spread, to explore the microbial strategies of iron acquisition in the ocean's bacterial community. We estimate iron metabolism strategies from metagenome gene content and investigate whether their prevalence varies with dissolved iron concentrations obtained from a biogeochemical model. We show significant quantitative and qualitative variations in iron metabolism pathways, with a higher proportion of iron metabolism genes in low iron environments. We found a striking difference between coastal and open ocean sites regarding Fe²⁺ versus Fe³⁺ uptake gene prevalence. We also show that non-specific siderophore uptake increases in low iron open ocean environments, suggesting bacteria may acquire iron from natural siderophore-like organic complexes. Despite the lack of knowledge of iron uptake mechanisms in most marine microorganisms, our approach provides insights into how the iron metabolic pathways of microbial communities may vary with seawater iron concentrations.