Signatures Of Selective Sweeps Research Articles

Population genomics of Plasmodium ovale species in sub-Saharan Africa.

Plasmodium ovale curtisi (Poc) and Plasmodium ovale wallikeri (Pow) are relapsing malaria parasites endemic to Africa and Asia that were previously thought to represent a single species. Amid increasing detection of ovale malaria in sub-Saharan Africa, we performed a population genomic study of both species across the continent. We conducted whole-genome sequencing of 25 isolates from Central and East Africa and analyzed them alongside 20 previously published African genomes. Isolates were predominantly monoclonal (43/45), with their genetic similarity aligning with geography. Pow showed lower average nucleotide diversity (1.8×10-4) across the genome compared to Poc (3.0×10-4) (p < 0.0001). Signatures of selective sweeps involving the dihydrofolate reductase gene were found in both species, as were signs of balancing selection at the merozoite surface protein 1 gene. Differences in the nucleotide diversity of Poc and Pow may reflect unique demographic history, even as similar selective forces facilitate their resilience to malaria control interventions.

Whole genome sequence analysis of population structure and insecticide resistance markers in Anopheles melas from the Bijagós Archipelago, Guinea-Bissau

BackgroundAnopheles melas is an understudied malaria vector with a potential role in malaria transmission on the Bijagós Archipelago of Guinea-Bissau. This study presents the first whole-genome sequencing and population genetic analysis for this species from the Bijagós. To our knowledge, this also represents the largest population genetic analysis using WGS data from non-pooled An. melas mosquitoes.MethodsWGS was conducted for 30 individual An. melas collected during the peak malaria transmission season in 2019 from six different islands on the Bijagós Archipelago. Bioinformatics tools were used to investigate the population structure and prevalence of insecticide resistance markers in this mosquito population.ResultsInsecticide resistance mutations associated with pyrethroid resistance in Anopheles gambiae s.s. from the Bijagós were absent in the An. melas population, and no signatures of selective sweeps were identified in insecticide resistance-associated genes. Analysis of structural variants identified a large duplication encompassing the cytochrome-P450 gene cyp9k1. Phylogenetic analysis using publicly available mitochondrial genomes indicated that An. melas from the Bijagós split into two phylogenetic groups because of differentiation on the mitochondrial genome attributed to the cytochrome C oxidase subunits COX I and COX II and the NADH dehydrogenase subunits 1, 4, 4L and 5.ConclusionsThis study identified an absence of insecticide-resistant SNPs common to An. gambiae in the An. melas population, but did identify structural variation over insecticide resistance-associated genes. Furthermore, this study presents novel insights into the population structure of this malaria vector using WGS analysis. Additional studies are required to further understand the role of this vector in malaria transmission.Graphical abstract

Genetic diversity during selective sweeps in non-recombining populations.

Selective sweeps, resulting from the spread of beneficial, neutral, or deleterious mutations through a population, shape patterns of genetic variation at linked neutral sites. While many theoretical, computational, and statistical advances have been made in understanding the genomic signatures of selective sweeps in recombining populations, substantially less is understood in populations with little/no recombination. We present a mathematical framework based on diffusion theory for obtaining the site frequency spectrum (SFS) at linked neutral sites immediately post and during the fixation of moderately or strongly beneficial mutations. We find that when a single hard sweep occurs, the SFS decays as 1/ x for low derived allele frequencies ( x ), similar to the neutral SFS at equilibrium, whereas at higher derived allele frequencies, it follows a 1/ x 2 power law. These power laws are universal in the sense that they are independent of the dominance and inbreeding coefficient, and also characterize the SFS during the sweep. Additionally, we find that the derived allele frequency where the SFS shifts from the 1/ x to 1/ x 2 law, is inversely proportional to the selection strength: thus under strong selection, the SFS follows the 1/ x 2 dependence for most allele frequencies, resembling a rapidly expanding neutral population. When clonal interference is pervasive, the SFS immediately post-fixation becomes U-shaped and is better explained by the equilibrium SFS of selected sites. Our results will be important in developing statistical methods to infer the timing and strength of recent selective sweeps in asexual populations, genomic regions that lack recombination, and clonally propagating tumor populations.

