Phylogenetic analyses of cyanobacterial genomes: quantification of horizontal gene transfer events

We are doing it wrong: Putting homology before phylogeny in cyanobacterial taxonomy

The rapid expansion of whole genome sequencing in bacterial taxonomy has revealed deep evolutionary relationships and speciation signals, but assembly methods often miss true nucleotide diversity in the ribosomal operons. Though it lacks sufficient phylogenetic signal at the species level, the 16S ribosomal RNA gene is still much used in bacterial taxonomy. In cyanobacterial taxonomy, comparisons of 16S-23S Internal Transcribed Spacer (ITS) regions are used to bridge this information gap. Although ITS rRNA region analyses are routinely being used to identify species, researchers often do not identify orthologous operons, which leads to improper comparisons. No method for delineating orthologous operon copies from paralogous ones has been established. A new method for recognizing orthologous ribosomal operons by quantifying the conserved paired nucleotides in a helical domain of the ITS, has been developed. The D1' Index quantifies differences in the ratio of pyrimidines to purines in paired nucleotide sequences of this helix. Comparing 111 operon sequences from 89 strains of Brasilonema, four orthologous operon types were identified. Plotting D1' Index values against the length of helices produced clear separation of orthologs. Most orthologous operons in this study were observed both with and without tRNA genes present. We hypothesize that genomic rearrangement, not gene duplication, is responsible for the variation among orthologs. This new method will allow cyanobacterial taxonomists to utilize ITS rRNA region data more correctly, preventing erroneous taxonomic hypotheses. Moreover, this work could assist genomicists in identifying and preserving evident sequence variability in ribosomal operons, which is an important proxy for evolution in prokaryotes.

Synechococcus elongatus Argonaute reduces natural transformation efficiency and provides immunity against exogenous plasmids

IMPORTANCE

S. elongatus is an important cyanobacterial model organism for the study of its prokaryotic circadian clock, photosynthesis, and other biological processes. It is also widely used for genetic engineering to produce renewable biochemicals. Our findings reveal an SeAgo-based defense mechanism in S. elongatus against the horizontal transfer of genetic material. We demonstrate that deletion of the ago gene facilitates genetic studies and genetic engineering of S. elongatus.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Structures of Get3d reveal a distinct architecture associated with the emergence of photosynthesis

Homologs of the protein Get3 have been identified in all domains yet remain to be fully characterized. In the eukaryotic cytoplasm, Get3 delivers tail-anchored (TA) integral membrane proteins, defined by a single transmembrane helix at their C terminus, to the endoplasmic reticulum. While most eukaryotes have a single Get3 gene, plants are notable for having multiple Get3 paralogs. Get3d is conserved across land plants and photosynthetic bacteria and includes a distinctive C-terminal α-crystallin domain. After tracing the evolutionary origin of Get3d, we solve the Arabidopsis thaliana Get3d crystal structure, identify its localization to the chloroplast, and provide evidence for a role in TA protein binding. The structure is identical to that of a cyanobacterial Get3 homolog, which is further refined here. Distinct features of Get3d include an incomplete active site, a "closed" conformation in the apo-state, and a hydrophobic chamber. Both homologs have ATPase activity and are capable of binding TA proteins, supporting a potential role in TA protein targeting. Get3d is first found with the development of photosynthesis and conserved across 1.2 billion years into the chloroplasts of higher plants across the evolution of photosynthesis suggesting a role in the homeostasis of photosynthetic machinery.

Current Methods for Recombination Detection in Bacteria

The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.

Genetic engineering of marine cyanophages reveals integration but not lysogeny in T7-like cyanophages

Marine cyanobacteria of the genera Synechococcus and Prochlorococcus are the most abundant photosynthetic organisms on earth, spanning vast regions of the oceans and contributing significantly to global primary production. Their viruses (cyanophages) greatly influence cyanobacterial ecology and evolution. Although many cyanophage genomes have been sequenced, insight into the functional role of cyanophage genes is limited by the lack of a cyanophage genetic engineering system. Here, we describe a simple, generalizable method for genetic engineering of cyanophages from multiple families, that we named REEP for REcombination, Enrichment and PCR screening. This method enables direct investigation of key cyanophage genes, and its simplicity makes it adaptable to other ecologically relevant host-virus systems. T7-like cyanophages often carry integrase genes and attachment sites, yet exhibit lytic infection dynamics. Here, using REEP, we investigated their ability to integrate and maintain a lysogenic life cycle. We found that these cyanophages integrate into the host genome and that the integrase and attachment site are required for integration. However, stable lysogens did not form. The frequency of integration was found to be low in both lab cultures and the oceans. These findings suggest that T7-like cyanophage integration is transient and is not part of a classical lysogenic cycle.

Phylum Gemmatimonadota and Its Role in the Environment

Bacteria are an important part of every ecosystem that they inhabit on Earth. Environmental microbiologists usually focus on a few dominant bacterial groups, neglecting less abundant ones, which collectively make up most of the microbial diversity. One of such less-studied phyla is Gemmatimonadota. Currently, the phylum contains only six cultured species. However, data from culture-independent studies indicate that members of Gemmatimonadota are common in diverse habitats. They are abundant in soils, where they seem to be frequently associated with plants and the rhizosphere. Moreover, Gemmatimonadota were found in aquatic environments, such as freshwaters, wastewater treatment plants, biofilms, and sediments. An important discovery was the identification of purple bacterial reaction centers and anoxygenic photosynthesis in this phylum, genes for which were likely acquired via horizontal gene transfer. So far, the capacity for anoxygenic photosynthesis has been described for two cultured species: Gemmatimonas phototrophica and Gemmatimonas groenlandica. Moreover, analyses of metagenome-assembled genomes indicate that it is also common in uncultured lineages of Gemmatimonadota. This review summarizes the current knowledge about this understudied bacterial phylum with an emphasis on its environmental distribution.

Thirtieth Anniversary of the Discovery of Laxaphycins. Intriguing Peptides Keeping a Part of Their Mystery

Lipopeptides are a class of compounds generally produced by microorganisms through hybrid biosynthetic pathways involving non-ribosomal peptide synthase and a polyketyl synthase. Cyanobacterial-produced laxaphycins are examples of this family of compounds that have expanded over the past three decades. These compounds benefit from technological advances helping in their synthesis and characterization, as well as in deciphering their biosynthesis. The present article attempts to summarize most of the articles that have been published on laxaphycins. The current knowledge on the ecological role of these complex sets of compounds will also be examined.

Systematic Detection of Large-Scale Multigene Horizontal Transfer in Prokaryotes

Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the “scale” of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.

Temperature Responses of Heterotrophic Bacteria in Co-culture With a Red Sea Synechococcus Strain

