African lungfish genome sheds light on the vertebrate water-to-land transition.

Lungfishes are the closest extant relatives of tetrapods and preserve ancestral traits linked with the water-to-land transition. However, their huge genome sizes have hindered understanding of this key transition in evolution. Here, we report a 40-Gb chromosome-level assembly of the African lungfish (Protopterus annectens) genome, which is the largest genome assembly ever reported and has a contig and chromosome N50 of 1.60 Mb and 2.81 Gb, respectively. The large size of the lungfish genome is due mainly to retrotransposons. Genes with ultra-long length show similar expression levels to other genes, indicating that lungfishes have evolved high transcription efficacy to keep gene expression balanced. Together with transcriptome and experimental data, we identified potential genes and regulatory elements related to such terrestrial adaptation traits as pulmonary surfactant, anxiolytic ability, pentadactyl limbs, and pharyngeal remodeling. Our results provide insights and key resources for understanding the evolutionary pathway leading from fishes to humans.

A Global Map of G Protein Signaling Regulation by RGS Proteins.

The control over the extent and timing of G protein signaling is provided by the regulator of G protein signaling (RGS) proteins that deactivate G protein α subunits (Gα). Mammalian genomes encode 20 canonical RGS and 16 Gα genes with key roles in physiology and disease. To understand the principles governing the selectivity of Gα regulation by RGS, we examine the catalytic activity of all canonical human RGS proteins and their selectivity for a complete set of Gα substrates using real-time kinetic measurements in living cells. The data reveal rules governing RGS-Gα recognition, the structural basis of its selectivity, and provide principles for engineering RGS proteins with defined selectivity. The study also explores the evolution of RGS-Gα selectivity through ancestral reconstruction and demonstrates how naturally occurring non-synonymous variants in RGS alter signaling. These results provide a blueprint for decoding signaling selectivity and advance our understanding of molecular recognition principles.

Ancestral ribonucleases back in motion for evolutionary-dynamics guided protein design.

The dynamics behavior of a protein is essential for its functionality. Here, Doucet et al. demonstrate how the evolutionary analysis of conformational pathways within a protein family serves to identify common core scaffolds that accommodate branch-specific functional regions controlled by flexibility switches, offering a model for evolutionary-dynamics based protein design.

Computationally guided AAV engineering for enhanced gene delivery.

Gene delivery vehicles based on adeno-associated viruses (AAVs) are enabling increasing success in human clinical trials, and they offer the promise of treating a broad spectrum of both genetic and non-genetic disorders. However, delivery efficiency and targeting must be improved to enable safe and effective therapies. In recent years, considerable effort has been invested in creating AAV variants with improved delivery, and computational approaches have been increasingly harnessed for AAV engineering. In this review, we discuss how computationally designed AAV libraries are enabling directed evolution. Specifically, we highlight approaches that harness sequences outputted by next-generation sequencing (NGS) coupled with machine learning (ML) to generate new functional AAV capsids and related regulatory elements, pushing the frontier of what vector engineering and gene therapy may achieve.

Metamorphic protein folding as evolutionary adaptation.

Metamorphic proteins switch reversibly between multiple distinct, stable structures, often with different functions. It was previously hypothesized that metamorphic proteins arose as intermediates in the evolution of a new fold - rare and transient exceptions to the 'one sequence, one fold' paradigm. However, as described herein, mounting evidence suggests that metamorphic folding is an adaptive feature, preserved and optimized over evolutionary time as exemplified by the NusG family and the chemokine XCL1. Analysis of extant protein families and resurrected protein ancestors demonstrates that large regions of sequence space are compatible with metamorphic folding. As a category that enhances biological fitness, metamorphic proteins are likely to employ fold switching to perform important biological functions and may be more common than previously thought.

Substitution of a single non-coding nucleotide upstream of TMEM216 causes non-syndromic retinitis pigmentosa and is associated with reduced TMEM216 expression.

