Recent Selective Sweep Research Articles

Low nucleotide diversity of the Plasmodium falciparum AP2-EXP2 gene among clinical samples from Ghana.

PfAP2-EXP2 is located within chromosome 6 of Plasmodium falciparum recently identified to be undergoing an extensive selective sweep in West African isolates. The gene encoding this transcription factor, PfAP2-EXP2, is essential and thus likely subject to purifying selection that limits variants in the parasite population despite its genomic location. 72 Plasmodium falciparum field samples and 801 clinical sequences from the Pf6 MalariaGEN dataset of Ghanaian origin, were integrated and analysed. A total of 14 single nucleotide variants of which 5 were missense variants, were identified after quality checks and filtering. Except for one, all identified variants were rare among the clinical samples obtained in this study (Minor allelic frequency < 0.01). Further results revealed a considerably low dN/dS value (0.208) suggesting the presence of purifying selection. Further, all the mutant amino acids were wildtype residues in AP2-EXP2 orthologous proteins-tentatively suggesting a genus-level conservation of amino acid residues. Computational analysis and predictions corroborated these findings. Despite the recent extensive selective sweep within chromosome 6 of West African isolates, PfAP2-EXP2 of Ghanaian origin exhibits low nucleotide diversity and very low dN/dS consistent with purifying selection acting to maintain the function of an essential gene. The conservation of AP2-EXP2 is an important factor that makes it a potential drug target.

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.

Two decades of molecular surveillance in Senegal reveal rapid changes in known drug resistance mutations over time

BackgroundDrug resistance in Plasmodium falciparum is a major threat to malaria control efforts. Pathogen genomic surveillance could be invaluable for monitoring current and emerging parasite drug resistance.MethodsData from two decades (2000–2020) of continuous molecular surveillance of P. falciparum parasites from Senegal were retrospectively examined to assess historical changes in malaria drug resistance mutations. Several known drug resistance markers and their surrounding haplotypes were profiled using a combination of single nucleotide polymorphism (SNP) molecular surveillance and whole genome sequence based population genomics.ResultsThis dataset was used to track temporal changes in drug resistance markers whose timing correspond to historically significant events such as the withdrawal of chloroquine (CQ) and the introduction of sulfadoxine-pyrimethamine (SP) in 2003. Changes in the mutation frequency at Pfcrt K76T and Pfdhps A437G coinciding with the 2014 introduction of seasonal malaria chemoprevention (SMC) in Senegal were observed. In 2014, the frequency of Pfcrt K76T increased while the frequency of Pfdhps A437G declined. Haplotype-based analyses of Pfcrt K76T showed that this rapid increase was due to a recent selective sweep that started after 2014.Discussion (Conclusion)The rapid increase in Pfcrt K76T is troubling and could be a sign of emerging amodiaquine (AQ) resistance in Senegal. Emerging AQ resistance may threaten the future clinical efficacy of artesunate-amodiaquine (ASAQ) and AQ-dependent SMC chemoprevention. These results highlight the potential of molecular surveillance for detecting rapid changes in parasite populations and stress the need to monitor the effectiveness of AQ as a partner drug for artemisinin-based combination therapy (ACT) and for chemoprevention.

Subterranean termites raise the alarm when their anti-fungal weapon falters.

Termicin is an anti-fungal defensin that is disseminated from termite salivary glands. The peptide appears to be critical for the elimination with mutual grooming (allogrooming) of pathogenic spores (conidia) that have attached to the insect cuticle. There has been a recent selective sweep for an advantageous variant of this peptide in the subterranean termite Reticulitermes flavipes. We tested the anti-mycotic activity of a recombinant termicin corresponding with this variant against the conidia of different Metarhizium fungal isolates from soil close to foraging R. flavipes workers. Termicin was most effective against isolates that had previously been shown to elicit a relatively weak alarm response, as indicated by brief bouts of rapid longitudinal oscillatory movement (LOM). These isolates that elicited weak alarm were also the deadliest apparently because the survival of termites exposed to the fungus depends on a strong social immune response (LOMs and allogrooming). The selective pressure for a single termicin variant may have been driven by the most dangerous isolates that elicit a weak behavioral response. The correlation between termicin anti-fungal activity and LOM suggests that pathogen-associated molecular patterns that affect termite recognition of conidial contamination and the onset of elevated allogrooming also affect the vulnerability of conidia to the disruption of their cell membranes by termicin.