Genome structural variants shape adaptive success of an invasive urban malaria vector Anopheles stephensi.

Global changes are associated with the emergence of several invasive species. However, the genomic determinants of the adaptive success of an invasive species in a new environment remain poorly understood. Genomic structural variants (SVs), consisting of copy number variants, play an important role in adaptation. SVs often cause large adaptive shifts in ecologically important traits, which makes SVs compelling candidates for driving rapid adaptations to environmental changes, which is critical to invasive success. To address this problem, we investigated the role SVs play in the adaptive success of Anopheles stephensi , a primary vector of urban malaria in South Asia and an invasive malaria vector in several South Asian islands and Africa. We collected whole genome sequencing data from 115 mosquitoes from invasive island populations and four locations from mainland India, an ancestral range for the species. We identified 2,988 duplication copy number variants and 16,038 deletions in these strains, with ∼50% overlapping genes. SVs are enriched in genomic regions with signatures of selective sweeps in the mainland and invasive island populations, implying a putative adaptive role of SVs. Nearly all high-frequency SVs, including the candidate adaptive variants, in the invasive island populations are present on the mainland, suggesting a major contribution of existing variation to the success of the island populations. Among the candidate adaptive SVs, three duplications involving toxin-resistance genes evolved, likely due to the widespread application of insecticides in India since the 1950s. We also identify two SVs associated with the adaptation of An. stephensi larvae to brackish water in the island and two coastal mainland populations, where the mutations likely originated. Our results suggest that existing SVs play a vital role in the evolutionary success of An. stephensi in new environmental conditions.

Signatures of selective sweeps in continuous-space populations.

Selective sweeps describe the process by which an adaptive mutation arises and rapidly fixes in the population, thereby removing genetic variation in its genomic vicinity. The expected signatures of selective sweeps are relatively well understood in panmictic population models, yet natural populations often extend across larger geographic ranges where individuals are more likely to mate with those born nearby. To investigate how such spatial population structure can affect sweep dynamics and signatures, we simulated selective sweeps in populations inhabiting a two-dimensional continuous landscape. The maximum dispersal distance of offspring from their parents can be varied in our simulations from an essentially panmictic population to scenarios with increasingly limited dispersal. We find that in low-dispersal populations, adaptive mutations spread more slowly than in panmictic ones, while recombination becomes less effective at breaking up genetic linkage around the sweep locus. Together, these factors result in a trough of reduced genetic diversity around the sweep locus that looks very similar across dispersal rates. We also find that the site frequency spectrum around hard sweeps in low-dispersal populations becomes enriched for intermediate-frequency variants, making these sweeps appear softer than they are. Furthermore, haplotype heterozygosity at the sweep locus tends to be elevated in low-dispersal scenarios as compared to panmixia, contrary to what we observe in neutral scenarios without sweeps. The haplotype patterns generated by these hard sweeps in low-dispersal populations can resemble soft sweeps from standing genetic variation that arose from substantially older alleles. Our results highlight the need for better accounting for spatial population structure when making inferences about selective sweeps.

Determining the genetic basis of ginsenosides variation during ginseng domestication by evolutionary transcriptomics

