Genetic hitchhiking versus background selection: the controversy and its implications

Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole-genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Enrichment of hard sweeps on the X chromosome compared to autosomes in six Drosophila species

The X chromosome, being hemizygous in males, is exposed one-third of the time increasing the visibility of new mutations to natural selection, potentially leading to different evolutionary dynamics than autosomes. Recently, we found an enrichment of hard selective sweeps over soft selective sweeps on the X chromosome relative to the autosomes in a North American population of Drosophila melanogaster. To understand whether this enrichment is a universal feature of evolution on the X chromosome, we analyze diversity patterns across 6 commonly studied Drosophila species. We find an increased proportion of regions with steep reductions in diversity and elevated homozygosity on the X chromosome compared to autosomes. To assess if these signatures are consistent with positive selection, we simulate a wide variety of evolutionary scenarios spanning variations in demography, mutation rate, recombination rate, background selection, hard sweeps, and soft sweeps and find that the diversity patterns observed on the X are most consistent with hard sweeps. Our findings highlight the importance of sex chromosomes in driving evolutionary processes and suggest that hard sweeps have played a significant role in shaping diversity patterns on the X chromosome across multiple Drosophila species.

Whole-Genome Sequencing for Identifying Candidate Genes Related to the Special Phenotypes of the Taihu Dianzi Pigeon

The Taihu Dianzi pigeon is a breed native to China, and its special piebalding, crest, and polydactyly phenotypes are the result of artificial and natural selection. Here, we analyzed the genetic differences among three kinds of pigeons with different phenotypes at the genomic level. A selective sweep was conducted based on the fixation index (FST) and nucleotide diversity (π) ratio, and the results revealed that MC1R was related to the formation of the distinctive piebalding of the Taihu Dianzi pigeon. Combined with the results of genome-wide association studies, we identified candidate genes associated with the crest (SMYD and STOX2) and polydactyly (SLC52A3 and ANGPT4). The candidate genes identified in this study and their variants may be useful for understanding the genetic mechanism underlying the special phenotypes of the Taihu Dianzi pigeon. This study provides new insights into the genetic factors that may influence the formation of the special piebalding, crest, and polydactyly characteristics in pigeons.

Genomic Architecture Predicts Tree Topology, Population Structuring, and Demographic History in Amazonian Birds

Geographic barriers are frequently invoked to explain genetic structuring across the landscape. However, inferences on the spatial and temporal origins of population variation have been largely limited to evolutionary neutral models, ignoring the potential role of natural selection and intrinsic genomic processes known as genomic architecture in producing heterogeneity in differentiation across the genome. To test how variation in genomic characteristics (e.g. recombination rate) impacts our ability to reconstruct general patterns of differentiation between species that cooccur across geographic barriers, we sequenced the whole genomes of multiple bird populations that are distributed across rivers in southeastern Amazonia. We found that phylogenetic relationships within species and demographic parameters varied across the genome in predictable ways. Genetic diversity was positively associated with recombination rate and negatively associated with species tree support. Gene flow was less pervasive in genomic regions of low recombination, making these windows more likely to retain patterns of population structuring that matched the species tree. We further found that approximately a third of the genome showed evidence of selective sweeps and linked selection, skewing genome-wide estimates of effective population sizes and gene flow between populations toward lower values. In sum, we showed that the effects of intrinsic genomic characteristics and selection can be disentangled from neutral processes to elucidate spatial patterns of population differentiation.

Enrichment of hard sweeps on the X chromosome compared to autosomes in six Drosophila species

The X chromosome, being hemizygous in males, is exposed one third of the time increasing the visibility of new mutations to natural selection, potentially leading to different evolutionary dynamics than autosomes. Recently, we found an enrichment of hard selective sweeps over soft selective sweeps on the X chromosome relative to the autosomes in a North American population of Drosophila melanogaster. To understand whether this enrichment is a universal feature of evolution on the X chromosome, we analyze diversity patterns across six commonly studied Drosophila species. We find an increased proportion of regions with steep reductions in diversity and elevated homozygosity on the X chromosome compared to autosomes. To assess if these signatures are consistent with positive selection, we simulate a wide variety of evolutionary scenarios spanning variations in demography, mutation rate, recombination rate, background selection, hard sweeps, and soft sweeps, and find that the diversity patterns observed on the X are most consistent with hard sweeps. Our findings highlight the importance of sex chromosomes in driving evolutionary processes and suggest that hard sweeps have played a significant role in shaping diversity patterns on the X chromosome across multiple Drosophila species.

Genomic signatures of selection, local adaptation and production type characterisation of East Adriatic sheep breeds

BACKGROUND

The importance of sheep breeding in the Mediterranean part of the eastern Adriatic has a long tradition since its arrival during the Neolithic migrations. Sheep production system is extensive and generally carried out in traditional systems without intensive systematic breeding programmes for high uniform trait production (carcass, wool and milk yield). Therefore, eight indigenous Croatian sheep breeds from eastern Adriatic treated here as metapopulation (EAS), are generally considered as multipurpose breeds (milk, meat and wool), not specialised for a particular type of production, but known for their robustness and resistance to certain environmental conditions. Our objective was to identify genomic regions and genes that exhibit patterns of positive selection signatures, decipher their biological and productive functionality, and provide a "genomic" characterization of EAS adaptation and determine its production type.

RESULTS

We identified positive selection signatures in EAS using several methods based on reduced local variation, linkage disequilibrium and site frequency spectrum (eROHi, iHS, nSL and CLR). Our analyses identified numerous genomic regions and genes (e.g., desmosomal cadherin and desmoglein gene families) associated with environmental adaptation and economically important traits. Most candidate genes were related to meat/production and health/immune response traits, while some of the candidate genes discovered were important for domestication and evolutionary processes (e.g., HOXa gene family and FSIP2). These results were also confirmed by GO and QTL enrichment analysis.

CONCLUSIONS

Our results contribute to a better understanding of the unique adaptive genetic architecture of EAS and define its productive type, ultimately providing a new opportunity for future breeding programmes. At the same time, the numerous genes identified will improve our understanding of ruminant (sheep) robustness and resistance in the harsh and specific Mediterranean environment.

Analysis of the Selection Signal of the Tibetan Black Chicken Genome Based on Whole-Genome Sequencing

BACKGROUND

The Tibetan chicken has adapted well to high altitudes genetically after its long-term habitation in the plateau. In this study, we analyzed the selection signal of Tibetan black chickens (TBCs) and discovered genes associated with the characteristics of germplasm.

METHODS

Whole-genome sequencing (WGS) was used to identify the single-nucleotide polymorphism (SNP) markers and genetic structures in the genome of Tibetan black chickens. Further, we performed a comparative population genomics analysis between the genomic data obtained in this present study and the genomic data for five wild red jungle fowls (RJFs) accessed from the NCBI database (GenBank accession number PRJNA241474). Thereafter, the Fst and Pi selections were used to identify genes under positive selection in the Tibetan black chicken genome.

RESULTS

