Genomic-assisted breeding for climate-smart coffee

Coffee is a universal beverage that drives a multi-industry market on a global basis. Today, the sustainability of coffee production is threatened by accelerated climate changes. In this work, we propose the implementation of genomic-assisted breeding for climate-smart coffee in Coffea canephora. This species is adapted to higher temperatures and is more resilient to biotic and abiotic stresses. After evaluating two populations, over multiple harvests, and under severe drought weather condition, we dissected the genetic architecture of yield, disease resistance, and quality-related traits. By integrating genome-wide association studies and diallel analyses, our contribution is four-fold: (i) we identified a set of molecular markers with major effects associated with disease resistance and post-harvest traits, while yield and plant architecture presented a polygenic background; (ii) we demonstrated the relevance of nonadditive gene actions and projected hybrid vigor when genotypes from different geographically botanical groups are crossed; (iii) we computed medium-to-large heritability values for most of the traits, representing potential for fast genetic progress; and (iv) we provided a first step toward implementing molecular breeding to accelerate improvements in C. canephora. Altogether, this work is a blueprint for how quantitative genetics and genomics can assist coffee breeding and support the supply chain in the face of the current global changes.

Understanding the joint evolution of dispersal and host specialisation using phytophagous arthropods as a model group

Theory generally predicts that host specialisation and dispersal should evolve jointly. Indeed, many models predict that specialists should be poor dispersers to avoid landing on unsuitable hosts while generalists will have high dispersal abilities. Phytophagous arthropods are an excellent group to test this prediction, given extensive variation in their host range and dispersal abilities. Here, we explore the degree to which the empirical literature on this group is in accordance with theoretical predictions. We first briefly outline the theoretical reasons to expect such a correlation. We then report empirical studies that measured both dispersal and the degree of specialisation in phytophagous arthropods. We find a correlation between dispersal and levels of specialisation in some studies, but with wide variation in this result. We then review theoretical attributes of species and environment that may blur this correlation, namely environmental grain, temporal heterogeneity, habitat selection, genetic architecture, and coevolution between plants and herbivores. We argue that theoretical models fail to account for important aspects, such as phenotypic plasticity and the impact of selective forces stemming from other biotic interactions, on both dispersal and specialisation. Next, we review empirical caveats in the study of this interplay. We find that studies use different measures of both dispersal and specialisation, hampering comparisons. Moreover, several studies do not provide independent measures of these two traits. Finally, variation in these traits may occur at scales that are not being considered. We conclude that this correlation is likely not to be expected from large-scale comparative analyses as it is highly context dependent and should not be considered in isolation from the factors that modulate it, such as environmental scale and heterogeneity, intrinsic traits or biotic interactions. A stronger crosstalk between theoretical and empirical studies is needed to understand better the prevalence and basis of the correlation between dispersal and specialisation.

DNA methylation differences between stick insect ecotypes

Epigenetic mechanisms, such as DNA methylation, can influence gene regulation and affect phenotypic variation, raising the possibility that they contribute to ecological adaptation. Beginning to address this issue requires high-resolution sequencing studies of natural populations to pinpoint epigenetic regions of potential ecological and evolutionary significance. However, such studies are still relatively uncommon, especially in insects, and are mainly restricted to a few model organisms. Here, we characterize patterns of DNA methylation for natural populations of Timema cristinae adapted to two host plant species (i.e. ecotypes). By integrating results from sequencing of whole transcriptomes, genomes and methylomes, we investigate whether environmental, host and genetic differences of these stick insects are associated with methylation levels of cytosine nucleotides in the CpG context. We report an overall genome-wide methylation level for T. cristinae of ~14%, with methylation being enriched in gene bodies and impoverished in repetitive elements. Genome-wide DNA methylation variation was strongly positively correlated with genetic distance (relatedness), but also exhibited significant host-plant effects. Using methylome-environment association analysis, we pinpointed specific genomic regions that are differentially methylated between ecotypes, with these regions being enriched for genes with functions in membrane processes. The observed association between methylation variation and genetic relatedness, and with the ecologically important variable of host plant, suggests a potential role for epigenetic modification in T. cristinae adaptation. To substantiate such adaptive significance, future studies could test whether methylation can be transmitted across generations and the extent to which it responds to experimental manipulation in field and laboratory studies.