Genome analysis of individuals affected by retinitis pigmentosa (RP) identified two rare nucleotide substitutions at the same genomic location on chromosome 11 (g.61392563 [GRCh38]), 69 base pairs upstream of the start codon of the ciliopathy gene TMEM216 (c.-69G>A, c.-69G>T [GenBank: NM_001173991.3]), in individuals of South Asian and African ancestry, respectively. Genotypes included 71 homozygotes and 3 mixed heterozygotes in trans with a predicted loss-of-function allele. Haplotype analysis showed single-nucleotide variants (SNVs) common across families, suggesting ancestral alleles within the two distinct ethnic populations. Clinical phenotype analysis of 62 available individuals from 49 families indicated a similar clinical presentation with night blindness in the first decade and progressive peripheral field loss thereafter. No evident systemic ciliopathy features were noted. Functional characterization of these variants by luciferase reporter gene assay showed reduced promotor activity. Nanopore sequencing confirmed the lower transcription of the TMEM216 c.-69G>T allele in blood-derived RNA from a heterozygous carrier, and reduced expression was further recapitulated by qPCR, using both leukocytes-derived RNA of c.-69G>T homozygotes and total RNA from genome-edited hTERT-RPE1 cells carrying homozygous TMEM216 c.-69G>A. In conclusion, these variants explain a significant proportion of unsolved cases, specifically in individuals of African ancestry, suggesting that reduced TMEM216 expression might lead to abnormal ciliogenesis and photoreceptor degeneration.

Amoebozoan testate amoebae illuminate the diversity of heterotrophs and the complexity of ecosystems throughout geological time.

Heterotrophic protists are vital in Earth's ecosystems, influencing carbon and nutrient cycles and occupying key positions in food webs as microbial predators. Fossils and molecular data suggest the emergence of predatory microeukaryotes and the transition to a eukaryote-rich marine environment by 800 million years ago (Ma). Neoproterozoic vase-shaped microfossils (VSMs) linked to Arcellinida testate amoebae represent the oldest evidence of heterotrophic microeukaryotes. This study explores the phylogenetic relationship and divergence times of modern Arcellinida and related taxa using a relaxed molecular clock approach. We estimate the origin of nodes leading to extant members of the Arcellinida Order to have happened during the latest Mesoproterozoic and Neoproterozoic (1054 to 661 Ma), while the divergence of extant infraorders postdates the Silurian. Our results demonstrate that at least one major heterotrophic eukaryote lineage originated during the Neoproterozoic. A putative radiation of eukaryotic groups (e.g., Arcellinida) during the early-Neoproterozoic sustained by favorable ecological and environmental conditions may have contributed to eukaryotic life endurance during the Cryogenian severe ice ages. Moreover, we infer that Arcellinida most likely already inhabited terrestrial habitats during the Neoproterozoic, coexisting with terrestrial Fungi and green algae, before land plant radiation. The most recent extant Arcellinida groups diverged during the Silurian Period, alongside other taxa within Fungi and flowering plants. These findings shed light on heterotrophic microeukaryotes' evolutionary history and ecological significance in Earth's ecosystems, using testate amoebae as a proxy.

Evolution and diversification of carboxylesterase-like [4+2] cyclases in aspidosperma and iboga alkaloid biosynthesis.

Monoterpene indole alkaloids (MIAs) are a large and diverse class of plant natural products, and their biosynthetic construction has been a subject of intensive study for many years. The enzymatic basis for the production of aspidosperma and iboga alkaloids, which are produced exclusively by members of the Apocynaceae plant family, has recently been discovered. Three carboxylesterase (CXE)-like enzymes from Catharanthus roseus and Tabernanthe iboga catalyze regio- and enantiodivergent [4+2] cycloaddition reactions to generate the aspidosperma (tabersonine synthase, TS) and iboga (coronaridine synthase, CorS; catharanthine synthase, CS) scaffolds from a common biosynthetic intermediate. Here, we use a combined phylogenetic and biochemical approach to investigate the evolution and functional diversification of these cyclase enzymes. Through ancestral sequence reconstruction, we provide evidence for initial evolution of TS from an ancestral CXE followed by emergence of CorS in two separate lineages, leading in turn to CS exclusively in the Catharanthus genus. This progression from aspidosperma to iboga alkaloid biosynthesis is consistent with the chemotaxonomic distribution of these MIAs. We subsequently generate and test a panel of chimeras based on the ancestral cyclases to probe the molecular basis for differential cyclization activity. Finally, we show through partial heterologous reconstitution of tabersonine biosynthesis using non-pathway enzymes how aspidosperma alkaloids could have first appeared as "underground metabolites" via recruitment of promiscuous enzymes from common protein families. Our results provide insight into the evolution of biosynthetic enzymes and how new secondary metabolic pathways can emerge through small but important sequence changes following co-option of preexisting enzymatic functions.

