Mitochondrial (mt) and nuclear-encoded proteins are integrated in aerobic respiration, requiring co-functionality among gene products from fundamentally different genomes. Different evolutionary rates, inheritance mechanisms, and selection pressures set the stage for incompatibilities between interacting products of the two genomes. The mitonuclear coevolution hypothesis posits that incompatibilities may be avoided if evolution in one genome selects for complementary changes in interacting genes encoded by the other genome. Nuclear compensation, in which deleterious mtDNA changes are offset by compensatory nuclear changes, is often invoked as the primary mechanism for mitonuclear coevolution. Yet, direct evidence supporting nuclear compensation is rare. Here, we used data from 58 mammalian species representing eight orders to show strong correlations between evolutionary rates of mt and nuclear-encoded mt-targeted (N-mt) proteins, but not between mt and non-mt-targeted nuclear proteins, providing strong support for mitonuclear coevolution across mammals. N-mt genes with direct mt interactions also showed the strongest correlations. Although most N-mt genes had elevated dN/dS ratios compared to mt genes (as predicted under nuclear compensation), N-mt sites in close contact with mt proteins were not overrepresented for signs of positive selection compared to noncontact N-mt sites (contrary to predictions of nuclear compensation). Furthermore, temporal patterns of N-mt and mt amino acid substitutions did not support predictions of nuclear compensation, even in positively selected, functionally important residues with direct mitonuclear contacts. Overall, our results strongly support mitonuclear coevolution across ∼170 million years of mammalian evolution but fail to support nuclear compensation as the major mode of mitonuclear coevolution.

Ancestral sequence reconstruction (ASR) infers predicted ancestral states for sites within sequences and can constrain the functions and properties of ancestors of extant protein families. Here, we compare the likely sequences of inferred nitrogenase ancestors to extant nitrogenase sequence diversity. We show that the most-likely combinations of ancestral states for key substrate channel residues are not represented in extant sequence space, and rarely found within a more broadly defined physiochemical space-supporting that the earliest ancestors of extant nitrogenases likely had alternative substrate channel composition. These differences may indicate differing environmental selection pressures acting on nitrogenase substrate specificity in ancient environments. These results highlight ASR's potential as an in silico tool for developing hypotheses about ancestral enzyme functions, as well as improving hypothesis testing through more targeted in vitro and in vivo experiments.

The inference of genome rearrangement events has been extensively studied, as they play a major role in molecular evolution. However, probabilistic evolutionary models that explicitly imitate the evolutionary dynamics of such events, as well as methods to infer model parameters, are yet to be fully utilized. Here, we developed a probabilistic approach to infer genome rearrangement rate parameters using an Approximate Bayesian Computation (ABC) framework. We developed two genome rearrangement models, a basic model, which accounts for genomic changes in gene order, and a more sophisticated one which also accounts for changes in chromosome number. We characterized the ABC inference accuracy using simulations and applied our methodology to both prokaryotic and eukaryotic empirical datasets. Knowledge of genome-rearrangement rates can help elucidate their role in evolution as well as help simulate genomes with evolutionary dynamics that reflect empirical genomes.

Gene-by-environment interactions play a crucial role in horizontal gene transfer by affecting how the transferred genes alter host fitness. However, how the environment modulates the fitness effect of transferred genes has not been tested systematically in an experimental study. We adapted a high-throughput technique for obtaining very precise estimates of bacterial fitness, in order to measure the fitness effects of 44 orthologs transferred from Salmonella Typhimurium to Escherichia coli in six physiologically relevant environments. We found that the fitness effects of individual genes were highly dependent on the environment, while the distributions of fitness effects across genes were not, with all tested environments resulting in distributions of same shape and spread. Furthermore, the extent to which the fitness effects of a gene varied between environments depended on the average fitness effect of that gene across all environments, with nearly neutral and nearly lethal genes having more consistent fitness effects across all environments compared to deleterious genes. Put together, our results reveal the unpredictable nature of how environmental conditions impact the fitness effects of each individual gene. At the same time, distributions of fitness effects across environments exhibit consistent features, pointing to the generalizability of factors that shape horizontal gene transfer of orthologous genes.

