Genome reduction as the dominant mode of evolution

Adaptive gene loss in the common bean pan-genome during range expansion and domestication

The common bean (Phaseolus vulgaris L.) is a crucial legume crop and an ideal evolutionary model to study adaptive diversity in wild and domesticated populations. Here, we present a common bean pan-genome based on five high-quality genomes and whole-genome reads representing 339 genotypes. It reveals ~234 Mb of additional sequences containing 6,905 protein-coding genes missing from the reference, constituting 49% of all presence/absence variants (PAVs). More non-synonymous mutations are found in PAVs than core genes, probably reflecting the lower effective population size of PAVs and fitness advantages due to the purging effect of gene loss. Our results suggest pan-genome shrinkage occurred during wild range expansion. Selection signatures provide evidence that partial or complete gene loss was a key adaptive genetic change in common bean populations with major implications for plant adaptation. The pan-genome is a valuable resource for food legume research and breeding for climate change mitigation and sustainable agriculture.

Evolutionary plasticity and functional repurposing of the essential metabolic enzyme MoeA

MoeA, or gephyrin in higher eukaryotes, is crucial for molybdenum cofactor biosynthesis required in redox reactions. Gephyrin is a moonlighting protein also involved in postsynaptic receptor clustering, a feature thought to be a recent evolutionary trait. We showed previously that a repurposed copy of MoeA (Glp) is involved in bacterial cell division. To investigate how MoeA acquired multifunctionality, we used phylogenetic inference and protein structure analyses to understand the diversity and evolutionary history of MoeA. Glp-expressing Bacteria have at least two copies of the gene, and our analysis suggests that Glp has lost its enzymatic role. In Archaea we identified an ancestral duplication where one of the paralogs might bind tungsten instead of molybdenum. In eukaryotes, the acquisition of the moonlighting activity of gephyrin comprised three major events: first, MoeA was obtained from Bacteria by early eukaryotes, second, MogA was fused to the N-terminus of MoeA in the ancestor of opisthokonts, and finally, it acquired the function of anchoring GlyR receptors in neurons. Our results support the functional versatility and adaptive nature of the MoeA scaffold, which has been repurposed independently both in eukaryotes and bacteria to carry out analogous functions in network organization at the cell membrane.

Genome-Wide Identification, Characterisation, and Evolution of the Transcription Factor WRKY in Grapevine (Vitis vinifera): New View and Update

WRKYs are a multigenic family of transcription factors that are plant-specific and involved in the regulation of plant development and various stress response processes. However, the evolution of WRKY genes is not fully understood. This family has also been incompletely studied in grapevine, and WRKY genes have been named with different numbers in different studies, leading to great confusion. In this work, 62 Vitis vinifera WRKY genes were identified based on six genomes of different cultivars. All WRKY genes were numbered according to their chromosomal location, and a complete revision of the numbering was performed. Amino acid variability between different cultivars was assessed for the first time and was greater than 5% for some WRKYs. According to the gene structure, all WRKYs could be divided into two groups: more exons/long length and fewer exons/short length. For the first time, some chimeric WRKY genes were found in grapevine, which may play a specific role in the regulation of different processes: VvWRKY17 (an N-terminal signal peptide region followed by a non-cytoplasmic domain) and VvWRKY61 (Frigida-like domain). Five phylogenetic clades A-E were revealed and correlated with the WRKY groups (I, II, III). The evolution of WRKY was studied, and we proposed a WRKY evolution model where there were two dynamic phases of complexity and simplification in the evolution of WRKY.

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. In this manuscript, we present a genome for a laboratory-tractable Entomophthora muscae isolate that infects fruit flies. Our E. muscae assembly is 1.03 Gb, consists of 7810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with recent datasets from entomophthoralean fungi, we show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity. In addition, we find that E. muscae and its closest allies possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Finally, we offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex that support the need for taxonomic revision within this group. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

Bacterial sensor evolved by decreasing complexity

Bacterial receptors feed into multiple signal transduction pathways that regulate a variety of cellular processes including gene expression, second messenger levels and motility. Receptors are typically activated by signal binding to ligand binding domains (LBD). Cache domains are omnipresent LBDs found in bacteria, archaea, and eukaryotes, including humans. They form the predominant family of extracytosolic bacterial LBDs and were identified in all major receptor types. Cache domains are composed of either a single (sCache) or a double (dCache) structural module. The functional relevance of bimodular LBDs remains poorly understood. Here, we identify the PacF chemoreceptor in the phytopathogen Pectobacterium atrosepticum that recognizes formate at the membrane distal module of its dCache domain, triggering chemoattraction. We further demonstrate that a family of formate-specific sCache domains has evolved from a dCache domain, exemplified by PacF, by losing the membrane proximal module. By solving high-resolution structures of two family members in complex with formate, we show that the molecular basis for formate binding at sCache and dCache domains is highly similar, despite their low sequence identity. The apparent loss of the membrane proximal module may be related to the observation that dCache domains bind ligands typically at the membrane distal module, whereas the membrane proximal module is not involved in signal sensing. This work advances our understanding of signal sensing in bacterial receptors and suggests that evolution by reducing complexity may be a common trend shaping their diversity.

Significance

Many bacterial receptors contain multi-modular sensing domains indicative of complex sensory processes. The presence of more than one sensing module likely permits the integration of multiple signals, although, the molecular detail and functional relevance for these complex sensors remain poorly understood. Bimodular sensory domains are likely to have arisen from the fusion or duplication of monomodular domains. Evolution by increasing complexity is generally believed to be a dominant force. Here we reveal the opposite - how a monomodular sensing domain has evolved from a bimodular one. Our findings will thus motivate research to establish whether evolution by decreasing complexity is typical of other sensory domains.

From cyanobacteria and cyanophages to chloroplasts: the fate of the genomes of oxyphototrophs and the genes encoding photosystem II proteins

Chloroplasts are the result of endosymbiosis of cyanobacterial organisms with proto-eukaryotes. The psbA, psbD and psbO genes are present in all oxyphototrophs and encode the D1/D2 proteins of photosystem II (PSII) and PsbO, respectively. PsbO is a peripheral protein that stabilizes the O2-evolving complex in PSII. Of these genes, psbA and psbD remained in the chloroplastic genome, while psbO was transferred to the nucleus. The genomes of selected cyanobacteria, chloroplasts and cyanophages carrying psbA and psbD, respectively, were analysed. The highest density of genes and coding sequences (CDSs) was estimated for the genomes of cyanophages, cyanobacteria and chloroplasts. The synonymous mutation rate (rS) of psbA and psbD in chloroplasts remained almost unchanged and is lower than that of psbO. The results indicate that the decreasing genome size in chloroplasts is more similar to the genome reduction observed in contemporary endosymbiotic organisms than in streamlined genomes of free-living cyanobacteria. The rS of atpA, which encodes the α-subunit of ATP synthase in chloroplasts, suggests that psbA and psbD, and to a lesser extent psbO, are ancient and conservative and arose early in the evolution of oxygenic photosynthesis. The role of cyanophages in the evolution of oxyphototrophs and chloroplastic genomes is discussed.