Two-speed genomes of Epichloe fungal pathogens show contrasting signatures of selection between species and across populations.

Antagonistic selection between pathogens and their hosts can drive rapid evolutionary change and leave distinct molecular footprints of past and ongoing selection in the genomes of the interacting species. Despite an increasing availability of tools able to identify signatures of selection, the genetic mechanisms underlying coevolutionary interactions and the specific genes involved are still poorly understood, especially in heterogeneous natural environments. We searched the genomes of two species of Epichloe plant pathogen for evidence of recent selection. The Epichloe genus includes highly host-specific species that can sterilize their grass hosts. We performed selection scans using genome-wide SNP data from seven natural populations of two co-occurring Epichloe sibling species specialized on different hosts. We found evidence of recent (and ongoing) selective sweeps across the genome in both species. However, selective sweeps were more abundant in the species with a larger effective population size. Sweep regions often overlapped with highly polymorphic AT-rich regions supporting the role of these genome compartments in adaptive evolution. Although most loci under selection were specific to individual populations, we could also identify several candidate genes targeted by selection in sweep regions shared among populations. The genes encoded small secreted proteins typical of fungal effectors and cell wall-degrading enzymes. By investigating the genomic signatures of selection across multiple populations and species, this study contributes to our understanding of complex adaptive processes in natural plant pathogen systems.

Do chromosome rearrangements fix by genetic drift or natural selection? Insights from Brenthis butterflies.

Large-scale chromosome rearrangements, such as fissions and fusions, are a common feature of eukaryote evolution. They can have considerable influence on the evolution of populations, yet it remains unclear exactly how rearrangements become established and eventually fix. Rearrangements could fix by genetic drift if they are weakly deleterious or neutral, or they may instead be favoured by positive natural selection. Here, we compare genome assemblies of three closely related Brenthis butterfly species and characterize a complex history of fission and fusion rearrangements. An inferred demographic history of these species suggests that rearrangements became fixed in populations with large long-term effective size (Ne ), consistent with rearrangements being selectively neutral or only very weakly underdominant. Using a recently developed analytic framework for characterizing hard selective sweeps, we find that chromosome fusions are not enriched for evidence of past sweeps compared to other regions of the genome. Nonetheless, we do infer a strong and recent selective sweep around one chromosome fusion in the B. daphne genome. Our results suggest that rearrangements in these species likely have weak absolute fitness effects and fix by genetic drift. However, one putative selective sweep raises the possibility that natural selection may sometimes play a role in the fixation of chromosome fusions.

Phylogeography of the Asian green mussel, Perna viridis: pronounced regional genetic differentiation is explained by major freshwater input into the Bay of Bengal

The Asian green mussel (AGM), Perna viridis, is a mussel of ecological and economic importance throughout much of the tropical and subtropical regions of Asia, from China and Korea in the east to Oman on the Arabian Peninsula in the west. We collated published DNA sequence data and then analysed 467 bp of cytochrome c oxidase I gene (COI) sequence variation from 649 mussels (576 from SE Asia and 73 from India) to test for regional differentiation. Analysis of molecular variation, a haplotype network and a neighbour-joining tree all revealed significant differentiation between mussels from India and those from SE Asia. We observed a greater number of segregating sites (haplotypes) and private haplotypes than expected for Indian mussels compared to SE Asian mussels, based on the respective sample sizes. Tests of neutrality revealed population expansion or a recent selective sweep for only the mussels from India. We suggest that the differentiation of the Indian from the SE Asian mussels is explained by the pronounced seasonal input of freshwater into the northern Bay of Bengal area, which acts as a (semi-permeable) barrier to gene flow between the two regions (i.e., east versus west of the Bay). This suggestion is based on the oceanography of the region and is consistent with the biogeography and reports of genetic breaks in a range of taxa in this region. Further targetted sampling of AGMs from the east coast of India, Bangladesh and Myanmar is required to test this hypothesis, and additional sampling of AGMs from west of India (e.g., Pakistan, Iran, the Arabian Peninsula) will also be informative. Finally, this study is based on the analysis of a single marker (COI) because this is all that is presently available from published sources but we note that the application of new molecular markers such as single nucleotide polymorphisms to newly collected mussel samples will greatly advance our understanding of the AGM genetic discontinuity and its age, and help test its geographical location and its nature (e.g., a steep/shallow cline or a stepped cline).

