Convergent Amino Acid Research Articles

Comparative genomics sheds new light on the convergent evolution of infrared vision in snakes.

Infrared vision is a highly specialized sensory system that evolved independently in three clades of snakes. Apparently, convergent evolution occurred in the transient receptor potential ankyrin 1 (TRPA1) proteins of infrared-sensing snakes. However, this gene can only explain how infrared signals are received, and not the transduction and processing of those signals. We sequenced the genome of Xenopeltis unicolor, a key outgroup species of pythons, and performed a genome-wide analysis of convergence between two clades of infrared-sensing snakes. Our results revealed pervasive molecular adaptation in pathways associated with neural development and other functions, with parallel selection on loci associated with trigeminal nerve structural organization. In addition, we found evidence of convergent amino acid substitutions in a set of genes, including TRPA1 and TRPM2. The analysis also identified convergent accelerated evolution in non-coding elements near 12 genes involved in facial nerve structural organization and optic nerve development. Thus, convergent evolution occurred across multiple dimensions of infrared vision in vipers and pythons, as well as amino acid substitutions, non-coding elements, genes and functions. These changes enabled independent groups of snakes to develop and use infrared vision.

Comparative genomics reveals convergent signals associated with the high metabolism and longevity in birds and bats.

Birds and bats have long lifespans relative to their body size compared with non-flying animals. However, the genomic basis associated with longer lifespan of flying species despite their higher metabolism was unclear. In this study, we hypothesized that genes involved in the regulation of metabolism and lifespan changed with the acquisition of flight and searched for genes that show specific evolutionary patterns in flying species. As a result, we identified several genes that show different evolutionary rates in bird and bat lineages. Genes in pathways involved in lifespan regulation were conserved in birds, while they evolved at an accelerated rate in bats. We also searched for genes in which convergent amino acid substitutions occurred in birds and bats and found such substitutions in genes involved in cancer, reactive oxygen species control and immunity. Our study revealed genomic changes associated with the acquisition of flight in birds and bats and suggested that multiple genes involved in the regulation of lifespan and metabolism support both high metabolism and longevity in flying species.

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage‐specific molecular evolution in deep‐sea adaptation

