Signatures Of Local Adaptation Research Articles

Complex patterns of genetic structure in the sea cucumber Holothuria (Metriatyla) scabra from the Philippines: implications for aquaculture and fishery management

The sandfish Holothuria (Metriatyla) scabra, is a high-value tropical sea cucumber harvested from wild stocks for over four centuries in multi-species fisheries across its Indo-Pacific distribution, for the global bêche-de-mer (BDM) trade. Within Southeast Asia, the Philippines is an important centre of the BDM trade, however overharvesting and largely open fishery management have resulted in declining catch volumes. Sandfish mariculture has been developed to supplement BDM supply and assist restocking efforts; however, it is heavily reliant on wild populations for broodstock supply. Consequently, to inform fishery, mariculture, germplasm and translocation management policies for both wild and captive resources, a high-resolution genomic audit of 16 wild sandfish populations was conducted, employing a proven genotyping-by-sequencing approach for this species (DArTseq). Genomic data (8,266 selectively-neutral and 117 putatively-adaptive SNPs) were used to assess fine-scale genetic structure, diversity, relatedness, population connectivity and local adaptation at both broad (biogeographic region) and local (within-biogeographic region) scales. An independent hydrodynamic particle dispersal model was also used to assess population connectivity. The overall pattern of population differentiation at the country level for H. scabra in the Philippines is complex, with nine genetic stocks and respective management units delineated across 5 biogeographic regions: (1) Celebes Sea, (2) North and (3) South Philippine Seas, (4) South China and Internal Seas and (5) Sulu Sea. Genetic connectivity is highest within proximate marine biogeographic regions (mean Fst=0.016), with greater separation evident between geographically distant sites (Fst range=0.041–0.045). Signatures of local adaptation were detected among six biogeographic regions, with genetic bottlenecks at 5 sites, particularly within historically heavily-exploited locations in the western and central Philippines. Genetic structure is influenced by geographic distance, larval dispersal capacity, species-specific larval development and settlement attributes, variable ocean current-mediated gene flow, source and sink location geography and habitat heterogeneity across the archipelago. Data reported here will inform accurate and sustainable fishery regulation, conservation of genetic diversity, direct broodstock sourcing for mariculture and guide restocking interventions across the Philippines.

Shared xerophytic genes and their re-use in local adaptation to aridity in the desert plant Gymnocarpos przewalskii.

In order to thrive and survive, plant species need to combine stability in the long term and rapid response to environmental challenges in the short term. The former would be reflected by parallel or convergent adaptation across species, and the latter by pronounced local adaptation among populations of the same species. In the present study, we generated a high-quality genome and re-sequenced 177 individuals for Gymnocarpos przewalskii, an important desert plant species from North-West China, to detect local adaptation. We first focus on ancient adaptation to aridity at the molecular level by comparing the genomic data of 15 species that vary in their ability to withstand aridity. We found that a total of 118 genes were shared across xerophytic species but absent from non-xerophytic species. Of the 65 found in G. przewalskii, 63 were under purifying selection and two under positive selection. We then focused on local adaptation. Up to 20% of the G. przewalskii genome showed signatures of local adaptation to aridity during population divergence. Thirteen of the selected shared xerophytic genes were reused in local adaptation after population differentiation. Hence, only about 20% of the genes shared and specific to xerophytic species and associated with adaptation to aridity were later recruited for local adaptation in G. przewalskii.

Signatures of adaptation at key insecticide resistance loci in Anopheles gambiae in Southern Ghana revealed by reduced-coverage WGS

Resistance to insecticides and adaptation to a diverse range of environments present challenges to Anopheles gambiae s.l. mosquito control efforts in sub-Saharan Africa. Whole-genome-sequencing is often employed for identifying the genomic basis underlying adaptation in Anopheles, but remains expensive for large-scale surveys. Reduced coverage whole-genome-sequencing can identify regions of the genome involved in adaptation at a lower cost, but is currently untested in Anopheles mosquitoes. Here, we use reduced coverage WGS to investigate population genetic structure and identify signatures of local adaptation in Anopheles mosquitoes across southern Ghana. In contrast to previous analyses, we find no structuring by ecoregion, with Anopheles coluzzii and Anopheles gambiae populations largely displaying the hallmarks of large, unstructured populations. However, we find signatures of selection at insecticide resistance loci that appear ubiquitous across ecoregions in An. coluzzii, and strongest in forest ecoregions in An. gambiae. Our study highlights resistance candidate genes in this region, and validates reduced coverage WGS, potentially to very low coverage levels, for population genomics and exploratory surveys for adaptation in Anopheles taxa.