Whole-genome Analyses Reveal Past Population Fluctuations and Low Genetic Diversities of the North Pacific Albatrosses.

Throughout the Plio-Pleistocene, climate change has impacted tropical marine ecosystems substantially, with even more severe impacts predicted in the Anthropocene. Although many studies have clarified demographic histories of seabirds in polar regions, the history of keystone seabirds of the tropics is unclear, despite the prominence of albatrosses (Diomedeidae, Procellariiformes) as the largest and most threatened group of oceanic seabirds. To understand the impact of climate change on tropical albatrosses, we investigated the evolutionary and demographic histories of all four North Pacific albatrosses and their prey using whole-genome analyses. We report a striking concordance in demographic histories among the four species, with a notable dip in effective population size at the beginning of the Pleistocene and a population expansion in the Last Glacial Period when sea levels were low, which resulted in increased potential coastal breeding sites. Abundance of the black-footed albatross dropped again during the Last Glacial Maximum, potentially linked to climate-driven loss of breeding sites and concordant genome-derived decreases in its major prey. We find very low genome-wide (π < 0.001) and adaptative genetic diversities across the albatrosses, with genes of the major histocompatibility complex close to monomorphic. We also identify recent selective sweeps at genes associated with hyperosmotic adaptation, longevity, and cognition and memory. Our study has shed light on the evolutionary and demographic histories of the largest tropical oceanic seabirds and provides evidence for their large population fluctuations and alarmingly low genetic diversities.

The expansion of agriculture has shaped the recent evolutionary history of a specialized squash pollinator

The expansion of agriculture is responsible for the mass conversion of biologically diverse natural environments into managed agroecosystems dominated by a handful of genetically homogeneous crop species. Agricultural ecosystems typically have very different abiotic and ecological conditions from those they replaced and create potential niches for those species that are able to exploit the abundant resources offered by crop plants. While there are well-studied examples of crop pests that have adapted into novel agricultural niches, the impact of agricultural intensification on the evolution of crop mutualists such as pollinators is poorly understood. We combined genealogical inference from genomic data with archaeological records to demonstrate that the Holocene demographic history of a wild specialist pollinator of Cucurbita (pumpkins, squashes, and gourds) has been profoundly impacted by the history of agricultural expansion in North America. Populations of the squash bee Eucera pruinosa experienced rapid growth in areas where agriculture intensified within the past 1,000 y, suggesting that the cultivation of Cucurbita in North America has increased the amount of floral resources available to these bees. In addition, we found that roughly 20% of this bee species' genome shows signatures of recent selective sweeps. These signatures are overwhelmingly concentrated in populations from eastern North America where squash bees were historically able to colonize novel environments due to human cultivation of Cucurbita pepo and now exclusively inhabit agricultural niches. These results suggest that the widespread cultivation of crops can prompt adaptation in wild pollinators through the distinct ecological conditions imposed by agricultural environments.

Population Genomics of Plasmodium falciparum and Malaria Control Implications in Abidjan (Cote d’ivoire)

Introduction: The onset of Plasmodium falciparum (P. falciparum) resistance to antimalarial drugs requires the careful surveillance of African parasite populations. Genomic tools are implemented to detect evolutionary changes that could impact malaria control and elimination strategies. Here, we evaluate the genome-wide pattern of selection and sequence variation in P. falciparum populations in Abidjan, Côte d’Ivoire. Methods: The study was conducted in three localities of Abidjan from 2013 to 2014. We collected 70 blood samples after receiving written informed consent from patients above two years of age. After extracting P. falciparum and human DNA from isolates, we performed whole-genome sequencing and used population genomics approaches to investigate the genetic diversity and complexity of infections and identify loci under positive directional selection. Results: We observed an excess of rare variants in the population, showing a clear mutation process in the isolates. Moderate Fst estimates (0.3) was detected for surfin, an immune invasion gene family. Seven iHS regions that had at least two SNPs with a score &gt; 3.2 were identified. These regions code for genes that have been under strong directional selection. Two of these genes were the chloroquine resistance transporter (crt) on chromosome 7 and the dihydropteroate reductase (dhps) on chromosome 8. Our analyses showed that a recent selective sweep occurred for the erythrocyte membrane protein (Pfemp1). Conclusion: Our analyses identified genes under selective drug pressure and balancing selection on protective immune-specific genes. These findings demonstrate the effectiveness of genomics analyses to follow the evolution of malaria parasites and adopt appropriate strategies to eliminate malaria in Côte d’Ivoire.

