RAD‐sequencing for estimating genomic relatedness matrix‐based heritability in the wild: A case study in roe deer

Estimating the evolutionary potential of quantitative traits and reliably predicting responses to selection in wild populations are important challenges in evolutionary biology. The genomic revolution has opened up opportunities for measuring relatedness among individuals with precision, enabling pedigree‐free estimation of trait heritabilities in wild populations. However, until now, most quantitative genetic studies based on a genomic relatedness matrix (GRM) have focused on long‐term monitored populations for which traditional pedigrees were also available, and have often had access to knowledge of genome sequence and variability. Here, we investigated the potential of RAD‐sequencing for estimating heritability in a free‐ranging roe deer (Capreolous capreolus) population for which no prior genomic resources were available. We propose a step‐by‐step analytical framework to optimize the quality and quantity of the genomic data and explore the impact of the single nucleotide polymorphism (SNP) calling and filtering processes on the GRM structure and GRM‐based heritability estimates. As expected, our results show that sequence coverage strongly affects the number of recovered loci, the genotyping error rate and the amount of missing data. Ultimately, this had little effect on heritability estimates and their standard errors, provided that the GRM was built from a minimum number of loci (above 7,000). Genomic relatedness matrix‐based heritability estimates thus appear robust to a moderate level of genotyping errors in the SNP data set. We also showed that quality filters, such as the removal of low‐frequency variants, affect the relatedness structure of the GRM, generating lower h² estimates. Our work illustrates the huge potential of RAD‐sequencing for estimating GRM‐based heritability in virtually any natural population.     

Germline variants in oculocutaneous albinism genes and predisposition to familial cutaneous melanoma

Approximately 1%–2% of cutaneous melanoma (CM) is classified as strongly familial. We sought to investigate unexplained CM predisposition in families negative for the known susceptibility genes using next‐generation sequencing of affected individuals. Segregation of germline variants of interest within families was assessed by Sanger sequencing. Several heterozygous variants in oculocutaneous albinism (OCA) genes: TYR, OCA2, TYRP1 and SLC45A2, were present in our CM cohort. OCA is a group of autosomal recessive genetic disorders, resulting in pigmentation defects of the eyes, hair and skin. Missense variants classified as pathogenic for OCA were present in multiple families and some fully segregated with CM. The functionally compromised TYR p.T373K variant was present in three unrelated families. In OCA2, known pathogenic variants: p.V443I and p.N489D, were present in three families and one family, respectively. We identified a likely pathogenic SLC45A2 frameshift variant that fully segregated with CM in a family of four cases. Another four‐case family harboured cosegregating variants (p.A24T and p.R153C) of uncertain functional significance in TYRP1. We conclude that rare, heterozygous variants in OCA genes confer moderate risk for CM.     

Selection for barriers between honey bees and a devastating parasite

Rivaling pesticides and a dearth of flowers, the parasitic mite Varroa destructor presents a tremendous threat to western honey bees, Apis mellifera. A longstanding, but minor, pest for the Asian honey bee Apis cerana, these obligate bee parasites feast on developing and adult A. mellifera across several continents. Varroa reproduction is limited to a short window when developing bee pupae are concealed in wax cells. Mated females target developing bees just before pupation and then have about one day to initiate reproduction, eventually laying one male and up to several female offspring. Female mites often fail to reproduce at all, instead waiting in cells until their bee host finishes development and then hitching dangerous rides on a succession of adult bees for up to several weeks, before scouting for a new host pupa. In this issue of Molecular Ecology, Conlon et al. (2019) have explored mite reproductive success via a clever and thought‐provoking association study. In so doing, they have identified a protein whose actions could be integral to the dance between bees and their mite parasites.     

Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex

Chromosomal inversions can play an important role in adaptation, but the mechanism of their action in many natural populations remains unclear. An inversion could suppress recombination between locally beneficial alleles, thereby preventing maladaptive reshuffling with less‐fit, migrant alleles. The recombination suppression hypothesis has gained much theoretical support but empirical tests are lacking. Here, we evaluated the evolutionary history and phenotypic effects of a chromosomal inversion which differentiates annual and perennial forms of Mimulus guttatus. We found that perennials likely possess the derived orientation of the inversion. In addition, this perennial orientation occurs in a second perennial species, M. decorus, where it is strongly associated with life history differences between co‐occurring M. decorus and annual M. guttatus. One prediction of the recombination suppression hypothesis is that loci contributing to local adaptation will predate the inversion. To test whether the loci influencing perenniality pre‐date this inversion, we mapped QTLs for life history traits that differ between annual M. guttatus and a more distantly related, collinear perennial species, M. tilingii. Consistent with the recombination suppression hypothesis, we found that this region is associated with life history in the absence of the inversion, and this association can be broken into at least two QTLs. However, the absolute phenotypic effect of the LG8 inversion region on life history is weaker in M. tilingii than in perennials which possess the inversion. Thus, while we find support for the recombination suppression hypothesis, the contribution of this inversion to life history divergence in this group is likely complex.     

Putting the genes into community genetics

As part of the long‐term fusion of evolutionary biology and ecology (Ford, 1964), the field of community genetics has made tremendous progress in describing the impacts of plant genetic variation on community and ecosystem processes. In the “genes‐to‐ecosystems” framework (Whitham et al., 2003), genetically based traits of plant species have ecological consequences, but previous studies have not identified specific plant genes responsible for community phenotypes. The study by Barker et al. (2019) in this issue of Molecular Ecology uses an impressive common garden experiment of trembling aspen (Figure 1) to test for the genetic basis of tree traits that shape the insect community composition. Using a Genome‐Wide Association Study (GWAS), they found that genomic regions associated with phytochemical traits best explain variation in herbivore community composition, and identified specific genes associated with different types of leaf‐modifying herbivores and ants. This is one of the first studies to identify candidate genes underlying the heritable plant traits that explain patterns of insect biodiversity.     

Genomic evidence of genetic variation with pleiotropic effects on caterpillar fitness and plant traits in a model legume

Plant–insect interactions are ubiquitous, and have been studied intensely because of their relevance to damage and pollination in agricultural plants, and to the ecology and evolution of biodiversity. Variation within species can affect the outcome of these interactions. Specific genes and chemicals that mediate these interactions have been identified, but genome‐ or metabolome‐scale studies might be necessary to better understand the ecological and evolutionary consequences of intraspecific variation for plant–insect interactions. Here, we present such a study. Specifically, we assess the consequences of genome‐wide genetic variation in the model plant Medicago truncatula for Lycaeides melissa caterpillar growth and survival (larval performance). Using a rearing experiment and a whole‐genome SNP data set (>5 million SNPs), we found that polygenic variation in M. truncatula explains 9%–41% of the observed variation in caterpillar growth and survival. Genetic correlations among caterpillar performance and other plant traits, including structural defences and some anonymous chemical features, suggest that multiple M. truncatula alleles have pleiotropic effects on plant traits and caterpillar performance (or that substantial linkage disequilibrium exists among distinct loci affecting subsets of these traits). A moderate proportion of the genetic effect of M. truncatula alleles on L. melissa performance can be explained by the effect of these alleles on the plant traits we measured, especially leaf toughness. Taken together, our results show that intraspecific genetic variation in M. truncatula has a substantial effect on the successful development of L. melissa caterpillars (i.e., on a plant–insect interaction), and further point toward traits potentially mediating this genetic effect.     

Unravelling the invasion history of the Asian tiger mosquito in Europe

