The repatterning of eukaryotic genomes by random genetic drift

Genomic Gigantism is not Associated with Reduced Selection Efficiency in Neotropical Salamanders

Genome size variation in eukaryotes has myriad effects on organismal biology from the genomic to whole-organism level. Large genome size may be associated with lower selection efficiency because lower effective population sizes allow fixation of deleterious mutations via genetic drift, increasing genome size and decreasing selection efficiency. Because of a hypothesized negative relationship between genome size and recombination rate per base pair, increased genome size could also increase the effect of linked selection in the genome, decreasing the efficiency with which natural selection can fix or remove mutations. We used a transcriptomic dataset of 15 and a subset of six Neotropical salamander species ranging in genome size from 12 to 87 pg to study the relationship between genome size and efficiency of selection. We estimated dN/dS of salamanders with small and large genomes and tested for relaxation of selection in the larger genomes. Contrary to our expectations, we did not find a significant relationship between genome size and selection efficiency or strong evidence for higher dN/dS values in species with larger genomes for either species set. We also found little evidence for relaxation of selection in species with larger genomes. A positive correlation between genome size and range size (a proxy of population size) in this group disagrees with predictions of stronger drift in species with larger genomes. Our results highlight the complex interactions between the many forces shaping genomic variation in organisms with genomic gigantism.

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.

Parallel gene size and isoform expansion of ancient neuronal genes

How nervous systems evolved is a central question in biology. A diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning hundreds of thousands of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. Although many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and, in many cases, the emergence of neurons as dedicated cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes potentially driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection, as demonstrated by low dN/dS ratios, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

Decoding populations in the ocean microbiome

Understanding the characteristics and structure of populations is fundamental to comprehending ecosystem processes and evolutionary adaptations. While the study of animal and plant populations has spanned a few centuries, microbial populations have been under scientific scrutiny for a considerably shorter period. In the ocean, analyzing the genetic composition of microbial populations and their adaptations to multiple niches can yield important insights into ecosystem function and the microbiome's response to global change. However, microbial populations have remained elusive to the scientific community due to the challenges associated with isolating microorganisms in the laboratory. Today, advancements in large-scale metagenomics and metatranscriptomics facilitate the investigation of populations from many uncultured microbial species directly from their habitats. The knowledge acquired thus far reveals substantial genetic diversity among various microbial species, showcasing distinct patterns of population differentiation and adaptations, and highlighting the significant role of selection in structuring populations. In the coming years, population genomics is expected to significantly increase our understanding of the architecture and functioning of the ocean microbiome, providing insights into its vulnerability or resilience in the face of ongoing global change. Video Abstract.

Genome evolution and divergence in cis-regulatory architecture is associated with condition-responsive development in horned dung beetles

Phenotypic plasticity is thought to be an important driver of diversification and adaptation to environmental variation, yet the genomic mechanisms mediating plastic trait development and evolution remain poorly understood. The Scarabaeinae, or true dung beetles, are a species-rich clade of insects recognized for their highly diversified nutrition-responsive development including that of cephalic horns-evolutionarily novel, secondary sexual weapons that exhibit remarkable intra- and interspecific variation. Here, we investigate the evolutionary basis for horns as well as other key dung beetle traits via comparative genomic and developmental assays. We begin by presenting chromosome-level genome assemblies of three dung beetle species in the tribe Onthophagini (> 2500 extant species) including Onthophagus taurus, O. sagittarius, and Digitonthophagus gazella. Comparing these assemblies to those of seven other species across the order Coleoptera identifies evolutionary changes in coding sequence associated with metabolic regulation of plasticity and metamorphosis. We then contrast chromatin accessibility in developing head horn tissues of high- and low-nutrition O. taurus males and females and identify distinct cis-regulatory architectures underlying nutrition- compared to sex-responsive development, including a large proportion of recently evolved regulatory elements sensitive to horn morph determination. Binding motifs of known and new candidate transcription factors are enriched in these nutrition-responsive open chromatin regions. Our work highlights the importance of chromatin state regulation in mediating the development and evolution of plastic traits, demonstrates gene networks are highly evolvable transducers of environmental and genetic signals, and provides new reference-quality genomes for three species that will bolster future developmental, ecological, and evolutionary studies of this insect group.

The Concept of Neuromuscular Repatterning in Dancers: A Systematic Review

Repatterning is a term that can be used in different fields, including genetics, molecular biology, neurology, psychology, or rehabilitation. Our aim is to identify the key concept of neuromuscular repatterning in somatic training programmes for dancers. A systematic search of eight databases was conducted using the PRISMA (Preferred Reporting Items for Systematic Review and Meta-Analysis) guidelines. The Quality Assessment Tool for Quantitative Studies and the Oxford Levels of Evidence scales were used. The search yielded 1218 results, of which 5 met the inclusion criteria. Five studies (n = 5) were related to psychosomatic health (n = 5), two studies highlighted integration and inter-articular connectivity in movement (n = 2), four studies investigated the neurological component of alignment and efficiency in dance practice (n = 4), and two studies investigated self-confidence (n = 2). Five studies (n = 5) used imagery based on the anatomical and physiological experience of body systems as the main analytical method. Four studies (n = 4) used developmental movement through Bartenieff fundamentals as the main technique for this methodology. Developmental movement and imagery are two methodologies strongly connected to the concept of neuromuscular repatterning in somatic training programmes for dancers. The acquisition of further quantitative experimental or quasi-experimental studies is warranted to better define the level of improvement or impact of neuromuscular repatterning in dancers.

Ancient origins of complex neuronal genes

How nervous systems evolved is a central question in biology. An increasing diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning millions of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. While many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and in many cases the emergence of dedicated neuronal cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection as demonstrated by low dN/dS scores, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

Nascent evolution of recombination rate differences as a consequence of chromosomal rearrangements

