RELAX: detecting relaxed selection in a phylogenetic framework

Natural resistance to cancers in long-lived mammals: genomic mechanisms and experimental evidence to explain Peto's paradox

Long-lived mammals are reported to have rare or considerably fewer instances of spontaneous tumors, suggesting they might have evolved specific or convergent mechanisms of cancer resistance to extend lifespan; however, the underlying mechanisms remain insufficiently explored. Here, we conducted comparative analysis across 60 mammalian genomes to investigate the genomic features associated with natural cancer resistance. We identified 296 strongly selected genes unique to long-lived species and associated with immune response, DNA repair, and cancer, which might contribute to cancer resistance and lifespan extension in long-lived species. Further, 229 convergent cancer-related genes were detected in the four extremely long-lived species and in-vitro assays confirmed a convergent mutation of LZTS1, shared by bowhead whales and naked mole rats, could suppress cancer development. Importantly, 16 genes were significantly related to both body weight and cancer, defined as candidate genes of Peto's paradox. Of them, the YAP1 gene, harboring the A214S mutation, was identified as a key gene that upregulated tumor suppression genes by localizing to the cytoplasm, which might prohibit cancer development in the large and long-lived cetaceans. These findings provide novel insights into the molecular mechanisms underlying natural cancer resistance in long-lived mammals and the biological basis of Peto's paradox.

Repeated loss of function at HD mating-type genes and of recombination in anther-smut fungi

Basidiomycete fungi typically have two mating-type loci controlling mating compatibility, HD and PR, residing on different chromosomes. Loss-of-function in mating compatibility has been reported at the PR genes in a few heterothallic basidiomycetes, but not for the HD genes. In Microbotryum anther-smut fungi, there have been repeated linkage events between the HD and PR loci through chromosome fusions, leading to non-recombining regions. Here, we found that two sister Microbotryum species parasitizing Dianthus plants, M. superbum and M. shykoffianum, as well as the distantly related M. scorzonerae, have their HD and PR loci on different chromosomes, but with the PR chromosome fused with a part of the ancestral HD chromosome. In addition, recombination suppression has extended stepwise, generating evolutionary strata. In all three species, the HD genes lost their function in mating compatibility, natural diploid strains being often homozygous at the HD locus. Strains could be homozygous for a disrupted HD2 gene, that was hardly expressed during mating. Mating tests confirmed that a single genetic factor controlled mating compatibility and that haploid strains with identical HD alleles could mate and produce hyphae. This study shows that a unifactorial mating-type determinism can evolve, repeatedly, from a bifactorial system, by different mechanisms.

Positive selection and relaxed purifying selection contribute to rapid evolution of sex-biased genes in green seaweed Ulva

BACKGROUND

The evolution of differences in gamete size and number between sexes is a cornerstone of sexual selection theories. The green macroalga Ulva, with incipient anisogamy and parthenogenetic gametes, provides a unique system to investigate theoretical predictions regarding the evolutionary pressures that drive the transition from isogamy to anisogamy, particularly in relation to gamete size differentiation and sexual selection. Its minimal gamete dimorphism and facultative parthenogenesis enable a rare window into early evolutionary steps toward anisogamy.

RESULTS

By analyzing the expression profiles of sex-biased genes (SBGs) during gametogenesis, we found that SBGs evolve faster than unbiased genes, driven by higher rates of non-synonymous substitution (dN), indicating that SBGs are under stronger selective pressures. Mating type minus-biased genes (mt-BGs) exhibit higher dN/dS values than mating type plus-biased genes (mt+BGs), suggesting stronger selective pressures on mt-BGs, although this difference was not statistically significant (P = 0.08). Using branch-site and RELAX models, we found positive selection and relaxed purifying selection acting on a significant proportion of SBGs, particularly those associated with flagella function.

CONCLUSIONS

This study highlights the selective pressures shaping anisogamy and provides insights into the molecular mechanisms underlying its evolution. The faster evolution of SBGs, particularly mt-BGs, and the positive selection on genes associated with motility, such as those related to flagella function, suggest the importance of enhanced gamete motility in the transition to anisogamy. These findings contribute to our understanding of sexual selection and the evolutionary forces that drive the differentiation of gamete size and number between sexes.

Convergent Isobilateral Leaves Increase the Risk for Mangroves Facing Human-Induced Rapid Environmental Changes

Understanding plant adaptations in extreme environments is crucial, as these adaptations often confer advantages for survival. However, a significant gap exists regarding the genetic mechanisms underlying these adaptations and their responses to human-induced rapid environmental change (HIREC). This study addresses the question of whether genetic convergence occurs among plants with similar adaptive features, specifically focusing on isobilateral leaves in mangrove species. Here, we analyse the genetic convergence of isobilateral leaves in mangroves that have independently adapted to coastal intertidal zones. Our findings reveal that genetic convergence is evident in gene families involved in leaf abaxial and adaxial development, with strong selection pressures identified in photosynthesis and leaf polarity pathways. Despite these adaptations, mangrove species with isobilateral leaves occupy narrower ecological niches and face diminishing suitable habitat areas projected under various HIREC scenarios. These results indicate that while convergent traits enhance local adaptation, they may also increase vulnerability to ongoing environmental changes. This research provides valuable insight into the interplay between genetic adaptation and environmental resilience, underscoring the necessity for targeted biodiversity conservation strategies that safeguard specific adaptive traits amid rapid environmental shifts.

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.

Multiple Horizontal Transfers of Immune Genes Between Distantly Related Teleost Fishes

Horizontal gene transfer (HGT) is less frequent in eukaryotes than in prokaryotes, yet can have strong functional implications and was proposed as a causal factor for major adaptations in several eukaryotic lineages. Most cases of eukaryote HGT reported to date are inter-domain transfers, and few studies have investigated eukaryote-to-eukaryote HGTs. Here, we performed a large-scale survey of HGT among 242 species of ray-finned fishes. We found multiple lines of evidence supporting 19 teleost-to-teleost HGT events that involve 17 different genes in 11 teleost fish orders. The genes involved in these transfers show lower synonymous divergence than expected under vertical transmission, their phylogeny is inconsistent with that of teleost fishes, and they occur at non-syntenic positions in donor and recipient lineages. The distribution of HGT events in the teleost tree is heterogenous, with 8 of the 19 transfers occurring between the same two orders (Osmeriformes and Clupeiformes). Though we favor a scenario involving multiple HGT events, future work should evaluate whether hybridization between species belonging to different teleost orders may generate HGT-like patterns. Besides the previously reported transfer of an antifreeze protein, most transferred genes play roles in immunity or are pore-forming proteins, suggesting that such genes may be more likely than others to confer a strong selective advantage to the recipient species. Overall, our work shows that teleost-to-teleost HGT has occurred on multiple occasions, and it will be worth further quantifying these transfers and evaluating their impact on teleost evolution as more genomes are sequenced.

Genomic signatures associated with the evolutionary loss of egg yolk in parasitoid wasps

Trait regression and loss have occurred repeatedly in numerous lineages in response to environmental changes. In parasitoid wasps, a megadiverse group of hymenopteran insects, yolk protein reduction or loss has been observed in many species, likely linked to the transition from ectoparasitism to endoparasitism. However, the genetic basis of this trait and the impact of its loss on genome evolution remain poorly understood. Here, we performed a comparative genomic analysis of 64 hymenopteran insects. The conserved insect yolk protein gene vitellogenin (Vg) underwent five independent loss events in four families, involving 23 of the analyzed endoparasitoid species. Whole-genome alignment suggested that Vg loss occurred during genome rearrangement events. Analysis of Vg receptor gene (VgR) loss, selection, and structural variation in lineages lacking Vg demonstrated functional biases in the patterns of gene loss. The ectoparasitism to endoparasitism transition did not appear to be the primary driver of Vg loss or the subsequent VgR evolution. Together, these findings reveal the genomic changes underlying a unique trait loss in parasitoid wasps. More broadly, this study enhances our understanding of yolk protein loss evolution outside the class Mammalia, highlighting a potential evolutionary trend arising from the availability of an alternative nutrient source for embryonic development.

Detectable episodic positive selection in the virion strand a-strain Maize streak virus genes may have a role in its host adaptation

Maize streak virus (MSV) has four genes: cp, encoding the coat protein; mp, the movement protein; and repA and rep, encoding two distinct replication-associated proteins from an alternatively spliced transcript. These genes play roles in encapsidation, movement, replication, and interactions with the external environment, making them prone to stimuli-driven molecular adaptation. We accomplished selection studies on publicly available curated, recombination-free, complete coding sequences for representative A-strain maize streak virus (MSV-A) cp and mp genes. We found evidence of gene-wide selection in these two MSV genes at specific sites within the genes (cp 1.23% and mp 0.99%). Positively selected sites have amino acids that are 60% hydrophilic and 40% hydrophobic in nature. We found significant evidence of positive selection at branches (cp: 0.76 and mp:1.66%) representing the diversity of MSV-A-strain in South Africa, which is related to the MSV-A-matA isolate (GenBank accession number: AF329881), well disseminated and adapted to the maize plant in sub-Saharan Africa. In the mp gene, selection significantly intensified for the overall diversities of the MSV-A sequences and those more related to the MSV-Mat-A isolate. These findings reveal that despite predominantly undergoing non-diversifying selection, the detectable diversifying positive selection observed in these genes may play a major role in MSV-A host adaptive evolution, ensuring sufficient pathogenicity for onward transmission without killing the host.

