Patterns of genetic diversity and gene flow associated with an aridity gradient in populations of common mole-rats, Cryptomys hottentotus hottentotus.

Genetic adaptation is the change of a population towards a phenotype that best fits the present ecological conditions in which it inhabits. As environmental conditions change, allele frequencies shift, resulting in different populations of the same species possessing genetic variation and divergent phenotypes. Cooperatively breeding common mole-rats (Cryptomys hottentotus hottentotus) inhabit environments along an aridity gradient in South Africa, which provides an opportunity for local genetic adaptations to occur. Using one mitochondrial gene (cytochrome b) and 3,540 SNP loci across the whole genome, we determined the phylogenetic relationship, population structure and genetic diversity of five populations of C. h. hottentotus located along an aridity gradient. Mitochondrial data identified population-specific clades that were less distinct in the two mesic populations, potentially indicating historical or recent gene flow, or the retention of ancestral haplotypes. Arid and semi-arid populations formed a distinct cluster from the non-arid populations. Genetic diversity and gene flow were higher in arid-dwelling individuals, suggesting greater connectivity and interactions between colonies in arid regions in comparison to mesic ones. Using an Aridity Index, we determined that isolation by environment, rather than isolation by geographical distance, best explains the genetic distance between the populations. Further analyses using target loci may determine if there are differing underlying genetic adaptations among populations of C. h. hottentotus. These analyses could help unravel population differences in response to environmental factors within a sub-species of bathyergid mole-rat and determine the adaptive capacity of this small non-migratory subterranean rodent species in response to aridification in the face of climate change.

Limited Parallelism in Genetic Adaptation to Brackish Water Bodies in European Sprat and Atlantic Herring.

The European sprat is a small plankton-feeding clupeid present in the northeastern Atlantic Ocean, in the Mediterranean Sea, and in the brackish Baltic Sea and Black Sea. This species is the target of a major fishery and, therefore, an accurate characterization of its genetic population structure is crucial to delineate proper stock assessments that aid ensuring the fishery's sustainability. Here, we present (i) a draft genome assembly, (ii) pooled whole genome sequencing of 19 population samples covering most of the species' distribution range, and (iii) the design and test of a single nucleotide polymorphism (SNP)-chip resource and use this to validate the population structure inferred from pooled sequencing. These approaches revealed, using the populations sampled here, three major groups of European sprat: Oceanic, Coastal, and Brackish with limited differentiation within groups even over wide geographical stretches. Genetic structure is largely driven by six large putative inversions that differentiate Oceanic and Brackish sprats, while Coastal populations display intermediate frequencies of haplotypes at each locus. Interestingly, populations from the Baltic and the Black Seas share similar frequencies of haplotypes at these putative inversions despite their distant geographic location. The closely related clupeids European sprat and Atlantic herring both show genetic adaptation to the brackish Baltic Sea, providing an opportunity to explore the extent of genetic parallelism. This analysis revealed limited parallelism because out of 125 independent loci detected in the Atlantic herring, three showed sharp signals of selection that overlapped between the two species and contained single genes such as PRLRA, which encodes the receptor for prolactin, a freshwater-adapting hormone in euryhaline species, and THRB, a receptor for thyroid hormones, important both for metabolic regulation and the development of red cone photoreceptors.

Sex-biased regulatory changes in the placenta of native highlanders contribute to adaptive fetal development.

Compared with lowlander migrants, native Tibetans have a higher reproductive success at high altitude though the underlying mechanism remains unclear. Here, we compared the transcriptome and histology of full-term placentas between native Tibetans and Han migrants. We found that the placental trophoblast shows the largest expression divergence between Tibetans and Han, and Tibetans show decreased immune response and endoplasmic reticulum stress. Remarkably, we detected a sex-biased expression divergence, where the male-infant placentas show a greater between-population difference than the female-infant placentas. The umbilical cord plays a key role in the sex-biased expression divergence, which is associated with the higher birth weight of the male newborns of Tibetans. We also identified adaptive histological changes in the male-infant placentas of Tibetans, including larger umbilical artery wall and umbilical artery intima and media, and fewer syncytial knots. These findings provide valuable insights into the sex-biased adaptation of human populations, with significant implications for medical and genetic studies of human reproduction.

