Genomes in flux: the evolution of archaeal and proteobacterial gene content

Corynebacterium pseudotuberculosis: Whole genome sequencing reveals unforeseen and relevant genetic diversity in this pathogen

Corynebacterium pseudotuberculosis (CPS) is an important bacterial animal pathogen. CPS causes chronic, debilitating and currently incurable infectious diseases affecting a wide range of livestock and wild herbivores including camelids worldwide. Belonging to the Corynebacterium diphtheriae complex, this pathogen can also infect humans. The classical characterization of CPS is typically based on the testing of nitrate reductase activity, separating the two biovars Equi and Ovis. However, more refined resolutions are required to unravel routes of infection. This was realized in our study by generating and analyzing whole genome sequencing (WGS) data. Using newly created core genome multilocus sequence typing (cgMLST) profiles we were the first to discover isolates grouping in a cluster adjacent to clusters formed by CPS biovar Equi isolates. This novel cluster includes CPS isolates from alpacas, llamas, camels and dromedaries, which are characterized by a lack of nitrate reductase activity as encountered in biovar Ovis. This is of special interest for molecular epidemiology. Nevertheless, these isolates bear the genes of the nitrate locus, which are characteristic of biovar Equi isolates. However, sequence analysis of the genes narG and narH of the nitrate locus revealed indels leading to frameshifts and inactivity of the enzymes involved in nitrate reduction. Interestingly, one CPS isolate originating from another lama with an insertion in the MFS transporter (narT) is adjacent to a cluster formed by ovine CPS isolates biovar Equi. Based on this knowledge, the combination of biochemical and PCR based molecular biological nitrate reductase detection can be used for a fast and uncomplicated classification of isolates in routine diagnostics in order to check the origin of camelid CPS isolates. Further analysis revealed that partial sequencing of the ABC transporter substrate binding protein (CP258_RS07935) is a powerful tool to assign the biovars and the novel genomovar.

RNA-guided RNA silencing by an Asgard archaeal Argonaute

Argonaute proteins are the central effectors of RNA-guided RNA silencing pathways in eukaryotes, playing crucial roles in gene repression and defense against viruses and transposons. Eukaryotic Argonautes are subdivided into two clades: AGOs generally facilitate miRNA- or siRNA-mediated silencing, while PIWIs generally facilitate piRNA-mediated silencing. It is currently unclear when and how Argonaute-based RNA silencing mechanisms arose and diverged during the emergence and early evolution of eukaryotes. Here, we show that in Asgard archaea, the closest prokaryotic relatives of eukaryotes, an evolutionary expansion of Argonaute proteins took place. In particular, a deep-branching PIWI protein (HrAgo1) encoded by the genome of the Lokiarchaeon 'Candidatus Harpocratesius repetitus' shares a common origin with eukaryotic PIWI proteins. Contrasting known prokaryotic Argonautes that use single-stranded DNA as guides and/or targets, HrAgo1 mediates RNA-guided RNA cleavage, and facilitates gene silencing when expressed in human cells and supplied with miRNA precursors. A cryo-EM structure of HrAgo1, combined with quantitative single-molecule experiments, reveals that the protein displays structural features and target-binding modes that are a mix of those of eukaryotic AGO and PIWI proteins. Thus, this deep-branching archaeal PIWI may have retained an ancestral molecular architecture that preceded the functional and mechanistic divergence of eukaryotic AGOs and PIWIs.

Identifying and tracking mobile elements in evolving compost communities yields insights into the nanobiome

Microbial evolution is driven by rapid changes in gene content mediated by horizontal gene transfer (HGT). While mobile genetic elements (MGEs) are important drivers of gene flux, the nanobiome-the zoo of Darwinian replicators that depend on microbial hosts-remains poorly characterised. New approaches are necessary to increase our understanding beyond MGEs shaping individual populations, towards their impacts on complex microbial communities. A bioinformatic pipeline (xenoseq) was developed to cross-compare metagenomic samples from microbial consortia evolving in parallel, aimed at identifying MGE dissemination, which was applied to compost communities which underwent periodic mixing of MGEs. We show that xenoseq can distinguish movement of MGEs from demographic changes in community composition that otherwise confounds identification, and furthermore demonstrate the discovery of various unexpected entities. Of particular interest was a nanobacterium of the candidate phylum radiation (CPR) which is closely related to a species identified in groundwater ecosystems (Candidatus Saccharibacterium), and appears to have a parasitic lifestyle. We also highlight another prolific mobile element, a 313 kb plasmid hosted by a Cellvibrio lineage. The host was predicted to be capable of nitrogen fixation, and acquisition of the plasmid coincides with increased ammonia production. Taken together, our data show that new experimental strategies combined with bioinformatic analyses of metagenomic data stand to provide insight into the nanobiome as a driver of microbial community evolution.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

In vitro functional characterization predicts the impact of bacterial root endophytes on plant growth

Utilizing beneficial microbes for crop improvement is one strategy to achieve sustainable agriculture. However, identifying microbial isolates that promote crop growth is challenging, in part because using bacterial taxonomy to predict an isolate's effect on plant growth may not be reliable. The overall aim of this work was to determine whether in vitro functional traits of bacteria were predictive of their in planta impact. We isolated 183 bacterial endophytes from field-grown roots of two tomato species, Solanum lycopersicum and S. pimpinellifolium. Sixty isolates were screened for six in vitro functional traits: auxin production, siderophore production, phosphate solubilization, antagonism to a soilborne pathogen, and the presence of two antimicrobial metabolite synthesis genes. Hierarchical clustering of the isolates based on the in vitro functional traits identified several groups of isolates sharing similar traits. We called these groups 'functional groups'. To understand how in vitro functional traits of bacteria relate to their impact on plants, we inoculated three isolates from each of the functional groups on tomato seedlings. Isolates within the same functional group promoted plant growth at similar levels, regardless of their host origin or taxonomy. Together, our results demonstrate the importance of examining root endophyte functions for improving crop production.

Retrieving a disrupted gene encoding phospholipase A for fibre enhancement in allotetraploid cultivated cotton

After polyploidization originated from one interspecific hybridization event in Gossypium, Gossypium barbadense evolved to produce extra-long staple fibres than Gossypium hirsutum (Upland cotton), which produces a higher fibre yield. The genomic diversity between G. barbadense and G. hirsutum thus provides a genetic basis for fibre trait variation. Recently, rapid accumulation of gene disruption or deleterious mutation was reported in allotetraploid cotton genomes, with unknown impacts on fibre traits. Here, we identified gene disruptions in allotetraploid G. hirsutum (18.14%) and G. barbadense (17.38%) through comparison with their presumed diploid progenitors. Relative to conserved genes, these disrupted genes exhibited faster evolution rate, lower expression level and altered gene co-expression networks. Within a module regulating fibre elongation, a hub gene experienced gene disruption in G. hirsutum after polyploidization, with a 2-bp deletion in the coding region of GhNPLA1D introducing early termination of translation. This deletion was observed in all of the 34 G. hirsutum landraces and 36 G. hirsutum cultivars, but not in 96% of 57 G. barbadense accessions. Retrieving the disrupted gene GhNPLA1D using its homoeolog GhNPLA1A achieved longer fibre length in G. hirsutum. Further enzyme activity and lipids analysis confirmed that GhNPLA1A encodes a typical phospholipase A and promotes cotton fibre elongation via elevating intracellular levels of linolenic acid and 34:3 phosphatidylinositol. Our work opens a strategy for identifying disrupted genes and retrieving their functions in ways that can provide valuable resources for accelerating fibre trait enhancement in cotton breeding.

Nutrition or nature: using elementary flux modes to disentangle the complex forces shaping prokaryote pan-genomes