De Novo Creation of Two Novel Spliceosomal Introns of RECG1 by Intronization of Formerly Exonic Sequences in Orchidaceae

Spliceosomal introns are a key characteristic of eukaryotic genes. However, the origins and mechanisms of new spliceosomal introns remain elusive, and definitive case studies documenting intron creation are still limited. This study examined the RECG1 genes of 49 land plants, including 21 orchids and 28 non-orchid species. Sequence comparison revealed that the fourth intron of Gastrodia and Platanthera (Orchidaceae) is a newly gained spliceosomal intron, originating from the intronization of former exonic sequences. This intronization event was accompanied by the creation of novel recognizable GT/AG splice sites. In contrast, other orchid species lack the corresponding splice sites in the counterpart regions. Moreover, the secondary and tertiary protein structures implied that the intronization events do not affect the protein function. Given the diverse trophic modes of the two genera, we infer that relaxed selection may have contributed to the fluidity of gene structures. This study provides a typical example of de novo lineage-specific intron creation via intronization in orchids supported by multiple lines of evidence, and the two intronization events occurred independently in the same gene. This research enhances our understanding of gene evolution in orchids and provides valuable insights that may assist the annotation of structurally complex genes.

Stronger Evidence for Relaxed Selection Than Adaptive Evolution in High-elevation Animal mtDNA

Mitochondrial (mt) genes are the subject of many adaptive hypotheses due to the key role of mitochondria in energy production and metabolism. One widespread adaptive hypothesis is that selection imposed by life at high elevation leads to the rapid fixation of beneficial alleles in mtDNA, reflected in the increased rates of mtDNA evolution documented in many high-elevation species. However, the assumption that fast mtDNA evolution is caused by positive selection, rather than relaxed purifying selection, has rarely been tested. Here, we calculated the dN/dS ratio, a metric of nonsynonymous substitution bias, and explicitly tested for relaxed selection in the mtDNA of over 700 species of terrestrial vertebrates, freshwater fishes, and arthropods, with information on elevation and latitudinal range limits, range sizes, and body sizes. We confirmed that mitochondrial genomes of high-elevation taxa have slightly higher dN/dS ratios compared to low-elevation relatives. High-elevation species tend to have smaller ranges, which predict higher dN/dS ratios and more relaxed selection across species and clades, while absolute elevation and latitude do not predict higher dN/dS. We also find a positive relationship between body mass and dN/dS, supporting a role for small effective population size leading to relaxed selection. We conclude that higher mt dN/dS among high-elevation species is more likely to reflect relaxed selection due to smaller ranges and reduced effective population size than adaptation to the environment. Our results highlight the importance of rigorously testing adaptive stories against non-adaptive alternative hypotheses, especially in mt genomes.

Comparative mitogenomic analysis of Chinese cavefish Triplophysa (Cypriniformes: Nemacheilidae): novel gene tandem duplication and evolutionary implications

BACKGROUND

Cavefish exhibit significant morphological changes that result in trade-offs in metabolic requirements and energy utilization in perpetual darkness. As cellular "powerhouses", mitochondria play crucial roles in energy metabolism, suggesting that mitochondrial genes have likely experienced selective pressures during cavefish evolution.

RESULTS

This study presents the first assembly of the complete mitogenome of Triplophysa yangi, a typical cavefish species in China. The mitogenome is 17,068 bp long, marking the longest recorded for the genus Triplophysa, and includes 13 protein-coding genes (PCGs), 2 rRNAs, 25 tRNAs, and a noncoding control region. An ~ 500 bp insertion between ND2 and WANCY regions was observed, comprising a large intact tandem repeat unit (A'-N'-OL'-C') flanked by two unannotated sequences (U1/U2). The evolutionary origin of this repeat unit may involve either in situ duplication events with subsequent functional divergence-where neofunctionalization, subfunctionalization, or pseudogenization drove differential mutation rates between paralogs-or alternatively, horizontal acquisition from exogenous genetic material that became functionally integrated into the ancestral T. yangi mitogenome through co-option mechanisms. Phylogenetic analyses revealed two major clades within Triplophysa-epigean and hypogean lineages-consistent with previous classifications, while cave-restricted species exhibited signs of parallel evolution within the hypogean lineage. Selective pressure analysis indicated that the hypogean lineage (cave-dwelling groups, II & III) have a significantly increased ratio of nonsynonymous to synonymous substitution rates (ω) compared to the epigean lineage (surface-dwelling group, I), suggesting a combination of adaptive selection and relaxed functional constraints in cave-dwelling species.

CONCLUSIONS

The duplication of tRNAs in T. yangi and the potential positive selection sites identified in Triplophysa cavefish further indicated adaptive evolution in mitochondrial PCGs in response to extreme subterranean conditions.

The genomic consequences and persistence of sociality in spiders

In cooperatively breeding social animals, a few individuals account for all reproduction. In some taxa, sociality is accompanied by a transition from outcrossing to inbreeding. In concert, these traits reduce effective population size, potentially rendering transitions to sociality "evolutionarily dead-ends." We addressed this hypothesis in a comparative genomic study in spiders, in which sociality has evolved independently at least 23 times, but social branches are recent and short. We present genomic evidence for the evolutionary dead-end hypothesis in a spider genus with three independent transitions to sociality. We assembled and annotated high-quality, chromosome-level reference genomes from three pairs of closely related social and subsocial Stegodyphus species. We timed the divergence between the social and subsocial species pairs to be from 1.3 million to 1.8 million years. Social evolution in spiders involves a shift from outcrossing to inbreeding and from an equal to a female-biased sex ratio, causing severe reductions in effective population size and decreased efficacy of selection. We show that transitions to sociality only had full effect on purifying selection at 119, 260, and 279 kya, respectively, and follow similar convergent trajectories of progressive loss of diversity and shifts to an increasingly female-biased sex ratio. This almost deterministic genomic response to sociality may explain why social spider lineages do not persist. What causes species extinction is not clear, but either could be selfish meiotic drive eliminating the production of males or could be an inability to retain genome integrity in the face of extremely reduced efficacy of selection.

Plastid-nuclear coevolution of ribosomal protein genes in papilionoid legumes

In plants, cellular function is orchestrated by three distinct genomes located within the nucleus, mitochondrion, and plastid. These genomes are interdependent, requiring tightly coordinated maintenance and expression. Plastids host several multisubunit protein complexes encoded by both the plastid and nuclear genomes. To investigate plastid-nuclear coevolution, this study focused on plastid ribosomal protein genes that are encoded by both plastid and nuclear genomes from 50 taxa across 15 of the 22 early branching major clades of the legume subfamily Papilionoideae. Comparative analysis of substitution rates was conducted across five gene sets: nuclear-encoded plastid-targeted ribosomal protein genes (NuCpRP), nuclear-encoded cytosol-targeted ribosomal genes (NuCyRP), other nuclear-encoded plastid-targeted genes that are not involved in ribosomes (NuCpOT), plastid-encoded ribosomal protein genes (CpRP) and plastid-encoded photosynthesis genes (CpPS).1 Elevated nonsynonymous substitution rates (dN) and ratios of nonsynonymous to synonymous substitution rates (dN/dS; ω) were observed in both CpRP and NuCpRP compared to the other gene sets. Significant differences in dN for CpRP and NuCpRP were found between the papilionoid 50-kb inversion clade and other legumes. Using coevolution statistics and evolutionary rate covariation, strong signals of cytonuclear coevolution were identified, where nonsynonymous substitutions in CpRP and NuCpRP genes co-occur along the same branches of the Papilionoideae phylogeny. Increased ω in a few CpRP genes was due to intensified positive selection whereas most of the CpRP and NuCpRP increased ω was caused by relaxed purifying selection. This pattern not only underscores the role of cytonuclear incompatibility in driving speciation but also highlights its constraints on the genetic enhancement of papilionoid crop species.

The unexpected loss of the 'hunger hormone' ghrelin in true passerines: a game changer in migration physiology

Migratory birds must accumulate large amounts of fat prior to migration to sustain long flights. In passerines, the small body size limits the amount of energy stores that can be transported, and therefore birds undergo cycles of extreme fattening and rapid exhaustion of reserves. Research on these physiological adaptations was rattled by the discovery that birds have lost the main vertebrate regulator of fat deposition, leptin. Recent studies have thus focused on ghrelin, known as 'hunger hormone', a peptide secreted by the gastrointestinal tract to regulate, e.g. food intake and body mass in vertebrates. Studies on domestic species showed that, in birds, ghrelin has effects opposite to those described in mammals such as inhibiting instead of promoting food intake. Furthermore, recent studies have shown that ghrelin administration influences migratory behaviour in passerine birds. Using comparative genomics and immunoaffinity chromatography, we show that ghrelin has been lost in Eupasseres after the basic split from Acanthisitti about 50 Ma. We found that the ghrelin receptor is still conserved in passerines. The maintenance of a functional receptor system suggests that in Eupasserines, another ligand has replaced ghrelin, perhaps to bypass the feedback system that would hinder the large pre-migratory accumulation of subcutaneous fat.