Strong, Recent Selective Sweeps Reshape Genetic Diversity in Freshwater Bivalve Megalonaias nervosa.

Freshwater Unionid bivalves have recently faced ecological upheaval through pollution, barriers to dispersal, harvesting, and changes in fish-host prevalence. Currently, over 70% of species in North America are threatened, endangered or extinct. To characterize the genetic response to recent selective pressures, we collected population genetic data for one successful bivalve species, Megalonaias nervosa. We identify megabase-sized regions that are nearly monomorphic across the population, signals of strong, recent selection reshaping diversity across 73 Mb total. These signatures of selection are greater than is commonly seen in population genetic models. We observe 102 duplicate genes with high dN/dS on terminal branches among regions with sweeps, suggesting that gene duplication is a causative mechanism of recent adaptation in M. nervosa. Genes in sweeps reflect functional classes important for Unionid survival, including anticoagulation genes important for fish host parasitization, detox genes, mitochondria management, and shell formation. We identify sweeps in regions with no known functional impacts, suggesting mechanisms of adaptation that deserve greater attention in future work on species survival. In contrast, polymorphic transposable elements (TEs) appear to be detrimental and underrepresented among regions with sweeps. TE site frequency spectra are skewed toward singleton variants, and TEs among regions with sweeps are present at low frequency. Our work suggests that duplicate genes are an essential source of genetic novelty that has helped this species succeed in environments where others have struggled. These results suggest that gene duplications deserve greater attention in non-model population genomics, especially in species that have recently faced sudden environmental challenges.

Multilocus phylogeography, population genetics and niche evolution of Australian brown and black-tailed treecreepers (Aves:Climacteris)

AbstractThe Carpentarian barrier across north-eastern Australia is a major biogeographic barrier and a generator of biodiversity within the Australian Monsoonal Tropics. Here we present a continent-wide analysis of mitochondrial (control region) and autosomal (14 anonymous loci) sequence and indel variation and niche modelling of brown and black-tailed treecreepers (Climacteris picumnus and Climacteris melanurus), a clade with a classic distribution on either side of the Carpentarian barrier. mtDNA control region sequences exhibited reciprocal monophyly and strong differentiation (Fst = 0.91), and revealed a signature of a recent selective sweep in C. picumnus. A variety of tests support an isolation-with-migration model of divergence, albeit with low levels of gene flow across the Carpentarian barrier and a divergence time between species of ~1.7–2.8 Mya. Palaeoecological niche models show that both range size as measured by available habitat and estimated historical population sizes of both species declined in the past ~600 kyr and that the area of interspecific range overlap was never historically large, perhaps decreasing opportunities for extensive gene flow. The relatively long divergence time and low opportunity for gene flow may have facilitated speciation more so than in other co-distributed bird taxa across the Australian Monsoonal Tropics.

Vectorial capacity and TEP1 genotypes of Anopheles gambiae sensu lato mosquitoes on the Kenyan coast

