Detecting selection using extended haplotype homozygosity (EHH)-based statistics in unphased or unpolarized data

Parallel Signatures of Diet Adaptation in the Invasive Common Myna Genome

Invasive species offer uniquely replicated model systems to study rapid adaptation. The common myna (Acridotheres tristis) has been introduced to over a dozen countries and is classified as one of the most invasive birds in the world. Their multiple invasions provide an opportunity to identify repeated adaptation, as invasive populations originated from multiple source populations. We compared whole-genome resequencing data from 80 individuals from four native and seven invasive populations, representing two independent introduction pathways. Results from two different selection scan methods were combined and identified a strongly selected region on chromosome 8 that spans two copies of AMY2A, part of the alpha-amylase gene family, a putative ncRNA and an insertion-deletion structural variant (SV) that contains an ERVK transposable element (TE). Outlier SNPs and the SV are polymorphic in native populations, but fixed or close-to-fixed in the two invasive pathways, with the fixation of the same alleles in two independent lineages providing evidence for parallel selection on standing variation. Intriguingly, the second copy of AMY2A has a non-conservative missense mutation at a phylogenetically conserved site. This mutation, alongside variation in the SV, TE and ncRNA, provide possible routes for changes to protein function or expression. AMY2A has been associated with human commensalism in house sparrows, and genes in this family have been linked to adaptation to high-starch diets in humans and dogs. This study illustrates the value of replicated analyses within and across species to understand rapid adaptation at the molecular level.

Spread of yellow-bill-color alleles favored by selection in the long-tailed finch hybrid system

Carotenoid pigments produce the yellow and red colors of birds and other vertebrates. Despite their importance in social signaling and sexual selection, our understanding of how carotenoid ornamentation evolves in nature remains limited. Here, we examine the long-tailed finch Poephila acuticauda, an Australian songbird with a yellow-billed western subspecies acuticauda and a red-billed eastern subspecies hecki, which hybridize where their ranges overlap. We found that yellow bills can be explained by the loss of C(4)-oxidation, thus preventing yellow dietary carotenoids from being converted to red. Combining linked-read genomic sequencing and reflectance spectrophotometry measurements of bill color collected from wild-sampled finches and laboratory crosses, we identify four loci that together explain 53% of variance in this trait. The two loci of largest effect contain the genes CYP2J19, an essential enzyme for producing red carotenoids, and TTC39B, an enhancer of carotenoid metabolism. A paucity of protein-coding changes and an enrichment of associated upstream variants suggest that the loss of C(4)-oxidation results from cis-regulatory evolution. Evolutionary genealogy reconstruction indicates that the red-billed phenotype is ancestral and that yellow alleles at CYP2J19 and TTC39B first arose and fixed in acuticauda approximately 100 kya. Yellow alleles subsequently introgressed into hecki less than 5 kya. Across all color loci, acuticauda-derived variants show evidence of selective sweeps, implying that yellow bill coloration has been favored by natural selection. Our study illustrates how evolutionary transitions between yellow and red coloration can be achieved by successive selective events acting on regulatory changes at a few interacting genes.

Assessment of genome complementarity in three beef-on-dairy crossbreds reveals sire-specific effects on production traits with comparable rates of genomic inbreeding reduction

BACKGROUND

Crossbreeding beef bulls with dairy cows can improve the economic value and fitness of calves not entering dairy production owing to increased meat yield and heterosis. However, outcrossing might reduce the dosage of alleles that confer local adaptation or result in a higher risk of dystocia due to increased calf size. Given the clear phenotypic differences between beef breeds, the varying phylogenetic distances between beef and dairy breeds, and the genomic variations within breeds, the attainable economic and fitness gains of calves will strongly depend on the selection of sires for crossing. Thus, the aim of this study was to assess genome complementarity between Angus (AAN), Limousin (LIM), or Simmental (SIM) beef bulls and Brown Swiss (BSW) dairy cows by quantifying genomic inbreeding reduction in F1 crosses and identifying genes potentially under BSW-specific selection that might be affected by outcrossing.

RESULTS