Reshuffling of genetic variation occurs both by independent assortment of chromosomes and by homologous recombination. Such reshuffling can generate novel allele combinations and break linkage between advantageous and deleterious variants which increases both the potential and the efficacy of natural selection. Here we used high-density linkage maps to characterize global and regional recombination rate variation in two populations of the wood white butterfly (Leptidea sinapis) that differ considerably in their karyotype as a consequence of at least 27 chromosome fissions and fusions. The recombination data were compared to estimates of genetic diversity and measures of selection to assess the relationship between chromosomal rearrangements, crossing over, maintenance of genetic diversity and adaptation. Our data show that the recombination rate is influenced by both chromosome size and number, but that the difference in the number of crossovers between karyotypes is reduced as a consequence of a higher frequency of double crossovers in larger chromosomes. As expected from effects of selection on linked sites, we observed an overall positive association between recombination rate and genetic diversity in both populations. Our results also revealed a significant effect of chromosomal rearrangements on the rate of intergenic diversity change between populations, but limited effects on polymorphisms in coding sequence. We conclude that chromosomal rearrangements can have considerable effects on the recombination landscape and consequently influence both maintenance of genetic diversity and efficiency of selection in natural populations.

Intracellular common gardens reveal niche differentiation in transposable element community during bacterial adaptive evolution

The distribution and abundance of transposable elements across the tree of life have significantly shaped the evolution of cellular organisms, but the underlying mechanisms shaping these ecological patterns remain elusive. Here we establish a "common garden" approach to study causal ecological interactions between a xenogeneic conditional lethal sacB gene and the community of transposable insertion sequences (ISs) in a multipartite prokaryote genome. Xenogeneic sacB of low, medium, or high GC content was individually inserted into three replicons of a model bacterium Sinorhizobium fredii, and exhibited replicon- and GC-dependent variation in genetic stability. This variation was largely attributable to multidimensional niche differentiation for IS community members. The transposition efficiency of major active ISs depended on the nucleoid-associated xenogeneic silencer MucR. Experimentally eliminating insertion activity of specific ISs by deleting MucR strongly demonstrated a dominant role of niche differentiation among ISs. This intracellular common garden approach in the experimental evolution context allows not only for evaluating genetic stability of natural and synthetic xenogeneic genes of different sequence signatures in host cells but also for tracking and testing causal relationships in unifying ecological principles in genome ecology.

Evolutionary scaling of maximum growth rate with organism size

Data from nearly 1000 species reveal the upper bound to rates of biomass production achievable by natural selection across the Tree of Life. For heterotrophs, maximum growth rates scale positively with organism size in bacteria but negatively in eukaryotes, whereas for phototrophs, the scaling is negligible for cyanobacteria and weakly negative for eukaryotes. These results have significant implications for understanding the bioenergetic consequences of the transition from prokaryotes to eukaryotes, and of the expansion of some groups of the latter into multicellularity. The magnitudes of the scaling coefficients for eukaryotes are significantly lower than expected under any proposed physical-constraint model. Supported by genomic, bioenergetic, and population-genetic data and theory, an alternative hypothesis for the observed negative scaling in eukaryotes postulates that growth-diminishing mutations with small effects passively accumulate with increasing organism size as a consequence of associated increases in the power of random genetic drift. In contrast, conditional on the structural and functional features of ribosomes, natural selection has been able to promote bacteria with the fastest possible growth rates, implying minimal conflicts with both bioenergetic constraints and random genetic drift. If this extension of the drift-barrier hypothesis is correct, the interpretations of comparative studies of biological traits that have traditionally ignored differences in population-genetic environments will require revisiting.

The Linkage-Disequilibrium and Recombinational Landscape in Daphnia pulex

By revealing the influence of recombinational activity beyond what can be achieved with controlled crosses, measures of linkage disequilibrium (LD) in natural populations provide a powerful means of defining the recombinational landscape within which genes evolve. In one of the most comprehensive studies of this sort ever performed, involving whole-genome analyses on nearly 1,000 individuals of the cyclically parthenogenetic microcrustacean Daphnia pulex, the data suggest a relatively uniform pattern of recombination across the genome. Patterns of LD are quite consistent among populations; average rates of recombination are quite similar for all chromosomes; and although some chromosomal regions have elevated recombination rates, the degree of inflation is not large, and the overall spatial pattern of recombination is close to the random expectation. Contrary to expectations for models in which crossing-over is the primary mechanism of recombination, and consistent with data for other species, the distance-dependent pattern of LD indicates excessively high levels at both short and long distances and unexpectedly low levels of decay at long distances, suggesting significant roles for factors such as nonindependent mutation, population subdivision, and recombination mechanisms unassociated with crossing over. These observations raise issues regarding the classical LD equilibrium model widely applied in population genetics to infer recombination rates across various length scales on chromosomes.

Evaluation of genome size and phylogenetic relationships of the Saccharum complex species

"Saccharum complex" is a hypothetical group of species, which is supposed to be involved in the origin of modern sugarcane, and displays large genomes and complex chromosomal alterations. The utilization of restricted parents in breeding programs of modern cultivated sugarcane has resulted in a genetic blockage, which controlled its improvement because of the limited genetic diversity. The use of wild relatives is an effective way to broaden the genetic composition of cultivated sugarcane. Due to the infrequent characterization of genomes, the potential of wild relatives is diffused in improving the cultivated sugarcane. To characterize the genomes of the wild relatives, the genome size and phylogenetic relationships among eight species, including Saccharum spontaneum, Erianthus arundinaceus, E. fulvus, E. rockii, Narenga porphyrocoma, Miscanthus floridulus, Eulalia quadrinervis, and M. sinensis were evaluated based on flow cytometry, genome surveys, K-mer analysis, chloroplast genome sequencing, and whole-genome SNPs analysis. We observed highly heterozygous genomes of S. spontaneum, E. rockii, and E. arundinaceus and the highly repetitive genome of E. fulvus. The genomes of Eulalia quadrinervis, N. porphyrocoma, M. sinensis, and M. floridulus were highly complex. Phylogenetic results of the two approaches were dissimilar, however, both indicate E. fulvus displayed closer relationships to Miscanthus and Saccharum than other species of Saccharum complex. Eulalia quadrinervis was more closely related to M. floridulus than M. sinensis; E. arundinaceus differ significantly from Miscanthus, Narenga, and Saccharum, but was relatively close to Erianthus. We proved the point of E. rockii and E. fulvus should not be classified as one genus, and E. fulvus should be classified as the Saccharum genus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03338-5.

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.

Characterising genome architectures using genome decomposition analysis