BackgroundMalaria remains one of the most important infectious diseases in sub-Saharan Africa, responsible for approximately 228 million cases and 602,000 deaths in 2020. In this region, malaria transmission is driven mainly by mosquitoes of the Anopheles gambiae and, more recently, Anopheles funestus complex. The gains made in malaria control are threatened by insecticide resistance and behavioural plasticity among these vectors. This, therefore, calls for the development of alternative approaches such as malaria transmission-blocking vaccines or gene drive systems. The thioester-containing protein 1 (TEP1) gene, which mediates the killing of Plasmodium falciparum in the mosquito midgut, has recently been identified as a promising target for gene drive systems. Here we investigated the frequency and distribution of TEP1 alleles in wild-caught malaria vectors on the Kenyan coast.MethodsMosquitoes were collected using CDC light traps both indoors and outdoors from 20 houses in Garithe village, along the Kenyan coast. The mosquitoes were dissected, and the different parts were used to determine their species, blood meal source, and sporozoite status. The data were analysed and visualised using the R (v 4.0.1) and STATA (v 17.0).ResultsA total of 18,802 mosquitoes were collected, consisting of 77.8% (n = 14,631) Culex spp., 21.4% (n = 4026) An. gambiae sensu lato, 0.4% (n = 67) An. funestus, and 0.4% (n = 78) other Anopheles (An. coustani, An. pharoensis, and An. pretoriensis). Mosquitoes collected were predominantly exophilic, with the outdoor catches being higher across all the species: Culex spp. 93% (IRR = 11.6, 95% Cl [5.9–22.9] P < 0.001), An. gambiae s.l. 92% (IRR = 7.2, 95% Cl [3.6–14.5]; P < 0.001), An. funestus 91% (IRR = 10.3, 95% Cl [3.3–32.3]; P < 0.001). A subset of randomly selected An. gambiae s.l. (n = 518) was identified by polymerase chain reaction (PCR), among which 77.2% were An. merus, 22% were An. arabiensis, and the rest were not identified. We were also keen on identifying and describing the TEP1 genotypes of these mosquitoes, especially the *R3/R3 allele that was identified recently in the study area. We identified the following genotypes among An. merus: *R2/R2, *R3/R3, *R3/S2, *S1/S1, and *S2/S2. Among An. arabiensis, we identified *R2/R2, *S1/S1, and *S2/S2. Tests on haplotype diversity showed that the most diverse allele was TEP1*S1, followed by TEP1*R2. Tajima’s D values were positive for TEP1*S1, indicating that there is a balancing selection, negative for TEP1*R2, indicating there is a recent selective sweep, and as for TEP1*R3, there was no evidence of selection. Phylogenetic analysis showed two distinct clades: refractory and susceptible alleles.ConclusionsWe find that the malaria vectors An. gambiae s.l. and An. funestus are predominantly exophilic. TEP1 genotyping for An. merus revealed five allelic combinations, namely *R2/R2, *R3/R3, *R3/S2, *S1/S1 and *S2/S2, while in An. arabiensis we only identified three allelic combinations: *R2/R2, *S1/S1, and *S2/S2. The TEP1*R3 allele was restricted to only An. merus among these sympatric mosquito species, and we find that there is no evidence of recombination or selection in this allele.Graphical

Diverse selection pressures shaping the genetic architecture of behçet disease susceptibility.

Behçet disease (BD) is a polygenic, multifactorial, multisystem inflammatory condition with unknown etiology. Global distribution of BD is geographically structured, highest prevalence observed among East Asian, Middle Eastern, and Mediterranean populations. Although adaptive selection on a few BD susceptibility loci is speculated, a thorough evolutionary analysis on the genetic architecture of BD is lacking. We aimed to understand whether increased BD risk in the human populations with high prevalence is due to past selection on BD associated genes. We performed population genetics analyses with East Asian (high BD prevalence), European (low/very low BD prevalence), and African (very low/no BD prevalence) populations. Comparison of ancestral and derived alleles’ frequencies versus their reported susceptible or protective effect on BD showed both derived and ancestral alleles are associated with increased BD risk. Variants showing higher risk to and more significant association with BD had smaller allele frequency differences, and showed less population differentiation compared to variants that showed smaller risk and less significant association with BD. Results suggest BD alleles are not unique to East Asians but are also found in other world populations at appreciable frequencies, and argue against selection favoring these variants only in populations with high BD prevalence. BD associated gene analyses showed similar evolutionary histories driven by neutral processes for many genes or balancing selection for HLA (Human Leukocyte Antigen) genes in all three populations studied. However, nucleotide diversity in several HLA region genes was much higher in East Asians suggesting selection for high nucleotide and haplotype diversity in East Asians. Recent selective sweep for genes involved in antigen recognition, peptide processing, immune and cellular differentiation regulation was observed only in East Asians. We conclude that the evolutionary processes shaping the genetic diversity in BD risk genes are diverse, and elucidating the underlying specific selection mechanisms is complex. Several of the genes examined in this study are risk factors (such as ERAP1, IL23R, HLA-G) for other inflammatory diseases. Thus, our conclusions are not only limited to BD but may have broader implications for other inflammatory diseases.

Identification of a rapidly-spreading triple mutant for high-level metabolic insecticide resistance in Anopheles gambiae provides a real-time molecular diagnostic for antimalarial intervention deployment.

