Indigenous cattle account for approximately 80% of Uganda's cattle population. These animals are well adapted to the country's ten agroecological zones and are mainly kept under pastoral and agropastoral systems. Unlike commercial breeds, they thrive on low-quality feeds, while tolerating major tropical diseases and parasites including tsetse flies, ticks, and vector-borne infections. Whole-genome sequence (WGS) analysis offers opportunities to uncover genomic regions underlying these adaptations and to trace the genetic footprints of long-term breeding decisions taken by cattle keepers. In this study, WGS data from 95 animals representing six indigenous cattle populations (Ankole, Karamojong, Nganda10, Nganda17, Nkedi, and Ntuku) were analyzed to identify genomic regions under putative selection. Two complementary approaches were applied: enumeration of the µ-statistic in RAiSD and runs of homozygosity (ROH) analysis. RAiSD identified population-level signals, while conserved ROH regions were defined using breed-specific SNP-incidence thresholds. The two methods identified 803 and 49 candidate genes respectively. The top genes identified included SLC37A1 (BTA1), CHCHD3 (BTA4), and RAB3GAP1 (BTA2) detected by RAiSD, and IL26 (BTA5), FBXL7 (BTA20), and HSPA9 (BTA7) contained in ROH. Furthermore, the regions harbored 107 novel genes (92 detected by RAiSD and 15 by ROH), corresponding to 255 quantitative trait loci. The identified genes under putative selection are associated with economically important traits including adaptation to tropical environments, resistance to parasites and diseases, and other farmer-preferred characteristics. These findings provide insights into the genetic basis of adaptation, selection and production in Ugandan indigenous cattle, supporting conservation and breeding strategies to enhance resilience and productivity.

Understanding the genomic basis of diversification is a central goal in evolutionary biology. In recent years, the development and use of pangenomes, a genomic representation of multiple individuals within a lineage (a set of related populations, subspecies, ecotypes, or species), has enabled researchers to differentiate between DNA sequences shared by all individuals of a given lineage (core regions) from those present only in some individuals (accessory or variable regions). Differentiating between core and accessory regions has highlighted a key limitation of relying on a single reference genome: It captures the genetic code of only one individual and this biases genomic analyses and our understanding of diversification. Here, we propose that by identifying genes associated with both core and accessory regions, we can deepen our understanding of the processes underlying diversification. We suggest that analyzing pangenomes and accessory regions will provide deeper insights into diversification, hybridization, and the genetic basis of adaptation and speciation.

ABSTRACTThe use of alternative methods to control cyst nematode populations has accelerated since the ban of chemical nematicides in Europe. The resistant QTL GpaVvrn, derived from the wild species Solanum vernei, is widely present in resistant European potato cultivars and provides strong protection against Globodera pallida populations although a risk of resistance breakdown has already been demonstrated in both experimental evolution studies and field populations. The wild relative S. sparsipilum, harbouring the resistant QTL GpaVspl, would be an interesting alternative source of resistance to control virulent G. pallida. The goal of the present study was to understand the genomics of adaptation of the nematode to these two colinear resistant QTLs. Starting with two natural populations, an experimental evolution approach allowed, after 10 generations on resistant potato genotypes, selecting independent nematode lineages adapted to each QTL. These virulent lineages were analysed through a combination of phenotyping and genome scans approaches. Phenotyping enabled the quantification of virulence levels and confirmed resistance breakdowns. Pool‐Seq whole genome sequencing followed by genome scan analyses identified genomic regions under selection, potentially involved in the adaptive mechanisms to each resistance factor. Candidate genes within these regions provided insights into the genetic basis of adaptation, revealing effectors known to suppress plant immunity. As genome scans highlighted distinct genomic regions for the adaptation to both resistant factors, we were able to predict and phenotypically confirm the absence of cross‐virulence between nematode lineages evolving on GpaVvrn and GpaVspl. These findings have significant implications for the design of effective and sustainable resistance management strategies.