BACKGROUND

Microbial pan-genomes are shaped by a complex combination of stochastic and deterministic forces. Even closely related genomes exhibit extensive variation in their gene content. Understanding what drives this variation requires exploring the interactions of gene products with each other and with the organism's external environment. However, to date, conceptual models of pan-genome dynamics often represent genes as independent units and provide limited information about their mechanistic interactions.

RESULTS

We simulated the stochastic process of gene-loss using the pooled genome-scale metabolic reaction networks of 46 taxonomically diverse bacterial and archaeal families as proxies for their pan-genomes. The frequency by which reactions are retained in functional networks when stochastic gene loss is simulated in diverse environments allowed us to disentangle the metabolic reactions whose presence depends on the metabolite composition of the external environment (constrained by "nutrition") from those that are independent of the environment (constrained by "nature"). By comparing the frequency of reactions from the first group with their observed frequencies in bacterial and archaeal families, we predicted the metabolic niches that shaped the genomic composition of these lineages. Moreover, we found that the lineages that were shaped by a more diverse metabolic niche also occur in more diverse biomes as assessed by global environmental sequencing datasets.

CONCLUSION

We introduce a computational framework for analyzing and interpreting pan-reactomes that provides novel insights into the ecological and evolutionary drivers of pan-genome dynamics.

Realistic Gene Transfer to Gene Duplication Ratios Identify Different Roots in the Bacterial Phylogeny Using a Tree Reconciliation Method

The rooting of phylogenetic trees permits important inferences about ancestral states and the polarity of evolutionary events. Recently, methods that reconcile discordance between gene-trees and species-trees—tree reconciliation methods—are becoming increasingly popular for rooting species trees. Rooting via reconciliation requires values for a particular parameter, the gene transfer to gene duplication ratio (T:D), which in current practice is estimated on the fly from discordances observed in the trees. To date, the accuracy of T:D estimates obtained by reconciliation analyses has not been compared to T:D estimates obtained by independent means, hence the effect of T:D upon inferences of species tree roots is altogether unexplored. Here we investigated the issue in detail by performing tree reconciliations of more than 10,000 gene trees under a variety of T:D ratios for two phylogenetic cases: a bacterial (prokaryotic) tree with 265 species and a fungal-metazoan (eukaryotic) tree with 31 species. We show that the T:D ratios automatically estimated by a current tree reconciliation method, ALE, generate virtually identical T:D ratios across bacterial genes and fungal-metazoan genes. The T:D ratios estimated by ALE differ 10- to 100-fold from robust, ALE-independent estimates from real data. More important is our finding that the root inferences using ALE in both datasets are strongly dependent upon T:D. Using more realistic T:D ratios, the number of roots inferred by ALE consistently increases and, in some cases, clearly incorrect roots are inferred. Furthermore, our analyses reveal that gene duplications have a far greater impact on ALE’s preferences for phylogenetic root placement than gene transfers or gene losses do. Overall, we show that obtaining reliable species tree roots with ALE is only possible when gene duplications are abundant in the data and the number of falsely inferred gene duplications is low. Finding a sufficient sample of true gene duplications for rooting species trees critically depends on the T:D ratios used in the analyses. T:D ratios, while being important parameters of genome evolution in their own right, affect the root inferences with tree reconciliations to an unanticipated degree.

SEQUENCE-FREE PHYLOGENETICS WITH MASS SPECTROMETRY

An alternative, more rapid, sequence-free approach to build phylogenetic trees has been conceived and implemented. Molecular phylogenetics has continued to mostly focus on improvement in tree construction based on gene sequence alignments. Here protein-based phylogenies are constructed using numerical data sets ("phylonumerics") representing the masses of peptide segments recorded in a mass mapping experiment. This truly sequence-free method requires no gene sequences, nor their alignment, to build the trees affording a considerable time and cost-saving to conventional phylogenetics methods. The approach also calculates single point amino acid mutations from a comparison of mass pairs from different maps in the data set and displays these at branch nodes across the tree together with their frequency. Studies of the consecutive, and near-consecutive, ancestral and descendant mutations across interconnected branches of a mass tree allow putative adaptive, epistatic, and compensatory mutations to be identified in order to investigate mechanisms associated with evolutionary processes and pathways. A side-by-side comparison of this sequence-free approach and conventional gene sequence phylogenetics is discussed.

Genetic and phylogenetic analysis of dissimilatory iodate-reducing bacteria identifies potential niches across the world's oceans

Iodine is oxidized and reduced as part of a biogeochemical cycle that is especially pronounced in the oceans, where the element naturally concentrates. The use of oxidized iodine in the form of iodate (IO3-) as an electron acceptor by microorganisms is poorly understood. Here, we outline genetic, physiological, and ecological models for dissimilatory IO3- reduction to iodide (I-) by a novel estuarine bacterium, Denitromonas sp. IR-12. Our results show that dissimilatory iodate reduction (DIR) by strain IR-12 is molybdenum-dependent and requires an IO3- reductase (idrA) and likely other genes in a mobile cluster with a conserved association across known and predicted DIR microorganisms (DIRM). Based on genetic and physiological data, we propose a model where three molecules of IO3- are likely reduced to three molecules of hypoiodous acid (HIO), which rapidly disproportionate into one molecule of IO3- and two molecules of iodide (I-), in a respiratory pathway that provides an energy yield equivalent to that of nitrate or perchlorate respiration. Consistent with the ecological niche expected of such a metabolism, idrA is enriched in the metagenome sequence databases of marine sites with a specific biogeochemical signature (high concentrations of nitrate and phosphate) and diminished oxygen. Taken together, these data suggest that DIRM help explain the disequilibrium of the IO3-:I- concentration ratio above oxygen-minimum zones and support a widespread iodine redox cycle mediated by microbiology.

Gene Duplications Are At Least 50 Times Less Frequent than Gene Transfers in Prokaryotic Genomes

The contribution of gene duplications to the evolution of eukaryotic genomes is well studied. By contrast, studies of gene duplications in prokaryotes are scarce and generally limited to a handful of genes or careful analysis of a few prokaryotic lineages. Systematic broad-scale studies of prokaryotic genomes that sample available data are lacking, leaving gaps in our understanding of the contribution of gene duplications as a source of genetic novelty in the prokaryotic world. Here, we report conservative and robust estimates for the frequency of recent gene duplications within prokaryotic genomes relative to recent lateral gene transfer (LGT), as mechanisms to generate multiple copies of related sequences in the same genome. We obtain our estimates by focusing on evolutionarily recent events among 5,655 prokaryotic genomes, thereby avoiding vagaries of deep phylogenetic inference and confounding effects of ancient events and differential loss. We find that recent, genome-specific gene duplications are at least 50 times less frequent and probably 100 times less frequent than recent, genome-specific, gene acquisitions via LGT. The frequency of gene duplications varies across lineages and functional categories. The findings improve our understanding of genome evolution in prokaryotes and have far-reaching implications for evolutionary models that entail LGT to gene duplications ratio as a parameter.

The metabolic footprint of Clostridia and Erysipelotrichia reveals their role in depleting sugar alcohols in the cecum

BACKGROUND

The catabolic activity of the microbiota contributes to health by aiding in nutrition, immune education, and niche protection against pathogens. However, the nutrients consumed by common taxa within the gut microbiota remain incompletely understood.

METHODS

Here we combined microbiota profiling with an un-targeted metabolomics approach to determine whether depletion of small metabolites in the cecum of mice correlated with the presence of specific bacterial taxa. Causality was investigated by engrafting germ-free or antibiotic-treated mice with complex or defined microbial communities.

RESULTS