Disruptive selection and the evolution of discrete color morphs in Timema stick insects

A major unresolved issue in biology is why phenotypic and genetic variation is sometimes continuous, yet other times packaged into discrete units of diversity, such as morphs, ecotypes, and species. In theory, ecological discontinuities can impose strong disruptive selection that promotes the evolution of discrete forms, but direct tests of this hypothesis are lacking. Here, we show that Timema stick insects exhibit genetically determined color morphs that range from weakly to strongly discontinuous. Color data from nature and a manipulative field experiment demonstrate that greater morph differentiation is associated with shifts from host plants exhibiting more continuous color variation to those exhibiting greater coloration distance between green leaves and brown stems, the latter of which generates strong disruptive selection. Our results show how ecological factors can promote discrete variation, and we further present results on how this can have variable effects on the genetic differentiation that promotes speciation.

Genomic prediction in the wild: A case study in Soay sheep

Genomic prediction, the technique whereby an individual's genetic component of their phenotype is estimated from its genome, has revolutionised animal and plant breeding and medical genetics. However, despite being first introduced nearly two decades ago, it has hardly been adopted by the evolutionary genetics community studying wild organisms. Here, genomic prediction is performed on eight traits in a wild population of Soay sheep. The population has been the focus of a >30 year evolutionary ecology study and there is already considerable understanding of the genetic architecture of the focal Mendelian and quantitative traits. We show that the accuracy of genomic prediction is high for all traits, but especially those with loci of large effect segregating. Five different methods are compared, and the two methods that can accommodate zero-effect and large-effect loci in the same model tend to perform best. If the accuracy of genomic prediction is similar in other wild populations, then there is a real opportunity for pedigree-free molecular quantitative genetics research to be enabled in many more wild populations; currently the literature is dominated by studies that have required decades of field data collection to generate sufficiently deep pedigrees. Finally, some of the potential applications of genomic prediction in wild populations are discussed.

A polygenic architecture with habitat-dependent effects underlies ecological differentiation in Silene

Ecological differentiation can drive speciation but it is unclear how the genetic architecture of habitat-dependent fitness contributes to lineage divergence. We investigated the genetic architecture of cumulative flowering, a fitness component, in second-generation hybrids between Silene dioica and Silene latifolia transplanted into the natural habitat of each species. We used reduced-representation sequencing and Bayesian sparse linear mixed models (BSLMMs) to analyze the genetic control of cumulative flowering in each habitat. Our results point to a polygenic architecture of cumulative flowering. Allelic effects were mostly beneficial or deleterious in one habitat and neutral in the other. Positive-effect alleles often were derived from the native species, whereas negative-effect alleles, at other loci, tended to originate from the non-native species. We conclude that ecological differentiation is governed and maintained by many loci with small, habitat-dependent effects consistent with conditional neutrality. This pattern may result from differences in selection targets in the two habitats and from environmentally dependent deleterious load. Our results further suggest that selection for native alleles and against non-native alleles acts as a barrier to gene flow between species.

Genomic Traces of the Fruit Fly Anastrepha obliqua Associated with Its Polyphagous Nature

Anastrepha obliqua (Macquart) (Diptera: Tephritidae) is an important pest in the neotropical region. It is considered a polyphagous insect, meaning it infests plants of different taxonomic families and readily colonizes new host plants. The change to new hosts can lead to diversification and the formation of host races. Previous studies investigating the effect of host plants on population structure and selection in Anastrepha obliqua have focused on the use of data from the mitochondrial DNA sequence and microsatellite markers of nuclear DNA, and there are no analyses at the genomic level. To better understand this issue, we used a pooled restriction site-associated DNA sequencing (pooled RAD-seq) approach to assess genomic differentiation and population structure across sympatric populations of Anastrepha obliqua that infest three host plants-Spondias purpurea (red mombin), Mangifera indica (mango) of the family Anacardiaceae and Averrhoa carambola (carambola) of the family Oxalidaceae-in sympatric populations of the species Anastrepha obliqua of Inter-Andean Valley of the Cauca River in southwestern Colombia. Our results show genomic differentiation of populations from carambola compared to mango and red mombin populations, but the genetic structure was mainly established by geography rather than by the host plant. On the other hand, we identified 54 SNPs in 23 sequences significantly associated with the use of the host plant. Of these 23 sequences, we identified 17 candidate genes and nine protein families, of which four protein families are involved in the nutrition of these flies. Future studies should investigate the adaptive processes undergone by phytophagous insects in the Neotropics, using fruit flies as a model and state-of-the-art molecular tools.