Range expansions have been common in the history of most species. Serial founder effects and subsequent population growth at expansion fronts typically lead to a loss of genomic diversity along the expansion axis. A frequent consequence is the phenomenon of "gene surfing," where variants located near the expanding front can reach high frequencies or even fix in newly colonized territories. Although gene surfing events have been characterized thoroughly for a specific locus, their effects on linked genomic regions and the overall patterns of genomic diversity have been little investigated. In this study, we simulated the evolution of whole genomes during several types of 1D and 2D range expansions differing by the extent of migration, founder events, and recombination rates. We focused on the characterization of local dips of diversity, or "troughs," taken as a proxy for surfing events. We find that, for a given recombination rate, once we consider the amount of diversity lost since the beginning of the expansion, it is possible to predict the initial evolution of trough density and their average width irrespective of the expansion condition. Furthermore, when recombination rates vary across the genome, we find that troughs are over-represented in regions of low recombination. Therefore, range expansions can leave local and global genomic signatures often interpreted as evidence of past selective events. Given the generality of our results, they could be used as a null model for species having gone through recent expansions, and thus be helpful to correctly interpret many evolutionary biology studies.

Topologically associating domains (TADs) are thought to play an important role in preventing gene misexpression by spatially constraining enhancer-promoter contacts. The deleterious nature of gene misexpression implies that TADs should, therefore, be conserved among related species. Several early studies comparing chromosome conformation between species reported high levels of TAD conservation; however, more recent studies have questioned these results. Furthermore, recent work suggests that TAD reorganization is not associated with extensive changes in gene expression. Here, we investigate the evolutionary conservation of TADs among 11 species of Drosophila. We use Hi-C data to identify TADs in each species and employ a comparative phylogenetic approach to derive empirical estimates of the rate of TAD evolution. Surprisingly, we find that TADs evolve rapidly. However, we also find that the rate of evolution depends on the chromatin state of the TAD, with TADs enriched for developmentally regulated chromatin evolving significantly slower than TADs enriched for broadly expressed, active chromatin. We also find that, after controlling for differences in chromatin state, highly conserved TADs do not exhibit higher levels of gene expression constraint. These results suggest that, in general, most TADs evolve rapidly and their divergence is not associated with widespread changes in gene expression. However, higher levels of evolutionary conservation and gene expression constraints in TADs enriched for developmentally regulated chromatin suggest that these TAD subtypes may be more important for regulating gene expression, likely due to the larger number of long-distance enhancer-promoter contacts associated with developmental genes.

S* is a widely used statistic for detecting archaic admixture from population genetic data. Previous studies used freezing-archer to apply S*, which is only directly applicable to the specific case of Neanderthal and Denisovan introgression in Papuans. Here, we implemented sstar for a more general purpose. Compared with several tools, including SPrime, SkovHMM, and ArchaicSeeker2.0, for detecting introgressed fragments with simulations, our results suggest that sstar is robust to differences in demographic models, including ghost introgression and two-source introgression. We believe sstar will be a useful tool for detecting introgressed fragments in various scenarios and in non-human species.

The selection of the optimal substitution model of molecular evolution imposes a high computational burden for long sequence alignments in phylogenomics. We discovered that the analysis of multiple tiny subsamples of site patterns from a full sequence alignment recovers the correct optimal substitution model when sites in the subsample are upsampled to match the total number of sites in the full alignment. The computational costs of maximum-likelihood analyses are reduced by orders of magnitude in the subsample-upsample (SU) approach because the upsampled alignment contains only a small fraction of all site patterns. We present an adaptive protocol, ModelTamer, that implements the new SU approach and automatically selects subsamples to estimate optimal models reliably. ModelTamer selects models hundreds to thousands of times faster than the full data analysis while needing megabytes rather than gigabytes of computer memory.