We noted that a depletion of Clostridia and Erysipelotrichia from the gut microbiota triggered by antibiotic treatment was associated with an increase in the cecal concentration of sugar acids and sugar alcohols (polyols). Notably, when we inoculated germ-free mice with a defined microbial community of 14 Clostridia and 3 Erysipelotrichia isolates, we observed the inverse, with a marked decrease in the concentrations of sugar acids and polyols in cecal contents. The carbohydrate footprint produced by the defined microbial community was similar to that observed in gnotobiotic mice receiving a cecal microbiota transplant from conventional mice. Supplementation with sorbitol, a polyol used as artificial sweetener, increased cecal sorbitol concentrations in antibiotic-treated mice, which was abrogated after inoculation with a Clostridia isolate able to grow on sorbitol in vitro.

CONCLUSIONS

We conclude that consumption of sugar alcohols by Clostridia and Erysipelotrichia species depletes these metabolites from the intestinal lumen during homeostasis. Video abstract.

Mirage: estimation of ancestral gene-copy numbers by considering different evolutionary patterns among gene families

Motivation

Reconstruction of gene copy number evolution is an essential approach for understanding how complex biological systems have been organized. Although various models have been proposed for gene copy number evolution, existing evolutionary models have not appropriately addressed the fact that different gene families can have very different gene gain/loss rates.

Results

In this study, we developed Mirage (MIxtuRe model for Ancestral Genome Estimation), which allows different gene families to have flexible gene gain/loss rates. Mirage can use three models for formulating heterogeneous evolution among gene families: the discretized Γ model, probability distribution-free model and pattern mixture (PM) model. Simulation analysis showed that Mirage can accurately estimate heterogeneous gene gain/loss rates and reconstruct gene-content evolutionary history. Application to empirical datasets demonstrated that the PM model fits genome data from various taxonomic groups better than the other heterogeneous models. Using Mirage, we revealed that metabolic function-related gene families displayed frequent gene gains and losses in all taxa investigated.

Availability and implementation

The source code of Mirage is freely available at https://github.com/fukunagatsu/Mirage.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Evolution of Microbial Genomics: Conceptual Shifts over a Quarter Century

Prokaryote genomics started in earnest in 1995, with the complete sequences of two small bacterial genomes, those of Haemophilus influenzae and Mycoplasma genitalium. During the next quarter century, the prokaryote genome database has been growing exponentially, with no saturation in sight. For most of these 25 years, genome sequencing remained limited to cultivable microbes. Together with next-generation sequencing methods, advances in metagenomics and single-cell genomics have lifted this limitation, providing for an increasingly unbiased characterization of the global prokaryote diversity. Advances in computational genomics followed the progress of genome sequencing, even if occasionally lagging behind. Several major new branches of bacteria and archaea were discovered, including Asgard archaea, the apparent closest relatives of eukaryotes and expansive groups of bacteria and archaea with small genomes thought to be symbionts of other prokaryotes. Comparative analysis of numerous prokaryote genomes spanning a wide range of evolutionary distances changed the conceptual foundations of microbiology, supplanting the notion of species genomes with fixed gene sets with that of dynamic pangenomes and the notion of a single Tree of Life (ToL) with a statistical tree-like trend among individual gene trees. Strides were also made towards a theory and quantitative laws of prokaryote genome evolution.

Evolution, Expression Profile, Regulatory Mechanism, and Functional Verification of EBP-Like Gene in Cholesterol Biosynthetic Process in Chickens (Gallus Gallus)

The emopamil binding protein (EBP) is an important enzyme participating in the final steps of cholesterol biosynthesis in mammals. A predictive gene EBP-like, which encodes the protein with a high identity to human EBP, was found in chicken genome. No regulatory mechanisms and biological functions of EBP-like have been characterized in chickens. In the present study, the coding sequence of EBP-like was cloned, the phylogenetic trees of EBP/EBP-like were constructed and the genomic synteny of EBP-like was analyzed. The regulatory mechanism of EBP-like were explored with in vivo and in vitro experiments. The biological functions of EBP-like in liver cholesterol biosynthetic were examined by using gain- or loss-of-function strategies. The results showed that chicken EBP-like gene was originated from a common ancestral with Japanese quail EBP gene, and was relatively conservative with EBP gene among different species. The EBP-like gene was highly expressed in liver, its expression level was significantly increased in peak-laying stage, and was upregulated by estrogen. Inhibition of the EBP-like mRNA expression could restrain the expressions of EBP-like downstream genes (SC5D, DHCR24, and DHCR7) in the cholesterol synthetic pathway, therefore downregulate the liver intracellular T-CHO level. In conclusion, as substitute of EBP gene in chickens, EBP-like plays a vital role in the process of chicken liver cholesterol synthesis. This research provides a basis for revealing the molecular regulatory mechanism of cholesterol synthesis in birds, contributes insights into the improvement of the growth and development, laying performance and egg quality in poultry.

Bacterial dispersal and drift drive microbiome diversity patterns within a population of feral hindgut fermenters

Studies of microbiome variation in wildlife often emphasize host physiology and diet as proximate selective pressures acting on host-associated microbiota. In contrast, microbial dispersal and ecological drift are more rarely considered. Using amplicon sequencing, we characterized the bacterial microbiome of adult female (n = 86) Sable Island horses (Nova Scotia, Canada) as part of a detailed individual-based study of this feral population. Using data on sampling date, horse location, age, parental status, and local habitat variables, we contrasted the ability of spatiotemporal, life history, and environmental factors to explain microbiome diversity among Sable Island horses. We extended inferences made from these analyses with both phylogeny-informed and phylogeny-independent null modelling approaches to identify deviations from stochastic expectations. Phylogeny-informed diversity measures were correlated with spatial and local habitat variables, but null modelling results suggested that heterogeneity in ecological drift, rather than differential selective pressures acting on the microbiome, was responsible for these correlations. Conversely, phylogeny-independent diversity measures were best explained by host spatial and social structure, suggesting that taxonomic composition of the microbiome was shaped most strongly by bacterial dispersal. Parental status was important but correlated with measures of β-dispersion rather than β-diversity (mares without foals had lower alpha diversity and more variable microbiomes than mares with foals). Our results suggest that between host microbiome variation within the Sable Island horse population is driven more strongly by bacterial dispersal and ecological drift than by differential selective pressures. These results emphasize the need to consider alternative ecological processes in the study of microbiomes.

Slightly beneficial genes are retained by bacteria evolving DNA uptake despite selfish elements

Horizontal gene transfer (HGT) and gene loss result in rapid changes in the gene content of bacteria. While HGT aids bacteria to adapt to new environments, it also carries risks such as selfish genetic elements (SGEs). Here, we use modelling to study how HGT of slightly beneficial genes impacts growth rates of bacterial populations, and if bacterial collectives can evolve to take up DNA despite selfish elements. We find four classes of slightly beneficial genes: indispensable, enrichable, rescuable, and unrescuable genes. Rescuable genes - genes with small fitness benefits that are lost from the population without HGT - can be collectively retained by a community that engages in costly HGT. While this 'gene-sharing' cannot evolve in well-mixed cultures, it does evolve in a spatial population like a biofilm. Despite enabling infection by harmful SGEs, the uptake of foreign DNA is evolutionarily maintained by the hosts, explaining the coexistence of bacteria and SGEs.

The complex phylogenetic relationships of a 4mC/6mA DNA methyltransferase in prokaryotes