Genome architecture describes how genes and other features are arranged in genomes. These arrangements reflect the evolutionary pressures on genomes and underlie biological processes such as chromosomal segregation and the regulation of gene expression. We present a new tool called Genome Decomposition Analysis (GDA) that characterises genome architectures and acts as an accessible approach for discovering hidden features of a genome assembly. With the imminent deluge of high-quality genome assemblies from projects such as the Darwin Tree of Life and the Earth BioGenome Project, GDA has been designed to facilitate their exploration and the discovery of novel genome biology. We highlight the effectiveness of our approach in characterising the genome architectures of single-celled eukaryotic parasites from the phylum Apicomplexa and show that it scales well to large genomes.

Ciliate Research. From Myth to Trendsetting Science

This special issue of the Journal of Eukaryotic Microbiology (JEM) summarizes achievements obtained by generations of researchers with ciliates in widely different disciplines. In fact, ciliates range among the first cells seen under the microscope centuries ago. Their beauty made them an object of scientia amabilis and their manifold reactions made them attractive for college experiments and finally challenged causal analyses at the cellular level. Some of this work was honored by a Nobel Prize. Some observations yielded a baseline for additional novel discoveries, occasionally facilitated by specific properties of some ciliates. This also offers some advantage in the exploration of closely related parasites (malaria). Articles contributed here by colleagues from all over the world encompass a broad spectrum of ciliate life, from genetics to evolution, from molecular cell biology to ecology, from intercellular signaling to epigenetics etc. This introductory chapter, largely based on my personal perception, aims at integrating work presented in this special issue of JEM into a broader historical context up to current research.

Non-Darwinian Molecular Biology

With the discovery of the double helical structure of DNA, a shift occurred in how biologists investigated questions surrounding cellular processes, such as protein synthesis. Instead of viewing biological activity through the lens of chemical reactions, this new field used biological information to gain a new profound view of how biological systems work. Molecular biologists asked new types of questions that would have been inconceivable to the older generation of researchers, such as how cellular machineries convert inherited biological information into functional molecules like proteins. This new focus on biological information also gave molecular biologists a way to link their findings to concepts developed by genetics and the modern synthesis. However, by the late 1960s this all changed. Elevated rates of mutation, unsustainable genetic loads, and high levels of variation in populations, challenged Darwinian evolution, a central tenant of the modern synthesis, where adaptation was the main driver of evolutionary change. Building on these findings, Motoo Kimura advanced the neutral theory of molecular evolution, which advocates that selection in multicellular eukaryotes is weak and that most genomic changes are neutral and due to random drift. This was further elaborated by Jack King and Thomas Jukes, in their paper “Non-Darwinian Evolution”, where they pointed out that the observed changes seen in proteins and the types of polymorphisms observed in populations only become understandable when we take into account biochemistry and Kimura’s new theory. Fifty years later, most molecular biologists remain unaware of these fundamental advances. Their adaptionist viewpoint fails to explain data collected from new powerful technologies which can detect exceedingly rare biochemical events. For example, high throughput sequencing routinely detects RNA transcripts being produced from almost the entire genome yet are present less than one copy per thousand cells and appear to lack any function. Molecular biologists must now reincorporate ideas from classical biochemistry and absorb modern concepts from molecular evolution, to craft a new lens through which they can evaluate the functionality of transcriptional units, and make sense of our messy, intricate, and complicated genome.

Satellitome comparison of two oedipodine grasshoppers highlights the contingent nature of satellite DNA evolution

BACKGROUND

The full catalog of satellite DNA (satDNA) within a same genome constitutes the satellitome. The Library Hypothesis predicts that satDNA in relative species reflects that in their common ancestor, but the evolutionary mechanisms and pathways of satDNA evolution have never been analyzed for full satellitomes. We compare here the satellitomes of two Oedipodine grasshoppers (Locusta migratoria and Oedaleus decorus) which shared their most recent common ancestor about 22.8 Ma ago.

RESULTS

We found that about one third of their satDNA families (near 60 in every species) showed sequence homology and were grouped into 12 orthologous superfamilies. The turnover rate of consensus sequences was extremely variable among the 20 orthologous family pairs analyzed in both species. The satDNAs shared by both species showed poor association with sequence signatures and motives frequently argued as functional, except for short inverted repeats allowing short dyad symmetries and non-B DNA conformations. Orthologous satDNAs frequently showed different FISH patterns at both intra- and interspecific levels. We defined indices of homogenization and degeneration and quantified the level of incomplete library sorting between species.

CONCLUSIONS

Our analyses revealed that satDNA degenerates through point mutation and homogenizes through partial turnovers caused by massive tandem duplications (the so-called satDNA amplification). Remarkably, satDNA amplification increases homogenization, at intragenomic level, and diversification between species, thus constituting the basis for concerted evolution. We suggest a model of satDNA evolution by means of recursive cycles of amplification and degeneration, leading to mostly contingent evolutionary pathways where concerted evolution emerges promptly after lineages split.

Environmental stress leads to genome streamlining in a widely distributed species of soil bacteria

Bacteria have highly flexible pangenomes, which are thought to facilitate evolutionary responses to environmental change, but the impacts of environmental stress on pangenome evolution remain unclear. Using a landscape pangenomics approach, I demonstrate that environmental stress leads to consistent, continuous reduction in genome content along four environmental stress gradients (acidity, aridity, heat, salinity) in naturally occurring populations of Bradyrhizobium diazoefficiens (widespread soil-dwelling plant mutualists). Using gene-level network and duplication functional traits to predict accessory gene distributions across environments, genes predicted to be superfluous are more likely lost in high stress, while genes with multi-functional roles are more likely retained. Genes with higher probabilities of being lost with stress contain significantly higher proportions of codons under strong purifying and positive selection. Gene loss is widespread across the entire genome, with high gene-retention hotspots in close spatial proximity to core genes, suggesting Bradyrhizobium has evolved to cluster essential-function genes (accessory genes with multifunctional roles and core genes) in discrete genomic regions, which may stabilise viability during genomic decay. In conclusion, pangenome evolution through genome streamlining are important evolutionary responses to environmental change. This raises questions about impacts of genome streamlining on the adaptive capacity of bacterial populations facing rapid environmental change.