Weakened tanning ability is an important mechanism for evolutionary skin lightening in East Asians.

The evolution of light-skin pigmentation among Eurasians is considered as an adaptation to the high-latitude environments. East Asians are ideal populations for studying skin color evolution because of the complex environment of East Asia. Here, we report a strong selection signal for the pigmentation gene phenylalanine hydroxylase (PAH) in light-skinned Han Chinese individuals. The intron mutation rs10778203 in PAH is enriched in East Asians and is significantly associated with skin color of the back of the hand in Han Chinese males (P < 0.05). In vitro luciferase and transcription factor binding assays show that the ancestral allele of rs10778203 could bind to SMAD2 and has a significant enhancer activity for PAH. However, the derived T allele (the major allele in East Asians) of rs10778203 decreases the binding activity of transcription factors and enhancer activity. Meanwhile, the derived T allele of rs10778203 shows a weaker ultraviolet radiation response in A375 cells and zebrafish embryos. Furthermore, rs10778203 decreases melanin production in transgenic zebrafish embryos after ultraviolet B (UVB) treatment. Collectively, PAH is a potential pigmentation gene that regulates skin tanning ability. Natural selection has enriched the adaptive allele, resulting in weakened tanning ability in East Asians, suggesting a unique genetic mechanism for evolutionary skin lightening in East Asians.

Exploring Adaptive Phenotypes for the Human Calcium-Sensing Receptor Polymorphism R990G.

Rainforest hunter-gatherers from Southeast Asia are characterized by specific morphological features including a particularly dark skin color (D), short stature (S), woolly hair (W), and the presence of steatopygia (S)-fat accumulation localized in the hips (DSWS phenotype). Based on previous evidence in the Andamanese population, we first characterized signatures of adaptive natural selection around the calcium-sensing receptor gene in Southeast Asian rainforest groups presenting the DSWS phenotype and identified the R990G substitution (rs1042636) as a putative adaptive variant for experimental follow-up. Although the calcium-sensing receptor has a critical role in calcium homeostasis by directly regulating the parathyroid hormone secretion, it is expressed in different tissues and has been described to be involved in many biological functions. Previous works have also characterized the R990G substitution as an activating polymorphism of the calcium-sensing receptor associated with hypocalcemia. Therefore, we generated a knock-in mouse for this substitution and investigated organismal phenotypes that could have become adaptive in rainforest hunter-gatherers from Southeast Asia. Interestingly, we found that mouse homozygous for the derived allele show not only lower serum calcium concentration but also greater body weight and fat accumulation, probably because of enhanced preadipocyte differentiation and lipolysis impairment resulting from the calcium-sensing receptor activation mediated by R990G. We speculate that such differential features in humans could have facilitated the survival of hunter-gatherer groups during periods of nutritional stress in the challenging conditions of the Southeast Asian tropical rainforests.

Assessing genome-wide adaptations associated with range expansion in the pink rice borer, Sesamia inferens.

Understanding the genetic basis of adaptive evolution following habitat expansion can have important implications for pest management. The pink rice borer (PRB), Sesamia inferens (Walker), is a destructive pest of rice that was historically restricted to regions south of 34° N latitude in China. However, with changes in global climate and farming practices, the distribution of this moth has progressively expanded, encompassing most regions in North China. Here, 3 highly differentiated subpopulations were discovered using high-quality single-nucleotide polymorphism and structural variant datasets across China, corresponding to northern, southern China regions, and the Yunnan-Guizhou Plateau, with significant patterns of isolation by geographic and environmental distances. Our estimates of evolutionary history indicate asymmetric migration with varying population sizes across the 3 subpopulations. Selective sweep analyses estimated strong selection at insect cuticle glycine-rich cuticular protein genes which are associated with enhanced desiccation adaptability in the northern group, and at the histone-lysine-N-methyltransferase gene associated with range expansion and local adaptation in the Shandong population. Our findings have significant implications for the development of effective strategies to control this pest.

