Structurally divergent and recurrently mutated regions of primate genomes.

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

Adaptive evolution of the enigmatic Takakia now facing climate change in Tibet.

The most extreme environments are the most vulnerable to transformation under a rapidly changing climate. These ecosystems harbor some of the most specialized species, which will likely suffer the highest extinction rates. We document the steepest temperature increase (2010-2021) on record at altitudes of above 4,000 m, triggering a decline of the relictual and highly adapted moss Takakia lepidozioides. Its de-novo-sequenced genome with 27,467 protein-coding genes includes distinct adaptations to abiotic stresses and comprises the largest number of fast-evolving genes under positive selection. The uplift of the study site in the last 65 million years has resulted in life-threatening UV-B radiation and drastically reduced temperatures, and we detected several of the molecular adaptations of Takakia to these environmental changes. Surprisingly, specific morphological features likely occurred earlier than 165 mya in much warmer environments. Following nearly 400 million years of evolution and resilience, this species is now facing extinction.

HGT is widespread in insects and contributes to male courtship in lepidopterans.

Horizontal gene transfer (HGT) is an important evolutionary force shaping prokaryotic and eukaryotic genomes. HGT-acquired genes have been sporadically reported in insects, a lineage containing >50% of animals. We systematically examined HGT in 218 high-quality genomes of diverse insects and found that they acquired 1,410 genes exhibiting diverse functions, including many not previously reported, via 741 distinct transfers from non-metazoan donors. Lepidopterans had the highest average number of HGT-acquired genes. HGT-acquired genes containing introns exhibited substantially higher expression levels than genes lacking introns, suggesting that intron gains were likely involved in HGT adaptation. Lastly, we used the CRISPR-Cas9 system to edit the prevalent unreported gene LOC105383139, which was transferred into the last common ancestor of moths and butterflies. In diamondback moths, males lacking LOC105383139 courted females significantly less. We conclude that HGT has been a major contributor to insect adaptation.

Genetic mechanisms of animal camouflage: an interdisciplinary perspective.

Camouflage is a classic example of a trait wherein animals respond to natural selection to avoid predation or attract prey. This unique phenomenon has attracted significant recent attention and the rapid development of integrative research methods is facilitating advances in our understanding of the in-depth genetic mechanisms of camouflage. In this review article, we revisit camouflage definitions and strategies and then we examine the underlying mechanisms of the two most common forms of camouflage, crypsis and masquerade, that have recently been elucidated using multiple approaches. We also discuss unresolved questions related to camouflage. Ultimately, we highlight the implications of camouflage for informing various key issues in ecology and evolution.

Finding genes and pathways that underlie coral adaptation.

Mass coral bleaching is one of the clearest threats of climate change to the persistence of marine biodiversity. Despite the negative impacts of bleaching on coral health and survival, some corals may be able to rapidly adapt to warming ocean temperatures. Thus, a significant focus in coral research is identifying the genes and pathways underlying coral heat adaptation. Here, we review state-of-the-art methods that may enable the discovery of heat-adaptive loci in corals and identify four main knowledge gaps. To fill these gaps, we describe an experimental approach combining seascape genomics with CRISPR/Cas9 gene editing to discover and validate heat-adaptive loci. Finally, we discuss how information on adaptive genotypes could be used in coral reef conservation and management strategies.

Widespread adaptive evolution in angiosperm photosystems provides insight into the evolution of photosystem II repair.

Oxygenic photosynthesis generates the initial energy source that fuels nearly all life on Earth. At the heart of the process are the photosystems, which are pigment binding multi-protein complexes that catalyse the first step of photochemical conversion of light energy into chemical energy. Here, we investigate the molecular evolution of the plastid-encoded photosystem subunits at single-residue resolution across 773 angiosperm species. We show that despite an extremely high level of conservation, 7% of residues in the photosystems, spanning all photosystem subunits, exhibit hallmarks of adaptive evolution. Through in silico modelling of these adaptive substitutions, we uncover the impact of these changes on the predicted properties of the photosystems, focussing on their effects on co-factor binding and inter-subunit interface formation. By analyzing these cohorts of changes, we reveal that evolution has repeatedly altered the interaction between photosystem II and its D1 subunit in a manner that is predicted to reduce the energetic barrier for D1 turnover and photosystem repair. Together, these results provide insight into the trajectory of photosystem adaptation during angiosperm evolution.