Dynamics of natural selection preceding human viral epidemics and pandemics

Using a phylogenetic framework to characterize natural selection, we investigate the hypothesis that zoonotic viruses require adaptation prior to zoonosis to sustain human-to-human transmission. Examining the zoonotic emergence of Ebola virus, Marburg virus, influenza A virus, SARS-CoV, and SARS-CoV-2, we find no evidence of a change in the intensity of natural selection immediately prior to a host switch, compared with typical selection within reservoir hosts. We conclude that extensive pre-zoonotic adaptation is not necessary for human-to-human transmission of zoonotic viruses. In contrast, the reemergence of H1N1 influenza A virus in 1977 showed a change in selection, consistent with the hypothesis of passage in a laboratory setting prior to its reintroduction into the human population, purportedly during a vaccine trial. Holistic phylogenetic analysis of selection regimes can be used to detect evolutionary signals of host switching or laboratory passage, providing insight into the circumstances of past and future viral emergence.

Evolution of UCP1 Gene and Its Significance to Temperature Adaptation in Rodents

Adaptive thermogenesis comprises shivering thermogenesis dependent on skeletal muscles and non-shivering thermogenesis (NST) mediated by uncoupling protein 1 (UCP1). Although the thermogenic function of UCP1 was adopted early in some placental mammals, positive selection predominantly occurred in the ancestral branches of small-bodied species. Some previous studies have revealed that rodents living in northern or high mountain regions adapt to cold environments by increasing NST, whereas those living in tropical and subtropical regions that are not exposed to cold stress express low concentrations of UCP1, indicating that UCP1 may have evolved to adapt to ambient temperatures. In this study, we explored the evolution of UCP1 and its significance to temperature adaptation by performing detailed evolutionary and statistical analyses on 64 rodents with known genomes. As a result, a total of 71 UCP1 gene sequences were obtained, including 47 intact genes, 22 partial genes, and 2 pseudogenes. Further, 47 intact genes and 3 previously published intact UCP1 genes were incorporated into evolutionary analyses, and correlation analyses between evolutionary rate and ambient temperatures (including average annual temperature, maximum temperature, and minimum temperature) of the rodent survives were conducted. The results show that UCP1 is under purifying selection (ω = 0.11), and among rodents with intact UCP1 sequences, Urocitellus parryii and Dicrostonyx groenlandicus-the two species with the lowest ambient temperatures among the rodents used here-have higher evolutionary rates than others. In the statistical analyses, in addition to ambient temperatures, body weight and weight at birth were also taken into account since weight was previously proposed to be linked to UCP1 evolution. The results showed that after controlling for the phylogenetic effect, the maximum temperature was significantly negatively correlated with the evolutionary rate of UCP1, whereas weight did not have a relationship with UCP1 evolutionary rate. Consequently, it is suggested that ambient temperature can drive the evolution of rodent UCP1, thereby enhancing NST adaptation to cold stress.

Plastome evolution in Santalales involves relaxed selection prior to loss of ndh genes and major boundary shifts of the inverted repeat

BACKGROUND AND AIMS

Biological aspects of haustorial parasitism have significant effects on the configuration of the plastid genome. Approximately half the diversity of haustorial parasites belongs to the order Santalales, where a clearer picture of plastome evolution in relation to parasitism is starting to emerge. However, in previous studies of plastome evolution there is still a notable under-representation of members from non-parasitic and deep-branching hemiparasitic lineages, limiting evolutionary inference around the time of transition to a parasitic lifestyle. To expand taxon sampling relevant to this transition we therefore targeted three families of non-parasites (Erythropalaceae, Strombosiaceae and Coulaceae), two families of root-feeding hemiparasites (Ximeniaceae and Olacaceae) and two families of uncertain parasitic status (Aptandraceae and Octoknemaceae). With data from these lineages we aimed to explore plastome evolution in relation to the evolution of parasitism.

METHODS

From 29 new samples we sequenced and annotated plastomes and the nuclear ribosomal cistron. We examined phylogenetic patterns, plastome evolution, and patterns of relaxed or intensified selection in plastid genes. Available transcriptome data were analysed to investigate potential transfer of infA to the nuclear genome.

RESULTS

Phylogenetic relationships indicate a single functional loss of all plastid ndh genes (ndhA-K) in a clade formed by confirmed parasites and Aptandraceae, and the loss coincides with major size and boundary shifts of the inverted repeat (IR) region. Depending on an autotrophic or heterotrophic lifestyle in Aptandraceae, plastome changes are either correlated with or pre-date the evolution of parasitism. Phylogenetic patterns also indicate repeated loss of infA from the plastome, and based on the presence of transcribed sequences with presequences corresponding to thylakoid luminal transit peptides, we infer that the genes were transferred to the nuclear genome.

CONCLUSIONS

Except for the loss of the ndh complex, relatively few genes have been lost from the plastome in deep-branching root parasites in Santalales. Prior to loss of the ndh genes, they show signs of relaxed selection indicative of their dispensability. To firmly establish a potential correlation between ndh gene loss, plastome instability and evolution of parasitism, it is pertinent to refute or confirm a parasitic lifestyle in all Santalales clades.

Recombination and retroprocessing in broomrapes reveal a universal roadmap for mitochondrial evolution in heterotrophic plants

The altered life history strategies of heterotrophic organisms often leave a profound genetic footprint on energy metabolism related functions. In parasitic plants, the reliance on host-derived nutrients and loss of photosynthesis in holoparasites have led to highly degraded to absent plastid genomes, but its impact on mitochondrial genome (mitogenome) evolution has remained controversial. By examining mitogenomes from 45 Orobanchaceae species including three independent transitions to holoparasitism and key evolutionary intermediates, we identified measurable and predictable genetic alterations in genomic shuffling, RNA editing, and intracellular (IGT) and horizontal gene transfer (HGT) en route to a nonphotosynthetic lifestyle. In-depth comparative analyses revealed DNA recombination and repair processes, especially RNA-mediated retroprocessing, as significant drivers for genome structure evolution. In particular, we identified a novel RNA-mediated IGT and HGT mechanism, which has not been demonstrated in cross-species and inter-organelle transfers. Based on this, we propose a generalized dosage effect mechanism to explain the biased transferability of plastid DNA to mitochondria across green plants, especially in heterotrophic lineages like parasites and mycoheterotrophs. Evolutionary rates scaled with these genomic changes, but the direction and strength of selection varied substantially among genes and clades, resulting in high contingency in mitochondrial genome evolution. Finally, we describe a universal roadmap for mitochondrial evolution in heterotrophic plants where increased recombination and repair activities, rather than relaxed selection alone, lead to differentiated genome structure compared to free-living species.

Blue light restores functional circadian clocks in eyeless cave spiders

Evolution in profound darkness often leads to predictable, convergent traits, such as the loss of vision. Yet, the consequences of such repeated evolutionary experiments remain obscure, especially regarding fundamental regulatory behaviors like circadian rhythms. We studied circadian clocks of blind cave spiders and their sighted relatives. In the field, cave spiders exhibit low per expression and maintain constant activity levels. Curiously, their clocks are not permanently lost; exposure to monochromatic blue light restores both circadian gene expression and behavioral rhythms. Conversely, blocking blue light in sighted relatives induces an arrhythmic "cave phenotype." Our RNA interference experiments suggest that clock genes regulate the rhythmicity of the huddle response, establishing a link between circadian gene networks and this behavioral rhythm. We demonstrate that circadian regulation is readily toggled and may play a latent role, even in constant darkness. Overall, our study expands understanding of circadian clock variations and paves the way for future research on the maintenance of silent phenotypes.

Primary regulatory T cell activator FOXP3 is present across Amphibia

The overall structure of the immune system is highly conserved across jawed vertebrates, but characterization and description of the immune system is heavily biased toward mammals. One arm of the vertebrate immune system, the adaptive immune system, mounts pathogen-specific responses that tend to be robust and effective at clearing pathogens. This system requires selection against self-recognition and modulation of the immune response. One of the mechanisms of immune modulation is the presence of regulatory T cells that suppress other effector immune cells. Regulatory T cells and their primary activator forkhead box protein P3 (FOXP3) have been well characterized in mammalian models but unexplored in most other vertebrate taxa. Amphibians are a good focal group for the characterization of FOXP3 due to their phylogenetic position on the vertebrate tree of life, and their susceptibility to emerging pathogens. In this study, we mined available transcriptomic and genomic data to confirm the presence of FOXP3 across the amphibian tree of life. We find that FOXP3 is present in all major clades of amphibians. We also test whether selection on FOXP3 shows signatures of intensification among the three main clades of amphibians, which may reflect shifts in the stringency of natural selection on this gene. Our findings provide insights into the evolutionary history of the vertebrate immune system and confirm the conservation of vertebrate immune genes within amphibians.

Elucidating the evolutionary dynamics of parasitism in Cuscuta: in-depth phylogenetic reconstruction and extensive plastomes reduction

BACKGROUND