Transposable Elements Contribute to Genome Dynamics and Gene Expression Variation in the Fungal Plant Pathogen Verticillium dahliae

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization and adaptation. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as polymorphic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the total TE complement. As expected, the polyrmorphic TEs are enriched in the adaptive genomic regions. Besides, we observed an association of polymorphic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of adaptive genomic regions.

Reconstructing the History of Variation in Effective Population Size along Phylogenies

The nearly neutral theory predicts specific relations between effective population size (Ne) and patterns of divergence and polymorphism, which depend on the shape of the distribution of fitness effects (DFE) of new mutations. However, testing these relations is not straightforward, owing to the difficulty in estimating Ne. Here, we introduce an integrative framework allowing for an explicit reconstruction of the phylogenetic history of Ne, thus leading to a quantitative test of the nearly neutral theory and an estimation of the allometric scaling of the ratios of nonsynonymous over synonymous polymorphism (πN/πS) and divergence (dN/dS) with respect to Ne. As an illustration, we applied our method to primates, for which the nearly neutral predictions were mostly verified. Under a purely nearly neutral model with a constant DFE across species, we find that the variation in πN/πS and dN/dS as a function of Ne is too large to be compatible with current estimates of the DFE based on site frequency spectra. The reconstructed history of Ne shows a 10-fold variation across primates. The mutation rate per generation u, also reconstructed over the tree by the method, varies over a 3-fold range and is negatively correlated with Ne. As a result of these opposing trends for Ne and u, variation in πS is intermediate, primarily driven by Ne but substantially influenced by u. Altogether, our integrative framework provides a quantitative assessment of the role of Ne and u in modulating patterns of genetic variation, while giving a synthetic picture of their history over the clade.

Random Penetrance of Mutations Among Individuals: A New Type of Genetic Drift in Molecular Evolution

The determinative view of mutation penetrance is a fundamental assumption for the building of molecular evolutionary theory: individuals in the population with the same genotype have the same fitness effect. Since this view has been constantly challenged by experimental evidence, it is desirable to examine to what extent violation of this view could affect our understanding of molecular evolution. To this end, the author formulated a new theory of molecular evolution under a random model of penetrance: for any individual with the same mutational genotype, the coefficient of selection is a random variable. It follows that, in addition to the conventional N _e-genetic drift (N _e is the effective population size), the variance of penetrance among individuals (ε ²) represents a new type of genetic drift, coined by the ε ²-genetic drift. It has been demonstrated that these two genetic drifts together provided new insights on the nearly neutral evolution: the evolutionary rate is inversely related to the log-of-N _e when the ε ²-genetic drift is nontrivial. This log-of-N _e feature of ε ²-genetic drift did explain well why the d _N /d _S ratio (the nonsynonymous rate to the synonymous rate) in humans is only as twofold as that in mice, while the effective population size (N _e) of mice is about two-magnitude larger than that of humans. It was estimated that, for the first time, the variance of random penetrance in mammalian genes was approximately ε ²≈ 5.89 × 10^-3.

Exploring Potential Signals of Selection for Disordered Residues in Prokaryotic and Eukaryotic Proteins

Intrinsically disordered proteins (IDPs) are an important class of proteins in all domains of life for their functional importance. However, how nature has shaped the disorder potential of prokaryotic and eukaryotic proteins is still not clearly known. Randomly generated sequences are free of any selective constraints, thus these sequences are commonly used as null models. Considering different types of random protein models, here we seek to understand how the disorder potential of natural eukaryotic and prokaryotic proteins differs from random sequences. Comparing proteome-wide disorder content between real and random sequences of 12 model organisms, we noticed that eukaryotic proteins are enriched in disordered regions compared to random sequences, but in prokaryotes such regions are depleted. By analyzing the position-wise disorder profile, we show that there is a generally higher disorder near the N- and C-terminal regions of eukaryotic proteins as compared to the random models; however, either no or a weak such trend was found in prokaryotic proteins. Moreover, here we show that this preference is not caused by the amino acid or nucleotide composition at the respective sites. Instead, these regions were found to be endowed with a higher fraction of protein–protein binding sites, suggesting their functional importance. We discuss several possible explanations for this pattern, such as improving the efficiency of protein–protein interaction, ribosome movement during translation, and post-translational modification. However, further studies are needed to clearly understand the biophysical mechanisms causing the trend.

GC-content biases in protein-coding genes act as an "mRNA identity" feature for nuclear export

It has long been observed that human protein-coding genes have a particular distribution of GC-content: the 5' end of these genes has high GC-content while the 3' end has low GC-content. In 2012, it was proposed that this pattern of GC-content could act as an mRNA identity feature that would lead to it being better recognized by the cellular machinery to promote its nuclear export. In contrast, junk RNA, which largely lacks this feature, would be retained in the nucleus and targeted for decay. Now two recent papers have provided evidence that GC-content does promote the nuclear export of many mRNAs in human cells.

What We Don't Know About Diet-Breadth Evolution in Herbivorous Insects

Half a million species of herbivorous insects have been described. Most of them are diet specialists, using only a few plant species as hosts. Biologists suspect that their specificity is key to their diversity. But why do herbivorous insects tend to be diet specialists? In this review, we catalog a broad range of explanations. We review the evidence for each and suggest lines of research to obtain the evidence we lack. We then draw attention to a second major question, namely how changes in diet breadth affect the rest of a species’ biology. In particular, we know little about how changes in diet breadth feed back on genetic architecture, the population genetic environment, and other aspects of a species’ ecology. Knowing more about how generalists and specialists differ should go a long way toward sorting out potential explanations of specificity, and yield a deeper understanding of herbivorous insect diversity.

Genome streamlining in a minute herbivore that manipulates its host plant

The tomato russet mite, Aculops lycopersici, is among the smallest animals on earth. It is a worldwide pest on tomato and can potently suppress the host's natural resistance. We sequenced its genome, the first of an eriophyoid, and explored whether there are genomic features associated with the mite's minute size and lifestyle. At only 32.5 Mb, the genome is the smallest yet reported for any arthropod and, reminiscent of microbial eukaryotes, exceptionally streamlined. It has few transposable elements, tiny intergenic regions, and is remarkably intron-poor, as more than 80% of coding genes are intronless. Furthermore, in accordance with ecological specialization theory, this defense-suppressing herbivore has extremely reduced environmental response gene families such as those involved in chemoreception and detoxification. Other losses associate with this species' highly derived body plan. Our findings accelerate the understanding of evolutionary forces underpinning metazoan life at the limits of small physical and genome size.