Maize is a staple food of smallholder farmers living in highland regions up to 4,000 m above sea level worldwide. Mexican and South American highlands are two major highland maize growing regions, and population genetic data suggest the maize's adaptation to these regions occurred largely independently, providing a case study for convergent evolution. To better understand the mechanistic basis of highland adaptation, we crossed maize landraces from 108 highland and lowland sites of Mexico and South America with the inbred line B73 to produce F1 hybrids and grew them in both highland and lowland sites in Mexico. We identified thousands of genes with divergent expression between highland and lowland populations. Hundreds of these genes show patterns of convergent evolution between Mexico and South America. To dissect the genetic architecture of the divergent gene expression, we developed a novel allele-specific expression analysis pipeline to detect genes with divergent functional cis-regulatory variation between highland and lowland populations. We identified hundreds of genes with divergent cis-regulation between highland and lowland landrace alleles, with 20 in common between regions, further suggesting convergence in the genes underlying highland adaptation. Further analyses suggest multiple mechanisms contribute to this convergence in gene regulation. Although the vast majority of evolutionary changes associated with highland adaptation were region specific, our findings highlight an important role for convergence at the gene expression and gene regulation levels as well.

Karyotype refers to the configuration of the genome into a set of chromosomes. The karyotype difference between species is expected to impede various biological processes, such as chromosome segregation and meiotic chromosome pairing, potentially contributing to incompatibility. Karyotypes can rapidly change between closely related species and even among populations of the same species. However, the forces driving karyotype evolution are poorly understood. Here we describe a unique karyotype of a Drosophila melanogaster strain isolated from the Seychelles archipelago. This strain has lost the ribosomal DNA (rDNA) locus on the X chromosome. Because the Y chromosome is the only other rDNA-bearing chromosome, all females carry at least one Y chromosome as the source of rDNA. Interestingly, we found that the strain also carries a truncated Y chromosome (YS) that is stably maintained in the population despite its inability to support male fertility. Our modeling and cytological analysis suggest that the Y chromosome has a larger negative impact on female fitness than the YS chromosome. Moreover, we generated an independent strain that lacks X rDNA and has a karyotype of XXY females and XY males. This strain quickly evolved multiple karyotypes: two new truncated Y chromosomes (similar to YS), as well as two independent X chromosome fusions that contain the Y-derived rDNA fragment, eliminating females' dependence on the Y chromosome. Considering that Robertsonian fusions frequently occur at rDNA loci in humans, we propose that rDNA loci instability may be one of driving forces of karyotype evolution.

MicroRNAs (miRNAs) are fast evolving endogenous small RNAs that regulate organism function and behavior in both animals and plants. Although models for de novo miRNA biogenesis have been proposed, the genomic mechanisms driving swift diversification of the miRNA repertoires in plants remain elusive. Here, by comprehensively analyzing 21 phylogenetically representative plant species, ranging from green algae to angiosperms, we systematically identified de novo miRNA events associated with 8,649 miRNA loci. We found that 399 (4.6%), 466 (5.4%), and 1,402 (16.2%) miRNAs were derived from inverted gene duplication events, long terminal repeats of retrotransposons, and miniature inverted-repeat transposable elements (MITEs), respectively. Among the miRNAs of these origins, MITEs, especially those belonging to the Mutator, Tc1/Mariner, and PIF/Harbinger superfamilies, were the predominant genomic source for de novo miRNAs in the 15 examined angiosperms but not in the six non-angiosperms. Our data further illustrated a transposition-transcription process by which MITEs are converted into new miRNAs (termed MITE-miRNAs) whereby properly sized MITEs are transcribed and therefore become potential substrates for the miRNA processing machinery by transposing into introns of active genes. By analyzing the 58,038 putative target genes for the 8,095 miRNAs, we found that the target genes of MITE-miRNAs were preferentially associated with response to environmental stimuli such as temperature, suggesting that MITE-miRNAs are pertinent to plant adaptation. Collectively, these findings demonstrate that molecular conversion of MITEs is a genomic mechanism leading to rapid and continuous changes to the miRNA repertoires in angiosperm.