Interactions between autotrophic and heterotrophic bacteria are fundamental for marine biogeochemical cycling. How global warming will affect the dynamics of these essential microbial players is not fully understood. The aims of this study were to identify the major groups of heterotrophic bacteria present in a Synechococcus culture originally isolated from the Red Sea and assess their joint responses to experimental warming within the metabolic ecology framework. A co-culture of Synechococcus sp. RS9907 and their associated heterotrophic bacteria, after determining their taxonomic affiliation by 16S rRNA gene sequencing, was acclimated and maintained in the lab at different temperatures (24-34°C). The abundance and cellular properties of Synechococcus and the three dominant heterotrophic bacterial groups (pertaining to the genera Paracoccus, Marinobacter, and Muricauda) were monitored by flow cytometry. The activation energy of Synechococcus, which grew at 0.94-1.38 d-1, was very similar (0.34 ± 0.02 eV) to the value hypothesized by the metabolic theory of ecology (MTE) for autotrophs (0.32 eV), while the values of the three heterotrophic bacteria ranged from 0.16 to 1.15 eV and were negatively correlated with their corresponding specific growth rates (2.38-24.4 d-1). The corresponding carrying capacities did not always follow the inverse relationship with temperature predicted by MTE, nor did we observe a consistent response of bacterial cell size and temperature. Our results show that the responses to future ocean warming of autotrophic and heterotrophic bacteria in microbial consortia might not be well described by theoretical universal rules.

Adaptive responses of free-living and symbiotic microalgae to simulated future ocean conditions

Marine microalgae are a diverse group of microscopic eukaryotic and prokaryotic organisms capable of photosynthesis. They are important primary producers and carbon sinks but their physiology and persistence are severely affected by global climate change. Powerful experimental evolution technologies are being used to examine the potential of microalgae to respond adaptively to current and predicted future conditions, as well as to develop resources to facilitate species conservation and restoration of ecosystem functions. This review synthesizes findings and insights from experimental evolution studies of marine microalgae in response to elevated temperature and/or pCO₂ . Adaptation to these environmental conditions has been observed in many studies of marine dinoflagellates, diatoms and coccolithophores. An enhancement in traits such as growth and photo-physiological performance and an increase in upper thermal limit have been shown to be possible, although the extent and rate of change differ between microalgal taxa. Studies employing multiple monoclonal replicates showed variation in responses among replicates and revealed the stochasticity of mutations. The work to date is already providing valuable information on species' climate sensitivity or resilience to managers and policymakers but extrapolating these insights to ecosystem- and community-level impacts continues to be a challenge. We recommend future work should include in situ experiments, diurnal and seasonal fluctuations, multiple drivers and multiple starting genotypes. Fitness trade-offs, stable versus plastic responses and the genetic bases of the changes also need investigating, and the incorporation of genome resequencing into experimental designs will be invaluable.

Reacquisition of light-harvesting genes in a marine cyanobacterium confers a broader solar niche

The evolution of phenotypic plasticity, i.e., the environmental induction of alternative phenotypes by the same genotype, can be an important mechanism of biological diversification.^1,2 For example, an evolved increase in plasticity may promote ecological niche expansion as well as the innovation of novel traits;³ however, both the role of phenotypic plasticity in adaptive evolution and its underlying mechanisms are still poorly understood.^4,5 Here, we report that the Chlorophyll d-producing marine cyanobacterium Acaryochloris marina strain MBIC11017 has evolved greater photosynthetic plasticity by reacquiring light-harvesting genes via horizontal gene transfer. The genes, which had been lost by the A. marina ancestor, are involved in the production and degradation of the light-harvesting phycobiliprotein phycocyanin. A. marina MBIC11017 exhibits a high degree of wavelength-dependence in phycocyanin production, and this ability enables it to grow with yellow and green light wavelengths that are inaccessible to other A. marina. Consequently, this strain has a broader solar niche than its close relatives. We discuss the role of horizontal gene transfer for regaining a lost phenotype in light of Dollo's Law⁶ that the loss of a complex trait is irreversible.

Genetic, Genomics, and Responses to Stresses in Cyanobacteria: Biotechnological Implications

Cyanobacteria are widely-diverse, environmentally crucial photosynthetic prokaryotes of great interests for basic and applied science. Work to date has focused mostly on the three non-nitrogen fixing unicellular species Synechocystis PCC 6803, Synechococcus PCC 7942, and Synechococcus PCC 7002, which have been selected for their genetic and physiological interests summarized in this review. Extensive "omics" data sets have been generated, and genome-scale models (GSM) have been developed for the rational engineering of these cyanobacteria for biotechnological purposes. We presently discuss what should be done to improve our understanding of the genotype-phenotype relationships of these models and generate robust and predictive models of their metabolism. Furthermore, we also emphasize that because Synechocystis PCC 6803, Synechococcus PCC 7942, and Synechococcus PCC 7002 represent only a limited part of the wide biodiversity of cyanobacteria, other species distantly related to these three models, should be studied. Finally, we highlight the need to strengthen the communication between academic researchers, who know well cyanobacteria and can engineer them for biotechnological purposes, but have a limited access to large photobioreactors, and industrial partners who attempt to use natural or engineered cyanobacteria to produce interesting chemicals at reasonable costs, but may lack knowledge on cyanobacterial physiology and metabolism.

Genome Mining and Evolutionary Analysis Reveal Diverse Type III Polyketide Synthase Pathways in Cyanobacteria

Cyanobacteria are prolific producers of natural products, including polyketides and hybrid compounds thereof. Type III polyketide synthases (PKSs) are of particular interest, due to their wide substrate specificity and simple reaction mechanism, compared with both type I and type II PKSs. Surprisingly, only two type III PKS products, hierridins, and (7.7)paracyclophanes, have been isolated from cyanobacteria. Here, we report the mining of 517 cyanobacterial genomes for type III PKS biosynthesis gene clusters. Approximately 17% of the genomes analyzed encoded one or more type III PKSs. Together with already characterized type III PKSs, the phylogeny of this group of enzymes was investigated. Our analysis showed that type III PKSs in cyanobacteria evolved into three major lineages, including enzymes associated with 1) (7.7)paracyclophane-like biosynthesis gene clusters, 2) hierridin-like biosynthesis gene clusters, and 3) cytochrome b5 genes. The evolutionary history of these enzymes is complex, with some sequences partitioning primarily according to speciation and others putatively according to their reaction type. Protein modeling showed that cyanobacterial type III PKSs generally have a smaller active site cavity (mean = 109.035 ų) compared with enzymes from other organisms. The size of the active site did not correlate well with substrate size, however, the “Gatekeeper” amino acid residues within the active site were strongly correlated to enzyme phylogeny. Our study provides unprecedented insight into the distribution, diversity, and molecular evolution of cyanobacterial type III PKSs, which could facilitate the discovery, characterization, and exploitation of novel enzymes, biochemical pathways, and specialized metabolites from this biosynthetically talented clade of microorganisms.