Ginsenosides, the primary medicinally active substances in ginseng, are consistently used to treat a variety of diseases. Cultivated ginseng was domesticated from wild ginseng, for which the mode of domestication and the changes in the content of ginsenosides and their regulatory mechanisms during the domestication process were unclear. Here, transcriptome sequencing and metabolomic analysis were used to compare wild and cultivated ginsengs to gain insight into the genetic and ginsenosides content changes associated with domestication. Genetic structure analysis revealed that ginseng might have undergone independent domestication in China and Korea, respectively. A slightly reduced genetic diversity occurred in wild and cultivated ginsengs and ginseng was domesticated about 2000 years ago in China. 330 genes were identified as signatures of selective sweeps that occurred during ginseng domestication, including bZIP051, GARP222, HMGCS058, and UGTs225. The content of ginsenosides in leaves of cultivated ginseng was higher than leaves of wild ginseng, while substantial monomers content in wild ginseng was significantly higher than that in cultivated ginseng. The content of ginsenosides was higher in leaves than in roots and stems of both wild and cultivated ginsengs. Module-trait association analysis suggested abundant TFs, including MYB, AP2/ERF, NAC, and bZIP, were highly associated with expression level of ginsenosides synthetic genes and monomer content. The occurrence of alternative splicing events in genes was significantly reduced during ginseng domestication, with ginsenosides synthesis genes being particularly susceptible to alternative splicing. The variation of ginsenosides content was independently regulated by alternative splicing and transcription factors between cultivated and wild ginsengs. This study provided valuable genetic sources of ginsenosides that could contribute to ginseng breeding.

Whole-genome resequencing of major populations revealed domestication-related genes in yaks

BackgroundThe yak is a symbol of the Qinghai-Tibet Plateau and provides important basic resources for human life on the plateau. Domestic yaks have been subjected to strong artificial selection and environmental pressures over the long-term. Understanding the molecular mechanisms of phenotypic differences in yak populations can reveal key functional genes involved in the domestication process and improve genetic breeding.Material and methodHere, we re-sequenced 80 yaks (Maiwa, Yushu, and Huanhu populations) to identify single-nucleotide polymorphisms (SNPs) as genetic variants. After filtering and quality control, remaining SNPs were kept to identify the genome-wide regions of selective sweeps associated with domestic traits. The four methods (π, XPEHH, iHS, and XP-nSL) were used to detect the population genetic separation.ResultsBy comparing the differences in the population stratification, linkage disequilibrium decay rate, and characteristic selective sweep signals, we identified 203 putative selective regions of domestic traits, 45 of which were mapped to 27 known genes. They were clustered into 4 major GO biological process terms. All known genes were associated with seven major domestication traits, such as dwarfism (ANKRD28), milk (HECW1, HECW2, and OSBPL2), meat (SPATA5 and GRHL2), fertility (BTBD11 and ARFIP1), adaptation (NCKAP5, ANTXR1, LAMA5, OSBPL2, AOC2, and RYR2), growth (GRHL2, GRID2, SMARCAL1, and EPHB2), and the immune system (INPP5D and ADCYAP1R1).ConclusionsWe provided there is an obvious genetic different among domestic progress in these three yak populations. Our findings improve the understanding of the major genetic switches and domestic processes among yak populations.

Functional allele of a MATE gene selected during domestication modulates seed color in chickpea.

Seed color is one of the key target traits of domestication and artificial selection in chickpeas due to its implications on consumer preference and market value. The complex seed color trait has been well dissected in several crop species; however, the genetic mechanism underlying seed color variation in chickpea remains poorly understood. Here, we employed an integrated genomics strategy involving QTL mapping, high-density mapping, map-based cloning, association analysis, and molecular haplotyping in an inter-specific RIL mapping population, association panel, wild accessions, and introgression lines (ILs) of Cicer gene pool. This delineated a MATE gene, CaMATE23, encoding a Transparent Testa (TT) and its natural allele (8-bp insertion) and haplotype underlying a major QTL governing seed color on chickpea chromosome 4. Signatures of selective sweep and a strong purifying selection reflected that CaMATE23, especially its 8-bp insertion natural allelic variant, underwent selection during chickpea domestication. Functional investigations revealed that the 8-bp insertion containing the third cis-regulatory RY-motif element in the CaMATE23 promoter is critical for enhanced binding of CaFUSCA3 transcription factor, a key regulator of seed development and flavonoid biosynthesis, thereby affecting CaMATE23 expression and proanthocyanidin (PA) accumulation in the seed coat to impart varied seed color in chickpea. Consequently, overexpression of CaMATE23 in Arabidopsis tt12 mutant partially restored the seed color phenotype to brown pigmentation, ascertaining its functional role in PA accumulation in the seed coat. These findings shed new light on the seed color regulation and evolutionary history, and highlight the transcriptional regulation of CaMATE23 by CaFUSCA3 in modulating seed color in chickpea. The functionally relevant InDel variation, natural allele, and haplotype from CaMATE23 are vital for translational genomic research, including marker-assisted breeding, for developing chickpea cultivars with desirable seed color that appeal to consumers and meet global market demand.