The genomics of phenotypically differentiated Asellus aquaticus cave, surface stream and lake ecotypes

Organisms well suited for the study of ecotype formation have wide distribution ranges, where they adapt to multiple drastically different habitats repeatedly over space and time. Here we study such ecotypes in a Crustacean model, Asellus aquaticus, a commonly occurring isopod found in freshwater habitats as diverse as streams, caves and lakes. Previous studies focusing on cave vs. surface ecotypes have attributed depigmentation, eye loss and prolonged antennae to several south European cave systems. Likewise, surveys across multiple Swedish lakes have identified the presence of dark-pigmented "reed" and light-pigmented "stonewort" ecotypes, which can be found within the same lake. In this study, we sequenced the first draft genome of A. aquaticus, and subsequently use this to map reads and call variants in surface stream, cave and two lake ecotypes. In addition, the draft genome was combined with a RADseq approach to perform a quantitative trait locus (QTL) mapping study using a laboratory bred F₂ and F₄ cave × surface intercross. We identified genomic regions associated with body pigmentation, antennae length and body size. Furthermore, we compared genome-wide differentiation between natural populations and found several genes potentially associated with these habitats. The assessment of the cave QTL regions in the light-dark comparison of lake populations suggests that the regions associated with cave adaptation are also involved with genomic differentiation in the lake ecotypes. These demonstrate how troglomorphic adaptations can be used as a model for related ecotype formation.

Standing genetic variation in laboratory populations of insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae) for the evolution of resistance

Insecticides can exert strong selection on insect pest species, including those that vector diseases, and have led to rapid evolution of resistance. Despite such rapid evolution, relatively little is known about standing genetic variation for resistance in insecticide-susceptible populations of many species. To help fill this knowledge gap, we generated genotyping-by-sequencing data from insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis sand flies that survived or died from a sub-diagnostic exposure to either permethrin or malathion using a modified version of the Centers for Disease Control and Prevention bottle bioassay. Multi-locus genome-wide association mapping methods were used to quantify standing genetic variation for insecticide resistance in these populations and to identify specific alleles associated with insecticide survival. For each insecticide treatment, we estimated the proportion of the variation in survival explained by the genetic data (i.e., "chip" heritability) and the number and contribution of individual loci with measurable effects. For all treatments, survival to an insecticide exposure was heritable with a polygenic architecture. Both P. papatasi and L. longipalpis had alleles for survival that resided within many genes throughout their genomes. The implications for resistance conferred by many alleles, as well as inferences made about the utility of laboratory insecticide resistance association studies compared to field observations, are discussed.

Clair JB, Crepeau MW, Salzberg SL, Langley CH, Allen B, Neale DB. Dissecting the Polygenic Basis of Cold Adaptation Using Genome-Wide Association of Traits and Environmental Data in Douglas-fir

Understanding the genomic and environmental basis of cold adaptation is key to understand how plants survive and adapt to different environmental conditions across their natural range. Univariate and multivariate genome-wide association (GWAS) and genotype-environment association (GEA) analyses were used to test associations among genome-wide SNPs obtained from whole-genome resequencing, measures of growth, phenology, emergence, cold hardiness, and range-wide environmental variation in coastal Douglas-fir (Pseudotsuga menziesii). Results suggest a complex genomic architecture of cold adaptation, in which traits are either highly polygenic or controlled by both large and small effect genes. Newly discovered associations for cold adaptation in Douglas-fir included 130 genes involved in many important biological functions such as primary and secondary metabolism, growth and reproductive development, transcription regulation, stress and signaling, and DNA processes. These genes were related to growth, phenology and cold hardiness and strongly depend on variation in environmental variables such degree days below 0c, precipitation, elevation and distance from the coast. This study is a step forward in our understanding of the complex interconnection between environment and genomics and their role in cold-associated trait variation in boreal tree species, providing a baseline for the species' predictions under climate change.

