Identifying loci under positive selection in complex population histories

Leveraging graphical model techniques to study evolution on phylogenetic networks

The evolution of molecular and phenotypic traits is commonly modelled using Markov processes along a phylogeny. This phylogeny can be a tree, or a network if it includes reticulations, representing events such as hybridization or admixture. Computing the likelihood of data observed at the leaves is costly as the size and complexity of the phylogeny grows. Efficient algorithms exist for trees, but cannot be applied to networks. We show that a vast array of models for trait evolution along phylogenetic networks can be reformulated as graphical models, for which efficient belief propagation algorithms exist. We provide a brief review of belief propagation on general graphical models, then focus on linear Gaussian models for continuous traits. We show how belief propagation techniques can be applied for exact or approximate (but more scalable) likelihood and gradient calculations, and prove novel results for efficient parameter inference of some models. We highlight the possible fruitful interactions between graphical models and phylogenetic methods. For example, approximate likelihood approaches have the potential to greatly reduce computational costs for phylogenies with reticulations.This article is part of the theme issue '"A mathematical theory of evolution": phylogenetic models dating back 100 years'.

Identifying Significant SNPs of the Total Number of Piglets Born and Their Relationship with Leg Bumps in Pigs

The aim of this study was to identify genetic variants and pathways associated with the total number of piglets born and to investigate the potential negative consequences of the intensive selection for reproductive traits, particularly the formation of bumps on the legs of pigs. We used genome-wide association analysis and methods for identifying selection signatures. As a result, 47 SNPs were identified, localized in genes that play a significant role during sow pregnancy. These genes are involved in follicle growth and development (SGC), early embryonic development (CCDC3, LRRC8C, LRFN3, TNFRSF19), endometrial receptivity and implantation (NEBL), placentation, and embryonic development (ESRRG, GHRHR, TUSC3, NBAS). Several genes are associated with disorders of the nervous system and brain development (BCL11B, CDNF, ULK4, CC2D2A, KCNK2). Additionally, six SNPs are associated with the formation of bumps on the legs of pigs. These variants include intronic variants in the CCDC3, ULK4, and MINDY4 genes, as well as intergenic variants, regulatory region variants, and variants in the exons of non-coding transcripts. The results suggest important biological pathways and genetic variants associated with sow fertility and highlight the potential negative impacts on the health and physical condition of pigs.

Deep genetic substructure within bonobos

Establishing the genetic and geographic structure of populations is fundamental, both to understand their evolutionary past and preserve their future. Nevertheless, the patterns of genetic population structure are unknown for most endangered species. This is the case for bonobos (Pan paniscus), which, together with chimpanzees (Pan troglodytes), are humans' closest living relatives. Chimpanzees live across equatorial Africa and are classified into four subspecies,1 with some genetic population substructure even within subspecies. Conversely, bonobos live exclusively in the Democratic Republic of Congo and are considered a homogeneous group with low genetic diversity,2 despite some population structure inferred from mtDNA. Nevertheless, mtDNA aside, their genetic structure remains unknown, hampering our understanding of the species and conservation efforts. Mapping bonobo genetic diversity in space is, however, challenging because, being endangered, only non-invasive sampling is possible for wild individuals. Here, we jointly analyze the exomes and mtDNA from 20 wild-born bonobos, the whole genomes of 10 captive bonobos, and the mtDNA of 136 wild individuals. We identify three genetically distinct bonobo groups of inferred Central, Western, and Far-Western geographic origin within the bonobo range. We estimate the split time between the central and western populations to be ∼145,000 years ago and genetic differentiation to be in the order of that of the closest chimpanzee subspecies. Furthermore, our estimated long-term Ne for Far-West (∼3,000) is among the lowest estimated for any great ape lineage. Our results highlight the need to attend to the bonobo substructure, both in terms of research and conservation.

Adaptive Evolution of Two Distinct Adaptive Haplotypes of Neanderthal Origin at the Immunoglobulin Heavy-chain Locus in East Asian and European Populations