Low-pass sequencing data from 181 cows, 34 bulls, and 301 of their F1 progeny, and body weight and carcass composition measurements of 248 F1s were obtained. The high genomic inbreeding levels detected in the BSW cows were substantially reduced in the crossbreds, with only minor differences between the sire breeds. In the BSW cows, 585 candidate genes under selection were identified, overrepresenting genes associated with milk, meat and carcass, and production traits. Only a few genes were strongly differentiated at nonsynonymous variants between the BSW and beef breeds, including four tightly clustered genes (FAM184B, NCAPG, DCAF16, and LCORL) nearly fixed for alternate alleles in the BSW cows but mostly heterozygous or homozygous for the reference alleles in the AAN and LIM bulls. The alternate allele dosage at these genes significantly correlated with reduced carcass weight and protein mass in F1s.

CONCLUSION

Some of the few genes that were highly divergent between the BSW and beef breeds at nonsynonymous variants were likely under strong selection for reduced carcass weight in the BSW breed, potentially due to trade-offs between beef and dairy productions. As alleles with opposing effects still segregate in beef cattle, marker-assisted selection of mating pairs may be used to modulate the desired phenotypes and simultaneously decrease genomic inbreeding.

Selection signatures for high altitude adaptation in livestock: A review

High altitude adapted livestock species (cattle, yak, goat, sheep, and horse) has critical role in the human socioeconomic sphere and acts as good source of animal source products including milk, meat, and leather, among other things. These species sustain production and reproduction even in harsh environments on account of adaptation resulting from continued evolution of beneficial traits. Selection pressure leads to various adaptive strategies in livestock whose footprints are evident at the different genomic sites as the "Selection Signature". Scrutiny of these signatures provides us crucial insight into the evolutionary process and domestication of livestock adapted to diverse climatic conditions. These signatures have the potential to change the sphere of animal breeding and further usher the selection programmes in right direction. Technological revolution and recent strides made in genomic studies has opened the routes for the identification of selection signatures. Numerous statistical approaches and bioinformatics tools have been developed to detect the selection signature. Consequently, studies across years have identified candidate genes under selection region found associated with numerous traits which have a say in adaptation to high-altitude environment. This makes it pertinent to have a better understanding about the selection signature, the ways to identify and how to utilize them for betterment of livestock populations as well as farmers. This review takes a closer look into the general concept, various methodologies, and bioinformatics tools commonly employed in selection signature studies and summarize the results of recent selection signature studies related to high-altitude adaptation in various livestock species. This review will serve as an informative and useful insight for researchers and students in the field of animal breeding and evolutionary biology.

Signatures of selection analyses reveal genomic differences among three heavy pig breeds that constitute the genetic backbone of a dry-cured ham production system

The Italian pig farming industry is unique in its focus on raising heavy pigs primarily for the production of high-quality dry-cured hams. These products require pigs to be slaughtered at a live weight of around 170 kg at 9 months of age. The primary breeds used in this system are Italian Duroc, Italian Landrace, and Italian Large White which are crossed to produce lines that meet standard requirements. Over the past four decades, selection and breeding programmes for these breeds have been subjected to distinct selective pressures to highlight the characteristics of each breed. In this study, we investigated the genome of these breeds by analysing high-density single nucleotide polymorphism data from over 9 000 pigs to scan for signatures of selection using four different methods, two within breeds and two across breeds. This allowed to identify the genomic regions that differentiate these breeds as well as any relevant genes and biological terms. On a global scale, we found that the Italian Duroc breed exhibited a higher genetic differentiation from the Italian Landrace and Italian Large White breeds, with a pairwise FST value of 0.20 compared with the 0.13 between Italian Landrace and Italian Large White. This may reflect either their different origins or the different breeding goals, which are more similar for the Italian Landrace and Italian Large White breeds. Despite these genetic differences at a global level, few signatures of selection regions reached complete fixation, possibly due to challenges in detecting selection linked to quantitative polygenic traits. The differences among the three breeds are confirmed by the low level of overlap in the regions detected. Genetic enrichment analyses of the three breeds revealed pathways and genes related to various productive traits associated with growth and fat deposition. This may indicate a common selection direction aimed at enhancing specific production traits, though different biological mechanisms are likely targeted by the same directional selection in these three breeds. Therefore, these genes may play a critical role in determining the distinctive characteristics of Italian Duroc, Italian Landrace, and Italian Large White, and potentially influence the traits in crossbred pigs derived from them. Overall, the insights gained from this study will contribute to understanding how directional selection has shaped the genome of these heavy pig breeds and to better address selection strategies aimed at enhancing the meat processing industry linked with dry-cured ham production chains.

Search for signals of positive selection of circadian rhythm genes PER1, PER2, PER3 in different human populations