Pleiotropy of positive selection in ancient ACE2 suggests an alternative hypothesis for bat-specific adaptations to host coronaviruses.

Angiotensin-convertingenzyme 2 (ACE2) has dual functions, regulating cardiovascular physiology and serving as the receptor for coronaviruses. Bats, the only true flying mammals and natural viral reservoirs, have evolved positive alterations in traits related to both functions of ACE2. This suggests significant evolutionary changes in ACE2 during bat evolution. To test this hypothesis, we examine the selection pressure in ACE2 along the ancestral branch of all bats (AncBat-ACE2), where powered flight and bat-coronavirus coevolution occurred, and detect a positive selection signature. To assess the functional effects of positive selection, we resurrect AncBat-ACE2 and its mutant (AncBat-ACE2-mut) created by replacing the positively selected sites. Compared to AncBat-ACE2-mut, AncBat-ACE2 exhibits stronger enzymatic activity, enhances mice's performance in exercise fatigue, and shows lower affinity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our findings indicate the functional pleiotropy of positive selection in the ancient ACE2 of bats, providing an alternative hypothesis for the evolutionary origin of bats' defense against coronaviruses.

GPCR-MAPK signaling pathways underpin fitness trade-offs in whitefly.

Trade-offs between evolutionary gain and loss are prevalent in nature, yet their genetic basis is not well resolved. The evolution of insect resistance to insecticide is often associated with strong fitness costs; however, how the fitness trade-offs operates remains poorly understood. Here, we show that the mitogen-activated protein kinase (MAPK) pathway and its upstream and downstream actors underlie the fitness trade-offs associated with insecticide resistance in the whitefly Bemisia tabaci. Specifically, we find a key cytochrome P450 gene CYP6CM1, that confers neonicotinoids resistance to in B. tabaci, is regulated by the MAPKs p38 and ERK through their activation of the transcription factor cAMP-response element binding protein. However, phosphorylation of p38 and ERK also leads to the activation of the transcription repressor Cap "n" collar isoform C (CncC) that negatively regulates exuperantia (Ex), vasa (Va), and benign gonial cell neoplasm (Bg), key genes involved in oogenesis, leading to abnormal ovary growth and a reduction in female fecundity. We further demonstrate that the transmembrane G protein-coupled receptor (GPCR) neuropeptide FF receptor 2 (NPFF2) triggers the p38 and ERK pathways via phosphorylation. Additionally, a positive feedback loop between p38 and NPFF2 leads to the continuous activation of the MAPK pathways, thereby constitutively promoting neonicotinoids resistance but with a significant reproductive cost. Collectively, these findings provide fundamental insights into the role of cis-trans regulatory networks incurred by GPCR-MAPK signaling pathways in evolutionary trade-offs and applied knowledge that can inform the development of strategies for the sustainable pest control.

Immune cells use active tugging forces to distinguish affinity and accelerate evolution.

Cells are known to exert forces to sense their physical surroundings for guidance of motion and fate decisions. Here, we propose that cells might do mechanical work to drive their own evolution, taking inspiration from the adaptive immune system. Growing evidence indicates that immune B cells-capable of rapid Darwinian evolution-use cytoskeletal forces to actively extract antigens from other cells' surfaces. To elucidate the evolutionary significance of force usage, we develop a theory of tug-of-war antigen extraction that maps receptor binding characteristics to clonal reproductive fitness, revealing physical determinants of selection strength. This framework unifies mechanosensing and affinity-discrimination capabilities of evolving cells: Pulling against stiff antigen tethers enhances discrimination stringency at the expense of absolute extraction. As a consequence, active force usage can accelerate adaptation but may also cause extinction of cell populations, resulting in an optimal range of pulling strength that matches molecular rupture forces observed in cells. Our work suggests that nonequilibrium, physical extraction of environmental signals can make biological systems more evolvable at a moderate energy cost.