Cip1, a CDK regulator, determines heterothallic mating or homothallic selfing in a protist.

Mating type (sex) plays a crucial role in regulating sexual reproduction in most extant eukaryotes. One of the functions of mating types is ensuring self-incompatibility to some extent, thereby promoting genetic diversity. However, heterothallic mating is not always the best mating strategy. For example, in low-density populations or specific environments, such as parasitic ones, species may need to increase the ratio of potential mating partners. Consequently, many species allow homothallic selfing (i.e., self-fertility or intraclonal mating). Throughout the extensive evolutionary history of species, changes in environmental conditions have influenced mating strategies back and forth. However, the mechanisms through which mating-type recognition regulates sexual reproduction and the dynamics of mating strategy throughout evolution remain poorly understood. In this study, we show that the Cip1 protein is responsible for coupling sexual reproduction initiation to mating-type recognition in the protozoal eukaryote Tetrahymena thermophila. Deletion of the Cip1 protein leads to the loss of the selfing-avoidance function of mating-type recognition, resulting in selfing without mating-type recognition. Further experiments revealed that Cip1 is a regulatory subunit of the Cdk19-Cyc9 complex, which controls the initiation of sexual reproduction. These results reveal a mechanism that regulates the choice between mating and selfing. This mechanism also contributes to the debate about the ancestral state of sexual reproduction.

Tree-based QTL mapping with expected local genetic relatedness matrices.

Understanding the genetic basis of complex phenotypes is a central pursuit of genetics. Genome-wide association studies (GWASs) are a powerful way to find genetic loci associated with phenotypes. GWASs are widely and successfully used, but they face challenges related to the fact that variants are tested for association with a phenotype independently, whereas in reality variants at different sites are correlated because of their shared evolutionary history. One way to model this shared history is through the ancestral recombination graph (ARG), which encodes a series of local coalescent trees. Recent computational and methodological breakthroughs have made it feasible to estimate approximate ARGs from large-scale samples. Here, we explore the potential of an ARG-based approach to quantitative-trait locus (QTL) mapping, echoing existing variance-components approaches. We propose a framework that relies on the conditional expectation of a local genetic relatedness matrix (local eGRM) given the ARG. Simulations show that our method is especially beneficial for finding QTLs in the presence of allelic heterogeneity. By framing QTL mapping in terms of the estimated ARG, we can also facilitate the detection of QTLs in understudied populations. We use local eGRM to analyze two chromosomes containing known body size loci in a sample of Native Hawaiians. Our investigations can provide intuition about the benefits of using estimated ARGs in population- and statistical-genetic methods in general.

Reconstructing complex admixture history using a hierarchical model.

Various methods have been proposed to reconstruct admixture histories by analyzing the length of ancestral chromosomal tracts, such as estimating the admixture time and number of admixture events. However, available methods do not explicitly consider the complex admixture structure, which characterizes the joining and mixing patterns of different ancestral populations during the admixture process, and instead assume a simplified one-by-one sequential admixture model. In this study, we proposed a novel approach that considers the non-sequential admixture structure to reconstruct admixture histories. Specifically, we introduced a hierarchical admixture model that incorporated four ancestral populations and developed a new method, called HierarchyMix, which uses the length of ancestral tracts and the number of ancestry switches along genomes to reconstruct the four-way admixture history. By automatically selecting the optimal admixture model using the Bayesian information criterion principles, HierarchyMix effectively estimates the corresponding admixture parameters. Simulation studies confirmed the effectiveness and robustness of HierarchyMix. We also applied HierarchyMix to Uyghurs and Kazakhs, enabling us to reconstruct the admixture histories of Central Asians. Our results highlight the importance of considering complex admixture structures and demonstrate that HierarchyMix is a useful tool for analyzing complex admixture events.

A tale of two copies: Evolutionary trajectories of moth pheromone receptors.