Hypoxia Inducible Factor pathway proteins in high-altitude mammals.

Humans and other mammals inhabit hypoxic high-altitude locales. In many of these species, genes under positive selection include ones in the Hypoxia Inducible Factor (HIF) pathway. One is PHD2 (EGLN1), which encodes for a key oxygen sensor. Another is HIF2A (EPAS1), which encodes for a PHD2-regulated transcription factor. Recent studies have provided insights into mechanisms for these high-altitude alleles. These studies have (i) shown that selection can occur on nonconserved, unstructured regions of proteins, (ii) revealed that high altitude-associated amino acid substitutions can have differential effects on protein-protein interactions, (iii) provided evidence for convergent evolution by different molecular mechanisms, and (iv) suggested that mutations in different genes can complement one another to produce a set of adaptive phenotypes.

Genetic diet interactions of ACE: the increased hypertension predisposition in the Latin American population.

Hypertension is one of the primary risk factors associated with cardiovascular diseases (CVDs). It is a condition that affects people worldwide, and its prevalence is increasing due to several factors, such as lack of physical activity, population aging, and unhealthy diets. Notably, this increase has primarily occurred in low and middle-income countries (LMICs). In Latin America, approximately 40% of adults have been diagnosed with hypertension. Moreover, reports have shown that the Latin American genetic composition is highly diverse, and this genetic background can influence various biological processes, including disease predisposition and treatment effectiveness. Research has shown that Western dietary patterns, which include increased consumption of red meat, refined grains, sugar, and ultra-processed food, have spread across the globe, including Latin America, due to globalization processes. Furthermore, a higher than recommended sodium consumption, which has been associated with hypertension, has been identified across different regions, including Asia, Europe, America, Oceania, and Africa. In conclusion, hypertension is a multifactorial disease involving environmental and genetic factors. In Latin America, hypertension prevalence is increasing due to various factors, including age, the adoption of a "Westernized" diet, and potential genetic predisposition factors involving the ACE gene. Furthermore, identifying the genetic and molecular mechanisms of the disease, its association with diet, and how they interact is essential for the development of personalized treatments to increase its efficacy and reduce side effects.

Comparative Population Transcriptomics Provide New Insight into the Evolutionary History and Adaptive Potential of World Ocean Krill.

Genetic variation is instrumental for adaptation to changing environments but it is unclear how it is structured and contributes to adaptation in pelagic species lacking clear barriers to gene flow. Here, we applied comparative genomics to extensive transcriptome datasets from 20 krill species collected across the Atlantic, Indian, Pacific, and Southern Oceans. We compared genetic variation both within and between species to elucidate their evolutionary history and genomic bases of adaptation. We resolved phylogenetic interrelationships and uncovered genomic evidence to elevate the cryptic Euphausia similis var. armata into species. Levels of genetic variation and rates of adaptive protein evolution vary widely. Species endemic to the cold Southern Ocean, such as the Antarctic krill Euphausia superba, showed less genetic variation and lower evolutionary rates than other species. This could suggest a low adaptive potential to rapid climate change. We uncovered hundreds of candidate genes with signatures of adaptive evolution among Antarctic Euphausia but did not observe strong evidence of adaptive convergence with the predominantly Arctic Thysanoessa. We instead identified candidates for cold-adaptation that have also been detected in Antarctic fish, including genes that govern thermal reception such as TrpA1. Our results suggest parallel genetic responses to similar selection pressures across Antarctic taxa and provide new insights into the adaptive potential of important zooplankton already affected by climate change.

Rapid genetic adaptation to a novel ecosystem despite a large founder event.