The drivers of intraspecific trait variation and their implications for future tree productivity and survival.

Forests are facing unprecedented levels of stress from pest and disease outbreaks, disturbance, fragmentation, development, and a changing climate. These selective agents act to alter forest composition from regional to cellular levels. Thus, a central challenge for understanding how forests will be impacted by future change is how to integrate across scales of biology. Phenotype, or an observable trait, is the product of an individual's genes (G) and the environment in which an organism lives (E). To date, researchers have detailed how environment drives variation in tree phenotypes over long time periods (e.g., long-term ecological research sites [LTERs]) and across large spatial scales (e.g., flux network). In parallel, researchers have discovered the genes and pathways that govern phenotypes, finding high degrees of genetic control and signatures of local adaptation in many plant traits. However, the research in these two areas remain largely independent of each other, hindering our ability to generate accurate predictions of plant response to environment, an increasingly urgent need given threats to forest systems. I present the importance of both genes and environment in determining tree responses to climate stress. I highlight why the difference between G versus E in driving variation is critical for our understanding of climate responses, then propose means of accelerating research that examines G and E simultaneously by leveraging existing long-term, large-scale phenotypic data sets from ecological networks and adding newly affordable sequence (-omics) data to both drill down to find the genes and alleles influencing phenotypes and scale up to find how patterns of demography and local adaptation may influence future response to change.

Patterns of genetic variation and local adaptation of a native herbivore to a lethal invasive plant.

Understanding the evolutionary processes that influence fitness is critical to predicting species' responses to selection. Interactions among evolutionary processes including gene flow, drift and the strength of selection can lead to either local adaptation or maladaptation, especially in heterogenous landscapes. Populations experiencing novel environments or resources are ideal for understanding the mechanisms underlying adaptation or maladaptation, specifically in locally co-evolved interactions. We used the interaction between a native herbivore that oviposits on a patchily distributed introduced plant that in turn causes significant mortality to the larvae to test for signatures of local adaptation in areas where the two co-occurred. We used whole-genome sequencing to explore population structure, patterns of gene flow and signatures of local adaptation. We found signatures of local adaptation in response to the introduced plant in the absence of strong population structure with no genetic differentiation and low genetic variation. Additionally, we found localized allele frequency differences within a single population between habitats with and without the lethal plant, highlighting the effects of strong selection. Finally, we identified that selection was acting on larval ability to feed on the plant rather than on females' ability to avoid oviposition, thus uncovering the specific ontogenetic target of selection. Our work highlights the potential for adaptation to occur in a fine-grained landscape in the presence of gene flow and low genetic variation.

Mosaic of local adaptation between white clover and rhizobia along an urbanization gradient

Abstract Urbanization frequently alters the biotic and abiotic environmental context. Variation in both the biotic and abiotic environment can result in selection mosaics that lead to locally adapted phenotypes. Ecological changes to urban environments can result in evolutionary responses across a diverse array of plants and animals, but the effects of urbanization on local adaptation and coevolution between species remain largely unstudied. Using a common garden experiment of 30 populations and 1080 plants, we tested for local adaptation in the mutualism between white clover (Trifolium repens) and rhizobia (Rhizobium leguminosarum symbiovar trifolii) along an urbanization gradient. We asked: (Q1) are white clover and rhizobia locally adapted? (Q2) Does nitrogen (N) addition mediate the strength of local adaptation? (Q3) Is the strength of local adaptation in host plants related to the strength of local adaptation in rhizobia, as expected if white clover and rhizobia are coadapting? And (Q4) how does the degree of local adaptation vary with urbanization? We found evidence for local adaptation for both white clover and rhizobia, with stronger signatures of local adaptation for rhizobia than white clover (Q1). While N addition positively affected plant biomass and negatively affected nodule density, N addition did not consistently mediate patterns of local adaptation (Q2). The strength of local adaptation was positively related between white clover and rhizobia under N addition, suggesting that soil N mediates coadaptation between white clover and rhizobia (Q3). We did not find a consistent relationship between measures of urbanization and local adaptation (Q4). Additionally, fitness responses for white clover and rhizobia were unrelated to Rhizobium abundance and local adaptation was unrelated to the broader bacterial microbiome in the soil inoculants or root endosphere. Synthesis. Our results show a spatial mosaic in local adaptation, with stronger local adaptation for rhizobia than white clover. While urbanization does influence the ecology of plant–microbe interactions, our study suggests urbanization does not disrupt local adaptation and coevolution among mutualists.