AbstractAbyssal (3501–6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome‐level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal‐living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep‐sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep‐sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin‐containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep‐sea species shed new light on the common genetic innovations required for deep‐sea adaptation by the fishes.

Genomic signatures associated with recurrent scale loss in cyprinid fish.

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Convergent evolution in high-altitude and marine mammals: Molecular adaptations to pulmonary fibrosis and hypoxia.

High-altitude and marine mammals inhabit distinct ecosystems but share a common challenge: hypoxia. To survive in low-oxygen environments, these species have evolved similar phenotypic pulmonary adaptations, characterized by a high density of elastic fibers. In this study, we explored the molecular mechanisms underlying these adaptations, focusing on pulmonary fibrosis and hypoxia tolerance through comparative genomics and convergent evolution analyses. We observed significant expansions and contractions in certain gene families across both high-altitude and marine mammals, closely associated with processes involved in pulmonary fibrosis. Notably, members of the keratin gene family, such as KRT17 and KRT14, appear to be associated with the development of the dense elastic fiber phenotype observed in the lungs of hypoxia-tolerant mammals. Through selection pressure and amino acid substitution analyses, we identified multiple genes exhibiting convergent accelerated evolution, positive selection, and amino acid substitution in these species, associated with adaptation to hypoxic environments. Specifically, the convergent evolution of ZFP36L1, FN1, and NEDD9 was found to contribute to the high density of elastic fibers in the lungs of both high-altitude and marine mammals, facilitating their hypoxia tolerance. Additionally, we identified convergent amino acid substitutions and gene loss events associated with sperm development, differentiation, and spermatogenesis, such as amino acid substitutions in SLC26A3 and pseudogenization of CFAP47, as confirmed by PCR. These genetic alterations may be linked to changes in the reproductive capabilities of these animals. Overall, this study offers novel perspectives on the genetic and molecular adaptations of high-altitude and marine mammals to hypoxic environments, with a particular emphasis on pulmonary fibrosis.

Duplications of human longevity-associated genes across placental mammals.

Natural selection has shaped a wide range of lifespans across mammals, with a few long-lived species showing negligible signs of ageing. Approaches used to elucidate the genetic mechanisms underlying mammalian longevity usually involve phylogenetic selection tests on candidate genes, detections of convergent amino acid changes in long-lived lineages, analyses of differential gene expression between age cohorts or species, and measurements of age-related epigenetic changes. However, the link between gene duplication and evolution of mammalian longevity has not been widely investigated. Here, we explored the association between gene duplication and mammalian lifespan by analysing 287 human longevity-associated genes across 37 placental mammals. We estimated that the expansion rate of these genes is eight times higher than their contraction rate across these 37 species. Using phylogenetic approaches, we identified 43 genes whose duplication levels are significantly correlated with longevity quotients (FDR < 0.05). In particular, the strong correlation observed for four genes (CREBBP, PIK3R1, HELLS, FOXM1) appears to be driven mainly by their high duplication levels in two ageing extremists, the naked mole rat (Heterocephalus glaber) and the greater mouse-eared bat (Myotis myotis). Further sequence and expression analyses suggest that the gene PIK3R1 may have undergone a convergent duplication event, whereby the similar region of its coding sequence was independently duplicated multiple times in both of these long-lived species. Collectively, this study identified several candidate genes whose duplications may underlie the extreme longevity in mammals, and highlighted the potential role of gene duplication in the evolution of mammalian long lifespans.

CAAStools: a toolbox to identify and test Convergent Amino Acid Substitutions.

Coincidence of Convergent Amino Acid Substitutions (CAAS) with phenotypic convergences allow pinpointing genes and even individual mutations that are likely to be associated with trait variation within their phylogenetic context. Such findings can provide useful insights into the genetic architecture of complex phenotypes. Here we introduce CAAStools, a set of bioinformatics tools to identify and validate CAAS in orthologous protein alignments for pre-defined groups of species representing the phenotypic values targeted by the user. CAAStools source code is available at http://github.com/linudz/caastools, along with documentation and examples. Supplementary data are available at Bioinformatics online.

Signatures of Extreme Longevity: A Perspective from Bivalve Molecular Evolution.

Among Metazoa, bivalves have the highest lifespan disparity, ranging from 1 to 500+ years, making them an exceptional testing ground to understand mechanisms underlying aging and the evolution of extended longevity. Nevertheless, comparative molecular evolution has been an overlooked approach in this instance. Here, we leveraged transcriptomic resources spanning 30 bivalve species to unravel the signatures of convergent molecular evolution in four long-lived species: Margaritifera margaritifera, Elliptio complanata, Lampsilis siliquoidea, and Arctica islandica (the latter represents the longest-lived noncolonial metazoan known so far). We applied a comprehensive approach-which included inference of convergent dN/dS, convergent positive selection, and convergent amino acid substitution-with a strong focus on the reduction of false positives. Genes with convergent evolution in long-lived bivalves show more physical and functional interactions to each other than expected, suggesting that they are biologically connected; this interaction network is enriched in genes for which a role in longevity has been experimentally supported in other species. This suggests that genes in the network are involved in extended longevity in bivalves and, consequently, that the mechanisms underlying extended longevity are-at least partially-shared across Metazoa. Although we believe that an integration of different genes and pathways is required for the extended longevity phenotype, we highlight the potential central roles of genes involved in cell proliferation control, translational machinery, and response to hypoxia, in lifespan extension.

Convergent resistance to GABA receptor neurotoxins through plant-insect coevolution.

The molecular mechanisms of coevolution between plants and insects remain elusive. GABA receptors are targets of many neurotoxic terpenoids, which represent the most diverse array of natural products known. Over deep evolutionary time, as plant terpene synthases diversified in plants, so did plant terpenoid defence repertoires. Here we show that herbivorous insects and their predators evolved convergent amino acid changing substitutions in duplicated copies of the Resistance to dieldrin (Rdl) gene that encodes the GABA receptor, and that the evolution of duplicated Rdl and terpenoid-resistant GABA receptors is associated with the diversification of moths and butterflies. These same substitutions also evolved in pests exposed to synthetic insecticides that target the GABA receptor. We used in vivo genome editing in Drosophila melanogaster to evaluate the fitness effects of each putative resistance mutation and found that pleiotropy both facilitates and constrains the evolution of GABA receptor resistance. The same genetic changes that confer resistance to terpenoids across 300 Myr of insect evolution have re-evolved in response to synthetic analogues over one human lifespan.

Convergent molecular evolution of thermogenesis and circadian rhythm in Arctic ruminants.

The muskox and reindeer are the only ruminants that have evolved to survive in harsh Arctic environments. However, the genetic basis of this Arctic adaptation remains largely unclear. Here, we compared a de novo assembled muskox genome with reindeer and other ruminant genomes to identify convergent amino acid substitutions, rapidly evolving genes and positively selected genes among the two Arctic ruminants. We found these candidate genes were mainly involved in brown adipose tissue (BAT) thermogenesis and circadian rhythm. Furthermore, by integrating transcriptomic data from goat adipose tissues (white and brown), we demonstrated that muskox and reindeer may have evolved modulating mitochondrion, lipid metabolism and angiogenesis pathways to enhance BAT thermogenesis. In addition, results from co-immunoprecipitation experiments prove that convergent amino acid substitution of the angiogenesis-related gene hypoxia-inducible factor 2alpha (HIF2A), resulting in weakening of its interaction with prolyl hydroxylase domain-containing protein 2 (PHD2), may increase angiogenesis of BAT. Altogether, our work provides new insights into the molecular mechanisms involved in Arctic adaptation.

The evolution of the discrete multirenculate kidney in mammals from ecological and molecular perspectives.

Mammals have developed different kinds of renal structures during evolution, yet the origin of the renal structural phenotypes and the molecular mechanisms underlying their adaptive evolution remains unclear. Here, we reconstructed the ancestral state of the renal structures across mammals and found that the unilobar kidney was the ancestral character in mammals. The subsequent correlation analyses between renal phenotypes and life history traits revealed that species with a larger body or in aquatic habitats tend to have evolved discrete multirenculate kidneys. To explore the molecular convergent mechanisms among mammals with this most distinct renal structure, the discrete multirenculate kidney, we used 45 genes related to duplex/multiplex kidney diseases to compare the evolutions of species with discrete multirenculate kidneys and with other renal phenotypes. Twelve rapidly evolving genes that were functionally enriched in cilium assembly and centrosome were identified in species with discrete multirenculate kidneys, suggesting that these genes played key roles in the evolution of discrete multirenculate kidneys. In addition, positive selection was detected in six crucial genes which are mainly involved in epithelial tube morphogenesis and regulation of neurogenesis. Finally, 12 convergent amino acid substitutions, six of which are in crucial domain of proteins, were shared by two or more lineages with discrete multirenculate kidneys. These findings could provide some novel insights into the origin and evolution of renal structures across mammals and the pathogenesis of renal diseases in humans.

The adaptive evolution of Leuciscus waleckii in Lake Dali Nur and convergent evolution of Cypriniformes fishes inhabiting extremely alkaline environments.

Leuciscus waleckii is widely distributed in Northeast Asia and has high economic value. The population in Lake Dali Nur can adapt to extremely alkaline-saline water with bicarbonate over 50 mmol/L (pH 9.6), thus providing an exceptional model for exploring the mechanisms of adaptive evolution under extreme alkaline environments. Here, we assembled a high-quality chromosome-level reference genome for L. waleckii from Lake Dali Nur. Based on the resequencing of 85 individuals from divergent populations, the historical population size of L.waleckii in Lake Dali Nur dramatically expanded in a thousand years approximately 13,000 years ago, and experienced a cliff recession in the process of adapting to the alkaline environment of Lake Dali Nur approximately 6,000 years ago. Genome scans between freshwater and alkaline populations further revealed the significant selective sweep regions from Lake Dali Nur, which harbour a set of candidate genes involved in hypoxia tolerance, ion transport, acid-base regulation and nitrogen metabolism. 5 alkali population-specific nonsynonymous mutations were identified in CA15 gene copies. In addition, two sites with convergent amino acid mutation were detected in the RHCG-a gene among several alkali environment adapted Cypriniformes fish. Our findings provide comprehensive insight into the genomic mechanisms of L. waleckii and reveal their adaptative evolution under extreme alkaline environments.

Convergent Genomic Signatures of Cashmere Traits: Evidence for Natural and Artificial Selection.

Convergent evolution provides powerful opportunities to investigate the genetic basis of complex traits. The Tibetan antelope (Pantholops hodgsonii) and Siberian ibex (Capra sibirica) belong to different subfamilies in Bovidae, but both have evolved similar superfine cashmere characteristics to meet the cold temperature in plateau environments. The cashmere traits of cashmere goats underwent strong artificial selection, and some traces of domestication also remained in the genome. Hence, we investigated the convergent genomic signatures of cashmere traits between natural and artificial selection. We compared the patterns of convergent molecular evolution between Tibetan antelope and Siberian ibex by testing positively selected genes, rapidly evolving genes and convergent amino acid substitutions. In addition, we analyzed the selected genomic features of cashmere goats under artificial selection using whole-genome resequencing data, and skin transcriptome data of cashmere goats were also used to focus on the genes involved in regulating cashmere traits. We found that molecular convergent events were very rare, but natural and artificial selection genes were convergent enriched in similar functional pathways (e.g., ECM-receptor interaction, focal adhesion, PI3K-Akt signaling pathway) in a variety of gene sets. Type IV collagen family genes (COL4A2, COL4A4, COL4A5, COL6A5, COL6A6) and integrin family genes (ITGA2, ITGA4, ITGA9, ITGB8) may be important candidate genes for cashmere formation and development. Our results provide a comprehensive approach and perspective for exploring cashmere traits and offer a valuable reference for subsequent in-depth research on the molecular mechanisms regulating cashmere development and fineness.

Rubbing Salt in the Wound: Molecular Evolutionary Analysis of Pain-Related Genes Reveals the Pain Adaptation of Cetaceans in Seawater.

Pain, usually caused by a strong or disruptive stimulus, is an unpleasant sensation that serves as a warning to organisms. To adapt to extreme environments, some terrestrial animals have evolved to be inherently insensitive to pain. Cetaceans are known as supposedly indifferent to pain from soft tissue injury representatives of marine mammals. However, the molecular mechanisms that explain how cetaceans are adapted to pain in response to seawater environment remain unclear. Here, we performed a molecular evolutionary analysis of pain-related genes in selected representatives of cetaceans. ASIC4 gene was identified to be pseudogenized in all odontocetes (toothed whales) except from Physeter macrocephalus (sperm whales), and relaxed selection of this gene was detected in toothed whales with pseudogenized ASIC4. In addition, positive selection was detected in pain perception (i.e., ASIC3, ANO1, CCK, and SCN9A) and analgesia (i.e., ASIC3, ANO1, CCK, and SCN9A) genes among the examined cetaceans. In this study, potential convergent amino acid substitutions within predicted proteins were found among the examined cetaceans and other terrestrial mammals, inhabiting extreme environments (e.g., V441I of TRPV1 in cetaceans and naked mole rats). Moreover, specific amino acid substitutions within predicted sequences of several proteins were found in the studied representatives of cetaceans (e.g., F56L and D163A of ASIC3, E88G of GRK2, and F159L of OPRD1). Most of the substitutions were located within important functional domains of proteins, affecting their protein functions. The above evidence suggests that cetaceans might have undergone adaptive molecular evolution in pain-related genes through different evolutionary patterns to adapt to pain, resulting in greater sensitivity to pain and more effective analgesia. This study could have implications for diagnosis and treatment of human pain.

Epistatic Effects Between Amino Acid Insertions and Substitutions Mediate Toxin resistance of Vertebrate Na+,K+-ATPases.

The recurrent evolution of resistance to cardiotonic steroids (CTS) across diverse animals most frequently involves convergent amino acid substitutions in the H1-H2 extracellular loop of Na+,K+-ATPase (NKA). Previous work revealed that hystricognath rodents (e.g., chinchilla) and pterocliform birds (sandgrouse) have convergently evolved amino acid insertions in the H1-H2 loop, but their functional significance was not known. Using protein engineering, we show that these insertions have distinct effects on CTS resistance in homologs of each of the two species that strongly depend on intramolecular interactions with other residues. Removing the insertion in the chinchilla NKA unexpectedly increases CTS resistance and decreases NKA activity. In the sandgrouse NKA, the amino acid insertion and substitution Q111R both contribute to an augmented CTS resistance without compromising ATPase activity levels. Molecular docking simulations provide additional insight into the biophysical mechanisms responsible for the context-specific mutational effects on CTS insensitivity of the enzyme. Our results highlight the diversity of genetic substrates that underlie CTS insensitivity in vertebrate NKA and reveal how amino acid insertions can alter the phenotypic effects of point mutations at key sites in the same protein domain.

A high-quality genome of the dobsonfly Neoneuromus ignobilis reveals molecular convergences in aquatic insects

Neoneuromus ignobilis is an archaic holometabolous aquatic predatory insect. However, a lack of genomic resources hinders the use of whole genome sequencing to explore their genetic basis and molecular mechanisms for adaptive evolution. Here, we provided a high-contiguity, chromosome-level genome assembly of N. ignobilis using high coverage Nanopore and PacBio reads with the Hi-C technique. The final assembly is 480.67 MB in size, containing 12 telomere-ended pseudochromosomes with only 17 gaps. We compared 42 hexapod species genomes including six independent lineages comprising 11 aquatic insects, and found convergent expansions of long wavelength-sensitive and blue-sensitive opsins, thermal stress response TRP channels, and sulfotransferases in aquatic insects, which may be related to their aquatic adaptation. We also detected strong nonrandom signals of convergent amino acid substitutions in aquatic insects. Collectively, our comparative genomic analysis revealed the evidence of molecular convergences in aquatic insects during both gene family evolution and convergent amino acid substitutions.

Comparative Genomics Provides Insights into Adaptive Evolution in Tactile-Foraging Birds.

Tactile-foraging birds have evolved an enlarged principal sensory nucleus (PrV) but smaller brain regions related to the visual system, which reflects the difference in sensory dependence. The “trade-off” may exist between different senses in tactile foragers, as well as between corresponding sensory-processing areas in the brain. We explored the mechanism underlying the adaptive evolution of sensory systems in three tactile foragers (kiwi, mallard, and crested ibis). The results showed that olfaction-related genes in kiwi and mallard and hearing-related genes in crested ibis were expanded, indicating they may also have sensitive olfaction or hearing, respectively. However, some genes required for visual development were positively selected or had convergent amino acid substitutions in all three tactile branches, and it seems to show the possibility of visual degradation. In addition, we may provide a new visual-degradation candidate gene PDLIM1 who suffered dense convergent amino acid substitutions within the ZM domain. At last, two genes responsible for regulating the proliferation and differentiation of neuronal progenitor cells may play roles in determining the relative sizes of sensory areas in brain. This exploration offers insight into the relationship between specialized tactile-forging behavior and the evolution of sensory abilities and brain structures.

E-Volve: understanding the impact of mutations in SARS-CoV-2 variants spike protein on antibodies and ACE2 affinity through patterns of chemical interactions at protein interfaces.

BackgroundThe SARS-CoV-2 pandemic reverberated, posing health and social hygiene obstacles throughout the globe. Mutant lineages of the virus have concerned scientists because of convergent amino acid alterations, mainly on the viral spike protein. Studies have shown that mutants have diminished activity of neutralizing antibodies and enhanced affinity with its human cell receptor, the ACE2 protein.MethodsHence, for real-time measuring of the impacts caused by variant strains in such complexes, we implemented E-Volve, a tool designed to model a structure with a list of mutations requested by users and return analyses of the variant protein. As a proof of concept, we scrutinized the spike-antibody and spike-ACE2 complexes formed in the variants of concern, B.1.1.7 (Alpha), B.1.351 (Beta), and P.1 (Gamma), by using contact maps depicting the interactions made amid them, along with heat maps to quantify these major interactions.ResultsThe results found in this study depict the highly frequent interface changes made by the entire set of mutations, mainly conducted by N501Y and E484K. In the spike-Antibody complex, we have noticed alterations concerning electrostatic surface complementarity, breaching essential sites in the P17 and BD-368-2 antibodies. Alongside, the spike-ACE2 complex has presented new hydrophobic bonds.DiscussionMolecular dynamics simulations followed by Poisson-Boltzmann calculations corroborate the higher complementarity to the receptor and lower to the antibodies for the K417T/E484K/N501Y (Gamma) mutant compared to the wild-type strain, as pointed by E-Volve, as well as an intensification of this effect by changes at the protein conformational equilibrium in solution. A local disorder of the loop α1′/β1′, as well its possible effects on the affinity to the BD-368-2 antibody were also incorporated to the final conclusions after this analysis. Moreover, E-Volve can depict the main alterations in important biological structures, as shown in the SARS-CoV-2 complexes, marking a major step in the real-time tracking of the virus mutant lineages. E-Volve is available at http://bioinfo.dcc.ufmg.br/evolve.

Beyond RuBisCO: convergent molecular evolution of multiple chloroplast genes in C4 plants.

BackgroundThe recurrent evolution of the C4 photosynthetic pathway in angiosperms represents one of the most extraordinary examples of convergent evolution of a complex trait. Comparative genomic analyses have unveiled some of the molecular changes associated with the C4 pathway. For instance, several key enzymes involved in the transition from C3 to C4 photosynthesis have been found to share convergent amino acid replacements along C4 lineages. However, the extent of convergent replacements potentially associated with the emergence of C4 plants remains to be fully assessed. Here, we conducted an organelle-wide analysis to determine if convergent evolution occurred in multiple chloroplast proteins beside the well-known case of the large RuBisCO subunit encoded by the chloroplast gene rbcL.MethodsOur study was based on the comparative analysis of 43 C4 and 21 C3 grass species belonging to the PACMAD clade, a focal taxonomic group in many investigations of C4 evolution. We first used protein sequences of 67 orthologous chloroplast genes to build an accurate phylogeny of these species. Then, we inferred amino acid replacements along 13 C4 lineages and 9 C3 lineages using reconstructed protein sequences of their reference branches, corresponding to the branches containing the most recent common ancestors of C4-only clades and C3-only clades. Pairwise comparisons between reference branches allowed us to identify both convergent and non-convergent amino acid replacements between C4:C4, C3:C3 and C3:C4 lineages.ResultsThe reconstructed phylogenetic tree of 64 PACMAD grasses was characterized by strong supports in all nodes used for analyses of convergence. We identified 217 convergent replacements and 201 non-convergent replacements in 45/67 chloroplast proteins in both C4 and C3 reference branches. C4:C4 branches showed higher levels of convergent replacements than C3:C3 and C3:C4 branches. Furthermore, we found that more proteins shared unique convergent replacements in C4 lineages, with both RbcL and RpoC1 (the RNA polymerase beta’ subunit 1) showing a significantly higher convergent/non-convergent replacements ratio in C4 branches. Notably, more C4:C4 reference branches showed higher numbers of convergent vs. non-convergent replacements than C3:C3 and C3:C4 branches. Our results suggest that, in the PACMAD clade, C4 grasses experienced higher levels of molecular convergence than C3 species across multiple chloroplast genes. These findings have important implications for our understanding of the evolution of the C4 photosynthesis pathway.

Genomic basis of evolutionary adaptation in a warm-blooded fish

Few fishes have evolved elevated body temperatures compared with ambient temperatures, and only in opah (Lampris spp) is the entire body affected. To understand the molecular basis of endothermy, we analyzed the opah genome and identified 23 genes with convergent amino acid substitutions across fish, birds, and mammals, including slc8b1, which encodes the mitochondrial Na+/Ca2+ exchanger and is essential for heart function and metabolic heat production. Among endothermic fishes, 44 convergent genes with suggestive metabolic functions were identified, such as glrx3, encoding a crucial protein for hemoglobin maturation. Numerous genes involved in the production and retention of metabolic heat were also found to be under positive selection. Analyses of opah's unique inner-heat-producing pectoral muscle layer (PMI), an evolutionary key innovation, revealed that many proteins were co-opted from dorsal swimming muscles for thermogenesis and oxidative phosphorylation. Thus, the opah genome provides valuable resources and opportunities to uncover the genetic basis of thermal adaptations in fish.