An important goal of evolutionary genomics is to identify genomic regions whose substitution rates differ among lineages. For example, genomic regions experiencing accelerated molecular evolution in some lineages may provide insight into links between genotype to phenotype. Several comparative genomics methods have been developed to identify genomic accelerations between species, including a Bayesian method called PhyloAcc, which models shifts in substitution rate in multiple target lineages on a phylogeny. However, few methods consider the possibility of discordance between the trees of individual loci and the species tree due to incomplete lineage sorting, which might cause false positives. Here we present PhyloAcc-GT, which extends PhyloAcc by modeling gene tree heterogeneity. Given a species tree, we adopt the multispecies coalescent model as the prior distribution of gene trees, use Markov chain Monte Carlo (MCMC) for inference, and design novel MCMC moves to sample gene trees efficiently. Through extensive simulations, we show that PhyloAcc-GT outperforms PhyloAcc and other methods in identifying target-lineage-specific accelerations and detecting complex patterns of rate shifts, and is robust to specification of population size parameters. PhyloAcc-GT is usually more conservative than PhyloAcc in calling convergent rate shifts because it identifies more accelerations on ancestral than on terminal branches. We apply PhyloAcc-GT to two examples of convergent evolution: flightlessness in ratites and marine mammal adaptations, and show that PhyloAcc-GT is a robust tool to identify shifts in substitution rate associated with specific target lineages while accounting for incomplete lineage sorting.

Uric acid is the main means of nitrogen excretion in uricotelic vertebrates (birds and reptiles) and the end product of purine catabolism in humans and a few other mammals. While uricase is inactivated in mammals unable to degrade urate, the presence of orthologous genes without inactivating mutations in avian and reptilian genomes is unexplained. Here we show that the Gallus gallus gene we name cysteine-rich urate oxidase (CRUOX) encodes a functional protein representing a unique case of cysteine enrichment in the evolution of vertebrate orthologous genes. CRUOX retains the ability to catalyze urate oxidation to hydrogen peroxide and 5-hydroxyisourate (HIU), albeit with a 100-fold reduced efficiency. However, differently from all uricases hitherto characterized, it can also facilitate urate regeneration from HIU, a catalytic property which we propose depends on its enrichment in cysteine residues. X-ray structural analysis highlights differences in the active site compared to known orthologs and suggests a mechanism for cysteine-mediated self-aggregation under H2O2-oxidative conditions. Cysteine enrichment was concurrent with transition to uricotelism and a shift in gene expression from the liver to the skin where CRUOX is co-expressed with β-keratins. Therefore, the loss of urate degradation in amniotes has followed opposite evolutionary trajectories: while uricase has been eliminated by pseudogenization in some mammals, it has been repurposed as a redox-sensitive enzyme in the reptilian skin.

Mutation rate is a fundamental parameter in population genetics. Apart from being an important scaling parameter for demographic and phylogenetic inference, it allows one to understand at what rate new genetic diversity is generated and what the expected level of genetic diversity is in a population at equilibrium. However, except for well-established model organisms, accurate estimates of de novo mutation rates are available for a very limited number of organisms from the wild. We estimated mutation rates (µ) in two marine populations of the nine-spined stickleback (Pungitius pungitius) with the aid of several 2- and 3-generational family pedigrees, deep (>50×) whole-genome resequences and a high-quality reference genome. After stringent filtering, we discovered 308 germline mutations in 106 offspring translating to µ = 4.83 × 10-9 and µ = 4.29 × 10-9 per base per generation in the two populations, respectively. Up to 20% of the mutations were shared by full-sibs showing that the level of parental mosaicism was relatively high. Since the estimated µ was 3.1 times smaller than the commonly used substitution rate, recalibration with µ led to substantial increase in estimated divergence times between different stickleback species. Our estimates of the de novo mutation rate should provide a useful resource for research focused on fish population genetics and that of sticklebacks in particular.