Hemoglobin wonders: a fascinating gas transporter dive into molluscs.

Hemoglobin (Hb) has been identified in at least 14 molluscan taxa so far. Research spanning over 130 years on molluscan Hbs focuses on their genes, protein structures, functions, and evolution. Molluscan Hbs are categorized into single-, two-, and multiple-domain chains, including red blood cell, gill, and extracellular Hbs, based on the number of globin domains and their respective locations. These Hbs exhibit variation in assembly, ranging from monomeric and dimeric to higher-order multimeric forms. Typically, molluscan Hbs display moderately high oxygen affinity, weak cooperativity, and varying pH sensitivity. Hb's potential role in antimicrobial pathways could augment the immune defense of bivalves, which may be a complement to their lack of adaptive immunity. The role of Hb as a respiratory protein in bivalves likely originated from the substitution of hemocyanin. Molluscan Hbs demonstrate adaptive evolution in response to environmental changes via various strategies (e.g. increasing Hb types, multimerization, and amino acid residue substitutions at key sites), enhancing or altering functional properties for habitat adaptation. Concurrently, an increase in Hb assembly diversity, coupled with a downward trend in oxygen affinity, is observed during molluscan differentiation and evolution. Hb in Protobranchia, Heteroconchia, and Pteriomorphia bivalves originated from separate ancestors, with Protobranchia inheriting a relative ancient molluscan Hb gene. In bivalves, extracellular Hbs share a common origin, while gill Hbs likely emerged from convergent evolution. In summary, research on molluscan Hbs offers valuable insights into the origins, biological variations, and adaptive evolution of animal Hbs.

Aneuploidy Can Be an Evolutionary Diversion on the Path to Adaptation.

Aneuploidy is common in eukaryotes, often leading to decreased fitness. However, evidence from fungi and human tumur cells suggests that specific aneuploidies can be beneficial under stressful conditions and facilitate adaptation. In a previous evolutionary experiment with yeast, populations evolving under heat stress became aneuploid, only to later revert to euploidy after beneficial mutations accumulated. It was therefore suggested that aneuploidy is a "stepping stone" on the path to adaptation. Here, we test this hypothesis. We use Bayesian inference to fit an evolutionary model with both aneuploidy and mutation to the experimental results. We then predict the genotype frequency dynamics during the experiment, demonstrating that most of the evolved euploid population likely did not descend from aneuploid cells, but rather from the euploid wild-type population. Our model shows how the beneficial mutation supply-the product of population size and beneficial mutation rate-determines the evolutionary dynamics: with low supply, much of the evolved population descends from aneuploid cells; but with high supply, beneficial mutations are generated fast enough to outcompete aneuploidy due to its inherent fitness cost. Our results suggest that despite its potential fitness benefits under stress, aneuploidy can be an evolutionary "diversion" rather than a "stepping stone": it can delay, rather than facilitate, the adaptation of the population, and cells that become aneuploid may leave less descendants compared to cells that remain diploid.

Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus.

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

HaploSweep: Detecting and Distinguishing Recent Soft and Hard Selective Sweeps through Haplotype Structure.

