Adaptive Variation Research Articles

A distributed adaptive policy gradient method based on momentum for multi-agent reinforcement learning

Policy Gradient (PG) method is one of the most popular algorithms in Reinforcement Learning (RL). However, distributed adaptive variants of PG are rarely studied in multi-agent. For this reason, this paper proposes a distributed adaptive policy gradient algorithm (IS-DAPGM) incorporated with Adam-type updates and importance sampling technique. Furthermore, we also establish the theoretical convergence rate of O(1/T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {O}(1/\\sqrt{T})$$\\end{document}, where T represents the number of iterations, it can match the convergence rate of the state-of-the-art centralized policy gradient methods. In addition, many experiments are conducted in a multi-agent environment, which is a modification on the basis of Particle world environment. By comparing with some other distributed PG methods and changing the number of agents, we verify the performance of IS-DAPGM is more efficient than the existing methods.

Multilevel gene expression changes in lineages containing adaptive copy number variants.

Copy-number variants (CNVs) are an important class of recurrent variants that mediate adaptive evolution. While CNVs can increase the relative fitness of the organism, they can also incur a cost. We previously evolved populations of Saccharomyces cerevisiae over hundreds of generations in glutamine-limited (Gln-) chemostats and observed the recurrent evolution of CNVs at the GAP1 locus. To understand the role that expression plays in adaptation, both in relation to the adaptation of the organism to the selective condition, and as a consequence of the CNV, we measured the transcriptome, translatome, and proteome of 4 strains of evolved yeast, each with a unique CNV, and their ancestor in Gln- conditions. We find CNV-amplified genes correlate with higher RNA abundance; however, this effect is reduced at the level of the proteome, consistent with post-transcriptional dosage compensation. By normalizing each level of expression by the abundance of the preceding step we were able to identify widespread divergence in the efficiency of each step in the gene in the efficiency of each step in gene expression. Genes with significantly different translational efficiency were enriched for potential regulatory mechanisms including either upstream open reading frames, RNA binding sites for SSD1, or both. Genes with lower protein expression efficiency were enriched for genes encoding proteins in protein complexes. Taken together, our study reveals widespread changes in gene expression at multiple regulatory levels in lineages containing adaptive CNVs highlighting the diverse ways in which adaptive evolution shapes gene expression.

Ecological stress memory in wood architecture of two Neotropical hickory species from central-eastern Mexico

BackgroundDrought periods are major evolutionary triggers of wood anatomical adaptive variation in Lower Tropical Montane Cloud Forests tree species. We tested the influence of historical drought events on the effects of ecological stress memory on latewood width and xylem vessel traits in two relict hickory species (Carya palmeri and Carya myristiciformis) from central-eastern Mexico. We hypothesized that latewood width would decrease during historical drought years, establishing correlations between growth and water stress conditions, and that moisture deficit during past tree growth between successive drought events, would impact on wood anatomical features. We analyzed latewood anatomical traits that developed during historical drought and pre- and post-drought years in both species.ResultsWe found that repeated periods of hydric stress left climatic signatures for annual latewood growth and xylem vessel traits that are essential for hydric adaptation in tropical montane hickory species.ConclusionsOur results demonstrate the existence of cause‒effect relationships in wood anatomical architecture and highlight the ecological stress memory linked with historical drought events. Thus, combined time-series analysis of latewood width and xylem vessel traits is a powerful tool for understanding the ecological behavior of hickory species.

Adaptive orthogonal directional total variation with kernel regression for CT image denoising.

Low-dose computed tomography (CT) has been successful in reducing radiation exposure for patients. However, the use of reconstructions from sparse angle sampling in low-dose CT often leads to severe streak artifacts in the reconstructed images. In order to address this issue and preserve image edge details, this study proposes an adaptive orthogonal directional total variation method with kernel regression. The CT reconstructed images are initially processed through kernel regression to obtain the N-term Taylor series, which serves as a local representation of the regression function. By expanding the series to the second order, we obtain the desired estimate of the regression function and localized information on the first and second derivatives. To mitigate the noise impact on these derivatives, kernel regression is performed again to update the first and second derivatives. Subsequently, the original reconstructed image, its local approximation, and the updated derivatives are summed using a weighting scheme to derive the image used for calculating orientation information. For further removal of stripe artifacts, the study introduces the adaptive orthogonal directional total variation (AODTV) method, which denoises along both the edge direction and the normal direction, guided by the previously obtained orientation. Both simulation and real experiments have obtained good results. The results of two real experiments show that the proposed method has obtained PSNR values of 34.5408 dB and 29.4634 dB, which are 1.2392-5.9333 dB and 2.828-6.7995 dB higher than the contrast denoising algorithm, respectively, indicating that the proposed method has good denoising performance. The study demonstrates the effectiveness of the method in eliminating strip artifacts and preserving the fine details of the images.

Epigenetic diversity of genes with copy number variations among natural populations of the three-spined stickleback.

Duplicated genes provide the opportunity for evolutionary novelty and adaptive divergence. In many cases, having more gene copies increases gene expression, which might facilitate adaptation to stressful or novel environments. Conversely, overexpression or misexpression of duplicated genes can be detrimental and subject to negative selection. In this scenario, newly duplicate genes may evade purifying selection if they are epigenetically silenced, at least temporarily, leading them to persist in populations as copy number variations (CNVs). In animals and plants, younger gene duplicates tend to have higher levels of DNA methylation and lower levels of gene expression, suggesting epigenetic regulation could promote the retention of gene duplications via expression repression or silencing. Here, we test the hypothesis that DNA methylation variation coincides with young duplicate genes that are segregating as CNVs in six populations of the three-spined stickleback that span a salinity gradient from 4 to 30 PSU. Using reduced-representation bisulfite sequencing, we found DNA methylation and CNV differentiation outliers rarely overlapped. Whereas lineage-specific genes and young duplicates were found to be highly methylated, just two gene CNVs showed a significant association between promoter methylation level and copy number, suggesting that DNA methylation might not interact with CNVs in our dataset. If most new duplications are regulated for dosage by epigenetic mechanisms, our results do not support a strong contribution from DNA methylation soon after duplication. Instead, our results are consistent with a preference to duplicate genes that are already highly methylated.

High-intensity light promotes adaptive divergence of photosynthetic traits between sun-exposed and shaded populations in Saxifraga fortunei.

Light is essential for plants, and local populations exhibit adaptive photosynthetic traits depending on their habitats. Although plastic responses in morphological and/or physiological characteristics to different light intensities are well known, adaptive divergence with genetic variation remains to be explored. This study focused on Saxifraga fortunei (Saxifragaceae) growing in sun-exposed and shaded habitats. We measured the leaf anatomical structure and photosynthetic rate of plants grown in their natural habitats and in a common greenhouse (high- and low-intensity light experimental sites). To assess differences in ecophysiological tolerance to high-intensity light between the sun and shade types, we evaluated the level of photoinhibition of photosystem II and the leaf mortality rate under high-intensity light conditions. In addition, population genetic analysis was conducted to investigate phylogenetic origins. Clear phenotypic differences were found between the sun and shade types despite their recent phylogenetic origin. The leaf anatomical structure and photosynthetic rate showed plastic changes in response to growing conditions. Moreover, the sun type had a well-developed palisade parenchyma and a higher photosynthetic rate, which were genetically fixed, and a lower level of photoinhibition under high-intensity light. Our findings demonstrate that light intensity is a selective pressure that can rapidly promote phenotypic divergence between the sun and shade types. While phenotypic changes in multiple photosynthetic traits were plastic, genetic divergence in specific traits related to adaptation to high-intensity light would be fundamental for ecotypic divergence to different light regimes.

Adaptive infrared radiation modulator using reversible metal electrodeposition with silicon carbide-based fabry-perot cavity

Conventional materials that statically modulate radiation in the infrared (IR) region of the electromagnetic spectrum underpin the performance of many entrenched technologies. The development of their adaptive variants, whereby the IR radiation properties dynamically change as a response to external stimuli, has emerged as an important unfulfilled scientific challenge. Here, by combining the silicon carbide (SiC)-based Fabry-Perot (FP) cavity effect with the traditional reversible metal electrodeposition device (RMED), where the FP cavity is used as a color filter to emit a particular wavelength and the RMED is used to reversibly modulate IR radiation, we develop an adaptive IR radiation modulator with large and uniform IR emissivity tunability (Δε ≥ 0.579) in the 3–14 µm region. A specific structural color which is a function of the SiC thickness can be simultaneously displayed when the modulator is in the low-emissivity mode. By varying the SiC thickness from 100 to 180 nm, a full gamut of colors spanning the entire visible spectrum can be achieved. Additionally, the modulator can exhibit fantastic patterns consisting of different structural colors with simple designs achieved. We believe that our findings may open opportunities for camouflage in multispectral detection and many technologies related to IR radiation management.

Molecular evolution of the mammalian kinetochore complex.

Mammalian centromeres are satellite-rich chromatin domains that serve as sites for kinetochore complex assembly. Centromeres are highly variable in sequence and satellite organization across species, but the processes that govern the co-evolutionary dynamics between rapidly evolving centromeres and their associated kinetochore proteins remain poorly understood. Here, we pursue a course of phylogenetic analyses to investigate the molecular evolution of the complete kinetochore complex across primate and rodent species with divergent centromere repeat sequences and features. We show that many protein components of the core centromere associated network (CCAN) harbor signals of adaptive evolution, consistent with their intimate association with centromere satellite DNA and roles in the stability and recruitment of additional kinetochore proteins. Surprisingly, CCAN and outer kinetochore proteins exhibit comparable rates of adaptive divergence, suggesting that changes in centromere DNA can ripple across the kinetochore to drive adaptive protein evolution within distant domains of the complex. Our work further identifies kinetochore proteins subject to lineage-specific adaptive evolution, including rapidly evolving proteins in species with centromere satellites characterized by higher-order repeat structure and lacking CENP-B boxes. Thus, features of centromeric chromatin beyond the linear DNA sequence may drive selection on kinetochore proteins. Overall, our work spotlights adaptively evolving proteins with diverse centromere-associated functions, including centromere chromatin structure, kinetochore protein assembly, kinetochore-microtubule association, cohesion maintenance, and DNA damage response pathways. These adaptively evolving kinetochore protein candidates present compelling opportunities for future functional investigations exploring how their concerted changes with centromere DNA ensure the maintenance of genome stability.

Rolling down that mountain: microgeographical adaptive divergence during a fast population expansion along a steep environmental gradient in European beech

Forest tree populations harbour high genetic diversity thanks to large effective population sizes and strong gene flow, allowing them to diversify through adaptation to local environmental pressures within dispersal distance. Many tree populations also experienced historical demographic fluctuations, including spatial population contraction or expansions at various temporal scales, which may constrain their ability to adapt to environmental variations. Our aim is to investigate how recent contraction and expansion events interfere with local adaptation, by studying patterns of adaptive divergence between closely related stands undergoing environmentally contrasted conditions, and having or not recently expanded. To investigate genome-wide signatures of local adaptation while accounting for demography, we analysed divergence in a European beech population by testing pairwise differentiation among four tree stands at ~35k Single Nucleotide Polymorphisms from ~9k genomic regions. We applied three divergence outlier search methods resting on different assumptions and targeting either single SNPs or contiguous genomic regions, while accounting for the effect of population size variations on genetic divergence. We found 27 signals of selective signatures in 19 target regions. Putatively adaptive divergence involved all stand pairs. We retrieved signals both when comparing old-growth stands and recently colonised areas and when comparing stands within the old-growth area. Therefore, adaptive divergence processes have taken place both over short time spans, under strong environmental contrasts, and over short ecological gradients, in populations that have been stable in the long term. This suggests that standing genetic variation supports local, microgeographic divergence processes, which can maintain genetic diversity at the landscape level.

UV light and adaptive divergence of leaf physiology, anatomy, and ultrastructure drive heat stress tolerance in genetically distant grapevines.

The genetic basis of plant response to light and heat stresses had been unveiled, and different molecular mechanisms of leaf cell homeostasis to keep high physiological performances were recognized in grapevine varieties. However, the ability to develop heat stress tolerance strategies must be further elucidated since the morpho-anatomical and physiological traits involved may vary with genotype × environment combination, stress intensity, and duration. A 3-year experiment was conducted on potted plants of Sardinian red grapevine cultivars Cannonau (syn. Grenache) and Carignano (syn. Carignan), exposed to prolonged heat stress inside a UV-blocking greenhouse, either submitted to low daily UV-B doses of 4.63 kJ m-2 d-1 (+UV) or to 0 kJ m-2 d-1 (-UV), and compared to a control (C) exposed to solar radiation (4.05 kJ m-2 d-1 average UV-B dose). Irrigation was supplied to avoid water stress, and canopy light and thermal microclimate were monitored continuously. Heat stress exceeded one-third of the duration inside the greenhouse and 6% in C. In vivo spectroscopy, including leaf reflectance and fluorescence, allowed for characterizing different patterns of leaf traits and metabolites involved in oxidative stress protection. Cannonau showed lower stomatal conductance under C (200 mmol m-2 s-1) but more than twice the values inside the greenhouse (400 to 900 mmol m-2 s-1), where water use efficiency was reduced similarly in both varieties. Under severe heat stress and -UV, Cannonau showed a sharper decrease in primary photochemical activity and higher leaf pigment reflectance indexes and leaf mass area. UV-B increased the leaf pigments, especially in Carignano, and different leaf cell regulatory traits to prevent oxidative damage were observed in leaf cross-sections. Heat stress induced chloroplast swelling, plastoglobule diffusion, and the accumulation of secretion deposits in both varieties, aggravated in Cannonau -UV by cell vacuolation, membrane dilation, and diffused leaf blade spot swelling. Conversely, in Carignano UV-B, cell wall barriers and calcium oxalate crystals proliferated in mesophyll cells. These responses suggest an adaptive divergence among cultivars to prolonged heat stress and UV-B light. Further research on grapevine biodiversity, heat, and UV-B light interactions may give new insights on the extent of stress tolerance to improve viticulture adaptive strategies in climate change hotspots.

The Impact of Biologically Adaptive Sex Ratio Variation in Lamprey Population on Ecosystem Stability: Based on Lotka-Volterra Model

The Impact of Biologically Adaptive Sex Ratio Variation in Lamprey Population on Ecosystem Stability: Based on Lotka-Volterra Model

Ecological transcriptomics reveals stress response pathways of a ground-herb species in a waterlogging gradient of Amazonian riparian forests.

Environmental stress is a fundamental facet of life and a significant driver of natural selection in the wild. Gene expression diversity may facilitate adaptation to environmental changes, without necessary genetic change, but its role in adaptive divergence remains largely understudied in Neotropical systems. In Amazonian riparian forests, species distribution is predominantly influenced by species' waterlogging tolerance. The flooding gradient delineates distinct wetland forest types, shaping habitats and species characteristics. Here we investigated the molecular basis of environmental stress response in a tropical ground-herb species (Ischnosiphon puberulus) to environmental variation in Amazonian riparian forests. We compared environmental variables and gene expression profiles from individuals collected in two forest types: Igapó and Terra firme in the Amazonian riparian forests. Predictable seasonal flooding poses a significant challenge in Igapó compared to Terra firme environments, with the former presenting higher water column height and longer flooding duration. Our findings suggest that contrasting environmental conditions related to flooding regimes are important drivers of population genetic differentiation and differential gene expression in I. puberulus. Enriched gene ontology terms highlight associations with environmental stresses, such as defence response, water transport, phosphorylation, root development, response to auxin, salicylic acid and oxidative stress. By uncovering key environmental stress response pathways conserved across populations, I. puberulus offers novel genetic insights into the molecular basis of plant reactions to environmental constraints found in flooded areas of this highly biodiverse neotropical ecosystem.

Introgression of the gamete recognition molecule, bindin, in the sea urchin Diadema.

Hybridization is important in evolution, because it is a necessary (though not sufficient) step in the introgression of potentially adaptive variation between species. Bindin is a gamete recognition protein in echinoids and asteroids, capable of blocking cross-fertilization between species to varying degrees. Four species of the sea urchin genus Diadema are broadly sympatric in the Indo-Pacific: D. paucispinum, D. savignyi, D. clarki, and D. setosum. Data from three published studies, one of identification of hybrids through allozymes, one of the phylogeography of mitochondrial DNA, and one of the phylogeny of bindin, were combined to assess the degree of bindin introgression between these four species. I analyzed sequences of the ATPase 8 and ATPase 6 mitochondrial genes and of bindin, sampled throughout the species ranges, with an Isolation-Migration algorithm, IMa3. IMa3 uses a coalescent approach to produce Bayesian estimates of effective population sizes and gene flow between populations. The results showed that bindin alleles coalesce completely within the species bounds of D. clarki, and of D. setosum. The sister species D. paucispinum and D. savignyi, however, were estimated as having exchanged a bindin allele at an average of every one to two-and-a-half generations since they speciated from each other. As the allozyme study detected nine hybrids between three of these species in Okinawa (most of them between D. setosum and D. savignyi) in a single sample, hybrids between these species are produced, but bindin does not introgress. Therefore, bindin must not be efficient in blocking heterospecific fertilizations. Complete, or almost complete, reproductive isolation between species of Diadema must result from low hybrid fitness.

Adaptive convergence and divergence underpin the diversity of Asteraceae in a semi-arid lowland region

Asteraceae is the most diverse family primarily found in semi-arid areas and species from the same botanical family have comparable phenotypes, occupy semi-arid habitats, and have adaptive convergences to environmental stress. Such characteristics aid in understanding adaptive radiations of recent groups on a regional scale and their successful establishment in habitats with challenging climatic and edaphic conditions. We investigated the adaptive convergence of the family as well as the fragile balance between fine-scale variation in adaptive traits and general convergence towards stress-tolerance and adaptation to aridity, using thorough field surveys and trait-based ecological approaches. Data was collected between 2020 and 2022 using the quadrat approach encompassed across 60 sites and five habitat types. Multivariate statistical methodologies were applied to analyze plant traits and soil variables. The tribes Cichorieae, Cardueae, and Heliantheae emerged as dominant contributors to the flora, illustrating their adaptability to varied ecological conditions. Microphyll emerged as the dominant leaf type, and therophytes prevailed as the most abundant life form. Urban habitats exhibited heightened sensitivity to Asteraceae species, as reflected in diversity indices. Canonical Correspondence Analysis demonstrated the significant influence of soil moisture, pH, phosphorus, and soil saturation on species distribution, with indicator species analysis highlighting habitat-specific plant-environment relationships. The study further classified species into adaptive strategies using CSR (competitive strategy-stress tolerant strategy-ruderal strategy) categorization. While some species exhibited pronounced R or S strategies, a majority demonstrated mixed strategies, such as R/CR, R/CSR, CS/CSR, emphasizing the versatility of adaptation strategies within the Asteraceae family. The findings underscore the intricate interplay between environmental variables and Asteraceae phylogenetic and functional diversity, highlighting the significance of soil properties in influencing plant distribution in semi-arid ecosystems. This research contributes valuable insights into the ecological dynamics of Asteraceae, shedding light on their responses to diverse environmental conditions and providing a nuanced understanding of the factors influencing their spatial distribution.

Odd-Paired is Involved in Morphological Divergence of Snail-Feeding Beetles.

Body shape and size diversity and their evolutionary rates correlate with species richness at the macroevolutionary scale. However, the molecular genetic mechanisms underlying the morphological diversification across related species are poorly understood. In beetles, which account for one-fourth of the known species, adaptation to different trophic niches through morphological diversification appears to have contributed to species radiation. Here, we explored the key genes for the morphological divergence of the slender to stout body shape related to divergent feeding methods on large to small snails within the genus Carabus. We show that the zinc-finger transcription factor encoded by odd-paired (opa) controls morphological variation in the snail-feeding ground beetle Carabus blaptoides. Specifically, opa was identified as the gene underlying the slender to stout morphological difference between subspecies through genetic mapping and functional analysis via gene knockdown. Further analyses revealed that changes in opa cis-regulatory sequences likely contributed to the differences in body shape and size between C. blaptoides subspecies. Among opa cis-regulatory sequences, single nucleotide polymorphisms on the transcription factor binding sites may be associated with the morphological differences between C. blaptoides subspecies. opa was highly conserved in a wide range of taxa, especially in beetles. Therefore, opa may play an important role in adaptive morphological divergence in beetles.

Facilitated introgression from domestic goat into Alpine ibex at immune loci.

Hybridization can result in the transfer of adaptive genetic material from one species to another, known as adaptive introgression. Bottlenecked (and hence genetically depleted) species are expected to be particularly receptive to adaptive introgression, since introgression can introduce new or previously lost adaptive genetic variation. The Alpine ibex (Capra ibex), which recently recovered from near extinction, is known to hybridize with the domestic goat (Capra aegagrus hircus), and signals of introgression previously found at the major histocompatibility complex were suggested to potentially be adaptive. Here, we combine two ancient whole genomes of Alpine ibex with 29 modern Alpine ibex genomes and 31 genomes representing six related Capra species to investigate the genome-wide patterns of introgression and confirm the potential relevance of immune loci. We identified low rates of admixture in modern Alpine ibex through various F statistics and screening for putative introgressed tracts. Further results based on demographic modelling were consistent with introgression to have occurred during the last 300 years, coinciding with the known species bottleneck, and that in each generation, 1-2 out of 100 Alpine ibex had a domestic goat parent. The putatively introgressed haplotypes were enriched at immune-related genes, where the adaptive value of alternative alleles may give individuals with otherwise depleted genetic diversity a selective advantage. While interbreeding with domestic species is a prevalent issue in species conservation, in this specific case, it resulted in putative adaptive introgression. Our findings highlight the complex interplay between hybridization, adaptive evolution, and the potential risks and benefits associated with anthropogenic influences on wild species.

Genomics of hybrid parallel origin in Aquilegia ecalcarata

BackgroundThe parallel evolution of similar traits or species provides strong evidence for the role of natural selection in evolution. Traits or species that evolved repeatedly can be driven by separate de novo mutations or interspecific gene flow. Although parallel evolution has been reported in many studies, documented cases of parallel evolution caused by gene flow are scarce by comparison. Aquilegia ecalcarata and A. kansuensis belong to the genus of Aquilegia, and are the closest related sister species. Mutiple origins of A. ecalcarata have been reported in previous studies, but whether they have been driven by separate de novo mutations or gene flow remains unclear.ResultsIn this study, We conducted genomic analysis from 158 individuals of two repeatedly evolving pairs of A. ecalcarata and A. kansuensis. All samples were divided into two distinct clades with obvious geographical distribution based on phylogeny and population structure. Demographic modeling revealed that the origin of the A. ecalcarata in the Eastern of China was caused by gene flow, and the Eastern A. ecalcarata occurred following introgression from Western A. ecalcarata population. Analysis of Treemix and D-statistic also revealed that a strong signal of gene flow was detected from Western A. ecalcarata to Eastern A. ecalcarata. Genetic divergence and selective sweep analyses inferred parallel regions of genomic divergence and identified many candidate genes associated with ecologically adaptive divergence between species pair. Comparative analysis of parallel diverged regions and gene introgression confirms that gene flow contributed to the parallel evolution of A. ecalcarata.ConclusionsOur results further confirmed the multiple origins of A. ecalcarata and highlighted the roles of gene flow. These findings provide new evidence for parallel origin after hybridization as well as insights into the ecological adaptation mechanisms underlying the parallel origins of species.

Microbiome and epigenetic variation in wild fish with low genetic diversity

Non-genetic sources of phenotypic variation, such as the epigenome and the microbiome, could be important contributors to adaptive variation for species with low genetic diversity. However, little is known about the complex interaction between these factors and the genetic diversity of the host, particularly in wild populations. Here, we examine the skin microbiome composition of two closely-related mangrove killifish species with different mating systems (self-fertilising and outcrossing) under sympatric and allopatric conditions. This allows us to partition the influence of the genotype and the environment on their microbiome and (previously described) epigenetic profiles. We find the diversity and community composition of the skin microbiome are strongly shaped by the environment and, to a lesser extent, by species-specific influences. Heterozygosity and microbiome alpha diversity, but not epigenetic variation, are associated with the fluctuating asymmetry of traits related to performance (vision) and behaviour (aggression). Our study identifies that a proportion of the epigenetic diversity and microbiome differentiation is unrelated to genetic variation, and we find evidence for an associative relationship between microbiome and epigenetic diversity in these wild populations. This suggests that both mechanisms could potentially contribute to variation in species with low genetic diversity.

Genetic divergence and one-way gene flow influence contemporary evolution and ecology of a partially migratory fish.

Recent work has revealed the importance of contemporary evolution in shaping ecological outcomes. In particular, rapid evolutionary divergence between populations has been shown to impact the ecology of populations, communities, and ecosystems. While studies have focused largely on the role of adaptive divergence in generating ecologically important variation among populations, much less is known about the role of gene flow in shaping ecological outcomes. After divergence, populations may continue to interact through gene flow, which may influence evolutionary and ecological processes. Here, we investigate the role of gene flow in shaping the contemporary evolution and ecology of recently diverged populations of anadromous steelhead and resident rainbow trout (Oncorhynchus mykiss). Results show that resident rainbow trout introduced above waterfalls have diverged evolutionarily from downstream anadromous steelhead, which were the source of introductions. However, the movement of fish from above to below the waterfalls has facilitated gene flow, which has reshaped genetic and phenotypic variation in the anadromous source population. In particular, gene flow has led to an increased frequency of residency, which in turn has altered population density, size structure, and sex ratio. This result establishes gene flow as a contemporary evolutionary process that can have important ecological outcomes. From a management perspective, anadromous steelhead are generally regarded as a higher conservation priority than resident rainbow trout, even when found within the same watershed. Our results show that anadromous and resident O. mykiss populations may be connected via gene flow, with important ecological consequences. Such eco-evolutionary processes should be considered when managing recently diverged populations connected by gene flow.

Terminal regions of a protein are a hotspot for low complexity regions and selection.

Volatile low complexity regions (LCRs) are a novel source of adaptive variation, functional diversification and evolutionary novelty. An interplay of selection and mutation governs the composition and length of low complexity regions. High %GC and mutations provide length variability because of mechanisms like replication slippage. Owing to the complex dynamics between selection and mutation, we need a better understanding of their coexistence. Our findings underscore that positively selected sites (PSS) and low complexity regions prefer the terminal regions of genes, co-occurring in most Tetrapoda clades. We observed that positively selected sites within a gene have position-specific roles. Central-positively selected site genes primarily participate in defence responses, whereas terminal-positively selected site genes exhibit non-specific functions. Low complexity region-containing genes in the Tetrapoda clade exhibit a significantly higher %GC and lower ω (dN/dS: non-synonymous substitution rate/synonymous substitution rate) compared with genes without low complexity regions. This lower ω implies that despite providing rapid functional diversity, low complexity region-containing genes are subjected to intense purifying selection. Furthermore, we observe that low complexity regions consistently display ubiquitous prevalence at lower purity levels, but exhibit a preference for specific positions within a gene as the purity of the low complexity region stretch increases, implying a composition-dependent evolutionary role. Our findings collectively contribute to the understanding of how genetic diversity and adaptation are shaped by the interplay of selection and low complexity regions in the Tetrapoda clade.