The diversity of geographically distributed human populations shows considerable variation in external and internal traits of individuals. Such differences are largely attributed to genetic adaptation to various environmental influences, which include changes in climatic conditions, variations in sleep and wakefulness, dietary variations, and others. Whole-genome data from individuals of different populations make it possible to determine the specific genetic sites responsible for adaptations and to further understand the genetic structure underlying human adaptive characteristics. In this article, we searched for signals of single nucleotide polymorphisms (SNPs) under selection pressure in people of different populations. To identify selection signals in different population groups, the PER1, PER2 and PER3 genes that are involved in the coordination of thermogenic functions and regulation of circadian rhythms, which is directly reflected in the adaptive abilities of the organism, were investigated. Data were analyzed using publicly available data from the 1000 Genomes Project for 23 populations. The Extended Haplotype Homozygosity Score statistical method was chosen to search for traces of selection. The comparative analysis performed identified points subject to selection pressure. The SNPs were annotated through the GWAS catalog and manually by analyzing Internet resources. This study suggests that living conditions, climate, and other external factors directly influence the genetic structure of populations and vary across races and geographic locations. In addition, many of the selection variants in the PER1, PER2, PER3 genes appear to regulate biological processes that are associated with major modern diseases, including obesity, cancer, metabolic syndrome, bipolar personality disorder, depression, rheumatoid arthritis, diabetes mellitus, lupus erythematosus, stroke and Alzheimer's disease, making them extremely interesting targets for further research aimed at identifying the genetic causes of human disease.

Genomic Reconstruction of the Successful Establishment of a Feralized Bovine Population on the Subantarctic Island of Amsterdam

The feral cattle of the subantarctic island of Amsterdam provide an outstanding case study of a large mammalian population that was established by a handful of founders and thrived within a few generations in a seemingly inhospitable environment. Here, we investigated the genetic history and composition of this population using genotyping and sequencing data. Our inference showed an intense but brief founding bottleneck around the late 19th century and revealed contributions from European taurine and Indian Ocean Zebu in the founder ancestry. Comparative analysis of whole-genome sequences further revealed a moderate reduction in genetic diversity despite high levels of inbreeding. The brief and intense bottleneck was associated with high levels of drift, a flattening of the site frequency spectrum and a slight relaxation of purifying selection on mildly deleterious variants. Unlike some populations that have experienced prolonged reductions in effective population size, we did not observe any significant purging of highly deleterious variants. Interestingly, the population's success in the harsh environment can be attributed to preadaptation from their European taurine ancestry, suggesting no strong bioclimatic challenge, and also contradicting evidence for insular dwarfism. Genome scan for footprints of selection uncovered a majority of candidate genes related to nervous system function, likely reflecting rapid feralization driven by behavioral changes and complex social restructuring. The Amsterdam Island cattle offers valuable insights into rapid population establishment, feralization, and genetic adaptation in challenging environments. It also sheds light on the unique genetic legacies of feral populations, raising ethical questions according to conservation efforts.

Signatures of positive selection after the introduction of genomic selection in the Finnish Ayrshire population

The Finnish Ayrshire (FAY) belongs to the Nordic Red breeds and is characterized by high milk yield, high milk components, good fertility, and functional conformation. The FAY breeding program is based on genomic selection. Despite the benefits of selection on breeding values, autozygosity in the genome may increase due to selection, and increased autozygosity may cause inbreeding depression in selected traits. However, there is lack of studies concerning selection signatures in the FAY after genomic selection introduction. The aim of this study was to identify signatures of selection in FAY after the introduction of genomic selection. Genomic data included 45,834 SNPs. The genotyped animals were divided into 2 groups: animals born before genomic selection introduction (6,108 cows) and animals born after genomic selection introduction (47,361 cows). We identified the selection signatures using 3 complementary methods: 2 based on identification of selection signatures from runs of homozygosity (ROH) islands and one based on the decay of site-specific extended haplotype between populations at SNP sites (Rsb). In total, we identified 34 ROH islands on chromosomes 1, 3, 6, 8, 12-15, 17, 19, 22, and 26 in FAY animals born before genomic selection (between 1980 and 2011) and 30 ROH islands on chromosomes 1-3, 13-17, 22, and 25-26 in FAY animals born after genomic selection introduction (between 2015 and 2020). We additionally detected 22 ΔROH islands on chromosomes 2-3, 11, 13, 14, 16, 18, 20, and 25-26. Finally, a total of 31 Rsb regions on chromosomes 2, 3, 14, 18, 20, and 25 were identified. Based on the results, genomic selection has favored certain alleles and haplotypes on genomic regions related to traits relevant in the FAY breeding program: milk production, fertility, growth, beef production traits, and feed efficiency. Several genes related to these traits (e.g., PLA2G4A, MECR, CHUK, COX15, RICTOR, SHISA9, and SEMA4G) overlapped or partially overlapped the observed selection signature regions. The association of genotypes within these regions and their effects on traits relevant in the FAY breeding program should be studied and genetic regions undergoing selection monitored in the FAY population.