Following a duplication, the resulting paralogs tend to diverge. While mutation and natural selection can accelerate this process, they can also slow it. Here, we quantify the paralog homogenization that is caused by point mutations and interlocus gene conversion (IGC). Among 164 duplicated teleost genes, the median percentage of postduplication codon substitutions that arise from IGC rather than point mutation is estimated to be between 7% and 8%. By differentiating between the nonsynonymous codon substitutions that homogenize the protein sequences of paralogs and the nonhomogenizing nonsynonymous substitutions, we estimate the homogenizing nonsynonymous rates to be higher for 163 of the 164 teleost data sets as well as for all 14 data sets of duplicated yeast ribosomal protein-coding genes that we consider. For all 14 yeast data sets, the estimated homogenizing nonsynonymous rates exceed the synonymous rates.

Microbial strategies for resource use are an essential determinant of their fitness in complex habitats. When facing environments with multiple nutrients, microbes often use them sequentially according to a preference hierarchy, resulting in well-known patterns of diauxic growth. In theory, the evolutionary diversification of metabolic hierarchies could represent a mechanism supporting coexistence and biodiversity by enabling temporal segregation of niches. Despite this ecologically critical role, the extent to which substrate preference hierarchies can evolve and diversify remains largely unexplored. Here, we used genome-scale metabolic modeling to systematically explore the evolution of metabolic hierarchies across a vast space of metabolic network genotypes. We find that only a limited number of metabolic hierarchies can readily evolve, corresponding to the most commonly observed hierarchies in genome-derived models. We further show how the evolution of novel hierarchies is constrained by the architecture of central metabolism, which determines both the propensity to change ranks between pairs of substrates and the effect of specific reactions on hierarchy evolution. Our analysis sheds light on the genetic and mechanistic determinants of microbial metabolic hierarchies, opening new research avenues to understand their evolution, evolvability, and ecology.

The human brain is generally anatomically symmetrical, boasting mirror-like brain regions in the left and right hemispheres. Despite this symmetry, fine-scale structural asymmetries are prevalent and are believed to be responsible for distinct functional divisions within the brain. Prior studies propose that these asymmetric structures are predominantly primate specific or even unique to humans, suggesting that the genes contributing to the structural asymmetry of the human brain might have evolved recently. In our study, we identified approximately 1,500 traits associated with human brain asymmetry by collecting paired brain magnetic resonance imaging features from the UK Biobank. Each trait is measured in a specific region of one hemisphere and mirrored in the corresponding region of the other hemisphere. Conducting genome-wide association studies on these traits, we identified over 1,000 quantitative trait loci. Around these index single nucleotide polymorphisms, we found approximately 200 genes that are enriched in brain-related Gene Ontology terms and are predominantly upregulated in brain tissues. Interestingly, most of these genes are evolutionarily old, originating just prior to the emergence of Bilateria (bilaterally symmetrical animals) and Euteleostomi (bony vertebrates with a brain), at a significantly higher ratio than expected. Further analyses of these genes reveal a brain-specific upregulation in humans relative to other mammalian species. This suggests that the structural asymmetry of the human brain has been shaped by evolutionarily ancient genes that have assumed new functions over time.