DNA methyltransferases are proteins that modify DNA via attachment of methyl groups to nucleobases and are ubiquitous across the bacterial, archaeal, and eukaryotic domains of life. Here, we investigated the complex evolutionary history of the large and consequential 4mC/6mA DNA methyltransferase protein family using phylogenetic reconstruction of amino acid sequences. We present a well-supported phylogeny of this family based on systematic sampling of taxa across superphyla of bacteria and archaea. We compared the phylogeny to a current representation of the species tree of life and found that the 4mC/6mA methyltransferase family has a strikingly complex evolutionary history that likely began sometime after the last universal common ancestor of life diverged into the bacterial and archaeal lineages and probably involved many horizontal gene transfers within and between domains. Despite the complexity of its evolutionary history, we inferred that only one significant shift in molecular evolutionary rate characterizes the diversification of this protein family.

Phylogenetic Analysis of HIV-1 Genomes Based on the Position-Weighted K-mers Method

HIV-1 viruses, which are predominant in the family of HIV viruses, have strong pathogenicity and infectivity. They can evolve into many different variants in a very short time. In this study, we propose a new and effective alignment-free method for the phylogenetic analysis of HIV-1 viruses using complete genome sequences. Our method combines the position distribution information and the counts of the k-mers together. We also propose a metric to determine the optimal k value. We name our method the Position-Weighted k-mers (PWkmer) method. Validation and comparison with the Robinson-Foulds distance method and the modified bootstrap method on a benchmark dataset show that our method is reliable for the phylogenetic analysis of HIV-1 viruses. PWkmer can resolve within-group variations for different known subtypes of Group M of HIV-1 viruses. This method is simple and computationally fast for whole genome phylogenetic analysis.

Bipartite Network Analysis of Gene Sharings in the Microbial World

Extensive microbial gene flows affect how we understand virology, microbiology, medical sciences, genetic modification, and evolutionary biology. Phylogenies only provide a narrow view of these gene flows: plasmids and viruses, lacking core genes, cannot be attached to cellular life on phylogenetic trees. Yet viruses and plasmids have a major impact on cellular evolution, affecting both the gene content and the dynamics of microbial communities. Using bipartite graphs that connect up to 149,000 clusters of homologous genes with 8,217 related and unrelated genomes, we can in particular show patterns of gene sharing that do not map neatly with the organismal phylogeny. Homologous genes are recycled by lateral gene transfer, and multiple copies of homologous genes are carried by otherwise completely unrelated (and possibly nested) genomes, that is, viruses, plasmids and prokaryotes. When a homologous gene is present on at least one plasmid or virus and at least one chromosome, a process of "gene externalization," affected by a postprocessed selected functional bias, takes place, especially in Bacteria. Bipartite graphs give us a view of vertical and horizontal gene flow beyond classic taxonomy on a single very large, analytically tractable, graph that goes beyond the cellular Web of Life.

Genome rearrangements in Escherichia coli during de novo acquisition of resistance to a single antibiotic or two antibiotics successively

Background

The ability of bacteria to acquire resistance to antibiotics relies to a large extent on their capacity for genome modification. Prokaryotic genomes are highly plastic and can utilize horizontal gene transfer, point mutations, and gene deletions or amplifications to realize genome expansion and rearrangements. The contribution of point mutations to de novo acquisition of antibiotic resistance is well-established. In this study, the internal genome rearrangement of Escherichia coli during to de novo acquisition of antibiotic resistance was investigated using whole-genome sequencing.

Results

Cells were made resistant to one of the four antibiotics and subsequently to one of the three remaining. This way the initial genetic rearrangements could be documented together with the effects of an altered genetic background on subsequent development of resistance. A DNA fragment including ampC was amplified by a factor sometimes exceeding 100 as a result of exposure to amoxicillin. Excision of prophage e14 was observed in many samples with a double exposure history, but not in cells exposed to a single antibiotic, indicating that the activation of the SOS stress response alone, normally the trigger for excision, was not sufficient to cause excision of prophage e14. Partial deletion of clpS and clpA occurred in strains exposed to enrofloxacin and tetracycline. Other deletions were observed in some strains, but not in replicates with the exact same exposure history. Various insertion sequence transpositions correlated with exposure to specific antibiotics.

Conclusions

Many of the genome rearrangements have not been reported before to occur during resistance development. The observed correlation between genome rearrangements and specific antibiotic pressure, as well as their presence in independent replicates indicates that these events do not occur randomly. Taken together, the observed genome rearrangements illustrate the plasticity of the E. coli genome when exposed to antibiotic stress.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5353-y) contains supplementary material, which is available to authorized users.

Multilevel comparative bioinformatics to investigate evolutionary relationships and specificities in gene annotations: an example for tomato and grapevine

Background

“Omics” approaches may provide useful information for a deeper understanding of speciation events, diversification and function innovation. This can be achieved by investigating the molecular similarities at sequence level between species, allowing the definition of ortholog and paralog genes. However, the spreading of sequenced genome, often endowed with still preliminary annotations, requires suitable bioinformatics to be appropriately exploited in this framework.

Results

We presented here a multilevel comparative approach to investigate on genome evolutionary relationships and peculiarities of two fleshy fruit species of relevant agronomic interest, Solanum lycopersicum (tomato) and Vitis vinifera (grapevine). We defined 17,823 orthology relationships between tomato and grapevine reference gene annotations. The resulting orthologs are associated with the detected paralogs in each species, permitting the definition of gene networks, useful to investigate the different relationships. The reconciliation of the compared collections in terms of an updating of the functional descriptions was also exploited. All the results were made accessible in ComParaLogs, a dedicated bioinformatics platform available at http://biosrv.cab.unina.it/comparalogs/gene/search.

Conclusions

The aim of the work was to suggest a reliable approach to detect all similarities of gene loci between two species based on the integration of results from different levels of information, such as the gene, the transcript and the protein sequences, overcoming possible limits due to exclusive protein versus protein comparisons. This to define reliable ortholog and paralog genes, as well as species specific gene loci in the two species, overcoming limits due to the possible draft nature of preliminary gene annotations. Moreover, reconciled functional descriptions, as well as common or peculiar enzymatic classes and protein domains from tomato and grapevine, together with the definition of species-specific gene sets after the pairwise comparisons, contributed a comprehensive set of information useful to comparatively exploit the two species gene annotations and investigate on differences between species with climacteric and non-climacteric fruits. In addition, the definition of networks of ortholog genes and of associated paralogs, and the organization of web-based interfaces for the exploration of the results, defined a friendly computational bench-work in support of comparative analyses between two species.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2420-y) contains supplementary material, which is available to authorized users.

Pleiotropy, cooperation, and the social evolution of genetic architecture

Pleiotropy has been suggested as a novel mechanism for stabilising cooperation in bacteria and other microbes. The hypothesis is that linking cooperation with a trait that provides a personal (private) benefit can outweigh the cost of cooperation in situations when cooperation would not be favoured by mechanisms such as kin selection. We analysed the theoretical plausibility of this hypothesis, with analytical models and individual-based simulations. We found that (1) pleiotropy does not stabilise cooperation, unless the cooperative and private traits are linked via a genetic architecture that cannot evolve (mutational constraint); (2) if the genetic architecture is constrained in this way, then pleiotropy favours any type of trait and not especially cooperation; (3) if the genetic architecture can evolve, then pleiotropy does not favour cooperation; and (4) there are several alternative explanations for why traits may be linked, and causality can even be predicted in the opposite direction, with cooperation favouring pleiotropy. Our results suggest that pleiotropy could only explain cooperation under restrictive conditions and instead show how social evolution can shape the genetic architecture.

Four domains: The fundamental unicell and Post-Darwinian Cognition-Based Evolution