Cytonuclear coordination between biparental-nuclear genomes and uniparental-cytoplasmic organellar genomes in plants is often resolved by genetic and transcriptional cytonuclear responses. Whether this mechanism also acts in allopolyploid members of other kingdoms is not clear. Additionally, cytonuclear coordination of interleaved allopolyploid cells/individuals within the same population is underexplored. The yeast Saccharomyces pastorianus provides the opportunity to explore cytonuclear coevolution during different growth stages and from novel dimensions. Using S. pastorianus cells from multiple growth stages in the same environment, we show that nuclear mitochondria-targeted genes have undergone both asymmetric gene conversion and growth stage-specific biased expression favoring genes from the mitochondrial genome donor (Saccharomyces eubayanus). Our results suggest that cytonuclear coordination in allopolyploid lager yeast species entails an orchestrated and compensatory genetic and transcriptional evolutionary regulatory shift. The common as well as unique properties of cytonuclear coordination underlying allopolyploidy between unicellular yeasts and higher plants offers novel insights into mechanisms of cytonuclear evolution associated with allopolyploid speciation.

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

Proteins need to selectively interact with specific targets among a multitude of similar molecules in the cell. However, despite a firm physical understanding of binding interactions, we lack a general theory of how proteins evolve high specificity. Here, we present such a model that combines chemistry, mechanics, and genetics and explains how their interplay governs the evolution of specific protein-ligand interactions. The model shows that there are many routes to achieving molecular discrimination-by varying degrees of flexibility and shape/chemistry complementarity-but the key ingredient is precision. Harder discrimination tasks require more collective and precise coaction of structure, forces, and movements. Proteins can achieve this through correlated mutations extending far from a binding site, which fine-tune the localized interaction with the ligand. Thus, the solution of more complicated tasks is enabled by increasing the protein size, and proteins become more evolvable and robust when they are larger than the bare minimum required for discrimination. The model makes testable, specific predictions about the role of flexibility and shape mismatch in discrimination, and how evolution can independently tune affinity and specificity. Thus, the proposed theory of specific binding addresses the natural question of "why are proteins so big?". A possible answer is that molecular discrimination is often a hard task best performed by adding more layers to the protein.

Uncoupling protein 1 (UCP1) is an essential protein in the mitochondrial inner membrane that mediates nonshivering thermogenesis (NST) and plays an important role in thermoregulation and fat deposition. However, the relationship between the evolution of UCP1 and fat deposition in the blubber layer in cetaceans remains unclear. Here, frameshift mutations, premature termination, and relaxed selection pressure (ω = 0.9557, P < 0.05) were detected in UCP1 in cetaceans, suggesting that UCP1 was inactivated during cetacean evolution. By time estimation, it was found that the inactivation of UCP1 in cetaceans occurred between 53.1 and 50.2 Ma. However, combined with findings from immunohistochemical analysis of the blubber layer of the Yangtze finless porpoise and in vitro functional assays, a premature termination of cetacean UCP1 resulted in a reduction of UCP1-mediated NST capacity (about 50%) and lipolytic capacity (about 40%), both of which were beneficial to maintain blubber layer and body temperature without excessive fat consumption. This study provides new insights into the molecular mechanisms of the blubber thickening in cetaceans and highlights the importance of UCP1 attenuation in cetaceans for secondary aquatic adaptation.

Despite a century of genetic analysis, the evolutionary processes that have generated the patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic variation is partitioned within its putative ancestral range in Southern Africa remains unresolved. Here, we study patterns of population genetic structure, admixture, and the spatial structuring of candidate incompatibility alleles across a global sample, including 223 new accessions, predominantly from remote regions in Southern Africa. We identify nine major ancestries, six that primarily occur in Africa and one that has not been previously described. We find evidence for both contemporary and historical admixture between ancestries, with admixture rates varying both within and between continents. For example, while previous work has highlighted an admixture zone between broadly defined African and European ancestries in the Caribbean and southeastern USA, we identify West African ancestry as the most likely African contributor. Moreover, loci showing the strongest signal of introgression between West Africa and the Caribbean/southeastern USA include several genes relating to neurological development and male courtship behavior, in line with previous work showing shared mating behaviors between these regions. Finally, while we hypothesized that potential incompatibility loci may contribute to population genetic structure across the range of D. melanogaster; these loci are, on average, not highly differentiated between ancestries. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the partitioning of diversity across its range.

