Count: evolutionary analysis of phylogenetic profiles with parsimony and likelihood

From Root to Tips: Sporulation Evolution and Specialization in Bacillus subtilis and the Intestinal Pathogen Clostridioides difficile

Bacteria of the Firmicutes phylum are able to enter a developmental pathway that culminates with the formation of highly resistant, dormant endospores. Endospores allow environmental persistence, dissemination and for pathogens, are also infection vehicles. In both the model Bacillus subtilis, an aerobic organism, and in the intestinal pathogen Clostridioides difficile, an obligate anaerobe, sporulation mobilizes hundreds of genes. Their expression is coordinated between the forespore and the mother cell, the two cells that participate in the process, and is kept in close register with the course of morphogenesis. The evolutionary mechanisms by which sporulation emerged and evolved in these two species, and more broadly across Firmicutes, remain largely unknown. Here, we trace the origin and evolution of sporulation using the genes known to be involved in the process in B. subtilis and C. difficile, and estimating their gain-loss dynamics in a comprehensive bacterial macroevolutionary framework. We show that sporulation evolution was driven by two major gene gain events, the first at the base of the Firmicutes and the second at the base of the B. subtilis group and within the Peptostreptococcaceae family, which includes C. difficile. We also show that early and late sporulation regulons have been coevolving and that sporulation genes entail greater innovation in B. subtilis with many Bacilli lineage-restricted genes. In contrast, C. difficile more often recruits new sporulation genes by horizontal gene transfer, which reflects both its highly mobile genome, the complexity of the gut microbiota, and an adjustment of sporulation to the gut ecosystem.

Genomic characterization of Parengyodontium americanum sp. nov

Modern genome analysis and phylogenomic methods have increased the number of fungal species, as well as enhanced appreciation of the degree of diversity within the fungal kingdom. In this context, we describe a new Parengyodontium species, P. americanum, which is phylogenetically related to the opportunistic human fungal pathogen P. album. Five unusual fungal isolates were recovered from five unique and confirmed coccidioidomycosis patients, and these isolates were subsequently submitted to detailed molecular and morphological identification procedures to determine identity. Molecular and morphological diagnostic analyses showed that the isolates belong to the Cordycipitaceae. Subsequently, three representative genomes were sequenced and annotated, and a new species, P. americanum, was identified. Using various genomic analyses, gene family expansions related to novel compounds and potential for ability to grow in diverse habitats are predicted. A general description of the genomic composition of this newly described species and comparison of genome content with Beauveria bassiana, Isaria fumosorosea and Cordyceps militaris shows a shared core genome of 6371 genes, and 148 genes that appear to be specific for P. americanum. This work provides the framework for future investigations of this interesting fungal species.

Phylogenomic Analyses of Members of the Widespread Marine Heterotrophic Genus Pseudovibrio Suggest Distinct Evolutionary Trajectories and a Novel Genus, Polycladidibacter gen. nov

Bacteria belonging to the Pseudovibrio genus are widespread, metabolically versatile, and able to thrive as both free-living and host-associated organisms. Although more than 50 genomes are available, a comprehensive comparative genomics study to resolve taxonomic inconsistencies is currently missing. We analyzed all available genomes and used 552 core genes to perform a robust phylogenomic reconstruction. This in-depth analysis revealed the divergence of two monophyletic basal lineages of strains isolated from polyclad flatworm hosts, namely, Pseudovibrio hongkongensis and Pseudovibrio stylochi These strains have reduced genomes and lack sulfur-related metabolisms and major biosynthetic gene clusters, and their environmental distribution appears to be tightly associated with invertebrate hosts. We showed experimentally that the divergent strains are unable to utilize various sulfur compounds that, in contrast, can be utilized by the type strain Pseudovibrio denitrificans Our analyses suggest that the lineage leading to these two strains has been subject to relaxed purifying selection resulting in great gene loss. Overall genome relatedness indices (OGRI) indicate substantial differences between the divergent strains and the rest of the genus. While 16S rRNA gene analyses do not support the establishment of a different genus for the divergent strains, their substantial genomic, phylogenomic, and physiological differences strongly suggest a divergent evolutionary trajectory and the need for their reclassification. Therefore, we propose the novel genus Polycladidibacter gen. nov.IMPORTANCE The genus Pseudovibrio is commonly associated with marine invertebrates, which are essential for ocean health and marine nutrient cycling. Traditionally, the phylogeny of the genus has been based on 16S rRNA gene analysis. The use of the 16S rRNA gene or any other single marker gene for robust phylogenetic placement has recently been questioned. We used a large set of marker genes from all available Pseudovibrio genomes for in-depth phylogenomic analyses. We identified divergent monophyletic basal lineages within the Pseudovibrio genus, including two strains isolated from polyclad flatworms. These strains showed reduced sulfur metabolism and biosynthesis capacities. The phylogenomic analyses revealed distinct evolutionary trajectories and ecological adaptations that differentiate the divergent strains from the other Pseudovibrio members and suggest that they fall into a novel genus. Our data show the importance of widening the use of phylogenomics for better understanding bacterial physiology, phylogeny, and evolution.

Genome Evolution of the Obligate Endosymbiont Buchnera aphidicola

An evolutionary consequence of uniparentally transmitted symbiosis is degradation of symbiont genomes. We use the system of aphids and their maternally inherited obligate endosymbiont, Buchnera aphidicola, to explore the evolutionary process of genome degradation. We compared complete genome sequences for 39 Buchnera strains, including 23 newly sequenced symbiont genomes from diverse aphid hosts. We reconstructed the genome of the most recent shared Buchnera ancestor, which contained 616 protein-coding genes, and 39 RNA genes. The extent of subsequent gene loss varied across lineages, resulting in modern genomes ranging from 412 to 646 kb and containing 354-587 protein-coding genes. Loss events were highly nonrandom across loci. Genes involved in replication, transcription, translation, and amino acid biosynthesis are largely retained, whereas genes underlying ornithine biosynthesis, stress responses, and transcriptional regulation were lost repeatedly. Aside from losses, gene order is almost completely stable. The main exceptions involve movement between plasmid and chromosome locations of genes underlying tryptophan and leucine biosynthesis and supporting nutrition of aphid hosts. This set of complete genomes enabled tests for signatures of positive diversifying selection. Of 371 Buchnera genes tested, 29 genes show strong support for ongoing positive selection. These include genes encoding outer membrane porins that are expected to be involved in direct interactions with hosts. Collectively, these results indicate that extensive genome reduction occurred in the ancestral Buchnera prior to aphid diversification and that reduction has continued since, with losses greater in some lineages and for some loci.

Comparative Genomic Analysis of the Pheromone Receptor Class 1 Family (V1R) Reveals Extreme Complexity in Mouse Lemurs (Genus, Microcebus) and a Chromosomal Hotspot across Mammals

Sensory gene families are of special interest for both what they can tell us about molecular evolution and what they imply as mediators of social communication. The vomeronasal type-1 receptors (V1Rs) have often been hypothesized as playing a fundamental role in driving or maintaining species boundaries given their likely function as mediators of intraspecific mate choice, particularly in nocturnal mammals. Here, we employ a comparative genomic approach for revealing patterns of V1R evolution within primates, with a special focus on the small-bodied nocturnal mouse and dwarf lemurs of Madagascar (genera Microcebus and Cheirogaleus, respectively). By doubling the existing genomic resources for strepsirrhine primates (i.e. the lemurs and lorises), we find that the highly speciose and morphologically cryptic mouse lemurs have experienced an elaborate proliferation of V1Rs that we argue is functionally related to their capacity for rapid lineage diversification. Contrary to a previous study that found equivalent degrees of V1R diversity in diurnal and nocturnal lemurs, our study finds a strong correlation between nocturnality and V1R elaboration, with nocturnal lemurs showing elaborate V1R repertoires and diurnal lemurs showing less diverse repertoires. Recognized subfamilies among V1Rs show unique signatures of diversifying positive selection, as might be expected if they have each evolved to respond to specific stimuli. Furthermore, a detailed syntenic comparison of mouse lemurs with mouse (genus Mus) and other mammalian outgroups shows that orthologous mammalian subfamilies, predicted to be of ancient origin, tend to cluster in a densely populated region across syntenic chromosomes that we refer to as a V1R "hotspot."

Comparative Pangenomics of the Mammalian Gut Commensal Bifidobacterium longum