In this study, we report 21 ancient shotgun genomes from present-day Western Hungary, from previously understudied Late Copper Age Baden, and Bronze Age Somogyvár-Vinkovci, Kisapostag, and Encrusted Pottery archeological cultures (3,530-1,620 cal Bce). Our results indicate the presence of high steppe ancestry in the Somogyvár-Vinkovci culture. They were then replaced by the Kisapostag group, who exhibit an outstandingly high (up to ∼47%) Mesolithic hunter-gatherer ancestry, despite this component being thought to be highly diluted by the time of the Early Bronze Age. The Kisapostag population contributed the genetic basis for the succeeding community of the Encrusted Pottery culture. We also found an elevated hunter-gatherer component in a local Baden culture-associated individual, but no connections were proven to the Bronze Age individuals. The hunter-gatherer ancestry in Kisapostag is likely derived from two main sources, one from a Funnelbeaker or Globular Amphora culture-related population and one from a previously unrecognized source in Eastern Europe. We show that this ancestry not only appeared in various groups in Bronze Age Central Europe but also made contributions to Baltic populations. The social structure of Kisapostag and Encrusted Pottery cultures is patrilocal, similarly to most contemporaneous groups. Furthermore, we developed new methods and method standards for computational analyses of ancient DNA, implemented to our newly developed and freely available bioinformatic package. By analyzing clinical traits, we found carriers of aneuploidy and inheritable genetic diseases. Finally, based on genetic and anthropological data, we present here the first female facial reconstruction from the Bronze Age Carpathian Basin.

Polyploidy is recurrent across the tree of life and known as an evolutionary driving force in plant diversification and crop domestication. How polyploid plants adapt to various habitats has been a fundamental question that remained largely unanswered. Brassica napus is a major crop cultivated worldwide, resulting from allopolyploidy between unknown accessions of diploid B. rapa and B. oleracea. Here, we used whole-genome resequencing data of accessions representing the majority of morphotypes and ecotypes from the species B. rapa, B. oleracea, and B. napus to investigate the role of polyploidy during domestication. To do so, we first reconstructed the phylogenetic history of B. napus, which supported the hypothesis that the emergence of B. napus derived from the hybridization of European turnip of B. rapa and wild B. oleracea. These analyses also showed that morphotypes of swede and Siberian kale (used as vegetable and fodder) were domesticated before rapeseed (oil crop). We next observed that frequent interploidy introgressions from sympatric diploids were prominent throughout the domestication history of B. napus. Introgressed genomic regions were shown to increase the overall genetic diversity and tend to be localized in regions of high recombination. We detected numerous candidate adaptive introgressed regions and found evidence that some of the genes in these regions contributed to phenotypic diversification and adaptation of different morphotypes. Overall, our results shed light on the origin and domestication of B. napus and demonstrate interploidy introgression as an important mechanism that fuels rapid diversification in polyploid species.

Novel phenotypes are increasingly recognized to have evolved by co-option of conserved genes into new developmental contexts, yet the process by which co-opted genes modify existing developmental programs remains obscure. Here, we provide insight into this process by characterizing the role of co-opted doublesex in butterfly wing color pattern development. dsx is the master regulator of insect sex differentiation but has been co-opted to control the switch between discrete nonmimetic and mimetic patterns in Papilio alphenor and its relatives through the evolution of novel mimetic alleles. We found dynamic spatial and temporal expression pattern differences between mimetic and nonmimetic butterflies throughout wing development. A mimetic color pattern program is switched on by a pulse of dsx expression in early pupal development that causes acute and long-term differential gene expression, particularly in Wnt and Hedgehog signaling pathways. RNAi suggested opposing, novel roles for these pathways in mimetic pattern development. Importantly, Dsx co-option caused Engrailed, a primary target of Hedgehog signaling, to gain a novel expression domain early in pupal wing development that is propagated through mid-pupal development to specify novel mimetic patterns despite becoming decoupled from Dsx expression itself. Altogether, our findings provide multiple views into how co-opted genes can both cause and elicit changes to conserved networks and pathways to result in development of novel, adaptive phenotypes.