Geographically distinct populations can adapt to the temperature conditions of their local environment, leading to temperature-dependent fitness differences between populations. Consistent with local adaptation, phylogeographically distinct Caenorhabditis briggsae nematodes show distinct fitness responses to temperature. The genetic mechanisms underlying local adaptation, however, remain unresolved. To investigate the potential role of small noncoding RNAs in genotype-specific responses to temperature, we quantified small RNA expression using high-throughput sequencing of C. briggsae nematodes from tropical and temperate strain genotypes reared under three temperature conditions (14 °C, 20 °C, and 30 C). Strains representing both tropical and temperate regions showed significantly lower expression of PIWI-interacting RNAs (piRNAs) at high temperatures, primarily mapping to a large ∼7 Mb long piRNA cluster on chromosome IV. We also documented decreased expression of 22G-RNAs antisense to protein-coding genes and other genomic features at high rearing temperatures for the thermally-intolerant temperate strain genotype, but not for the tropical strain genotype. Reduced 22G-RNA expression was widespread along chromosomes and among feature types, indicative of a genome-wide response. Targets of the EGO-1/CSR-1 22G-RNA pathway were most strongly impacted compared with other 22G-RNA pathways, implicating the CSR-1 Argonaute and its RNA-dependent RNA polymerase EGO-1 in the genotype-dependent modulation of C. briggsae 22G-RNAs under chronic thermal stress. Our work suggests that gene regulation via small RNAs may be an important contributor to the evolution of local adaptations.

The appearance of genomic variations such as loss of heterozygosity (LOH) has a significant impact on phenotypic diversity observed in a population. Recent large-scale yeast population genomic surveys have shown a high frequency of these events in natural isolates and more particularly in polyploids. However, the frequency, extent, and spectrum of LOH in polyploid organisms have never been explored and are poorly characterized to date. Here, we accumulated 5,163 LOH events over 1,875 generations in 76 mutation accumulation (MA) lines comprising nine natural heterozygous diploid, triploid, and tetraploid natural S. cerevisiae isolates from different ecological and geographical origins. We found that the rate and spectrum of LOH are variable across ploidy levels. Of the total accumulated LOH events, 8.5%, 21%, and 70.5% of them were found in diploid, triploid, and tetraploid MA lines, respectively. Our results clearly show that the frequency of generated LOH events increases with ploidy level. In fact, the cumulative LOH rates were estimated to be 9.3 × 10-3, 2.2 × 10-2, and 8.4 × 10-2 events per division for diploids, triploids, and tetraploids, respectively. In addition, a clear bias toward the accumulation of interstitial and short LOH tracts is observed in triploids and tetraploids compared with diploids. The variation of the frequency and spectrum of LOH events across ploidy level could be related to the genomic instability, characterizing higher ploidy isolates.

Traits shared among distantly related lineages are indicators of common evolutionary constraints, at the ecological, physiological, or molecular level. Here, we show that the vertebral stripe, a cryptic color pattern, has evolved hundreds of times in the evolutionary history of anurans (frogs and toads) and is favored in terrestrial habitats. Using a genome-wide association study, we demonstrate that variation near the Agouti signaling protein gene (ASIP) is responsible for the different vertebral stripe phenotypes in the African grass frog Ptychadena robeensis. RNAseq and real-time quantitative PCR revealed that differential expression of the gene and an adjacent long non-coding RNA is linked to patterning in this species. Surprisingly, and although the stripe phenotypes are shared with closely related species, we found that the P. robeensis alleles are private to the species and unlikely to evolve under long-term balancing selection, thus indicating that the vertebral stripe phenotypes result from parallel evolution within the group. Our findings demonstrate that this cryptic color pattern evolved rapidly and recurrently in terrestrial anurans, and therefore constitutes an ideal system to study repeated evolution.