Pheromone communication is an essential component of reproductive isolation in animals. As such, evolution of pheromone signaling can be linked to speciation. For example, the evolution of sex pheromones is thought to have played a major role in the diversification of moths. In the crop pests Spodoptera littoralis and S. litura, the major component of the sex pheromone blend is (Z,E)-9,11-tetradecadienyl acetate, which is lacking in other Spodoptera species. It indicates that a major shift occurred in their common ancestor. It has been shown recently in S. littoralis that this compound is detected with high specificity by an atypical pheromone receptor, named SlitOR5. Here, we studied its evolutionary history through functional characterization of receptors from different Spodoptera species. SlitOR5 orthologs in S. exigua and S. frugiperda exhibited a broad tuning to several pheromone compounds. We evidenced a duplication of OR5 in a common ancestor of S. littoralis and S. litura and found that in these two species, one duplicate is also broadly tuned while the other is specific to (Z,E)-9,11-tetradecadienyl acetate. By using ancestral gene resurrection, we confirmed that this narrow tuning evolved only in one of the two copies issued from the OR5 duplication. Finally, we identified eight amino acid positions in the binding pocket of these receptors whose evolution has been responsible for narrowing the response spectrum to a single ligand. The evolution of OR5 is a clear case of subfunctionalization that could have had a determinant impact in the speciation process in Spodoptera species.

Role of Pfs47 in the dispersal of ancestral Plasmodium falciparum malaria through adaptation to different anopheline vectors.

Plasmodium falciparum malaria originated when Plasmodium praefalciparum, a gorilla malaria parasite transmitted by African sylvan anopheline mosquitoes, adapted to humans. Pfs47, a protein on the parasite surface mediates P. falciparum evasion of the mosquito immune system by interacting with a midgut receptor and is critical for Plasmodium adaptation to different anopheline species. Genetic analysis of 4,971 Pfs47 gene sequences from different continents revealed that Asia and Papua New Guinea harbor Pfs47 haplotypes more similar to its ortholog in P. praefalciparum at sites that determine vector compatibility, suggesting that ancestral P. falciparum readily adapted to Asian vectors. Consistent with this observation, Pfs47-receptor gene sequences from African sylvan malaria vectors, such as Anopheles moucheti and An. marshallii, were found to share greater similarity with those of Asian vectors than those of vectors of the African An. gambiae complex. Furthermore, experimental infections provide direct evidence that transformed P. falciparum parasites carrying Pfs47 orthologs of P. praefalciparum or P. reichenowi were more effective at evading the immune system of the Asian malaria vector An. dirus than An. gambiae. We propose that high compatibility of ancestral P. falciparum Pfs47 with the receptors of Asian vectors facilitated the early dispersal of human malaria to the Asian continent, without having to first adapt to sub-Saharan vectors of the An. gambiae complex.

Deciphering the evolution of flavin-dependent monooxygenase stereoselectivity using ancestral sequence reconstruction.

Controlling the selectivity of a reaction is critical for target-oriented synthesis. Accessing complementary selectivity profiles enables divergent synthetic strategies, but is challenging to achieve in biocatalytic reactions given enzymes' innate preferences of a single selectivity. Thus, it is critical to understand the structural features that control selectivity in biocatalytic reactions to achieve tunable selectivity. Here, we investigate the structural features that control the stereoselectivity in an oxidative dearomatization reaction that is key to making azaphilone natural products. Crystal structures of enantiocomplementary biocatalysts guided the development of multiple hypotheses centered on the structural features that control the stereochemical outcome of the reaction; however, in many cases, direct substitutions of active site residues in natural proteins led to inactive enzymes. Ancestral sequence reconstruction (ASR) and resurrection were employed as an alternative strategy to probe the impact of each residue on the stereochemical outcome of the dearomatization reaction. These studies suggest that two mechanisms are active in controlling the stereochemical outcome of the oxidative dearomatization reaction: one involving multiple active site residues in AzaH and the other dominated by a single Phe to Tyr switch in TropB and AfoD. Moreover, this study suggests that the flavin-dependent monooxygenases (FDMOs) adopt simple and flexible strategies to control stereoselectivity, which has led to stereocomplementary azaphilone natural products produced by fungi. This paradigm of combining ASR and resurrection with mutational and computational studies showcases sets of tools for understanding enzyme mechanisms and provides a solid foundation for future protein engineering efforts.

Ancestral sequence reconstruction dissects structural and functional differences among eosinophil ribonucleases.

