Characterising functionally important and ecologically meaningful genetic diversity using a candidate gene approach

Over the past two decades the fields of molecular ecology and population genetics have been dominated by the use of putatively neutral DNA markers, primarily to resolve spatio-temporal patterns of genetic variation to inform our understanding of population structure, gene flow and pedigree. Recent emphasis in comparative functional genomics, however, has fuelled a resurgence of interest in functionally important genetic variation that underpins phenotypic traits of adaptive or ecological significance. It may prove a major challenge to transfer genomics information from classical model species to examine functional diversity in non-model species in natural populations, but already multiple gene-targeted candidate loci with major effect on phenotype and fitness have been identified. Here we briefly describe some of the research strategies used for isolating and characterising functional genetic diversity at candidate gene-targeted loci, and illustrate the efficacy of some of these approaches using our own studies on red grouse (Lagopus lagopus scoticus). We then review how candidate gene markers have been used to: (1) quantify genetic diversity among populations to identify those depauperate in genetic diversity and requiring specific management action; (2) identify the strength and mode of selection operating on individuals within natural populations; and (3) understand direct mechanistic links between allelic variation at single genes and variance in individual fitness.     

solQTL: a tool for QTL analysis,visualization and linking to genomes at SGN database

Background  

A common approach to understanding the genetic basis of complex traits is through identification of associated quantitative trait loci (QTL). Fine mapping QTLs requires several generations of backcrosses and analysis of large populations, which is time-consuming and costly effort. Furthermore, as entire genomes are being sequenced and an increasing amount of genetic and expression data are being generated, a challenge remains: linking phenotypic variation to the underlying genomic variation. To identify candidate genes and understand the molecular basis underlying the phenotypic variation of traits, bioinformatic approaches are needed to exploit information such as genetic map, expression and whole genome sequence data of organisms in biological databases.     

Integrating candidate gene and quantitative genetic approaches to understand variation in timing of breeding in wild tit populations

Two commonly used techniques for estimating the effect of genes on traits in wild populations are the candidate gene approach and quantitative genetic analyses. However, whether these two approaches measure the same underlying processes remains unresolved. Here, we use these two methods to test whether they are alternative or complementary approaches to understanding genetic variation in the timing of reproduction - a key trait involved in adaptation to climate change - in wild tit populations. Our analyses of the candidate gene Clock show weak correlates with timing variables in blue tits, but no association in great tits, confirming earlier results. Quantitative genetic analyses revealed very low levels of both direct (female) and indirect (male) additive genetic variation in timing traits for both species, in contrast to previous studies on these traits, and much lower than generally assumed. Hence, neither method suggests strong genetic effects on the timing of breeding in birds, and further work should seek to assess the generality of these conclusions. We discuss how differences in the genetic control of traits, species life-history and confounding environmental variables may determine how useful integrating these two techniques is to understand the phenotypic variation in wild populations.     

Patterns of variation during adaptation in functionally linked loci

An understanding of the distribution of natural patterns of genetic variation is relevant to such fundamental biological fields as evolution and development. One recent approach to understanding such patterns has been to focus on the constraints that may arise as a function of the network or pathway context in which genes are embedded. Despite theoretical expectations of higher evolutionary constraint for genes encoding upstream versus downstream enzymes in metabolic pathways, empirical results have varied. Here we combine two complementary models from population genetics and enzyme kinetics to explore genetic variation as a function of pathway position when selection acts on whole-pathway flux. We are able to qualitatively reproduce empirically observed patterns of polymorphism and divergence and suggest that expectations should vary depending on the evolutionary trajectory of a population. Upstream genes are initially more polymorphic and diverge faster after an environmental change, while we see the opposite trend as the population approaches its fitness optimum.     

Identification of candidate loci for adaptive phenotypic plasticity in natural populations of spadefoot toads

Phenotypic plasticity allows organisms to alter their phenotype in direct response to changes in the environment. Despite growing recognition of plasticity's role in ecology and evolution, few studies have probed plasticity's molecular bases—especially using natural populations. We investigated the genetic basis of phenotypic plasticity in natural populations of spadefoot toads (Spea multiplicata). Spea tadpoles normally develop into an “omnivore” morph that is favored in long‐lasting, low‐density ponds. However, if tadpoles consume freshwater shrimp or other tadpoles, they can alternatively develop (via plasticity) into a “carnivore” morph that is favored in ephemeral, high‐density ponds. By combining natural variation in pond ecology and morph production with population genetic approaches, we identified candidate loci associated with each morph (carnivores vs. omnivores) and loci associated with adaptive phenotypic plasticity (adaptive vs. maladaptive morph choice). Our candidate morph loci mapped to two genes, whereas our candidate plasticity loci mapped to 14 genes. In both cases, the identified genes tended to have functions related to their putative role in spadefoot tadpole biology. Our results thereby form the basis for future studies into the molecular mechanisms that mediate plasticity in spadefoots. More generally, these results illustrate how diverse loci might mediate adaptive plasticity.     