What We Don't Know About Diet-Breadth Evolution in Herbivorous Insects

Half a million species of herbivorous insects have been described. Most of them are diet specialists, using only a few plant species as hosts. Biologists suspect that their specificity is key to their diversity. But why do herbivorous insects tend to be diet specialists? In this review, we catalog a broad range of explanations. We review the evidence for each and suggest lines of research to obtain the evidence we lack. We then draw attention to a second major question, namely how changes in diet breadth affect the rest of a species’ biology. In particular, we know little about how changes in diet breadth feed back on genetic architecture, the population genetic environment, and other aspects of a species’ ecology. Knowing more about how generalists and specialists differ should go a long way toward sorting out potential explanations of specificity, and yield a deeper understanding of herbivorous insect diversity.

Large-scale mutation in the evolution of a gene complex for cryptic coloration

The types of mutations affecting adaptation in the wild are only beginning to be understood. In particular, whether structural changes shape adaptation by suppressing recombination or by creating new mutations is unresolved. Here, we show that multiple linked but recombining loci underlie cryptic color morphs of Timema chumash stick insects. In a related species, these loci are found in a region of suppressed recombination, forming a supergene. However, in seven species of Timema, we found that a megabase-size "supermutation" has deleted color loci in green morphs. Moreover, we found that balancing selection likely contributes more to maintaining this mutation than does introgression. Our results show how suppressed recombination and large-scale mutation can help to package gene complexes into discrete units of diversity such as morphs, ecotypes, or species.

The target of selection matters: An established resistance-development-time negative genetic trade-off is not found when selecting on development time

Trade-offs are fundamental to evolutionary outcomes and play a central role in eco-evolutionary theory. They are often examined by experimentally selecting on one life-history trait and looking for negative correlations in other traits. For example, populations of the moth Plodia interpunctella selected to resist viral infection show a life-history cost with longer development times. However, we rarely examine whether the detection of such negative genetic correlations depends on the trait on which we select. Here, we examine a well-characterized negative genotypic trade-off between development time and resistance to viral infection in the moth Plodia interpunctella and test whether selection on a phenotype known to be a cost of resistance (longer development time) leads to the predicted correlated increase in resistance. If there is tight pleiotropic relationship between genes that determine development time and resistance underpinning this trade-off, we might expect increased resistance when we select on longer development time. However, we show that selecting for longer development time in this system selects for reduced resistance when compared to selection for shorter development time. This shows how phenotypes typically characterized by a trade-off can deviate from that trade-off relationship, and suggests little genetic linkage between the genes governing viral resistance and those that determine response to selection on the key life-history trait. Our results are important for both selection strategies in applied biological systems and for evolutionary modelling of host-parasite interactions.

The green-brown polymorphism of the club-legged grasshopper Gomphocerus sibiricus is heritable and appears genetically simple

Background

Local coexistence of distinct, genetically determined color morphs can be unstable and transitional. Stable, long-term coexistence requires some form of balancing selection to protect morphs from getting lost by directional selection or genetic drift. However, not all phenotypic polymorphism need to have a genetic basis. We here report on the genetic basis of two color polymorphisms in the club-legged grasshopper Gomphocerus sibiricus: a green-brown polymorphism that is phylogenetically and geographically widespread among orthopteran insects and a pied-brown pattern polymorphism that is shared among many gomphocerine grasshoppers.

Results

We found a remarkably clear outcome of matings within and between morph that suggest not only that the green-brown polymorphism is heritable in this species, but that results can be most parsimoniously explained by a single autosomal locus with two alleles in which the green allele is dominant over the brown allele. A few individuals did not match this pattern and suggest the existence of genetic modifiers and/or developmental phenocopies. We also show that the pied-brown polymorphism is highly heritable, although the evidence for the involvement of one or more loci is less clear-cut.

Conclusions