The genus Cuscuta L. (Convolvulaceae), commonly known as dodder, is a holoparasite plant that relies on host plants for nutrition, leading to significant genomic changes, particularly in plastomes. This dependency has led to significant reductions and modifications in their plastomes compared to autotrophic plants. In contrast to the well-conserved plastomes of photosynthetic plants, Cuscuta exhibits substantial genomic reductions reflecting the loss of photosynthetic functions and associated genes.

RESULT

This study examines eight plastomes within Cuscuta and reconstructs the phylogenetic relationships among 40 Cuscuta taxa using five other genera as an outgroup. The size of plastid genome varies significantly, with the smallest being 60 kb and the largest 121 kb, highlighting extensive genomic reduction. In special cases, the subgenera Cuscuta exhibit the loss of inverted repeats, distinguishing from them other subge within the Cuscuta genus. This reduction is most pronounced in genes related to photosynthesis, such as atp, pet, psa, psb, and ycf genes, particularly in the subg. Grammica (Lour.) Peter. The study also notes the frequent and independent loss of the plastid genes infA, rpl23, rpl32, rps15, and rps16 across various angiosperm lineages, often involving transfer to the nuclear genome. In parasitic plants like Cuscuta, the ndh genes, crucial for photosynthesis, are often lost. The study also highlights that in the subg. Grammica, the matK and rpo genes, along with trnR-ACG genes, are lost in parallel, indicating that these parasitic plants do not need matK and rpo genes after the loss of ndh genes for survival. Analysis of selective relaxation pressure on plastid genes shows a reductive trend, with genes such as atp, pet, psa, psb, rpo, and ycf progressively becoming pseudogenes over time, with housekeeping genes like rpl and rps expected to follow. However, the pseudogenization process is specific to the subg. Grammica, Pachystigma (Engelm.) Baker & C.H.Wright, and Cuscuta, rather than in the subg. Monogynella (Des Moul.) Peter, Engl. & Prantl (ancient clade species).

CONCLUSION

The study of Cuscuta plastomes reveals the profound impact of parasitism on genome evolution, highlighting the complex interplay of gene retention and loss through phylogenomic approaches. This research enriches our understanding of plant genome evolution and the intricate host-parasite relationships. It also sheds light on the evolutionary history and genomic adaptations of Cuscuta, illustrating the diverse strategies enabling subg. Grammica, Pachystigma, Cuscuta, and Monogynella thrive as parasitic species. These findings provide valuable insights into the molecular mechanisms underlying parasitism and its impact on plastid genome organization.

Whole-genome relaxed selection and molecular constraints in Triplophysa under adapted Qinghai-Tibetan Plateau

High-altitude environments are inhospitable, but Triplophysa, the largest taxon among the three major fish groups in the Qinghai-Tibetan Plateau (QTP), is an exception. However, the evolutionary profiling of the common ancestor and its contribution to the adaptation of existing QTP native species is unclear.We researched the comparative genomics of Triplophysa species and found that the genome-wide genes of Triplophysa and its ancestry have the characteristics of rapid evolution.Moreover, the rapid evolution of the ancestral genes was caused by relaxed selection. Natural selection analysis showed that more ancestral relaxed selection genes were under strongly purifying selection and showed higher expression in QTP endemic Triplophysa species.The change in natural selection might be associated with the adaptation to QTP. It should be noted that SPT5 homolog, DSIF elongation factor subunit (supt5h) experienced relaxed selection in common ancestral populations of Triplophysa but under purifying selection in extant species, which might be related to hypoxia adaptation of QTP. In summary, the extant species in different environments were used to infer the evolutionary profile of the common ancestor and to identify candidate genes based on changes in natural selection. Our work might provide new clues for understanding adaptation to extreme environments.

Caecilians maintain a functional long-wavelength-sensitive cone opsin gene despite signatures of relaxed selection and more than 200 million years of fossoriality

Visual systems are tuned to animals' ecologies, evolving in response to specific light environments and visual needs. Ecological transitions to fossorial lifestyles impose strong selective pressures favoring morphological adaptations for underground life, such as increased skull ossification and reduced eye protrusion. Fossoriality may simultaneously relax constraints on other aspects of vision leading to diminished visual capabilities. Caecilians (Gymnophiona)-specialized, fossorial amphibians-possess reduced eyes covered by skin or bone. For years, these traits, along with the presence of a single photoreceptor expressing one functional opsin gene, have been interpreted as evidence of limited visual capabilities, including an inability to focus or perceive color. Our results challenge these assumptions: we identified the long-wavelength-sensitive (LWS) opsin gene in 11 species of caecilians spanning 8 of 10 recognized families. Molecular evidence indicates that LWS is intact and transcribed in the eye of at least one species (Caecilia orientalis). Anatomical observations from five caecilian families indicate highly organized retinae even in families with vestigial eyes. While the presence of cone cells in our study species remains uncertain, a putatively functional LWS gene suggests that the visual capabilities of caecilians and the role of light perception in their ecology may be underestimated.

Homospermidine synthase evolution and the origin(s) of pyrrolizidine alkaloids in Apocynaceae

PREMISE

Enzymes that are encoded by paralogous genes and produce identical specialized metabolites in distantly related plant lineages are strong evidence of parallel phenotypic evolution. Inference of phenotypic homology for metabolites produced by orthologous genes is less straightforward, since orthologs may be recruited in parallel into novel pathways. In prior research on pyrrolizidine alkaloids (PAs), specialized metabolites of Apocynaceae, the evolution of homospermidine synthase (HSS), an enzyme of PA biosynthesis, was reconstructed and a single origin of PAs inferred because HSS enzymes of all known PA-producing Apocynaceae species are orthologous and descended from an ancestral enzyme with the motif (VXXXD) of an optimized HSS.

METHODS

We increased sampling, tested the effect of amino acid motif on HSS function, revisited motif evolution, and tested for selection to infer evolution of HSS function and its correlation with phenotype.

RESULTS

Some evidence supports a single origin of PAs: an IXXXD HSS-like gene, similar in function to VXXXD HSS, evolved in the shared ancestor of all PA-producing species; loss of HSS function occurred multiple times via pseudogenization and perhaps via evolution of an IXXXN motif. Other evidence indicates multiple origins: the VXXXD motif, highly correlated with the PA phenotype, evolved two or four times independently; the ancestral IXXXD gene was not under positive selection, while some VXXXD genes were; and substitutions at sites experiencing positive selection occurred on multiple branches in the HSS-like gene tree.

CONCLUSIONS

The complexity of the genotype-function-phenotype map confounds the inference of PA homology from HSS-like gene evolution in Apocynaceae.

Characterization of Complete Mitochondrial Genome of Badri Breed of Bos indicus (Bovidae: Bovinae): Selection Pressure and Comparative Analysis

High-altitude mammals are often subject to specific environmental obstacles, which exert selective pressure on their physiological and morphological traits, hence driving their evolutionary processes. It is anticipated that these circumstances will lead to the adaptive evolution of protein-coding genes (PCGs) in the mitochondrial genome, which play a crucial role in the oxidative phosphorylation system. In this study, we have generated the complete mitochondrial genome of the Badri breed of Bos indicus inhabiting a high-altitude environment to test the signatures of adaptive evolution on PCGs and their phylogenetic relationships. The complete mitogenome of the Badri breed is 16,339 bp and most tRNAs showed typical clover-leaf secondary structure with a few exceptions, like trnS1 and trnS2 without DHU arm and trnK without DHU loop. Comparative analysis of PCGs indicated that cox1 is the most conserved, while atp6 is the most variable gene. Moreover, the ratios of non-synonymous to synonymous substitution rates indicated the purifying selection (Ka/Ks < 1) in the protein-coding genes that shape the diversity in mitogenome of Bos indicus. Furthermore, Branch-site model (BSM) suggested that cox1, cox2, nad3, nad4L, and nad6 underwent stronger purifying selection (ω < 1) than other PCGs in 15 breeds of 4 species, including Badri. BSM also detected 10 positive sites in PCGs and one in 13 PCGs concatenated dataset. Additional analyses in Datamonkey indicated 11 positive sites and 23 purifying sites in the concatenated dataset, a relaxation of selection strength in nad3, and no evidence of episodic diversifying selection in any PCGs. Phylogeny revealed the sister relationship of the Badri with other breeds of Bos indicus as well as Bos frontalis (Gayal-2). The mitogenome of the Badri breed is an important genomic resource for conservation genetics of this species and also contributes to the understanding of the adaptive evolution of mitochondrial protein coding genes.

Hymenopteran-specific TRPA channel from the Texas leaf cutter ant (Atta texana) is heat and cold activated and expression correlates with environmental temperature