Multiple introductions are key features for the establishment and persistence of introduced species. However, little is known about the contribution of genetic admixture to the invasive potential of populations. To address this issue, we studied the recent invasion of the Asian tiger mosquito (Aedes albopictus) in Europe. Combining genome‐wide single nucleotide polymorphisms and historical knowledge using an approximate Bayesian computation framework, we reconstruct the colonization routes and establish the demographic dynamics of invasion. The colonization of Europe involved at least three independent introductions in Albania, North Italy and Central Italy that subsequently acted as dispersal centres throughout Europe. We show that the topology of human transportation networks shaped demographic histories with North Italy and Central Italy being the main dispersal centres in Europe. Introduction modalities conditioned the levels of genetic diversity in invading populations, and genetically diverse and admixed populations promoted more secondary introductions and have spread farther than single‐source invasions. This genomic study provides further crucial insights into a general understanding of the role of genetic diversity promoted by modern trade in driving biological invasions.     

Spatiotemporally explicit demographic modelling supports a joint effect of historical barriers to dispersal and contemporary landscape composition on structuring genomic variation in a red‐listed grasshopper

Inferring the processes underlying spatial patterns of genomic variation is fundamental to understand how organisms interact with landscape heterogeneity and to identify the factors determining species distributional shifts. Here, we use genomic data (restriction site‐associated DNA sequencing) to test biologically informed models representing historical and contemporary demographic scenarios of population connectivity for the Iberian cross‐backed grasshopper Dociostaurus hispanicus, a species with a narrow distribution that currently forms highly fragmented populations. All models incorporated biological aspects of the focal taxon that could hypothetically impact its geographical patterns of genomic variation, including (a) spatial configuration of impassable barriers to dispersal defined by topographic landscapes not occupied by the species; (b) distributional shifts resulting from the interaction between the species bioclimatic envelope and Pleistocene glacial cycles; and (c) contemporary distribution of suitable habitats after extensive land clearing for agriculture. Spatiotemporally explicit simulations under different scenarios considering these aspects and statistical evaluation of competing models within an Approximate Bayesian Computation framework supported spatial configuration of topographic barriers to dispersal and human‐driven habitat fragmentation as the main factors explaining the geographical distribution of genomic variation in the species, with no apparent impact of hypothetical distributional shifts linked to Pleistocene climatic oscillations. Collectively, this study supports that both historical (i.e., topographic barriers) and contemporary (i.e., anthropogenic habitat fragmentation) aspects of landscape composition have shaped major axes of genomic variation in the studied species and emphasizes the potential of model‐based approaches to gain insights into the temporal scale at which different processes impact the demography of natural populations.     

Heritabilities,social environment effects and genetic correlations of social behaviours in a cooperatively breeding vertebrate

Social animals interact frequently with conspecifics, and their behaviour is influenced by social context, environmental cues and the behaviours of interaction partners, allowing for adaptive, flexible adjustments to social encounters. This flexibility can be limited by part of the behavioural variation being genetically determined. Furthermore, behaviours can be genetically correlated, potentially constraining independent evolution. Understanding social behaviour thus requires carefully disentangling genetic, environmental, maternal and social sources of variations as well as the correlation structure between behaviours. Here, we assessed heritability, maternal, common environment and social effects of eight social behaviours in Neolamprologus pulcher, a cooperatively breeding cichlid. We bred wild‐caught fish in a paternal half‐sibling design and scored ability to defend a resource against conspecifics, to integrate into a group and the propensity to help defending the group territory (“helping behaviour”). We assessed genetic, social and phenotypic correlations within clusters of behaviours predicted to be functionally related, namely “competition,” “aggression,” “aggression‐sociability,” “integration” and “integration‐help.” Helping behaviour and two affiliative behaviours were heritable, whereas there was little evidence for a genetic basis in all other traits. Phenotypic social effects explained part of the variation in a sociable and a submissive behaviour, but there were no maternal or common environment effects. Genetic and phenotypic correlation within clusters was mostly positive. A group's social environment influenced covariances of social behaviours. Genetic correlations were similar in magnitude but usually exceeding the phenotypic ones, indicating that conclusions about the evolution of social behaviours in this species could be provisionally drawn from phenotypic data in cases where data for genetic analyses are unobtainable.