Whole genome analyses reveal weak signatures of population structure and environmentally associated local adaptation in an important North American pollinator, the bumble bee Bombus vosnesenskii.

Studies of species that experience environmental heterogeneity across their distributions have become an important tool for understanding mechanisms of adaptation and predicting responses to climate change. We examine population structure, demographic history and environmentally associated genomic variation in Bombus vosnesenskii, a common bumble bee in the western USA, using whole genome resequencing of populations distributed across a broad range of latitudes and elevations. We find that B. vosnesenskii exhibits minimal population structure and weak isolation by distance, confirming results from previous studies using other molecular marker types. Similarly, demographic analyses with Sequentially Markovian Coalescent models suggest that minimal population structure may have persisted since the last interglacial period, with genomes from different parts of the species range showing similar historical effective population size trajectories and relatively small fluctuations through time. Redundancy analysis revealed a small amount of genomic variation explained by bioclimatic variables. Environmental association analysis with latent factor mixed modelling (LFMM2) identified few outlier loci that were sparsely distributed throughout the genome and although a few putative signatures of selective sweeps were identified, none encompassed particularly large numbers of loci. Some outlier loci were in genes with known regulatory relationships, suggesting the possibility of weak selection, although compared with other species examined with similar approaches, evidence for extensive local adaptation signatures in the genome was relatively weak. Overall, results indicate B. vosnesenskii is an example of a generalist with a high degree of flexibility in its environmental requirements that may ultimately benefit the species under periods of climate change.

Chromosome-level genome and population genomics provide novel insights into adaptive divergence in allopatric Eleutheronema tetradactylum

Understanding the adaptive ecological divergence provides important information for revealing biodiversity generation and maintenance. Adaptive ecology divergence in populations occurs in various environments and locations, but its genetic underpinnings remain elusive. We generated a chromosome-level genome of Eleutheronema tetradactylum (~582 Mb) and re-sequenced 50 allopatric E. tetradactylum in two independent environmental axes in China and Thailand Coastal waters as well as 11 cultured relatives. A low level of whole genome-wide diversity explained their decreased adaptive potential in the wild environment. Demographic analysis showed evidence of historically high abundance followed by a continuous distinct decline, plus signs of recent inbreeding and accumulation of deleterious mutations. Extensive signals of selective sweeps with signs of local adaptation to environmental differentiation between China and Thailand at genes related to thermal and salinity adaptation were discovered, which might be the driving factors of the geographical divergence of E. tetradactylum. Many genes and pathways subjected to strong selection under artificial breeding were associated with fatty acids and immunity (ELOVL6L, MAPK, p53/NF-kB), likely contributing to the eventual adaptation of artificial selective breeding. Our comprehensive study provided crucial genetic information for E. tetradactylum, with implications for the further conservation efforts of this threatened and ecologically valuable fish.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii.

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km2 region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

Population Genetic Analysis of Six Chinese Indigenous Pig Meta-Populations Based on Geographically Isolated Regions.

The diversification of indigenous pig breeds in China has resulted from multiple climate, topographic, and human cultural influences. The numerous indigenous pig breeds can be geographically divided into six meta-populations; however, their genetic relationships, contributions to genetic diversity, and genetic signatures remain unclear. Whole-genome SNP data for 613 indigenous pigs from the six Chinese meta-populations were obtained and analyzed. Population genetic analyses confirmed significant genetic differentiation and a moderate mixture among the Chinese indigenous pig meta-populations. The North China (NC) meta-population had the largest contribution to genetic and allelic diversity. Evidence from selective sweep signatures revealed that genes related to fat deposition and heat stress response (EPAS1, NFE2L2, VPS13A, SPRY1, PLA2G4A, and UBE3D) were potentially involved in adaptations to cold and heat. These findings from population genetic analyses provide a better understanding of indigenous pig characteristics in different environments and a theoretical basis for future work on the conservation and breeding of Chinese indigenous pigs.