Complex Evolution of Light-Dependent Protochlorophyllide Oxidoreductases in Aerobic Anoxygenic Phototrophs: Origin, Phylogeny, and Function

Light-dependent protochlorophyllide oxidoreductase (LPOR) and dark-operative protochlorophyllide oxidoreductase are evolutionary and structurally distinct enzymes that are essential for the synthesis of (bacterio)chlorophyll, the primary pigment needed for both anoxygenic and oxygenic photosynthesis. In contrast to the long-held hypothesis that LPORs are only present in oxygenic phototrophs, we recently identified a functional LPOR in the aerobic anoxygenic phototrophic bacterium (AAPB) Dinoroseobacter shibae and attributed its presence to a single horizontal gene transfer event from cyanobacteria. Here, we provide evidence for the more widespread presence of genuine LPOR enzymes in AAPBs. An exhaustive bioinformatics search identified 36 putative LPORs outside of oxygenic phototrophic bacteria (cyanobacteria) with the majority being AAPBs. Using in vitro and in vivo assays, we show that the large majority of the tested AAPB enzymes are genuine LPORs. Solution structural analyses, performed for two of the AAPB LPORs, revealed a globally conserved structure when compared with a well-characterized cyanobacterial LPOR. Phylogenetic analyses suggest that LPORs were transferred not only from cyanobacteria but also subsequently between proteobacteria and from proteobacteria to Gemmatimonadetes. Our study thus provides another interesting example for the complex evolutionary processes that govern the evolution of bacteria, involving multiple horizontal gene transfer events that likely occurred at different time points and involved different donors.

Contribution of horizontal gene transfer to the functionality of microbial biofilm on a macroalgae

Horizontal gene transfer (HGT) is thought to be an important driving force for microbial evolution and niche adaptation and has been show in vitro to occur frequently in biofilm communities. However, the extent to which HGT takes place and what functions are being transferred in more complex and natural biofilm systems remains largely unknown. To address this issue, we investigated here HGT and enrichment of gene functions in the biofilm community of the common kelp (macroalgae) Ecklonia radiata in comparison to microbial communities in the surrounding seawater. We found that HGTs in the macroalgal biofilms were dominated by transfers between bacterial members of the same class or order and frequently involved genes for nutrient transport, sugar and phlorotannin degradation as well as stress responses, all functions that would be considered beneficial for bacteria living in this particular niche. HGT did not appear to be driven by mobile gene elements, indicating rather an involvement of unspecific DNA uptake (e.g. natural transformation). There was also a low overlap between the gene functions subject to HGT and those enriched in the biofilm community in comparison to planktonic community members. This indicates that much of the functionality required for bacteria to live in an E. radiata biofilm might be derived from vertical or environmental transmissions of symbionts. This study enhances our understanding of the relative role of evolutionary and ecological processes in driving community assembly and genomic diversity of biofilm communities.

Exploring syntenic conservation across genomes for phylogenetic studies of organisms subjected to horizontal gene transfers: A case study with Cyanobacteria and cyanolichens

Understanding the evolutionary history of symbiotic Cyanobacteria at a fine scale is essential to unveil patterns of associations with their hosts and factors driving their spatiotemporal interactions. As for bacteria in general, Horizontal Gene Transfers (HGT) are expected to be rampant throughout their evolution, which justified the use of single-locus phylogenies in macroevolutionary studies of these photoautotrophic bacteria. Genomic approaches have greatly increased the amount of molecular data available, but the selection of orthologous, congruent genes that are more likely to reflect bacterial macroevolutionary histories remains problematic. In this study, we developed a synteny-based approach and searched for Collinear Orthologous Regions (COR), under the assumption that genes that are present in the same order and orientation across a wide monophyletic clade are less likely to have undergone HGT. We searched sixteen reference Nostocales genomes and identified 99 genes, part of 28 COR comprising three to eight genes each. We then developed a bioinformatic pipeline, designed to minimize inter-genome contamination and processed twelve Nostoc-associated lichen metagenomes. This reduced our original dataset to 90 genes representing 25 COR, which were used to infer phylogenetic relationships within Nostocales and among lichenized Cyanobacteria. This dataset was narrowed down further to 71 genes representing 22 COR by selecting only genes part of one (largest) operon per COR. We found a relatively high level of congruence among trees derived from the 90-gene dataset, but congruence was only slightly higher among genes within a COR compared to genes across COR. However, topological congruence was significantly higher among the 71 genes part of one operon per COR. Nostocales phylogenies resulting from concatenation and species tree approaches based on the 90- and 71-gene datasets were highly congruent, but the most highly supported result was obtained when using synteny, collinearity, and operon information (i.e., 71-gene dataset) as gene selection criteria, which outperformed larger datasets with more genes.

Cyanobacterial blooms in wastewater treatment facilities: Significance and emerging monitoring strategies

Municipal wastewater treatment facilities (WWTFs) are prone to the proliferation of cyanobacterial species which thrive in stable, nutrient-rich environments. Dense cyanobacterial blooms frequently disrupt treatment processes and the supply of recycled water due to their production of extracellular polymeric substances, which hinder microfiltration, and toxins, which pose a health risk to end-users. A variety of methods are employed by water utilities for the identification and monitoring of cyanobacteria and their toxins in WWTFs, including microscopy, flow cytometry, ELISA, chemoanalytical methods, and more recently, molecular methods. Here we review the literature on the occurrence and significance of cyanobacterial blooms in WWTFs and discuss the pros and cons of the various strategies for monitoring these potentially hazardous events. Particular focus is directed towards next-generation metagenomic sequencing technologies for the development of site-specific cyanobacterial bloom management strategies. Long-term multi-omic observations will enable the identification of indicator species and the development of site-specific bloom dynamics models for the mitigation and management of cyanobacterial blooms in WWTFs. While emerging metagenomic tools could potentially provide deep insight into the diversity and flux of problematic cyanobacterial species in these systems, they should be considered a complement to, rather than a replacement of, quantitative chemoanalytical approaches.

A horizontally acquired expansin gene increases virulence of the emerging plant pathogen Erwinia tracheiphila