Essential gene disruptions reveal complex relationships between phenotypic robustness,pleiotropy, and fitness

The concept of robustness in biology has gained much attention recently, but a mechanistic understanding of how genetic networks regulate phenotypic variation has remained elusive. One approach to understand the genetic architecture of variability has been to analyze dispensable gene deletions in model organisms; however, the most important genes cannot be deleted. Here, we have utilized two systems in yeast whereby essential genes have been altered to reduce expression. Using high-throughput microscopy and image analysis, we have characterized a large number of morphological phenotypes, and their associated variation, for the majority of essential genes in yeast. Our results indicate that phenotypic robustness is more highly dependent upon the expression of essential genes than on the presence of dispensable genes. Morphological robustness appears to be a general property of a genotype that is closely related to pleiotropy. While the fitness profile across a range of expression levels is idiosyncratic to each gene, the global pattern indicates that there is a window in which phenotypic variation can be released before fitness effects are observable.     

An essay on the necessity and feasibility of conservation genomics

The basic premise of conservation genetics is that small populations may be genetically threatened. The two steps leading to this premise are: (1) due to prominent influence of random genetic drift and inbreeding allelic and genotypic diversity in small populations is expected to be low, and (2) low allelic diversity and high homozygosity are expected to lead to immediate fitness decreases (inbreeding depression) and a compromised potential for evolutionary adaptation. Conservation genetic research has been strongly stimulated by the application of neutral molecular markers like microsatellites and AFLPs. In general these marker studies have provided evidence for step 1. It is less evident how these markers may provide evidence for step 2. In this essay we argue that, in order to get detailed insight in step 2, adopting a conservation genomic approach, in which conservation genetics will use approaches from ecological and evolutionary functional genomics (ecogenomics), is both necessary and feasible. Conservation genomics is necessary for studying functional genomic variation as function of drift and inbreeding, for studying the mechanisms that relate low genetic variation to low fitness, for integrating environmental and genetic approaches to conservation biology, and for developing modern, fast monitoring tools. The rapid technical and financial developments in genomics currently make conservation genomics feasible, and will improve feasibility in the very near future even further. We therefore argue that conservation genomics personifies part of the near future of conservation genetics.     

Evolutionary adaptation to environmental stressors: a common response at the proteomic level

Mechanistic trade‐offs between traits under selection can shape and constrain evolutionary adaptation to environmental stressors. However, our knowledge of the quantitative and qualitative overlap in the molecular machinery among stress tolerance traits is highly restricted by the challenges of comparing and interpreting data between separate studies and laboratories, as well as to extrapolating between different levels of biological organization. We investigated the expression of the constitutive proteome (833 proteins) of 35 Drosophila melanogaster replicate populations artificially selected for increased resistance to six different environmental stressors. The evolved proteomes were significantly differentiated from replicated control lines. A targeted analysis of the constitutive proteomes revealed a regime‐specific selection response among heat‐shock proteins, which provides evidence that selection also adjusts the constitutive expression of these molecular chaperones. Although the selection response in some proteins was regime specific, the results were dominated by evidence for a “common stress response.” With the exception of high temperature survival, we found no evidence for negative correlations between environmental stress resistance traits, meaning that evolutionary adaptation is not constrained by mechanistic trade‐offs in regulation of functional important proteins. Instead, standing genetic variation and genetic trade‐offs outside regulatory domains likely constrain the evolutionary responses in natural populations.     

Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool‐Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using F_ST‐based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo‐climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo‐climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment.     

UNVEILing connections between genotype,phenotype, and fitness in natural populations

Understanding the links between genetic variation and fitness in natural populations is a central goal of evolutionary genetics. This monumental task spans the fields of classical and molecular genetics, population genetics, biochemistry, physiology, developmental biology, and ecology. Advances to our molecular and developmental toolkits are facilitating integrative approaches across these traditionally separate fields, providing a more complete picture of the genotype‐phenotype map in natural and non‐model systems. Here, we summarize research presented at the first annual symposium of the UNVEIL Network, an NSF‐funded collaboration between the University of Montana and the University of Nebraska, Lincoln, which took place from the 1st to the 3rd of June, 2018. We discuss how this body of work advances basic evolutionary science, what it implies for our ability to predict evolutionary change, and how it might inform novel conservation strategies.     

Next‐generation QTL mapping: crowdsourcing SNPs,without pedigrees