Rapid evolution of similar phenotypes in similar environments, giving rise to in situ parallel adaptation, is an important hallmark of ecological speciation. However, what appears to be in situ adaptation can also arise by dispersal of divergent lineages from elsewhere. We test whether two contrasting phenotypes repeatedly evolved in parallel, or have a single origin, in an archetypal example of ecological adaptive radiation: benthic-limnetic three-spined stickleback (Gasterosteus aculeatus) across species pair and solitary lakes in British Columbia. We identify two genomic clusters across freshwater populations, which differ in benthic-limnetic divergent phenotypic traits and separate benthic from limnetic individuals in species pair lakes. Phylogenetic reconstruction and niche evolution modeling both suggest a single evolutionary origin for each of these clusters. We detected strong phylogenetic signal in benthic-limnetic divergent traits, suggesting that they are ancestrally retained. Accounting for ancestral state retention, we identify local adaptation of body armor due to the presence of an intraguild predator, the sculpin (Cottus asper), and environmental effects of lake depth and pH on body size. Taken together, our results imply a predominant role for retention of ancestral characteristics in driving trait distribution, with further selection imposed on some traits by environmental factors.

Horizontal gene transfer (HGT) is a means of exchanging genetic material asexually. The process by which horizontally transferred genes are domesticated by the host genome is of great interest but is not well understood. In this study, we determined the telomere-to-telomere genome sequence of the wheat-infecting Pyricularia oryzae strain Br48. SNP analysis indicated that the Br48 strain is a hybrid of wheat- and Brachiaria-infecting strains by a sexual or parasexual cross. Comparative genomic analysis identified several megabase-scale "insertions" in the Br48 genome, some of which were possibly gained by HGT-related events from related species, such as P. pennisetigena or P. grisea. Notably, the mega-insertions often contained genes whose phylogeny is not congruent with the species phylogeny. Moreover, some of the genes have a close homolog even in distantly related organisms, such as basidiomycetes or prokaryotes, implying the involvement of multiple HGT events. Interestingly, the levels of the silent epigenetic marks H3K9me3 and H3K27me3 in a genomic region tended to be negatively correlated with the phylogenetic concordance of genes in the same region, suggesting that horizontally transferred DNA is preferentially targeted for epigenetic silencing. Indeed, the putative HGT-derived genes were activated when MoKmt6, the gene responsible for H3K27me3 modification, was deleted. Notably, these genes also tended to be up-regulated during infection, suggesting that they are now under host control and have contributed to establishing a fungal niche. In conclusion, this study suggests that epigenetic modifications have played an important role in the domestication of HGT-derived genes in the P. oryzae genome.

Nonadaptive hypotheses on the evolution of eukaryotic genome size predict an expansion when the process of purifying selection becomes weak. Accordingly, species with huge genomes, such as lungfish, are expected to show a genome-wide relaxation signature of selection compared with other organisms. However, few studies have empirically tested this prediction using genomic data in a comparative framework. Here, we show that 1) the newly assembled transcriptome of the Australian lungfish, Neoceratodus forsteri, is characterized by an excess of pervasive transcription, or transcriptional leakage, possibly due to suboptimal transcriptional control, and 2) a significant relaxation signature in coding genes in lungfish species compared with other vertebrates. Based on these observations, we propose that the largest known animal genomes evolved in a nearly neutral scenario where genome expansion is less efficiently constrained.

MADS-box transcription factors (TFs), among the first TFs extensively studied, exhibit a wide distribution across eukaryotes and play diverse functional roles. Varying by domain architecture, MADS-box TFs in land plants are categorized into Type I (M-type) and Type II (MIKC-type). Type I and II genes have been considered orthologous to the SRF and MEF2 genes in animals, respectively, presumably originating from a duplication before the divergence of eukaryotes. Here, we exploited the increasing availability of eukaryotic MADS-box sequences and reassessed their evolution. While supporting the ancient duplication giving rise to SRF- and MEF2-types, we found that Type I and II genes originated from the MEF2-type genes through another duplication in the most recent common ancestor (MRCA) of land plants. Protein structures predicted by AlphaFold2 and OmegaFold support our phylogenetic analyses, with plant Type I and II TFs resembling the MEF2-type structure, rather than SRFs. We hypothesize that the ancestral SRF-type TFs were lost in the MRCA of Archaeplastida (the kingdom Plantae sensu lato). The retained MEF2-type TFs acquired a Keratin-like domain and became MIKC-type before the divergence of Streptophyta. Subsequently in the MRCA of land plants, M-type TFs evolved from a duplicated MIKC-type precursor through loss of the Keratin-like domain, leading to the Type I clade. Both Type I and II TFs expanded and functionally differentiated in concert with the increasing complexity of land plant body architecture. The recruitment of these originally stress-responsive TFs into developmental programs, including those underlying reproduction, may have facilitated the adaptation to the terrestrial environment.