A total of 9,490,690 SNPs were identified in the Tibetan black chickens. In addition, the results from the gene ontology (GO) analysis showed that 732 genes of TBCs were enriched in a total of 210 GO terms with specific molecular functions such as regulation of cellular catabolic process, the MAPK signaling pathway, regulation of ion transport, growth, morphogenesis and lung alveolus development which may provide a better mechanism to facilitate oxygen transport and utilization in TBCs. Moreover, the results from the KEGG analysis showed that 732 genes of the TBCs were significantly enriched in the calcium signaling pathway, circadian entrainment (ADCY1, GNG7 and PER3), oxytocin signaling pathway and pathways of multiple neurodegeneration diseases. In addition, the CD86 antigen (CD86) was identified as a gene associated with the immune response in chickens. It was also revealed that genes such as TRIT1, HPCAL4, NT5C1A and HEYL were discovered under selection in Tibetan black chickens on chromosome 23. These genes may be related to the local adaptive characteristics of Tibetan black chickens, for instance, NT5C1A and HEYL may be involved in the high-altitude adaption of oxygen delivery in Tibetan black chickens.

CONCLUSIONS

In summary, we found that selection mainly affects the disease resistance and cold acclimatization of Tibetan black chickens. Hence, these results may provide important genetic information for the evolution and breeding of Tibetan black chickens.

Drivers of genomic landscapes of differentiation across a Populus divergence gradient

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation.

Timesweeper: accurately identifying selective sweeps using population genomic time series

Despite decades of research, identifying selective sweeps, the genomic footprints of positive selection, remains a core problem in population genetics. Of the myriad methods that have been developed to tackle this task, few are designed to leverage the potential of genomic time-series data. This is because in most population genetic studies of natural populations, only a single period of time can be sampled. Recent advancements in sequencing technology, including improvements in extracting and sequencing ancient DNA, have made repeated samplings of a population possible, allowing for more direct analysis of recent evolutionary dynamics. Serial sampling of organisms with shorter generation times has also become more feasible due to improvements in the cost and throughput of sequencing. With these advances in mind, here we present Timesweeper, a fast and accurate convolutional neural network-based tool for identifying selective sweeps in data consisting of multiple genomic samplings of a population over time. Timesweeper analyzes population genomic time-series data by first simulating training data under a demographic model appropriate for the data of interest, training a one-dimensional convolutional neural network on said simulations, and inferring which polymorphisms in this serialized data set were the direct target of a completed or ongoing selective sweep. We show that Timesweeper is accurate under multiple simulated demographic and sampling scenarios, identifies selected variants with high resolution, and estimates selection coefficients more accurately than existing methods. In sum, we show that more accurate inferences about natural selection are possible when genomic time-series data are available; such data will continue to proliferate in coming years due to both the sequencing of ancient samples and repeated samplings of extant populations with faster generation times, as well as experimentally evolved populations where time-series data are often generated. Methodological advances such as Timesweeper thus have the potential to help resolve the controversy over the role of positive selection in the genome. We provide Timesweeper as a Python package for use by the community.

Evaluation of the evolutionary genetics and population structure of Culex pipiens pallens in Shandong province, China based on knockdown resistance (kdr) mutations and the mtDNA-COI gene

BACKGROUND

Mosquitoes are important vectors for a range of diseases, contributing to high rates of morbidity and mortality in the human population. Culex pipiens pallens is dominant species of Culex mosquito in northern China and a major vector for both West Nile virus and Bancroftian filariasis. Insecticide application were largely applied to control the mosquito-mediated spread of these diseases, contributing to increasing rates of resistance in the mosquito population. The voltage-gated sodium channel (Vgsc) gene is the target site of pyrethroids, and mutations in this gene cause knockdown resistance (kdr). While these kdr mutations are known to be critical to pyrethroid resistance, their evolutionary origins remain poorly understood. Clarifying the origins of these mutations is potential to guide further vector control and disease prevention efforts. Accordingly, the present study was designed to study the evolutionary genetics of kdr mutations and their association with the population structure of Cx. p. pallens in Shandong province, China.

METHODS

Adult Culex females were collected from Shandong province and subjected to morphological identification under a dissection microscope. Genomic DNA were extracted from the collected mosquitoes, the Vgsc gene were amplified via PCR and sequenced to assess kdr allele frequencies, intron polymorphisms, and kdr codon evolution. In addition, population genetic diversity and related population characteristics were assessed by amplifying and sequencing the mitochondrial cytochrome C oxidase I (COI) gene.

RESULTS

Totally, 263 Cx. p. pallens specimens were used for DNA barcoding and sequencing analyses to assess kdr allele frequencies in nine Culex populations. The kdr codon L1014 in the Vgsc gene identified two non-synonymous mutations (L1014F and L1014S) in the analyzed population. These mutations were present in the eastern hilly area and west plain region of Shandong Province. However, only L1014F mutation was detected in the southern mountainous area and Dongying city of Shandong Province, where the mutation frequency was low. Compared to other cities, population in Qingdao revealed significant genetic differentiation. Spatial kdr mutation patterns are likely attributable to some combination of prolonged insecticide-mediated selection coupled with the genetic isolation of these mosquito populations.

CONCLUSIONS

These data suggest that multiple kdr alleles associated with insecticide resistance are present within the Cx. p. pallens populations of Shandong Province, China. The geographical distributions of kdr mutations in this province are likely that the result of prolonged and extensive insecticide application in agricultural contexts together with frequent mosquito population migrations. In contrast, the low-frequency kdr mutation detected in central Shandong Province populations may originate from the limited selection pressure in this area and the relative genetic isolation. Overall, the study compares the genetic patterns revealed by a functional gene with a neutral marker and demonstrates the combined impact of demographic and selection factors on population structure.

The genetic architecture of behavioral canalization

Behaviors are components of fitness and contribute to adaptive evolution. Behaviors represent the interactions of an organism with its environment, yet innate behaviors display robustness in the face of environmental change, which we refer to as 'behavioral canalization'. We hypothesize that positive selection of hub genes of genetic networks stabilizes the genetic architecture for innate behaviors by reducing variation in the expression of interconnected network genes. Robustness of these stabilized networks would be protected from deleterious mutations by purifying selection or suppressing epistasis. We propose that, together with newly emerging favorable mutations, epistatically suppressed mutations can generate a reservoir of cryptic genetic variation that could give rise to decanalization when genetic backgrounds or environmental conditions change to allow behavioral adaptation.

Interrogating the Roles of Mutation-Selection Balance, Heterozygote Advantage, and Linked Selection in Maintaining Recessive Lethal Variation in Natural Populations

For nearly a century, evolutionary biologists have observed chromosomes that cause lethality when made homozygous persisting at surprisingly high frequencies (>25%) in natural populations of many species. The evolutionary forces responsible for the maintenance of such detrimental mutations have been heavily debated-are some lethal mutations under balancing selection? We suggest that mutation-selection balance alone cannot explain lethal variation in nature and the possibility that other forces play a role. We review the potential that linked selection in particular may drive maintenance of lethal alleles through associative overdominance or linkage to beneficial mutations or by reducing effective population size. Over the past five decades, investigation into this mystery has tapered. During this time, key scientific advances have provided the ability to collect more accurate data and analyze them in new ways, making the underlying genetic bases and evolutionary forces of lethal alleles timely for study once more.