Selective footprints and genes relevant to cold adaptation and other phenotypic traits are unscrambled in the genomes of divergently selected chicken breeds

BackgroundThe genomes of worldwide poultry breeds divergently selected for performance and other phenotypic traits may also be affected by, and formed due to, past and current admixture events. Adaptation to diverse environments, including acclimation to harsh climatic conditions, has also left selection footprints in breed genomes.ResultsUsing the Chicken 50K_CobbCons SNP chip, we genotyped four divergently selected breeds: two aboriginal, cold tolerant Ushanka and Orloff Mille Fleur, one egg-type Russian White subjected to artificial selection for cold tolerance, and one meat-type White Cornish. Signals of selective sweeps were determined in the studied breeds using three methods: (1) assessment of runs of homozygosity islands, (2) FST based population differential analysis, and (3) haplotype differentiation analysis. Genomic regions of true selection signatures were identified by two or more methods or in two or more breeds. In these regions, we detected 540 prioritized candidate genes supplemented them with those that occurred in one breed using one statistic and were suggested in other studies. Amongst them, SOX5, ME3, ZNF536, WWP1, RIPK2, OSGIN2, DECR1, TPO, PPARGC1A, BDNF, MSTN, and beta-keratin genes can be especially mentioned as candidates for cold adaptation. Epigenetic factors may be involved in regulating some of these important genes (e.g., TPO and BDNF).ConclusionBased on a genome-wide scan, our findings can help dissect the genetic architecture underlying various phenotypic traits in chicken breeds. These include genes representing the sine qua non for adaptation to harsh environments. Cold tolerance in acclimated chicken breeds may be developed following one of few specific gene expression mechanisms or more than one overlapping response known in cold-exposed individuals, and this warrants further investigation.

Evolutionary responses of a reef-building coral to climate change at the end of the last glacial maximum.

Climate change threatens the survival of coral reefs on a global scale, primarily through mass bleaching and mortality as a result of marine heatwaves. While these short-term effects are clear, predicting the fate of coral reefs over the coming century is a major challenge. One way to understand the longer-term effect of rapid climate change is to examine the response of coral populations to past climate shifts. Coastal and shallow-water marine ecosystems such as coral reefs have been reshaped many times by sea-level changes during the Pleistocene, yet few studies have directly linked this with its consequences on population demographics, dispersal, and adaptation. Here we use powerful analytical techniques, afforded by haplotype-phased whole-genomes, to establish such links for the reef-building coral, Acropora digitifera. We show that three genetically distinct populations are present in northwestern Australia, and that their rapid divergence since the last glacial maximum (LGM) can be explained by a combination of founder-effects and restricted gene flow. Signatures of selective sweeps, too strong to be explained by demographic history, are present in all three populations and overlap with genes that show different patterns of functional enrichment between inshore and offshore habitats. In contrast to rapid divergence in the host, we find that photosymbiont communities are largely undifferentiated between corals from all three locations, spanning almost 1000 km, indicating that selection on host genes, and not acquisition of novel symbionts, has been the primary driver of adaptation for this species in northwestern Australia.

Comparative genomics of Sarcoptes scabiei provide new insights into adaptation to permanent parasitism and within-host species divergence.