The question of when adaptation involves genetic changes of large effect versus a polygenic response traces back to early debates around Darwin’s “Origin of Species” and remains unanswered today. While there are compelling reasons to expect polygenic adaptation to be common, direct evidence for it is still lacking. In turn, there are hundreds of examples of large effect adaptations across species, but it is unclear whether they are a common occurrence in any given species. Synthesizing the different lines of evidence is further complicated by differences in study designs, limitations and biases. Here, we reframe this long-standing question in terms of the trait under selection and ask how the genetic basis of adaptation is expected to depend on key properties of the genetic variation in the trait (i.e., the trait genetics) and on the changes in selection pressures that act on it (i.e., the “trait ecology”). To study this question, we consider a quantitative trait subject to stabilizing selection and model the response to selection when a population at mutation-selection-drift balance experiences a sudden shift in the optimal value. Using this model, we delimit how the contributions of large effect and polygenic changes to adaptation depend on the genetics and ecology of the trait, as well as other salient factors. This theory allows us to formulate testable predictions about when different modes of adaptation are expected and to outline a framework within which to interpret disparate sources of evidence about the genetic basis of adaptation.

Sugar kelp (Saccharina latissima) is an ecologically and increasingly economically important kelp, distributed from temperate to Arctic rocky shores. However, S. latissima is presently threatened by ongoing climate changes. Genetic variations have previously been identified across S. latissima populations. However, little is known regarding the genetic basis for adaptation and acclimation to different environmental conditions. In this study, a common garden experiment was performed with sporophytes originated from North-Norway (NN), Mid-Norway (MN), and South-Norway (SN), representing areas with highly different temperatures and photoperiods. Transcriptomic analyses revealed significant variation in the gene expression of cultures from North-Norway, associated with low temperature and long photoperiods, compared to Mid- and South-Norway. Differentially expressed genes included genes linked to photosynthesis, chlorophyll biosynthesis, and heat response, suggesting that they are directly involved in temperature and light adaptation. In addition, genes related to growth, metabolism, protein synthesis, and translation were upregulated in the NN genotype, providing evidence that the NN genotype is better adapted to low temperatures than the SN and MN genotypes. Significant variation in gene expression among populations found in this study is influenced by the environment, but genetic differentiation by origin seems to play a role as responses were population specific. This study provides a baseline for deeper insight into the local adaptation potential of S. latissima populations along the Norwegian Coast with implications for the conservation of natural populations.

Adaptation is a critical determinant of the diversification, persistence, and geographic range limits of species. Yet the genetic basis of adaptation is often unknown and potentially underpinned by a wide range of mutational types-from single nucleotide changes to large-scale alterations of chromosome structure. Copy number variation (CNV) is thought to be an important source of adaptive genetic variation, as indicated by decades of candidate gene studies that point to CNVs underlying rapid adaptation to strong selective pressures. Nevertheless, population-genomic studies of CNVs face unique logistical challenges not encountered by other forms of genetic variation. Consequently, few studies have systematically investigated the contributions of CNVs to adaptation at a genome-wide scale. We present a genome-wide analysis of CNV contributing to the adaptation of an invasive weed, Ambrosia artemisiifolia. CNVs show clear signatures of parallel local adaptation between North American (native) and European (invaded) ranges, implying widespread reuse of CNVs during adaptation to shared heterogeneous patterns of selection. We used a local principal component analysis (PCA) to genotype CNV regions in whole-genome sequences of samples collected over the last two centuries. We identified 16 large CNV regions of up to 11.85 megabases in length, eight of which show signals of rapid evolutionary change, with pronounced frequency shifts between historic and modern populations. Our results provide compelling genome-wide evidence that CNV underlies rapid adaptation over contemporary timescales of natural populations.