Bifidobacterium longum colonizes mammalian gastrointestinal tracts where it could metabolize host-indigestible oligosaccharides. Although B. longum strains are currently segregated into three subspecies that reflect common metabolic capacities and genetic similarity, heterogeneity within subspecies suggests that these taxonomic boundaries may not be completely resolved. To address this, the B. longum pangenome was analyzed from representative strains isolated from a diverse set of sources. As a result, the B. longum pangenome is open and contains almost 17,000 genes, with over 85% of genes found in ≤28 of 191 strains. B. longum genomes share a small core gene set of only ~500 genes, or ~3% of the total pangenome. Although the individual B. longum subspecies pangenomes share similar relative abundances of clusters of orthologous groups, strains show inter- and intrasubspecies differences with respect to carbohydrate utilization gene content and growth phenotypes.

Genome analysis of sponge symbiont 'Candidatus Halichondribacter symbioticus' shows genomic adaptation to a host-dependent lifestyle

The marine sponge Halichondria panicea inhabits coastal areas around the globe and is a widely studied sponge species in terms of its biology, yet the ecological functions of its dominant bacterial symbiont 'Candidatus Halichondribacter symbioticus' remain unknown. Here, we present the draft genome of 'Ca. H. symbioticus' HS1 (2.8 Mbp, ca. 87.6% genome coverage) recovered from the sponge metagenome of H. panicea in order to study functions and symbiotic interactions at the genome level. Functional genome comparison of HS1 against closely related free-living seawater bacteria revealed a reduction of genes associated with carbohydrate transport and transcription regulation, pointing towards a limited carbohydrate metabolism, and static transcriptional dynamics reminiscent of other bacterial symbionts. In addition, HS1 was enriched in sponge symbiont specific gene families related to host-symbiont interactions and defence. Similarity in the functional gene repertoire between HS1 and a phylogenetically more distant symbiont in the marine sponge Aplysina aerophoba, based on COG category distribution, suggest a convergent evolution of symbiont specific traits and general metabolic features. This warrants further investigation into convergent genomic evolution of symbionts across different sponge species and habitats.

Glutamicibacter mishrai sp. nov., isolated from the coral Favia veroni from Andaman Sea

A novel aerobic marine actinobacterium (strain S5-52^T) belonging to the genus Glutamicibacter was isolated from the coral Favia veroni sampled from the Andaman Sea, India. Cells are Gram stain positive and rod shaped. The DNA G+C content was 58.7 mol%. The major quinones were MK-8 and MK-9. The polar lipids were diphosphatidylglycerol, phosphatidylglycerol, glycolipid, trimannosyldiacylglycerol, phospholipid and dimannosylglyceride. The peptidoglycan type was A4α. Strain S5-52^T showed a maximum 16S rRNA similarity of 99.36% with Glutamicibacter halophytocola DSM 101718^T. The genome of strain S5-52^T was 3.57 Mb that contains 3274 protein coding sequences (CDS). DNA-DNA similarity and ANI values between S5-52^T and the reference strains were below 70% and 95-96%, respectively. Analysis of genomic reduction events in the evolutionary path from the LUCA (last universal common ancestor) to G. mishrai LMG 29155^T and G. halophytocola DSM 101718^T exhibit a number of genes involved in amino acid metabolism, cell wall biogenesis and replication, recombination and repair mechanism that reduced in both the species. Based on phenotypic, chemotaxonomic properties and comparative genomic studies, the strain S5-52^T is considered a novel species of the genus Glutamicibacter, for which the name Glutamicibacter mishrai sp. nov. is proposed. The type strain is S5-52^T (= KCTC 39846^T = LMG 29155^T).

Genome-wide Transcriptional Analysis of Tetrahymena thermophila Response to Exogenous Cholesterol

The ciliate Tetrahymena thermophila does not require sterols for growth and synthesizes pentacyclic triterpenoid alcohols, mainly tetrahymanol, as sterol surrogates. However, when sterols are present in the environment, T. thermophila efficiently incorporates and modifies them. These modifications consist of desaturation reactions at positions C5(6), C7(8), and C22(23), and de-ethylation at C24 of 29-carbon sterols (i.e. phytosterols). Three out of four of the enzymes involved in the sterol modification pathway have been previously identified. However, identification of the sterol C22 desaturase remained elusive, as did other basic aspects of this metabolism. To get more insights into this peculiar metabolism, we here perform a whole transcriptome analysis of T. thermophila in response to exogenous cholesterol. We found 356 T. thermophila genes to be differentially expressed after supplementation with cholesterol for 2 h. Among those that were upregulated, we found two genes belonging to the long spacing family of desaturases that we tentatively identified by RNAi analysis as sterol C22 desaturases. Additionally, we determined that the inhibition of tetrahymanol synthesis after supplementation with cholesterol occurs by a transcriptional downregulation of genes involved in squalene synthesis and cyclization. Finally, we identified several uncharacterized genes that are likely involved in sterols transport and signaling.

Epigenomic characterization of Clostridioides difficile finds a conserved DNA methyltransferase that mediates sporulation and pathogenesis

Clostridioides difficile is a leading cause of health care-associated infections. Although significant progress has been made in the understanding of its genome, the epigenome of C. difficile and its functional impact has not been systematically explored. Here, we performed a comprehensive DNA methylome analysis of C. difficile using 36 human isolates and observed great epigenomic diversity. We discovered an orphan DNA methyltransferase with a well-defined specificity whose corresponding gene is highly conserved across our dataset and in all ∼300 global C. difficile genomes examined. Inactivation of the methyltransferase gene negatively impacted sporulation, a key step in C. difficile disease transmission, consistently supported by multi-omics data, genetic experiments, and a mouse colonization model. Further experimental and transcriptomic analysis also suggested that epigenetic regulation is associated with cell length, biofilm formation, and host colonization. These findings provide a unique epigenetic dimension to characterize medically relevant biological processes in this critical pathogen. This work also provides a set of methods for comparative epigenomics and integrative analysis, which we expect to be broadly applicable to bacterial epigenomics studies.

Reconstruction of protein domain evolution using single-cell amplified genomes of uncultured choanoflagellates sheds light on the origin of animals

Understanding the origins of animal multicellularity is a fundamental biological question. Recent genome data have unravelled the role that co-option of pre-existing genes played in the origin of animals. However, there were also some important genetic novelties at the onset of Metazoa. To have a clear understanding of the specific genetic innovations and how they appeared, we need the broadest taxon sampling possible, especially among early-branching animals and their unicellular relatives. Here, we take advantage of single-cell genomics to expand our understanding of the genomic diversity of choanoflagellates, the sister-group to animals. With these genomes, we have performed an updated and taxon-rich reconstruction of protein evolution from the Last Eukaryotic Common Ancestor (LECA) to animals. Our novel data re-defines the origin of some genes previously thought to be metazoan-specific, like the POU transcription factor, which we show appeared earlier in evolution. Moreover, our data indicate that the acquisition of new genes at the stem of Metazoa was mainly driven by duplications and protein domain rearrangement processes at the stem of Metazoa. Furthermore, our analysis allowed us to reveal protein domains that are essential to the maintenance of animal multicellularity. Our analyses also demonstrate the utility of single-cell genomics from uncultured taxa to address evolutionary questions. This article is part of a discussion meeting issue 'Single cell ecology'.

Metabolic diversification of nitrogen-containing metabolites by the expression of a heterologous lysine decarboxylase gene in Arabidopsis

Lysine decarboxylase converts l-lysine to cadaverine as a branching point for the biosynthesis of plant Lys-derived alkaloids. Although cadaverine contributes towards the biosynthesis of Lys-derived alkaloids, its catabolism, including metabolic intermediates and the enzymes involved, is not known. Here, we generated transgenic Arabidopsis lines by expressing an exogenous lysine/ornithine decarboxylase gene from Lupinus angustifolius (La-L/ODC) and identified cadaverine-derived metabolites as the products of the emerged biosynthetic pathway. Through untargeted metabolic profiling, we observed the upregulation of polyamine metabolism, phenylpropanoid biosynthesis and the biosynthesis of several Lys-derived alkaloids in the transgenic lines. Moreover, we found several cadaverine-derived metabolites specifically detected in the transgenic lines compared with the non-transformed control. Among these, three specific metabolites were identified and confirmed as 5-aminopentanal, 5-aminopentanoate and δ-valerolactam. Cadaverine catabolism in a representative transgenic line (DC29) was traced by feeding stable isotope-labeled [α-¹⁵ N]- or [ε-¹⁵ N]-l-lysine. Our results show similar ¹⁵ N incorporation ratios from both isotopomers for the specific metabolite features identified, indicating that these metabolites were synthesized via the symmetric structure of cadaverine. We propose biosynthetic pathways for the metabolites on the basis of metabolite chemistry and enzymes known or identified through catalyzing specific biochemical reactions in this study. Our study shows that this pool of enzymes with promiscuous activities is the driving force for metabolite diversification in plants. Thus, this study not only provides valuable information for understanding the catabolic mechanism of cadaverine but also demonstrates that cadaverine accumulation is one of the factors to expand plant chemodiversity, which may lead to the emergence of Lys-derived alkaloid biosynthesis.