Evolutionarily conserved structural folds can give rise to diverse biological functions, yet predicting atomic-scale interactions that contribute to the emergence of novel activities within such folds remains challenging. Pancreatic-type ribonucleases illustrate this complexity, sharing a core structure that has evolved to accommodate varied functions. In this study, we used ancestral sequence reconstruction to probe evolutionary and molecular determinants that distinguish biological activities within eosinophil members of the RNase 2/3 subfamily. Our investigation unveils functional, structural, and dynamical behaviors that differentiate the evolved ancestral ribonuclease (AncRNase) from its contemporary eosinophil RNase orthologs. Leveraging the potential of ancestral reconstruction for protein engineering, we used AncRNase predictions to design a minimal 4-residue variant that transforms human RNase 2 into a chimeric enzyme endowed with the antimicrobial and cytotoxic activities of RNase 3 members. This work provides unique insights into mutational and evolutionary pathways governing structure, function, and conformational states within the eosinophil RNase subfamily, offering potential for targeted modulation of RNase-associated functions.

The haplotype-resolved genome of diploid Chrysanthemum indicum unveils new acacetin synthases genes and their evolutionary history.

Acacetin, a flavonoid compound, possesses a wide range of pharmacological effects, including antimicrobial, immune regulation, and anticancer effects. Some key steps in its biosynthetic pathway were largely unknown in flowering plants. Here, we present the first haplotype-resolved genome of Chrysanthemum indicum, whose dried flowers contain abundant flavonoids and have been utilized as traditional Chinese medicine. Various phylogenetic analyses revealed almost equal proportion of three tree topologies among three Chrysanthemum species (C. indicum, C. nankingense, and C. lavandulifolium), indicating that frequent gene flow among Chrysanthemum species or incomplete lineage sorting due to rapid speciation might contribute to conflict topologies. The expanded gene families in C. indicum were associated with oxidative functions. Through comprehensive candidate gene screening, we identified five flavonoid O-methyltransferase (FOMT) candidates, which were highly expressed in flowers and whose expressional levels were significantly correlated with the content of acacetin. Further experiments validated two FOMTs (CI02A009970 and CI03A006662) were capable of catalyzing the conversion of apigenin into acacetin, and these two genes are possibly responsible acacetin accumulation in disc florets and young leaves, respectively. Furthermore, combined analyses of ancestral chromosome reconstruction and phylogenetic trees revealed the distinct evolutionary fates of the two validated FOMT genes. Our study provides new insights into the biosynthetic pathway of flavonoid compounds in the Asteraceae family and offers a model for tracing the origin and evolutionary routes of single genes. These findings will facilitate in vitro biosynthetic production of flavonoid compounds through cellular and metabolic engineering and expedite molecular breeding of C. indicum cultivars.

Complex Genomes of Early Nucleocytoviruses Revealed by Ancient Origins of Viral Aminoacyl-tRNA Synthetases.

Aminoacyl-tRNA synthetases (aaRSs), also known as tRNA ligases, are essential enzymes in translation. Owing to their functional essentiality, these enzymes are conserved in all domains of life and used as informative markers to trace the evolutionary history of cellular organisms. Unlike cellular organisms, viruses generally lack aaRSs because of their obligate parasitic nature, but several large and giant DNA viruses in the phylum Nucleocytoviricota encode aaRSs in their genomes. The discovery of viral aaRSs led to the idea that the phylogenetic analysis of aaRSs can shed light on ancient viral evolution. However, conflicting results have been reported from previous phylogenetic studies: one posited that nucleocytoviruses recently acquired their aaRSs from their host eukaryotes, while another hypothesized that the viral aaRSs have ancient origins. Here, we investigated 4,168 nucleocytovirus genomes, including metagenome-assembled genomes (MAGs) derived from large-scale metagenomic studies. In total, we identified 780 viral aaRS sequences in 273 viral genomes. We generated and examined phylogenetic trees of these aaRSs with a large set of cellular sequences to trace evolutionary relationships between viral and cellular aaRSs. The analyses suggest that the origins of some viral aaRSs predate the last common eukaryotic ancestor. Inside viral aaRS clades, we identify intricate evolutionary trajectories of viral aaRSs with horizontal transfers, losses, and displacements. Overall, these results suggest that ancestral nucleocytoviruses already developed complex genomes with an expanded set of aaRSs in the proto-eukaryotic era.