Despite the increasing abundance of whole transcriptome data, few methods are available to analyze global gene expression across phylogenies. Here, we present a new software package (Computational Analysis of Gene Expression Evolution [CAGEE]) for inferring patterns of increases and decreases in gene expression across a phylogenetic tree, as well as the rate at which these changes occur. In contrast to previous methods that treat each gene independently, CAGEE can calculate genome-wide rates of gene expression, along with ancestral states for each gene. The statistical approach developed here makes it possible to infer lineage-specific shifts in rates of evolution across the genome, in addition to possible differences in rates among multiple tissues sampled from the same species. We demonstrate the accuracy and robustness of our method on simulated data and apply it to a data set of ovule gene expression collected from multiple self-compatible and self-incompatible species in the genus Solanum to test hypotheses about the evolutionary forces acting during mating system shifts. These comparisons allow us to highlight the power of CAGEE, demonstrating its utility for use in any empirical system and for the analysis of most morphological traits. Our software is available at https://github.com/hahnlab/CAGEE/.

Sequencing of reduced representation libraries enables genotyping of many individuals for population genomic studies. However, high amounts of DNA are required, and the method cannot be applied directly on single cells, preventing its use on most microbes. We developed and implemented the analysis of single amplified genomes followed by restriction-site-associated DNA sequencing to bypass labor-intensive culturing and to avoid culturing bias in population genomic studies of unicellular eukaryotes. This method thus opens the way for addressing important questions about the genetic diversity, gene flow, adaptation, dispersal, and biogeography of hitherto unexplored species.

The germline-restricted chromosome (GRC) is likely present in all songbird species but differs widely in size and gene content. This extra chromosome has been described as either a microchromosome with only limited basic gene content or a macrochromosome with enriched gene functions related to female gonad and embryo development. Here, we assembled, annotated, and characterized the first micro-GRC in the blue tit (Cyanistes caeruleus) using high-fidelity long-read sequencing data. Although some genes on the blue tit GRC show signals of pseudogenization, others potentially have important functions, either currently or in the past. We highlight the GRC gene paralog BMP15, which is among the highest expressed GRC genes both in blue tits and in zebra finches (Taeniopygia guttata) and is known to play a role in oocyte and follicular maturation in other vertebrates. The GRC genes of the blue tit are further enriched for functions related to the synaptonemal complex. We found a similar functional enrichment when analyzing published data on GRC genes from two nightingale species (Luscinia spp.). We hypothesize that these genes play a role in maintaining standard maternal inheritance or in recombining maternal and paternal GRCs during potential episodes of biparental inheritance.

Independently acquired envelope (env) genes from endogenous retroviruses have contributed to the placental trophoblast cell-cell fusion in therian mammals. Egg-laying mammals (monotremes) are an important sister clade for understanding mammalian placental evolution, but the env genes in their genomes have yet to be investigated. Here, env-derived open reading frames (env-ORFs) encoding more than 400 amino acid lengths were searched in the genomes of two monotremes: platypus and echidna. Only two env-ORFs were present in the platypus genome, whereas 121 env-ORFs were found in the echidna genome. The echidna env-ORFs were phylogenetically classified into seven groups named env-Tac1 to -Tac7. Among them, the env-Tac1 group contained only a single gene, and its amino acid sequence showed high similarity to those of the RD114/simian type D retroviruses. Using the pseudotyped virus assay, we demonstrated that the Env-Tac1 protein utilizes echidna sodium-dependent neutral amino acid transporter type 1 and 2 (ASCT1 and ASCT2) as entry receptors. Moreover, the Env-Tac1 protein caused cell-cell fusion in human 293T cells depending on the expression of ASCT1 and ASCT2. These results illustrate that fusogenic env genes are not restricted to placental mammals, providing insights into the evolution of retroviral genes and the placenta.