Retention of an Endosymbiont for the Production of a Single Molecule

Sap-feeding insects often maintain two or more nutritional endosymbionts that act in concert to produce compounds essential for insect survival. Many mealybugs have endosymbionts in a nested configuration: one or two bacterial species reside within the cytoplasm of another bacterium, and together, these bacteria have genomes that encode interdependent sets of genes needed to produce key nutritional molecules. Here, we show that the mealybug Pseudococcus viburni has three endosymbionts, one of which contributes only two unique genes that produce the host nutrition-related molecule chorismate. All three bacterial endosymbionts have tiny genomes, suggesting that they have been coevolving inside their insect host for millions of years.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Genomic Diversity in the Endosymbiotic Bacteria of Human Head Lice

Insects have repeatedly forged symbioses with heritable microbes, gaining novel traits. For the microbe, the transition to symbioses can lead to the degeneration of the symbiont's genome through transmission bottlenecks, isolation, and the loss of DNA repair enzymes. However, some insect-microbial symbioses have persisted for millions of years, suggesting that natural selection slows genetic drift and maintains functional consistency between symbiont populations. By sampling in multiple countries, we examine genomic diversity within a symbiont species, a heritable symbiotic bacterium found only in human head lice. We find that human head louse symbionts contain genetic diversity that appears to have arisen contemporaneously with the appearance of anatomically modern humans within Africa and/or during the colonization of Eurasia by humans. We predict that the observed genetic diversity underlies functional differences in extant symbiont lineages, through the inactivation of genes involved in symbiont membrane construction. Furthermore, we find evidence of additional gene losses prior to the appearance of modern humans, also impacting the symbiont membrane. From this, we conclude that symbiont genome degeneration is proceeding, via gene inactivation and subsequent loss, in human head louse symbionts, while genomic diversity is maintained. Collectively, our results provide a look into the genomic diversity within a single symbiont species and highlight the shared evolutionary history of humans, lice, and bacteria.

Chromosome-scale assembly of the streamlined picoeukaryote Picochlorum sp. SENEW3 genome reveals Rabl-like chromatin structure and potential for C4 photosynthesis

Genome sequencing and assembly of the photosynthetic picoeukaryotic Picochlorum sp. SENEW3 revealed a compact genome with a reduced gene set, few repetitive sequences, and an organized Rabl-like chromatin structure. Hi-C chromosome conformation capture revealed evidence of possible chromosomal translocations, as well as putative centromere locations. Maintenance of a relatively few selenoproteins, as compared to similarly sized marine picoprasinophytes Mamiellales, and broad halotolerance compared to others in Trebouxiophyceae, suggests evolutionary adaptation to variable salinity environments. Such adaptation may have driven size and genome minimization and have been enabled by the retention of a high number of membrane transporters. Identification of required pathway genes for both CAM and C4 photosynthetic carbon fixation, known to exist in the marine mamiellale pico-prasinophytes and seaweed Ulva, but few other chlorophyte species, further highlights the unique adaptations of this robust alga. This high-quality assembly provides a significant advance in the resources available for genomic investigations of this and other photosynthetic picoeukaryotes.

Macroevolutionary dynamics of gene family gain and loss along multicellular eukaryotic lineages

The gain and loss of genes fluctuate over evolutionary time in major eukaryotic clades. However, the full profile of these macroevolutionary trajectories is still missing. To give a more inclusive view on the changes in genome complexity across the tree of life, here we recovered the evolutionary dynamics of gene family gain and loss ranging from the ancestor of cellular organisms to 352 eukaryotic species. We show that in all considered lineages the gene family content follows a common evolutionary pattern, where the number of gene families reaches the highest value at a major evolutionary and ecological transition, and then gradually decreases towards extant organisms. This supports theoretical predictions and suggests that the genome complexity is often decoupled from commonly perceived organismal complexity. We conclude that simplification by gene family loss is a dominant force in Phanerozoic genomes of various lineages, probably underpinned by intense ecological specializations and functional outsourcing.

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. This is in part due to their large genome sizes and difficulty in obtaining sequenceable material. In this manuscript, we leveraged a recently-isolated, laboratory-tractable Entomophthora muscae isolate and improved long-read sequencing to obtain a largely-complete entomophthoralean genome. Our E. muscae assembly is 1.03 Gb, consists of 7,810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with other available (transcriptomic and genomic) datasets from entomophthoralean fungi, we provide new insight into the biology of these understudied pathogens. We offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex. Our findings suggest that substantial taxonomic revision is needed to define species within this group and we provide recommendations for differentiating strains and species in the context of the existing body of E. muscae scientific literature. We show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity, despite the presence of machinery to defend against transposable elements(RNAi). In addition, we find that E. muscae and its closest allies are enriched for M16A peptidases and possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We find that E. muscae has an expanded group of acid-trehalases, consistent with trehalose being the primary sugar component of fly (and insect) hemolymph. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

Complex Endosymbioses I: From Primary to Complex Plastids, Serial Endosymbiotic Events

A considerable part of the diversity of eukaryotic phototrophs consists of algae with plastids that evolved from endosymbioses between two eukaryotes. These complex plastids are characterized by a high number of envelope membranes (more than two) and some of them contain a residual nucleus of the endosymbiotic alga called a nucleomorph. Complex plastid-bearing algae are thus chimeric cell assemblies, eukaryotic symbionts living in a eukaryotic host. In contrast, the primary plastids of the Archaeplastida (plants, green algae, red algae, and glaucophytes) possibly evolved from a single endosymbiosis with a cyanobacterium and are surrounded by two membranes. Complex plastids have been acquired several times by unrelated groups of eukaryotic heterotrophic hosts, suggesting that complex plastids are somewhat easier to obtain than primary plastids. Evidence suggests that complex plastids arose twice independently in the green lineage (euglenophytes and chlorarachniophytes) through secondary endosymbiosis, and four times in the red lineage, first through secondary endosymbiosis in cryptophytes, then by higher-order events in stramenopiles, alveolates, and haptophytes. Engulfment of primary and complex plastid-containing algae by eukaryotic hosts (secondary, tertiary, and higher-order endosymbioses) is also responsible for numerous plastid replacements in dinoflagellates. Plastid endosymbiosis is accompanied by massive gene transfer from the endosymbiont to the host nucleus and cell adaptation of both endosymbiotic partners, which is related to the trophic switch to phototrophy and loss of autonomy of the endosymbiont. Such a process is essential for the metabolic integration and division control of the endosymbiont in the host. Although photosynthesis is the main advantage of acquiring plastids, loss of photosynthesis often occurs in algae with complex plastids. This chapter summarizes the essential knowledge of the acquisition, evolution, and function of complex plastids.