Immunoglobulins (Igs) have a crucial role in humoral immunity. Two recent studies have reported a high-frequency Neanderthal-introgressed haplotype throughout Eurasia and a high-frequency Neanderthal-introgressed haplotype specific to southern East Asia at the immunoglobulin heavy-chain (IGH) gene locus on chromosome 14q32.33. Surprisingly, we found the previously reported high-frequency Neanderthal-introgressed haplotype does not exist throughout Eurasia. Instead, our study identified two distinct high-frequency haplotypes of putative Neanderthal origin in East Asia and Europe, although they shared introgressed alleles. Notably, the alleles of putative Neanderthal origin reduced the expression of IGHG1 and increased the expression of IGHG2 and IGHG3 in various tissues. These putatively introgressed alleles also affected the production of IgG1 upon antigen stimulation and increased the risk of systemic lupus erythematosus. Additionally, the greatest genetic differentiation across the whole genome between southern and northern East Asians was observed for the East Asian haplotype of putative Neanderthal origin. The frequency decreased from southern to northern East Asia and correlated positively with the genome-wide proportion of southern East Asian ancestry, indicating that this putative positive selection likely occurred in the common ancestor of southern East Asian populations before the admixture with northern East Asian populations.

Genomic evidence for human-mediated introgressive hybridization and selection in the developed breed

BACKGROUND

The pig (Sus Scrofa) is one of the oldest domesticated livestock species that has undergone extensive improvement through modern breeding. European breeds have advantages in lean meat development and highly-productive body type, whereas Asian breeds possess extraordinary fat deposition and reproductive performance. Consequently, Eurasian breeds have been extensively used to develop modern commercial breeds for fast-growing and high prolificacy. However, limited by the sequencing technology, the genome architecture of some nascent developed breeds and the human-mediated impact on their genomes are still unknown.

RESULTS

Through whole-genome analysis of 178 individuals from an Asian locally developed pig breed, Beijing Black pig, and its two ancestors from two different continents, we found the pervasive inconsistent gene trees and species trees across the genome of Beijing Black pig, which suggests its introgressive hybrid origin. Interestingly, we discovered that this developed breed has more genetic relationships with European pigs and an unexpected introgression from Asian pigs to this breed, which indicated that human-mediated introgression could form the porcine genome architecture in a completely different type compared to native introgression. We identified 554 genomic regions occupied 63.30 Mb with signals of introgression from the Asian ancestry to Beijing Black pig, and the genes in these regions enriched in pathways associated with meat quality, fertility, and disease-resistant. Additionally, a proportion of 7.77% of genomic regions were recognized as regions that have been under selection. Moreover, combined with the results of a genome-wide association study for meat quality traits in the 1537 Beijing Black pig population, two important candidate genes related to meat quality traits were identified. DNAJC6 is related to intramuscular fat content and fat deposition, and RUFY4 is related to meat pH and tenderness.

CONCLUSIONS

Our research provides insight for analyzing the origins of nascent developed breeds and genome-wide selection remaining in the developed breeds mediated by humans during modern breeding.

Genome-wide scans for selective sweeps using convolutional neural networks

MOTIVATION

Recent methods for selective sweep detection cast the problem as a classification task and use summary statistics as features to capture region characteristics that are indicative of a selective sweep, thereby being sensitive to confounding factors. Furthermore, they are not designed to perform whole-genome scans or to estimate the extent of the genomic region that was affected by positive selection; both are required for identifying candidate genes and the time and strength of selection.

RESULTS