Genomic responses to climate: Understanding local adaptation in the Andean tree species Nothofagus pumilio and implications for a changing world

Societal Impact StatementForest trees tend to be strongly genetically adapted to their local environments, but climate change will probably subject trees to novel combinations of precipitation, temperature, and photoperiod. Local adaptation was investigated in the ecologically and economically important Patagonian tree species Nothofagus pumilio by characterizing its genetic diversity in relation to the varied environmental conditions across its range. These insights are useful for conservation and management decisions, for example by identifying suitable populations to establish seed source plantations for restoration and characterizing relationships with environmental drivers of selection to better understand how this species will respond to climate change.SummaryNothofagus pumilio is a foundation tree species that inhabits a 2000‐km‐long range in the southern Andes, a region with two perpendicular environmental gradients: temperature and photoperiod (North–South), and precipitation (West–East). We investigated local adaptation patterns by searching for relationships between environmental clines and signatures of adaptation in candidate genes related to stress response, growth, and phenology. Using a paired site sampling design within a landscape genome analysis, we analyzed 493 adult N. pumilio trees in 20 sampling sites across the species' latitudinal range. We screened 47,336 single nucleotide polymorphism (SNP) loci in 1632 contigs (i.e., coding regions along the genome). Population structure and genetic diversity analyses preceded four genome scan analyses using genetic and environmental data. Population structure and genetic diversity are mainly oriented along the latitude axis. Genome scans identified 445 outlier SNPs, which are loci showing signatures of selection. Temperature and photoperiod variables were associated with notably more outliers than precipitation. However, the most frequent biological functions among genes were water deprivation response and cold response, suggesting that stress response is comprised of complex and polygenic traits that are affected by many environmental variables. Our findings suggest that N. pumilio shows signatures of local adaptation to extant climate conditions, including temperature, photoperiod, and precipitation. However, climate change is likely to alter existing relationships among environmental conditions to which this species is currently adapted. These changes may have unpredictable consequences for the species' future survival, adaptation potential, and the people who depend upon these forests.

Signatures of local adaptation to current and future climate in phenology-related genes in natural populations of Quercus robur

BackgroundLocal adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change.ResultsHere, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 FST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change.ConclusionsThe model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

Divergent gametic thermal performance and floral warming across an elevation gradient.

Thermal environments vary widely across species ranges, establishing the potential for local adaptation of thermal performance optima and tolerance. In the absence of local adaptation, selection should favor mechanisms to meet thermal optima. Floral temperature is a major determinant of reproductive success in angiosperms, yet whether gametic thermal performance shows signatures of local adaptation across temperature gradients, and how variation in gametic thermal performance influences floral evolution, is unknown. We characterized flowering season temperatures for the forb, Argentina anserina, at extremes of a 1000 m elevation gradient and generated thermal performance curves for pollen and ovule performance in populations at each extreme. Thermal optima fell between mean and maximum intrafloral temperatures. However, cooler high-elevation populations had ~4 °C greater pollen thermal optima than warmer low-elevation populations, while tolerance breadths did not differ. We then tested whether plants at elevational extremes differentially warmed the floral microenvironment. High-elevation flowers warmed significantly more than low, bringing intrafloral temperatures nearer the pollen optima. A manipulative experiment demonstrated that stronger warming in high elevation was conferred by floral tissues. Elevational divergence in floral warming may be driven, in part, by selection on flowers to meet different thermal demands of the gametophytes.