Leaf cutting ants of the genus Atta cultivate fungal gardens, carefully modifying environmental conditions to maintain optimal temperature for fungal growth. Antennal nerves from Atta are highly temperature sensitive, but the underlying molecular sensor is unknown. Here, we utilize Atta texana (Texas leaf cutter ant) to investigate the molecular basis of ant temperature sensation and how it might have evolved as the range expanded northeast across Texas from ancestral populations in Mexico. We focus on transient receptor potential (TRP) channel genes, the best characterized temperature sensor proteins in animals. Atta texana antennae express 6 of 13 Hymenopteran TRP channel genes and sequences are under a mix of relaxed and intensified selection. In a behavioral assay, we find A. texana workers prefer 24 °C (range 21-26 °C) for fungal growth. There was no evidence of regulatory evolution across a temperature transect in Texas, but instead Hymenoptera-specific TRPA (HsTRPA) expression highly correlated with ambient temperature. When expressed in vitro, HsTRPA from A. texana is temperature activated with Q10 values exceeding 100 on initial exposure to temperatures above 33 °C. Surprisingly, HsTRPA also appears to be activated by cooling, and therefore to our knowledge, the first non-TRPA1 ortholog to be described with dual heat/cold activation and the first in any invertebrate.

Hierarchical architecture of neo-sex chromosomes and accelerated adaptive evolution in tortricid moths

Sex chromosomes can expand through fusion with autosomes, thereby acquiring unique evolutionary patterns. In butterflies and moths (Lepidoptera), these sex chromosome-autosome (SA) fusions occur relatively frequently, suggesting possible evolutionary advantages. Here, we investigated how SA fusion affects chromosome features and molecular evolution in leafroller moths (Lepidoptera: Tortricidae). Phylogenomic analysis showed that Tortricidae diverged ∼124 million years ago, accompanied by an SA fusion between the Merian elements M(20 + 17) and MZ. In contrast to partial autosomal fusions, the fused neo-Z Chromosome developed a hierarchical architecture, in which the three elements exhibit heterogeneous sequence features and evolutionary patterns. Specifically, the M17 part had a distinct base composition and chromatin domains. Unlike M20 and MZ, M17 was expressed at the same levels as autosomes in both sexes, compensating for the lost gene dosage in females. Concurrently, the SA fusion drove M17 as an evolutionary hotspot, accelerating the evolution of several genes related to ecological adaptation (e.g., ABCCs) and facilitating the divergence of closely related species, whereas the undercompensated M20 did not show such an effect. Thus, accelerated evolution under a novel pattern of dosage compensation may have favored the adaptive radiation of this group. This study demonstrates the association between a karyotype variant and adaptive evolution and explains the recurrent SA fusion in the Lepidoptera.

Dampened TLR2-mediated Inflammatory Signaling in Bats

Bats are considered natural hosts for numerous viruses. Their ability to carry viruses that cause severe diseases or even death in other mammals without falling ill themselves has attracted widespread research attention. Toll-like receptor 2 forms heterodimers with Toll-like receptor 1 or Toll-like receptor 6 on cell membranes, recognizing specific pathogen-associated molecular patterns and playing a key role in innate immune responses. Previous studies have shown that moderate Toll-like receptor 2-mediated immune signals aid in pathogen clearance, while excessive or inappropriate Toll-like receptor 2-mediated immune signals can cause self-damage. In this study, we observed that TLR2, unlike TLR1 or TLR6, has undergone relaxed selection in bats compared with other mammals, indicating a reduced functional constraint on TLR2 specifically in bats. Indeed, our cell-based functional assays demonstrated that the ability of Toll-like receptor 2 to bind with Toll-like receptor 1 or Toll-like receptor 6 was significantly reduced in bats, leading to dampened inflammatory signaling. We identified mutations unique to bats that were responsible for this observation. Additionally, we found that mutations at residues 375 and 376 of Toll-like receptor 2 in the common ancestor of bats also resulted in reduced inflammatory response, suggesting that this reduction occurred early in bat evolution. Together, our study reveals that the Toll-like receptor 2-mediated inflammatory response has been specifically dampened in bats, which may be one of the reasons why they could harbor many viruses without falling ill.

Functional differences between rodent and human PD-1 linked to evolutionary divergence

Mechanistic understanding of the inhibitory immunoreceptor PD-1 is largely based on mouse models, but human and mouse PD-1 share only 59.6% amino acid identity. Here, we found that human PD-1 is more inhibitory than mouse PD-1, owing to stronger interactions with the ligands PD-L1 and PD-L2 and more efficient recruitment of the effector phosphatase Shp2. In a mouse melanoma model with adoptively transferred T cells, humanization of a PD-1 intracellular domain disrupted the antitumor activity of CD8+ T cells and increased the magnitude of anti-PD-1 response. We identified a motif highly conserved across vertebrate PD-1 orthologs, absent in rodents, as a key determinant for differential Shp2 recruitment. Evolutionary analysis suggested that PD-1 underwent a rodent lineage-specific functional attenuation during evolution. Together, our study uncovers species-specific features of the PD-1 pathway, with implications for PD-1 evolution and differential anti-PD-(L)1 responses in mouse models and human patients.

Organellar phylogenomics at the epidendroid orchid base, with a focus on the mycoheterotrophic Wullschlaegelia

BACKGROUND AND AIMS

Heterotrophic plants have long been a challenge for systematists, exemplified by the base of the orchid subfamily Epidendroideae, which contains numerous mycoheterotrophic species.

METHODS

Here we address the utility of organellar genomes in resolving relationships at the epidendroid base, specifically employing models of heterotachy, or lineage-specific rate variation over time. We further conduct comparative analyses of plastid genome evolution in heterotrophs and structural variation in matK.

KEY RESULTS

We present the first complete plastid genomes (plastomes) of Wullschlaegelia, the sole genus of the tribe Wullschlaegelieae, revealing a highly reduced genome of 37 kb, which retains a fraction of the genes present in related autotrophs. Plastid phylogenomic analyses recovered a strongly supported clade composed exclusively of mycoheterotrophic species with long branches. We further analysed mitochondrial gene sets, which recovered similar relationships to those in other studies using nuclear data, but the placement of Wullschlaegelia remains uncertain. We conducted comparative plastome analyses among Wullschlaegelia and other heterotrophic orchids, revealing a suite of correlated substitutional and structural changes relative to autotrophic species. Lastly, we investigated evolutionary and structural variation in matK, which is retained in Wullschlaegelia and a few other 'late stage' heterotrophs and found evidence for structural conservation despite rapid substitution rates in both Wullschlaegelia and the leafless Gastrodia.

CONCLUSIONS

Our analyses reveal the limits of what the plastid genome can tell us on orchid relationships in this part of the tree, even when applying parameter-rich heterotachy models. Our study underscores the need for increased taxon sampling across all three genomes at the epidendroid base, and illustrates the need for further research on addressing heterotachy in phylogenomic analyses.

Convergent relaxation of molecular constraint in herbivores reveals the changing role of liver and kidney functions across mammalian diets

Mammalia comprises a great diversity of diet types and associated adaptations. An understanding of the genomic mechanisms underlying these adaptations may offer insights for improving human health. Comparative genomic studies of diet that employ taxonomically restricted analyses or simplified diet classifications may suffer reduced power to detect molecular convergence associated with diet evolution. Here, we use a quantitative carnivory score-indicative of the amount of animal protein in the diet-for 80 mammalian species to detect significant correlations between the relative evolutionary rates of genes and changes in diet. We have identified six genes-ACADSB, CLDN16, CPB1, PNLIP, SLC13A2, and SLC14A2-that experienced significant changes in evolutionary constraint alongside changes in carnivory score, becoming less constrained in lineages evolving more herbivorous diets. We further consider the biological functions associated with diet evolution and observe that pathways related to amino acid and lipid metabolism, biological oxidation, and small molecule transport experienced reduced purifying selection as lineages became more herbivorous. Liver and kidney functions show similar patterns of constraint with dietary change. Our results indicate that these functions are important for the consumption of animal matter and become less important with the evolution of increasing herbivory. So, genes expressed in these tissues experience a relaxation of evolutionary constraint in more herbivorous lineages.

Comparative genomics of the sexually transmitted parasite Trichomonas vaginalis reveals relaxed and convergent evolution and genes involved in spillover from birds to humans

Trichomonas vaginalis is the causative agent of the venereal disease trichomoniasis which infects men and women globally and is associated with serious outcomes during pregnancy and cancers of the human reproductive tract. Trichomonads parasitize a range of hosts in addition to humans including birds, livestock, and domesticated animals. Recent genetic analysis of trichomonads recovered from columbid birds has provided evidence that these parasite species undergo frequent host-switching, and that a current epoch spillover event from columbids likely gave rise to T. vaginalis in humans. We undertook a comparative evolutionary genomics study of seven trichomonad species, generating chromosome-scale reference genomes for T. vaginalis and its avian sister species Trichomonas stableri, and assemblies of five other species that infect birds and mammals. Human-infecting trichomonad lineages have undergone recent and convergent genome size expansions compared to their avian sister species, and the major contributor to their increased genome size is increased repeat expansions, especially multicopy gene families and transposable elements, with genetic drift likely a driver due to relaxed selection. Trichomonads have independently host-switched twice from birds to humans, and genes implicated in the transition to the human host include those associated with host tissue adherence and phagocytosis, extracellular vesicles, and CAZyme virulence factors.

Sperm competition intensity shapes divergence in both sperm morphology and reproductive genes across murine rodents