Why does the X chromosome lag behind autosomes in GWAS findings?

The X-chromosome is among the largest human chromosomes. It differs from autosomes by a number of important features including hemizygosity in males, an almost complete inactivation of one copy in females, and unique patterns of recombination. We used data from the Catalog of Published Genome Wide Association Studies to compare densities of the GWAS-detected SNPs on the X-chromosome and autosomes. The density of GWAS-detected SNPs on the X-chromosome is 6-fold lower compared to the density of the GWAS-detected SNPs on autosomes. Differences between the X-chromosome and autosomes cannot be explained by differences in the overall SNP density, lower X-chromosome coverage by genotyping platforms or low call rate of X-chromosomal SNPs. Similar differences in the density of GWAS-detected SNPs were found in female-only GWASs (e.g. ovarian cancer GWASs). We hypothesized that the lower density of GWAS-detected SNPs on the X-chromosome compared to autosomes is not a result of a methodological bias, e.g. differences in coverage or call rates, but has a real underlying biological reason-a lower density of functional SNPs on the X-chromosome versus autosomes. This hypothesis is supported by the observation that (i) the overall SNP density of X-chromosome is lower compared to the SNP density on autosomes and that (ii) the density of genic SNPs on the X-chromosome is lower compared to autosomes while densities of intergenic SNPs are similar.

Host-plant adaptation as a driver of incipient speciation in the fall armyworm (Spodoptera frugiperda)

Background

Divergent selection on host-plants is one of the main evolutionary forces driving ecological speciation in phytophagous insects. The ecological speciation might be challenging in the presence of gene flow and assortative mating because the direction of divergence is not necessarily the same between ecological selection (through host-plant adaptation) and assortative mating. The fall armyworm (FAW), a major lepidopteran pest species, is composed of two sympatric strains, corn and rice strains, named after two of their preferred host-plants. These two strains have been hypothesized to undergo incipient speciation, based on (i) several lines of evidence encompassing both pre- and post-zygotic reproductive isolation, and (ii) the presence of a substantial level of genetic differentiation. Even though the status of these two strains has been established a long time ago, it is still yet to be found whether these two strains indeed exhibit a marked level of genetic differentiation from a large number of genomic loci. Here, we analyzed whole genome sequences from 56 FAW individuals either collected from pasture grasses (a part of the favored host range of the rice strain) or corn to assess the role of host-plant adaptation in incipient speciation.

Results

Principal component analysis of whole genome data shows that the pattern of divergence in the fall armyworm is predominantly explained by the genetic differentiation associated with host-plants. The level of genetic differentiation between corn and rice strains is particularly marked in the Z chromosome. We identified one autosomal locus and two Z chromosome loci targeted by selective sweeps specific to rice strain and corn strain, respectively. The autosomal locus has both increased DXY and FST while the Z chromosome loci had decreased DXY and increased FST.

Conclusion

These results show that the FAW population structure is dominated by the genetic differentiation between corn and rice strains. This differentiation involves divergent selection targeting at least three loci, which include a locus potentially causing reproductive isolation. Taken together, these results suggest the evolutionary scenario that host-plant speciation is a driver of incipient speciation in the fall armyworm.

Evidence of selection in the uncoupling protein 1 gene region suggests local adaptation to solar irradiance in savannah monkeys (Chlorocebus spp.)

In the last 300 thousand years, the genus Chlorocebus expanded from equatorial Africa into the southernmost latitudes of the continent, where colder climate was a probable driver of natural selection. We investigated population-level genetic variation in the mitochondrial uncoupling protein 1 (UCP1) gene region-implicated in non-shivering thermogenesis (NST)-in 73 wild savannah monkeys from three taxa representing this southern expansion (Chlorocebus pygerythrus hilgerti, Chlorocebus cynosuros and Chlorocebus pygerythrus pygerythrus) ranging from Kenya to South Africa. We found 17 single nucleotide polymorphisms with extended haplotype homozygosity consistent with positive selective sweeps, 10 of which show no significant linkage disequilibrium with each other. Phylogenetic generalized least-squares modelling with ecological covariates suggest that most derived allele frequencies are significantly associated with solar irradiance and winter precipitation, rather than overall low temperatures. This selection and association with irradiance is demonstrated by a relatively isolated population in the southern coastal belt of South Africa. We suggest that sunbathing behaviours common to savannah monkeys, in combination with the strength of solar irradiance, may mediate adaptations to thermal stress via NST among savannah monkeys. The variants we discovered all lie in non-coding regions, some with previously documented regulatory functions, calling for further validation and research.

Ultra-rare and common genetic variant analysis converge to implicate negative selection and neuronal processes in the aetiology of schizophrenia

Both common and rare genetic variants (minor allele frequency >1% and <0.1% respectively) have been implicated in the aetiology of schizophrenia. In this study, we integrate single-cell gene expression data with publicly available Genome-Wide Association Study (GWAS) and exome sequenced data in order to investigate in parallel, the enrichment of common and (ultra-)rare variants related to schizophrenia in several functionally relevant gene-sets. Four types of gene-sets were constructed 1) protein-truncating variant (PTV)-intolerant (PI) genes 2) genes expressed in brain cell types and neurons ascertained from mouse and human brain tissue 3) genes defined by synaptic function and location and 4) intersection genes, i.e., PI genes that are expressed in the human and mouse brain cell gene-sets. We show that common as well as ultra-rare schizophrenia-associated variants are overrepresented in PI genes, in excitatory neurons from the prefrontal cortex and hippocampus, medium spiny neurons, and genes enriched for synaptic processes. We also observed stronger enrichment in the intersection genes. Our findings suggest that across the allele frequency spectrum, genes and genetic variants likely to be under stringent selection, and those expressed in particular brain cell types, are involved in the same biological pathways influencing the risk for schizophrenia.

Past and recent farming degrades aquatic insect genetic diversity

Recent declines in once-common species are triggering concern that an environmental crisis point has been reached. Yet, the lack of long abundance time series data for most species can make it difficult to attribute these changes to anthropogenic causes, and to separate them from normal cycles. Genetic diversity, on the other hand, is sensitive to past and recent environmental changes, and reflects a measure of a populations' potential to adapt to future stressors. Here, we consider whether patterns of genetic diversity among aquatic insects can be linked to historical and recent patterns of land use change. We collated mitochondrial cytochrome c oxidase subunit I (COI) variation for >700 aquatic insect species across the United States, where patterns of agricultural expansion and intensification have been documented since the 1800s. We found that genetic diversity was lowest in regions where cropland was historically (pre-1950) most extensive, suggesting a legacy of past environmental harm. Genetic diversity further declined where cropland has since expanded, even after accounting for climate and sampling effects. Notably though, genetic diversity also appeared to rebound where cropland has diminished. Our study suggests that genetic diversity at the community level can be a powerful tool to infer potential population declines and rebounds over longer time spans than is typically possible with ecological data. For the aquatic insects that we considered, patterns of land use many decades ago appear to have left long-lasting damage to genetic diversity that could threaten evolutionary responses to rapid global change.