A unicellular relative of animals generates a layer of polarized cells by actomyosin-dependent cellularization

In animals, cellularization of a coenocyte is a specialized form of cytokinesis that results in the formation of a polarized epithelium during early embryonic development. It is characterized by coordinated assembly of an actomyosin network, which drives inward membrane invaginations. However, whether coordinated cellularization driven by membrane invagination exists outside animals is not known. To that end, we investigate cellularization in the ichthyosporean Sphaeroforma arctica, a close unicellular relative of animals. We show that the process of cellularization involves coordinated inward plasma membrane invaginations dependent on an actomyosin network and reveal the temporal order of its assembly. This leads to the formation of a polarized layer of cells resembling an epithelium. We show that this stage is associated with tightly regulated transcriptional activation of genes involved in cell adhesion. Hereby we demonstrate the presence of a self-organized, clonally-generated, polarized layer of cells in a unicellular relative of animals.

Comparative genomics of Leishmania (Mundinia)

Background

Trypanosomatids of the genus Leishmania are parasites of mammals or reptiles transmitted by bloodsucking dipterans. Many species of these flagellates cause important human diseases with clinical symptoms ranging from skin sores to life-threatening damage of visceral organs. The genus Leishmania contains four subgenera: Leishmania, Sauroleishmania, Viannia, and Mundinia. The last subgenus has been established recently and remains understudied, although Mundinia contains human-infecting species. In addition, it is interesting from the evolutionary viewpoint, representing the earliest branch within the genus and possibly with a different type of vector. Here we analyzed the genomes of L. (M.) martiniquensis, L. (M.) enriettii and L. (M.) macropodum to better understand the biology and evolution of these parasites.

Results

All three genomes analyzed were approximately of the same size (~ 30 Mb) and similar to that of L. (Sauroleishmania) tarentolae, but smaller than those of the members of subgenera Leishmania and Viannia, or the genus Endotrypanum (~ 32 Mb). This difference was explained by domination of gene losses over gains and contractions over expansions at the Mundinia node, although only a few of these genes could be identified. The analysis predicts significant changes in the Mundinia cell surface architecture, with the most important ones relating to losses of LPG-modifying side chain galactosyltransferases and arabinosyltransferases, as well as β-amastins. Among other important changes were gene family contractions for the oxygen-sensing adenylate cyclases and FYVE zinc finger-containing proteins.

Conclusions

We suggest that adaptation of Mundinia to different vectors and hosts has led to alternative host-parasite relationships and, thereby, made some proteins redundant. Thus, the evolution of genomes in the genus Leishmania and, in particular, in the subgenus Mundinia was mainly shaped by host (or vector) switches.

Insights into the ecological roles and evolution of methyl-coenzyme M reductase-containing hot spring Archaea

Several recent studies have shown the presence of genes for the key enzyme associated with archaeal methane/alkane metabolism, methyl-coenzyme M reductase (Mcr), in metagenome-assembled genomes (MAGs) divergent to existing archaeal lineages. Here, we study the mcr-containing archaeal MAGs from several hot springs, which reveal further expansion in the diversity of archaeal organisms performing methane/alkane metabolism. Significantly, an MAG basal to organisms from the phylum Thaumarchaeota that contains mcr genes, but not those for ammonia oxidation or aerobic metabolism, is identified. Together, our phylogenetic analyses and ancestral state reconstructions suggest a mostly vertical evolution of mcrABG genes among methanogens and methanotrophs, along with frequent horizontal gene transfer of mcr genes between alkanotrophs. Analysis of all mcr-containing archaeal MAGs/genomes suggests a hydrothermal origin for these microorganisms based on optimal growth temperature predictions. These results also suggest methane/alkane oxidation or methanogenesis at high temperature likely existed in a common archaeal ancestor.

Methane metabolism by some lineages of Archaea contributes to the cycling of carbon on Earth. Here, the authors show high diversity of methyl-coenzyme M reductase (Mcr), a key enzyme associated with archaeal methane/alkane metabolism, in hot spring Archaea, and investigate their ecological roles and evolution.

Introducing Trait Networks to Elucidate the Fluidity of Organismal Evolution Using Palaeontological Data

Explaining the evolution of animals requires ecological, developmental, paleontological, and phylogenetic considerations because organismal traits are affected by complex evolutionary processes. Modeling a plurality of processes, operating at distinct time-scales on potentially interdependent traits, can benefit from approaches that are complementary treatments to phylogenetics. Here, we developed an inclusive network approach, implemented in the command line software ComponentGrapher, and analyzed trait co-occurrence of rhinocerotoid mammals. We identified stable, unstable, and pivotal traits, as well as traits contributing to complexes, that may follow to a common developmental regulation, that point to an early implementation of the postcranial Bauplan among rhinocerotoids. Strikingly, most identified traits are highly dissociable, used repeatedly in distinct combinations and in different taxa, which usually do not form clades. Therefore, the genes encoding these traits are likely recruited into novel gene regulation networks during the course of evolution. Our evo-systemic framework, generalizable to other evolved organizations, supports a pluralistic modeling of organismal evolution, including trees and networks.

Multiple Independent Origins of Apicomplexan-Like Parasites

The apicomplexans are a group of obligate animal pathogens that include Plasmodium (malaria), Toxoplasma (toxoplasmosis), and Cryptosporidium (cryptosporidiosis) [1]. They are an extremely diverse and specious group but are nevertheless united by a distinctive suite of cytoskeletal and secretory structures related to infection, called the apical complex, which is used to recognize and gain entry into animal host cells. The apicomplexans are also known to have evolved from free-living photosynthetic ancestors and retain a relict plastid (the apicoplast), which is non-photosynthetic but houses a number of other essential metabolic pathways [2]. Their closest relatives include a mix of both photosynthetic algae (chromerids) and non-photosynthetic microbial predators (colpodellids) [3]. Genomic analyses of these free-living relatives have revealed a great deal about how the alga-parasite transition may have taken place, as well as origins of parasitism more generally [4]. Here, we show that, despite the surprisingly complex origin of apicomplexans from algae, this transition actually occurred at least three times independently. Using single-cell genomics and transcriptomics from diverse uncultivated parasites, we find that two genera previously classified within the Apicomplexa, Piridium and Platyproteum, form separately branching lineages in phylogenomic analyses. Both retain cryptic plastids with genomic and metabolic features convergent with apicomplexans. These findings suggest a predilection in this lineage for both the convergent loss of photosynthesis and transition to parasitism, resulting in multiple lineages of superficially similar animal parasites.

Blast Fungal Genomes Show Frequent Chromosomal Changes, Gene Gains and Losses, and Effector Gene Turnover

Pyricularia is a fungal genus comprising several pathogenic species causing the blast disease in monocots. Pyricularia oryzae, the best-known species, infects rice, wheat, finger millet, and other crops. As past comparative and population genomics studies mainly focused on isolates of P. oryzae, the genomes of the other Pyricularia species have not been well explored. In this study, we obtained a chromosomal-level genome assembly of the finger millet isolate P. oryzae MZ5-1-6 and also highly contiguous assemblies of Pyricularia sp. LS, P. grisea, and P. pennisetigena. The differences in the genomic content of repetitive DNA sequences could largely explain the variation in genome size among these new genomes. Moreover, we found extensive gene gains and losses and structural changes among Pyricularia genomes, including a large interchromosomal translocation. We searched for homologs of known blast effectors across fungal taxa and found that most avirulence effectors are specific to Pyricularia, whereas many other effectors share homologs with distant fungal taxa. In particular, we discovered a novel effector family with metalloprotease activity, distinct from the well-known AVR-Pita family. We predicted 751 gene families containing putative effectors in 7 Pyricularia genomes and found that 60 of them showed differential expression in the P. oryzae MZ5-1-6 transcriptomes obtained under experimental conditions mimicking the pathogen infection process. In summary, this study increased our understanding of the structural, functional, and evolutionary genomics of the blast pathogen and identified new potential effector genes, providing useful data for developing crops with durable resistance.

The Role of Homology and Orthology in the Phylogenomic Analysis of Metazoan Gene Content