The dominance of beneficial mutations is a key evolutionary parameter affecting the rate and genetic basis of adaptation, yet it is notoriously difficult to estimate. A leading method to infer it is to compare the relative rates of adaptive substitution for X-linked and autosomal genes, which-according to a classic model by Charlesworth et al. (1987)-is a simple function of the dominance of new beneficial mutations. Recent evidence that rates of adaptive substitution are faster for X-linked genes implies, accordingly, that beneficial mutations are usually recessive. However, this conclusion is incompatible with leading theories of dominance, which predict that beneficial mutations tend to be dominant or overdominant with respect to fitness. To address this incompatibility, we use Fisher's geometric model to predict the distribution of fitness effects of new mutations and the relative rates of positively selected substitution on the X and autosomes. Previous predictions of faster-X theory emerge as a special case of our model in which the phenotypic effects of mutations are small relative to the distance to the phenotypic optimum. But as mutational effects become large relative to the optimum, we observe an elevated tempo of positively selected substitutions on the X relative to the autosomes across a broader range of dominance conditions, including those predicted by theories of dominance. Our results imply that, contrary to previous models, dominant and overdominant beneficial mutations can plausibly generate patterns of faster-X adaptation. We discuss resulting implications for genomic studies of adaptation and inferences of dominance.

Diamondback terrapins (Malaclemys terrapin centrata) exhibit strong environmental adaptability and live in both freshwater and saltwater. However, the genetic basis of this adaptability has not been the focus of research. In this study, we successfully constructed a ∼2.21-Gb chromosome-level genome assembly for M. t. centrata using high-coverage and high-depth genomic sequencing data generated on multiple platforms. The M. t. centrata genome contains 25 chromosomes and the scaffold N50 of ∼143.75 Mb, demonstrating high continuity and accuracy. In total, 53.82% of the genome assembly was composed of repetitive sequences, and 22435 protein-coding genes were predicted. Our phylogenetic analysis indicated that M. t. centrata was closely related to the red-eared slider turtle (Trachemys scripta elegans), with divergence approximately ∼23.6 million years ago (Mya) during the early Neogene period of the Cenozoic era. The population size of M. t. centrata decreased significantly over the past ∼14 Mya during the Cenozoic era. Comparative genomic analysis indicated that 36 gene families related to ion transport were expanded and several genes (AQP3, solute carrier subfamily, and potassium channel genes) underwent specific amino acid site mutations in the M. t. centrata genome. Changes to these ion transport-related genes may have contributed to the remarkable salinity adaptability of diamondback terrapin. The results of this study not only provide a high-quality reference genome for M. t. centrata but also elucidate the possible genetic basis for salinity adaptation in this species.

Papua New Guinea (PNG) hosts distinct environments mainly represented by the ecoregions of the Highlands and Lowlands that display increased altitude and a predominance of pathogens, respectively. Since its initial peopling approximately 50,000 years ago, inhabitants of these ecoregions might have differentially adapted to the environmental pressures exerted by each of them. However, the genetic basis of adaptation in populations from these areas remains understudied. Here, we investigated signals of positive selection in 62 highlanders and 43 lowlanders across 14 locations in the main island of PNG using whole-genome genotype data from the Oceanian Genome Variation Project (OGVP) and searched for signals of positive selection through population differentiation and haplotype-based selection scans. Additionally, we performed archaic ancestry estimation to detect selection signals in highlanders within introgressed regions of the genome. Among highland populations we identified candidate genes representing known biomarkers for mountain sickness (SAA4, SAA1, PRDX1, LDHA) as well as candidate genes of the Notch signaling pathway (PSEN1, NUMB, RBPJ, MAML3), a novel proposed pathway for high altitude adaptation in multiple organisms. We also identified candidate genes involved in oxidative stress, inflammation, and angiogenesis, processes inducible by hypoxia, as well as in components of the eye lens and the immune response. In contrast, candidate genes in the lowlands are mainly related to the immune response (HLA-DQB1, HLA-DQA2, TAAR6, TAAR9, TAAR8, RNASE4, RNASE6, ANG). Moreover, we find two candidate regions to be also enriched with archaic introgressed segments, suggesting that archaic admixture has played a role in the local adaptation of PNG populations.