Comparative analysis of cysteine proteases reveals gene family evolution of the group 1 allergens in astigmatic mites

BACKGROUND

Astigmatic mites contain potent allergens that can trigger IgE-mediated immune responses, leading to allergic diseases such as asthma, allergic rhinitis and atopic dermatitis. In house dust mites Dermatophagoides pteronyssinus and Dermatophagoides farinae, group 1 allergens (Der p 1 and Der f 1), characterized as papain-like cysteine proteases, have been defined as the major allergens that have high prevalence and potency. Previous studies of mite group 1 allergens mainly focused on identification, comparison of sequence and structure, as well as the investigation of cross-reactivity. To achieve a comprehensive view of mite group 1 allergens, we performed a comparative genomic analysis of all the cysteine proteases in six astigmatic mite species to elucidate the evolutionary relationships of group 1 allergens.

METHODS

Based on the high-quality and annotated genomes, all the cysteine proteases in six astigmatic mite species were identified by sequence homology search. The phylogenetic relationships, gene synteny and expression levels were revealed by bioinformatic tools. The allergenicity of recombinant cysteine proteases was evaluated by enzyme-linked immunosorbent assay.

RESULTS

Tandem duplication was revealed as the major feature of cysteine protease gene evolution in astigmatic mites. The high IgE-binding capacity and the significant expression level of the cysteine protease DP_007902.01 suggested its potential as a novel group 1 allergen of D. pteronyssinus. In addition, gene decay events were identified in the skin-burrowing parasitic mite Sarcoptes scabiei.

CONCLUSION

This comprehensive analysis provided insights into the evolution of cysteine proteases, as well as the component-resolved diagnosis of mite allergies.

Genome evolution in plants and the origins of innovation

Plant evolution has been characterised by a series of major novelties in their vegetative and reproductive traits that have led to greater complexity. Underpinning this diversification has been the evolution of the genome. When viewed at the scale of the plant kingdom, plant genome evolution has been punctuated by conspicuous instances of gene and whole-genome duplication, horizontal gene transfer and extensive gene loss. The periods of dynamic genome evolution often coincide with the evolution of key traits, demonstrating the coevolution of plant genomes and phenotypes at a macroevolutionary scale. Conventionally, plant complexity and diversity have been considered through the lens of gene duplication and the role of gene loss in plant evolution remains comparatively unexplored. However, in light of reductive evolution across multiple plant lineages, the association between gene loss and plant phenotypic diversity warrants greater attention.

Advances in genome sequencing reveal changes in gene content that contribute to arthropod macroevolution

Current sequencing technology allows for the relatively affordable generation of highly contiguous genomes. Technological advances have made it possible for researchers to investigate the consequences of diverse sorts of genomic variants, such as gene gain and loss. With the extraordinary number of high-quality genomes now available, we take stock of how these genomic variants impact phenotypic evolution. We take care to point out that the identification of genomic variants of interest is only the first step in understanding their impact. Painstaking lab or fieldwork is still required to establish causal relationships between genomic variants and phenotypic evolution. We focus mostly on arthropod research, as this phylum has an impressive degree of phenotypic diversity and is also the subject of much evolutionary genetics research. This article is intended to both highlight recent advances in the field and also to be a primer for learning about evolutionary genetics and genomics.

Lessons from the deep: mechanisms behind diversification of eukaryotic protein complexes

Genetic variation is the major mechanism behind adaptation and evolutionary change. As most proteins operate through interactions with other proteins, changes in protein complex composition and subunit sequence provide potentially new functions. Comparative genomics can reveal expansions, losses and sequence divergence within protein-coding genes, but in silico analysis cannot detect subunit substitutions or replacements of entire protein complexes. Insights into these fundamental evolutionary processes require broad and extensive comparative analyses, from both in silico and experimental evidence. Here, we combine data from both approaches and consider the gamut of possible protein complex compositional changes that arise during evolution, citing examples of complete conservation to partial and total replacement by functional analogues. We focus in part on complexes in trypanosomes as they represent one of the better studied non-animal/non-fungal lineages, but extend insights across the eukaryotes by extensive comparative genomic analysis. We argue that gene loss plays an important role in diversification of protein complexes and hence enhancement of eukaryotic diversity.

Shining the spotlight on the neglected: new high-quality genome assemblies as a gateway to understanding the evolution of Trypanosomatidae

BACKGROUND

Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group.

RESULTS

We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy.

CONCLUSIONS

The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

Genomes of fungi and relatives reveal delayed loss of ancestral gene families and evolution of key fungal traits

Fungi are ecologically important heterotrophs that have radiated into most niches on Earth and fulfil key ecological services. Despite intense interest in their origins, major genomic trends of their evolutionary route from a unicellular opisthokont ancestor to derived multicellular fungi remain poorly known. Here we provide a highly resolved genome-wide catalogue of gene family changes across fungal evolution inferred from the genomes of 123 fungi and relatives. We show that a dominant trend in early fungal evolution has been the gradual shedding of protist genes and the punctuated emergence of innovation by two main gene duplication events. We find that the gene content of non-Dikarya fungi resembles that of unicellular opisthokonts in many respects, owing to the conservation of protist genes in their genomes. The most rapidly duplicating gene groups included extracellular proteins and transcription factors, as well as ones linked to the coordination of nutrient uptake with growth, highlighting the transition to a sessile osmotrophic feeding strategy and subsequent lifestyle evolution as important elements of early fungal history. These results suggest that the genomes of pre-fungal ancestors evolved into the typical filamentous fungal genome by a combination of gradual gene loss, turnover and several large duplication events rather than by abrupt changes. Consequently, the taxonomically defined Fungi represents a genomically non-uniform assemblage of species.

Evolution of a minimal cell

Possessing only essential genes, a minimal cell can reveal mechanisms and processes that are critical for the persistence and stability of life1,2. Here we report on how an engineered minimal cell3,4 contends with the forces of evolution compared with the Mycoplasma mycoides non-minimal cell from which it was synthetically derived. Mutation rates were the highest among all reported bacteria, but were not affected by genome minimization. Genome streamlining was costly, leading to a decrease in fitness of greater than 50%, but this deficit was regained during 2,000 generations of evolution. Despite selection acting on distinct genetic targets, increases in the maximum growth rate of the synthetic cells were comparable. Moreover, when performance was assessed by relative fitness, the minimal cell evolved 39% faster than the non-minimal cell. The only apparent constraint involved the evolution of cell size. The size of the non-minimal cell increased by 80%, whereas the minimal cell remained the same. This pattern reflected epistatic effects of mutations in ftsZ, which encodes a tubulin-homologue protein that regulates cell division and morphology5,6. Our findings demonstrate that natural selection can rapidly increase the fitness of one of the simplest autonomously growing organisms. Understanding how species with small genomes overcome evolutionary challenges provides critical insights into the persistence of host-associated endosymbionts, the stability of streamlined chassis for biotechnology and the targeted refinement of synthetically engineered cells2,7-9.