It remains unclear how variation in the intensity of sperm competition shapes phenotypic and molecular evolution across clades. Mice and rats in the subfamily Murinae are a rapid radiation exhibiting incredible diversity in sperm morphology and production. We combined phenotypic and genomic data to perform phylogenetic comparisons of male reproductive traits and genes across 78 murine species. We identified several shifts towards smaller relative testes mass (RTM), presumably reflecting reduced sperm competition. Several sperm traits were associated with RTM, suggesting that mating system evolution selects for convergent suites of traits related to sperm competitive ability. We predicted that sperm competition would also drive more rapid molecular divergence in species with large testes. Contrary to this, we found that many spermatogenesis genes evolved more rapidly in species with smaller RTM due to relaxed purifying selection. While some reproductive genes evolved rapidly under recurrent positive selection, relaxed selection played a greater role in underlying rapid evolution in small testes species. Our work demonstrates that postcopulatory sexual selection can impose strong purifying selection shaping the evolution of male reproduction and that broad patterns of molecular evolution may help identify genes that contribute to male fertility.

Gene expansion in the hawkmoth Manduca sexta drives evolution of food-associated odorant receptors

In insects, odorant receptors (ORs) are required for the detection of most olfactory cues. We investigated the function of a clade of four duplicated ORs in the hawkmoth Manduca sexta and found that these paralogs encode broadly tuned receptors with overlapping but distinct response spectra. Two paralogs, which arose after divergence from a related lineage, show high sensitivity to floral esters released by a nectar-rich plant frequently visited by M. sexta. Functional imaging in mutant moths lacking one of the paralogs suggests that olfactory sensory neurons expressing this OR target a previously identified feeding-associated glomerulus in the primary olfactory center of the brain. However, only the response of this glomerulus to the single ligand unique to the now mutated OR disappeared, suggesting neuronal coexpression of the paralogs. Our results suggest a link between the studied OR expansion and enhanced detection of odors emitted by valuable nectar sources in M. sexta.

Genome dynamics across the evolutionary transition to endosymbiosis

Endosymbiosis-where a microbe lives and replicates within a host-is an important contributor to organismal function that has accelerated evolutionary innovations and catalyzed the evolution of complex life. The evolutionary processes associated with transitions to endosymbiosis, however, are poorly understood. Here, we leverage the wide diversity of host-associated lifestyles of the genus Arsenophonus to reveal the complex evolutionary processes that occur during the transition to a vertically transmitted endosymbiotic lifestyle from strains maintained solely by horizontal (infectious) transmission. We compared the genomes of 38 strains spanning diverse lifestyles from horizontally transmitted pathogens to obligate interdependent endosymbionts. Among culturable strains, we observed those with vertical transmission had larger genome sizes than closely related horizontally transmitting counterparts, consistent with evolutionary innovation and the rapid gain of new functions. Increased genome size was a consequence of prophage and plasmid acquisition, including a cargo of type III effectors, alongside the concomitant loss of CRISPR-Cas genome defense systems, enabling mobile genetic element expansion. Persistent endosymbiosis was also associated with loss of type VI secretion, which we hypothesize to be a consequence of reduced microbe-microbe competition. Thereafter, the transition to endosymbiosis with strict vertical inheritance was associated with the expected relaxation of purifying selection, gene pseudogenization, metabolic degradation, and genome reduction. We argue that reduced phage predation in endosymbiotic niches drives the loss of genome defense systems driving rapid genome expansion upon the adoption of endosymbiosis and vertical transmission. This remodeling enables rapid horizontal gene transfer-mediated evolutionary innovation and precedes the reductive evolution traditionally associated with adaptation to endosymbiosis.

Genomic Insights Into Genetic Basis of Evolutionary Conservatism and Innovation in Frogs

Examining closely related species evolving in similar environments offers valuable insights into the mechanisms driving phylogenetic conservatism and evolutionary lability. This can elucidate the intricate relationship between inheritance and environmental factors. Nonetheless, the precise genomic dynamics and molecular underpinnings of this process remain enigmatic. This study explores the evolutionary conservatism and adaptation exhibited by two closely related high-altitude frog species: Nanorana parkeri and N. pleskei. We assembled a high-quality genome for Tibetan N. pleskei and compared it to the genomes of N. parkeri and their lowland relatives. Our findings reveal that these two Tibetan frog species diverged approximately 16.6 million years ago, pointing to a possible ancestral colonization of high-elevation habitats. Following this colonization, significant adaptive evolution occurred in both coding and non-coding regions of the ancestral lineage. This evolution led to notable phenotypic alterations, as evidenced by the reduced body size. Also, due to purifying selection, most ancestral adaptive features persisted in descendant species, indicating a strong element of evolutionary conservatism. However, descendant species evolved novel adaptations to exacerbated environmental challenges in the Tibet Plateau, mainly related to hypoxia response. Furthermore, our analysis underscores the critical role of regulatory variations in descendant adaptive evolution. Notably, hub genes in networks, such as EGLN3, accumulated more variations in regulatory regions as they were transmitted from ancestors to descendants. In sum, our study sheds light on the profound and lasting impact of genetic heritage on species' adaptive evolution.

Adaptation During the Shift from Entomopathogen to Endosymbiont Is Accompanied by Gene Loss and Intensified Selection

Fungi have been found to be associated with many insect species, with some species transitioning to reside within insects as symbionts. However, the evolutionary pressures and genomic consequences associated with this transition are not well understood. Pathogenic fungi of the genus Ophiocordyceps have undergone multiple, independent transitions from pathogen to endosymbiont lifestyles, where they reside within the fatty tissues of infected soft-scale insects transgenerationally without killing their hosts. To gain an understanding of the genomic adaptations underlying this life history shift, long-read sequencing was utilized to assemble the genomes of both the soft-scale insect Parthenolecanium corni and its Ophiocordyceps endosymbiont from a single insect. Assembly and metagenomic-based binning produced a highly contiguous genome for Pa. corni and a chromosome-level assembly for the Ophiocordyceps endosymbiont. The endosymbiont genome was characterized by 524 gene loss events compared to free-living pathogenic Ophiocordyceps relatives, with predicted roles in hyphal growth, cell wall integrity, metabolism, gene regulation, and toxin production. Contrasting patterns of selection were observed between the nuclear and mitochondrial genomes specific to the endosymbiont lineage. Intensified selection was most frequently observed across orthologs in the nuclear genome, whereas selection on most mitochondrial genes was found to be relaxed. Scans for positive selection were enriched within the fatty acid metabolism pathway with endosymbiont specific selection within three adjacent enzymes catalyzing the conversion of acetoacetate to acetyl-coenzyme A, suggesting that the endosymbiont lineage is under selective pressure to effectively exploit the lipid rich environment of the insect fat bodies in which it is found.

Genomic signatures of sensory adaptation and evolution in pangolins

BACKGROUND

Pangolin is one of the most endangered mammals with many peculiar characteristics, yet the understanding of its sensory systems is still superficial. Studying the genomic basis of adaptation and evolution of pangolin's sensory system is expected to provide further potential assistance for their conservation in the future.

RESULTS

In this study, we performed a comprehensive comparative genomic analysis to explore the signature of sensory adaptation and evolution in pangolins. By comparing with the aardvark, Cape golden mole, and short-beaked echidna, 124 and 152 expanded gene families were detected in the genome of the Chinese and Malayan pangolins, respectively. The enrichment analyses showed olfactory-related genomic convergence among five concerned mammals. We found 769 and 733 intact OR genes, and 704 and 475 OR pseudogenes in the Chinese and Malayan pangolin species, respectively. Compared to other mammals, far more intact members of OR6 and OR14 were identified in pangolins, particularly for four genes with large copy numbers (OR6C2, OR14A2, OR14C36, and OR14L1). On the genome-wide scale, 1,523, 1,887, 1,110, and 2,732 genes were detected under positive selection (PSGs), intensified selection (ISGs), rapid evolution (REGs), and relaxed selection (RSGs) in pangolins. GO terms associated with visual perception were enriched in PSGs, ISGs, and REGs. Those related to rhythm and sound perception were enriched in both ISGs and REGs, ear development and morphogenesis were enriched in ISGs, and mechanical stimulus and temperature adaptation were enriched in RSGs. The convergence of two vision-related PSGs (OPN4 and ATXN7), with more than one parallel substituted site, was detected among five concerned mammals. Additionally, the absence of intact genes of PKD1L3, PKD2L1, and TAS1R2 and just six single-copy TAS2Rs (TAS2R1, TAS2R4, TAS2R7, TAS2R38, TAS2R40, and TAS2R46) were found in pangolins. Interestingly, we found two large insertions in TAS1R3, distributed in the N-terminal ectodomain, just in pangolins.

CONCLUSIONS

We found new features related to the adaptation and evolution of pangolin-specific sensory characteristics across the genome. These are expected to provide valuable and useful genome-wide genetic information for the future breeding and conservation of pangolins.

Convergent Evolution Associated with the Loss of Developmental Diapause May Promote Extended Lifespan in Bees