Genomic Signatures of Microgeographic Adaptation in Anopheles coluzzii Along an Anthropogenic Gradient in Gabon

Species distributed across heterogeneous environments often evolve locally adapted populations, but understanding how these persist in the presence of homogenizing gene flow remains puzzling. In Gabon, Anopheles coluzzii, a major African malaria mosquito is found along an ecological gradient, including a sylvatic population, away of any human presence. This study identifies into the genomic signatures of local adaptation in populations from distinct environments including the urban area of Libreville, and two proximate sites 10km apart in the La Lopé National Park (LLP), a village and its sylvatic neighborhood. Whole genome re-sequencing of 96 mosquitoes unveiled ∼ 5.7millions high-quality single nucleotide polymorphisms. Coalescent-based demographic analyses suggest an ∼ 8,000-year-old divergence between Libreville and La Lopé populations, followed by a secondary contact ( ∼ 4,000 ybp) resulting in asymmetric effective gene flow. The urban population displayed reduced effective size, evidence of inbreeding, and strong selection pressures for adaptation to urban settings, as suggested by the hard selective sweeps associated with genes involved in detoxification and insecticide resistance. In contrast, the two geographically proximate LLP populations showed larger effective sizes, and distinctive genomic differences in selective signals, notably soft-selective sweeps on the standing genetic variation. Although neutral loci and chromosomal inversions failed to discriminate between LLP populations, our findings support that microgeographic adaptation can swiftly emerge through selection on standing genetic variation despite high gene flow. This study contributes to the growing understanding of evolution of populations in heterogeneous environments amid ongoing gene flow and how major malaria mosquitoes adapt to human.

Significance

Anopheles coluzzii , a major African malaria vector, thrives from humid rainforests to dry savannahs and coastal areas. This ecological success is linked to its close association with domestic settings, with human playing significant roles in driving the recent urban evolution of this mosquito. Our research explores the assumption that these mosquitoes are strictly dependent on human habitats, by conducting whole-genome sequencing on An. coluzzii specimens from urban, rural, and sylvatic sites in Gabon. We found that urban mosquitoes show de novo genetic signatures of human-driven vector control, while rural and sylvatic mosquitoes exhibit distinctive genetic evidence of local adaptations derived from standing genetic variation. Understanding adaptation mechanisms of this mosquito is therefore crucial to predict evolution of vector control strategies.

Genomic Analyses Capture the Human-Induced Demographic Collapse and Recovery in a Wide-Ranging Cervid

The glacial cycles of the Quaternary heavily impacted species through successions of population contractions and expansions. Similarly, populations have been intensely shaped by human pressures such as unregulated hunting and land use changes. White-tailed and mule deer survived in different refugia through the Last Glacial Maximum, and their populations were severely reduced after the European colonization. Here, we analyzed 73 resequenced deer genomes from across their North American range to understand the consequences of climatic and anthropogenic pressures on deer demographic and adaptive history. We found strong signals of climate-induced vicariance and demographic decline; notably, multiple sequentially Markovian coalescent recovers a severe decline in mainland white-tailed deer effective population size (Ne) at the end of the Last Glacial Maximum. We found robust evidence for colonial overharvest in the form of a recent and dramatic drop in Ne in all analyzed populations. Historical census size and restocking data show a clear parallel to historical Ne estimates, and temporal Ne/Nc ratio shows patterns of conservation concern for mule deer. Signatures of selection highlight genes related to temperature, including a cold receptor previously highlighted in woolly mammoth. We also detected immune genes that we surmise reflect the changing land use patterns in North America. Our study provides a detailed picture of anthropogenic and climatic-induced decline in deer diversity and clues to understanding the conservation concerns of mule deer and the successful demographic recovery of white-tailed deer.

Comparative genomic analyses provide new insights into evolutionary history and conservation genomics of gorillas