Multiple Instances of Adaptive Evolution in Aquaporins of Amphibious Fishes

Aquaporins (AQPs) are a highly diverse family of transmembrane proteins involved in osmotic regulation that played an important role in the conquest of land by tetrapods. However, little is known about their possible implication in the acquisition of an amphibious lifestyle in actinopterygian fishes. Herein, we investigated the molecular evolution of AQPs in 22 amphibious actinopterygian fishes by assembling a comprehensive dataset that was used to (1) catalogue AQP paralog members and classes; (2) determine the gene family birth and death process; (3) test for positive selection in a phylogenetic framework; and (4) reconstruct structural protein models. We found evidence of adaptive evolution in 21 AQPs belonging to 5 different classes. Almost half of the tree branches and protein sites that were under positive selection were found in the AQP11 class. The detected sequence changes indicate modifications in molecular function and/or structure, which could be related to adaptation to an amphibious lifestyle. AQP11 orthologues appear to be the most promising candidates to have facilitated the processes of the water-to-land transition in amphibious fishes. Additionally, the signature of positive selection found in the AQP11b stem branch of the Gobiidae clade suggests a possible case of exaptation in this clade.

The Multifaceted Role of Homologous Recombination in a Fastidious Bacterial Plant Pathogen

Homologous recombination plays a key function in the evolution of bacterial genomes. Within Xylella fastidiosa, an emerging plant pathogen with increasing host and geographic ranges, it has been suggested that homologous recombination facilitates host switching, speciation, and the development of virulence. We used 340 whole-genome sequences to study the relationship between inter- and intrasubspecific homologous recombination, random mutation, and natural selection across individual X. fastidiosa genes. Individual gene orthologs were identified and aligned, and a maximum likelihood (ML) gene tree was generated. Each gene alignment and tree pair were then used to calculate gene-wide and branch-specific r/m values (relative effect of recombination to mutation), gene-wide and branch-site nonsynonymous over synonymous substitution rates (dN/dS values; episodic selection), and branch length (as a proxy for mutation rate). The relationships between these variables were evaluated at the global level (i.e., for all genes among and within a subspecies), among specific functional classes (i.e., COGs), and between pangenome components (i.e., accessory versus core genes). Our analysis showed that r/m varied widely among genes as well as across X. fastidiosa subspecies. While r/m and dN/dS values were positively correlated in some instances (e.g., core genes in X. fastidiosa subsp. fastidiosa and both core and accessory genes in X. fastidiosa subsp. multiplex), low correlation coefficients suggested no clear biological significance. Overall, our results indicate that, in addition to its adaptive role in certain genes, homologous recombination acts as a homogenizing and a neutral force across phylogenetic clades, gene functional groups, and pangenome components. IMPORTANCE There is ample evidence that homologous recombination occurs frequently in the economically important plant pathogen Xylella fastidiosa. Homologous recombination has been known to occur among sympatric subspecies and is associated with host-switching events and virulence-linked genes. As a consequence, is it generally assumed that recombinant events in X. fastidiosa are adaptive. This mindset influences expectations of how homologous recombination acts as an evolutionary force as well as how management strategies for X. fastidiosa diseases are determined. Yet, homologous recombination plays roles beyond that of a source for diversification and adaptation. Homologous recombination can act as a DNA repair mechanism, as a means to facilitate nucleotide compositional change, as a homogenization mechanism within populations, or even as a neutral force. Here, we provide a first assessment of long-held beliefs regarding the general role of recombination in adaptation for X. fastidiosa. We evaluate gene-specific variations in homologous recombination rate across three X. fastidiosa subspecies and its relationship to other evolutionary forces (e.g., natural selection, mutation, etc.). These data were used to assess the role of homologous recombination in X. fastidiosa evolution.

In silico analysis of sporulene biosynthesis pathway genes in the members of the class Bacilli

Sporulene, a pentacyclic triterpenoid, was discovered in Bacillus subtilis and is associated with bacterial endospores. However, the study was not further extended, leaving a trail of questions. One such question is what diversity of sporulenes exists among spore-forming members? Considering the sporulene biosynthesis pathway as a fundamental tool to survey the distribution of this molecule, a genome mining study was conducted. Mining for genes encoding putative proteins of sporulene biosynthesis pathway among the class Bacilli members revealed the presence of hepS, hepT, ytpB, and sqhC genes in the members of the family Bacillaceae, Caryophanaceae, Paenibacillaceae, and Sporolactobacillaceae. However, these genes were completely absent in the members of Staphylococcaceae, Lactobacillaceae, Aerococcaceae, Carnobacteriaceae, and Leuconostocaceae. Unlike other probable pathway related proteins, a conserved amino acid domain of putative terpenoid cyclase (YtpB) appeared deep-rooted among the genus Bacillus members. In-depth analysis showed the constant gene arrangement of hepS, hepT, ytpB, and sqhC genes in these members, there by demonstrating the conserved nature of sporulene biosynthesis pathway in the members of the genus Bacillus. Our study suggests confinement of the sporulene biosynthesis pathway to spore-forming members of the class Bacilli, majorly to the genus Bacillus.

Gene loss during a transition to multicellularity

Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited novel interactions, suggesting that gene loss could play a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest gene loss, in addition to gene gain and co-option, may be important for the evolution developmental complexity.

Linking prokaryotic genome size variation to metabolic potential and environment

While theories and models have appeared to explain genome size as a result of evolutionary processes, little work has shown that genome sizes carry ecological signatures. Our work delves into the ecological implications of microbial genome size variation in benthic and pelagic habitats across environmental gradients of the brackish Baltic Sea. While depth is significantly associated with genome size in benthic and pelagic brackish metagenomes, salinity is only correlated to genome size in benthic metagenomes. Overall, we confirm that prokaryotic genome sizes in Baltic sediments (3.47 Mbp) are significantly bigger than in the water column (2.96 Mbp). While benthic genomes have a higher number of functions than pelagic genomes, the smallest genomes coded for a higher number of module steps per Mbp for most of the functions irrespective of their environment. Some examples of this functions are amino acid metabolism and central carbohydrate metabolism. However, we observed that nitrogen metabolism was almost absent in pelagic genomes and was mostly present in benthic genomes. Finally, we also show that Bacteria inhabiting Baltic sediments and water column not only differ in taxonomy, but also in their metabolic potential, such as the Wood-Ljungdahl pathway or the presence of different hydrogenases. Our work shows how microbial genome size is linked to abiotic factors in the environment, metabolic potential and taxonomic identity of Bacteria and Archaea within aquatic ecosystems.

Shuffled ATG8 interacting motifs form an ancestral bridge between UFMylation and autophagy