Machine-Learning Prospects for Detecting Selection Signatures Using Population Genomics Data

Natural selection has been given a lot of attention because it relates to the adaptation of populations to their environments, both biotic and abiotic. An allele is selected when it is favored by natural selection. Consequently, the favored allele increases in frequency in the population and neighboring linked variation diminishes, causing so-called selective sweeps. A high-throughput genomic sequence allows one to disentangle the evolutionary forces at play in populations. With the development of high-throughput genome sequencing technologies, it has become easier to detect these selective sweeps/selection signatures. Various methods can be used to detect selective sweeps, from simple implementations using summary statistics to complex statistical approaches. One of the important problems of these statistical models is the potential to provide inaccurate results when their assumptions are violated. The use of machine learning (ML) in population genetics has been introduced as an alternative method of detecting selection by treating the problem of detecting selection signatures as a classification problem. Since the availability of population genomics data is increasing, researchers may incorporate ML into these statistical models to infer signatures of selection with higher predictive accuracy and better resolution. This article describes how ML can be used to aid in detecting and studying natural selection patterns using population genomic data.

Extreme variation in recombination rate and genetic diversity along the Sylvioidea neo-sex chromosome

Recombination strongly impacts sequence evolution by affecting the extent of linkage and the efficiency of selection. Here, we study recombination over the Z chromosome in great reed warblers (Acrocephalus arundinaceus) using pedigree‐based linkage mapping. This species has extended Z and W chromosomes (“neo‐sex chromosomes”) formed by a fusion between a part of chromosome 4A and the ancestral sex chromosomes, which provides a unique opportunity to assess recombination and sequence evolution in sex‐linked regions of different ages. We assembled an 87.54 Mbp and 90.19 cM large Z with a small pseudoautosomal region (0.89 Mbp) at one end and the fused Chr4A‐part at the other end of the chromosome. A prominent feature in our data was an extreme variation in male recombination rate along Z with high values at both chromosome ends, but an apparent lack of recombination over a substantial central section, covering 78% of the chromosome. The nonrecombining region showed a drastic loss of genetic diversity and accumulation of repeats compared to the recombining parts. Thus, our data emphasize a key role of recombination in affecting local levels of polymorphism. Nonetheless, the evolutionary rate of genes (dN/dS) did not differ between high and low recombining regions, suggesting that the efficiency of selection on protein‐coding sequences can be maintained also at very low levels of recombination. Finally, the Chr4A‐derived part showed a similar recombination rate as the part of the ancestral Z that did recombine, but its sequence characteristics reflected both its previous autosomal, and current Z‐linked, recombination patterns.

Selection signatures in melanocortin-1 receptor gene of turkeys (Meleagris gallopavo) raised in hot humid tropics

Feather colours are used by avian species for defense, adaptation and signaling. Melanocortin-1 receptor (MC1R) gene is one of the genes responsible for feather colour. This study identified selection signatures in MC1R gene of Nigerian indigenous turkeys (NIT) using British United turkeys (BUT) as control breed to investigate the evolutionary processes that have shaped NIT with various feather colours. Complete MC1R gene of 146 NIT (76 males and 70 females) and 32 BUT (18 males and 14 females) were sequenced. Transition/transversion and codon usage biases were predicted using MEGA v6 software. The selective force acting on the gene was predicted using HyPhy software. The FST values were estimated using Arlequin v3.5. The highest transition/transversion bias was predicted for white BUT (1.00) while the lowest was predicted for black NIT (0.50). Negative dN-dS values, indicative of purifying selection, were observed in MC1R gene of all the turkeys. The highest pairwise FST was observed between the MC1R gene of white BUT and black NIT while the least was observed between lavender NIT and white NIT. No recombination event was observed in black NIT and white BUT. The relative synonymous codon usage was the same among different colours for some codons. Presence of purifying selection in MC1R gene of all the turkeys with different feather colours confirms that the gene plays role in many biological processes such as feather colouration, behaviour, pain perception, immunity, growth and adaptation. The results also suggested that the genetic mechanisms generating different feather colours in turkeys are conserved.

Molecular Evolutionary Rate Predicts Intraspecific Genetic Polymorphism and Species-Specific Selection

It is unknown what determines genetic diversity and how genetic diversity is associated with various biological traits. In this work, we provide insight into these issues. By comparing genetic variation of 14,671 mammalian gene trees with thousands of individual human, chimpanzee, gorilla, mouse, and dog/wolf genomes, we found that intraspecific genetic diversity can be predicted by long-term molecular evolutionary rates rather than de novo mutation rates. This relationship was established during the early stage of mammalian evolution. Moreover, we developed a method to detect fluctuations of species-specific selection on genes based on the deviations of intraspecific genetic diversity predicted from long-term rates. We showed that the evolution of epithelial cells, rather than connective tissue, mainly contributed to morphological evolution of different species. For humans, evolution of the immune system and selective sweeps caused by infectious diseases are the most representative examples of adaptive evolution.

Whole-genome resequencing to investigate the determinants of the multi-lumbar vertebrae trait in sheep

Multi-lumbar vertebrae trait is a beneficial mutation that can significantly improve livestock meat production. However, the genetic basis of the multi-lumbar vertebrae in sheep is still unclear. Here, we analysed the number of lumbar vertebrae of Duolang sheep and found three different traits of lumbar vertebrae number. Compared with the normal sheep, the length and weight of animal carcass from the multi-lumbar vertebrae sheep increased by 2.21 cm and 0.78 kg, respectively. We performed high-throughput genome resequencing on multi-lumbar vertebrae (n = 18) and normal (n = 11) Duolang sheep and obtained a total of more than 528.87 GB data. We found that the most significantly selective region were located in the 49.68-49.74 MB of chromosome 4 by selective-sweep analysis. We annotated this region and found that it contains SFRP4 which is known to regulate bone development. We further used the PCR-SSCP technology to detect the single nucleotide polymorphism (SNP) of the putative candidate SFRP4 and found that the two SNPs (rs600370085:C > T and rs415133338: A > G) of this gene were significantly associated with the multi-lumbar vertebrae of Duolang sheep. Our study indicates that the SFRP4 may be a potential major gene that affects the number of lumbar vertebrae in Duolang sheep, and has the potential to be utilized for sheep breeding in the future.

Bottleneck-associated changes in the genomic landscape of genetic diversity in wild lynx populations