Resolving the relationships of animals (Metazoa) is crucial to our understanding of the origin of key traits such as muscles, guts, and nerves. However, a broadly accepted metazoan consensus phylogeny has yet to emerge. In part, this is because the genomes of deeply diverging and fast-evolving lineages may undergo significant gene turnover, reducing the number of orthologs shared with related phyla. This can limit the usefulness of traditional phylogenetic methods that rely on alignments of orthologous sequences. Phylogenetic analysis of gene content has the potential to circumvent this orthology requirement, with binary presence/absence of homologous gene families representing a source of phylogenetically informative characters. Applying binary substitution models to the gene content of 26 complete animal genomes, we demonstrate that patterns of gene conservation differ markedly depending on whether gene families are defined by orthology or homology, that is, whether paralogs are excluded or included. We conclude that the placement of some deeply diverging lineages may exceed the limit of resolution afforded by the current methods based on comparisons of orthologous protein sequences, and novel approaches are required to fully capture the evolutionary signal from genes within genomes.

Evolution of the Drosophila melanogaster Chromatin Landscape and Its Associated Proteins

In the nucleus of eukaryotic cells, genomic DNA associates with numerous protein complexes and RNAs, forming the chromatin landscape. Through a genome-wide study of chromatin-associated proteins in Drosophila cells, five major chromatin types were identified as a refinement of the traditional binary division into hetero- and euchromatin. These five types were given color names in reference to the Greek word chroma. They are defined by distinct but overlapping combinations of proteins and differ in biological and biochemical properties, including transcriptional activity, replication timing, and histone modifications. In this work, we assess the evolutionary relationships of chromatin-associated proteins and present an integrated view of the evolution and conservation of the fruit fly Drosophila melanogaster chromatin landscape. We combine homology prediction across a wide range of species with gene age inference methods to determine the origin of each chromatin-associated protein. This provides insight into the evolution of the different chromatin types. Our results indicate that for the euchromatic types, YELLOW and RED, young associated proteins are more specialized than old ones; and for genes found in either chromatin type, intron/exon structure is lineage-specific. Next, we provide evidence that a subset of GREEN-associated proteins is involved in a centromere drive in D. melanogaster. Our results on BLUE chromatin support the hypothesis that the emergence of Polycomb Group proteins is linked to eukaryotic multicellularity. In light of these results, we discuss how the regulatory complexification of chromatin links to the origins of eukaryotic multicellularity.

Genomic evidence of the illumination response mechanism and evolutionary history of magnetotactic bacteria within the Rhodospirillaceae family

BACKGROUND

Magnetotactic bacteria (MTB) are ubiquitous in natural aquatic environments. MTB can produce intracellular magnetic particles, navigate along geomagnetic field, and respond to light. However, the potential mechanism by which MTB respond to illumination and their evolutionary relationship with photosynthetic bacteria remain elusive.

RESULTS

We utilized genomes of the well-sequenced genus Magnetospirillum, including the newly sequenced MTB strain Magnetospirillum sp. XM-1 to perform a comprehensive genomic comparison with phototrophic bacteria within the family Rhodospirillaceae regarding the illumination response mechanism. First, photoreceptor genes were identified in the genomes of both MTB and phototrophic bacteria in the Rhodospirillaceae family, but no photosynthesis genes were found in the MTB genomes. Most of the photoreceptor genes in the MTB genomes from this family encode phytochrome-domain photoreceptors that likely induce red/far-red light phototaxis. Second, illumination also causes damage within the cell, and in Rhodospirillaceae, both MTB and phototrophic bacteria possess complex but similar sets of response and repair genes, such as oxidative stress response, iron homeostasis and DNA repair system genes. Lastly, phylogenomic analysis showed that MTB cluster closely with phototrophic bacteria in this family. One photoheterotrophic genus, Phaeospirillum, clustered within and displays high genomic similarity with Magnetospirillum. Moreover, the phylogenetic tree topologies of magnetosome synthesis genes in MTB and photosynthesis genes in phototrophic bacteria from the Rhodospirillaceae family were reasonably congruent with the phylogenomic tree, suggesting that these two traits were most likely vertically transferred during the evolution of their lineages.

CONCLUSION

Our new genomic data indicate that MTB and phototrophic bacteria within the family Rhodospirillaceae possess diversified photoreceptors that may be responsible for phototaxis. Their genomes also contain comprehensive stress response genes to mediate the negative effects caused by illumination. Based on phylogenetic studies, most of MTB and phototrophic bacteria in the Rhodospirillaceae family evolved vertically with magnetosome synthesis and photosynthesis genes. The ancestor of Rhodospirillaceae was likely a magnetotactic phototrophic bacteria, however, gain or loss of magnetotaxis and phototrophic abilities might have occurred during the evolution of ancestral Rhodospirillaceae lineages.

Microevolution and Adaptive Strategy of Psychrophilic Species Flavobacterium bomense sp. nov. Isolated From Glaciers

Numerous mountain glaciers located on the Tibetan Plateau are inhabited by abundant microorganisms. The microorganisms on the glacier surface are exposed to the cold, barren, and high-ultraviolet radiation environments. Although the microbial community composition on glaciers has been revealed by high-throughput sequencing, little is known about the microevolution and adaptive strategy of certain bacterial populations. In this study, we used a polyphasic approach to determine the taxonomic status of 11 psychrophilic Flavobacterium strains isolated from glaciers on the Tibetan Plateau and performed a comparative genomic analysis. The phylogenetic tree based on the concatenated single-copy gene sequences showed the 11 strains clustered together, forming a distinct and novel clade in the genus Flavobacterium. The average nucleotide identity (ANI) values among these strains were higher than 96%. However, the values much lower than 90% between them and related species indicated that they represent a novel species and the name Flavobacterium bomense sp. nov. is proposed. The core and accessory genomes of strains in this new Flavobacterium species showed diverse distinct patterns of gene content and metabolism pathway. In order to infer the driving evolutionary forces of the core genomes, homologous recombination was found to contribute twice as much to nucleotide substitutions as mutations. A series of genes encoding proteins with known or predicted roles in cold adaptation were found in their genomes, for example, cold-shock protein, proteorhodopsin, osmoprotection, and membrane-related proteins. A comparative analysis of the group with optimal growth temperature (OGT) ≤ 20°C and the group with OGT > 20°C of the 32 Flavobacterium type strains and 11 new strains revealed multiple amino acid substitutions, including the decrease of the proline and glutamine content and the increase of the methionine and isoleucine content in the group with OGT ≤ 20°C, which may contribute to increased protein flexibility at low temperatures. Thus, this study discovered a novel Flavobacterium species in glaciers, which has high intraspecific diversity and multiple adaptation mechanisms that enable them to cope and thrive in extreme habitats.

Genome-centric metagenomics resolves microbial diversity and prevalent truncated denitrification pathways in a denitrifying PAO-enriched bioprocess

Denitrification is the stepwise microbial reduction of nitrate or nitrite (NO₂^-) to nitrogen gas via the obligate intermediates nitric oxide (NO) and nitrous oxide (N₂O). Substantial N₂O accumulation has been reported in denitrifying enhanced biological phosphorus removal (EBPR) bioreactors enriched in denitrifying polyphosphate accumulating organisms (DPAOs), but little is known about underlying mechanisms for N₂O generation, prevalence of complete versus truncated denitrification pathways, or the impact of NO₂^- feed on DPAO-enriched consortia. To address this knowledge gap, we employed genome-resolved metagenomics to investigate nitrogen transformation potential in a NO₂^- fed denitrifying EBPR bioreactor enriched in Candidatus Accumulibacter and prone to N₂O accumulation. Our analysis yielded 41 near-complete metagenome-assembled genomes (MAGs), including two co-occurring Accumulibacter strains affiliated with clades IA and IC (the first published genome from this clade) and 39 non-PAO flanking bacterial genomes. The dominant Accumulibacter clade IA encoded genes for complete denitrification, while the lower abundance Accumulibacter clade IC harbored all denitrification genes except for a canonical respiratory NO reductase. Analysis of the 39 non-PAO MAGs revealed a high prevalence of taxa harboring an incomplete denitrification pathway. Of the 27 MAGs harboring capacity for at least one step in the denitrification pathway, 10 were putative N₂O producers lacking N₂O reductase, 16 were putative N₂O reducers that lacked at least one upstream denitrification gene, and only one harbored a complete denitrification pathway. We also documented increasing abundance over the course of reactor operation of putative N₂O producers. Our results suggest that the unusually high levels of N₂O production observed in this Accumulibacter-enriched consortium are linked in part to the selection for non-PAO flanking microorganisms with truncated denitrification pathways.

Rapid functional divergence after small-scale gene duplication in grasses

BACKGROUND