Adaptation from standing genetic variation is an important process underlying evolution in natural populations, but we rarely get the opportunity to observe the dynamics of fitness and genomic changes in real time. Here, we used experimental evolution and Pool-Seq to track the phenotypic and genomic changes of genetically diverse asexual populations of the yeast Saccharomyces cerevisiae in four environments with different fitness costs. We found that populations rapidly and in parallel increased in fitness in stressful environments. In contrast, allele frequencies showed a range of trajectories, with some populations fixing all their ancestral variation in <30 generations and others maintaining diversity across hundreds of generations. We detected parallelism at the genomic level (involving genes, pathways, and aneuploidies) within and between environments, with idiosyncratic changes recurring in the environments with higher stress. In particular, we observed a tendency of becoming haploid-like in one environment, whereas the populations of another environment showed low overall parallelism driven by standing genetic variation despite high selective pressure. This work highlights the interplay between standing genetic variation and the influx of de novo mutations in populations adapting to a range of selective pressures with different underlying trait architectures, advancing our understanding of the constraints and drivers of adaptation.

Type II DNA topoisomerases regulate topology by double-stranded DNA cleavage and ligation. The TopoVI family of DNA topoisomerase, first identified and biochemically characterized in Archaea, represents, with TopoVIII and mini-A, the type IIB family. TopoVI has several intriguing features in terms of function and evolution. TopoVI has been identified in some eukaryotes, and a global view is lacking to understand its evolutionary pattern. In addition, in eukaryotes, the two TopoVI subunits (TopoVIA and TopoVIB) have been duplicated and have evolved to give rise to Spo11 and TopoVIBL, forming TopoVI-like (TopoVIL), a complex essential for generating DNA breaks that initiate homologous recombination during meiosis. TopoVIL is essential for sexual reproduction. How the TopoVI subunits have evolved to ensure this meiotic function is unclear. Here, we investigated the phylogenetic conservation of TopoVI and TopoVIL. We demonstrate that BIN4 and RHL1, potentially interacting with TopoVIB, have co-evolved with TopoVI. Based on model structures, this observation supports the hypothesis for a role of TopoVI in decatenation of replicated chromatids and predicts that in eukaryotes the TopoVI catalytic complex includes BIN4 and RHL1. For TopoVIL, the phylogenetic analysis of Spo11, which is highly conserved among Eukarya, highlighted a eukaryal-specific N-terminal domain that may be important for its regulation. Conversely, TopoVIBL was poorly conserved, giving rise to ATP hydrolysis-mutated or -truncated protein variants, or was undetected in some species. This remarkable plasticity of TopoVIBL provides important information for the activity and function of TopoVIL during meiosis.

Neotropical Heliconius butterflies are well known for their intricate behaviors and multiple instances of incipient speciation. Chemosensing plays a fundamental role in the life history of these groups of butterflies and in the establishment of reproductive isolation. However, chemical communication involves synergistic sensory and accessory functions, and it remains challenging to investigate the molecular mechanisms underlying behavioral differences. Here, we examine the gene expression profiles and genomic divergence of three sensory tissues (antennae, legs, and mouthparts) between sexes (females and males) and life stages (different adult stages) in two hybridizing butterflies, Heliconius melpomene and Heliconius cydno. By integrating comparative transcriptomic and population genomic approaches, we found evidence of widespread gene expression divergence, supporting a crucial role of sensory tissues in the establishment of species barriers. We also show that sensory diversification increases in a manner consistent with evolutionary divergence based on comparison with the more distantly related species Heliconius charithonia. The findings of our study strongly support the unique chemosensory function of antennae in all three species, the importance of the Z chromosome in interspecific divergence, and the nonnegligible role of nonchemosensory genes in the divergence of chemosensory tissues. Collectively, our results provide a genome-wide illustration of diversification in the chemosensory system under incomplete reproductive isolation, revealing strong molecular separation in the early stage of speciation. Here, we provide a unique perspective and relevant view of the genetic architecture (sensory and accessory functions) of chemosensing beyond the classic chemosensory gene families, leading to a better understanding of the magnitude and complexity of molecular changes in sensory tissues that contribute to the establishment of reproductive isolation and speciation.