Contemporary research supports the viewpoint that self-referential cognition is the proper definition of life. From that initiating platform, a cohesive alternative evolutionary narrative distinct from standard Neodarwinism can be presented. Cognition-Based Evolution contends that biological variation is a product of a self-reinforcing information cycle that derives from self-referential attachment to biological information space-time with its attendant ambiguities. That information cycle is embodied through obligatory linkages among energy, biological information, and communication. Successive reiterations of the information cycle enact the informational architectures of the basic unicellular forms. From that base, inter-domain and cell-cell communications enable genetic and cellular variations through self-referential natural informational engineering and cellular niche construction. Holobionts are the exclusive endpoints of that self-referential cellular engineering as obligatory multicellular combinations of the essential Four Domains: Prokaryota, Archaea, Eukaryota and the Virome. Therefore, it is advocated that these Four Domains represent the perpetual object of the living circumstance rather than the visible macroorganic forms. In consequence, biology and its evolutionary development can be appraised as the continual defense of instantiated cellular self-reference. As the survival of cells is as dependent upon limitations and boundaries as upon any freedom of action, it is proposed that selection represents only one of many forms of cellular constraint that sustain self-referential integrity.

Phenotypic and genotypic characterisation of an unique indigenous hypersaline unicellular cyanobacterium, Euhalothece sp.nov

A novel halotolerant species of cyanobacterium of the order Chroococcales was isolated from hypersaline estuary in Kwa-Zulu Natal, South Africa. A comprehensive polyphasic approach viz., cell morphology, pigment composition and complete genome sequence analysis was conducted to elucidate the taxonomic position of the isolated strain. The blue-green oval to rod-shaped cells were 14-18 μm in size, and contained a high amount of phycocyanin pigments. The strain was moderate thermotolerant/alkalitolerant halophile with the optimum conditions for growth at 35 °C, pH 8.5 and 120 g/l of NaCl. Based on 16S rRNA gene sequence phylogeny, the strain was related to members of the 'Euhalothece' subcluster (99%). The whole genome sequence was determined, and the annotated genes showed a 90% sequence similarity to the gas-vacuolate, spindle-shaped Dactylococcopsis salina PCC 8305. The size of the genome was determined to be 5,113,178 bp and contained 4332 protein-coding genes and 69 RNA genes with a G + C content of 46.7%. Genes encoding osmoregulation, oxidative stress, heat shock, persister cells, and UV-absorbing secondary metabolites, among others, were identified. Based on the phylogenetic analysis of the 16S rRNA gene sequences, physiological data, pigment compositions and genomic data, the strain is considered to represent a novel species of Euhalothece.

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.

Evolution of Proteins and Proteomes: A Phylogenetics Approach

The study of evolutionary relationships among protein sequences was one of the first applications of bioinformatics. Since then, and accompanying the wealth of biological data produced by genome sequencing and other high-throughput techniques, the use of bioinformatics in general and phylogenetics in particular has been gaining ground in the study of protein and proteome evolution. Nowadays, the use of phylogenetics is instrumental not only to infer the evolutionary relationships among species and their genome sequences, but also to reconstruct ancestral states of proteins and proteomes and hence trace the paths followed by evolution. Here I survey recent progress in the elucidation of mechanisms of protein and proteome evolution in which phylogenetics has played a determinant role.

Lack of the PGA exopolysaccharide in Salmonella as an adaptive trait for survival in the host

Many bacteria build biofilm matrices using a conserved exopolysaccharide named PGA or PNAG (poly-β-1,6-N-acetyl-D-glucosamine). Interestingly, while E. coli and other members of the family Enterobacteriaceae encode the pgaABCD operon responsible for PGA synthesis, Salmonella lacks it. The evolutionary force driving this difference remains to be determined. Here, we report that Salmonella lost the pgaABCD operon after the divergence of Salmonella and Citrobacter clades, and previous to the diversification of the currently sequenced Salmonella strains. Reconstitution of the PGA machinery endows Salmonella with the capacity to produce PGA in a cyclic dimeric GMP (c-di-GMP) dependent manner. Outside the host, the PGA polysaccharide does not seem to provide any significant benefit to Salmonella: resistance against chlorine treatment, ultraviolet light irradiation, heavy metal stress and phage infection remained the same as in a strain producing cellulose, the main biofilm exopolysaccharide naturally produced by Salmonella. In contrast, PGA production proved to be deleterious to Salmonella survival inside the host, since it increased susceptibility to bile salts and oxidative stress, and hindered the capacity of S. Enteritidis to survive inside macrophages and to colonize extraintestinal organs, including the gallbladder. Altogether, our observations indicate that PGA is an antivirulence factor whose loss may have been a necessary event during Salmonella speciation to permit survival inside the host.

During bacterial evolution, specific traits that optimize the organism’s fitness are selected. The production of exopolysaccharides is widespread among bacteria in which they play a protective shielding role as main constituents of biofilms. In contrast to closely related siblings, Salmonella has lost the capacity to produce the exopolysaccharide PGA. Our study reveals that Salmonella lost pga genes, and that the driving force for such a loss may have been the detrimental impact that PGA has during Salmonella invasion of internal organs where it augments the susceptibility to bile salts and oxygen radicals, reducing bacterial survival inside macrophages and rendering Salmonella avirulent. These results suggest that gene-loss has played an important role during Salmonella evolution.

MimiLook: A Phylogenetic Workflow for Detection of Gene Acquisition in Major Orthologous Groups of Megavirales

With the inclusion of new members, understanding about evolutionary mechanisms and processes by which members of the proposed order, Megavirales, have evolved has become a key area of interest. The central role of gene acquisition has been shown in previous studies. However, the major drawback in gene acquisition studies is the focus on few MV families or putative families with large variation in their genetic structure. Thus, here we have tried to develop a methodology by which we can detect horizontal gene transfers (HGTs), taking into consideration orthologous groups of distantly related Megavirale families. Here, we report an automated workflow MimiLook, prepared as a Perl command line program, that deduces orthologous groups (OGs) from ORFomes of Megavirales and constructs phylogenetic trees by performing alignment generation, alignment editing and protein-protein BLAST (BLASTP) searching across the National Center for Biotechnology Information (NCBI) non-redundant (nr) protein sequence database. Finally, this tool detects statistically validated events of gene acquisitions with the help of the T-REX algorithm by comparing individual gene tree with NCBI species tree. In between the steps, the workflow decides about handling paralogs, filtering outputs, identifying Megavirale specific OGs, detection of HGTs, along with retrieval of information about those OGs that are monophyletic with organisms from cellular domains of life. By implementing MimiLook, we noticed that nine percent of Megavirale gene families (i.e., OGs) have been acquired by HGT, 80% OGs were Megaviralespecific and eight percent were found to be sharing common ancestry with members of cellular domains (Eukaryote, Bacteria, Archaea, Phages or other viruses) and three percent were ambivalent. The results are briefly discussed to emphasize methodology. Also, MimiLook is relevant for detecting evolutionary scenarios in other targeted phyla with user defined modifications. It can be accessed at following link 10.6084/m9.figshare.4653622.

Reconstruction of the evolution of microbial defense systems

Background

Evolution of bacterial and archaeal genomes is a highly dynamic process that involves intensive loss of genes as well as gene gain via horizontal transfer, with a lesser contribution from gene duplication. The rates of these processes can be estimated by comparing genomes that are linked by an evolutionary tree. These estimated rates of genome dynamics events substantially differ for different functional classes of genes. The genes involved in defense against viruses and other invading DNA are among those that are gained and lost at the highest rates.

Results