Parallel shifts in flight-height associated with altitude across incipient Heliconius species.

Vertical gradients in microclimate, resource availability, and interspecific interactions are thought to underly stratification patterns in tropical insect communities. However, only a few studies have explored the adaptive significance of vertical space use during the early stages of reproductive isolation. We analysed flight-height variation across speciation events in Heliconius butterflies, representing parallel colonizations of high-altitude forest. We measured flight-height in wild H. erato venus and H. chestertonii, parapatric lowland and mountain specialists, respectively, and found that H. chestertonii consistently flies at a lower height. By comparing our data to previously published results for the ecologically equivalent H. e. cyrbia (lowland) and H. himera (high altitude), we found that the species flying closest to the ground are those that recently colonized high-altitude forests. We show that these repeated trends largely result from shared patterns of ecological selection producing parallel trait-shifts in H. himera and H. chestertonii. Although our results imply a signature of local adaptation, we did not find an association between resource distribution and flight-height in H. e. venus and H. chestertonii. We discuss how this pattern may be explained by variations in forest structure and microclimate. Overall, our findings underscore the importance of behavioural adjustments during early divergence mediated by altitude-shifts.

Divergent selection in a Mediterranean pine on local spatial scales

Abstract The effects of selection on an organism's genome are hard to detect on small spatial scales, as gene flow can swamp signatures of local adaptation. Therefore, most genome scans to detect signatures of environmental selection are performed on large spatial scales; however, divergent selection on the local scale (e.g. between contrasting soil conditions) has also been demonstrated, in particular for herbaceous plants. Here, we hypothesised that in topographically complex landscapes, microenvironment variability is strong enough to leave a selective footprint in the genomes of long‐lived organisms. To test this, we investigated paired south‐ versus north‐facing Pinus pinaster stands on the local scale, with trees growing in close vicinity (≤820 m distance between paired south‐ and north‐facing stands), in a Mediterranean mountain area. While trees on north‐facing slopes experience less radiation, trees on south‐facing slopes suffer from especially harsh conditions, particularly during the dry summer season. Two outlier analyses consistently revealed five putatively adaptive loci (out of 4034), in candidate genes two of which encoded non‐synonymous substitutions. Additionally, one locus showed consistent allele frequency differences in all three stand pairs indicating divergent selection despite high gene flow on the local scale. Permutation tests demonstrated that our findings were robust. Functional annotation of these candidate genes revealed biological functions related to abiotic stress response, such as water availability, in other plant species. Synthesis. Our study highlights how divergent selection in heterogeneous microenvironments shapes and maintains the functional genetic variation within populations of long‐lived forest tree species, being the first to focus on adaptive genetic divergence between south‐ and north‐facing slopes within continuous forest stands. This is especially relevant in the current context of climate change, as this variation is at the base of plant population responses to future climate.

Genetic architecture of the response of Arabidopsis thaliana to a native plant-growth-promoting bacterial strain.

By improving plant nutrition and alleviating abiotic and biotic stresses, plant growth-promoting bacteria (PGPB) can help to develop eco-friendly and sustainable agricultural practices. Besides climatic conditions, soil conditions, and microbe-microbe interactions, the host genotype influences the effectiveness of PGPB. Yet, most GWAS conducted to characterize the genetic architecture of response to PGPB are based on non-native interactions between a host plant and PGPB strains isolated from the belowground compartment of other plants. In this study, a GWAS was set up under in vitro conditions to describe the genetic architecture of the response of Arabidopsis thaliana to the PGPB Pseudomonas siliginis, by inoculating seeds of 162 natural accessions from the southwest of France with one strain isolated from the leaf compartment in the same geographical region. Strong genetic variation of plant growth response to this native PGPB was observed at a regional scale, with the strain having a positive effect on the vegetative growth of small plants and a negative effect on the vegetative growth of large plants. The polygenic genetic architecture underlying this negative trade-off showed suggestive signatures of local adaptation. The main eco-evolutionary relevant candidate genes are involved in seed and root development.

Genomic signals of local adaptation across climatically heterogenous habitats in an invasive tropical fruit fly (Bactrocera tryoni)