Studies of insecticide resistance provide insights into the capacity of populations to show rapid evolutionary responses to contemporary selection. Malaria control remains heavily dependent on pyrethroid insecticides, primarily in long lasting insecticidal nets (LLINs). Resistance in the major malaria vectors has increased in concert with the expansion of LLIN distributions. Identifying genetic mechanisms underlying high-level resistance is crucial for the development and deployment of resistance-breaking tools. Using the Anopheles gambiae 1000 genomes (Ag1000g) data we identified a very recent selective sweep in mosquitoes from Uganda which localized to a cluster of cytochrome P450 genes. Further interrogation revealed a haplotype involving a trio of mutations, a nonsynonymous point mutation in Cyp6p4 (I236M), an upstream insertion of a partial Zanzibar-like transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae from the Kenya-Uganda border, with stepwise replacement of the double-mutant (Zanzibar-like TE and Cyp6p4-236 M) with the triple-mutant haplotype (including Cyp6aa1 duplication), which has spread into the Democratic Republic of Congo and Tanzania. The triple-mutant haplotype is strongly associated with increased expression of genes able to metabolize pyrethroids and is strongly predictive of resistance to pyrethroids most notably deltamethrin. Importantly, there was increased mortality in mosquitoes carrying the triple-mutation when exposed to nets cotreated with the synergist piperonyl butoxide (PBO). Frequencies of the triple-mutant haplotype remain spatially variable within countries, suggesting an effective marker system to guide deployment decisions for limited supplies of PBO-pyrethroid cotreated LLINs across African countries.

A comparative study of population genetic structure reveals patterns consistent with selection at functional microsatellites in common sunflower.

Microsatellites, also known as short tandem repeats (STRs), have long been considered non-functional, neutrally evolving regions of the genome. Recent findings suggest that they can function as drivers of rapid adaptive evolution. Previous work on the common sunflower identified 479 transcribed microsatellites where allele length significantly correlates with gene expression (eSTRs) in a stepwise manner. Here, a population genetic approach is used to test whether eSTR allele length variation is under selection. Genotypic variation among and within populations at 13 eSTRs was compared with that at 19 anonymous microsatellites in 672 individuals from 17 natural populations of sunflower from across a cline running from Saskatchewan to Oklahoma (distance of approximately 1600km). Expected heterozygosity, allelic richness, and allelic diversity were significantly lower at eSTRs, a pattern consistent with higher relative rates of purifying selection. Further, an analysis of variation in microsatellite allele lengths (lnRV), and heterozygosities (lnRH), indicate recent selective sweeps at the eSTRs. Mean microsatellite allele lengths at four eSTRs within populations are significantly correlated with latitude consistent with the predictions of the tuning-knob model which predicts stepwise relationships between microsatellite allele length and phenotypes. This finding suggests that shorter or longer alleles at eSTRs may be favored in climatic extremes. Collectively, our results imply that eSTRs are likely under selection and that they may be playing a role in facilitating local adaptation across a well-defined cline in the common sunflower.

Epistatic selection on a selfish Segregation Distorter supergene - drive, recombination, and genetic load.

Meiotic drive supergenes are complexes of alleles at linked loci that together subvert Mendelian segregation resulting in preferential transmission. In males, the most common mechanism of drive involves the disruption of sperm bearing one of a pair of alternative alleles. While at least two loci are important for male drive-the driver and the target-linked modifiers can enhance drive, creating selection pressure to suppress recombination. In this work, we investigate the evolution and genomic consequences of an autosomal, multilocus, male meiotic drive system, Segregation Distorter (SD) in the fruit fly, Drosophila melanogaster. In African populations, the predominant SD chromosome variant, SD-Mal, is characterized by two overlapping, paracentric inversions on chromosome arm 2R and nearly perfect (~100%) transmission. We study the SD-Mal system in detail, exploring its components, chromosomal structure, and evolutionary history. Our findings reveal a recent chromosome-scale selective sweep mediated by strong epistatic selection for haplotypes carrying Sd, the main driving allele, and one or more factors within the double inversion. While most SD-Mal chromosomes are homozygous lethal, SD-Mal haplotypes can recombine with other, complementing haplotypes via crossing over, and with wildtype chromosomes via gene conversion. SD-Mal chromosomes have nevertheless accumulated lethal mutations, excess non-synonymous mutations, and excess transposable element insertions. Therefore, SD-Mal haplotypes evolve as a small, semi-isolated subpopulation with a history of strong selection. These results may explain the evolutionary turnover of SD haplotypes in different populations around the world and have implications for supergene evolution broadly.