Gene duplication has played an important role in the evolution and domestication of flowering plants. Yet little is known about how plant duplicate genes evolve and are retained over long timescales, particularly those arising from small-scale duplication (SSD) rather than whole-genome duplication (WGD) events.

RESULTS

We address this question in the Poaceae (grass) family by analyzing gene expression data from nine tissues of Brachypodium distachyon, Oryza sativa japonica (rice), and Sorghum bicolor (sorghum). Consistent with theoretical predictions, expression profiles of most grass genes are conserved after SSD, suggesting that functional conservation is the primary outcome of SSD in grasses. However, we also uncover support for widespread functional divergence, much of which occurs asymmetrically via the process of neofunctionalization. Moreover, neofunctionalization preferentially targets younger (child) duplicate gene copies, is associated with RNA-mediated duplication, and occurs quickly after duplication. Further analysis reveals that functional divergence of SSD-derived genes is positively correlated with both sequence divergence and tissue specificity in all three grass species, and particularly with anther expression in B. distachyon.

CONCLUSIONS

Our results suggest that SSD-derived grass genes often undergo rapid functional divergence that may be driven by natural selection on male-specific phenotypes. These observations are consistent with those in several animal species, suggesting that duplicate genes take similar evolutionary trajectories in plants and animals.

Recombination and intraspecific polymorphism for the presence and absence of entire chromosomes in mitochondrial genomes

Although mitochondrial genomes are typically thought of as single circular molecules, these genomes are fragmented into multiple chromosomes in many eukaryotes, raising intriguing questions about inheritance and (in)stability of mtDNA in such systems. A previous comparison of mitochondrial genomes from two different individuals of the angiosperm species Silene noctiflora found variation in the presence of entire mitochondrial chromosomes. Here, we expand on this work with a geographically diverse sampling of 25 S. noctiflora populations and the closely related species S. turkestanica and S. undulata. Using a combination of deep sequencing and PCR-based screening for the presence of 22 different mitochondrial chromosomes, we found extensive variation in the complement of chromosomes across individuals. Much of this variation could be attributed to recent chromosome loss events, suggesting that the massively expanded and fragmented mitochondrial genomes of S. noctiflora may have entered a phase of genome reduction in which they are losing entire chromosomes at a rapid rate. Sequence analysis of mitochondrial and plastid genomes revealed genealogical differences both between these organelles and within the mitochondrial genome, indicating a history of recombination. Evidence that recombination has generated novel combinations of alleles was more frequent between loci on different mitochondrial chromosomes than it was within chromosomes. Therefore, the fragmentation of mitochondrial genomes and the assortment of chromosomes during mitochondrial inheritance appears to have contributed to a history of sexual-like recombination in the mtDNA of this species.

The phytopathogenic nature of Dickeya aquatica 174/2 and the dynamic early evolution of Dickeya pathogenicity

Dickeya is a genus of phytopathogenic enterobacterales causing soft rot in a variety of plants (e.g. potato, chicory, maize). Among the species affiliated to this genus, Dickeya aquatica, described in 2014, remained particularly mysterious because it had no known host. Furthermore, while D. aquatica was proposed to represent a deep-branching species among Dickeya genus, its precise phylogenetic position remained elusive. Here, we report the complete genome sequence of the D. aquatica type strain 174/2. We demonstrate the affinity of D. aquatica strain 174/2 for acidic fruits such as tomato and cucumber and show that exposure of this bacterium to acidic pH induces twitching motility. An in-depth phylogenomic analysis of all available Dickeya proteomes pinpoints D. aquatica as the second deepest branching lineage within this genus and reclassifies two lineages that likely correspond to new genomospecies (gs.): Dickeya gs. poaceaephila (Dickeya sp NCPPB 569) and Dickeya gs. undicola (Dickeya sp 2B12), together with a new putative genus, tentatively named Prodigiosinella. Finally, from comparative analyses of Dickeya proteomes, we infer the complex evolutionary history of this genus, paving the way to study the adaptive patterns and processes of Dickeya to different environmental niches and hosts. In particular, we hypothesize that the lack of xylanases and xylose degradation pathways in D. aquatica could reflect adaptation to aquatic charophyte hosts which, in contrast to land plants, do not contain xyloglucans.

Dissecting the Evolutionary Development of the Species Bifidobacterium animalis through Comparative Genomics Analyses

Bifidobacteria are members of the gut microbiota of animals, including mammals, birds, and social insects. In this study, we analyzed and determined the pangenome of Bifidobacterium animalis species, encompassing B. animalis subsp. animalis and the B. animalis subsp. lactis taxon, which is one of the most intensely exploited probiotic bifidobacterial species. In order to reveal differences within the B. animalis species, detailed comparative genomics and phylogenomics analyses were performed, indicating that these two subspecies recently arose through divergent evolutionary events. A subspecies-specific core genome was identified for both B. animalis subspecies, revealing the existence of subspecies-defining genes involved in carbohydrate metabolism. Notably, these in silico analyses coupled with carbohydrate profiling assays suggest genetic adaptations toward a distinct glycan milieu for each member of the B. animalis subspecies, resulting in a divergent evolutionary development of the two subspecies.IMPORTANCE The majority of characterized B. animalis strains have been isolated from human fecal samples. In order to explore genome variability within this species, we isolated 15 novel strains from the gastrointestinal tracts of different animals, including mammals and birds. The present study allowed us to reconstruct the pangenome of this taxon, including the genome contents of 56 B. animalis strains. Through careful assessment of subspecies-specific core genes of the B. animalis subsp. animalis/lactis taxon, we identified genes encoding enzymes involved in carbohydrate transport and metabolism, while unveiling specific gene acquisition and loss events that caused the evolutionary emergence of these two subspecies.

Population Genomics Insights into Adaptive Evolution and Ecological Differentiation in Streptomycetes

Deciphering the genomic variation that represents microevolutionary processes toward species divergence is key to understanding microbial speciation, which has long been under debate. Streptomycetes are filamentous bacteria that are ubiquitous in nature and the richest source of antibiotics; however, their speciation processes remain unknown. To tackle this issue, we performed a comprehensive population genomics analysis on Streptomyces albidoflavus residing in different habitats and with a worldwide distribution and identified and characterized the foundational changes within the species. We detected three well-defined phylogenomic clades, of which clades I and III mainly contained free-living (soil/marine) and insect-associated strains, respectively, and clade II had a mixed origin. By performing genome-wide association studies (GWAS), we identified a number of genetic variants associated with free-living or entomic (denoting or relating to insects) habitats in both the accessory and core genomes. These variants contributed collectively to the population structure and had annotated or confirmed functions that likely facilitate differential adaptation of the species. In addition, we detected higher levels of homologous recombination within each clade and in the free-living group than within the whole species and in the entomic group. A subset of the insect-associated strains (clade III) showed a relatively independent evolutionary trajectory with more symbiosis-favorable genes but little genetic interchange with the other lineages. Our results demonstrate that ecological adaptation promotes genetic differentiation in S. albidoflavus, suggesting a model of ecological speciation with gene flow in streptomycetes.IMPORTANCE Species are the fundamental units of ecology and evolution, and speciation leads to the astounding diversity of life on Earth. Studying speciation is thus of great significance to understand, protect, and exploit biodiversity, but it is a challenge in the microbial world. In this study, using population genomics, we placed Streptomyces albidoflavus strains in a spectrum of speciation and showed that the genetic differences between phylogenomic clusters evolved mainly by environmental selection and gene-specific sweeps. These findings highlight the role of ecology in structuring recombining bacterial species, making a step toward a deeper understanding of microbial speciation. Our results also raise concerns of an underrated microbial diversity at the intraspecies level, which can be utilized for mining of ecologically relevant natural products.

The 4-Celled Tetrabaena socialis Nuclear Genome Reveals the Essential Components for Genetic Control of Cell Number at the Origin of Multicellularity in the Volvocine Lineage

Multicellularity is the premier example of a major evolutionary transition in individuality and was a foundational event in the evolution of macroscopic biodiversity. The volvocine chlorophyte lineage is well suited for studying this process. Extant members span unicellular, simple colonial, and obligate multicellular taxa with germ-soma differentiation. Here, we report the nuclear genome sequence of one of the most morphologically simple organisms in this lineage-the 4-celled colonial Tetrabaena socialis and compare this to the three other complete volvocine nuclear genomes. Using conservative estimates of gene family expansions a minimal set of expanded gene families was identified that associate with the origin of multicellularity. These families are rich in genes related to developmental processes. A subset of these families is lineage specific, which suggests that at a genomic level the evolution of multicellularity also includes lineage-specific molecular developments. Multiple points of evidence associate modifications to the ubiquitin proteasomal pathway (UPP) with the beginning of coloniality. Genes undergoing positive or accelerating selection in the multicellular volvocines were found to be enriched in components of the UPP and gene families gained at the origin of multicellularity include components of the UPP. A defining feature of colonial/multicellular life cycles is the genetic control of cell number. The genomic data presented here, which includes diversification of cell cycle genes and modifications to the UPP, align the genetic components with the evolution of this trait.