Local adaptation plays a key role in the successful establishment of pest populations in new environments by enabling them to tolerate novel biotic and abiotic conditions experienced outside their native range. However, the genomic underpinnings of such adaptive responses remain unclear, especially for agriculturally important pests. We investigated population genomic signatures in the tropical/subtropical Queensland fruit fly, Bactrocera tryoni, which has an expanded range encompassing temperate and arid zones in Australia, and tropical zones in the Pacific Islands. Using reduced representation sequencing data from 28 populations, we detected allele frequency shifts associated with the native/invasive status of populations and identified environmental factors that have likely driven population differentiation. We also determined that precipitation, temperature, and geographic variables explain allelic shifts across the distribution range of B. tryoni. We found spatial heterogeneity in signatures of local adaptation across various climatic conditions in invaded areas. Specifically, disjunct invasive populations in the tropical Pacific Islands and arid zones of Australia were characterised by multiple significantly differentiated single nucleotide polymorphisms (SNPs), some of which were associated with genes with well-understood function in environmental stress (e.g., heat and desiccation) response. However, invasive populations in southeast Australian temperate zones showed higher gene flow with the native range and lacked a strong local adaptive signal. These results suggest that population connectivity with the native range has differentially affected local adaptive patterns in different invasive populations. Overall, our findings provide insights into the evolutionary underpinnings of invasion success of an important horticultural pest in climatically distinct environments.

Heritable variation in thermal profiles is associated with reproductive success in the world's largest bird.

Organisms inhabiting extreme thermal environments, such as desert birds, have evolved spectacular adaptations to thermoregulate during hot and cold conditions. However, our knowledge of selection for thermoregulation and the potential for evolutionary responses is limited, particularly for large organisms experiencing extreme temperature fluctuations. Here we use thermal imaging to quantify selection and genetic variation in thermoregulation in ostriches (Struthio camelus), the world's largest bird species that is experiencing increasingly volatile temperatures. We found that females who are better at regulating their head temperatures ("thermoregulatory capacity") had higher egg-laying rates under hotter conditions. Thermoregulatory capacity was both heritable and showed signatures of local adaptation: females originating from more unpredictable climates were better at regulating their head temperatures in response to temperature fluctuations. Together these results reveal that past and present evolutionary processes have shaped genetic variation in thermoregulatory capacity, which appears to protect critical organs, such as the brain, from extreme temperatures during reproduction.

High-performing plastic clones best explain the spread of yellow monkeyflower from lowland to higher elevation areas in New Zealand.

The relative contribution of adaptation and phenotypic plasticity can vary between core and edge populations, with implications for invasive success. We investigated the spread of the invasive yellow monkeyflower, Erythranthe gutatta in New Zealand, where it is spreading from lowland agricultural land into high-elevation conservation areas. We investigated the extent of phenotypic variation among clones from across the South Island, looked for adaptation and compared degrees of plasticity among lowland core versus montane range-edge populations. We grew 34 clones and measured their vegetative and floral traits in two common gardens, one in the core range at 9 m a.s.l. and one near the range-edge at 560 m a.s.l. Observed trait variation was explained by a combination of genotypic diversity (as identified through common gardens) and high phenotypic plasticity. We found a subtle signature of local adaptation to lowland habitats but all clones were plastic and able to survive and reproduce in both gardens. In the range-edge garden, above-ground biomass was on average almost double and stolon length almost half that of the same clone in the core garden. Clones from low-elevation sites showed higher plasticity on average than those from higher elevation sites. The highest performing clones in the core garden were also top performers in the range-edge garden. These results suggest some highly fit general-purpose genotypes, possibly pre-adapted to New Zealand montane conditions, best explains the spread of E. gutatta from lowland to higher elevation areas.

Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production.

Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.

Multi-locus genomic signatures of local adaptation to snow across the landscape in California populations of a willow leaf beetle.