Genome sequencing is a powerful tool to understand species evolutionary history, uncover genes under selection, which could be informative of local adaptation, and infer measures of genetic diversity, inbreeding and mutational load that could be used to inform conservation efforts. Gorillas, critically endangered primates, have received considerable attention and with the recently sequenced Bwindi mountain gorilla population, genomic data is now available from all gorilla subspecies and both mountain gorilla populations. Here, we reanalysed this rich dataset with a focus on evolutionary history, local adaptation and genomic parameters relevant for conservation. We estimate a recent split between western and eastern gorillas of 150,000-180,000 years ago, with gene flow around 20,000 years ago, primarily between the Cross River and Grauer's gorilla subspecies. This gene flow event likely obscures evolutionary relationships within eastern gorillas: after excluding putatively introgressed genomic regions, we uncover a sister relationship between Virunga mountain gorillas and Grauer's gorillas to the exclusion of Bwindi mountain gorillas. This makes mountain gorillas paraphyletic. Eastern gorillas are less genetically diverse and more inbred than western gorillas, yet we detected lower genetic load in the eastern species. Analyses of indels fit remarkably well with differences in genetic diversity across gorilla taxa as recovered with nucleotide diversity measures. We also identified genes under selection and unique gene variants specific for each gorilla subspecies, encoding, among others, traits involved in immunity, diet, muscular development, hair morphology and behavior. The presence of this functional variation suggests that the subspecies may be locally adapted. In conclusion, using extensive genomic resources we provide a comprehensive overview of gorilla genomic diversity, including a so-far understudied Bwindi mountain gorilla population, identify putative genes involved in local adaptation, and detect population-specific gene flow across gorilla species.

selscan 2.0: scanning for sweeps in unphased data

SUMMARY

Several popular haplotype-based statistics for identifying recent or ongoing positive selection in genomes require knowledge of haplotype phase. Here, we provide an update to selscan which implements a re-definition of these statistics for use in unphased data.

AVAILABILITY AND IMPLEMENTATION

Source code and binaries are freely available at https://github.com/szpiech/selscan, implemented in C/C++, and supported on Linux, Windows, and MacOS.

Genomic selection pressure discovery using site-frequency spectrum and reduced local variability statistics in Pakistani Dera-Din-Panah goat

BACKGROUND

Population geneticists have long sought to comprehend various selection traces accumulated in the goat genome due to natural or human driven artificial selection through breeding practices, which led the wild animals to domestication, so understanding evolutionary process may helpful to utilize the full genetic potential of goat genome.

METHODS AND RESULTS