Erwinia tracheiphila is a bacterial plant pathogen that causes a fatal wilt infection in some cucurbit crop plants. Wilt symptoms are thought to be caused by systemic bacterial colonization through xylem that impedes sap flow. However, the genetic determinants of within-plant movement are unknown for this pathogen species. Here, we find that E. tracheiphila has horizontally acquired an operon with a microbial expansin (exlx) gene adjacent to a glycoside hydrolase family 5 (gh5) gene. Plant inoculation experiments with deletion mutants in the individual genes (Δexlx and Δgh5) and the full operon (Δexlx-gh5) resulted in decreased severity of wilt symptoms, decreased mortality rate, and impaired systemic colonization compared to the Wt strain. Co-inoculation experiments with Wt and Δexlx-gh5 rescued the movement defect of the mutant strain, suggesting that expansin and GH5 function extracellularly. Together, these results show that expansin-GH5 contributes to systemic movement through xylem, leading to rapid wilt symptom development and higher rates of plant death. The presence of expansin genes in diverse species of bacterial and fungal wilt-inducing pathogens suggests that microbial expansin proteins may be an under-appreciated virulence factor for many pathogen species.

A New Phylogenomic Approach For Quantifying Horizontal Gene Transfer Trends in Prokaryotes

It is well established nowadays that among prokaryotes, various families of orthologous genes exhibit conflicting evolutionary history. A prime factor for this conflict is horizontal gene transfer (HGT) - the transfer of genetic material not via vertical descent. Thus, the prevalence of HGT is challenging the meaningfulness of the classical Tree of Life concept. Here we present a comprehensive study of HGT representing the entire prokaryotic world. We mainly rely on a novel analytic approach for analyzing an aggregate of gene histories, by means of the quartet plurality distribution (QPD) that we develop. Through the analysis of real and simulated data, QPD is used to reveal evidence of a barrier against HGT, separating the archaea from the bacteria and making HGT between the two domains, in general, quite rare. In contrast, bacteria's confined HGT is substantially more frequent than archaea's. Our approach also reveals that despite intensive HGT, a strong tree-like signal can be extracted, corroborating several previous works. Thus, QPD, which enables one to analytically combine information from an aggregate of gene trees, can be used for understanding patterns and rates of HGT in prokaryotes, as well as for validating or refuting models of horizontal genetic transfers and evolution in general.

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.

Assessment of the Chemical Diversity and Potential Toxicity of Benthic Cyanobacterial Blooms in the Lagoon of Moorea Island (French Polynesia)

In the last decades, an apparent increase in the frequency of benthic cyanobacterial blooms has occurred in coral reefs and tropical lagoons, possibly in part because of global change and anthropogenic activities. In the frame of the survey of marine benthic cyanobacteria proliferating in the lagoon of Moorea Island (French Polynesia), 15 blooms were collected, mainly involving three species—Anabaena sp.1, Lyngbya majuscula and Hydrocoleum majus-B. Their chemical fingerprints, obtained through high performance liquid chromatography combined with UV detection and mass spectrometry (HPLC-UV-MS) analyses, revealed a high extent of species-specificity. The chemical profile of Anabaena sp.1 was characterized by three major cyclic lipopeptides of the laxaphycin family, whereas the one of L. majuscula was characterized by a complex mixture including tiahuramides, trungapeptins and serinol-derived malyngamides. Toxicity screening analyses conducted on these cyanobacterial samples using Artemia salina and mouse neuroblastoma cell-based (CBA-N2a) cytotoxic assays failed to show any toxicity to a degree that would merit risk assessment with regard to public health. However, the apparently increasing presence of blooms of Lyngbya, Hydrocoleum, Anabaena or other benthic cyanobacteria on coral reefs in French Polynesia encourages the implementation of ad hoc monitoring programs for the surveillance of their proliferation and potential assessment of associated hazards.

The circadian clock and darkness control natural competence in cyanobacteria

The cyanobacterium Synechococcus elongatus is a model organism for the study of circadian rhythms. It is naturally competent for transformation—that is, it takes up DNA from the environment, but the underlying mechanisms are unclear. Here, we use a genome-wide screen to identify genes required for natural transformation in S. elongatus, including genes encoding a conserved Type IV pilus, genes known to be associated with competence in other bacteria, and others. Pilus biogenesis occurs daily in the morning, while natural transformation is maximal when the onset of darkness coincides with the dusk circadian peak. Thus, the competence state in cyanobacteria is regulated by the circadian clock and can adapt to seasonal changes of day length.

The cyanobacterium Synechococcus elongatus is a model organism for the study of circadian rhythms, and is naturally competent for transformation. Here, Taton et al. identify genes required for natural transformation in this organism, and show that the coincidence of circadian dusk and darkness regulates the competence state in different day lengths.

Networks Consolidate the Core Concepts of Evolution by Natural Selection

Microbiology has unraveled rich evidence of ongoing reticulate evolutionary processes and complex interactions both within and between cells. These phenomena feature real biological networks, which can logically be analyzed using network-based tools. It is thus not surprising that network sciences, a field independent from evolutionary biology and microbiology, have recently pervasively infused their methods into both fields. Importantly, network tools bring forward observations enhancing the understanding of three core evolutionary concepts: variation, fitness, and heredity. Consequently, our work shows how network sciences can enhance evolutionary theory by explaining the evolution by natural selection of a broad diversity of units of selection, while updating the popular figure of Darwin's tree of life with a comprehensive sketch of the networks of evolution.

First-tier detection of intragenomic 16S rRNA gene variation in culturable endophytic bacteria from cacao seeds

Background

Intragenomic variability in 16S rDNA is a limiting factor for taxonomic and diversity characterization of Bacteria, and studies on its occurrence in natural/environmental populations are scarce. In this work, direct DNA amplicon sequencing coupled with frequent-cutter restriction analysis allowed detection of intragenomic 16S rDNA variation in culturable endophytic bacteria from cacao seeds in a fast and attractive manner.

Methods

Total genomic DNA from 65 bacterial strains was extracted and the 16S rDNA hyper variable V5-V9 regions were amplified for enzyme digestion and direct Sanger-type sequencing. The resulting electropherograms were visually inspected and compared to the corresponding AluI-restriction profiles, as well as to complete genome sequences in databases. Restriction analysis were employed to substitute the need of amplicon cloning and re-sequencing. A specifically improved polyacrylamide-gradient electrophoresis allowed to resolve 5-bp differences in restriction fragment sizes. Chi-square analysis on 2 × 2 contingency table tested for the independence between the 'number of AluI bands' and 'type of eletropherogram'.

Results

Two types of electropherograms were obtained: unique template, with single peaks per base (clean chromatograms), and heterogeneous template, with various levels of multiple peaks per base (mixed chromatograms). Statistics revealed significant interaction between number of restriction fragments and type of electropherogram for the same amplicons: clean or mixed ones associated to ≤5 or ≥6 bands, respectively. The mixed-template pattern combined with the AluI-restriction profiles indicated a high proportion of 49% of the culturable endophytes from a tropical environment showing evidence of intragenomic 16S rDNA heterogeneity.