Diapause has long been proposed to play a significant role in the evolution of eusociality in Hymenoptera. Recent studies have shown that shifts in the diapause stage precede social evolution in wasps and bees; however, the genomic basis remains unknown. Given the overlap in molecular pathways that regulate diapause and lifespan, we hypothesized that the evolutionary loss of developmental diapause may lead to extended lifespan among adults, which is a prerequisite for the evolution of eusociality. To test whether the loss of prepupal diapause is followed by genomic changes associated with lifespan extension, we compared 27 bee genomes with or without prepupal diapause. Our results point to several potential mechanisms for lifespan extension in species lacking prepupal diapause, including the loss of the growth hormone PTTH and its receptor TORSO, along with convergent selection in genes known to regulate lifespan in animals. Specifically, we observed purifying selection of prolongevity genes and relaxed selection of antilongevity genes within the IIS/TOR pathway in species that have lost prepupal diapause. Changes in selection pressures on this pathway may lead to the evolution of new phenotypes, such as lifespan extension and altered responses to nutritional signals that are crucial for social evolution.

Molecular evolution and genetic diversity of defective chorion 1 in Anastrepha fraterculus and Anastrepha obliqua (Diptera, Tephritidae)

The family Tephritidae comprises numerous fruit fly species, some of which are economically significant, such as several in the genus Anastrepha. Most pest species in this genus belong to the fraterculus group, characterized by closely related species that are difficult to differentiate due to recent divergence and gene flow. Identifying genetic markers for their study is paramount for understanding the group's evolution and eventual phytosanitary control. Because there is variation in eggshell morphology among species in the genus, the study of the rapidly evolving defective chorion 1 (dec-1) gene, which is crucial for chorion formation and reproduction, could provide relevant information for Anastrepha differentiation. We compared transcriptome sequences of dec-1 from two of the most important pest species in the genus, Anastrepha fraterculus and Anastrepha obliqua to dec-1 sequences from Anastrepha ludens, which was used for structure prediction. Furthermore, we amplified a conserved exon across populations of these species. These data revealed three alternative transcripts in A. fraterculus and A. obliqua, consistent with patterns found in other Tephritidae; we obtained orthologous sequences for these other tephritids from NCBI to investigate patterns of selection affecting this gene at different hierarchical levels using different methods. These analyses show a general pattern of purifying selection across the whole gene and throughout its history at different hierarchical levels, from populations to more distantly related species. That notwithstanding, we still found evidence of positive and episodic diversifying selection at different levels. Different parts of the gene have shown distinct evolutionary rates, which were associated with the diverse proproteins produced by posttranslational changes of DEC-1, with proproteins that are incorporated in the chorion earlier in egg formation being in general more conserved than others that are incorporated later. This correlation appears more evident in certain lineages, including the branch that separates Anastrepha, as well as other internal branches that differentiate species within the genus. Our data showed that this gene shows remarkable variation across its different exons, which has proven to be informative at different evolutionary levels. These changes hold promise not only for studying differentiation in Anastrepha but also for the eventual management of selected pest species.

Extensive longevity and DNA virus-driven adaptation in nearctic Myotis bats

The genus Myotis is one of the largest clades of bats, and exhibits some of the most extreme variation in lifespans among mammals alongside unique adaptations to viral tolerance and immune defense. To study the evolution of longevity-associated traits and infectious disease, we generated near-complete genome assemblies and cell lines for 8 closely related species of Myotis. Using genome-wide screens of positive selection, analyses of structural variation, and functional experiments in primary cell lines, we identify new patterns of adaptation contributing to longevity, cancer resistance, and viral interactions in bats. We find that Myotis bats have some of the most significant variation in cancer risk across mammals and demonstrate a unique DNA damage response in primary cells of the long-lived M. lucifugus. We also find evidence of abundant adaptation in response to DNA viruses - but not RNA viruses - in Myotis and other bats in sharp contrast with other mammals, potentially contributing to the role of bats as reservoirs of zoonoses. Together, our results demonstrate how genomics and primary cells derived from diverse taxa uncover the molecular bases of extreme adaptations in non-model organisms.

Phylogenomic Signatures of a Lineage of Vesicular Stomatitis Indiana Virus Circulating During the 2019-2020 Epidemic in the United States

For the first time, we describe phylogenomic signatures of an epidemic lineage of vesicular stomatitis Indiana virus (VSIV). We applied multiple evolutionary analyses to a dataset of 87 full-length genome sequences representing the circulation of an epidemic VSIV lineage in the US between 2019 and 2020. Based on phylogenetic analyses, we predicted the ancestral relationship of this lineage with a specific group of isolates circulating in the endemic zone of Chiapas, Mexico. Subsequently, our findings indicate that the lineage diversified into at least four different subpopulations during its circulation in the US. We identified single nucleotide polymorphisms (SNPs) that differentiate viral subpopulations and assessed their potential relevance using comparative phylogenetic methods, highlighting the preponderance of synonymous mutations during the differentiation of these populations. Purifying selection was the main evolutionary force favoring the conservation of this epidemic phenotype, with P and G genes as the main drivers of the evolution of this lineage. Our analyses identified multiple codon sites under positive selection and the association of these sites with specific functional domains at P, M, G, and L proteins. Based on ancestral reconstruction analyses, we showed the potential relevance of some of the sites identified under positive selection to the adaptation of the epidemic lineage at the population level. Finally, using a representative group of viruses from Colorado, we established a positive correlation between genetic and geographical distances, suggesting that positive selection on specific codon positions might have favored the adaptation of different subpopulations to circulation in specific geographical settings. Collectively, our study reveals the complex dynamics that accompany the evolution of an epidemic lineage of VSIV in nature. Our analytical framework provides a model for conducting future evolutionary analyses. The ultimate goal is to support the implementation of an early warning system for vesicular stomatitis virus in the US, enabling early detection of epidemic precursors from Mexico.

Molecular evolution of toothed whale genes reveals adaptations to echolocating in different environments

BACKGROUND

Echolocation was a key development in toothed whale evolution, enabling their adaptation and diversification across various environments. Previous bioacoustic and morphological studies suggest that environmental pressures have influenced the evolution of echolocation in toothed whales. This hypothesis demands further investigation, especially regarding the molecular mechanisms involved in the adaptive radiation of toothed whales across multiple habitats. Here we show that the coding sequences of four hearing genes involved in echolocation (CDH23, prestin, TMC1, and CLDN14) have different signatures of molecular evolution among riverine, coastal, and oceanic dolphins, suggesting that the evolutionary constraints of these habitats shaped the underlying genetic diversity of the toothed whale sonar.

RESULTS

Our comparative analysis across 37 odontocete species revealed patterns of accelerated evolution within coastal and riverine lineages, supporting the hypothesis that shallow habitats pose specific selective pressures to sonar propagation, which are not found in the deep ocean. All toothed whales with genes evolving under positive selection are shallow coastal species, including three species that have recently diverged from freshwater lineages (Cephalorhynchus commersonii, Sotalia guianensis, and Orcaella heinsohni - CDH23), and three species that operate specialized Narrow Band High Frequency (NBHF) Sonars (Phocoena sinus - prestin, Neophocaena phocaenoides - TMC1 and Cephalorhynchus commersonii - CDH23). For river dolphins and deep-diving toothed whales, we found signatures of positive selection and molecular convergence affecting specific sites on CDH23, TMC1, and prestin. Positively selected sites (PSS) were different in number, identity, and substitution rates (dN/dS) across riverine, coastal, and oceanic toothed whales.

CONCLUSION

Here we shed light on potential molecular mechanisms underlying the diversification of toothed whale echolocation. Our results suggest that toothed whale hearing genes changed under different selective pressures in coastal, riverine, and oceanic environments.

Accelerated differentiation of neo-W nuclear-encoded mitochondrial genes between two climate-associated bird lineages signals potential co-evolution with mitogenomes

There is considerable evidence for mitochondrial-nuclear co-adaptation as a key evolutionary driver. Hypotheses regarding the roles of sex-linkage have emphasized Z-linked nuclear genes with mitochondrial function (N-mt genes), whereas it remains contentious whether the perfect co-inheritance of W genes with mitogenomes could hinder or facilitate co-adaptation. Young (neo-) sex chromosomes that possess relatively many N-mt genes compared to older chromosomes provide unprecedented hypothesis-testing opportunities. Eastern Yellow Robin (EYR) lineages in coastal and inland habitats with different climates are diverged in mitogenomes, and in a ~ 15.4 Mb nuclear region enriched with N-mt genes, in contrast with otherwise-similar nuclear genomes. This nuclear region maps to passerine chromosome 1A, previously found to be neo-sex in the inland EYR genome. To compare sex-linked Chr1A-derived genes between lineages, we assembled and annotated the coastal EYR genome. We found that: (i) the coastal lineage shares a similar neo-sex system with the inland lineage, (ii) neo-W and neo-Z N-mt genes are not more diverged between lineages than are comparable non-N-mt genes, and showed little evidence for broad positive selection, (iii) however, W-linked N-mt genes are more diverged between lineages than are their Z-linked gametologs. The latter effect was ~7 times stronger for N-mt than non-N-mt genes, suggesting that W-linked N-mt genes might have diverged between lineages under environmental selection through co-evolution with mitogenomes. Finally, we identify a candidate gene driver for divergent selection, NDUFA12. Our data represent a rare example suggesting a possible role for W-associated mitochondrial-nuclear interactions in climate-associated adaptation and lineage differentiation.