The genus Pristimantis diversified in the tropical Andes mountains and is the most speciose genus of terrestrial vertebrates. Pristimantis are notable among frogs in that they thrive at high elevations (>2,000 m) and are direct developers without a tadpole stage. Despite their ecological significance, little is known about the genetic and physiological traits enabling their success. We conducted transcriptomic analysis on seven Pristimantis species sampled across elevations in the Ecuadorean Andes to explore three hypotheses for their success: (i) unique genes are under selection relative to all other frogs, (ii) common selection occurs across all direct developers, or (iii) common selection occurs across all high-elevation frog clades. Comparative analysis with 34 frog species revealed unique positive selection in Pristimantis genes related to aerobic respiration, hemostasis, signaling, cellular transportation of proteins and ions, and immunity. Additionally, we detected positive selection across all direct developers for genes associated with oxygenase activity and metal ion binding. While many genes under selection in Pristimantis were not positively selected in other high-elevation frog species, we identified some shared genes and pathways linked to lipid metabolism, innate immunity, and cellular redox processes. We observed more positive selection in duplicated- versus single-copy genes, while relaxed purifying selection was prevalent in single-copy genes. Notably, copy number of an innate immunity complement gene was positively correlated with Pristimantis species elevation. Our findings contribute novel insights into the genetic basis of adaptation in Pristimantis and provide a foundation for future studies on the evolutionary mechanisms leading to direct development and coping with high elevations.

Understanding how evolutionary factors related to climate adaptation and human selection have influenced the genetic architecture of domesticated animals is of great interest in biology. In the current study, by using 304 whole genomes from different geographical regions (including Europe, north Africa, Southwest Asia, east Asia, west Africa, south Asia, east Africa, Australia and Turkey), We evaluate global sheep population dynamics in terms of genetic variation and population structure. We further conducted comparative population analysis to study the genetic underpinnings of climate adaption to local environments and also morphological traits. In order to identify genomic signals under selection, we applied fixation index (FST) and also nucleotide diversity (θπ) statistical measurements. Our results revealed several candidate genes on different chromosomes under selection for local climate adaptation (e.g. HOXC12, HOXC13, IRF1, FGD2 and GNAQ), body size (PDGFA, HMGA2, PDE3A) and also morphological related traits (RXFP2). The discovered candidate genes may offer newel insights into genetic underpinning of regional adaptation and commercially significant features in local sheep.

Most traits are polygenic, and the contributing loci can be identified by genome-wide association studies. The genetic basis of adaptation (adaptive architecture) is, however, difficult to characterize. Here, we propose to study the adaptive architecture of traits by monitoring the evolution of their phenotypic variance during adaptation to a new environment in well-defined laboratory conditions. Extensive computer simulations show that the evolution of phenotypic variance in a replicated experimental evolution setting can distinguish between oligogenic and polygenic adaptive architectures. We compared gene expression variance in male Drosophila simulans before and after 100 generations of adaptation to a novel hot environment. The variance change in gene expression was indistinguishable for genes with and without a significant change in mean expression after 100 generations of evolution. We suggest that the majority of adaptive gene expression evolution can be explained by a polygenic architecture. We propose that tracking the evolution of phenotypic variance across generations can provide an approach to characterize the adaptive architecture.