UFMylation involves the covalent modification of substrate proteins with UFM1 (Ubiquitin-fold modifier 1) and is important for maintaining ER homeostasis. Stalled translation triggers the UFMylation of ER-bound ribosomes and activates C53-mediated autophagy to clear toxic polypeptides. C53 contains noncanonical shuffled ATG8-interacting motifs (sAIMs) that are essential for ATG8 interaction and autophagy initiation. However, the mechanistic basis of sAIM-mediated ATG8 interaction remains unknown. Here, we show that C53 and sAIMs are conserved across eukaryotes but secondarily lost in fungi and various algal lineages. Biochemical assays showed that the unicellular alga Chlamydomonas reinhardtii has a functional UFMylation pathway, refuting the assumption that UFMylation is linked to multicellularity. Comparative structural analyses revealed that both UFM1 and ATG8 bind sAIMs in C53, but in a distinct way. Conversion of sAIMs into canonical AIMs impaired binding of C53 to UFM1, while strengthening ATG8 binding. Increased ATG8 binding led to the autoactivation of the C53 pathway and sensitization of Arabidopsis thaliana to ER stress. Altogether, our findings reveal an ancestral role of sAIMs in UFMylation-dependent fine-tuning of C53-mediated autophagy activation.

Highly Dynamic Gene Family Evolution Suggests Changing Roles for PON Genes Within Metazoa

Change in gene family size has been shown to facilitate adaptation to different selective pressures. This includes gene duplication to increase dosage or diversification of enzymatic substrates and gene deletion due to relaxed selection. We recently found that the PON1 gene, an enzyme with arylesterase and lactonase activity, was lost repeatedly in different aquatic mammalian lineages, suggesting that the PON gene family is responsive to environmental change. We further investigated if these fluctuations in gene family size were restricted to mammals and approximately when this gene family was expanded within mammals. Using 112 metazoan protein models, we explored the evolutionary history of the PON family to characterize the dynamic evolution of this gene family. We found that there have been multiple, independent expansion events in tardigrades, cephalochordates, and echinoderms. In addition, there have been partial gene loss events in monotremes and sea cucumbers and what appears to be complete loss in arthropods, urochordates, platyhelminths, ctenophores, and placozoans. In addition, we show the mammalian expansion to three PON paralogs occurred in the ancestor of all mammals after the divergence of sauropsida but before the divergence of monotremes from therians. We also provide evidence of a novel PON expansion within the brushtail possum. In the face of repeated expansions and deletions in the context of changing environments, we suggest a range of selective pressures, including pathogen infection and mitigation of oxidative damage, are likely influencing the diversification of this dynamic gene family across metazoa.

Bacterial genome reductions: Tools, applications, and challenges

Bacterial cells are widely used to produce value-added products due to their versatility, ease of manipulation, and the abundance of genome engineering tools. However, the efficiency of producing these desired biomolecules is often hindered by the cells’ own metabolism, genetic instability, and the toxicity of the product. To overcome these challenges, genome reductions have been performed, making strains with the potential of serving as chassis for downstream applications. Here we review the current technologies that enable the design and construction of such reduced-genome bacteria as well as the challenges that limit their assembly and applicability. While genomic reductions have shown improvement of many cellular characteristics, a major challenge still exists in constructing these cells efficiently and rapidly. Computational tools have been created in attempts at minimizing the time needed to design these organisms, but gaps still exist in modelling these reductions in silico. Genomic reductions are a promising avenue for improving the production of value-added products, constructing chassis cells, and for uncovering cellular function but are currently limited by their time-consuming construction methods. With improvements to and the creation of novel genome editing tools and in silico models, these approaches could be combined to expedite this process and create more streamlined and efficient cell factories.

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

Background

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

Results

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

Conclusion

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

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02052-3.

The origin of RNA interference: Adaptive or neutral evolution?

The origin of RNA interference (RNAi) is usually explained by a defense-based hypothesis, in which RNAi evolved as a defense against transposable elements (TEs) and RNA viruses and was already present in the last eukaryotic common ancestor (LECA). However, since RNA antisense regulation and double-stranded RNAs (dsRNAs) are ancient and widespread phenomena, the origin of defensive RNAi should have occurred in parallel with its regulative functions to avoid imbalances in gene regulation. Thus, we propose a neutral evolutionary hypothesis for the origin of RNAi in which qualitative system drift from a prokaryotic antisense RNA gene regulation mechanism leads to the formation of RNAi through constructive neutral evolution (CNE). We argue that RNAi was already present in the ancestor of LECA before the need for a new defense system arose and that its presence helped to shape eukaryotic genomic architecture and stability.

Where does RNA interference come from? This Essay describes a new step-by-step evolutionary model of how RNA interference might have originated in early eukaryotes through neutral events from the molecular machinery present in prokaryotes.

Construction and Analysis of the Complete Genome Sequence of Leprosy Agent Mycobacterium lepromatosis

Leprosy is caused by Mycobacterium leprae and Mycobacterium lepromatosis. We report construction and analyses of the complete genome sequence of M. lepromatosis FJ924. The genome contained 3,271,694 nucleotides to encode 1,789 functional genes and 1,564 pseudogenes. It shared 1,420 genes and 885 pseudogenes (71.4%) with M. leprae but differed in 1,281 genes and pseudogenes (28.6%). In phylogeny, the leprosy bacilli started from a most recent common ancestor (MRCA) that diverged ~30 million years ago (Mya) from environmental organism Mycobacterium haemophilum. The MRCA then underwent reductive evolution with pseudogenization, gene loss, and chromosomal rearrangements. Analysis of the shared pseudogenes estimated the pseudogenization event ~14 Mya, shortly before species bifurcation. Afterwards, genomic changes occurred to lesser extent in each species. Like M. leprae, four major types of highly repetitive sequences were detected in M. lepromatosis, contributing to chromosomal rearrangements within and after MRCA. Variations in genes and copy numbers were noted, such as three copies of the gene encoding bifunctional diguanylate cyclase/phosphodiesterase in M. lepromatosis, but single copy in M. leprae; 6 genes encoding the TetR family transcriptional regulators in M. lepromatosis, but 11 such genes in M. leprae; presence of hemW gene in M. lepromatosis, but absence in M. leprae; and others. These variations likely aid unique pathogenesis, such as diffuse lepromatous leprosy associated with M. lepromatosis, while the shared genomic features should explain the common pathogenesis of dermatitis and neuritis in leprosy. Together, these findings and the genomic data of M. lepromatosis may facilitate future research and care for leprosy. IMPORTANCE Leprosy is a dreaded infection that still affects millions of people worldwide. Mycobacterium lepromatosis is a recently recognized cause in addition to the well-known Mycobacterium leprae. M. lepromatosis is likely specific for diffuse lepromatous leprosy, a severe form of the infection and endemic in Mexico. This study constructed and annotated the complete genome sequence of M. lepromatosis FJ924 and performed comparative genomic analyses with related mycobacteria. The results afford new and refined insights into the genome size, gene repertoire, pseudogenes, phylogenomic relationship, genome organization and plasticity, process and timing of reductive evolution, and genetic and proteomic basis for pathogenesis. The availability of the complete M. lepromatosis genome may prove to be useful for future research and care for the infection.