Biases in ARG-Based Inference of Historical Population Size in Populations Experiencing Selection.

Inferring the demographic history of populations provides fundamental insights into species dynamics and is essential for developing a null model to accurately study selective processes. However, background selection and selective sweeps can produce genomic signatures at linked sites that mimic or mask signals associated with historical population size change. While the theoretical biases introduced by the linked effects of selection have been well established, it is unclear whether ancestral recombination graph (ARG)-based approaches to demographic inference in typical empirical analyses are susceptible to misinference due to these effects. To address this, we developed highly realistic forward simulations of human and Drosophila melanogaster populations, including empirically estimated variability of gene density, mutation rates, recombination rates, purifying, and positive selection, across different historical demographic scenarios, to broadly assess the impact of selection on demographic inference using a genealogy-based approach. Our results indicate that the linked effects of selection minimally impact demographic inference for human populations, although it could cause misinference in populations with similar genome architecture and population parameters experiencing more frequent recurrent sweeps. We found that accurate demographic inference of D. melanogaster populations by ARG-based methods is compromised by the presence of pervasive background selection alone, leading to spurious inferences of recent population expansion, which may be further worsened by recurrent sweeps, depending on the proportion and strength of beneficial mutations. Caution and additional testing with species-specific simulations are needed when inferring population history with non-human populations using ARG-based approaches to avoid misinference due to the linked effects of selection.

Please Mind the Gap: Indel-Aware Parsimony for Fast and Accurate Ancestral Sequence Reconstruction and Multiple Sequence Alignment Including Long Indels.

Despite having important biological implications, insertion, and deletion (indel) events are often disregarded or mishandled during phylogenetic inference. In multiple sequence alignment, indels are represented as gaps and are estimated without considering the distinct evolutionary history of insertions and deletions. Consequently, indels are usually excluded from subsequent inference steps, such as ancestral sequence reconstruction and phylogenetic tree search. Here, we introduce indel-aware parsimony (indelMaP), a novel way to treat gaps under the parsimony criterion by considering insertions and deletions as separate evolutionary events and accounting for long indels. By identifying the precise location of an evolutionary event on the tree, we can separate overlapping indel events and use affine gap penalties for long indel modeling. Our indel-aware approach harnesses the phylogenetic signal from indels, including them into all inference stages. Validation and comparison to state-of-the-art inference tools on simulated data show that indelMaP is most suitable for densely sampled datasets with closely to moderately related sequences, where it can reach alignment quality comparable to probabilistic methods and accurately infer ancestral sequences, including indel patterns. Due to its remarkable speed, our method is well suited for epidemiological datasets, eliminating the need for downsampling and enabling the exploitation of the additional information provided by dense taxonomic sampling. Moreover, indelMaP offers new insights into the indel patterns of biologically significant sequences and advances our understanding of genetic variability by considering gaps as crucial evolutionary signals rather than mere artefacts.

Deciphering Structural Traits for Thermal and Kinetic Stability across Protein Family Evolution through Ancestral Sequence Reconstruction.

Natural proteins are frequently marginally stable, and an increase in environmental temperature can easily lead to unfolding. As a result, protein engineering to improve protein stability is an area of intensive research. Nonetheless, since there is usually a high degree of structural homology between proteins from thermophilic organisms and their mesophilic counterparts, the identification of structural determinants for thermoadaptation is challenging. Moreover, in many cases, it has become clear that the success of stabilization strategies is often dependent on the evolutionary history of a protein family. In the last few years, the use of ancestral sequence reconstruction (ASR) as a tool for elucidation of the evolutionary history of functional traits of a protein family has gained strength. Here, we used ASR to trace the evolutionary pathways between mesophilic and thermophilic kinases that participate in the biosynthetic pathway of vitamin B1 in bacteria. By combining biophysics approaches, X-ray crystallography, and molecular dynamics simulations, we found that the thermal stability of these enzymes correlates with their kinetic stability, where the highest thermal/kinetic stability is given by an increase in small hydrophobic amino acids that allow a higher number of interatomic hydrophobic contacts, making this type of interaction the main support for stability in this protein architecture. The results highlight the potential benefits of using ASR to explore the evolutionary history of protein sequence and structure to identify traits responsible for the kinetic and thermal stability of any protein architecture.