Demographic bottlenecks generally reduce genetic diversity through more intense genetic drift, but their net effect may vary along the genome due to the random nature of genetic drift and to local effects of recombination, mutation, and selection. Here, we analyzed the changes in genetic diversity following a bottleneck by comparing whole-genome diversity patterns in populations with and without severe recent documented declines of Iberian (Lynx pardinus, n = 31) and Eurasian lynx (Lynx lynx, n = 29). As expected, overall genomic diversity correlated negatively with bottleneck intensity and/or duration. Correlations of genetic diversity with divergence, chromosome size, gene or functional site content, GC content, or recombination were observed in nonbottlenecked populations, but were weaker in bottlenecked populations. Also, functional features under intense purifying selection and the X chromosome showed an increase in the observed density of variants, even resulting in higher θ W diversity than in nonbottlenecked populations. Increased diversity seems to be related to both a higher mutational input in those regions creating a large collection of low-frequency variants, a few of which increase in frequency during the bottleneck to the point they become detectable with our limited sample, and the reduced efficacy of purifying selection, which affects not only protein structure and function but also the regulation of gene expression. The results of this study alert to the possible reduction of fitness and adaptive potential associated with the genomic erosion in regulatory elements. Further, the detection of a gain of diversity in ultra-conserved elements can be used as a sensitive and easy-to-apply signature of genetic erosion in wild populations.

Population Genomics on the Fly: Recent Advances in Drosophila

Drosophila melanogaster, a small dipteran of African origin, represents one of the best-studied model organisms. Early work in this system has uniquely shed light on the basic principles of genetics and resulted in a versatile collection of genetic tools that allow to uncover mechanistic links between genotype and phenotype. Moreover, given its worldwide distribution in diverse habitats and its moderate genome-size, Drosophila has proven very powerful for population genetics inference and was one of the first eukaryotes whose genome was fully sequenced. In this book chapter, we provide a brief historical overview of research in Drosophila and then focus on recent advances during the genomic era. After describing different types and sources of genomic data, we discuss mechanisms of neutral evolution including the demographic history of Drosophila and the effects of recombination and biased gene conversion. Then, we review recent advances in detecting genome-wide signals of selection, such as soft and hard selective sweeps. We further provide a brief introduction to background selection, selection of noncoding DNA and codon usage and focus on the role of structural variants, such as transposable elements and chromosomal inversions, during the adaptive process. Finally, we discuss how genomic data helps to dissect neutral and adaptive evolutionary mechanisms that shape genetic and phenotypic variation in natural populations along environmental gradients. In summary, this book chapter serves as a starting point to Drosophila population genomics and provides an introduction to the system and an overview to data sources, important population genetic concepts and recent advances in the field.

Divergence and Gene Flow Between Fusarium subglutinans and F. temperatum Isolated from Maize in Argentina

Fusarium subglutinans and F. temperatum are two important fungal pathogens of maize whose distinctness as separate species has been difficult to assess. We isolated strains of these species from commercial and native maize varieties in Argentina and sequenced >28,000 loci to estimate genetic variation in the sample. Our objectives were to measure genetic divergence between the species, infer demographic parameters related to their split, and describe the population structure of the sample. When analyzed together, over 30% of each species' polymorphic sites (>2,500 sites) segregate as polymorphisms in the other. Demographic modeling confirmed the species split predated maize domestication, but subsequent between-species gene flow has occurred, with gene flow from F. subglutinans into F. temperatum greater than gene flow in the reverse direction. In F. subglutinans, little evidence exists for substructure or recent selective sweeps, but there is evidence for limited sexual reproduction. In F. temperatum, there is clear evidence for population substructure and signals of abundant recent selective sweeps, with sexual reproduction probably less common than in F. subglutinans. Both genetic variation and the relative number of polymorphisms shared between species increase near the telomeres of all 12 chromosomes, where genes related to plant-pathogen interactions often are located. Our results suggest that species boundaries between closely related Fusarium species can be semipermeable and merit further study. Such semipermeability could facilitate unanticipated genetic exchange between species and enable quicker permanent responses to changes in the agro-ecosystem, e.g., pathogen-resistant host varieties, new chemical and biological control agents, and agronomic practices.

Hidden paths to endless forms most wonderful: Complexity of bacterial motility shapes diversification of latent phenotypes

BACKGROUND

Evolution in one selective environment often latently generates phenotypic change that is manifested only later in different environments, but the complexity of behavior important to fitness in the original environment might influence the character of such latent-phenotype evolution. Using Myxococcus xanthus, a bacterium possessing two motility systems differing in effectiveness on hard vs. soft surfaces, we test (i) whether and how evolution while swarming on one surface-the selective surface-latently alters motility on the alternative surface type and (ii) whether patterns of such latent-phenotype evolution depend on the complexity of ancestral motility, specific ancestral motility genotypes and/or the selective surface of evolution. We analysze an experiment in which populations established from three ancestral genotypes-one with both motility systems intact and two others with one system debilitated-evolved while swarming across either hard or soft agar in six evolutionary treatments. We then compare motility-phenotype patterns across selective vs. alternative surface types.

RESULTS

Latent motility evolution was pervasive but varied in character as a function of the presence of one or two functional motility systems and, for some individual-treatment comparisons, the specific ancestral genotype and/or selective surface. Swarming rates on alternative vs. selective surfaces were positively correlated generally among populations with one functional motility system but not among those with two. This suggests that opportunities for pleiotropy and epistasis generated by increased genetic complexity underlying behavior can alter the character of latent-phenotype evolution. No tradeoff between motility performance across surface types was detected in the dual-system treatments, even after adaptation on a surface on which one motility system dominates strongly over the other in driving movement, but latent-phenotype evolution was instead idiosyncratic in these treatments. We further find that the magnitude of stochastic diversification at alternative-surface swarming among replicate populations greatly exceeded diversification of selective-surface swarming within some treatments and varied across treatments.

CONCLUSION

Collectively, our results suggest that increases in the genetic and mechanistic complexity of behavior can increase the complexity of latent-phenotype evolution outcomes and illustrate that diversification manifested during evolution in one environment can be augmented greatly by diversification of latent phenotypes manifested later.

Whole-genome sequencing of Tarim red deer (Cervus elaphus yarkandensis) reveals demographic history and adaptations to an arid-desert environment

Background

The initiation of desert conditions in the Tarim Basin in China since the late Miocene has led to the significant genetic structuring of local organisms. Tarim Red Deer (Cervus elaphus yarkandensis, TRD) have adapted to the harsh environmental conditions in this basin, including high solar radiation and temperature, aridity, and poor nutritional conditions. However, the underlying genetic basis of this adaptation is poorly understood.

Results

We sequenced the whole genomes of 13 TRD individuals, conducted comparative genomic analyses, and estimated demographic fluctuation. The ∂a∂i model estimated that the TRD and Tule elk (Cervus canadensis nannodes) populations diverged approximately 0.98 Mya. Analyses revealed a substantial influence of the Earth’s climate on the effective population size of TRD, associated with glacial advances and retreat, and human activities likely underlie a recent serious decline in population. A marked bottleneck may have profoundly affected the genetic diversity of TRD populations. We detected a set of candidate genes, pathways, and GO categories related to oxidative stress, water reabsorption, immune regulation, energy metabolism, eye protection, heat stress, respiratory system adaptation, prevention of high blood pressure, and DNA damage and repair that may directly or indirectly be involved in the adaptation of TRD to an arid-desert environment.

Conclusions

Our analyses highlight the role of historical global climates in the population dynamics of TRD. In light of ongoing global warming and the increasing incidence of droughts, our study offers insights into the genomic adaptations of animals, especially TRD, to extreme arid-desert environments and provides a valuable resource for future research on conservation design and biological adaptations to environmental change.