Ancestry Specific variation in neuropsychological disorders among the South Asian population

The enormous genetic diversity in South Asia resulting from a long and complex admixture history resulted in the emergence of variation in various traits and variations in disease susceptibility. Neuropsychological disorders are one such example that shows variation at the population level. In this study, we aimed at understanding the variation in neuropsychological disorders at the population level among South Asian populations by curating, comparing and contrasting single nucleotide polymorphisms (SNPs), known to be associated with the same. Whole-genome data comprising of 1662 South Asians, belonging to 241 distinct populations were obtained from the database of Dr. David Reich, Harvard Medical School, USA. Principal Component Analysis (PCA) revealed that the Ancestral Tibeto Burman (ATB) genomes form a distinct and distinguishable cluster for the SNPs known to be associated with neuropsychological disorders. Identical By Descent (IBD) analysis showed that out of the top seven populations in terms of IBD sharing, six are from Southern India indicating that these populations may have undergone a recent selective sweep for these SNPs. Further, out of the top ten genomes, according to the number of genomes fixed for the minor alleles, seven were from Southern India. Furthermore, several indigenous populations from South India depicted high F values (&gt;0.25) for SNPs associated with neuropsychological disorders, indicating higher susceptibility for neuropsychological disorders among these South Indian populations. Interestingly, we found that most of the SNPs, fixed for the alternative alleles, were also found to be fixed among the ancient genomes from Indus Valley Civilization (IVC), indicating that these SNPs likely got transmitted to various modern-day South Indian populations from IVC.

Phylogeography and evolutionary dynamism of marine gastropods from the Lord Howe Island Group

Context There has been considerable interest in the biogeography of the marine fauna of the Lord Howe Island Group (LHIG) lying ~580 km from continental Australia. Most phylogenetic studies of the group’s marine taxa have included few specimens, limiting their phylogeographic perspective. Aims This investigation was conducted to enhance understanding of the phylogeographic relationships of the LHIG, especially with temperate Australia. Methods The investigation used mitochondrial DNA sequences of the muricid subfamily Ergalataxinae, especially Tenguella, which has a broad distribution in the Indo-West Pacific, and two other gastropod genera, Cellana and Bembicium, each proposed to have endemic species in the LHIG. Key results Most species of LHIG Ergalataxinae have recent genetic connections with tropical landmasses, but Tenguella marginalba is closely related to temperate mainland conspecifics. The endemic taxa Cellana howensis and Bembicium flavescens respectively, have tropical and southern temperate affinities. Conclusions The most recent common ancestors (MCRAs) of LHIG specimens of each of B. flavescens, C. howensis and T. marginalba all date from long after the archipelago’s emergence. Implications The studied species may have arrived only recently on the archipelago. Alternatively, their estimated MCRAs may be due to recent demographic bottlenecks or selective sweeps influenced by Pleistocene environmental changes.

PopHumanVar: an interactive application for the functional characterization and prioritization of adaptive genomic variants in humans.

Adaptive challenges that humans faced as they expanded across the globe left specific molecular footprints that can be decoded in our today's genomes. Different sets of metrics are used to identify genomic regions that have undergone selection. However, there are fewer methods capable of pinpointing the allele ultimately responsible for this selection. Here, we present PopHumanVar, an interactive online application that is designed to facilitate the exploration and thorough analysis of candidate genomic regions by integrating both functional and population genomics data currently available. PopHumanVar generates useful summary reports of prioritized variants that are putatively causal of recent selective sweeps. It compiles data and graphically represents different layers of information, including natural selection statistics, as well as functional annotations and genealogical estimations of variant age, for biallelic single nucleotide variants (SNVs) of the 1000 Genomes Project phase 3. Specifically, PopHumanVar amasses SNV-based information from GEVA, SnpEFF, GWAS Catalog, ClinVar, RegulomeDB and DisGeNET databases, as well as accurate estimations of iHS, nSL and iSAFE statistics. Notably, PopHumanVar can successfully identify known causal variants of frequently reported candidate selection regions, including EDAR in East-Asians, ACKR1 (DARC) in Africans and LCT/MCM6 in Europeans. PopHumanVar is open and freely available at https://pophumanvar.uab.cat.