Diversity and evolution of cytochromes P450 in stramenopiles

MAIN CONCLUSION

Comparative genomic analysis of cytochromes P450 revealed high diversification and dynamic changes in stramenopiles, associated with transcriptional responsiveness to various environmental stimuli. Comparative genomic and molecular evolution approaches were used to characterize cytochromes P450 (P450) diversity in stramenopiles. Phylogenetic analysis pointed to a high diversity of P450 in stramenopiles and identified three major clans. The CYP51 and CYP97 clans were present in brown algae, diatoms and Nannochloropsis gaditana, whereas the CYP5014 clan mainly includes oomycetes. Gene gain and loss patterns revealed that six CYP families-CYP51, CYP97, CYP5160, CYP5021, CYP5022, and CYP5165-predated the split of brown algae and diatoms. After they diverged, diatoms gained more CYP families, especially in the cold-adapted species Fragilariopsis cylindrus, in which eight new CYP families were found. Selection analysis revealed that the expanded CYP51 family in the brown alga Cladosiphon okamuranus exhibited a more relaxed selection constraint compared with those of other brown algae and diatoms. Our RNA-seq data further evidenced that most of P450s in Saccharina japonica are highly expressed in large sporophytes, which could potentially promote the large kelp formation in this developmental stage. A survey of Ectocarpus siliculosus and diatom transcriptomes showed that many P450s are responsive to stress, nutrient limitation or light quality, suggesting pivotal roles in detoxification or metabolic processes under adverse environmental conditions. The information provided in this study will be helpful in designing functional experiments and interpreting P450 roles in this particular lineage.

Evolutionary dynamics of the chromatophore genome in three photosynthetic Paulinella species

The thecate amoeba Paulinella is a valuable model for understanding plastid organellogenesis because this lineage has independently gained plastids (termed chromatophores) of alpha-cyanobacterial provenance. Plastid primary endosymbiosis in Paulinella occurred relatively recently (90–140 million years ago, Mya), whereas the origin of the canonical Archaeplastida plastid occurred >1,500 Mya. Therefore, these two events provide independent perspectives on plastid formation on vastly different timescales. Here we generated the complete chromatophore genome sequence from P. longichromatophora (979,356 bp, GC-content = 38.8%, 915 predicted genes) and P. micropora NZ27 (977,190 bp, GC-content = 39.9%, 911 predicted genes) and compared these data to that from existing chromatophore genomes. Our analysis suggests that when a basal split occurred among photosynthetic Paulinella species ca. 60 Mya, only 35% of the ancestral orthologous gene families from the cyanobacterial endosymbiont remained in chromatophore DNA. Following major gene losses during the early stages of endosymbiosis, this process slowed down significantly, resulting in a conserved gene content across extant taxa. Chromatophore genes faced relaxed selection when compared to homologs in free-living alpha-cyanobacteria, likely reflecting the homogeneous intracellular environment of the Paulinella host. Comparison of nucleotide substitution and insertion/deletion events among different P. micropora strains demonstrates that increases in AT-content and genome reduction are ongoing and dynamic processes in chromatophore evolution.

Evolution of substrate-specific gene expression and RNA editing in brown rot wood-decaying fungi

Fungi that decay wood have characteristic associations with certain tree species, but the mechanistic bases for these associations are poorly understood. We studied substrate-specific gene expression and RNA editing in six species of wood-decaying fungi from the 'Antrodia clade' (Polyporales, Agaricomycetes) on three different wood substrates (pine, spruce, and aspen) in submerged cultures. We identified dozens to hundreds of substrate-biased genes (i.e., genes that are significantly upregulated in one substrate relative to the other two substrates) in each species, and these biased genes are correlated with their host ranges. Evolution of substrate-biased genes is associated with gene family expansion, gain and loss of genes, and variation in cis- and trans- regulatory elements, rather than changes in protein coding sequences. We also demonstrated widespread RNA editing events in the Antrodia clade, which differ from those observed in the Ascomycota in their distribution, substitution types, and the genomic environment. Moreover, we found that substrates could affect editing positions and frequency, including editing events occurring in mRNA transcribed from wood-decay-related genes. This work shows the extent to which gene expression and RNA editing differ among species and substrates, and provides clues into mechanisms by which wood-decaying fungi may adapt to different hosts.

Pan-genome analyses of model fungal species

The concept of the species ‘pan-genome’, the union of ‘core’ conserved genes and all ‘accessory’ non-conserved genes across all strains of a species, was first proposed in prokaryotes to account for intraspecific variability. Species pan-genomes have been extensively studied in prokaryotes, but evidence of species pan-genomes has also been demonstrated in eukaryotes such as plants and fungi. Using a previously published methodology based on sequence homology and conserved microsynteny, in addition to bespoke pipelines, we have investigated the pan-genomes of four model fungal species: Saccharomyces cerevisiae, Candida albicans, Cryptococcus neoformans var. grubii and Aspergillus fumigatus. Between 80 and 90 % of gene models per strain in each of these species are core genes that are highly conserved across all strains of that species, many of which are involved in housekeeping and conserved survival processes. In many of these species, the remaining ‘accessory’ gene models are clustered within subterminal regions and may be involved in pathogenesis and antimicrobial resistance. Analysis of the ancestry of species core and accessory genomes suggests that fungal pan-genomes evolve by strain-level innovations such as gene duplication as opposed to wide-scale horizontal gene transfer. Our findings lend further supporting evidence to the existence of species pan-genomes in eukaryote taxa.

Trophic Specialization Results in Genomic Reduction in Free-Living Marine Idiomarina Bacteria

The streamlining hypothesis is usually used to explain the genomic reduction events in free-living bacteria like SAR11. However, we find that the genomic reduction phenomenon in the bacterial genus Idiomarina is different from that in SAR11. Therefore, we propose a new hypothesis to explain genomic reduction in this genus based on trophic specialization that could result in genomic reduction, which would be not uncommon in nature. Not only can the trophic specialization hypothesis explain the genomic reduction in the genus Idiomarina, but it also sheds new light on our understanding of the genomic reduction processes in other free-living bacterial lineages.

The streamlining hypothesis is generally used to explain the genomic reduction events related to the small genome size of free-living bacteria like marine bacteria SAR11. However, our current understanding of the correlation between bacterial genome size and environmental adaptation relies on too few species. It is still unclear whether there are other paths leading to genomic reduction in free-living bacteria. The genome size of marine free-living bacteria of the genus Idiomarina belonging to the order Alteromonadales (Gammaproteobacteria) is much smaller than the size of related genomes from bacteria in the same order. Comparative genomic and physiological analyses showed that the genomic reduction pattern in this genus is different from that of the classical SAR11 lineage. Genomic reduction reconstruction and substrate utilization profile showed that Idiomarina spp. lost a large number of genes related to carbohydrate utilization, and instead they specialized on using proteinaceous resources. Here we propose a new hypothesis to explain genomic reduction in this genus; we propose that trophic specialization increasing the metabolic efficiency for using one kind of substrate but reducing the substrate utilization spectrum could result in bacterial genomic reduction, which would be not uncommon in nature. This hypothesis was further tested in another free-living genus, Kangiella, which also shows dramatic genomic reduction. These findings highlight that trophic specialization is potentially an important path leading to genomic reduction in some marine free-living bacteria, which is distinct from the classical lineages like SAR11.

Insights into the Evolution of the New World Diploid Cottons (Gossypium, Subgenus Houzingenia) Based on Genome Sequencing

We employed phylogenomic methods to study molecular evolutionary processes and phylogeny in the geographically widely dispersed New World diploid cottons (Gossypium, subg. Houzingenia). Whole genome resequencing data (average of 33× genomic coverage) were generated to reassess the phylogenetic history of the subgenus and provide a temporal framework for its diversification. Phylogenetic analyses indicate that the subgenus likely originated following transoceanic dispersal from Africa about 6.6 Ma, but that nearly all of the biodiversity evolved following rapid diversification in the mid-Pleistocene (0.5-2.0 Ma), with multiple long-distance dispersals required to account for range expansion to Arizona, the Galapagos Islands, and Peru. Comparative analyses of cpDNAversus nuclear data indicate that this history was accompanied by several clear cases of interspecific introgression. Repetitive DNAs contribute roughly half of the total 880 Mb genome, but most transposable element families are relatively old and stable among species. In the genic fraction, pairwise synonymous mutation rates average 1% per Myr, with nonsynonymous changes being about seven times less frequent. Over 1.1 million indels were detected and phylogenetically polarized, revealing a 2-fold bias toward deletions over small insertions. We suggest that this genome down-sizing bias counteracts genome size growth by TE amplification and insertions, and helps explain the relatively small genomes that are restricted to this subgenus. Compared with the rate of nucleotide substitution, the rate of indel occurrence is much lower averaging about 17 nucleotide substitutions per indel event.