Identifying soft selective sweeps using genomic data is a challenging yet crucial task in population genetics. In this study, we present HaploSweep, a novel method for detecting and categorizing soft and hard selective sweeps based on haplotype structure. Through simulations spanning a broad range of selection intensities, softness levels, and demographic histories, we demonstrate that HaploSweep outperforms iHS, nSL, and H12 in detecting soft sweeps. HaploSweep achieves high classification accuracy-0.9247 for CHB, 0.9484 for CEU, and 0.9829 YRI-when applied to simulations in line with the human Out-of-Africa demographic model. We also observe that the classification accuracy remains consistently robust across different demographic models. Additionally, we introduce a refined method to accurately distinguish soft shoulders adjacent to hard sweeps from soft sweeps. Application of HaploSweep to genomic data of CHB, CEU, and YRI populations from the 1000 genomes project has led to the discovery of several new genes that bear strong evidence of population-specific soft sweeps (HRNR, AMBRA1, CBFA2T2, DYNC2H1, and RANBP2 etc.), with prevalent associations to immune functions and metabolic processes. The validated performance of HaploSweep, demonstrated through both simulated and real data, underscores its potential as a valuable tool for detecting and comprehending the role of soft sweeps in adaptive evolution.

Adaptive Evolution of Pseudomonas aeruginosa in Human Airways Shows Phenotypic Convergence Despite Diverse Patterns of Genomic Changes.

Selective forces in the environment drive bacterial adaptation to novel niches, choosing the fitter variants in the population. However, in dynamic and changing environments, the evolutionary processes controlling bacterial adaptation are difficult to monitor. Here, we follow 9 people with cystic fibrosis chronically infected with Pseudomonas aeruginosa, as a proxy for bacterial adaptation. We identify and describe the bacterial changes and evolution occurring between 15 and 35 yr of within-host evolution. We combine whole-genome sequencing, RNA sequencing, and metabolomics and compare the evolutionary trajectories directed by the adaptation of 4 different P. aeruginosa lineages to the lung. Our data suggest divergent evolution at the genomic level for most of the genes, with signs of convergent evolution with respect to the acquisition of mutations in regulatory genes, which drive the transcriptional and metabolomic program at late time of evolution. Metabolomics further confirmed convergent adaptive phenotypic evolution as documented by the reduction of the quorum-sensing molecules acyl-homoserine lactone, phenazines, and rhamnolipids (except for quinolones). The modulation of the quorum-sensing repertoire suggests that similar selective forces characterize at late times of evolution independent of the patient. Collectively, our data suggest that similar environments and similar P. aeruginosa populations in the patients at prolonged time of infection are associated with an overall reduction of virulence-associated features and phenotypic convergence.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory.

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Genome structural variation in human evolution.

Structural variation (SV) is a large difference (typically >100 bp) in the genomic structure of two genomes and includes both copy number variation and variation that does not change copy number of a genomic region, such as an inversion. Improved reference genomes, combined with widespread genome sequencing using short-read sequencing technology, and increasingly using long-read sequencing, have reignited interest in SV. Recent large-scale studies and functional focused analyses have highlighted the role of SV in human evolution. In this review, we highlight human-specific SVs involved in changes in the brain, population-specific SVs that affect response to the environment, including adaptation to diet and infectious diseases, and summarise the contribution of archaic hominin admixture to present-day human SV.

Envelope domain III E324, E351, and E380 mutations lever adaptive evolution of DENV-1 genotype I.

Dengue virus (DENV) gains genetic mutations during continuous transmission and evolution, making the virus more adaptive and virulent. The clade of DENV-1 genotype I has expanded and become the predominant genotype in Asia and the Pacific areas, but the underlying mechanisms are unclear. A combined analysis of nonsynonymous mutations in domain III of the envelope protein and their biological effects on virus pathogenesis and transmission was evaluated. Phylogenetic analyses found three nonsynonymous mutations (V324I, V351L, and V380I) in domain III of the envelope protein, which emerged in 1970s-1990s and stably inherited and expanded in contemporary strains after 2000. We generated reverse-mutated viruses (I324V, L351V, and I380V) based on an infectious clone of an epidemic DENV-1 strain (NIID02-20), and the results suggested that the infectivity of the contemporary epidemic virus (wild type, WT) has increased compared to the reverse mutant viruses in mammalian hosts but not mosquito vectors. The WT virus showed a higher binding affinity to host cells and increased virion stability. In addition, weaker immunogenicity and higher resistance to neutralizing antibodies of the WT virus indicated a trend of immune escape. The data suggested that nonsynonymous mutations of the E protein (V324I, V351L, and V380I) promote infectivity and immune evasion of DENV-1 genotype I, which may facilitate its onward transmission on a global scale. IMPORTANCE: We provide evidence that minor sequence variation among dengue virus (DENV) strains can result in increased adaptability and virulence, impacting both the biology of the virus and the antiviral immune response. The genetic mutations of DENV-1 gained during continuous transmission and evolution will offer new clues for the design of novel vaccines against flaviviruses.