For many molecular ecologists, the mantra and mission of the field of ecological genomics could be encapsulated by the phrase ‘to find the genes that matter’ (Mitchell‐Olds 2001 ; Rockman 2012 ). This phrase of course refers to the early hope and current increasing success in the search for genes whose variation underlies phenotypic variation and fitness in natural populations. In the years since the modern incarnation of the field of ecological genomics, many would agree that the low‐hanging fruit has, at least in principle, been plucked: we now have several elegant examples of genes whose variation influences key adaptive traits in natural populations, and these examples have revealed important insights into the architecture of adaptive variation (Hoekstra et al. 2006 ; Shapiro et al. 2009 ; Chan et al. 2010 ). But how well will these early examples, often involving single genes of large effect on discrete or near‐discrete phenotypes, represent the dynamics of adaptive change for the totality of phenotypes in nature? Will traits exhibiting continuous rather than discrete variation in natural populations have as simple a genetic basis as these early examples suggest (Prasad et al. 2012 ; Rockman 2012 )? Two papers in this issue (Robinson et al. 2013 ; Santure et al. 2013 ) not only suggest answers to these questions but also provide useful extensions of statistical approaches for ecological geneticists to study the genetics of continuous variation in nature. Together these papers, by the same research groups studying evolution in a natural population of Great Tits (Parus major), provide a glimpse of what we should expect as the field begins to dissect the genetic basis of what is arguably the most common type of variation in nature, and how genome‐wide surveys of variation can be applied to natural populations without pedigrees.     

Admixture and the organization of genetic diversity in a butterfly species complex revealed through common and rare genetic variants

Detailed information about the geographic distribution of genetic and genomic variation is necessary to better understand the organization and structure of biological diversity. In particular, spatial isolation within species and hybridization between them can blur species boundaries and create evolutionary relationships that are inconsistent with a strictly bifurcating tree model. Here, we analyse genome‐wide DNA sequence and genetic ancestry variation in Lycaeides butterflies to quantify the effects of admixture and spatial isolation on how biological diversity is organized in this group. We document geographically widespread and pervasive historical admixture, with more restricted recent hybridization. This includes evidence supporting previously known and unknown instances of admixture. The genome composition of admixed individuals varies much more among than within populations, and tree‐ and genetic ancestry‐based analyses indicate that multiple distinct admixed lineages or populations exist. We find that most genetic variants in Lycaeides are rare (minor allele frequency <0.5%). Because the spatial and taxonomic distributions of alleles reflect demographic and selective processes since mutation, rare alleles, which are presumably younger than common alleles, were spatially and taxonomically restricted compared with common variants. Thus, we show patterns of genetic variation in this group are multifaceted, and we argue that this complexity challenges simplistic notions concerning the organization of biological diversity into discrete, easily delineated and hierarchically structured entities.     

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation.     

CONVERGENCE AND DIVERGENCE DURING THE ADAPTATION TO SIMILAR ENVIRONMENTS BY AN AUSTRALIAN GROUNDSEL

Adaptation to replicate environments is often achieved through similar phenotypic solutions. Whether selection also produces convergent genomic changes in these situations remains largely unknown. The variable groundsel, Senecio lautus, is an excellent system to investigate the genetic underpinnings of convergent evolution, because morphologically similar forms of these plants have adapted to the same environments along the coast of Australia. We compared range‐wide patterns of genomic divergence in natural populations of this plant and searched for regions putatively affected by natural selection. Our results indicate that environmental adaptation followed complex genetic trajectories, affecting multiple loci, implying both the parallel recruitment of the same alleles and the divergence of completely different genomic regions across geography. An analysis of the biological functions of candidate genes suggests that adaptation to coastal environments may have occurred through the recruitment of different genes participating in similar processes. The relatively low genetic convergence that characterizes the parallel evolution of S. lautus forms suggests that evolution is more constrained at higher levels of biological organization.     

GENOMIC BASIS OF AGING AND LIFE‐HISTORY EVOLUTION IN DROSOPHILA MELANOGASTER

Natural diversity in aging and other life‐history patterns is a hallmark of organismal variation. Related species, populations, and individuals within populations show genetically based variation in life span and other aspects of age‐related performance. Population differences are especially informative because these differences can be large relative to within‐population variation and because they occur in organisms with otherwise similar genomes. We used experimental evolution to produce populations divergent for life span and late‐age fertility and then used deep genome sequencing to detect sequence variants with nucleotide‐level resolution. Several genes and genome regions showed strong signatures of selection, and the same regions were implicated in independent comparisons, suggesting that the same alleles were selected in replicate lines. Genes related to oogenesis, immunity, and protein degradation were implicated as important modifiers of late‐life performance. Expression profiling and functional annotation narrowed the list of strong candidate genes to 38, most of which are novel candidates for regulating aging. Life span and early age fecundity were negatively correlated among populations; therefore, the alleles we identified also are candidate regulators of a major life‐history trade‐off. More generally, we argue that hitchhiking mapping can be a powerful tool for uncovering the molecular bases of quantitative genetic variation.