Evolution of the Large Nucleocytoplasmic DNA Viruses of Eukaryotes and Convergent Origins of Viral Gigantism

The Nucleocytoplasmic Large DNA Viruses (NCLDV) of eukaryotes (proposed order "Megavirales") comprise an expansive group of eukaryotic viruses that consists of the families Poxviridae, Asfarviridae, Iridoviridae, Ascoviridae, Phycodnaviridae, Marseilleviridae, Pithoviridae, and Mimiviridae, as well as Pandoraviruses, Molliviruses, and Faustoviruses that so far remain unaccounted by the official virus taxonomy. All these viruses have double-stranded DNA genomes that range in size from about 100 kilobases (kb) to more than 2.5 megabases. The viruses with genomes larger than 500kb are informally considered "giant," and the largest giant viruses surpass numerous bacteria and archaea in both particle and genome size. The discovery of giant viruses has been highly unexpected and has changed the perception of viral size and complexity, and even, arguably, the entire concept of a virus. Given that giant viruses encode multiple proteins that are universal among cellular life forms and are components of the translation system, the quintessential cellular molecular machinery, attempts have been made to incorporate these viruses in the evolutionary tree of cellular life. Moreover, evolutionary scenarios of the origin of giant viruses from a fourth, supposedly extinct domain of cellular life have been proposed. However, despite all the differences in the genome size and gene repertoire, the NCLDV can be confidently defined as monophyletic group, on the strength of the presence of about 40 genes that can be traced back to their last common ancestor. Using several most strongly conserved genes from this ancestral set, a well-resolved phylogenetic tree of the NCLDV was built and employed as the scaffold to reconstruct the history of gene gain and loss throughout the course of the evolution of this group of viruses. This reconstruction reveals extremely dynamic evolution that involved extensive gene gain and loss in many groups of viruses and indicates that giant viruses emerged independently in several clades of the NCLDV. Thus, these giants of the virus world evolved repeatedly from smaller and simpler viruses, rather than from a fourth domain of cellular life, and captured numerous genes, including those for translation system components, from eukaryotes, along with some bacterial genes. Even deeper evolutionary reconstructions reveal apparent links between the NCLDV and smaller viruses of eukaryotes, such as adenoviruses, and ultimately, derive all these viruses from tailless bacteriophages.

Genetic exchanges are more frequent in bacteria encoding capsules

Capsules allow bacteria to colonize novel environments, to withstand numerous stresses, and to resist antibiotics. Yet, even though genetic exchanges with other cells should be adaptive under such circumstances, it has been suggested that capsules lower the rates of homologous recombination and horizontal gene transfer. We analysed over one hundred pan-genomes and thousands of bacterial genomes for the evidence of an association between genetic exchanges (or lack thereof) and the presence of a capsule system. We found that bacteria encoding capsules have larger pan-genomes, higher rates of horizontal gene transfer, and higher rates of homologous recombination in their core genomes. Accordingly, genomes encoding capsules have more plasmids, conjugative elements, transposases, prophages, and integrons. Furthermore, capsular loci are frequent in plasmids, and can be found in prophages. These results are valid for Bacteria, independently of their ability to be naturally transformable. Since we have shown previously that capsules are commonly present in nosocomial pathogens, we analysed their co-occurrence with antibiotic resistance genes. Genomes encoding capsules have more antibiotic resistance genes, especially those encoding efflux pumps, and they constitute the majority of the most worrisome nosocomial bacteria. We conclude that bacteria with capsule systems are more genetically diverse and have fast-evolving gene repertoires, which may further contribute to their success in colonizing novel niches such as humans under antibiotic therapy.

Previous works showed that bacteria encoding capsules are better colonizers and are dominant in most environments suggesting a positive role for capsules in the genetic diversification of bacteria. Yet, it has been repeatedly suggested, based almost exclusively studies in few model species, that such bacteria are less diverse and engage in fewer genetic exchanges. Here, we reverse the current paradigm and show that bacteria encoding capsules have larger and more diverse gene repertoires, which change faster by horizontal gene transfer and recombination. Our study alters the traditional view of the capsule as a barrier to gene flow and raises novel questions about the role of capsules in bacterial adaptation.

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

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

Factors driving effective population size and pan-genome evolution in bacteria

Background

Knowledge of population-level processes is essential to understanding the efficacy of selection operating within a species. However, attempts at estimating effective population sizes (Ne) are particularly challenging in bacteria due to their extremely large census populations sizes, varying rates of recombination and arbitrary species boundaries.

Results

In this study, we estimated Ne for 153 species (152 bacteria and one archaeon) defined under a common framework and found that ecological lifestyle and growth rate were major predictors of Ne; and that contrary to theoretical expectations, Ne was unaffected by recombination rate. Additionally, we found that Ne shapes the evolution and diversity of total gene repertoires of prokaryotic species.

Conclusion

Together, these results point to a new model of genome architecture evolution in prokaryotes, in which pan-genome sizes, not individual genome sizes, are governed by drift-barrier evolution.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1272-4) contains supplementary material, which is available to authorized users.

Genome-wide analysis of horizontally acquired genes in the genus Mycobacterium

Horizontal gene transfer (HGT) was attributed as a major driving force for the innovation and evolution of prokaryotic genomes. Previously, multiple research endeavors were undertaken to decipher HGT in different bacterial lineages. The genus Mycobacterium houses some of the most deadly human pathogens; however, the impact of HGT in Mycobacterium has never been addressed in a systematic way. Previous initiatives to explore the genomic imprints of HGTs in Mycobacterium were focused on few selected species, specifically among the members of Mycobacterium tuberculosis complex. Considering the recent availability of a large number of genomes, the current study was initiated to decipher the probable events of HGTs among 109 completely sequenced Mycobacterium species. Our comprehensive phylogenetic analysis with more than 9,000 families of Mycobacterium proteins allowed us to list several instances of gene transfers spread across the Mycobacterium phylogeny. Moreover, by examining the topology of gene phylogenies here, we identified the species most likely to donate and receive these genes and provided a detailed overview of the putative functions these genes may be involved in. Our study suggested that horizontally acquired foreign genes had played an enduring role in the evolution of Mycobacterium genomes and have contributed to their metabolic versatility and pathogenicity.

Evolutionary Genomics of Metchnikovella incurvata (Metchnikovellidae): An Early Branching Microsporidium

Metchnikovellids are highly specialized hyperparasites, which infect and reproduce inside gregarines (Apicomplexa) inhabiting marine invertebrates. Their phylogenetic affiliation was under constant discussion until recently, when analysis of the first near-complete metchnikovellid genome, that of Amphiamblys sp., placed it in a basal position with respect to most Microsporidia. Microsporidia are a highly diversified lineage of extremely reduced parasites related to Rozellida (Rozellosporidia = Rozellomycota = Cryptomycota) within the Holomycota clade of Opisthokonta. By sequencing DNA from a single-isolated infected gregarine cell we obtained an almost complete genome of a second metchnikovellid species, and the first one of a taxonomically described and well-documented species, Metchnikovella incurvata. Our phylogenomic analyses show that, despite being considerably divergent from each other, M. incurvata forms a monophyletic group with Amphiamplys sp., and confirm that metchnikovellids are one of the deep branches of Microsporidia. Comparative genomic analysis demonstrates that, like most Microsporidia, metchnikovellids lack mitochondrial genes involved in energy transduction and are thus incapable of synthesizing their own ATP via mitochondrial oxidative phosphorylation. They also lack the horizontally acquired ATP transporters widespread in most Microsporidia. We hypothesize that a family of mitochondrial carrier proteins evolved to transport ATP from the host into the metchnikovellid cell. We observe the progressive reduction of genes involved in DNA repair pathways along the evolutionary path of Microsporidia, which might explain, at least partly, the extremely high evolutionary rate of the most derived species. Our data also suggest that genome reduction and acquisition of novel genes co-occurred during the adaptation of Microsporidia to their hosts.