Sarcoptes scabiei is the causative agent of a highly contagious skin disease in humans and more than 100 mammals. Here, we reported the first chromosome-level reference genome of S. scabiei isolated from rabbits, with a contig N50 size of 5.92Mb, a total assembled length of 57.30Mb, ∼12.65% repetitive sequences and 9333 predicted protein-coding genes. The phylogenetic tree based on 1338 shared high-confidence single-copy orthologous genes estimated that the mammalian ectoparasite S. scabiei and the plant-feeding mite Tetranychus urticae separated ∼340 million years ago. Both neighbour-joining tree and principal component analysis of 20 mite populations isolated from four hosts (humans, pigs, dogs and rabbits) distributed in three countries (China, Australia and the USA) consistently supported genetic subdivisions according to host species rather than geographical location. The demographic history of S. scabiei reconstructed by multiple sequentially Markovian coalescent analysis suggested that S. scabiei isolated from rabbits, humans, dogs and pigs diverged ∼5000 years ago. Investigation of the homeobox (Hox) genes revealed that S. scabiei contains 8 of 10 canonical Hox genes that are present in the arthropod ancestor, and the absence of the Abd-A gene may correlate with the long gap between their front and back legs. Comparative genomics demonstrated that genes specific to scabies mites were mainly enriched in nutrition digestive systems, whereas genes in the families that involved detoxification (cytochrome P450, carboxyl/cholinesterases and the ATP-binding cassette transporter C group) were extremely contracted compared with that of other mites analysed in this study. Selective sweep analysis of mite populations from various pairs of two out of the four host species revealed that the strongest selective sweep signals were mainly enriched in cysteine-type peptidase activity and apoptosis. The results provided clues for the mechanisms of S. scabiei adaptation to a permanent parasitic lifestyle and knowledge that would enable further control of this highly contagious skin disease.

Identifying Candidate Genes for Short Gestation Length Trait in Chinese Qingping Pigs by Whole-Genome Resequencing and RNA Sequencing.

Gestation length is a complex polygenic trait that affects pig fetal development. The Qingping (QP) pig, a Chinese native black pig breed, is characterized by short gestation length. However, the genetic architecture of short gestation length is still not clear. The present study aimed to explore the genetic architecture of short gestation length in QP pigs. In this study, selective sweep analyses were performed to detect selective sweep signatures for short gestation length traits between 100 QP pigs and 219 pigs from 15 other breeds. In addition, differentially expressed genes for the short gestation length between QP pigs and Large White pigs were detected by RNA sequencing. Comparing candidate genes from these methods with known genes for preterm birth in the database, we obtained 111 candidate genes that were known preterm birth genes. Prioritizing other candidate genes, 839 novel prioritized candidate genes were found to have significant functional similarity to preterm birth genes. In particular, we highlighted EGFR, which was the most prioritized novel candidate relative to preterm birth genes. Experimental validations in placental and porcine trophectoderm cells suggest that EGFR is highly expressed in the QP pigs with short gestation length and could regulate the NF-κΒ pathway and downstream expression of PTGS2. These findings comprehensively identified candidate genes for short gestation length trait at the genomic and transcriptomic levels. These candidate genes provide an important new resource for further investigation and genetic improvement of gestation length.

Strong neutral sweeps occurring during a population contraction.

A strong reduction in diversity around a specific locus is often interpreted as a recent rapid fixation of a positively selected allele, a phenomenon called a selective sweep. Rapid fixation of neutral variants can however lead to a similar reduction in local diversity, especially when the population experiences changes in population size, e.g. bottlenecks or range expansions. The fact that demographic processes can lead to signals of nucleotide diversity very similar to signals of selective sweeps is at the core of an ongoing discussion about the roles of demography and natural selection in shaping patterns of neutral variation. Here, we quantitatively investigate the shape of such neutral valleys of diversity under a simple model of a single population size change, and we compare it to signals of a selective sweep. We analytically describe the expected shape of such "neutral sweeps" and show that selective sweep valleys of diversity are, for the same fixation time, wider than neutral valleys. On the other hand, it is always possible to parametrize our model to find a neutral valley that has the same width as a given selected valley. Our findings provide further insight into how simple demographic models can create valleys of genetic diversity similar to those attributed to positive selection.

Selective Sweeps and Polygenic Adaptation Drive Local Adaptation along Moisture and Temperature Gradients in Natural Populations of Coast Redwood and Giant Sequoia.