Introduced and invasive species make excellent natural experiments for investigating rapid evolution. Here, we describe the effects of genetic drift and rapid genetic adaptation in pink salmon (Oncorhynchus gorbuscha) that were accidentally introduced to the Great Lakes via a single introduction event 31 generations ago. Using whole-genome resequencing for 134 fish spanning five sample groups across the native and introduced range, we estimate that the source population's effective population size was 146,886 at the time of introduction, whereas the founding population's effective population size was just 72-a 2040-fold decrease. As expected with a severe founder event, we show reductions in genome-wide measures of genetic diversity, specifically a 37.7% reduction in the number of SNPs and an 8.2% reduction in observed heterozygosity. Despite this decline in genetic diversity, we provide evidence for putative selection at 47 loci across multiple chromosomes in the introduced populations, including missense variants in genes associated with circadian rhythm, immunological response and maturation, which match expected or known phenotypic changes in the Great Lakes. For one of these genes, we use a species-specific agent-based model to rule out genetic drift and conclude our results support a strong response to selection occurring in a period gene (per2) that plays a predominant role in determining an organism's daily clock, matching large day length differences experienced by introduced salmon during important phenological periods. Together, these results inform how populations might evolve rapidly to new environments, even with a small pool of standing genetic variation.

Population Genomics for Insect Conservation.

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

Stay in shape: Assessing the adaptive potential of shell morphology and its sensitivity to temperature in the invasive New Zealand mud snail Potamopyrgus antipodarum through phenotypic plasticity and natural selection in Europe.

Climate change may force organisms to adapt genetically or plastically to new environmental conditions. Invasive species show remarkable potential for rapid adaptation. The ovoviviparous New Zealand mud snail (NZMS), Potamopyrgus antipodarum, has successfully established across Europe with two clonally reproducing mitochondrial lineages since its arrival in the first half of the 19th century. Its remarkable variation in shell morphology was shown to be fitness relevant. We investigated the effects of temperature on shell morphology across 11 populations from Germany and the Iberian Peninsula in a common garden across three temperatures. We analyzed size and shape using geometric morphometrics. For both, we compared reaction norms and estimated heritabilities. For size, the interaction of temperature and haplotype explained about 50% of the total variance. We also observed more genotype by environment interactions indicating a higher degree of population differentiation than in shape. Across the three temperatures, size followed the expectations of the temperature-size rule, with individuals growing larger in cold environments. Changes in shape may have compensated for changes in size affecting space for brooding embryos. Heritability estimates were relatively high. As indicated by the very low coefficients of variation for clonal repeatability (CV A ), they can probably not be compared in absolute terms. However, they showed some sensitivity to temperature, in haplotype t more so than in z, which was only found in Portugal. The low CV A values indicate that genetic variation among European populations is still restricted with a low potential to react to selection. A considerable fraction of the genetic variation was due to differences between the clonal lineages. The NZMS has apparently not been long enough in Europe to accumulate significant genetic variation relevant for morphological adaptation. As temperature is obviously not the sole factor influencing shell morphology, their interaction will probably not be a factor limiting population persistence under a warming climate in Europe.

The blunted vascular endothelial growth factor-A (VEGF-A) response to high-altitude hypoxia and genetic variants in the promoter region of the VEGFA gene in Sherpa highlanders.