We employed a stochastic birth-and-death model to obtain maximum likelihood estimates of the rates of gain and loss of defense genes in 35 groups of closely related bacterial genomes and one group of archaeal genomes. We find that on average, the defense genes experience 1.4 fold higher flux than the rest of microbial genes. This excessive flux of defense genes over the genomic mean is consistent across diverse microbial groups. The few exceptions include intracellular parasites with small, degraded genomes that possess few defense systems which are more stable than in other microbes. Generally, defense genes follow the previously established pattern of genome dynamics, with gene family loss being about 3 times more common than gain and an order of magnitude more common than expansion or contraction of gene families. Case by case analysis of the evolutionary dynamics of defense genes indicates frequent multiple events in the same locus and widespread involvement of mobile elements in the gain and loss of defense genes.

Conclusions

Evolution of microbial defense systems is highly dynamic but, notwithstanding the host-parasite arms race, generally follows the same trends that have been established for the rest of the genes. Apart from the paucity and the low flux of defense genes in parasitic bacteria with deteriorating genomes, there is no clear connection between the evolutionary regime of defense systems and microbial life style.

Electronic supplementary material

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

Bacterial intra-species gene loss occurs in a largely clocklike manner mostly within a pool of less conserved and constrained genes

Gene loss is a major contributor to the evolution of bacterial gene content. Gene loss may occur as a result of shifts in environment leading to changes in the intensity and/or directionality of selection applied for the maintenance of specific genes. Gene loss may also occur in a more neutral manner, when gene functions are lost that were not subject to strong selection to be maintained, irrespective of changes to environment. Here, we used a pangenome-based approach to investigate patterns of gene loss across 15 bacterial species. We demonstrate that gene loss tends to occur mostly within a pool of genes that are less constrained within species, even in those strains from which they are not lost, and less conserved across bacterial species. Our results indicate that shifts in selection, resulting from shifts in environment are not required to explain the majority of gene loss events occurring within a diverse collection of bacterial species. Caution should therefore be taken when attributing differences in gene content to differences in environment.

Phylogenetic signal in phenotypic traits related to carbon source assimilation and chemical sensitivity in Acinetobacter species

A common belief is that the phylogeny of bacteria may reflect molecular functions and phenotypic characteristics, pointing towards phylogenetic conservatism of traits. Here, we tested this hypothesis for a large set of Acinetobacter strains. Members of the genus Acinetobacter are widespread in nature, demonstrate a high metabolic diversity and are resistant to several environmental stressors. Notably, some species are known to cause opportunistic human infections. A total of 133 strains belonging to 33 species with validly published names, two genomic species and species of an as-yet unknown taxonomic status were analyzed using the GENIII technology of Biolog, which allows high-throughput phenotyping. We estimated the strength and significance of the phylogenetic signal of each trait across phylogenetic reconstructions based on partial RNA polymerase subunit B (rpoB) and core genome sequences. Secondly, we tested whether phylogenetic distance was a good predictor of trait differentiation by Mantel test analysis. And finally, evolutionary model fitting was used to determine if the data for each phenotypic character was consistent with a phylogenetic or an essentially random model of trait distribution. Our data revealed that some key phenotypic traits related to substrate assimilation and chemical sensitivity are linked to the phylogenetic placement of Acinetobacter species. The strongest phylogenetic signals found were for utilization of different carbon sources such as some organic acids, amino acids and sugars, thus suggesting that in the diversification of Acinetobacter carbon source assimilation has had a relevant role. Future work should be aimed to clarify how such traits have shaped the remarkable ability of this bacterial group to dominate in a wide variety of habitats.

Horizontal gene transfer: essentiality and evolvability in prokaryotes, and roles in evolutionary transitions

The wide spread of gene exchange and loss in the prokaryotic world has prompted the concept of ‘lateral genomics’ to the point of an outright denial of the relevance of phylogenetic trees for evolution. However, the pronounced coherence congruence of the topologies of numerous gene trees, particularly those for (nearly) universal genes, translates into the notion of a statistical tree of life (STOL), which reflects a central trend of vertical evolution. The STOL can be employed as a framework for reconstruction of the evolutionary processes in the prokaryotic world. Quantitatively, however, horizontal gene transfer (HGT) dominates microbial evolution, with the rate of gene gain and loss being comparable to the rate of point mutations and much greater than the duplication rate. Theoretical models of evolution suggest that HGT is essential for the survival of microbial populations that otherwise deteriorate due to the Muller’s ratchet effect. Apparently, at least some bacteria and archaea evolved dedicated vehicles for gene transfer that evolved from selfish elements such as plasmids and viruses. Recent phylogenomic analyses suggest that episodes of massive HGT were pivotal for the emergence of major groups of organisms such as multiple archaeal phyla as well as eukaryotes. Similar analyses appear to indicate that, in addition to donating hundreds of genes to the emerging eukaryotic lineage, mitochondrial endosymbiosis severely curtailed HGT. These results shed new light on the routes of evolutionary transitions, but caution is due given the inherent uncertainty of deep phylogenies.

Metagenomics and Bioinformatics in Microbial Ecology: Current Status and Beyond

Metagenomic approaches are now commonly used in microbial ecology to study microbial communities in more detail, including many strains that cannot be cultivated in the laboratory. Bioinformatic analyses make it possible to mine huge metagenomic datasets and discover general patterns that govern microbial ecosystems. However, the findings of typical metagenomic and bioinformatic analyses still do not completely describe the ecology and evolution of microbes in their environments. Most analyses still depend on straightforward sequence similarity searches against reference databases. We herein review the current state of metagenomics and bioinformatics in microbial ecology and discuss future directions for the field. New techniques will allow us to go beyond routine analyses and broaden our knowledge of microbial ecosystems. We need to enrich reference databases, promote platforms that enable meta- or comprehensive analyses of diverse metagenomic datasets, devise methods that utilize long-read sequence information, and develop more powerful bioinformatic methods to analyze data from diverse perspectives.

Functional profiling of cyanobacterial genomes and its role in ecological adaptations

With the availability of complete genome sequences of many cyanobacterial species, it is becoming feasible to study the broad prospective of the environmental adaptation and the overall changes at transcriptional and translational level in these organisms. In the evolutionary phase, niche-specific competitive forces have resulted in specific features of the cyanobacterial genomes. In this study, functional composition of the 84 different cyanobacterial genomes and their adaptations to different environments was examined by identifying the genomic composition for specific cellular processes, which reflect their genomic functional profile and ecological adaptation. It was identified that among cyanobacterial genomes, metabolic genes have major share over other categories and differentiation of genomic functional profile was observed for the species inhabiting different habitats. The cyanobacteria of freshwater and other habitats accumulate large number of poorly characterized genes. Strain specific functions were also reported in many cyanobacterial members, of which an important feature was the occurrence of phage-related sequences. From this study, it can be speculated that habitat is one of the major factors in giving the shape of functional composition of cyanobacterial genomes towards their ecological adaptations.

Reconstruction of ancestral chromosome architecture and gene repertoire reveals principles of genome evolution in a model yeast genus

Reconstructing genome history is complex but necessary to reveal quantitative principles governing genome evolution. Such reconstruction requires recapitulating into a single evolutionary framework the evolution of genome architecture and gene repertoire. Here, we reconstructed the genome history of the genus Lachancea that appeared to cover a continuous evolutionary range from closely related to more diverged yeast species. Our approach integrated the generation of a high-quality genome data set; the development of AnChro, a new algorithm for reconstructing ancestral genome architecture; and a comprehensive analysis of gene repertoire evolution. We found that the ancestral genome of the genus Lachancea contained eight chromosomes and about 5173 protein-coding genes. Moreover, we characterized 24 horizontal gene transfers and 159 putative gene creation events that punctuated species diversification. We retraced all chromosomal rearrangements, including gene losses, gene duplications, chromosomal inversions and translocations at single gene resolution. Gene duplications outnumbered losses and balanced rearrangements with 1503, 929, and 423 events, respectively. Gene content variations between extant species are mainly driven by differential gene losses, while gene duplications remained globally constant in all lineages. Remarkably, we discovered that balanced chromosomal rearrangements could be responsible for up to 14% of all gene losses by disrupting genes at their breakpoints. Finally, we found that nonsynonymous substitutions reached fixation at a coordinated pace with chromosomal inversions, translocations, and duplications, but not deletions. Overall, we provide a granular view of genome evolution within an entire eukaryotic genus, linking gene content, chromosome rearrangements, and protein divergence into a single evolutionary framework.