Background Selection Does Not Mimic the Patterns of Genetic Diversity Produced by Selective Sweeps

It is increasingly evident that natural selection plays a prominent role in shaping patterns of diversity across the genome. The most commonly studied modes of natural selection are positive selection and negative selection, which refer to directional selection for and against derived mutations, respectively. Positive selection can result in hitchhiking events, in which a beneficial allele rapidly replaces all others in the population, creating a valley of diversity around the selected site along with characteristic skews in allele frequencies and linkage disequilibrium among linked neutral polymorphisms. Similarly, negative selection reduces variation not only at selected sites but also at linked sites, a phenomenon called background selection (BGS). Thus, discriminating between these two forces may be difficult, and one might expect efforts to detect hitchhiking to produce an excess of false positives in regions affected by BGS. Here, we examine the similarity between BGS and hitchhiking models via simulation. First, we show that BGS may somewhat resemble hitchhiking in simplistic scenarios in which a region constrained by negative selection is flanked by large stretches of unconstrained sites, echoing previous results. However, this scenario does not mirror the actual spatial arrangement of selected sites across the genome. By performing forward simulations under more realistic scenarios of BGS, modeling the locations of protein-coding and conserved noncoding DNA in real genomes, we show that the spatial patterns of variation produced by BGS rarely mimic those of hitchhiking events. Indeed, BGS is not substantially more likely than neutrality to produce false signatures of hitchhiking. This holds for simulations modeled after both humans and Drosophila, and for several different demographic histories. These results demonstrate that appropriately designed scans for hitchhiking need not consider BGS's impact on false-positive rates. However, we do find evidence that BGS increases the false-negative rate for hitchhiking, an observation that demands further investigation.

Genome-Wide Association Study Reveals Novel Candidate Genes Associated with Productivity and Disease Resistance to Moniliophthora spp. in Cacao (Theobroma cacao L.)

Cacao (Theobroma cacao L.), the source of chocolate, is one of the most important commodity products worldwide that helps improve the economic livelihood of farmers. Diseases like frosty pod rot caused by Moniliophthora roreri and witches’ broom caused by Moniliophthora perniciosa limit the cacao productivity, this can be solved by using resistant varieties. In the current study, we sequenced 229 cacao accessions using genotyping-by-sequencing to examine the genetic diversity and population structure employing 9,003 and 8,131 single nucleotide polymorphisms recovered by mapping against two cacao genomes (Criollo B97-61/B2 v2 and Matina 1-6 v1.1). In the phenotypic evaluation, three promising accessions for productivity and 10 with good tolerance to the frosty pod rot and witches’ broom diseases were found. A genome-wide association study was performed on 102 accessions, discovering two genes associated with productivity and seven to disease resistance. The results enriched the knowledge of the genetic regions associated with important cacao traits that can have significant implications for conservation and breeding strategies like marker-assisted selection.

Toward an Evolutionarily Appropriate Null Model: Jointly Inferring Demography and Purifying Selection

The question of the relative evolutionary roles of adaptive and nonadaptive processes has been a central debate in population genetics for nearly a century. While advances have been made in the theoretical development of the underlying models, and statistical methods for estimating their parameters from large-scale genomic data, a framework for an appropriate null model remains elusive. A model incorporating evolutionary processes known to be in constant operation, genetic drift (as modulated by the demographic history of the population) and purifying selection, is lacking. Without such a null model, the role of adaptive processes in shaping within- and between-population variation may not be accurately assessed. Here, we investigate how population size changes and the strength of purifying selection affect patterns of variation at "neutral" sites near functional genomic components. We propose a novel statistical framework for jointly inferring the contribution of the relevant selective and demographic parameters. By means of extensive performance analyses, we quantify the utility of the approach, identify the most important statistics for parameter estimation, and compare the results with existing methods. Finally, we reanalyze genome-wide population-level data from a Zambian population of Drosophila melanogaster, and find that it has experienced a much slower rate of population growth than was inferred when the effects of purifying selection were neglected. Our approach represents an appropriate null model, against which the effects of positive selection can be assessed.

The genetic basis of hybrid male sterility in sympatric Primulina species

BACKGROUND

Sympatric sister species provide an opportunity to investigate the genetic mechanisms and evolutionary forces that maintain species boundaries. The persistence of morphologically and genetically distinct populations in sympatry can only occur if some degree of reproductive isolation exists. A pair of sympatric sister species of Primulina (P. depressa and P. danxiaensis) was used to explore the genetic architecture of hybrid male sterility.

RESULTS

We mapped one major- and seven minor-effect quantitative trait loci (QTLs) that underlie pollen fertility rate (PFR). These loci jointly explained 55.4% of the phenotypic variation in the F2 population. A Bateson-Dobzhansky-Muller (BDM) model involving three loci was observed in this system. We found genotypic correlations between hybrid male sterility and flower morphology, consistent with the weak but significant phenotypic correlations between PFR and floral traits.

CONCLUSIONS

Hybrid male sterility in Primulina is controlled by a polygenic genetic basis with a complex pattern. The genetic incompatibility involves a three-locus BDM model. Hybrid male sterility is genetically correlated with floral morphology and divergence hitchhiking may occur between them.

Detecting Positive Selection in Populations Using Genetic Data

High-throughput genomic sequencing allows to disentangle the evolutionary forces acting in populations. Among evolutionary forces, positive selection has received a lot of attention because it is related to the adaptation of populations in their environments, both biotic and abiotic. Positive selection, also known as Darwinian selection, occurs when an allele is favored by natural selection. The frequency of the favored allele increases in the population and, due to genetic hitchhiking, neighboring linked variation diminishes, creating so-called selective sweeps. Such a process leaves traces in genomes that can be detected in a future time point. Detecting traces of positive selection in genomes is achieved by searching for signatures introduced by selective sweeps, such as regions of reduced variation, a specific shift of the site frequency spectrum, and particular linkage disequilibrium (LD) patterns in the region. A variety of approaches can be used for detecting selective sweeps, ranging from simple implementations that compute summary statistics to more advanced statistical approaches, e.g., Bayesian approaches, maximum-likelihood-based methods, and machine learning methods. In this chapter, we discuss selective sweep detection methodologies on the basis of their capacity to analyze whole genomes or just subgenomic regions, and on the specific polymorphism patterns they exploit as selective sweep signatures. We also summarize the results of comparisons among five open-source software releases (SweeD, SweepFinder, SweepFinder2, OmegaPlus, and RAiSD) regarding sensitivity, specificity, and execution times. Furthermore, we test and discuss machine learning methods and present a thorough performance analysis. In equilibrium neutral models or mild bottlenecks, most methods are able to detect selective sweeps accurately. Methods and tools that rely on linkage disequilibrium (LD) rather than single SNPs exhibit higher true positive rates than the site frequency spectrum (SFS)-based methods under the model of a single sweep or recurrent hitchhiking. However, their false positive rate is elevated when a misspecified demographic model is used to build the distribution of the statistic under the null hypothesis. Both LD and SFS-based approaches suffer from decreased accuracy on localizing the true target of selection in bottleneck scenarios. Furthermore, we present an extensive analysis of the effects of gene flow on selective sweep detection, a problem that has been understudied in selective sweep literature.