Conclusion

The approach presented here was useful for a rapid, first-tier detection of intragenomic variation in culturable isolates, which can be applied in studies of other natural populations; a preliminary view of intragenomic heterogeneity levels can complement culture-dependent and -independent methods. Consequences of these findings in taxonomic and diversity studies in complex bacterial communities are discussed.

Quartet-Based Computations of Internode Certainty Provide Robust Measures of Phylogenetic Incongruence

Incongruence, or topological conflict, is prevalent in genome-scale data sets. Internode certainty (IC) and related measures were recently introduced to explicitly quantify the level of incongruence of a given internal branch among a set of phylogenetic trees and complement regular branch support measures (e.g., bootstrap, posterior probability) that instead assess the statistical confidence of inference. Since most phylogenomic studies contain data partitions (e.g., genes) with missing taxa and IC scores stem from the frequencies of bipartitions (or splits) on a set of trees, IC score calculation typically requires adjusting the frequencies of bipartitions from these partial gene trees. However, when the proportion of missing taxa is high, the scores yielded by current approaches that adjust bipartition frequencies in partial gene trees differ substantially from each other and tend to be overestimates. To overcome these issues, we developed three new IC measures based on the frequencies of quartets, which naturally apply to both complete and partial trees. Comparison of our new quartet-based measures to previous bipartition-based measures on simulated data shows that: (1) on complete data sets, both quartet-based and bipartition-based measures yield very similar IC scores; (2) IC scores of quartet-based measures on a given data set with and without missing taxa are more similar than the scores of bipartition-based measures; and (3) quartet-based measures are more robust to the absence of phylogenetic signal and errors in phylogenetic inference than bipartition-based measures. Additionally, the analysis of an empirical mammalian phylogenomic data set using our quartet-based measures reveals the presence of substantial levels of incongruence for numerous internal branches. An efficient open-source implementation of these quartet-based measures is freely available in the program QuartetScores (https://github.com/lutteropp/QuartetScores).

Life at Home and on the Roam: Genomic Adaptions Reflect the Dual Lifestyle of an Intracellular, Facultative Symbiont

"Candidatus Synechococcus feldmannii" is a facultative intracellular symbiont of the Atlanto-Mediterranean sponge Petrosia ficiformis. Genomic information of sponge-associated cyanobacteria derives thus far from the obligate and extracellular symbiont "Candidatus Synechococcus spongiarum." Here we utilized a differential methylation-based approach for bacterial DNA enrichment combined with metagenomics to obtain the first draft genomes of "Ca. Synechococcus feldmannii." By comparative genomics, we revealed that some genomic features (e.g., iron transport mediated by siderophores, eukaryotic-like proteins, and defense mechanisms, like CRISPR-Cas [clustered regularly interspaced short palindromic repeats-associated proteins]) are unique to both symbiont types and absent or rare in the genomes of taxonomically related free-living cyanobacteria. These genomic features likely enable life under the conditions found inside the sponge host. Interestingly, there are many genomic features that are shared by "Ca. Synechococcus feldmannii" and free-living cyanobacteria, while they are absent in the obligate symbiont "Ca. Synechococcus spongiarum." These include genes related to cell surface structures, genetic regulation, and responses to environmental stress, as well as the composition of photosynthetic genes and DNA metabolism. We speculate that the presence of these genes confers on "Ca. Synechococcus feldmannii" its facultative nature (i.e., the ability to respond to a less stable environment when free-living). Our comparative analysis revealed that distinct genomic features depend on the nature of the symbiotic interaction: facultative and intracellular versus obligate and extracellular. IMPORTANCE Given the evolutionary position of sponges as one of the earliest phyla to depart from the metazoan stem lineage, studies on their distinct and exceptionally diverse microbial communities should yield a better understanding of the origin of animal-bacterium interactions. While genomes of several extracellular sponge symbionts have been published, the intracellular symbionts have, so far, been elusive. Here we compare the genomes of two unicellular cyanobacterial sponge symbionts that share an ancestor but followed different evolutionary paths-one became intracellular and the other extracellular. Counterintuitively, the intracellular cyanobacteria are facultative, while the extracellular ones are obligate. By sequencing the genomes of the intracellular cyanobacteria and comparing them to the genomes of the extracellular symbionts and related free-living cyanobacteria, we show how three different cyanobacterial lifestyles are reflected by adaptive genomic features.

A 16S rDNA PCR-based theoretical to actual delta approach on culturable mock communities revealed severe losses of diversity information

BACKGROUND

Subunits of ribosomal RNA genes (rDNAs) characterized by PCR-based protocols have been the proxy for studies in microbial taxonomy, phylogenetics, evolution and ecology. However, relevant factors have shown to interfere in the experimental outputs in a variety of systems. In this work, a 'theoretical' to 'actual' delta approach was applied to data on culturable mock bacterial communities (MBCs) to study the levels of losses in operational taxonomic units (OTUs) detectability. Computational and lab-bench strategies based on 16S rDNA amplification by 799F and U1492R primers were employed, using a fingerprinting method with highly improved detectability of fragments as a case-study tool. MBCs were of two major types: in silico MBCs, assembled with database-retrieved sequences, and in vitro MBCs, with AluI digestions of PCR data generated from culturable endophytes isolated from cacao trees.

RESULTS

Interfering factors for the 16 s rDNA amplifications, such as the type of template, direct and nested PCR, proportion of chloroplast DNA from a tropical plant source (Virola officinalis), and biased-amplification by the primers resulted in altered bacterial 16S rDNA amplification, both on MBCs and V. officinalis leaf-extracted DNA. For the theoretical data, the maximum number of fragments for in silico and in vitro cuts were not significantly different from each other. Primers' preferences for certain sequences were detected, depending on the MBCs' composition prior to PCR. The results indicated overall losses from 2.3 up to 8.2 times in the number of OTUs detected from actual AluI digestions of MBCs when compared to in silico and in vitro theoretical data.

CONCLUSIONS

Due to all those effects, the final amplification profile of the bacterial community assembled was remarkably simplified when compared to the expected number of detectable fragments known to be present in the MBC. From these findings, the scope of hypotheses generation and conclusions from experiments based on PCR amplifications of bacterial communities was discussed.

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.

Cyanobacterial phylogenetic analysis based on phylogenomics approaches render evolutionary diversification and adaptation: an overview of representative orders

Phylogenetic studies based on a definite set of marker genes usually reconstruct evolutionary relationships among the prokaryotic species. Based on specific target sequences, such studies represent variations and allow identification of similarities or dissimilarities in organisms. With the advent of completely sequenced genomes and accumulation of information on whole prokaryotic genomes, phylogenetic reconstructions should be considered more reliable if they are ideally based on entire genomes to resolve phylogenetic interest. We applied phylogenomics approaches taking into account completely sequenced cyanobacterial genomes to reconstruct underlying species that represented major taxonomic classes and belonged to distinctly different habitats (freshwater, marine, soils, and rocks). We did not rely on describing phylogeny of all representative class of cyanobacterial species on the basis of only ribosomal gene, 16S rDNA gene. In contrast, we analyzed combined molecular marker and phylogenomics approaches (genome alignment, gene content and gene order, composition vector and protein domain content) for accurately inferring phylogenetic relationship of species. We have shown that this approach reflects the impact of evolution on the organisms and considers connects with the ecological adaptation in cyanobacteria in different habitats. Analysis revealed that the members from marine habitat occupy different profile than those from freshwater. Impact of GC content and genomic repetitiveness over the diversification of cyanobacterial species and their possible role in adaptation was also reflected. Members occupying similar habitats cover more evolutionary distance together and also evolve various strategies for adaptation and survival either through genomic repetitiveness or preferences for genes of particular functions or modified GC content. Genomes undergo different changes for their adaptation in diverse habitats.

Cycad Coralloid Roots Contain Bacterial Communities Including Cyanobacteria and Caulobacter spp. That Encode Niche-Specific Biosynthetic Gene Clusters

Cycads are the only early seed plants that have evolved a specialized root to host endophytic bacteria that fix nitrogen. To provide evolutionary and functional insights into this million-year old symbiosis, we investigate endophytic bacterial sub-communities isolated from coralloid roots of species from Dioon (Zamiaceae) sampled from their natural habitats. We employed a sub-community co-culture experimental strategy to reveal both predominant and rare bacteria, which were characterized using phylogenomics and detailed metabolic annotation. Diazotrophic plant endophytes, including Bradyrhizobium, Burkholderia, Mesorhizobium, Rhizobium, and Nostoc species, dominated the epiphyte-free sub-communities. Draft genomes of six cyanobacteria species were obtained after shotgun metagenomics of selected sub-communities. These data were used for whole-genome inferences that suggest two Dioon-specific monophyletic groups, and a level of specialization characteristic of co-evolved symbiotic relationships. Furthermore, the genomes of these cyanobacteria were found to encode unique biosynthetic gene clusters, predicted to direct the synthesis of specialized metabolites, mainly involving peptides. After combining genome mining with detection of pigment emissions using multiphoton excitation fluorescence microscopy, we also show that Caulobacter species co-exist with cyanobacteria, and may interact with them by means of a novel indigoidine-like specialized metabolite. We provide an unprecedented view of the composition of the cycad coralloid root, including phylogenetic and functional patterns mediated by specialized metabolites that may be important for the evolution of ancient symbiotic adaptations.

A New Quartet-Based Statistical Method for Comparing Sets of Gene Trees Is Developed Using a Generalized Hoeffding Inequality

Extracting the strength of the tree signal that is encompassed by a collection of gene trees is an exceptionally challenging problem in phylogenomics. Often, this problem not only involves the construction of individual phylogenies based on different genes, which may be a difficult endeavor on its own, but is also exacerbated by many factors that create conflicts between the evolutionary histories of different gene families, such as duplications or losses of genes; hybridization events; incomplete lineage sorting; and horizontal gene transfer, the latter two play central roles in the evolution of eukaryotes and prokaryotes, respectively. In this work, we tackle the aforementioned problem by focusing on quartet trees, which are the most basic unit of information in the context of unrooted phylogenies. In the first part, we show how a theorem of Janson that generalizes the classical Hoeffding inequality can be used to develop a statistical test involving quartets. In the second part, we study real and simulated data using this theoretical advancement, thus demonstrating how the significance of the differences between sets of quartets can be assessed. Our results are particularly intriguing since they nonstandardly require the analysis of dependent random variables.

A phylogenomic study quantifies competing mechanisms for pseudogenization in prokaryotes-The Mycobacterium leprae case

BACKGROUND

Pseudogenes are non-functional sequences in the genome with homologous sequences that are functional (i.e. genes). They are abundant in eukaryotes where they have been extensively investigated, while in prokaryotes they are significantly scarcer and less well studied. Here we conduct a comprehensive analysis of the evolution of orthologs of Mycobacterium leprae pseudogenes in prokaryotes. The leprosy pathogen M. leprae is of particular interest since it contains an unusually large number of pseudogenes, comprising approximately 40% of its entire genome. The analysis is conducted in both broad and narrow phylogenetic ranges.

RESULTS

We have developed an informatics-based approach to characterize the evolution of pseudogenes. This approach combines tools from phylogenomics, genomics, and transcriptomics. The results we obtain are used to assess the contributions of two mechanisms for pseudogene formation: failed horizontal gene transfer events and disruption of native genes.

CONCLUSIONS

We conclude that, although it was reported that in most bacteria the former is most likely responsible for the majority of pseudogenization events, in mycobacteria, and in particular in M. leprae with its exceptionally high pseudogene numbers, the latter predominates. We believe that our study sheds new light on the evolution of pseudogenes in bacteria, by utilizing new methodologies that are applied to the unusually abundant M. leprae pseudogenes and their orthologs.

Reconstruction of real and simulated phylogenies based on quartet plurality inference

BACKGROUND

Deciphering the history of life on Earth has long been regarded as one of the most central tasks in biology. In past years, widespread discordance between the evolutionary histories of different groups of orthologous genes of prokaryotes have been revealed, primarily due to horizontal gene transfers (HGTs). Nonetheless, evidence that support a strong tree-like signal of evolution have been uncovered, despite the presence of HGT events. Therefore, a challenging task is to distill this tree-like signal from the noise induced by all sources of non-tree-like events.

RESULTS

In this work we tackle this question, using real and simulated data. We first tighten a recent related theoretical result in this field. In a simulation study, we infer individual quartet topologies, and then use the inferred quartets to reconstruct simulated species trees. We demonstrate that accurate tree reconstruction is feasible despite surprisingly high rates of HGT. In a real data study, we construct phylogenies of two sets of prokaryotes, and show that our tree reconstruction scheme is comparable with (and complementary better than) other commonly used methods.

CONCLUSIONS

Using a blend of theoretical and empirical investigations, our study proves the feasibility of accurate quartet-based phylogenetic reconstruction, the vast impact of HGT events notwithstanding.

Consensus assessment of the contamination level of publicly available cyanobacterial genomes

Publicly available genomes are crucial for phylogenetic and metagenomic studies, in which contaminating sequences can be the cause of major problems. This issue is expected to be especially important for Cyanobacteria because axenic strains are notoriously difficult to obtain and keep in culture. Yet, despite their great scientific interest, no data are currently available concerning the quality of publicly available cyanobacterial genomes. As reliably detecting contaminants is a complex task, we designed a pipeline combining six methods in a consensus strategy to assess the contamination level of 440 genome assemblies of Cyanobacteria. Two methods are based on published reference databases of ribosomal genes (SSU rRNA 16S and ribosomal proteins), one is indirectly based on a reference database of marker genes (CheckM), and three are based on complete genome analysis. Among those genome-wide methods, Kraken and DIAMOND blastx share the same reference database that we derived from Ensembl Bacteria, whereas CONCOCT does not require any reference database, instead relying on differences in DNA tetramer frequencies. Given that all the six methods appear to have their own strengths and limitations, we used the consensus of their rankings to infer that >5% of cyanobacterial genome assemblies are highly contaminated by foreign DNA (i.e., contaminants were detected by 5 or 6 methods). Our results will help researchers to check the quality of publicly available genomic data before use in their own analyses. Moreover, we argue that journals should make mandatory the submission of raw read data along with genome assemblies in order to facilitate the detection of contaminants in sequence databases.

Genomic Changes Associated with the Evolutionary Transitions of Nostoc to a Plant Symbiont

Cyanobacteria belonging to the genus Nostoc comprise free-living strains and also facultative plant symbionts. Symbiotic strains can enter into symbiosis with taxonomically diverse range of host plants. Little is known about genomic changes associated with evolutionary transition of Nostoc from free-living to plant symbiont. Here, we compared the genomes derived from 11 symbiotic Nostoc strains isolated from different host plants and infer phylogenetic relationships between strains. Phylogenetic reconstructions of 89 Nostocales showed that symbiotic Nostoc strains with a broad host range, entering epiphytic and intracellular or extracellular endophytic interactions, form a monophyletic clade indicating a common evolutionary history. A polyphyletic origin was found for Nostoc strains which enter only extracellular symbioses, and inference of transfer events implied that this trait was likely acquired several times in the evolution of the Nostocales. Symbiotic Nostoc strains showed enriched functions in transport and metabolism of organic sulfur, chemotaxis and motility, as well as the uptake of phosphate, branched-chain amino acids, and ammonium. The genomes of the intracellular clade differ from that of other Nostoc strains, with a gain/enrichment of genes encoding proteins to generate l-methionine from sulfite and pathways for the degradation of the plant metabolites vanillin and vanillate, and of the macromolecule xylan present in plant cell walls. These compounds could function as C-sources for members of the intracellular clade. Molecular clock analysis indicated that the intracellular clade emerged ca. 600 Ma, suggesting that intracellular Nostoc symbioses predate the origin of land plants and the emergence of their extant hosts.

Phylogeny and salt-tolerance of freshwater Nostocales strains: Contribution to their systematics and evolution

Phylogenetic relationships among heterocytous genera (the Nostocales order) have been profoundly modified since the use of polyphasic approaches that include molecular data. There is nonetheless still ample scope for improving phylogenetic delineations of genera with broad ecological distributions, particularly by integrating specimens from specific or up-to-now poorly sampled habitats. In this context, we studied 36 new isolates belonging to Chrysosporum, Dolichospermum, Anabaena, Anabaenopsis, and Cylindrospermopsis from freshwater ecosystems of Burkina-Faso, Senegal, and Mayotte Island. Studying strains from these habitats is of particular interest as we suspected different range of salt variations during underwent periods of drought in small ponds and lakes. Such salt variation may cause different adaptation to salinity. We then undertook a polyphasic approach, combining molecular phylogenies, morphological analyses, and physiological measurements of tolerance to salinity. Molecular phylogenies of 117 Nostocales sequences showed that the 36 studied strains were distributed in seven lineages: Dolichospermum, Chrysosporum, Cylindrospermopsis/Raphidiopsis, Anabaenopsis, Anabaena sphaerica var tenuis/Sphaerospermopsis, and two independent Anabaena sphaerica lineages. Physiological data were congruent with molecular results supporting the separation into seven lineages. In an evolutionary context, salinity tolerance can be used as an integrative marker to reinforce the delineation of some cyanobacterial lineages. The history of this physiological trait contributes to a better understanding of processes leading to the divergence of cyanobacteria. In this study, most of the cyanobacterial strains isolated from freshwater environments were salt-tolerant, thus suggesting this trait constituted an ancestral trait of the heterocytous cyanobacteria and that it was probably lost two times secondarily and independently in the ancestor of Dolichospermum and of Cylindrospermopsis.

The Performance of Two Supertree Schemes Compared Using Synthetic and Real Data Quartet Input

Despite impressive advancements in technological and theoretical tools, construction of phylogenetic (evolutionary) trees is still a challenging task. The availability of enormous quantities of molecular data has made large-scale phylogenetic reconstruction involving thousands of species, a more viable goal. For this goal, separate trees over different, overlapping subsets of species, representing histories of various markers of these species, are collected. These trees, typically with conflicting signals, are subsequently combined into a single tree over the full set, an operation denoted as supertree construction. The amalgamation of such trees into a single tree lies at the heart of many tasks in phylogenetics, yet remains a daunting endeavor, especially in light of conflicting signals. In this work, we study the performance of matrix representation with parsimony (MRP), the most widely used supertree method to date, when confronted with quartet trees. Quartet trees are the most basic informational unit when amalgamation of unrooted trees is attempted, and they remain relevant in more general settings even though standard supertree methods are not necessarily confined to quartets. This study involves both real and simulated data, and the effects of several parameters on the results are evaluated, revealing a number of anomalies associated with MRP. We show that these anomalies are surmountable when using a recently introduced supertree method, weighted quartet MaxCut (wQMC).

Toxic genes present a unique phylogenetic signature

Horizontal gene transfer (HGT) is a major part of the evolution of Archaea and Bacteria, to the extent that the validity of the Tree of Life concept for prokaryotes has been seriously questioned. The patterns and routes of HGT remain a subject of intense study and debate. It was discovered that while several genes exhibit rampant HGT across the whole prokaryotic tree of life, others are lethal to certain organisms and therefore cannot be successfully transferred to them. We distinguish between these two classes of genes and show analytically that genes found to be toxic to a specific species (E. coli) also resist HGT in general. Several tools we employ show evidence to support that claim. One of those tools is the quartet plurality distribution (QPD), a mathematical tool that measures tendency to HGT over a large set of genes and species. When aggregated over a collection of genes, it can reveal important properties of this collection. We conclude that evidence of toxicity of certain genes to a wide variety of prokaryotes are revealed using the new tool of quartet plurality distribution.

Recurrent horizontal transfer of arsenite methyltransferase genes facilitated adaptation of life to arsenic

The toxic metalloid arsenic has been environmentally ubiquitous since life first arose nearly four billion years ago and presents a challenge for the survival of all living organisms. Its bioavailability has varied dramatically over the history of life on Earth. As life spread, biogeochemical and climate changes cyclically increased and decreased bioavailable arsenic. To elucidate the history of arsenic adaptation across the tree of life, we reconstructed the phylogeny of the arsM gene that encodes the As(III) S-adenosylmethionine (SAM) methyltransferase. Our results suggest that life successfully moved into arsenic-rich environments in the late Archean Eon and Proterozoic Eon, respectively, by the spread of arsM genes. The arsM genes of bacterial origin have been transferred to other kingdoms of life on at least six occasions, and the resulting domesticated arsM genes promoted adaptation to environmental arsenic. These results allow us to peer into the history of arsenic adaptation of life on our planet and imply that dissemination of genes encoding diverse adaptive functions to toxic chemicals permit adaptation to changes in concentrations of environmental toxins over evolutionary history.

A Metagenomic Approach to Cyanobacterial Genomics

Cyanobacteria, or oxyphotobacteria, are primary producers that establish ecological interactions with a wide variety of organisms. Although their associations with eukaryotes have received most attention, interactions with bacterial and archaeal symbionts have also been occurring for billions of years. Due to these associations, obtaining axenic cultures of cyanobacteria is usually difficult, and most isolation efforts result in unicyanobacterial cultures containing a number of associated microbes, hence composing a microbial consortium. With rising numbers of cyanobacterial blooms due to climate change, demand for genomic evaluations of these microorganisms is increasing. However, standard genomic techniques call for the sequencing of axenic cultures, an approach that not only adds months or even years for culture purification, but also appears to be impossible for some cyanobacteria, which is reflected in the relatively low number of publicly available genomic sequences of this phylum. Under the framework of metagenomics, on the other hand, cumbersome techniques for achieving axenic growth can be circumvented and individual genomes can be successfully obtained from microbial consortia. This review focuses on approaches for the genomic and metagenomic assessment of non-axenic cyanobacterial cultures that bypass requirements for axenity. These methods enable researchers to achieve faster and less costly genomic characterizations of cyanobacterial strains and raise additional information about their associated microorganisms. While non-axenic cultures may have been previously frowned upon in cyanobacteriology, latest advancements in metagenomics have provided new possibilities for in vitro studies of oxyphotobacteria, renewing the value of microbial consortia as a reliable and functional resource for the rapid assessment of bloom-forming cyanobacteria.

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.

Fast algorithm for the reconciliation of gene trees and LGT networks

In phylogenomics, reconciliations aim at explaining the discrepancies between the evolutionary histories of genes and species. Several reconciliation models are available when the evolution of the species of interest is modelled via phylogenetic trees; the most commonly used are the DL model, accounting for duplications and losses in gene evolution and yielding polynomially-solvable problems, and the DTL model, which also accounts for gene transfers and implies NP-hard problems. However, when dealing with non-tree-like evolutionary events such as hybridisations, phylogenetic networks - and not phylogenetic trees - should be used to model species evolution. Reconciliation models involving phylogenetic networks are still at their early days. In this paper, we propose a new reconciliation model in which the evolution of species is modelled by a special kind of phylogenetic networks - the LGT networks. Our model considers duplications, losses and transfers of genes, but restricts transfers to happen through some specific arcs of the network, called secondary arcs. Moreover, we provide a polynomial algorithm to compute the most parsimonious reconciliation between a gene tree and an LGT network under this model. Our method, when combined with quartet decomposition methods to detect putative "highways" of transfers, permits to refine their analyses by allowing to examine the two possible directions of a highway and even consider combinations of highways.

A comprehensive phylogenetic analysis of copper transporting P1B ATPases from bacteria of the Rhizobiales order uncovers multiplicity, diversity and novel taxonomic subtypes

The ubiquitous cytoplasmic membrane copper transporting P_1B‐1 and P_1B‐3‐type ATPases pump out Cu⁺ and Cu²⁺, respectively, to prevent cytoplasmic accumulation and avoid toxicity. The presence of five copies of Cu‐ATPases in the symbiotic nitrogen‐fixing bacteria Sinorhizobium meliloti is remarkable; it is the largest number of Cu⁺‐transporters in a bacterial genome reported to date. Since the prevalence of multiple Cu‐ATPases in members of the Rhizobiales order is unknown, we performed an in silico analysis to understand the occurrence, diversity and evolution of Cu⁺‐ATPases in members of the Rhizobiales order. Multiple copies of Cu‐ATPase coding genes (2–8) were detected in 45 of the 53 analyzed genomes. The diversity inferred from a maximum‐likelihood (ML) phylogenetic analysis classified Cu‐ATPases into four monophyletic groups. Each group contained additional subtypes, based on the presence of conserved motifs. This novel phylogeny redefines the current classification, where they are divided into two subtypes (P_1B‐1 and P_1B‐3). Horizontal gene transfer (HGT) as well as the evolutionary dynamic of plasmid‐borne genes may have played an important role in the functional diversification of Cu‐ATPases. Homologous cytoplasmic and periplasmic Cu⁺‐chaperones, CopZ, and CusF, that integrate a CopZ‐CopA‐CusF tripartite efflux system in gamma‐proteobacteria and archeae, were found in 19 of the 53 surveyed genomes of the Rhizobiales. This result strongly suggests a high divergence of CopZ and CusF homologs, or the existence of unexplored proteins involved in cellular copper transport.

Robust Inference of Genetic Exchange Communities from Microbial Genomes Using TF-IDF

Bacteria and archaea can exchange genetic material across lineages through processes of lateral genetic transfer (LGT). Collectively, these exchange relationships can be modeled as a network and analyzed using concepts from graph theory. In particular, densely connected regions within an LGT network have been defined as genetic exchange communities (GECs). However, it has been problematic to construct networks in which edges solely represent LGT. Here we apply term frequency-inverse document frequency (TF-IDF), an alignment-free method originating from document analysis, to infer regions of lateral origin in bacterial genomes. We examine four empirical datasets of different size (number of genomes) and phyletic breadth, varying a key parameter (word length k) within bounds established in previous work. We map the inferred lateral regions to genes in recipient genomes, and construct networks in which the nodes are groups of genomes, and the edges natively represent LGT. We then extract maximum and maximal cliques (i.e., GECs) from these graphs, and identify nodes that belong to GECs across a wide range of k. Most surviving lateral transfer has happened within these GECs. Using Gene Ontology enrichment tests we demonstrate that biological processes associated with metabolism, regulation and transport are often over-represented among the genes affected by LGT within these communities. These enrichments are largely robust to change of k.