Origin and evolution of the blue light receptor cryptochromes (CRY1/2) in aquatic angiosperms.

Cryptochromes (CRYs), which are responsible for sensing blue light in plants, play a critical role in regulating blue light signals and circadian rhythms. However, their functions extend beyond light detection, as they also aid plants in adapting to stress and potentially other regulatory mechanisms. Aquatic angiosperms, which independently evolved from various angiosperm lineages, have developed specific adaptations to unique light qualities and environmental stressors found in aquatic habitats compared to terrestrial ones. It was hypothesized that the sequences and regulatory networks of angiosperm CRY1/2 underwent adaptive evolution in different aquatic angiosperm lineages. To test this hypothesis, we compiled comprehensive datasets consisting of 55 green plant genomes (including 37 angiosperm genomes), 80 angiosperm transcriptomes, and 4 angiosperm expression networks. Through comparative analysis, we found that CRY1 originated from a common ancestor of seed plants, whereas CRY2 originated from a common ancestor of land plants. In angiosperms, the CRY1/2 sequences of aquatic lineages exhibited positive selection, and the conserved valine-proline (VP) motif of CRY2 showed a convergent loss in two aquatic species. Co-expressed genes associated with blue light receptors (CRY) showed adaptations to aquatic environments, specifically in relation to flooding and osmotic stress. These discoveries shed light on the adaptive evolution of CRY1/2, encompassing their origins, sequences, and regulatory networks. Furthermore, these results provide valuable insights for investigating the uncharacterized functions and regulatory pathways of CRY and offer potential targets for enhancing growth and adaptation in agricultural plants.

Functional genomics analysis reveals the evolutionary adaptation and demographic history of pygmy lorises.

Lorises are a group of globally threatened strepsirrhine primates that exhibit many unusual physiological and behavioral features, including a low metabolic rate, slow movement, and hibernation. Here, we assembled a chromosome-level genome sequence of the pygmy loris (Xanthonycticebus pygmaeus) and resequenced whole genomes from 50 pygmy lorises and 6 Bengal slow lorises (Nycticebus bengalensis). We found that many gene families involved in detoxification have been specifically expanded in the pygmy loris, including the GSTA gene family, with many newly derived copies functioning specifically in the liver. We detected many genes displaying evolutionary convergence between pygmy loris and koala, including PITRM1. Significant decreases in PITRM1 enzymatic activity in these two species may have contributed to their characteristic low rate of metabolism. We also detected many evolutionarily convergent genes and positively selected genes in the pygmy loris that are involved in muscle development. Functional assays demonstrated the decreased ability of one positively selected gene, MYOF, to up-regulate the fast-type muscle fiber, consistent with the lower proportion of fast-twitch muscle fibers in the pygmy loris. The protein product of another positively selected gene in the pygmy loris, PER2, exhibited weaker binding to the key circadian core protein CRY, a finding that may be related to this species' unusual circadian rhythm. Finally, population genomics analysis revealed that these two extant loris species, which coexist in the same habitat, have exhibited an inverse relationship in terms of their demography over the past 1 million years, implying strong interspecies competition after speciation.

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)6, Ge] and Gossypium stephensii [(AD)7, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)1, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates-including phenotypic differentiation, genetic isolation, and genetic convergence-that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.