Males, Outcrossing, and Sexual Selection in Caenorhabditis Nematodes

Males of Caenorhabditis elegans provide a crucial practical tool in the laboratory, but, as the rarer and more finicky sex, have not enjoyed the same depth of research attention as hermaphrodites. Males, however, have attracted the attention of evolutionary biologists who are exploiting the C. elegans system to test longstanding hypotheses about sexual selection, sexual conflict, transitions in reproductive mode, and genome evolution, as well as to make new discoveries about Caenorhabditis organismal biology. Here, we review the evolutionary concepts and data informed by study of males of C. elegans and other Caenorhabditis We give special attention to the important role of sperm cells as a mediator of inter-male competition and male-female conflict that has led to drastic trait divergence across species, despite exceptional phenotypic conservation in many other morphological features. We discuss the evolutionary forces important in the origins of reproductive mode transitions from males being common (gonochorism: females and males) to rare (androdioecy: hermaphrodites and males) and the factors that modulate male frequency in extant androdioecious populations, including the potential influence of selective interference, host-pathogen coevolution, and mutation accumulation. Further, we summarize the consequences of males being common vs rare for adaptation and for trait divergence, trait degradation, and trait dimorphism between the sexes, as well as for molecular evolution of the genome, at both micro-evolutionary and macro-evolutionary timescales. We conclude that C. elegans male biology remains underexploited and that future studies leveraging its extensive experimental resources are poised to discover novel biology and to inform profound questions about animal function and evolution.

Whole-Genome Resequencing Reveals Loci Associated With Thoracic Vertebrae Number in Sheep

The number of vertebrae, especially thoracic vertebrae, is an important economic trait that may influence carcass length and meat production in animals. However, the genetic basis of vertebrae number in sheep is still poorly understood. To detect the candidate genes, 400 increased number of thoracic vertebrae (T14L6) and 200 normal (T13L6) Kazakh sheep were collected. We generated and sequenced 60 pools of genomic DNA (each pool prepared by mixing genomic DNA from 10 sheep with the same thoracic traits), with an average depth of coverage of 25.65×. We identified a total of 42,075,402 SNPs and 11 putatively selected genomic regions, including the VRTN gene and the HoxA gene family that regulate vertebral development. The most prominent areas of selective elimination were located in a region of chromosome 7, including VRTN, which regulates spinal development and morphology. Further investigation indicated that the expression level of the VRTN gene during fetal development was significantly higher in sheep with more thoracic vertebrae than in those with a normal number of thoracic vertebrae. A genome-wide comparison between sheep with increased and normal numbers of thoracic vertebrae showed that the VRTN gene is the major selection locus for the number of thoracic vertebrae in sheep and has the potential to be utilized in sheep breeding in the future.

Understanding the Factors That Shape Patterns of Nucleotide Diversity in the House Mouse Genome

A major goal of population genetics has been to determine the extent by which selection at linked sites influences patterns of neutral nucleotide diversity in the genome. Multiple lines of evidence suggest that diversity is influenced by both positive and negative selection. For example, in many species there are troughs in diversity surrounding functional genomic elements, consistent with the action of either background selection (BGS) or selective sweeps. In this study, we investigated the causes of the diversity troughs that are observed in the wild house mouse genome. Using the unfolded site frequency spectrum, we estimated the strength and frequencies of deleterious and advantageous mutations occurring in different functional elements in the genome. We then used these estimates to parameterize forward-in-time simulations of chromosomes, using realistic distributions of functional elements and recombination rate variation in order to determine whether selection at linked sites can explain the observed patterns of nucleotide diversity. The simulations suggest that BGS alone cannot explain the dips in diversity around either exons or conserved noncoding elements. A combination of BGS and selective sweeps produces deeper dips in diversity than BGS alone, but the inferred parameters of selection cannot fully explain the patterns observed in the genome. Our results provide evidence of sweeps shaping patterns of nucleotide diversity across the mouse genome and also suggest that infrequent, strongly advantageous mutations play an important role in this. The limitations of using the unfolded site frequency spectrum for inferring the frequency and effects of advantageous mutations are discussed.

The Effects on Neutral Variability of Recurrent Selective Sweeps and Background Selection

Levels of variability and rates of adaptive evolution may be affected by hitchhiking, the effect of selection on evolution at linked sites. Hitchhiking can be caused either by "selective sweeps" or by background selection, involving the spread of new favorable alleles or the elimination of deleterious mutations, respectively. Recent analyses of population genomic data have fitted models where both these processes act simultaneously, to infer the parameters of selection. Here, we investigate the consequences of relaxing a key assumption of some of these studies, that the time occupied by a selective sweep is negligible compared with the neutral coalescent time. We derive a new expression for the expected level of neutral variability in the presence of recurrent selective sweeps and background selection. We also derive approximate integral expressions for the effects of recurrent selective sweeps. The accuracy of the theoretical predictions was tested against multilocus simulations, with selection, recombination, and mutation parameters that are realistic for Drosophila melanogaster In the presence of crossing over, there is approximate agreement between the theoretical and simulation results. We show that the observed relationships between the rate of crossing over, and the level of synonymous site diversity and rate of adaptive evolution in Drosophila are probably mainly caused by background selection, whereas selective sweeps and population size changes are needed to produce the observed distortions of the site frequency spectrum.

Low genetic variation is associated with low mutation rate in the giant duckweed

Mutation rate and effective population size (N_e) jointly determine intraspecific genetic diversity, but the role of mutation rate is often ignored. Here we investigate genetic diversity, spontaneous mutation rate and N_e in the giant duckweed (Spirodela polyrhiza). Despite its large census population size, whole-genome sequencing of 68 globally sampled individuals reveals extremely low intraspecific genetic diversity. Assessed under natural conditions, the genome-wide spontaneous mutation rate is at least seven times lower than estimates made for other multicellular eukaryotes, whereas N_e is large. These results demonstrate that low genetic diversity can be associated with large-N_e species, where selection can reduce mutation rates to very low levels. This study also highlights that accurate estimates of mutation rate can help to explain seemingly unexpected patterns of genome-wide variation.

While the role of effective population size (N_e) in explaining variation in genetic diversity has received much attention, the role of spontaneous mutation rate is largely ignored. Here, Xu et al. show that giant duckweed has a high N_e yet low genetic diversity, likely due to its low mutation rate.

How Linked Selection Shapes the Diversity Landscape in Ficedula Flycatchers