Adaptive evolution of complex innovations through stepwise metabolic niche expansion

A central challenge in evolutionary biology concerns the mechanisms by which complex metabolic innovations requiring multiple mutations arise. Here, we propose that metabolic innovations accessible through the addition of a single reaction serve as stepping stones towards the later establishment of complex metabolic features in another environment. We demonstrate the feasibility of this hypothesis through three complementary analyses. First, using genome-scale metabolic modelling, we show that complex metabolic innovations in Escherichia coli can arise via changing nutrient conditions. Second, using phylogenetic approaches, we demonstrate that the acquisition patterns of complex metabolic pathways during the evolutionary history of bacterial genomes support the hypothesis. Third, we show how adaptation of laboratory populations of E. coli to one carbon source facilitates the later adaptation to another carbon source. Our work demonstrates how complex innovations can evolve through series of adaptive steps without the need to invoke non-adaptive processes.

Conserved Units of Co-Expression in Bacterial Genomes: An Evolutionary Insight into Transcriptional Regulation

Genome-wide measurements of transcriptional activity in bacteria indicate that the transcription of successive genes is strongly correlated beyond the scale of operons. Here, we analyze hundreds of bacterial genomes to identify supra-operonic segments of genes that are proximal in a large number of genomes. We show that these synteny segments correspond to genomic units of strong transcriptional co-expression. Structurally, the segments contain operons with specific relative orientations (co-directional or divergent) and nucleoid-associated proteins are found to bind at their boundaries. Functionally, operons inside a same segment are highly co-expressed even in the apparent absence of regulatory factors at their promoter regions. Remote operons along DNA can also be co-expressed if their corresponding segments share a transcriptional or sigma factor, without requiring these factors to bind directly to the promoters of the operons. As evidence that these results apply across the bacterial kingdom, we demonstrate them both in the Gram-negative bacterium Escherichia coli and in the Gram-positive bacterium Bacillus subtilis. The underlying process that we propose involves only RNA-polymerases and DNA: it implies that the transcription of an operon mechanically enhances the transcription of adjacent operons. In support of a primary role of this regulation by facilitated co-transcription, we show that the transcription en bloc of successive operons as a result of transcriptional read-through is strongly and specifically enhanced in synteny segments. Finally, our analysis indicates that facilitated co-transcription may be evolutionary primitive and may apply beyond bacteria.

One step beyond a ribosome: The ancient anaerobic core

Life arose in a world without oxygen and the first organisms were anaerobes. Here we investigate the gene repertoire of the prokaryote common ancestor, estimating which genes it contained and to which lineages of modern prokaryotes it was most similar in terms of gene content. Using a phylogenetic approach we found that among trees for all 8779 protein families shared between 134 archaea and 1847 bacterial genomes, only 1045 have sequences from at least two bacterial and two archaeal groups and retain the ancestral archaeal–bacterial split. Among those, the genes shared by anaerobes were identified as candidate genes for the prokaryote common ancestor, which lived in anaerobic environments. We find that these anaerobic prokaryote common ancestor genes are today most frequently distributed among methanogens and clostridia, strict anaerobes that live from low free energy changes near the thermodynamic limit of life. The anaerobic families encompass genes for bifunctional acetyl-CoA-synthase/CO-dehydrogenase, heterodisulfide reductase subunits C and A, ferredoxins, and several subunits of the Mrp-antiporter/hydrogenase family, in addition to numerous S-adenosyl methionine (SAM) dependent methyltransferases. The data indicate a major role for methyl groups in the metabolism of the prokaryote common ancestor. The data furthermore indicate that the prokaryote ancestor possessed a rotor stator ATP synthase, but lacked cytochromes and quinones as well as identifiable redox-dependent ion pumping complexes. The prokaryote ancestor did possess, however, an Mrp-type H⁺/Na⁺ antiporter complex, capable of transducing geochemical pH gradients into biologically more stable Na⁺-gradients. The findings implicate a hydrothermal, autotrophic, and methyl-dependent origin of life. This article is part of a Special Issue entitled ‘EBEC 2016: 19th European Bioenergetics Conference, Riva del Garda, Italy, July 2–6, 2016’, edited by Prof. Paolo Bernardi.

•

Life arose without oxygen, the universal ancestor (Luca) was an anaerobe.

•

We used phylogenetic and physiological criteria to identify genes present in Luca.

•

An ancient core of 65 metabolic genes shed light on Luca's anaerobic lifestyle.

•

Ancient core genes are most widespread among modern methanogens and clostridia.

•

The data implicate a major role for methyl groups in Luca's anaerobic metabolism.

Origin and Evolution of Rickettsial Plasmids

Background

Rickettsia species are strictly intracellular bacteria that have undergone a reductive genomic evolution. Despite their allopatric lifestyle, almost half of the 26 currently validated Rickettsia species have plasmids. In order to study the origin, evolutionary history and putative roles of rickettsial plasmids, we investigated the evolutionary processes that have shaped 20 plasmids belonging to 11 species, using comparative genomics and phylogenetic analysis between rickettsial, microbial and non-microbial genomes.

Results

Plasmids were differentially present among Rickettsia species. The 11 species had 1 to 4 plasmid (s) with a size ranging from 12 kb to 83 kb. We reconstructed pRICO, the last common ancestor of the current rickettsial plasmids. pRICO was vertically inherited mainly from Rickettsia/Orientia chromosomes and diverged vertically into a single or multiple plasmid(s) in each species. These plasmids also underwent a reductive evolution by progressive gene loss, similar to that observed in rickettsial chromosomes, possibly leading to cryptic plasmids or complete plasmid loss. Moreover, rickettsial plasmids exhibited ORFans, recent gene duplications and evidence of horizontal gene transfer events with rickettsial and non-rickettsial genomes mainly from the α/γ-proteobacteria lineages. Genes related to maintenance and plasticity of plasmids, and to adaptation and resistance to stress mostly evolved under vertical and/or horizontal processes. Those involved in nucleotide/carbohydrate transport and metabolism were under the influence of vertical evolution only, whereas genes involved in cell wall/membrane/envelope biogenesis, cycle control, amino acid/lipid/coenzyme and secondary metabolites biosynthesis, transport and metabolism underwent mainly horizontal transfer events.

Conclusion

Rickettsial plasmids had a complex evolution, starting with a vertical inheritance followed by a reductive evolution associated with increased complexity via horizontal gene transfer as well as gene duplication and genesis. The plasmids are plastic and mosaic structures that may play biological roles similar to or distinct from their co-residing chromosomes in an obligate intracellular lifestyle.

In Silico Gene-Level Evolution Explains Microbial Population Diversity through Differential Gene Mobility