Landscape, climate, and culture can all structure human populations, but few existing methods are designed to simultaneously disentangle among a large number of variables in explaining genetic patterns. We developed a machine learning method for identifying the variables which best explain migration rates, as measured by the coalescent-based program MAPS that uses shared identical by descent tracts to infer spatial migration across a region of interest. We applied our method to 30 human populations in eastern Africa with high-density single nucleotide polymorphism array data. The remarkable diversity of ethnicities, languages, and environments in this region offers a unique opportunity to explore the variables that shape migration and genetic structure. We explored more than 20 spatial variables relating to landscape, climate, and presence of tsetse flies. The full model explained ∼40% of the variance in migration rate over the past 56 generations. Precipitation, minimum temperature of the coldest month, and elevation were the variables with the highest impact. Among the three groups of tsetse flies, the most impactful was fusca which transmits livestock trypanosomiasis. We also tested for adaptation to high elevation among Ethiopian populations. We did not identify well-known genes related to high elevation, but we did find signatures of positive selection related to metabolism and disease. We conclude that the environment has influenced the migration and adaptation of human populations in eastern Africa; the remaining variance in structure is likely due in part to cultural or other factors not captured in our model.

The data available for reconstructing molecular phylogenies have become wildly disparate. Phylogenomic studies can generate data for thousands of genetic markers for dozens of species, but for hundreds of other taxa, data may be available from only a few genes. Can these two types of data be integrated to combine the advantages of both, addressing the relationships of hundreds of species with thousands of genes? Here, we show that this is possible, using data from frogs. We generated a phylogenomic data set for 138 ingroup species and 3,784 nuclear markers (ultraconserved elements [UCEs]), including new UCE data from 70 species. We also assembled a supermatrix data set, including data from 97% of frog genera (441 total), with 1-307 genes per taxon. We then produced a combined phylogenomic-supermatrix data set (a "gigamatrix") containing 441 ingroup taxa and 4,091 markers but with 86% missing data overall. Likelihood analysis of the gigamatrix yielded a generally well-supported tree among families, largely consistent with trees from the phylogenomic data alone. All terminal taxa were placed in the expected families, even though 42.5% of these taxa each had >99.5% missing data and 70.2% had >90% missing data. Our results show that missing data need not be an impediment to successfully combining very large phylogenomic and supermatrix data sets, and they open the door to new studies that simultaneously maximize sampling of genes and taxa.

The formation of new genes during evolution is an important motor of functional innovation, but the rate at which new genes originate and the likelihood that they persist over longer evolutionary periods are still poorly understood questions. Two important mechanisms by which new genes arise are gene duplication and de novo formation from a previously noncoding sequence. Does the mechanism of formation influence the evolutionary trajectories of the genes? Proteins arisen by gene duplication retain the sequence and structural properties of the parental protein, and thus they may be relatively stable. Instead, de novo originated proteins are often species specific and thought to be more evolutionary labile. Despite these differences, here we show that both types of genes share a number of similarities, including low sequence constraints in their initial evolutionary phases, high turnover rates at the species level, and comparable persistence rates in deeper branchers, in both yeast and flies. In addition, we show that putative de novo proteins have an excess of substitutions between charged amino acids compared with the neutral expectation, which is reflected in the rapid loss of their initial highly basic character. The study supports high evolutionary dynamics of different kinds of new genes at the species level, in sharp contrast with the stability observed at later stages.