Functional Long Non-coding RNAs Evolve from Junk Transcripts

Transcriptome studies reveal pervasive transcription of complex genomes, such as those of mammals. Despite popular arguments for functionality of most, if not all, of these transcripts, genome-wide analysis of selective constraints indicates that most of the produced RNA are junk. However, junk is not garbage. On the contrary, junk transcripts provide the raw material for the evolution of diverse long non-coding (lnc) RNAs by non-adaptive mechanisms, such as constructive neutral evolution. The generation of many novel functional entities, such as lncRNAs, that fuels organismal complexity does not seem to be driven by strong positive selection. Rather, the weak selection regime that dominates the evolution of most multicellular eukaryotes provides ample material for functional innovation with relatively little adaptation involved.

How Much Does Ne Vary Among Species?

Genetic drift is an important evolutionary force of strength inversely proportional to Ne , the effective population size. The impact of drift on genome diversity and evolution is known to vary among species, but quantifying this effect is a difficult task. Here we assess the magnitude of variation in drift power among species of animals via its effect on the mutation load - which implies also inferring the distribution of fitness effects of deleterious mutations. To this aim, we analyze the nonsynonymous (amino-acid changing) and synonymous (amino-acid conservative) allele frequency spectra in a large sample of metazoan species, with a focus on the primates vs. fruit flies contrast. We show that a Gamma model of the distribution of fitness effects is not suitable due to strong differences in estimated shape parameters among taxa, while adding a class of lethal mutations essentially solves the problem. Using the Gamma + lethal model and assuming that the mean deleterious effects of nonsynonymous mutations is shared among species, we estimate that the power of drift varies by a factor of at least 500 between large-Ne and small-Ne species of animals, i.e., an order of magnitude more than the among-species variation in genetic diversity. Our results are relevant to Lewontin's paradox while further questioning the meaning of the Ne parameter in population genomics.

Intron and gene size expansion during nervous system evolution

Background

The evolutionary radiation of animals was accompanied by extensive expansion of gene and genome sizes, increased isoform diversity, and complexity of regulation.

Results

Here we show that the longest genes are enriched for expression in neuronal tissues of diverse vertebrates and of invertebrates. Additionally, we show that neuronal gene size expansion occurred predominantly through net gains in intron size, with a positional bias toward the 5′ end of each gene.

Conclusions

We find that intron and gene size expansion is a feature of many genes whose expression is enriched in nervous systems. We speculate that unique attributes of neurons may subject neuronal genes to evolutionary forces favoring net size expansion. This process could be associated with tissue-specific constraints on gene function and/or the evolution of increasingly complex gene regulation in nervous systems.

The evolutionary scaling of cellular traits imposed by the drift barrier

Owing to internal homeostatic mechanisms, cellular traits may experience long periods of stable selective pressures, during which the stochastic forces of drift and mutation conspire to generate variation. However, even in the face of invariant selection, the drift barrier defined by the genetic effective population size, which is negatively associated with organism size, can have a substantial influence on the location and dispersion of the long-term steady-state distribution of mean phenotypes. In addition, for multilocus traits, the multiplicity of alternative, functionally equivalent states can draw mean phenotypes away from selective optima, even in the absence of mutation bias. Using a framework for traits with an additive genetic basis, it is shown that 1) optimal phenotypic states may be only rarely achieved; 2) gradients of mean phenotypes with respect to organism size (i.e., allometric relationships) are likely to be molded by differences in the power of random genetic drift across the tree of life; and 3) for any particular set of population-genetic conditions, significant variation in mean phenotypes may exist among lineages exposed to identical selection pressures. These results provide a potentially useful framework for understanding numerous aspects of cellular diversification and illustrate the risks of interpreting such variation in a purely adaptive framework.

Mitonuclear Compensatory Coevolution

In bilaterian animals, the mitochondrial genome is small, haploid, does not typically recombine, and is subject to accumulation of deleterious alleles via Muller's ratchet. These basic features of the genomic architecture present a paradox: mutational erosion of these genomes should lead to decline in mitochondrial function over time, yet no such decline is observed. Compensatory coevolution, whereby the nuclear genome evolves to compensate for the deleterious alleles in the mitochondrial genome, presents a potential solution to the paradox of Muller's ratchet without loss of function. Here, I review different proposed forms of mitonuclear compensatory coevolution. Empirical evidence from diverse eukaryotic taxa supports the mitonuclear compensatory coevolution hypothesis, but the ubiquity and importance of such compensatory coevolution remains a topic of debate.

Characterizing lineage-specific evolution and the processes driving genomic diversification in chordates

Background

Understanding the origins of genome content has long been a goal of molecular evolution and comparative genomics. By examining genome evolution through the guise of lineage-specific evolution, it is possible to make inferences about the evolutionary events that have given rise to species-specific diversification. Here we characterize the evolutionary trends found in chordate species using The Adaptive Evolution Database (TAED). TAED is a database of phylogenetically indexed gene families designed to detect episodes of directional or diversifying selection across chordates. Gene families within the database have been assessed for lineage-specific estimates of dN/dS and have been reconciled to the chordate species to identify retained duplicates. Gene families have also been mapped to the functional pathways and amino acid changes which occurred on high dN/dS lineages have been mapped to protein structures.

Results

An analysis of this exhaustive database has enabled a characterization of the processes of lineage-specific diversification in chordates. A pathway level enrichment analysis of TAED determined that pathways most commonly found to have elevated rates of evolution included those involved in metabolism, immunity, and cell signaling. An analysis of protein fold presence on proteins, after normalizing for frequency in the database, found common folds such as Rossmann folds, Jelly Roll folds, and TIM barrels were overrepresented on proteins most likely to undergo directional selection. A set of gene families which experience increased numbers of duplications within short evolutionary times are associated with pathways involved in metabolism, olfactory reception, and signaling. An analysis of protein secondary structure indicated more relaxed constraint in β-sheets and stronger constraint on alpha Helices, amidst a general preference for substitutions at exposed sites. Lastly a detailed analysis of the ornithine decarboxylase gene family, a key enzyme in the pathway for polyamine synthesis, revealed lineage-specific evolution along the lineage leading to Cetacea through rapid sequence evolution in a duplicate gene with amino acid substitutions causing active site rearrangement.