A comparative genomic analysis of lichen-forming fungi reveals new insights into fungal lifestyles

Lichen-forming fungi are mutualistic symbionts of green algae or cyanobacteria. We report the comparative analysis of six genomes of lichen-forming fungi in classes Eurotiomycetes and Lecanoromycetes to identify genomic information related to their symbiotic lifestyle. The lichen-forming fungi exhibited genome reduction via the loss of dispensable genes encoding plant-cell-wall-degrading enzymes, sugar transporters, and transcription factors. The loss of these genes reflects the symbiotic biology of lichens, such as the absence of pectin in the algal cell wall and obtaining specific sugars from photosynthetic partners. The lichens also gained many lineage- and species-specific genes, including those encoding small secreted proteins. These genes are primarily induced during the early stage of lichen symbiosis, indicating their significant roles in the establishment of lichen symbiosis.Our findings provide comprehensive genomic information for six lichen-forming fungi and novel insights into lichen biology and the evolution of symbiosis.

Human follicular mites: Ectoparasites becoming symbionts

Most humans carry mites in the hair follicles of their skin for their entire lives. Follicular mites are the only metazoans tha continuously live on humans. We propose that Demodex folliculorum (Acari) represents a transitional stage from a host-injuring obligate parasite to an obligate symbiont. Here, we describe the profound impact of this transition on the genome and physiology of the mite. Genome sequencing revealed that the permanent host association of D. folliculorum led to an extensive genome reduction through relaxed selection and genetic drift, resulting in the smallest number of protein-coding genes yet identified among panarthropods. Confocal microscopy revealed that this gene loss coincided with an extreme reduction in the number of cells. Single uninucleate muscle cells are sufficient to operate each of the three segments that form each walking leg. While it has been assumed that the reduction of the cell number in parasites starts early in development, we identified a greater total number of cells in the last developmental stage (nymph) than in the terminal adult stage, suggesting that reduction starts at the adult or ultimate stage of development. This is the first evolutionary step in an arthropod species adopting a reductive, parasitic or endosymbiotic lifestyle. Somatic nuclei show underreplication at the diploid stage. Novel eye structures or photoreceptors as well as a unique human host melatonin-guided day/night rhythm are proposed for the first time. The loss of DNA repair genes coupled with extreme endogamy might have set this mite species on an evolutionary dead-end trajectory.

Old genes in new places: A taxon-rich analysis of interdomain lateral gene transfer events

Vertical inheritance is foundational to Darwinian evolution, but fails to explain major innovations such as the rapid spread of antibiotic resistance among bacteria and the origin of photosynthesis in eukaryotes. While lateral gene transfer (LGT) is recognized as an evolutionary force in prokaryotes, the role of LGT in eukaryotic evolution is less clear. With the exception of the transfer of genes from organelles to the nucleus, a process termed endosymbiotic gene transfer (EGT), the extent of interdomain transfer from prokaryotes to eukaryotes is highly debated. A common critique of studies of interdomain LGT is the reliance on the topology of single-gene trees that attempt to estimate more than one billion years of evolution. We take a more conservative approach by identifying cases in which a single clade of eukaryotes is found in an otherwise prokaryotic gene tree (i.e. exclusive presence). Starting with a taxon-rich dataset of over 13,600 gene families and passing data through several rounds of curation, we identify and categorize the function of 306 interdomain LGT events into diverse eukaryotes, including 189 putative EGTs, 52 LGTs into Opisthokonta (i.e. animals, fungi and their microbial relatives), and 42 LGTs nearly exclusive to anaerobic eukaryotes. To assess differential gene loss as an explanation for exclusive presence, we compare branch lengths within each LGT tree to a set of vertically-inherited genes subsampled to mimic gene loss (i.e. with the same taxonomic sampling) and consistently find shorter relative distance between eukaryotes and prokaryotes in LGT trees, a pattern inconsistent with gene loss. Our methods provide a framework for future studies of interdomain LGT and move the field closer to an understanding of how best to model the evolutionary history of eukaryotes.

Genome evolution of a nonparasitic secondary heterotroph, the diatom Nitzschia putrida

Secondary loss of photosynthesis is observed across almost all plastid-bearing branches of the eukaryotic tree of life. However, genome-based insights into the transition from a phototroph into a secondary heterotroph have so far only been revealed for parasitic species. Free-living organisms can yield unique insights into the evolutionary consequence of the loss of photosynthesis, as the parasitic lifestyle requires specific adaptations to host environments. Here, we report on the diploid genome of the free-living diatom Nitzschia putrida (35 Mbp), a nonphotosynthetic osmotroph whose photosynthetic relatives contribute ca. 40% of net oceanic primary production. Comparative analyses with photosynthetic diatoms and heterotrophic algae with parasitic lifestyle revealed that a combination of gene loss, the accumulation of genes involved in organic carbon degradation, a unique secretome, and the rapid divergence of conserved gene families involved in cell wall and extracellular metabolism appear to have facilitated the lifestyle of a free-living secondary heterotroph.

A synthetic 'essentialome' for axenic culturing of 'Candidatus Liberibacter asiaticus'

OBJECTIVE

'Candidatus Liberibacter asiaticus' (CLas) is associated with the devastating citrus 'greening' disease. All attempts to achieve axenic growth and complete Koch's postulates with CLas have failed to date, at best yielding complex cocultures with very low CLas titers detectable only by PCR. Reductive genome evolution has rendered all pathogenic 'Ca. Liberibacter' spp. deficient in multiple key biosynthetic, metabolic and structural pathways that are highly unlikely to be rescued in vitro by media supplementation alone. By contrast, Liberibacter crescens (Lcr) is axenically cultured and its genome is both syntenic and highly similar to CLas. Our objective is to achieve replicative axenic growth of CLas via addition of missing culturability-related Lcr genes.

RESULTS

Bioinformatic analyses identified 405 unique ORFs in Lcr but missing (or truncated) in all 24 sequenced CLas strains. Site-directed mutagenesis confirmed and extended published EZ-Tn5 mutagenesis data, allowing elimination of 310 of these 405 genes as nonessential, leaving 95 experimentally validated Lcr genes as essential for CLas growth in axenic culture. Experimental conditions for conjugation of large GFP-expressing plasmids from Escherichia coli to Lcr were successfully established for the first time, providing a practical method for transfer of large groups of 'essential' Lcr genes to CLas.

An extremely streamlined macronuclear genome in the free-living protozoan Fabrea salina