Overall, our data demonstrate that the two polymorphisms are heritable in the club-legged grasshopper and appear genetically simple, at least with respect to green morphs. The results are consistent with the idea that the synthesis or transport of a pigment involved in the production of green coloration (likely biliverdin) is lost by homozygosity for loss-of-function alleles in brown individuals. The apparently simple genetic architecture of the green-brown polymorphism offer potential for studying balancing selection in the field and for genetic mapping in this species.

Host plant-related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi

Elucidating the mechanisms and conditions facilitating the formation of biodiversity are central topics in evolutionary biology. A growing number of studies imply that divergent ecological selection may often play a critical role in speciation by counteracting the homogenising effects of gene flow. Several examples involve phytophagous insects, where divergent selection pressures associated with host plant shifts may generate reproductive isolation, promoting speciation. Here, we use ddRADseq to assess the population structure and to test for host‐related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi (L., 1758) (Diptera: Tephritidae). This tephritid is distributed throughout Europe and western Asia, and has adapted to two different genera of host plants, Prunus spp. (cherries) and Lonicera spp. (honeysuckle). Our data imply that geographic distance and geomorphic barriers serve as the primary factors shaping genetic population structure across the species range. Locally, however, flies genetically cluster according to host plant, with consistent allele frequency differences displayed by a subset of loci between Prunus and Lonicera flies across four sites surveyed in Germany and Norway. These 17 loci display significantly higher F_ST values between host plants than others. They also showed high levels of linkage disequilibrium within and between Prunus and Lonicera flies, supporting host‐related selection and reduced gene flow. Our findings support the existence of sympatric host races in R. cerasi embedded within broader patterns of geographic variation in the fly, similar to the related apple maggot, Rhagoletis pomonella, in North America.

Admixture mapping in interspecific Populus hybrids identifies classes of genomic architectures for phytochemical, morphological and growth traits

The genomic architecture of functionally important traits is key to understanding the maintenance of reproductive barriers and trait differences when divergent populations or species hybridize. We conducted a genome-wide association study (GWAS) to study trait architecture in natural hybrids of two ecologically divergent Populus species. We genotyped 472 seedlings from a natural hybrid zone of Populus alba and Populus tremula for genome-wide markers from reduced representation sequencing, phenotyped the plants in common gardens for 46 phytochemical (phenylpropanoid), morphological and growth traits, and used a Bayesian polygenic model for mapping. We detected three classes of genomic architectures: traits with finite, detectable associations of genetic loci with phenotypic variation in addition to highly polygenic heritability; traits with indications for polygenic heritability only; and traits with no detectable heritability. For the first class, we identified genome regions with plausible candidate genes for phenylpropanoid biosynthesis or its regulation, including MYB transcription factors and glycosyl transferases. GWAS in natural, recombinant hybrids represent a promising step towards resolving the genomic architecture of phenotypic traits in long-lived species. This facilitates the fine-mapping and subsequent functional characterization of genes and networks causing differences in hybrid performance and fitness.

The Evolutionary History of Nebraska Deer Mice: Local Adaptation in the Face of Strong Gene Flow

The interplay of gene flow, genetic drift, and local selective pressure is a dynamic process that has been well studied from a theoretical perspective over the last century. Wright and Haldane laid the foundation for expectations under an island-continent model, demonstrating that an island-specific beneficial allele may be maintained locally if the selection coefficient is larger than the rate of migration of the ancestral allele from the continent. Subsequent extensions of this model have provided considerably more insight. Yet, connecting theoretical results with empirical data has proven challenging, owing to a lack of information on the relationship between genotype, phenotype, and fitness. Here, we examine the demographic and selective history of deer mice in and around the Nebraska Sand Hills, a system in which variation at the Agouti locus affects cryptic coloration that in turn affects the survival of mice in their local habitat. We first genotyped 250 individuals from 11 sites along a transect spanning the Sand Hills at 660,000 single nucleotide polymorphisms across the genome. Using these genomic data, we found that deer mice first colonized the Sand Hills following the last glacial period. Subsequent high rates of gene flow have served to homogenize the majority of the genome between populations on and off the Sand Hills, with the exception of the Agouti pigmentation locus. Furthermore, mutations at this locus are strongly associated with the pigment traits that are strongly correlated with local soil coloration and thus responsible for cryptic coloration.