Background: Sherpa highlanders demonstrate extraordinary tolerance to hypoxia at high altitudes, which may be achieved by mechanisms promoting microcirculatory blood flow and capillary density at high altitudes for restoring oxygen supply to tissues. Vascular endothelial growth factors (VEGFs) are important signaling proteins involved in vasculogenesis and angiogenesis which are stimulated by hypoxia. We hypothesize that the VEGF-A, the major member of the VEGF family, and the gene encoding VEGF-A (VEGFA) play a part in the adaptation to high-altitude hypoxia in Sherpa highlanders. Methods: Fifty-one Sherpa highlanders in Namche Bazaar village at a high altitude of 3,440 meters (m) above sea level and 76 non-Sherpa lowlanders in Kathmandu city at 1,300 m in Nepal were recruited for the study. Venous blood was sampled to obtain plasma and extract DNA from each subject. The plasma VEGF-A concentrations were measured and five single-nucleotide polymorphisms (SNPs, rs699947, rs833061, rs1570360, rs2010963, and rs3025039) in the VEGFA were genotyped. The VEGF-A levels and allelic frequencies of the SNPs were compared between the two populations. Results: A significant difference in oxygen saturation (SpO2) was observed between the two ethnic groups locating at different elevations (93.7 ± 0.2% in Sherpas at 3,440 m vs. 96.7 ± 0.2% in non-Sherpas at 1,300 m, P < 0.05). The plasma VEGF-A concentration in the Sherpas at high altitude was on the same level as that in the non-Sherpas at low altitude (262.8 ± 17.9 pg/ml vs. 266.8 ± 21.8 pg/ml, P = 0.88). This result suggested that the plasma VEGF-A concentration in Sherpa highlanders was stable despite a high-altitude hypoxic stimulus and that therefore the Sherpas exhibited a phenotype of blunted response to hypoxic stress. Moreover, the allele frequencies of the SNPs rs699947, rs833061, and rs2010963 in the promoter region of the VEGFA were different between the Sherpa highlanders and non-Sherpa lowlanders (corrected P values = 3.30 ×10-5, 4.95 ×10-4, and 1.19 ×10-7, respectively). Conclusions: Sherpa highlanders exhibited a blunted VEGF-A response to hypoxia at high altitudes, which was speculated to be associated with the distinctive genetic variations of the SNPs and haplotype in the promoter region of VEGFA in Sherpa highlanders.

Whole genome resequencing reveals signatures of rapid selection in a virus-affected commercial fishery.

Infectious diseases are recognized as one of the greatest global threats to biodiversity and ecosystem functioning. Consequently, there is a growing urgency to understand the speed at which adaptive phenotypes can evolve and spread in natural populations to inform future management. Here we provide evidence of rapid genomic changes in wild Australian blacklip abalone (Haliotis rubra) following a major population crash associated with an infectious disease. Genome scans on H. rubra were performed using pooled whole genome resequencing data from commercial fishing stocks varying in historical exposure to haliotid herpesvirus-1 (HaHV-1). Approximately 25,000 single nucleotide polymorphism loci associated with virus exposure were identified, many of which mapped to genes known to contribute to HaHV-1 immunity in the New Zealand paua (Haliotis iris) and herpesvirus response pathways in haliotids and other animal systems. These findings indicate genetic changes across a single generation in H. rubra fishing stocks decimated by HaHV-1, with stock recovery potentially determined by rapid evolutionary changes leading to virus resistance. This is a novel example of apparently rapid adaptation in natural populations of a nonmodel marine organism, highlighting the pace at which selection can potentially act to counter disease in wildlife communities.

Complementary genomic and epigenomic adaptation to environmental heterogeneity.

While adaptation is commonly thought to result from selection on DNA sequence-based variation, recent studies have highlighted an analogous epigenetic component as well. However, the relative roles of these mechanisms in facilitating population persistence under environmental heterogeneity remain unclear. To address the underlying genetic and epigenetic mechanisms and their relationship during environmental adaptation, we screened the genomes and epigenomes of nine global populations of a predominately sessile marine invasive tunicate, Botryllus schlosseri. We detected clear population differentiation at the genetic and epigenetic levels. Patterns of genetic and epigenetic structure were significantly influenced by local environmental variables. Among these variables, minimum annual sea surface temperature was identified as the top explanatory variable for both genetic and epigenetic variation. However, patterns of population structure driven by genetic and epigenetic variation were somewhat distinct, suggesting possible autonomy of epigenetic variation. We found both shared and specific genes and biological pathways among genetic and epigenetic loci associated with environmental factors, consistent with complementary and independent contributions of genetic and epigenetic variation to environmental adaptation in this system. Collectively, these mechanisms may facilitate population persistence under environmental change and sustain successful invasions across novel environments.

Rapid adaptation of Brachionus dorcas (Rotifera) to tetracycline antibiotic stress.