Organisms living in mountains contend with extreme climatic conditions, including short growing seasons and long winters with extensive snow cover. Anthropogenic climate change is driving unprecedented, rapid warming of montane regions across the globe, resulting in reduced winter snowpack. Loss of snow as a thermal buffer may have serious consequences for animals overwintering in soil, yet little is known about how variability in snowpack acts as a selective agent in montane ecosystems. Here, we examine genomic variation in California populations of the leaf beetle Chrysomela aeneicollis, an emerging natural model system for understanding how organisms respond to climate change. We used a genotype-environment association approach to identify genomic signatures of local adaptation to microclimate in populations from three montane regions with variable snowpack and a coastal region with no snow. We found that both winter-associated environmental variation and geographical distance contribute to overall genomic variation across the landscape. We identified non-synonymous variation in novel candidate loci associated with cytoskeletal function, ion transport and membrane stability, cellular processes associated with cold tolerance in other insects. These findings provide intriguing evidence that variation in snowpack imposes selective gradients in montane ecosystems.

The meaning of wild: Genetic and adaptive consequences from large-scale releases of domestic mallards

The translocation of individuals around the world is leading to rising incidences of anthropogenic hybridization, particularly between domestic and wild congeners. We apply a landscape genomics approach for thousands of mallard (Anas platyrhynchos) samples across continental and island populations to determine the result of over a century of supplementation practices. We establish that a single domestic game-farm mallard breed is the source for contemporary release programs in Eurasia and North America, as well as for established feral populations in New Zealand and Hawaii. In particular, we identify central Europe and eastern North America as epicenters of ongoing anthropogenic hybridization, and conclude that the release of game-farm mallards continues to affect the genetic integrity of wild mallards. Conversely, self-sustaining feral populations in New Zealand and Hawaii not only show strong differentiation from their original stock, but also signatures of local adaptation occurring in less than a half-century since game-farm mallard releases have ceased. We conclude that ‘wild’ is not singular, and that even feral populations are capable of responding to natural processes. Although considered paradoxical to biological conservation, understanding the capacity for wildness among feral and feral admixed populations in human landscapes is critical as such interactions increase in the Anthropocene.

Genomic analysis of the emergent aquatic plant Sparganium stoloniferum provides insights into its clonality, local adaptation and demographic history.

Clonal propagation and extensive dispersal of seeds and asexual propagules are two important features of aquatic plants that help them adapt to aquatic environments. Accurate measurements of clonality and effective clonal dispersal are essential for understanding the evolution of aquatic plants. Here, we first assembled a high-quality chromosome-level genome of a widespread emergent aquatic plant Sparganium stoloniferum to provide a reference for its population genomic study. We then performed high-depth resequencing of 173 individuals from 20 populations covering different basins across its range in China. Population genomic analyses revealed three genetic lineages reflecting the northeast (NE), southwest (SW) and northwest (NW) of its geographical distribution. The NE lineage diverged in the middle Pleistocene while the SW and NW lineages diverged until about 2400 years ago. Clonal relationship analyses identified nine populations as monoclonal population. Dispersal of vegetative propagules was identified between five populations covering three basins in the NE lineage, and dispersal distance was up to 1041 km, indicating high dispersibility in emergent aquatic plant species. We also identified lineage-specific positively selected genes that are likely to be involved in adaptations to saline wetlands and high-altitude environments. Our findings accurately measure the clonality, determine the dispersal range and frequency of vegetative propagules, and detect genetic signatures of local adaptation in a widespread emergent aquatic plant species, providing new perspectives on the evolution of aquatic plants.

Human subsistence and signatures of selection on chemosensory genes

Chemosensation (olfaction, taste) is essential for detecting and assessing foods, such that dietary shifts elicit evolutionary changes in vertebrate chemosensory genes. The transition from hunting and gathering to agriculture dramatically altered how humans acquire food. Recent genetic and linguistic studies suggest agriculture may have precipitated olfactory degeneration. Here, we explore the effects of subsistence behaviors on olfactory (OR) and taste (TASR) receptor genes among rainforest foragers and neighboring agriculturalists in Africa and Southeast Asia. We analyze 378 functional OR and 26 functional TASR genes in 133 individuals across populations in Uganda (Twa, Sua, BaKiga) and the Philippines (Agta, Mamanwa, Manobo) with differing subsistence histories. We find no evidence of relaxed selection on chemosensory genes in agricultural populations. However, we identify subsistence-related signatures of local adaptation on chemosensory genes within each geographic region. Our results highlight the importance of culture, subsistence economy, and drift in human chemosensory perception.