Microbial communities can show astonishing ecological and phylogenetic diversity. What is the role of pervasive horizontal gene transfer (HGT) in shaping this diversity in the presence of clonally expanding “killer strains”? Does HGT of antibiotic production and resistance genes erase phylogenetic structure? To answer these questions, we study a spatial eco-evolutionary model of prokaryotes, inspired by recent findings on antagonistic interactions in Vibrionaceae populations. We find toxin genes evolve to be highly mobile, whereas resistance genes minimize mobility. This differential gene mobility is a requirement to maintain a diverse and dynamic ecosystem. The resistance gene repertoire acts as a core genome that corresponds to the phylogeny of cells, whereas toxin genes do not follow this phylogeny and have a patchy distribution. We also show that interstrain HGT makes the emergent phylogenetic structure robust to selective sweeps. Finally, in this evolved ecosystem we observe antagonistic interactions between, rather than within, spatially structure subpopulations, as has been previously observed for prokaryotes in soils and oceans. In contrast to ascribing the diversification and evolution of microbial communities to clonal dynamics, we show that multilevel evolution can elegantly explain the observed phylogenetic structure and ecosystem diversity.

Horizontal acquisition of toxic alkaloid synthesis in a clade of plant associated fungi

•

Loline and ergot alkaloids are fundamental for Clavicipitaceae.

•

These secondary metabolism clusters are found in other, far related, species.

•

Phylogenetics establish a close relationship between clusters of different species.

•

Detection of large amounts of gene loss supports the HGT events.

Clavicipitaceae is a fungal group that comprises species that closely interact with plants as pathogens, parasites or symbionts. A key factor in these interactions is the ability of these fungi to synthesize toxic alkaloid compounds that contribute to the protection of the plant host against herbivores. Some of these compounds such as ergot alkaloids are toxic to humans and have caused important epidemics throughout history. The gene clusters encoding the proteins responsible for the synthesis of ergot alkaloids and lolines in Clavicipitaceae have been elucidated. Notably, homologs to these gene clusters can be found in distantly related species such as Aspergillus fumigatus and Penicillium expansum, which diverged from Clavicipitaceae more than 400 million years ago. We here use a phylogenetic approach to analyze the evolution of these gene clusters. We found that the gene clusters conferring the ability to synthesize ergot alkaloids and loline emerged first in Eurotiomycetes and were then likely transferred horizontally to Clavicipitaceae. Horizontal gene transfer is known to play a role in shaping the distribution of secondary metabolism clusters across distantly related fungal species. We propose that HGT events have played an important role in the capability of Clavicipitaceae to produce two key secondary metabolites that have enhanced the ability of these species to protect their plant hosts, therefore favoring their interactions.

Cave microbial community composition in oceanic islands: disentangling the effect of different colored mats in diversity patterns of Azorean lava caves

Processes determining diversity and composition of bacterial communities in island volcanic caves are still poorly understood. Here, we characterized colored microbial mats in 14 volcanic caves from two oceanic islands of the Azores using 16S rRNA gene sequences. Factors determining community diversity (α) and composition (β) were explored, namely colored mats, caves and islands, as well as environmental and chemical characteristics of caves. Additive partitioning of diversity using OTU occurrence showed a greater influence of β-diversity between islands and caves that may relate to differences in rare OTUs (singletons and doubletons) across scales. In contrast, Shannon diversity partitioning revealed the importance of the lowest hierarchical level (α diversity, colored mat), suggesting a dominance of cosmopolitan OTUs (>1%) in most samples. Cosmopolitan OTUs included members involved in nitrogen cycling, supporting the importance of this process in Azorean caves. Environmental and chemical conditions in caves did not show any significant relationship to OTU diversity and composition. The absence of clear differences between mat colors and across scales may be explained by (1) the geological youth of the cave system (cave communities have not had enough time to diverge) or/and (2) community convergence, as the result of selection pressure in extreme environments.

Gene Acquisitions from Bacteria at the Origins of Major Archaeal Clades Are Vastly Overestimated

In a recent article, Nelson-Sathi et al. (NS) report that the origins of major archaeal lineages (MAL) correspond to massive group-specific gene acquisitions via HGT from bacteria (Nelson-Sathi et al. 2015. Origins of major archaeal clades correspond to gene acquisitions from bacteria. Nature 517(7532):77-80.). If correct, this would have fundamental implications for the process of diversification in microbes. However, a reexamination of these data and results shows that the methodology used by NS systematically inflates the number of genes acquired at the root of each MAL, and incorrectly assumes bacterial origins for these genes. A reanalysis of their data with appropriate phylogenetic models accounting for the dynamics of gene gain and loss between lineages supports the continuous acquisition of genes over long periods in the evolution of Archaea.

Complementarity among plant growth promoting traits in rhizospheric bacterial communities promotes plant growth

An assessment of roles of rhizospheric microbial diversity in plant growth is helpful in understanding plant-microbe interactions. Using random combinations of rhizospheric bacterial species at different richness levels, we analysed the contribution of species richness, compositions, interactions and identity on soil microbial respiration and plant biomass. We showed that bacterial inoculation in plant rhizosphere enhanced microbial respiration and plant biomass with complementary relationships among bacterial species. Plant growth was found to increase linearly with inoculation of rhizospheric bacterial communities with increasing levels of species or plant growth promoting trait diversity. However, inoculation of diverse bacterial communities having single plant growth promoting trait, i.e., nitrogen fixation could not enhance plant growth over inoculation of single bacteria. Our results indicate that bacterial diversity in rhizosphere affect ecosystem functioning through complementary relationship among plant growth promoting traits and may play significant roles in delivering microbial services to plants.

Event inference in multidomain families with phylogenetic reconciliation

BACKGROUND

Reconstructing evolution provides valuable insights into the processes of gene evolution and function. However, while there have been great advances in algorithms and software to reconstruct the history of gene families, these tools do not model the domain shuffling events (domain duplication, insertion, transfer, and deletion) that drive the evolution of multidomain protein families. Protein evolution through domain shuffling events allows for rapid exploration of functions by introducing new combinations of existing folds. This powerful mechanism was key to some significant evolutionary innovations, such as multicellularity and the vertebrate immune system. A method for reconstructing this important evolutionary process is urgently needed.

RESULTS

Here, we introduce a novel, event-based framework for studying multidomain evolution by reconciling a domain tree with a gene tree, with additional information provided by the species tree. In the context of this framework, we present the first reconciliation algorithms to infer domain shuffling events, while addressing the challenges inherent in the inference of evolution across three levels of organization.

CONCLUSIONS

We apply these methods to the evolution of domains in the Membrane associated Guanylate Kinase family. These case studies reveal a more vivid and detailed evolutionary history than previously provided. Our algorithms have been implemented in software, freely available at http://www.cs.cmu.edu/˜durand/Notung.

Evolutionary history influences the salinity preference of bacterial taxa in wetland soils

Salinity is a major driver of bacterial community composition across the globe. Despite growing recognition that different bacterial species are present or active at different salinities, the mechanisms by which salinity structures community composition remain unclear. We tested the hypothesis that these patterns reflect ecological coherence in the salinity preferences of phylogenetic groups using a reciprocal transplant experiment of fresh- and saltwater wetland soils. The salinity of both the origin and host environments affected community composition (16S rRNA gene sequences) and activity (CO2 and CH4 production, and extracellular enzyme activity). These changes in community composition and activity rates were strongly correlated, which suggests the effect of environment on function could be mediated, at least in part, by microbial community composition. Based on their distribution across treatments, each phylotype was categorized as having a salinity preference (freshwater, saltwater, or none) and phylogenetic analyses revealed a significant influence of evolutionary history on these groupings. This finding was corroborated by examining the salinity preferences of high-level taxonomic groups. For instance, we found that the majority of α- and γ-proteobacteria in these wetland soils preferred saltwater, while many β-proteobacteria prefer freshwater. Overall, our results indicate the effect of salinity on bacterial community composition results from phylogenetically-clustered salinity preferences.