We present ASDEC (https://github.com/pephco/ASDEC), a neural-network-based framework that can scan whole genomes for selective sweeps. ASDEC achieves similar classification performance to other convolutional neural network-based classifiers that rely on summary statistics, but it is trained 10× faster and classifies genomic regions 5× faster by inferring region characteristics from the raw sequence data directly. Deploying ASDEC for genomic scans achieved up to 15.2× higher sensitivity, 19.4× higher success rates, and 4× higher detection accuracy than state-of-the-art methods. We used ASDEC to scan human chromosome 1 of the Yoruba population (1000Genomes project), identifying nine known candidate genes.

A graph-based genome and pan-genome variation of the model plant Setaria

Setaria italica (foxtail millet), a founder crop of East Asian agriculture, is a model plant for C4 photosynthesis and developing approaches to adaptive breeding across multiple climates. Here we established the Setaria pan-genome by assembling 110 representative genomes from a worldwide collection. The pan-genome is composed of 73,528 gene families, of which 23.8%, 42.9%, 29.4% and 3.9% are core, soft core, dispensable and private genes, respectively; 202,884 nonredundant structural variants were also detected. The characterization of pan-genomic variants suggests their importance during foxtail millet domestication and improvement, as exemplified by the identification of the yield gene SiGW3, where a 366-bp presence/absence promoter variant accompanies gene expression variation. We developed a graph-based genome and performed large-scale genetic studies for 68 traits across 13 environments, identifying potential genes for millet improvement at different geographic sites. These can be used in marker-assisted breeding, genomic selection and genome editing to accelerate crop improvement under different climatic conditions.

Sweepstakes reproductive success via pervasive and recurrent selective sweeps

Highly fecund natural populations characterized by high early mortality abound, yet our knowledge about their recruitment dynamics is somewhat rudimentary. This knowledge gap has implications for our understanding of genetic variation, population connectivity, local adaptation, and the resilience of highly fecund populations. The concept of sweepstakes reproductive success, which posits a considerable variance and skew in individual reproductive output, is key to understanding the distribution of individual reproductive success. However, it still needs to be determined whether highly fecund organisms reproduce through sweepstakes and, if they do, the relative roles of neutral and selective sweepstakes. Here, we use coalescent-based statistical analysis of population genomic data to show that selective sweepstakes likely explain recruitment dynamics in the highly fecund Atlantic cod. We show that the Kingman coalescent (modelling no sweepstakes) and the Xi-Beta coalescent (modelling random sweepstakes), including complex demography and background selection, do not provide an adequate fit for the data. The Durrett-Schweinsberg coalescent, in which selective sweepstakes result from recurrent and pervasive selective sweeps of new mutations, offers greater explanatory power. Our results show that models of sweepstakes reproduction and multiple-merger coalescents are relevant and necessary for understanding genetic diversity in highly fecund natural populations. These findings have fundamental implications for understanding the recruitment variation of fish stocks and general evolutionary genomics of high-fecundity organisms.

Locus-level antagonistic selection shaped the polygenic architecture of human complex diseases

BACKGROUND

We aimed to evaluate the potential role of antagonistic selection in polygenic diseases: if one variant increases the risk of one disease and decreases the risk of another disease, the signals of genetic risk elimination by natural selection will be distorted, which leads to a higher frequency of risk alleles.

METHODS

We applied local genetic correlations and transcriptome-wide association studies to identify genomic loci and genes adversely associated with at least two diseases. Then, we used different population genetic metrics to measure the signals of natural selection for these loci and genes.

RESULTS

First, we identified 2120 cases of antagonistic pleiotropy (negative local genetic correlation) among 87 diseases in 716 genomic loci (antagonistic loci). Next, by comparing with non-antagonistic loci, we observed that antagonistic loci explained an excess proportion of disease heritability (median 6%), showed enhanced signals of balancing selection, and reduced signals of directional polygenic adaptation. Then, at the gene expression level, we identified 31,991 cases of antagonistic pleiotropy among 98 diseases at 4368 genes. However, evidence of altered signals of selection pressure and heritability distribution at the gene expression level is limited.

CONCLUSION

We conclude that antagonistic pleiotropy is widespread among human polygenic diseases, and it has distorted the evolutionary signal and genetic architecture of diseases at the locus level.

The immunogenetic impact of European colonization in the Americas

The introduction of pathogens originating from Eurasia into the Americas during early European contact has been associated with high mortality rates among Indigenous peoples, likely contributing to their historical and precipitous population decline. However, the biological impacts of imported infectious diseases and resulting epidemics, especially in terms of pathogenic effects on the Indigenous immunity, remain poorly understood and highly contentious to this day. Here, we examine multidisciplinary evidence underpinning colonization-related immune genetic change, providing contextualization from anthropological studies, paleomicrobiological evidence of contrasting host-pathogen coevolutionary histories, and the timings of disease emergence. We further summarize current studies examining genetic signals reflecting post-contact Indigenous population bottlenecks, admixture with European and other populations, and the putative effects of natural selection, with a focus on ancient DNA studies and immunity-related findings. Considering current genetic evidence, together with a population genetics theoretical approach, we show that post-contact Indigenous immune adaptation, possibly influenced by selection exerted by introduced pathogens, is highly complex and likely to be affected by multifactorial causes. Disentangling putative adaptive signals from those of genetic drift thus remains a significant challenge, highlighting the need for the implementation of population genetic approaches that model the short time spans and complex demographic histories under consideration. This review adds to current understandings of post-contact immunity evolution in Indigenous peoples of America, with important implications for bettering our understanding of human adaptation in the face of emerging infectious diseases.

Disentangling Signatures of Selection Before and After European Colonization in Latin Americans

Throughout human evolutionary history, large-scale migrations have led to intermixing (i.e., admixture) between previously separated human groups. Although classical and recent work have shown that studying admixture can yield novel historical insights, the extent to which this process contributed to adaptation remains underexplored. Here, we introduce a novel statistical model, specific to admixed populations, that identifies loci under selection while determining whether the selection likely occurred post-admixture or prior to admixture in one of the ancestral source populations. Through extensive simulations, we show that this method is able to detect selection, even in recently formed admixed populations, and to accurately differentiate between selection occurring in the ancestral or admixed population. We apply this method to genome-wide SNP data of ∼4,000 individuals in five admixed Latin American cohorts from Brazil, Chile, Colombia, Mexico, and Peru. Our approach replicates previous reports of selection in the human leukocyte antigen region that are consistent with selection post-admixture. We also report novel signals of selection in genomic regions spanning 47 genes, reinforcing many of these signals with an alternative, commonly used local-ancestry-inference approach. These signals include several genes involved in immunity, which may reflect responses to endemic pathogens of the Americas and to the challenge of infectious disease brought by European contact. In addition, some of the strongest signals inferred to be under selection in the Native American ancestral groups of modern Latin Americans overlap with genes implicated in energy metabolism phenotypes, plausibly reflecting adaptations to novel dietary sources available in the Americas.

The genomic signatures of natural selection in admixed human populations

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

A review of inbreeding depression in dairy cattle: current status, emerging control strategies, and future prospects

Dairy cattle breeding has historically focused on relatively small numbers of elite bulls as sires of sons. In recent years, even if generation intervals were reduced and more diverse sires of sons could have been selected, genomic selection has not fundamentally changed the fact that a large number of individuals are being analyzed. However, a relatively small number of elite bulls are still siring those animals. Therefore inbreeding-derived negative consequences in the gene pool have brought concern. The detrimental effects of non-additive genetic changes such as inbreeding depression and dominance have been widely disseminated while seriously affecting bioeconomically important parameters because of an antagonistic relationship between dairy production and reproductive traits. Therefore, the estimation of benefits and limitations of inbreeding and variance of the selection response deserves to be evaluated and discussed to preserve genetic variability, a significant concern in the selection of individuals for reproduction and production. Short-term strategies for genetic merit improvement through modern breeding programs have severely lowered high-producing dairy cattle fertility potential. Since the current selection programs potentially increase long-term costs, genetic diversity has decreased globally as a consequence. Therefore, a greater understanding of the potential that selection programs have for supporting long-term genetic sustainability and genetic diversity among dairy cattle populations should be prioritized in managing farm profitability. The present review provides a broad approach to current inbreeding-derived problems, identifying critical points to be solved and possible alternative strategies to control selection against homozygous haplotypes while maintaining sustained selection pressure. Moreover, this manuscript explores future perspectives, emphasizing theoretical applications and critical points, and strategies to avoid the adverse effects of inbreeding in dairy cattle. Finally, this review provides an overview of challenges that will soon require multidisciplinary approaches to managing dairy cattle populations, intending to combine increases in productive trait phenotypes with improvements in reproductive, health, welfare, linear conformation, and adaptability traits into the foreseeable future.

Detecting Selection in Multiple Populations by Modeling Ancestral Admixture Components

One of the most powerful and commonly used approaches for detecting local adaptation in the genome is the identification of extreme allele frequency differences between populations. In this article, we present a new maximum likelihood method for finding regions under positive selection. It is based on a Gaussian approximation to allele frequency changes and it incorporates admixture between populations. The method can analyze multiple populations simultaneously and retains power to detect selection signatures specific to ancestry components that are not representative of any extant populations. Using simulated data, we compare our method to related approaches, and show that it is orders of magnitude faster than the state-of-the-art, while retaining similar or higher power for most simulation scenarios. We also apply it to human genomic data and identify loci with extreme genetic differentiation between major geographic groups. Many of the genes identified are previously known selected loci relating to hair pigmentation and morphology, skin, and eye pigmentation. We also identify new candidate regions, including various selected loci in the Native American component of admixed Mexican-Americans. These involve diverse biological functions, such as immunity, fat distribution, food intake, vision, and hair development.

Sporadic occurrence of recent selective sweeps from standing variation in humans as revealed by an approximate Bayesian computation approach

During their dispersals over the last 100,000 years, modern humans have been exposed to a large variety of environments, resulting in genetic adaptation. While genome-wide scans for the footprints of positive Darwinian selection have increased knowledge of genes and functions potentially involved in human local adaptation, they have globally produced evidence of a limited contribution of selective sweeps in humans. Conversely, studies based on machine learning algorithms suggest that recent sweeps from standing variation are widespread in humans, an observation that has been recently questioned. Here, we sought to formally quantify the number of recent selective sweeps in humans, by leveraging approximate Bayesian computation and whole-genome sequence data. Our computer simulations revealed suitable ABC estimations, regardless of the frequency of the selected alleles at the onset of selection and the completion of sweeps. Under a model of recent selection from standing variation, we inferred that an average of 68 (from 56 to 79) and 140 (from 94 to 198) sweeps occurred over the last 100,000 years of human history, in African and Eurasian populations, respectively. The former estimation is compatible with human adaptation rates estimated since divergence with chimps, and reveals numbers of sweeps per generation per site in the range of values estimated in Drosophila. Our results confirm the rarity of selective sweeps in humans and show a low contribution of sweeps from standing variation to recent human adaptation.

Genome analysis traces regional dispersal of rice in Taiwan and Southeast Asia

The dispersal of rice (Oryza sativa) following domestication influenced massive social and cultural changes across South, East, and Southeast Asia. The history of dispersal across islands of Southeast Asia, and the role of Taiwan and the Austronesian expansion in this process remain largely unresolved. Here, we reconstructed the routes of dispersal of O. sativa ssp. japonica rice through Taiwan and the northern Philippines using whole-genome re-sequencing of indigenous rice landraces coupled with archaeological and paleoclimate data. Our results indicate that japonica rice found in the northern Philippines diverged from Indonesian landraces as early as 3500 BP. In contrast, rice cultivated by the indigenous peoples of the Taiwanese mountains has complex origins. It comprises two distinct populations, each best explained as a result of admixture between temperate japonica that presumably came from northeast Asia, and tropical japonica from the northern Philippines and mainland Southeast Asia respectively. We find that the temperate japonica component of these indigenous Taiwan populations diverged from northeast Asia subpopulations at about 2600 BP, while gene flow from the northern Philippines occurred before ∼1300 years BP. This coincides with a period of intensified trade established across the South China Sea. Finally, we find evidence for positive selection acting on distinct genomic regions in different rice subpopulations, indicating local adaptation associated with the spread of japonica rice.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.

An overview of current population genomics methods for the analysis of whole-genome resequencing data in eukaryotes

Characterizing the population history of a species and identifying loci underlying local adaptation is crucial in functional ecology, evolutionary biology, conservation and agronomy. The constant improvement of high-throughput sequencing techniques has facilitated the production of whole genome data in a wide range of species. Population genomics now provides tools to better integrate selection into a historical framework, and take into account selection when reconstructing demographic history. However, this improvement has come with a profusion of analytical tools that can confuse and discourage users. Such confusion limits the amount of information effectively retrieved from complex genomic data sets, and impairs the diffusion of the most recent analytical tools into fields such as conservation biology. It may also lead to redundancy among methods. To address these isssues, we propose an overview of more than 100 state-of-the-art methods that can deal with whole genome data. We summarize the strategies they use to infer demographic history and selection, and discuss some of their limitations. A website listing these methods is available at www.methodspopgen.com.

The timing of human adaptation from Neanderthal introgression

Admixture has the potential to facilitate adaptation by providing alleles that are immediately adaptive in a new environment or by simply increasing the long-term reservoir of genetic diversity for future adaptation. A growing number of cases of adaptive introgression are being identified in species across the tree of life, however the timing of selection, and therefore the importance of the different evolutionary roles of admixture, is typically unknown. Here, we investigate the spatio-temporal history of selection favoring Neanderthal-introgressed alleles in modern human populations. Using both ancient and present-day samples of modern humans, we integrate the known demographic history of populations, namely population divergence and migration, with tests for selection. We model how a sweep placed along different branches of an admixture graph acts to modify the variance and covariance in neutral allele frequencies among populations at linked loci. Using a method based on this model of allele frequencies, we study previously identified cases of adaptive Neanderthal introgression. From these, we identify cases in which Neanderthal-introgressed alleles were quickly beneficial and other cases in which they persisted at low frequency for some time. For some of the alleles that persisted at low frequency, we show that selection likely independently favored them later on in geographically separated populations. Our work highlights how admixture with ancient hominins has contributed to modern human adaptation and contextualizes observed levels of Neanderthal ancestry in present-day and ancient samples.

Evolutionary genetics of skin pigmentation in African populations

Skin color is a highly heritable human trait, and global variation in skin pigmentation has been shaped by natural selection, migration and admixture. Ethnically diverse African populations harbor extremely high levels of genetic and phenotypic diversity, and skin pigmentation varies widely across Africa. Recent genome-wide genetic studies of skin pigmentation in African populations have advanced our understanding of pigmentation biology and human evolutionary history. For example, novel roles in skin pigmentation for loci near MFSD12 and DDB1 have recently been identified in African populations. However, due to an underrepresentation of Africans in human genetic studies, there is still much to learn about the evolutionary genetics of skin pigmentation. Here, we summarize recent progress in skin pigmentation genetics in Africans and discuss the importance of including more ethnically diverse African populations in future genetic studies. In addition, we discuss methods for functional validation of adaptive variants related to skin pigmentation.

A Whole-Genome Scan for Association with Invasion Success in the Fruit Fly Drosophila suzukii Using Contrasts of Allele Frequencies Corrected for Population Structure

Evidence is accumulating that evolutionary changes are not only common during biological invasions but may also contribute directly to invasion success. The genomic basis of such changes is still largely unexplored. Yet, understanding the genomic response to invasion may help to predict the conditions under which invasiveness can be enhanced or suppressed. Here, we characterized the genome response of the spotted wing drosophila Drosophila suzukii during the worldwide invasion of this pest insect species, by conducting a genome-wide association study to identify genes involved in adaptive processes during invasion. Genomic data from 22 population samples were analyzed to detect genetic variants associated with the status (invasive versus native) of the sampled populations based on a newly developed statistic, we called C2, that contrasts allele frequencies corrected for population structure. We evaluated this new statistical framework using simulated data sets and implemented it in an upgraded version of the program BayPass. We identified a relatively small set of single-nucleotide polymorphisms that show a highly significant association with the invasive status of D. suzukii populations. In particular, two genes, RhoGEF64C and cpo, contained single-nucleotide polymorphisms significantly associated with the invasive status in the two separate main invasion routes of D. suzukii. Our methodological approaches can be applied to any other invasive species, and more generally to any evolutionary model for species characterized by nonequilibrium demographic conditions for which binary covariables of interest can be defined at the population level.

Genomic Approaches Reveal Pleiotropic Effects in Crossbred Beef Cattle

Carcass and meat quality are two important attributes for the beef industry because they drive profitability and consumer demand. These traits are of even greater importance in crossbred cattle used in subtropical and tropical regions for their superior adaptability because they tend to underperform compared to their purebred counterparts. Many of these traits are challenging and expensive to measure and unavailable until late in life or after the animal is harvested, hence unrealistic to improve through traditional phenotypic selection, but perfect candidates for genomic selection. Before genomic selection can be implemented in crossbred populations, it is important to explore if pleiotropic effects exist between carcass and meat quality traits. Therefore, the objective of this study was to identify genomic regions with pleiotropic effects on carcass and meat quality traits in a multibreed Angus-Brahman population that included purebred and crossbred animals. Data included phenotypes for 10 carcass and meat quality traits from 2,384 steers, of which 1,038 were genotyped with the GGP Bovine F-250. Single-trait genome-wide association studies were first used to investigate the relevance of direct additive genetic effects on each carcass, sensory and visual meat quality traits. A second analysis for each trait included all other phenotypes as covariates to correct for direct causal effects from identified genomic regions with pure direct effects on the trait under analysis. Five genomic windows on chromosomes BTA5, BTA7, BTA18, and BTA29 explained more than 1% of additive genetic variance of two or more traits. Moreover, three suggestive pleiotropic regions were identified on BTA10 and BTA19. The 317 genes uncovered in pleiotropic regions included anchoring and cytoskeletal proteins, key players in cell growth, muscle development, lipid metabolism and fat deposition, and important factors in muscle proteolysis. A functional analysis of these genes revealed GO terms directly related to carcass quality, meat quality, and tenderness in beef cattle, including calcium-related processes, cell signaling, and modulation of cell-cell adhesion. These results contribute with novel information about the complex genetic architecture and pleiotropic effects of carcass and meat quality traits in crossbred beef cattle.

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

Widespread haploid-biased gene expression enables sperm-level natural selection

Sperm are haploid but must be functionally equivalent to distribute alleles equally among progeny. Accordingly, gene products are shared through spermatid cytoplasmic bridges that erase phenotypic differences between individual haploid sperm. Here, we show that a large class of mammalian genes are not completely shared across these bridges. We call these genes "genoinformative markers" (GIMs) and show that a subset can act as selfish genetic elements that spread alleles unevenly through murine, bovine, and human populations. We identify evolutionary pressure to avoid conflict between sperm and somatic function as GIMs are enriched for testis-specific gene expression, paralogs, and isoforms. Therefore, GIMs and sperm-level natural selection may help to explain why testis gene expression patterns are an outlier relative to all other tissues.

The Evolutionary and Historical Foundation of the Modern Horse: Lessons from Ancient Genomics

The domestication of the horse some 5,500 years ago followed those of dogs, sheep, goats, cattle, and pigs by ∼2,500-10,000 years. By providing fast transportation and transforming warfare, the horse had an impact on human history with no equivalent in the animal kingdom. Even though the equine sport industry has considerable economic value today, the evolutionary history underlying the emergence of the modern domestic horse remains contentious. In the last decade, novel sequencing technologies have revolutionized our capacity to sequence the complete genome of organisms, including from archaeological remains. Applied to horses, these technologies have provided unprecedented levels of information and have considerably changed models of horse domestication. This review illustrates how ancient DNA, especially ancient genomes, has inspired researchers to rethink the process by which horses were first domesticated and then diversified into a variety of breeds showing a range of traits that are useful to humans.

Inference of natural selection from ancient DNA

Evolutionary processes, including selection, can be indirectly inferred based on patterns of genomic variation among contemporary populations or species. However, this often requires unrealistic assumptions of ancestral demography and selective regimes. Sequencing ancient DNA from temporally spaced samples can inform about past selection processes, as time series data allow direct quantification of population parameters collected before, during, and after genetic changes driven by selection. In this Comment and Opinion, we advocate for the inclusion of temporal sampling and the generation of paleogenomic datasets in evolutionary biology, and highlight some of the recent advances that have yet to be broadly applied by evolutionary biologists. In doing so, we consider the expected signatures of balancing, purifying, and positive selection in time series data, and detail how this can advance our understanding of the chronology and tempo of genomic change driven by selection. However, we also recognize the limitations of such data, which can suffer from postmortem damage, fragmentation, low coverage, and typically low sample size. We therefore highlight the many assumptions and considerations associated with analyzing paleogenomic data and the assumptions associated with analytical methods.

Antagonistic Pleiotropy in Human Disease

Between the 1930s and 1950s, scientists developed key principles of population genetics to try and explain the aging process. Almost a century later, these aging theories, including antagonistic pleiotropy and mutation accumulation, have been experimentally validated in animals. Although the theories have been much harder to test in humans despite research dating back to the 1970s, recent research is closing this evidence gap. Here we examine the strength of evidence for antagonistic pleiotropy in humans, one of the leading evolutionary explanations for the retention of genetic risk variation for non-communicable diseases. We discuss the analytical tools and types of data that are used to test for patterns of antagonistic pleiotropy and provide a primer of evolutionary theory on types of selection as a guide for understanding this mechanism and how it may manifest in other diseases. We find an abundance of non-experimental evidence for antagonistic pleiotropy in many diseases. In some cases, several studies have independently found corroborating evidence for this mechanism in the same or related sets of diseases including cancer and neurodegenerative diseases. Recent studies also suggest antagonistic pleiotropy may be involved in cardiovascular disease and diabetes. There are also compelling examples of disease risk variants that confer fitness benefits ranging from resistance to other diseases or survival in extreme environments. This provides increasingly strong support for the theory that antagonistic pleiotropic variants have enabled improved fitness but have been traded for higher burden of disease later in life. Future research in this field is required to better understand how this mechanism influences contemporary disease and possible consequences for their treatment.

Ancient Genomes Reveal Unexpected Horse Domestication and Management Dynamics

The horse was essential to past human societies but became a recreational animal during the twentieth century as the world became increasingly mechanized. As the author reviews here, recent studies of ancient genomes have revisited the understanding of horse domestication, from the very early stages to the most modern developments. They have uncovered several extinct lineages roaming the far ends of Eurasia some 4000 years ago. They have shown that the domestic horse has been significantly reshaped during the last millennium and experienced a sharp decline in genetic diversity within the last two centuries. At a time when no truly wild horses exist any longer, this calls for enhanced conservation in all endangered populations. These include the Przewalski's horse native to Mongolia, and the many local breeds side-lined by the modern agenda, but yet representing the living heritage of over five millennia of horse breeding.

Indications of strong adaptive population genetic structure in albacore tuna (Thunnus alalunga) in the southwest and central Pacific Ocean

Albacore tuna (Thunnus alalunga) has a distinctly complex life history in which juveniles and adults separate geographically but at times inhabit the same spaces sequentially. The species also migrates long distances and presumably experiences varied regimes of physical stress over a lifetime. There are, therefore, many opportunities for population structure to arise based on stochastic differences or environmental factors that promote local adaptation. However, with the extent of mobility consistently demonstrated by tagged individuals, there is also a strong argument for panmixia within an ocean basin. It is important to confirm such assumptions from a population genetics standpoint for this species in particular because albacore is one of the principal market tuna species that sustains massive global fisheries and yet is also a slow-growing temperate tuna. Consequently, we used 1,837 neutral SNP loci and 89 loci under potential selection to analyze population genetic structure among five sample groups collected from the western and central South Pacific. We found no evidence to challenge panmixia at neutral loci, but strong indications of structuring at adaptive loci. One population sample, from French Polynesia in 2004, was particularly differentiated. Unfortunately, the current study cannot infer whether the divergence is geographic or temporal, or possibly caused by sample distribution. We encourage future studies to include potentially adaptive loci and to continue fine-scale observations within an ocean basin, and not to assume genome-wide panmixia.