SNP-skimming: A fast approach to map loci generating quantitative variation in natural populations

Genome-wide association mapping (GWAS) is a method to estimate the contribution of segregating genetic loci to trait variation. A major challenge for applying GWAS to nonmodel species has been generating dense genome-wide markers that satisfy the key requirement that marker data are error-free. Here, we present an approach to map loci within natural populations using inexpensive shallow genome sequencing. This "SNP-skimming" approach involves two steps: an initial genome-wide scan to identify putative targets followed by deep sequencing for confirmation of targeted loci. We apply our method to a test data set of floral dimension variation in the plant Penstemon virgatus, a member of a genus that has experienced dynamic floral adaptation that reflects repeated transitions in primary pollinator. The ability to detect SNPs that generate phenotypic variation depends on population genetic factors such as population allele frequency, effect size and epistasis, as well as sampling effects contingent on missing data and genotype uncertainty. However, both simulations and the Penstemon data suggest that the most significant tests from the initial SNP skim are likely to be true positives-loci with subtle but significant quantitative effects on phenotype. We discuss the promise and limitations of this method and consider optimal experimental design for a given sequencing effort. Simulations demonstrate that sampling a larger number of individual at the expense of average read depth per individual maximizes the power to detect loci.

Water availability drives signatures of local adaptation in whitebark pine (Pinus albicaulis Engelm.) across fine spatial scales of the Lake Tahoe Basin, USA

Patterns of local adaptation at fine spatial scales are central to understanding how evolution proceeds, and are essential to the effective management of economically and ecologically important forest tree species. Here, we employ single and multilocus analyses of genetic data (n = 116 231 SNPs) to describe signatures of fine-scale adaptation within eight whitebark pine (Pinus albicaulis Engelm.) populations across the local extent of the environmentally heterogeneous Lake Tahoe Basin, USA. We show that despite highly shared genetic variation (F_ST  = 0.0069), there is strong evidence for adaptation to the rain shadow experienced across the eastern Sierra Nevada. Specifically, we build upon evidence from a common garden study and find that allele frequencies of loci associated with four phenotypes (mean = 236 SNPs), 18 environmental variables (mean = 99 SNPs), and those detected through genetic differentiation (n = 110 SNPs) exhibit significantly higher signals of selection (covariance of allele frequencies) than could be expected to arise, given the data. We also provide evidence that this covariance tracks environmental measures related to soil water availability through subtle allele frequency shifts across populations. Our results replicate empirical support for theoretical expectations of local adaptation for populations exhibiting strong gene flow and high selective pressures and suggest that ongoing adaptation of many P. albicaulis populations within the Lake Tahoe Basin will not be constrained by the lack of genetic variation. Even so, some populations exhibit low levels of heritability for the traits presumed to be related to fitness. These instances could be used to prioritize management to maintain adaptive potential. Overall, we suggest that established practices regarding whitebark pine conservation be maintained, with the additional context of fine-scale adaptation.

Novel Candidate Genes Underlying Extreme Trophic Specialization in Caribbean Pupfishes

The genetic changes responsible for evolutionary transitions from generalist to specialist phenotypes are poorly understood. Here we examine the genetic basis of craniofacial traits enabling novel trophic specialization in a sympatric radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas. This recent radiation consists of a generalist species and two novel specialists: a small-jawed "snail-eater" and a large-jawed "scale-eater." We genotyped 12 million single nucleotide polymorphisms (SNPs) by whole-genome resequencing of 37 individuals of all three species from nine populations and integrated genome-wide divergence scans with association mapping to identify divergent regions containing putatively causal SNPs affecting jaw size-the most rapidly diversifying trait in this radiation. A mere 22 fixed variants accompanied extreme ecological divergence between generalist and scale-eater species. We identified 31 regions (20 kb) containing variants fixed between specialists that were significantly associated with variation in jaw size which contained 11 genes annotated for skeletal system effects and 18 novel candidate genes never previously associated with craniofacial phenotypes. Six of these 31 regions showed robust signs of hard selective sweeps after accounting for demographic history. Our data are consistent with predictions based on quantitative genetic models of adaptation, suggesting that the effect sizes of regions influencing jaw phenotypes are positively correlated with distance between fitness peaks on a complex adaptive landscape.