Conclusion

Episodes of lineage-specific evolution are frequent throughout chordate species. Both duplication and directional selection have played large roles in the evolution of the phylum. TAED is a powerful tool for facilitating this understanding of lineage-specific evolution.

Within-Population Genome Size Variation is Mediated by Multiple Genomic Elements That Segregate Independently during Meiosis

Within-species variation in genome size has been documented in many animals and plants. Despite its importance for understanding eukaryotic genome diversity, there is only sparse knowledge about how individual-level processes mediate genome size variation in populations. Here, we study a natural population of the rotifer Brachionus asplanchnoidis whose members differ up to 1.9-fold in diploid genome size, but were still able to interbreed and produce viable offspring. We show that genome size is highly heritable and can be artificially selected up or down, but not below a certain basal diploid genome size for this species. Analyses of segregation patterns in haploid males reveal that large genomic elements (several megabases in size) provide the substrate of genome size variation. These elements, and their segregation patterns, explain the generation of new genome size variants, the short-term evolutionary potential of genome size change in populations, and some seemingly paradoxical patterns, like an increase in genome size variation among highly inbred lines. Our study suggests that a conceptual model involving only two variables, 1) a basal genome size of the population, and 2) a vector containing information on additional elements that may increase genome size in this population (size, number, and meiotic segregation behavior), can effectively address most scenarios of short-term evolutionary change of genome size in a population.

Transposable elements contribute to fungal genes and impact fungal lifestyle

The last decade brought a still growing experimental evidence of mobilome impact on host's gene expression. We systematically analysed genomic location of transposable elements (TEs) in 625 publicly available fungal genomes from the NCBI database in order to explore their potential roles in genome evolution and correlation with species' lifestyle. We found that non-autonomous TEs and remnant copies are evenly distributed across genomes. In consequence, they also massively overlap with regions annotated as genes, which suggests a great contribution of TE-derived sequences to host's coding genome. Younger and potentially active TEs cluster with one another away from genic regions. This non-randomness is a sign of either selection against insertion of TEs in gene proximity or target site preference among some types of TEs. Proteins encoded by genes with old transposable elements insertions have significantly less repeat and protein-protein interaction motifs but are richer in enzymatic domains. However, genes only proximal to TEs do not display any functional enrichment. Our findings show that adaptive cases of TE insertion remain a marginal phenomenon, and the overwhelming majority of TEs are evolving neutrally. Eventually, animal-related and pathogenic fungi have more TEs inserted into genes than fungi with other lifestyles. This is the first systematic, kingdom-wide study concerning mobile elements and their genomic neighbourhood. The obtained results should inspire further research concerning the roles TEs played in evolution and how they shape the life we know today.

Multivariate analysis of genomic variables, effective population size, and mutation rate

Objective

The relationship between genomic variables (genome size, gene number, intron size, and intron number) and evolutionary forces has two implications. First, they help to unravel the mechanism underlying genome evolution. Second, they provide a solution to the debate over discrepancy between genome size variation and organismal complexity. Previously, a clear correlation between genomic variables and effective population size and mutation rate (Neu) led to an important hypothesis to consider random genetic drift as a major evolutionary force during evolution of genome size and complexity. But recent reports also support natural selection as the leading evolutionary force. As such, the debate remains unresolved.

Results

Here, we used a multivariate method to explore the relationship between genomic variables and Neu in order to understand the evolution of genome. Previously reported patterns between genomic variables and Neu were not observed in our multivariate study. We found only one association between intron number and Neu, but no relationships were observed between genome size, intron size, gene number, and Neu, suggesting that Neu of the organisms solely does not influence genome evolution. We, therefore, concluded that Neu influences intron evolution, while it may not be the only force that provides mechanistic insights into genome evolution and complexity.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4097-3) contains supplementary material, which is available to authorized users.

Comparative genomics of the major parasitic worms

Parasitic nematodes (roundworms) and platyhelminths (flatworms) cause debilitating chronic infections of humans and animals, decimate crop production and are a major impediment to socioeconomic development. Here we report a broad comparative study of 81 genomes of parasitic and non-parasitic worms. We have identified gene family births and hundreds of expanded gene families at key nodes in the phylogeny that are relevant to parasitism. Examples include gene families that modulate host immune responses, enable parasite migration though host tissues or allow the parasite to feed. We reveal extensive lineage-specific differences in core metabolism and protein families historically targeted for drug development. From an in silico screen, we have identified and prioritized new potential drug targets and compounds for testing. This comparative genomics resource provides a much-needed boost for the research community to understand and combat parasitic worms.

Genome size evolution in the Archaea

What determines variation in genome size, gene content and genetic diversity at the broadest scales across the tree of life? Much of the existing work contrasts eukaryotes with prokaryotes, the latter represented mainly by Bacteria. But any general theory of genome evolution must also account for the Archaea, a diverse and ecologically important group of prokaryotes that represent one of the primary domains of cellular life. Here, we survey the extant diversity of Bacteria and Archaea, and ask whether the general principles of genome evolution deduced from the study of Bacteria and eukaryotes also apply to the archaeal domain. Although Bacteria and Archaea share a common prokaryotic genome architecture, the extant diversity of Bacteria appears to be much higher than that of Archaea. Compared with Archaea, Bacteria also show much greater genome-level specialisation to specific ecological niches, including parasitism and endosymbiosis. The reasons for these differences in long-term diversification rates are unclear, but might be related to fundamental differences in informational processing machineries and cell biological features that may favour archaeal diversification in harsher or more energy-limited environments. Finally, phylogenomic analyses suggest that the first Archaea were anaerobic autotrophs that evolved on the early Earth.

Manifold Routes to a Nucleus