Although natural populations can rapidly adapt to selection pressures, the fitness consequences of selection are controversial. In this study, a selection experiment was conducted with replicate populations of Brachionus dorcas that were exposed to two sublethal concentrations (26.8 and 78.3 mg/L) of oxytetracycline (OTC), followed by two common garden experiments (population growth and life table experiments). During the 102-day (approximately 36 asexual generations) selection experiment, a markedly increased growth rate but a significantly decreased mictic ratio over time in the populations exposed to OTC when compared to the control populations suggested that the former adapted to the selection pressures and that a trade-off exists between asexual and sexual reproduction. The high and stable population growth rates after 90 days of OTC selection illustrate an example of evolutionary rescue. After 102 days of selection, OTC-selected populations showed higher population growth rates than the control populations when exposed to OTC, indicating significantly increased tolerance. OTC-selected populations showed a lower average growth rate, longer average generation time and life expectancy at hatching, and higher average net reproduction rate and proportion of mictic offspring than the control populations in the absence of OTC, which indicate that OTC selection results in two fitness costs and three fitness gains and that the effect of OTC selection on fitness differs with the measured fitness variables. Both the evolutionary potential of populations under the stress of higher concentrations of OTC and the fitness costs and gains of selection in the absence of OTC indicate that past exposures to pollutants cannot be neglected when evaluating the effects of current stressors on natural populations.

The genomic signatures of natural selection in admixed human populations.

Admixture has been a pervasive phenomenon in human history, extensively shaping the patterns of population genetic diversity. There is increasing evidence to suggest that admixture can also facilitate genetic adaptation to local environments, i.e., admixed populations acquire beneficial mutations from source populations, a process that we refer to as "adaptive admixture." However, the role of adaptive admixture in human evolution and the power to detect it remain poorly characterized. Here, we use extensive computer simulations to evaluate the power of several neutrality statistics to detect natural selection in the admixed population, assuming multiple admixture scenarios. We show that statistics based on admixture proportions, Fadm and LAD, show high power to detect mutations that are beneficial in the admixed population, whereas other statistics, including iHS and FST, falsely detect neutral mutations that have been selected in the source populations only. By combining Fadm and LAD into a single, powerful statistic, we scanned the genomes of 15 worldwide, admixed populations for signatures of adaptive admixture. We confirm that lactase persistence and resistance to malaria have been under adaptive admixture in West Africans and in Malagasy, North Africans, and South Asians, respectively. Our approach also uncovers other cases of adaptive admixture, including APOL1 in Fulani nomads and PKN2 in East Indonesians, involved in resistance to infection and metabolism, respectively. Collectively, our study provides evidence that adaptive admixture has occurred in human populations whose genetic history is characterized by periods of isolation and spatial expansions resulting in increased gene flow.

Persistence and genetic adaptation of Pseudomonas aeruginosa in patients with chronic obstructive pulmonary disease.

OBJECTIVES: It is unclear whether recurrent sputum culture with Pseudomonas aeruginosa from patients with chronic obstructive pulmonary disease (COPD) is caused by intermittent airway carriage by different P. aeruginosa lineages or persistent carriage by the same lineage, and whether lineages genetically adapt during carriage. METHODS: Whole-genome sequencing was performed for P. aeruginosa isolates sampled longitudinally from sputum cultures in patients with COPD who were enrolled in an ongoing randomized controlled trial (clinicaltrials.gov: NCT03262142). RESULTS: A total of 153 P. aeruginosa isolates were sequenced for 23 patients during 365 days of follow-up. Recurrent presence of P. aeruginosa was seen in 19 patients (83%) and was caused by persistence of the same clonal lineage in all but one patient. We identified 38 genes mutated in parallel in two or more lineages, suggesting positive selection for adaptive mutations. Mutational enrichment analysis revealed genes important in antibiotic resistance and chronic infections to be more frequently mutated. DISCUSSION: Recurrent P. aeruginosa was common and carried for a prolonged time after initial detection in the airways of patients with COPD. Recurrence was caused by persistence of the same clonal lineage and was associated with genetic adaptation. Trial data on possible clinical benefits of attempting antibiotic eradication of P. aeruginosa in COPD are warranted.