Ciliated protists are among the oldest unicellular organisms with a heterotrophic lifestyle and share a common ancestor with Plantae. Unlike any other eukaryotes, there are two distinct nuclei in ciliates with separate germline and somatic cell functions. Here, we assembled a near-complete macronuclear genome of Fabrea salina, which belongs to one of the oldest clades of ciliates. Its extremely minimized genome (18.35 Mb) is the smallest among all free-living heterotrophic eukaryotes and exhibits typical streamlined genomic features, including high gene density, tiny introns, and shrinkage of gene paralogs. Gene families involved in hypersaline stress resistance, DNA replication proteins, and mitochondrial biogenesis are expanded, and the accumulation of phosphatidic acid may play an important role in resistance to high osmotic pressure. We further investigated the morphological and transcriptomic changes in the macronucleus during sexual reproduction and highlighted the potential contribution of macronuclear residuals to this process. We believe that the minimized genome generated in this study provides novel insights into the genome streamlining theory and will be an ideal model to study the evolution of eukaryotic heterotrophs.

Genomic evolution and adaptation of arthropod-associated Rickettsia

Rickettsia species are endosymbionts hosted by arthropods and are known to cause mild to fatal diseases in humans. Here, we analyse the evolution and diversity of 34 Rickettsia species using a pangenomic meta-analysis (80 genomes/41 plasmids). Phylogenomic trees showed that Rickettsia spp. diverged into two Spotted Fever groups, a Typhus group, a Canadensis group and a Bellii group, and may have inherited their plasmids from an ancestral plasmid that persisted in some strains or may have been lost by others. The results suggested that the ancestors of Rickettsia spp. might have infected Acari and/or Insecta and probably diverged by persisting inside and/or switching hosts. Pangenomic analysis revealed that the Rickettsia genus evolved through a strong interplay between genome degradation/reduction and/or expansion leading to possible distinct adaptive trajectories. The genus mainly shared evolutionary relationships with α-proteobacteria, and also with γ/β/δ-proteobacteria, cytophagia, actinobacteria, cyanobacteria, chlamydiia and viruses, suggesting lateral exchanges of several critical genes. These evolutionary processes have probably been orchestrated by an abundance of mobile genetic elements, especially in the Spotted Fever and Bellii groups. In this study, we provided a global evolutionary genomic view of the intracellular Rickettsia that may help our understanding of their diversity, adaptation and fitness.

Nuclear genome of a pedinophyte pinpoints genomic innovation and streamlining in the green algae

The genomic diversity underpinning high ecological and species diversity in the green algae (Chlorophyta) remains little known. Here, we aimed to track genome evolution in the Chlorophyta, focusing on loss and gain of homologous genes, and lineage-specific innovations of the core Chlorophyta. We generated a high-quality nuclear genome for pedinophyte YPF701, a sister lineage to others in the core Chlorophyta and incorporated this genome in a comparative analysis with 25 other genomes from diverse Viridiplantae taxa. The nuclear genome of pedinophyte YPF701 has an intermediate size and gene number between those of most prasinophytes and the remainder of the core Chlorophyta. Our results suggest positive selection for genome streamlining in the Pedinophyceae, independent from genome minimisation observed among prasinophyte lineages. Genome expansion was predicted along the branch leading to the UTC clade (classes Ulvophyceae, Trebouxiophyceae and Chlorophyceae) after divergence from their last common ancestor with pedinophytes, with genomic novelty implicated in a range of basic biological functions. Results emphasise multiple independent signals of genome minimisation within the Chlorophyta, as well as the genomic novelty arising before diversification in the UTC clade, which may underpin the success of this species-rich clade in a diversity of habitats.

Enhanced Loss of Retinoic Acid Network Genes in Xenopus laevis Achieves a Tighter Signal Regulation

Retinoic acid (RA) is a major regulatory signal during embryogenesis produced from vitamin A (retinol) by an extensive, autoregulating metabolic and signaling network to prevent fluctuations that result in developmental malformations. Xenopus laevis is an allotetraploid hybrid frog species whose genome includes L (long) and S (short) chromosomes from the originating species. Evolutionarily, the X. laevis subgenomes have been losing either L or S homoeologs in about 43% of genes to generate singletons. In the RA network, out of the 47 genes, about 47% have lost one of the homoeologs, like the genome average. Interestingly, RA metabolism genes from storage (retinyl esters) to retinaldehyde production exhibit enhanced gene loss with 75% singletons out of 28 genes. The effect of this gene loss on RA signaling autoregulation was studied. Employing transient RA manipulations, homoeolog gene pairs were identified in which one homoeolog exhibits enhanced responses or looser regulation than the other, while in other pairs both homoeologs exhibit similar RA responses. CRISPR/Cas9 targeting of individual homoeologs to reduce their activity supports the hypothesis where the RA metabolic network gene loss results in tighter network regulation and more efficient RA robustness responses to overcome complex regulation conditions.

Viroids and Viroid-like Circular RNAs: Do They Descend from Primordial Replicators?

Viroids are a unique class of plant pathogens that consist of small circular RNA molecules, between 220 and 450 nucleotides in size. Viroids encode no proteins and are the smallest known infectious agents. Viroids replicate via the rolling circle mechanism, producing multimeric intermediates which are cleaved to unit length either by ribozymes formed from both polarities of the viroid genomic RNA or by coopted host RNAses. Many viroid-like small circular RNAs are satellites of plant RNA viruses. Ribozyviruses, represented by human hepatitis delta virus, are larger viroid-like circular RNAs that additionally encode the viral nucleocapsid protein. It has been proposed that viroids are direct descendants of primordial RNA replicons that were present in the hypothetical RNA world. We argue, however, that much later origin of viroids, possibly, from recently discovered mobile genetic elements known as retrozymes, is a far more parsimonious evolutionary scenario. Nevertheless, viroids and viroid-like circular RNAs are minimal replicators that are likely to be close to the theoretical lower limit of replicator size and arguably comprise the paradigm for replicator emergence. Thus, although viroid-like replicators are unlikely to be direct descendants of primordial RNA replicators, the study of the diversity and evolution of these ultimate genetic parasites can yield insights into the earliest stages of the evolution of life.

OUP accepted manuscript

Operons are a hallmark of the genomic and regulatory architecture of prokaryotes. However, the mechanism by which two genes placed far apart gradually come close and form operons remains to be elucidated. Here, we propose a new model of the origin of operons: Mobile genetic elements called insertion sequences can facilitate the formation of operons by consecutive insertion–deletion–excision reactions. This mechanism barely leaves traces of insertion sequences and thus difficult to detect in nature. In this study, as a proof-of-concept, we reproducibly demonstrated operon formation in the laboratory. The insertion sequence IS3 and the insertion sequence excision enhancer are genes found in a broad range of bacterial species. We introduced these genes into insertion sequence-less Escherichia coli and found that, supporting our hypothesis, the activity of the two genes altered the expression of genes surrounding IS3, closed a 2.7 kb gap between a pair of genes, and formed new operons. This study shows how insertion sequences can facilitate the rapid formation of operons through locally increasing the structural mutation rates and highlights how coevolution with mobile elements may shape the organization of prokaryotic genomes and gene regulation.