It is widely assumed that there is a clear distinction between eukaryotes, with cell nuclei, and prokaryotes, which lack nuclei. This suggests the evolution of nuclear compartmentation is a singular event. However, emerging knowledge of the diversity of bacterial internal cell structures suggests the picture may not be as black-and-white as previously thought. For instance, some members of the bacterial PVC superphylum appear to have nucleus-like compartmentation, where transcription and translation are physically separated, and some jumbophages have recently been shown to create nucleus-like structures within their Pseudomonad hosts. Moreover, there is also tantalizing metagenomic identification of new Archaea that carry homologs of genes associated with internal cell membrane structure in eukaryotes. All these cases invite comparison with eukaryote cell biology. While the bacterial cases of genetic compartmentation are likely convergent, and thus viewed by many as not germane to the question of eukaryote origins, we argue here that, in addressing the broader question of the evolution of compartmentation, other instances are at least as important: they provide us with a point of comparison which is critical for a more general understanding of both the conditions favoring the emergence of intracellular compartmentation of DNA and the evolutionary consequences of such cellular architecture. Finally, we consider three classes of explanation for the emergence of compartmentation: physical protection, crosstalk avoidance and nonadaptive origins.

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.

Macroevolutionary shift in the size of amphibian genomes and the role of life history and climate

The evolution and great diversity of genome size has been of long-standing interest to biologists, but has seldom been investigated on a broad phylogenetic scale. Here we present a comparative quantitative analysis of factors shaping genome size evolution in amphibians, the extant class of vertebrates with the largest variation in genome size. We find that amphibian genomes have undergone saltations in size, although these are rare and the evolutionary history of genome size in amphibians has otherwise been one of gradual, time-dependent variation (that is, Brownian motion). This macroevolutionary homogeneity is remarkable given the evolutionary and ecological diversity of most other aspects of the natural history of amphibians. Contrary to previous claims, we find no evidence for associations between life cycle complexity and genome size despite the high diversity of reproductive modes and the multiple events of independent evolution of divergent life cycles in the group. Climate (temperature and humidity) affects genome size indirectly, at least in frogs, as a consequence of its effect on premetamorphic developmental period, although directionality of the relationship between developmental period and genome size is not unequivocal.

Interkingdom microbial consortia mechanisms to guide biotechnological applications

Microbial consortia are capable of surviving diverse conditions through the formation of synergistic population-level structures, such as stromatolites, microbial mats and biofilms. Biotechnological applications are poised to capitalize on these unique interactions. However, current artificial co-cultures constructed for societal benefits, including biosynthesis, agriculture and bioremediation, face many challenges to perform as well as natural consortia. Interkingdom microbial consortia tend to be more robust and have higher productivity compared with monocultures and intrakingdom consortia, but the control and design of these diverse artificial consortia have received limited attention. Further, feasible research techniques and instrumentation for comprehensive mechanistic insights have only recently been established for interkingdom microbial communities. Here, we review these recent advances in technology and our current understanding of microbial interaction mechanisms involved in sustaining or developing interkingdom consortia for biotechnological applications. Some of the interactions among members from different kingdoms follow similar mechanisms observed for intrakingdom microbial consortia. However, unique interactions in interkingdom consortia, including endosymbiosis or interkingdom-specific cell-cell interactions, provide improved mitigation to external stresses and inhibitory compounds. Furthermore, antagonistic interactions among interkingdom species can promote fitness, diversification and adaptation, along with the production of beneficial metabolites and enzymes for society. Lastly, we shed light on future research directions to develop study methods at the level of metabolites, genes and meta-omics. These potential research methods could lead to the control and utilization of highly diverse microbial communities.

Linking autoimmunity to the origin of the adaptive immune system

In jawed vertebrates, the adaptive immune system (AIS) cooperates with the innate immune system (IIS) to protect hosts from infections. Although targeting non-self-components, the AIS also generates self-reactive antibodies which, when inadequately counter-selected, can give rise to autoimmune diseases (ADs). ADs are on the rise in western countries. Why haven’t ADs been eliminated during the evolution of a ∼500 million-year old system? And why have they become more frequent in recent decades? Self-recognition is an attribute of the phylogenetically more ancient IIS and empirical data compellingly show that some self-reactive antibodies, which are classifiable as elements of the IIS rather then the AIS, may protect from (rather than cause) ADs. Here, we propose that the IIS’s self-recognition system originally fathered the AIS and, as a consequence of this relationship, its activity is dampened in hygienic environments. Rather than a mere breakdown or failure of the mechanisms of self-tolerance, ADs might thus arise from architectural constraints.

Polymorphic duplicate genes and persistent non-coding sequences reveal heterogeneous patterns of mitochondrial DNA loss in salamanders

Background

Mitochondria are the site of the citric acid cycle and oxidative phosphorylation (OXPHOS). In metazoans, the mitochondrial genome is a small, circular molecule averaging 16.5 kb in length. Despite evolutionarily conserved gene content, metazoan mitochondrial genomes show a diversity of gene orders most commonly explained by the duplication-random loss (DRL) model. In the DRL model, (1) a sequence of genes is duplicated in tandem, (2) one paralog sustains a loss-of-function mutation, resulting in selection to retain the other copy, and (3) the non-functional paralog is eventually deleted from the genome. Despite its apparent role in generating mitochondrial gene order diversity, little is known about the tempo and mode of random gene loss after duplication events. Here, we determine mitochondrial gene order across the salamander genus Aneides, which was previously shown to include at least two DRL-mediated rearrangement events. We then analyze these gene orders in a phylogenetic context to reveal patterns of DNA loss after mitochondrial gene duplication.

Results

We identified two separate duplication events that resulted in mitochondrial gene rearrangements in Aneides; one occurred at the base of the clade tens of millions of years ago, while the other occurred much more recently (i.e. within a single species), resulting in gene order polymorphism and paralogs that are readily identifiable. We demonstrate that near-complete removal of duplicate rRNA genes has occurred since the recent duplication event, whereas duplicate protein-coding genes persist as pseudogenes and duplicate tRNAs persist as functionally intact paralogs. In addition, we show that non-coding DNA duplicated at the base of the clade has persisted across species for tens of millions of years.

Conclusions

The evolutionary history of the mitochondrial genome, from its inception as a bacterial endosymbiont, includes massive genomic reduction. Consistent with this overall trend, selection for efficiency of mitochondrial replication and transcription has been hypothesized to favor elimination of extra sequence. Our results, however, suggest that there may be no strong disadvantage to extraneous sequences in salamander mitochondrial genomes, although duplicate rRNA genes may be deleterious.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4358-2) contains supplementary material, which is available to authorized users.