A vast body of studies is available that describe age-dependent gene expression in relation to aging in a number of different model species. These data were obtained from animals kept in conditions with reduced environmental challenges, abundant food, and deprivation of natural sensory stimulation. Here, we compared wild- and captive aging in the short-lived turquoise killifish (Nothobranchius furzeri). These fish inhabit temporary ponds in the African savannah. When the ponds are flooded, eggs hatch synchronously, enabling a precise timing of their individual and population age. We collected the brains of wild fish of different ages and quantified the global age-dependent regulation of transcripts using RNAseq. A major difference between captive and wild populations is that wild populations had unlimited access to food and hence grew to larger sizes and reached asymptotic size more rapidly, enabling the analysis of age-dependent gene expression without the confounding effect of adult brain growth. We found that the majority of differentially expressed genes show the same direction of regulation in wild and captive populations. However, a number of genes were regulated in opposite direction. Genes downregulated in the wild and upregulated in captivity were enriched for terms related to neuronal communication. Genes upregulated in the wild and downregulated in captive conditions were enriched in terms related to DNA replication. Finally, the rate of age-dependent gene regulation was higher in wild animals, suggesting a phenomenon of accelerated aging.

Key innovations enable access to new adaptive zones and are often linked to increased species diversification. As such, innovations have attracted much attention, yet their concrete consequences on the subsequent evolutionary trajectory and diversification of the bearing lineages remain unclear. Water striders and relatives (Hemiptera: Heteroptera: Gerromorpha) represent a monophyletic lineage of insects that transitioned to live on the water-air interface and that diversified to occupy ponds, puddles, streams, mangroves and even oceans. This lineage offers an excellent model to study the patterns and processes underlying species diversification following the conquest of new adaptive zones. However, such studies require a reliable and comprehensive phylogeny of the infraorder. Based on whole transcriptomic datasets of 97 species and fossil records, we reconstructed a new phylogeny of the Gerromorpha that resolved inconsistencies and uncovered strong support for previously unknown relationships between some important taxa. We then used this phylogeny to reconstruct the ancestral state of a set of adaptations associated with water surface invasion (fluid locomotion, dispersal and transition to saline waters) and sexual dimorphism. Our results uncovered important patterns and dynamics of phenotypic evolution, revealing how the initial event of water surface invasion enabled multiple subsequent transitions to new adaptive zones on the water surfaces. This phylogeny and the associated transcriptomic datasets constitute highly valuable resources, making Gerromorpha an attractive model lineage to study phenotypic evolution.

To investigate correlation between the ameloblastin (Ambn) amino acid sequence and the emergence of prismatic enamel, a notable event in the evolution of ectodermal hard tissues, we analyzed Ambn sequences of 53 species for which enamel microstructures have been previously reported. We found that a potential amphipathic helix (AH) within the sequence encoded by Exon 5 of Ambn appeared in species with prismatic enamel, with a few exceptions. We studied this correlation by investigating synthetic peptides from different species. A blue shift in fluorescence spectroscopy suggested that the peptides derived from mammalian Ambn interacted with liposomes. A downward shift at 222 nm in circular dichroism spectroscopy of the peptides in the presence of liposomes suggested that the peptides of mammals with prismatic enamel underwent a transition from disordered to helical structure. The peptides of species without prismatic enamel did not show similar secondary structural changes in the presence of liposomes. Peptides of mammals with prismatic enamel caused liposome leakage and inhibited LS8 and ALC cell spreading regulated by full-length Ambn. RT-PCR showed that AH is involved in Ambn's regulation of cell polarization genes: Vangl2, Vangl1, Prickle1, ROCK1, ROCK2, and Par3. Our comprehensive sequence analysis clearly demonstrates that AH motif is closely related to the emergence of enamel prismatic structure, providing insight into the evolution of complex enamel microstructure. We speculate that the AH motif evolved in mammals to interact with cell membrane, triggering signaling pathways required for specific changes in cell morphology associated with the formation of enamel prismatic structure.