Dissecting the genomic basis of local adaptation is a major goal in evolutionary biology and conservation science. Rapid changes in the climate pose significant challenges to the survival of natural populations, and the genomic basis of long-generation plant species is still poorly understood. Here, we investigated genome-wide climate adaptation in giant sequoia and coast redwood, two iconic and ecologically important tree species. We used a combination of univariate and multivariate genotype–environment association methods and a selective sweep analysis using non-overlapping sliding windows. We identified genomic regions of potential adaptive importance, showing strong associations to moisture variables and mean annual temperature. Our results found a complex architecture of climate adaptation in the species, with genomic regions showing signatures of selective sweeps, polygenic adaptation, or a combination of both, suggesting recent or ongoing climate adaptation along moisture and temperature gradients in giant sequoia and coast redwood. The results of this study provide a first step toward identifying genomic regions of adaptive significance in the species and will provide information to guide management and conservation strategies that seek to maximize adaptive potential in the face of climate change.

Linked selection shapes the landscape of genomic variation in three oak species.

Natural selection shapes genome-wide patterns of diversity within species and divergence between species. However, quantifying the efficacy of selection and elucidating the relative importance of different types of selection in shaping genomic variation remain challenging. We sequenced whole genomes of 101 individuals of three closely related oak species to track the divergence history, and to dissect the impacts of selective sweeps and background selection on patterns of genomic variation. We estimated that the three species diverged around the late Neogene and experienced a bottleneck during the Pleistocene. We detected genomic regions with elevated relative differentiation ('FST -islands'). Population genetic inferences from the site frequency spectrum and ancestral recombination graph indicated that FST -islands were formed by selective sweeps. We also found extensive positive selection; the fixation of adaptive mutations and reduction neutral diversity around substitutions generated a signature of selective sweeps. Prevalent negative selection and background selection have reduced genetic diversity in both genic and intergenic regions, and contributed substantially to the baseline variation in genetic diversity. Our results demonstrate the importance of linked selection in shaping genomic variation, and illustrate how the extent and strength of different selection models vary across the genome.

Haplotype Variations and Evolutionary Analysis of the Granule-Bound Starch Synthase I Gene in the Korean World Rice Collection

Granule-bound starch synthase I (GBSSI) is responsible for Waxy gene encoding the, which is involved in the amylose synthesis step of starch biosynthesis. We investigated the genotypic and haplotypic variations of GBSSI (Os06g0133000) gene, including its evolutionary relatedness in the nucleotide sequence level using single-nucleotide polymorphisms (SNPs), indels, and structural variations (SVs) from 475 Korean World Rice Collection (KRICE_CORE), which comprised 54 wild rice and 421 cultivated represented by 6 ecotypes (temperate japonica, indica, tropical japonica, aus, aromatic, and admixture) or in another way by 3 varietal types (landrace, weedy, and bred). The results revealed that 27 of 59 haplotypes indicated a total of 12 functional SNPs (fSNPs), identifying 9 novel fSNPs. According to the identified novel fSNPs, we classified the entire rice collection into three groups: cultivated, wild, and mixed (cultivated and wild) rice. Five novel fSNPs were localized in wild rice: four G/A fSNPs in exons 2, 9, and 12 and one T/C fSNP in exon 13. We also identified the three previously reported fSNPs, namely, a G/A fSNP (exon 4), an A/C fSNP (exon 6), and a C/T fSNP (exon 10), which were observed only in cultivated rice, whereas an A/G fSNP (exon 4) was observed exclusively in wild rice. All-against-all comparison of four varietal types or six ecotypes of cultivated rice with wild rice showed that the GBSSI diversity was higher only in wild rice (π = 0.0056). The diversity reduction in cultivated rice can be useful to encompass the origin of this gene GBSSI during its evolution. Significant deviations of positive (wild and indica under balancing selection) and negative (temperate and tropical japonica under purifying selection) Tajima's D values from a neutral model can be informative about the selective sweeps of GBSSI genome insights. Despite the estimation of the differences in population structure and principal component analysis (PCA) between wild and subdivided cultivated subgroups, an inbreeding effect was quantified by FST statistic, signifying the genetic relatedness of GBSSI. Our findings of a novel wild fSNPS can be applicable for future breeding of waxy rice varieties. Furthermore, the signatures of selective sweep can also be of informative into further deeper insights during domestication.