The genomic basis of eco-evolutionary dynamics

Recent recognition that ecological and evolutionary processes can operate on similar timescales has led to a rapid increase in theoretical and empirical studies on eco-evolutionary dynamics. Progress in the fields of evolutionary biology, genomics and ecology is greatly enhancing our understanding of rapid adaptive processes, the predictability of adaptation and the genetics of ecologically important traits. However, progress in these fields has proceeded largely independently of one another. In an attempt to better integrate these fields, the centre for 'Adaptation to a Changing Environment' organized a conference entitled 'The genomic basis of eco-evolutionary change' and brought together experts in ecological genomics and eco-evolutionary dynamics. In this review, we use the work of the invited speakers to summarize eco-evolutionary dynamics and discuss how they are relevant for understanding and predicting responses to contemporary environmental change. Then, we show how recent advances in genomics are contributing to our understanding of eco-evolutionary dynamics. Finally, we highlight the gaps in our understanding of eco-evolutionary dynamics and recommend future avenues of research in eco-evolutionary dynamics.

Assessing genotype-phenotype associations in three dorsal colour morphs in the meadow spittlebug Philaenus spumarius (L.) (Hemiptera: Aphrophoridae) using genomic and transcriptomic resources

Background

Colour polymorphisms are common among animal species. When combined with genetic and ecological data, these polymorphisms can be excellent systems in which to understand adaptation and the molecular changes underlying phenotypic evolution. The meadow spittlebug, Philaenus spumarius (L.) (Hemiptera, Aphrophoridae), a widespread insect species in the Holarctic region, exhibits a striking dorsal colour/pattern balanced polymorphism. Although experimental crosses have revealed the Mendelian inheritance of this trait, its genetic basis remains unknown. In this study we aimed to identify candidate genomic regions associated with the colour balanced polymorphism in this species.

Results

By using restriction site-associated DNA (RAD) sequencing we were able to obtain a set of 1,837 markers across 33 individuals to test for associations with three dorsal colour phenotypes (typicus, marginellus, and trilineatus). Single and multi-association analyses identified a total of 60 SNPs associated with dorsal colour morphs. The genome size of P. spumarius was estimated by flow cytometry, revealing a 5.3 Gb genome, amongst the largest found in insects. A partial genome assembly, representing 24% of the total size, and an 81.4 Mb transcriptome, were also obtained. From the SNPs found to be associated with colour, 35% aligned to the genome and 10% to the transcriptome. Our data suggested that major loci, consisting of multi-genomic regions, may be involved in dorsal colour variation among the three dorsal colour morphs analysed. However, no homology was found between the associated loci and candidate genes known to be responsible for coloration pattern in other insect species. The associated markers showed stronger differentiation of the trilineatus colour phenotype, which has been shown previously to be more differentiated in several life-history and physiological characteristics as well. It is possible that colour variation and these traits are linked in a complex genetic architecture.

Conclusions

The loci detected to have an association with colour and the genomic and transcriptomic resources developed here constitute a basis for further research on the genetic basis of colour pattern in the meadow spittlebug P. spumarius.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0455-5) contains supplementary material, which is available to authorized users.

Observational evidence that maladaptive gene flow reduces patch occupancy in a wild insect metapopulation

Theory predicts that dispersal throughout metapopulations has a variety of consequences for the abundance and distribution of species. Immigration is predicted to increase abundance and habitat patch occupancy, but gene flow can have both positive and negative demographic consequences. Here, we address the eco-evolutionary effects of dispersal in a wild metapopulation of the stick insect Timema cristinae, which exhibits variable degrees of local adaptation throughout a heterogeneous habitat patch network of two host-plant species. To disentangle the ecological and evolutionary contributions of dispersal to habitat patch occupancy and abundance, we contrasted the effects of connectivity to populations inhabiting conspecific host plants and those inhabiting the alternate host plant. Both types of connectivity should increase patch occupancy and abundance through increased immigration and sharing of beneficial alleles through gene flow. However, connectivity to populations inhabiting the alternate host-plant species may uniquely cause maladaptive gene flow that counters the positive demographic effects of immigration. Supporting these predictions, we find the relationship between patch occupancy and alternate-host connectivity to be significantly smaller in slope than the relationship between patch occupancy and conspecific-host connectivity. Our findings illustrate the ecological and evolutionary roles of dispersal in driving the distribution and abundance of species.