Expansion of MHC-IIB Has Constrained the Evolution of MHC-IIA in Passerines

The major histocompatibility complex (MHC) is central in adaptive immunity, with the highly polymorphic MHC genes encoding antigen-presenting molecules. Two MHC class II (MHC-II) loci, DA1 and DA2, predate the radiation of extant birds and persist throughout much of the avian phylogeny. Within each locus, the MHC-II molecules are encoded by A-genes (DAA) and B-genes (DAB), which are arranged in A-B dyads. However, in passerines (order Passeriformes), the DA2 locus has been lost, and the ancestral A-B dyad at the DA1 locus has been replaced by a putatively single A-gene (DAA1) and an array of highly polymorphic B-genes (DAB1). In this study, we genotyped the DAA1 gene of 15 passerine species and confirmed that passerines possess just one copy of DAA1. We then compared selection patterns in DAA1 between passerines and nonpasserines and found that exon 2, which encodes the antigen-presenting domain, has been subject to weaker positive selection and stronger negative selection in passerines compared with nonpasserines. Additional comparisons showed that the patterns of selection in the passerine DAA1 gene are unlikely to be related to the loss of the DA2 locus. Instead, our findings suggest that the expansion of DAB1 (MHC-IIB) has imposed an evolutionary constraint on the passerine DAA1 (MHC-IIA) gene. We speculate that this constraint may be the result of each DAA1 chain forming heterodimers with many different DAB1 chains.

Bank vole genomics links determinate and indeterminate growth of teeth

BACKGROUND

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Bank and prairie voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars while retaining similar size and shape, providing alternative models for studying roots.

RESULTS

We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. Bulk transcriptomics comparisons of embryonic molar development between bank voles and mice demonstrated overall conservation of gene expression levels, with species-specific differences corresponding to the accelerated and more extensive patterning of the vole molar. We leverage convergent evolution of unrooted molars across the clade to examine changes that may underlie the evolution of unrooted molars. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation.

CONCLUSIONS

Our results support ongoing evolution of dental genes across Glires and identify candidate genes for mechanistic studies of root formation. Comparative research using the bank vole as a model species can reveal the complex evolutionary background of convergent evolution for ever-growing molars.

Revisiting the functional annotation of TriTryp using sequence similarity tools

Trypanosomatids are the causative agents of deadly diseases in humans and livestock. Given the high phylogenetic distance of trypanosomatids from model organisms, these organisms have ample unannotated genes. Manual functional annotation is time-consuming, highlighting the importance of automated functional annotation tools. The development of automated functional tools is a hot research topic, and multiple tools have been developed for the task. PANNZER2 is an automated functional annotation tool that merely relies on the sequence similarity of the query to the annotated proteins. We tried PANNZER2 on Trypanosoma brucei, the most studied organism among trypanosomatids, to see if it could improve our knowledge of the functions of the genes. Even with the availability of automated annotation tools like InterPro2GO in databases such as TriTrypDB, PANNZER2 has made surprisingly confident predictions for some hypothetical proteins in T. brucei. In this study, we identify gaps in such annotations because of not employing pairwise sequence alignment tools in TriTrypDB's automated annotation process. Our findings demonstrate that even the use of stringent cutoffs can successfully annotate a significant number of proteins. Additionally, we discovered that adjusting the open reading frames in certain genes leads to sequences with increased sequence signature coverage-characterized by the length covered by at least one sequence signature-compared to the original sequences. This enhanced sequence signature coverage suggests these genomic fragments could be pseudogenes. To facilitate further exploration, we developed a script to help identify potential pseudogenes within an organism's genome, offering researchers a new tool for genomic analysis and understanding. We extended all our analysis to Trypanosoma cruzi and Leishmania major to assess the impact of this approach across different species. Our study demonstrates that by utilizing pairwise sequence similarity alignment, even with stringent cutoffs, we can attribute 2986, 3953, and 3798 new GO terms to the genomes of T. brucei, T. cruzi, and L. major. Additionally, we found that 210, 239, and 29 genes exhibit increased sequence signature coverage following frame correction, suggesting the presence of pseudogenes.

Rethinking asexuality: the enigmatic case of functional sexual genes in Lepraria (Stereocaulaceae)

BACKGROUND

The ubiquity of sex across eukaryotes, given its high costs, strongly suggests it is evolutionarily advantageous. Asexual lineages can avoid, for example, the risks and energetic costs of recombination, but suffer short-term reductions in adaptive potential and long-term damage to genome integrity. Despite these costs, lichenized fungi have frequently evolved asexual reproduction, likely because it allows the retention of symbiotic algae across generations. The lichenized fungal genus Lepraria is thought to be exclusively asexual, while its sister genus Stereocaulon completes a sexual reproductive cycle. A comparison of sister sexual and asexual clades should shed light on the evolution of asexuality in lichens in general, as well as the apparent long-term maintenance of asexuality in Lepraria, specifically.

RESULTS

In this study, we assembled and annotated representative long-read genomes from the putatively asexual Lepraria genus and its sexual sister genus Stereocaulon, and added short-read assemblies from an additional 22 individuals across both genera. Comparative genomic analyses revealed that both genera were heterothallic, with intact mating-type loci of both idiomorphs present across each genus. Additionally, we identified and assessed 29 genes involved in meiosis and mitosis and 45 genes that contribute to formation of fungal sexual reproductive structures (ascomata). All genes were present and appeared functional in nearly all Lepraria, and we failed to identify a general pattern of relaxation of selection on these genes across the Lepraria lineage. Together, these results suggest that Lepraria may be capable of sexual reproduction, including mate recognition, meiosis, and production of ascomata.

CONCLUSIONS

Despite apparent maintenance of machinery essential for fungal sex, over 200 years of careful observations by lichenologists have produced no evidence of canonical sexual reproduction in Lepraria. We suggest that Lepraria may have instead evolved a form of parasexual reproduction, perhaps by repurposing MAT and meiosis-specific genes. This may, in turn, allow these lichenized fungi to avoid long-term consequences of asexuality, while maintaining the benefit of an unbroken bond with their algal symbionts.

Evidence of Intraspecific Adaptive Variation in the American Pika (Ochotona princeps) on a Continental Scale Using a Target Enrichment and Mitochondrial Genome Skimming Approach

Montane landscapes present an array of abiotic challenges that drive adaptive evolution amongst organisms. These adaptations can promote habitat specialisation, which may heighten the risk of extirpation from environmental change. For example, higher metabolic rates in an endothermic species may contribute to heightened cold tolerance, whilst simultaneously limiting heat tolerance. Here, using the climate-sensitive American pika (Ochotona princeps), we test for evidence of intraspecific adaptive variation amongst environmental gradients across the Intermountain West of North America. We leveraged results from previous studies on pika adaptation to generate a custom nuclear target enrichment design to sequence several hundred candidate genes related to cold, hypoxia and dietary detoxification. We also applied a 'genome skimming' approach to sequence mitochondrial DNA. Using genotype-environment association tests, we identified rare genomic variants associated with elevation and temperature variation amongst populations. Amongst mitochondrial genes, we identified intraspecific variation in selective signals and significant changes to the amino acid property equilibrium constant, which may relate to electron transport chain efficiency. These results illustrate a complex dynamic of adaptive variation amongst O. princeps where lineages and populations have adapted to unique regional conditions. Some of the clearest signals of selection were in a genetic lineage that includes pikas of the Great Basin region, which is also where recent localised extirpations have taken place and highlights the risk of losing adaptive alleles during environmental change.

Novel Sex-Specific Genes and Diverse Interspecific Expression in the Antennal Transcriptomes of Ithomiine Butterflies

The olfactory sense is crucial for organisms, facilitating environmental recognition and interindividual communication. Ithomiini butterflies exemplify this importance not only because they rely strongly on olfactory cues for both inter- and intra-sexual behaviors, but also because they show convergent evolution of specialized structures within the antennal lobe, called macroglomerular complexes (MGCs). These structures, widely absent in butterflies, are present in moths where they enable heightened sensitivity to, and integration of, information from various types of pheromones. In this study, we investigate chemosensory evolution across six Ithomiini species and identify possible links between expression profiles and neuroanatomical. To enable this, we sequenced four new high-quality genome assemblies and six sex-specific antennal transcriptomes for three of these species with different MGC morphologies. With extensive genomic analyses, we found that the expression of antennal transcriptomes across species exhibit profound divergence, and identified highly expressed ORs, which we hypothesize may be associated to MGCs, as highly expressed ORs are absent in Methona, an Ithomiini lineage which also lacks MGCs. More broadly, we show how antennal sexual dimorphism is prevalent in both chemosensory genes and non-chemosensory genes, with possible relevance for behavior. As an example, we show how lipid-related genes exhibit consistent sexual dimorphism, potentially linked to lipid transport or host selection. In this study, we investigate the antennal chemosensory adaptations, suggesting a link between genetic diversity, ecological specialization, and sensory perception with the convergent evolution of MCGs. Insights into chemosensory gene evolution, expression patterns, and potential functional implications enhance our knowledge of sensory adaptations and sexual dimorphisms in butterflies, laying the foundation for future investigations into the genetic drivers of insect behavior, adaptation, and speciation.