Environmental Determination of Indigenous Bifidobacteria of the Human Intestine

The environmental determination of indigenous (constantly present) bifidobacteria of the human large intestine is considered in this review. Environmental determination (from the Latin determinere, "I determine") is understood as a set of natural phenomena of a habitat (biotope) that determine the role of indigenous microorganisms in the microbiocenosis. Using the symbiotic approach, an attempt is made to identify the environmental conditions for the habitat of bifidobacteria and their physiological effects in the microsymbiocenosis. The features of indigenous bifidobacteria in terms of their nature have been established: evolutionary-genetic (phylogenetic remoteness, genome conservation, metabolic specialization), biochemical (lysozyme resistance, constitutive acetate production), and physiological (microbial "friend-foe" identification, immunoregulation), which are important in adaptation (persistence) and the provision of mutualistic effects and stability of the bifidoflora in the population.

Three topological features of regulatory networks control life-essential and specialized subsystems

Gene regulatory networks (GRNs) play key roles in development, phenotype plasticity, and evolution. Although graph theory has been used to explore GRNs, associations amongst topological features, transcription factors (TFs), and systems essentiality are poorly understood. Here we sought the relationship amongst the main GRN topological features that influence the control of essential and specific subsystems. We found that the K_nn, page rank, and degree are the most relevant GRN features: the ones are conserved along the evolution and are also relevant in pluripotent cells. Interestingly, life-essential subsystems are governed mainly by TFs with intermediary K_nn and high page rank or degree, whereas specialized subsystems are mainly regulated by TFs with low K_nn. Hence, we suggest that the high probability of TFs be toured by a random signal, and the high probability of the signal propagation to target genes ensures the life-essential subsystems' robustness. Gene/genome duplication is the main evolutionary process to rise K_nn as the most relevant feature. Herein, we shed light on unexplored topological GRN features to assess how they are related to subsystems and how the duplications shaped the regulatory systems along the evolution. The classification model generated can be found here: https://github.com/ivanrwolf/NoC/ .

The evolutionary emergence of land plants

There can be no doubt that early land plant evolution transformed the planet but, until recently, how and when this was achieved was unclear. Coincidence in the first appearance of land plant fossils and formative shifts in atmospheric oxygen and CO2 are an artefact of the paucity of earlier terrestrial rocks. Disentangling the timing of land plant bodyplan assembly and its impact on global biogeochemical cycles has been precluded by uncertainty concerning the relationships of bryophytes to one another and to the tracheophytes, as well as the timescale over which these events unfolded. New genome and transcriptome sequencing projects, combined with the application of sophisticated phylogenomic modelling methods, have yielded increasing support for the Setaphyta clade of liverworts and mosses, within monophyletic bryophytes. We consider the evolution of anatomy, genes, genomes and of development within this phylogenetic context, concluding that many vascular plant (tracheophytes) novelties were already present in a comparatively complex last common ancestor of living land plants (embryophytes). Molecular clock analyses indicate that embryophytes emerged in a mid-Cambrian to early Ordovician interval, compatible with hypotheses on their role as geoengineers, precipitating early Palaeozoic glaciations.

Unique phage-bacterium interplay in sponge holobionts from the southern Okinawa Trough hydrothermal vent

Deep-sea hydrothermal vents harbour diverse and abundant animals and their symbiotic microorganisms, which together comprise holobionts. The interplay between bacterial members of holobionts and their viruses (phages) is important for maintaining these symbiotic systems; however, phage-bacterium interactions in deep-sea vent holobionts are not well understood. Marine sponges serve as good models for such studies and are used to unveil phage-bacterium interplay via metagenomic analysis. In three demosponges from deep-sea hydrothermal vent fields in the southern Okinawa Trough, the genomes of a diverse array of symbiotic bacteria, including 10 bacterial phyla, were found to lack intact prophages. Genes related to diverse anti-viral defence systems, for example, the restriction-modification and toxin-antitoxin systems, were abundant in the bacterial communities. We also detected phage genes that could complement or compensate host bacterial metabolism, indicating beneficial roles of phage infection. Our findings provide insight into phage-bacterium interplay in sponges from deep-sea hydrothermal vents.

Evolutionary Dynamics of the Spindle Assembly Checkpoint in Eukaryotes

The tremendous diversity in eukaryotic life forms can ultimately be traced back to evolutionary modifications at the level of molecular networks. Deep understanding of these modifications will not only explain cellular diversity, but will also uncover different ways to execute similar processes and expose the evolutionary 'rules' that shape the molecular networks. Here, we review the evolutionary dynamics of the spindle assembly checkpoint (SAC), a signaling network that guards fidelity of chromosome segregation. We illustrate how the interpretation of divergent SAC systems in eukaryotic species is facilitated by combining detailed molecular knowledge of the SAC and extensive comparative genome analyses. Ultimately, expanding this to other core cellular systems and experimentally interrogating such systems in organisms from all major lineages may start outlining the routes to and eventual manifestation of the cellular diversity of eukaryotic life.

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

Motivation

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

Results

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

Availability and implementation

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

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Metabolic adaptations to marine environments: molecular diversity and evolution of ovothiol biosynthesis in Bacteria

Ovothiols are sulfur-containing amino acids synthesized by marine invertebrates, protozoans, and bacteria. They act as pleiotropic molecules in signaling and protection against oxidative stress. The discovery of ovothiol biosynthetic enzymes, sulfoxide synthase OvoA and β-lyase OvoB, paves the way for a systematic investigation of ovothiol distribution and molecular diversification in nature. In this work, we conducted genomic and metagenomics data mining to investigate the distribution and diversification of ovothiol biosynthetic enzymes in Bacteria. We identified the bacteria endowed with this secondary metabolic pathway, described their taxonomy, habitat and biotic interactions in order to provide insight into their adaptation to specific environments. We report that OvoA and OvoB are mostly encountered in marine aerobic Proteobacteria, some of them establishing symbiotic or parasitic relationships with other organisms. We identified a horizontal gene transfer event of OvoB from Bacteroidetes living in symbiosis with Hydrozoa. Our search within the Ocean Gene Atlas revealed the occurrence of ovothiol biosynthetic genes in Proteobacteria living in a wide range of pelagic and highly oxygenated environments. Finally, we tracked the evolutionary history of ovothiol biosynthesis from marine bacteria to unicellular eukaryotes and metazoans. Our analysis provides new conceptual elements to unravel the evolutionary and ecological significance of ovothiol biosynthesis.