There is an increasing awareness that selection affecting linked neutral sites strongly influences on how diversity is distributed across the genome. In particular, linked selection is likely involved in the formation of heterogenous landscapes of genetic diversity, including genomic regions with locally reduced effective population sizes that manifest as dips in diversity, and "islands" of differentiation between closely related populations or species. Linked selection can be in the form of background selection or selective sweeps, and a long-standing quest in population genetics has been to unveil the relative importance of these processes. Here, we analyzed the theoretically expected reduction of diversity caused by linked selection in the collared flycatcher (Ficedula albicollis) genome and compared this with population genomic data on the distribution of diversity across the flycatcher genome. By incorporating data on recombination rate variation and the density of target sites for selection (including both protein-coding genes and conserved noncoding elements), we found that background selection can explain most of the observed baseline variation in genetic diversity. However, positive selection was necessary to explain the pronounced local diversity dips in the collared flycatcher genome. We confirmed our analytical findings by comprehensive simulations. Therefore, our study demonstrates that even though both background selection and selective sweeps contribute to the heterogeneous diversity landscape seen in this avian system, they play different roles in shaping it.

Selective Sweeps

For almost 20 years, many inference methods have been developed to detect selective sweeps and localize the targets of directional selection in the genome. These methods are based on population genetic models that describe the effect of a beneficial allele (e.g., a new mutation) on linked neutral variation (driven by directional selection from a single copy to fixation). Here, I discuss these models, ranging from selective sweeps in a panmictic population of constant size to evolutionary traffic when simultaneous sweeps at multiple loci interfere, and emphasize the important role of demography and population structure in data analysis. In the past 10 years, soft sweeps that may arise after an environmental change from directional selection on standing variation have become a focus of population genetic research. In contrast to selective sweeps, they are caused by beneficial alleles that were neutrally segregating in a population before the environmental change or were present at a mutation-selection balance in appreciable frequency.

The importance of the Neutral Theory in 1968 and 50 years on: A response to Kern and Hahn 2018

A recent article reassessing the Neutral Theory of Molecular Evolution claims that it is no longer as important as is widely believed. The authors argue that "the neutral theory was supported by unreliable theoretical and empirical evidence from the beginning, and that in light of modern, genome-scale data, we can firmly reject its universality." Claiming that "the neutral theory has been overwhelmingly rejected," they propose instead that natural selection is the major force shaping both between-species divergence and within-species variation. Although this is probably a minority view, it is important to evaluate such claims carefully in the context of current knowledge, as inaccuracies can sometimes morph into an accepted narrative for those not familiar with the underlying science. We here critically examine and ultimately reject Kern and Hahn's arguments and assessment, and instead propose that it is now abundantly clear that the foundational ideas presented five decades ago by Kimura and Ohta are indeed correct.

Background selection and biased gene conversion affect more than 95% of the human genome and bias demographic inferences

Disentangling the effect on genomic diversity of natural selection from that of demography is notoriously difficult, but necessary to properly reconstruct the history of species. Here, we use high-quality human genomic data to show that purifying selection at linked sites (i.e. background selection, BGS) and GC-biased gene conversion (gBGC) together affect as much as 95% of the variants of our genome. We find that the magnitude and relative importance of BGS and gBGC are largely determined by variation in recombination rate and base composition. Importantly, synonymous sites and non-transcribed regions are also affected, albeit to different degrees. Their use for demographic inference can lead to strong biases. However, by conditioning on genomic regions with recombination rates above 1.5 cM/Mb and mutation types (C↔G, A↔T), we identify a set of SNPs that is mostly unaffected by BGS or gBGC, and that avoids these biases in the reconstruction of human history.

RAiSD detects positive selection based on multiple signatures of a selective sweep and SNP vectors

Selective sweeps leave distinct signatures locally in genomes, enabling the detection of loci that have undergone recent positive selection. Multiple signatures of a selective sweep are known, yet each neutrality test only identifies a single signature. We present RAiSD (Raised Accuracy in Sweep Detection), an open-source software that implements a novel, to our knowledge, and parameter-free detection mechanism that relies on multiple signatures of a selective sweep via the enumeration of SNP vectors. RAiSD achieves higher sensitivity and accuracy than the current state of the art, while the computational complexity is greatly reduced, allowing up to 1000 times faster processing than widely used tools, and negligible memory requirements.

Nikolaos Alachiotis and Pavlos Pavlidis present RAiSD, a computational method for identifying multiple signatures of selective sweeps using single nucleotide polymorphism vectors. They show that RAiSD has higher sensitivity and accuracy with reduced computational complexity than current methods.

Human demographic history has amplified the effects of background selection across the genome

Natural populations often grow, shrink, and migrate over time. Such demographic processes can affect genome-wide levels of genetic diversity. Additionally, genetic variation in functional regions of the genome can be altered by natural selection, which drives adaptive mutations to higher frequencies or purges deleterious ones. Such selective processes affect not only the sites directly under selection but also nearby neutral variation through genetic linkage via processes referred to as genetic hitchhiking in the context of positive selection and background selection (BGS) in the context of purifying selection. While there is extensive literature examining the consequences of selection at linked sites at demographic equilibrium, less is known about how non-equilibrium demographic processes influence the effects of hitchhiking and BGS. Utilizing a global sample of human whole-genome sequences from the Thousand Genomes Project and extensive simulations, we investigate how non-equilibrium demographic processes magnify and dampen the consequences of selection at linked sites across the human genome. When binning the genome by inferred strength of BGS, we observe that, compared to Africans, non-African populations have experienced larger proportional decreases in neutral genetic diversity in strong BGS regions. We replicate these findings in admixed populations by showing that non-African ancestral components of the genome have also been affected more severely in these regions. We attribute these differences to the strong, sustained/recurrent population bottlenecks that non-Africans experienced as they migrated out of Africa and throughout the globe. Furthermore, we observe a strong correlation between FST and the inferred strength of BGS, suggesting a stronger rate of genetic drift. Forward simulations of human demographic history with a model of BGS support these observations. Our results show that non-equilibrium demography significantly alters the consequences of selection at linked sites and support the need for more work investigating the dynamic process of multiple evolutionary forces operating in concert.

Concurrent action of purifying selection and gene conversion results in extreme conservation of the major stress-inducible Hsp70 genes in mammals

Several evolutionary mechanisms alter the fate of mutations and genes within populations based on their exhibited functional effects. To understand the underlying mechanisms involved in the evolution of the cellular stress response, a very conserved mechanism in the course of organismal evolution, we studied the patterns of natural genetic variation and functional consequences of polymorphisms of two stress-inducible Hsp70 genes. These genes, HSPA1A and HSPA1B, are major orchestrators of the cellular stress response and are associated with several human diseases. Our phylogenetic analyses revealed that the duplication of HSPA1A and HSPA1B originated in a lineage proceeding to placental mammals, and henceforth they remained in conserved synteny. Additionally, analyses of synonymous and non-synonymous changes suggest that purifying selection shaped the HSPA1 gene diversification, while gene conversion resulted in high sequence conservation within species. In the human HSPA1-cluster, the vast majority of mutations are synonymous and specific genic regions are devoid of mutations. Furthermore, functional characterization of several human polymorphisms revealed subtle differences in HSPA1A stability and intracellular localization. Collectively, the observable patterns of HSPA1A-1B variation describe an evolutionary pattern, in which purifying selection and gene conversion act simultaneously and conserve a major orchestrator of the cellular stress response.

Genomewide patterns of variation in genetic diversity are shared among populations, species and higher-order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e ) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage-specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e . Utilizing population-level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population-scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST , PBS, d_xy ) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.