Comparative Genomics of an Unusual Biogeographic Disjunction in the Cotton Tribe (Gossypieae) Yields Insights into Genome Downsizing

Long-distance insular dispersal is associated with divergence and speciation because of founder effects and strong genetic drift. The cotton tribe (Gossypieae) has experienced multiple transoceanic dispersals, generating an aggregate geographic range that encompasses much of the tropics and subtropics worldwide. Two genera in the Gossypieae, Kokia and Gossypioides, exhibit a remarkable geographic disjunction, being restricted to the Hawaiian Islands and Madagascar/East Africa, respectively. We assembled and use de novo genome sequences to address questions regarding the divergence of these two genera from each other and from their sister-group, Gossypium. In addition, we explore processes underlying the genome downsizing that characterizes Kokia and Gossypioides relative to other genera in the tribe. Using 13,000 gene orthologs and synonymous substitution rates, we show that the two disjuncts last shared a common ancestor ∼5 Ma, or half as long ago as their divergence from Gossypium. We report relative stasis in the transposable element fraction. In comparison to Gossypium, there is loss of ∼30% of the gene content in the two disjunct genera and a history of genome-wide accumulation of deletions. In both genera, there is a genome-wide bias toward deletions over insertions, and the number of gene losses exceeds the number of gains by ∼2- to 4-fold. The genomic analyses presented here elucidate genomic consequences of the demographic and biogeographic history of these closest relatives of Gossypium, and enhance their value as phylogenetic outgroups.

Strong phylogenetic inertia on genome size and transposable element content among 26 species of flies

While the evolutionary mechanisms driving eukaryote genome size evolution are still debated, repeated element content appears to be crucial. Here, we reconstructed the phylogeny and identified repeats in the genome of 26 Drosophila exhibiting a twofold variation in genome size. The content in transposable elements (TEs) is highly correlated to genome size evolution among these closely related species. We detected a strong phylogenetic signal on the evolution of both genome size and TE content, and a genome contraction in the Drosophila melanogaster subgroup.

Emerging and evolving concepts in gene essentiality

Gene essentiality is a founding concept of genetics with important implications in both fundamental and applied research. Multiple screens have been performed over the years in bacteria, yeasts, animals and more recently in human cells to identify essential genes. A mounting body of evidence suggests that gene essentiality, rather than being a static and binary property, is both context dependent and evolvable in all kingdoms of life. This concept of a non-absolute nature of gene essentiality changes our fundamental understanding of essential biological processes and could directly affect future treatment strategies for cancer and infectious diseases.

Evolution of intrinsic disorder in eukaryotic proteins

Conformational flexibility conferred though regions of intrinsic structural disorder allows proteins to behave as dynamic molecules. While it is well-known that intrinsically disordered regions can undergo disorder-to-order transitions in real-time as part of their function, we also are beginning to learn more about the dynamics of disorder-to-order transitions along evolutionary time-scales. Intrinsically disordered regions endow proteins with functional promiscuity, which is further enhanced by the ability of some of these regions to undergo real-time disorder-to-order transitions. Disorder content affects gene retention after whole genome duplication, but it is not necessarily conserved. Altered patterns of disorder resulting from evolutionary disorder-to-order transitions indicate that disorder evolves to modify function through refining stability, regulation, and interactions. Here, we review the evolution of intrinsically disordered regions in eukaryotic proteins. We discuss the interplay between secondary structure and disorder on evolutionary time-scales, the importance of disorder for eukaryotic proteome expansion and functional divergence, and the evolutionary dynamics of disorder.

The Evolution of Bacterial Genome Architecture

The genome architecture of bacteria and eukaryotes evolves in opposite directions when subject to genetic drift, a difference that can be ascribed to the fact that bacteria exhibit a mutational bias that deletes superfluous sequences, whereas eukaryotes are biased toward large insertions. Expansion of eukaryotic genomes occurs through the addition of non-functional sequences, such as repetitive sequences and transposable elements, whereas variation in bacterial genome size is largely due to the acquisition and loss of functional accessory genes. These properties create the situation in which eukaryotes with very similar numbers of genes can have vastly different genome sizes, while in bacteria, gene number scales linearly with genome size. Some bacterial genomes, however, particularly those of species that undergo bottlenecks due to recent association with hosts, accumulate pseudogenes and mobile elements, conferring them a low gene content relative to their genome size. These non-functional sequences are gradually eroded and eliminated after long-term association with hosts, with the result that obligate symbionts have the smallest genomes of any cellular organism. The architecture of bacterial genomes is shaped by complex and diverse processes, but for most bacterial species, genome size is governed by a non-adaptive process, i.e., genetic drift coupled with a mutational bias toward deletions. Thus, bacteria with small effective population sizes typically have the smallest genomes. Some marine bacteria counter this near-universal trend: despite having immense population sizes, selection, not drift, acts to reduce genome size in response to metabolic constraints in their nutrient-limited environment.

Large-Scale Comparative Analysis Reveals the Mechanisms Driving Plastomic Compaction, Reduction, and Inversions in Conifers II (Cupressophytes)

Conifers II (cupressophytes), comprising about 400 tree species in five families, are the most diverse group of living gymnosperms. Their plastid genomes (plastomes) are highly variable in size and organization, but such variation has never been systematically studied. In this study, we assessed the potential mechanisms underlying the evolution of cupressophyte plastomes. We analyzed the plastomes of 24 representative genera in all of the five cupressophyte families, focusing on their variation in size, noncoding DNA content, and nucleotide substitution rates. Using a tree-based method, we further inferred the ancestral plastomic organizations of internal nodes and evaluated the inversions across the evolutionary history of cupressophytes. Our data showed that variation in plastome size is statistically associated with the dynamics of noncoding DNA content, which results in different degrees of plastomic compactness among the cupressophyte families. The degrees of plastomic inversions also vary among the families, with the number of inversions per genus ranging from 0 in Araucariaceae to 1.27 in Cupressaceae. In addition, we demonstrated that synonymous substitution rates are significantly correlated with plastome size as well as degree of inversions. These data suggest that in cupressophytes, mutation rates play a critical role in driving the evolution of plastomic size while plastomic inversions evolve in a neutral manner.