Comprehensive Genome-Wide Classification Reveals That Many Plant-Specific Transcription Factors Evolved in Streptophyte Algae

Plant genomes encode many lineage-specific, unique transcription factors. Expansion of such gene families has been previously found to coincide with the evolution of morphological complexity, although comparative analyses have been hampered by severe sampling bias. Here, we make use of the recently increased availability of plant genomes. We have updated and expanded previous rule sets for domain-based classification of transcription associated proteins (TAPs), comprising transcription factors and transcriptional regulators. The genome-wide annotation of these protein families has been analyzed and made available via the novel TAPscan web interface. We find that many TAP families previously thought to be specific for land plants actually evolved in streptophyte (charophyte) algae; 26 out of 36 TAP family gains are inferred to have occurred in the common ancestor of the Streptophyta (uniting the land plants—Embryophyta—with their closest algal relatives). In contrast, expansions of TAP families were found to occur throughout streptophyte evolution. 17 out of 76 expansion events were found to be common to all land plants and thus probably evolved concomitant with the water-to-land-transition.

Whole Genome Sequencing of the Pirarucu (Arapaima gigas) Supports Independent Emergence of Major Teleost Clades

The Pirarucu (Arapaima gigas) is one of the world's largest freshwater fishes and member of the superorder Osteoglossomorpha (bonytongues), one of the oldest lineages of ray-finned fishes. This species is an obligate air-breather found in the basin of the Amazon River with an attractive potential for aquaculture. Its phylogenetic position among bony fishes makes the Pirarucu a relevant subject for evolutionary studies of early teleost diversification. Here, we present, for the first time, a draft genome version of the A. gigas genome, providing useful information for further functional and evolutionary studies. The A. gigas genome was assembled with 103-Gb raw reads sequenced in an Illumina platform. The final draft genome assembly was ∼661 Mb, with a contig N50 equal to 51.23 kb and scaffold N50 of 668 kb. Repeat sequences accounted for 21.69% of the whole genome, and a total of 24,655 protein-coding genes were predicted from the genome assembly, with an average of nine exons per gene. Phylogenomic analysis based on 24 fish species supported the postulation that Osteoglossomorpha and Elopomorpha (eels, tarpons, and bonefishes) are sister groups, both forming a sister lineage with respect to Clupeocephala (remaining teleosts). Divergence time estimations suggested that Osteoglossomorpha and Elopomorpha lineages emerged independently in a period of ∼30 Myr in the Jurassic. The draft genome of A. gigas provides a valuable genetic resource for further investigations of evolutionary studies and may also offer a valuable data for economic applications.

Genomic inference of the metabolism and evolution of the archaeal phylum Aigarchaeota

Microbes of the phylum Aigarchaeota are widely distributed in geothermal environments, but their physiological and ecological roles are poorly understood. Here we analyze six Aigarchaeota metagenomic bins from two circumneutral hot springs in Tengchong, China, to reveal that they are either strict or facultative anaerobes, and most are chemolithotrophs that can perform sulfide oxidation. Applying comparative genomics to the Thaumarchaeota and Aigarchaeota, we find that they both originated from thermal habitats, sharing 1154 genes with their common ancestor. Horizontal gene transfer played a crucial role in shaping genetic diversity of Aigarchaeota and led to functional partitioning and ecological divergence among sympatric microbes, as several key functional innovations were endowed by Bacteria, including dissimilatory sulfite reduction and possibly carbon monoxide oxidation. Our study expands our knowledge of the possible ecological roles of the Aigarchaeota and clarifies their evolutionary relationship to their sister lineage Thaumarchaeota.

Genomes of two archaeal endosymbionts show convergent adaptations to an intracellular lifestyle

Endosymbiosis is a widespread phenomenon in the microbial world and can be based on diverse interactions between endosymbiont and host cell. The vast majority of the known endosymbiotic interactions involve bacteria that have invaded eukaryotic host cells. However, methanogenic archaea have been found to thrive in anaerobic, hydrogenosome-containing protists and it was suggested that this symbiosis is based on the transfer of hydrogen. Here, we used culture-independent genomics approaches to sequence the genomes of two distantly related methanogenic endosymbionts that have been acquired in two independent events by closely related anaerobic ciliate hosts Nyctotherus ovalis and Metopus contortus, respectively. The sequences obtained were then validated as originating from the ciliate endosymbionts by in situ probing experiments. Comparative analyses of these genomes and their closest free-living counterparts reveal that the genomes of both endosymbionts are in an early stage of adaptation towards endosymbiosis as evidenced by the large number of genes undergoing pseudogenization. For instance, the observed loss of genes involved in amino acid biosynthesis in both endosymbiont genomes indicates that the endosymbionts rely on their hosts for obtaining several essential nutrients. Furthermore, the endosymbionts appear to have gained significant amounts of genes of potentially secreted proteins, providing targets for future studies aiming to elucidate possible mechanisms underpinning host-interactions. Altogether, our results provide the first genomic insights into prokaryotic endosymbioses from the archaeal domain of life.

A well-resolved fern nuclear phylogeny reveals the evolution history of numerous transcription factor families

Ferns account for 80% of nonflowering vascular plant species and are the sister lineage of seed plants. Recent molecular phylogenetics have greatly advanced understanding of fern tree of life, but relationships among some major lineages remain unclear. To better resolve the phylogenetic relationships of ferns, we generated transcriptomes from 125 ferns and two lycophytes, with three additional public datasets, to represent all 11 orders and 85% of families of ferns. Our nuclear phylogeny provides strong supports for the monophyly of all four subclasses and nearly all orders and families, and for relationships among these lineages. The only exception is Gleicheniales, which was highly supported as being paraphyletic with Dipteridaceae sister to a clade with Gleicheniaceae + Hymenophyllales. In addition, new and strongly supported phylogenetic relationships are found for suborders and families in Polypodiales. We provide the first dated fern phylogenomic tree using many nuclear genes from a large majority of families, with an estimate for separation of the ancestors of ferns and seed plants in early Devonian at ∼400 Mya and subsequent gradual divergences of fern orders from ∼380 to 200 Mya. Moreover, the newly obtained fern phylogeny provides a framework for gene family analyses, which indicate that the vast majority of transcription factor families found in seed plants were already present in the common ancestor of extant vascular plants. In addition, fern transcription factor genes show similar duplication patterns to those in seed plants, with some showing stable copy number and others displaying independent expansions in both ferns and seed plants. This study provides a robust phylogenetic and gene family evolution framework, as well as rich molecular resources for understanding the morphological and functional evolution in ferns.

Fern genomes elucidate land plant evolution and cyanobacterial symbioses

Ferns are the closest sister group to all seed plants, yet little is known about their genomes other than that they are generally colossal. Here, we report on the genomes of Azolla filiculoides and Salvinia cucullata (Salviniales) and present evidence for episodic whole-genome duplication in ferns-one at the base of 'core leptosporangiates' and one specific to Azolla. One fern-specific gene that we identified, recently shown to confer high insect resistance, seems to have been derived from bacteria through horizontal gene transfer. Azolla coexists in a unique symbiosis with N2-fixing cyanobacteria, and we demonstrate a clear pattern of cospeciation between the two partners. Furthermore, the Azolla genome lacks genes that are common to arbuscular mycorrhizal and root nodule symbioses, and we identify several putative transporter genes specific to Azolla-cyanobacterial symbiosis. These genomic resources will help in exploring the biotechnological potential of Azolla and address fundamental questions in the evolution of plant life.

Parallels between experimental and natural evolution of legume symbionts

The emergence of symbiotic interactions has been studied using population genomics in nature and experimental evolution in the laboratory, but the parallels between these processes remain unknown. Here we compare the emergence of rhizobia after the horizontal transfer of a symbiotic plasmid in natural populations of Cupriavidus taiwanensis, over 10 MY ago, with the experimental evolution of symbiotic Ralstonia solanacearum for a few hundred generations. In spite of major differences in terms of time span, environment, genetic background, and phenotypic achievement, both processes resulted in rapid genetic diversification dominated by purifying selection. We observe no adaptation in the plasmid carrying the genes responsible for the ecological transition. Instead, adaptation was associated with positive selection in a set of genes that led to the co-option of the same quorum-sensing system in both processes. Our results provide evidence for similarities in experimental and natural evolutionary transitions and highlight the potential of comparisons between both processes to understand symbiogenesis.

It is unclear if experimental evolution is a good model for natural processes. Here, Clerissi et al. find parallels between the evolution of symbiosis in rhizobia after horizontal transfer of a plasmid over 10 million years ago and experimentally evolved symbionts.

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.