As a step forward to pinpoint the selection signals in Pakistani Dera-Din-Panah (DDP) goat, whole-genome pooled sequencing (n = 12) was performed, and 618,236,192 clean paired-end reads were mapped against ARS1 reference goat assembly. Five different selection signature statistics were applied using four site-frequency spectrum (SFS) methods (Tajima's D ([Formula: see text]), Fay and Wu's H ([Formula: see text]), Zeng's E ([Formula: see text]), [Formula: see text]) and one reduced local variability approach named pooled heterozygosity ([Formula: see text]). The under-selection regions were annotated with significant threshold values of [Formula: see text]≥4.7, [Formula: see text]≥6, [Formula: see text]≥2.5, Pool-HMM ≥ 12, and [Formula: see text]≥5 that resulted in accumulative 364 candidate gene hits. The highest genomic selection signals were observed on Chr. 4, 6, 10, 12, 15, 16, 18, 20, and 27 and harbor ADAMTS6, CWC27, RELN, MYCBP2, FGF14, STIM1, CFAP74, GNB1, CALML6, TMEM52, FAM149A, NADK, MMP23B, OPN3, FH, MFHAS1, KLKB1, RRM1, KMO, SPEF2, F11, KIT, KMO, ERI1, ATP8B4, and RHOG genes. Next, the validation of our captured genomic hits was also performed by more than one applied statistics which harbor meat production, immunity, and reproduction associated genes to strengthen our hypothesis of under-selection traits in this Pakistani goat breed. Furthermore, common candidate genes captured by more than one statistical method were subjected to gene ontology and KEGG pathway analysis to get insights of particular biological processes associated with this goat breed.

CONCLUSION

Current perception of genomic architecture of DDP goat provides a better understanding to improve its genetic potential and other economically important traits of medium to large body size, milk, and fiber production by updating the genomic insight driven breeding strategies to boost the livestock and agriculture-based economy of the country.

Exploring genome-wide differentiation and signatures of selection in Italian and North American Holstein populations

Among Italian dairy cattle, the Holstein is the most reared breed for the production of Parmigiano Reggiano protected designation of origin cheese, which represents one of the most renowned products in the entire Italian dairy industry. In this work, we used a medium-density genome-wide data set consisting of 79,464 imputed SNPs to study the genetic structure of Italian Holstein breed, including the population reared in the area of Parmigiano Reggiano cheese production, and assessing its distinctiveness from the North American population. Multidimensional scaling and ADMIXTURE approaches were used to explore the genetic structure among populations. We also investigated putative genomic regions under selection among these 3 populations by combining 4 different statistical methods based either on allele frequencies (single marker and window-based) or extended haplotype homozygosity (EHH; standardized log-ratio of integrated EHH and cross-population EHH). The genetic structure results allowed us to clearly distinguish the 3 Holstein populations; however, the most remarkable difference was observed between Italian and North American stock. Selection signature analyses identified several significant SNPs falling within or closer to genes with known roles in several traits such as milk quality, resistance to disease, and fertility. In particular, a total of 22 genes related to milk production have been identified using the 2 allele frequency approaches. Among these, a convergent signal has been found in the VPS8 gene which resulted to be involved in milk traits, whereas other genes (CYP7B1, KSR2, C4A, LIPE, DCDC1, GPR20, and ST3GAL1) resulted to be associated with quantitative trait loci related to milk yield and composition in terms of fat and protein percentage. In contrast, a total of 7 genomic regions were identified combining the results of standardized log-ratio of integrated EHH and cross-population EHH. In these regions candidate genes for milk traits were also identified. Moreover, this was also confirmed by the enrichment analyses in which we found that the majority of the significantly enriched quantitative trait loci were linked to milk traits, whereas the gene ontology and pathway enrichment analysis pointed to molecular functions and biological processes involved in AA transmembrane transport and methane metabolism pathway. This study provides information on the genetic structure of the examined populations, showing that they are distinguishable from each other. Furthermore, the selection signature analyses can be considered as a starting point for future studies in the identification of causal mutations and consequent implementation of more practical application.

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

The genomic signature of wild-to-crop introgression during the domestication of scarlet runner bean (Phaseolus coccineus L.)

The scarlet runner bean (Phaseolus coccineus) is one of the five domesticated Phaseolus species. It is cultivated in small-scale agriculture in the highlands of Mesoamerica for its dry seeds and immature pods, and unlike the other domesticated beans, P. coccineus is an open-pollinated legume. Contrasting with its close relative, the common bean, few studies focusing on its domestication history have been conducted. Demographic bottlenecks associated with domestication might reduce genetic diversity and facilitate the accumulation of deleterious mutations. Conversely, introgression from wild relatives could be a source of variation. Using Genotyping by Sequencing data (79,286 single-nucleotide variants) from 237 cultivated and wild samples, we evaluated the demographic history of traditional varieties from different regions of Mexico and looked for evidence of introgression between sympatric wild and cultivated populations. Traditional varieties have high levels of diversity, even though there is evidence of a severe initial genetic bottleneck followed by a population expansion. Introgression from wild to domesticated populations was detected, which might contribute to the recovery of the genetic variation. Introgression has occurred at different times: constantly in the center of Mexico; recently in the North West; and anciently in the South. Several factors are acting together to increase and maintain genetic diversity in P. coccineus cultivars, such as demographic expansion and introgression. Wild relatives represent a valuable genetic resource and have played a key role in scarlet runner bean evolution via introgression into traditional varieties.

Inferring Signatures of Positive Selection in Whole-Genome Sequencing Data: An Overview of Haplotype-Based Methods

Signatures of positive selection in the genome are a characteristic mark of adaptation that can reveal an ongoing, recent, or ancient response to environmental change throughout the evolution of a population. New sources of food, climate conditions, and exposure to pathogens are only some of the possible sources of selective pressure, and the rise of advantageous genetic variants is a crucial determinant of survival and reproduction. In this context, the ability to detect these signatures of selection may pinpoint genetic variants that are responsible for a significant change in gene regulation, gene expression, or protein synthesis, structure, and function. This review focuses on statistical methods that take advantage of linkage disequilibrium and haplotype determination to reveal signatures of positive selection in whole-genome sequencing data, showing that they emerge from different descriptions of the same underlying event. Moreover, considerations are provided around the application of these statistics to different species, their suitability for ancient DNA, and the usefulness of discovering variants under selection for biomedicine and public health in an evolutionary medicine framework.