A genomic perspective on the generation and maintenance of genetic diversity in herbivorous insects

Understanding the processes that generate and maintain genetic variation within populations is a central goal in evolutionary biology. Theory predicts that some of this variation is maintained as a consequence of adapting to variable habitats. Studies in herbivorous insects have played a key role in confirming this prediction. Here, we highlight theoretical and conceptual models for the maintenance of genetic diversity in herbivorous insects, empirical genomic studies testing these models, and pressing questions within the realm of evolutionary and functional genomic studies. To address key gaps, we propose an integrative approach combining population genomic scans for adaptation, genome-wide characterization of targets of selection through experimental manipulations, mapping the genetic architecture of traits influencing fitness, and functional studies. We also stress the importance of studying the maintenance of genetic variation across biological scales-from variation within populations to divergence among populations-to form a comprehensive view of adaptation in herbivorous insects.

Genome-Wide Association Studies of Grain Yield Components in Diverse Sorghum Germplasm

Grain yield and its primary determinants, grain number and weight, are important traits in cereal crops that have been well studied; however, the genetic basis of and interactions between these traits remain poorly understood. Characterization of grain yield per primary panicle (YPP), grain number per primary panicle (GNP), and 1000-grain weight (TGW) in sorghum [ (L.) Moench], a hardy C cereal with a genome size of ∼730 Mb, was implemented in a diversity panel containing 390 accessions. These accessions were genotyped to obtain 268,830 single-nucleotide polymorphisms (SNPs). Genome-wide association studies (GWAS) were performed to identify loci associated with each grain yield component and understand the genetic interactions between these traits. Genome-wide association studies identified associations across the genome with YPP, GNP, and TGW that were located within previously mapped sorghum QTL for panicle weight, grain yield, and seed size, respectively. There were no significant associations between GNP and TGW that were within 100 kb, much greater than the average linkage disequilibrium (LD) in sorghum. The identification of nonoverlapping loci for grain number and weight suggests these traits may be manipulated independently to increase the grain yield of sorghum. Following GWAS, genomic regions surrounding each associated SNP were mined for candidate genes. Previously published expression data indicated several TGW candidate genes, including an ethylene receptor homolog, were primarily expressed within developing seed tissues to support GWAS. Furthermore, maize ( L.) homologs of identified TGW candidates were differentially expressed within the seed between small- and large-kernel lines from a segregating maize population.

Whole-genome resequencing of extreme phenotypes in collared flycatchers highlights the difficulty of detecting quantitative trait loci in natural populations

Dissecting the genetic basis of phenotypic variation in natural populations is a long-standing goal in evolutionary biology. One open question is whether quantitative traits are determined only by large numbers of genes with small effects, or whether variation also exists in large-effect loci. We conducted genomewide association analyses of forehead patch size (a sexually selected trait) on 81 whole-genome-resequenced male collared flycatchers with extreme phenotypes, and on 415 males sampled independent of patch size and genotyped with a 50K SNP chip. No SNPs were genomewide statistically significantly associated with patch size. Simulation-based power analyses suggest that the power to detect large-effect loci responsible for 10% of phenotypic variance was <0.5 in the genome resequencing analysis, and <0.1 in the SNP chip analysis. Reducing the recombination by two-thirds relative to collared flycatchers modestly increased power. Tripling sample size increased power to >0.8 for resequencing of extreme phenotypes (N = 243), but power remained <0.2 for the 50K SNP chip analysis (N = 1245). At least 1 million SNPs were necessary to achieve power >0.8 when analysing 415 randomly sampled phenotypes. However, power of the 50K SNP chip to detect large-effect loci was nearly 0.8 in simulations with a small effective population size of 1500. These results suggest that reliably detecting large-effect trait loci in large natural populations will often require thousands of individuals and near complete sampling of the genome. Encouragingly, far fewer individuals and loci will often be sufficient to reliably detect large-effect loci in small populations with widespread strong linkage disequilibrium.