Abstract Insects are essential not only for ecosystem functioning and food security but also comprise some of the world's most destructive invasive species. Therefore, both insect declines and invasions raise major conservation concerns globally and call for respective conservation or mitigation measures. However, studies of insects are hampered by intrinsic biological features of these organisms, which include extreme population fluctuations, a huge diversity of ecological strategies, and common cryptic species. Population genetics provides a large toolkit to adequately accommodate those features, thereby enabling researchers to inform and monitor the efficacy of conservation and mitigation programs. Here, we provide an overview of the molecular and analytical methods that are relevant to insect conservation or mitigation and highlight the challenges involved in their implementation. We detail how and why temporal changes in genetic diversity, population structure and migration, and the genetic basis of adaptation should be taken into account to inform insect management programs. Finally, we review the barriers to the broad adoption of population genetics in insect research and provide guidelines to facilitate the use of these methods by stakeholders. Overall, this review provides theoretical and practical guidelines for implementing population genetics in both insect conservation and control.

Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST ), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genetic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families – genes directly mediating interactions with plant chemical defenses – underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many lineages are ancient (>150 million years ago [mya]), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several non-herbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza have among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on plants (bitter or electrophilic phytotoxins) or their ancestral diet (yeast and fruit volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight strong gene candidates that have also been linked to other dietary transitions in Drosophila.

The genetic basis of adaptation is driven by both selection and the spectrum of available mutations. Given that the rate of mutation is not uniformly distributed across the genome and varies depending on the environment, understanding the signatures of selection across the genome is aided by first establishing what the expectations of genetic change are from mutation. To determine the interaction between salt stress, selection, and mutation across the genome, we compared mutations observed in a selection experiment for salt tolerance in Chlamydomonas reinhardtii to those observed in mutation accumulation (MA) experiments with and without salt exposure. MA lines evolved under salt stress had a single-nucleotide mutation rate of , similar to that of MA lines under standard conditions ( ). However, we found that salt stress led to an increased rate of indel mutations, but that many of these mutations were removed under selection. Finally, lines adapted to salt also showed excess clustering of mutations in the genome and the co-expression network, suggesting a role for positive selection in retaining mutations in particular compartments of the genome during the evolution of salt tolerance. Our study shows that characterizing mutation rates and spectra expected under stress helps disentangle the effects of environment and selection duringadaptation.

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

Does evolution proceed in small steps or large leaps? How repeatable is evolution? How constrained is the evolutionary process? Answering these long-standing questions in evolutionary biology is indispensable for both understanding how extant biodiversity has evolved and predicting how organisms and ecosystems will respond to changing environments in the future. Understanding the genetic basis of phenotypic diversification and speciation in natural populations is key to properly answering these questions. The leap forward in genome sequencing technologies has made it increasingly easier to not only investigate the genetic architecture but also identify the variant sites underlying adaptation and speciation in natural populations. Furthermore, recent advances in genome editing technologies are making it possible to investigate the functions of each candidate gene in organisms from natural populations. In this article, we discuss how these recent technological advances enable the analysis of causative genes and mutations and how such analysis can help answer long-standing evolutionary biology questions.This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

At the birth of supergenes, the genomic landscape is dramatically re-organized leading to pronounced differences in phenotypes and increased intrasexual diversity. Two of the best-studied supergenes in vertebrates are arguably the inversion polymorphisms on chromosomes 2 and 11 in the white-throated sparrow (Zonotrichia albicollis) and the ruff (Calidris pugnax), respectively. In both species, regions of suppressed recombination determine plumage coloration and social behavioural phenotypes. Despite the apparent lack of gene overlap between these two supergenes, in both cases the alternative phenotypes seem to be driven largely by alterations in steroid hormone pathways. Here, we explore the interplay between genomic architecture and steroid-related genes. Due to the highly pleiotropic effects of steroid-related genes and their universal involvement in social behaviour and transcriptomic regulation, processes favouring their linkage are likely to have substantial effects on the evolution of behavioural phenotypes, individual fitness, and life-history strategies. We propose that inversion-related differentiation and regulatory changes in steroid-related genes lie at the core of phenotypic differentiation in both of these interesting species. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Identifying the genetic basis of adaptation is a central goal of evolutionary biology. However, identifying genes and mutations affecting fitness remains challenging because a large number of traits and variants can influence fitness. Selected phenotypes can also be difficult to know a priori, complicating top-down genetic approaches for trait mapping that involve crosses or genome-wide association studies. In such cases, experimental genetic approaches, where one maps fitness directly and attempts to infer the traits involved afterwards, can be valuable. Here, we re-analyse data from a transplant experiment involving Timema stick insects, where five physically clustered single-nucleotide polymorphisms associated with cryptic body coloration were shown to interact to affect survival. Our analysis covers a larger genomic region than past work and revealed a locus previously not identified as associated with survival. This locus resides near a gene, Punch (Pu), involved in pteridine pigments production, implying that it could be associated with an unmeasured coloration trait. However, by combining previous and newly obtained phenotypic data, we show that this trait is not eye or body coloration. We discuss the implications of our results for the discovery of traits, genes and mutations associated with fitness in other systems, as well as for supergene evolution. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.