Genetic by environment interactions affect plant–soil linkages

The role of plant intraspecific variation in plant–soil linkages is poorly understood, especially in the context of natural environmental variation, but has important implications in evolutionary ecology. We utilized three 18‐ to 21‐year‐old common gardens across an elevational gradient, planted with replicates of five Populus angustifolia genotypes each, to address the hypothesis that tree genotype (G), environment (E), and G × E interactions would affect soil carbon and nitrogen dynamics beneath individual trees. We found that soil nitrogen and carbon varied by over 50% and 62%, respectively, across all common garden environments. We found that plant leaf litter (but not root) traits vary by genotype and environment while soil nutrient pools demonstrated genotype, environment, and sometimes G × E interactions, while process rates (net N mineralization and net nitrification) demonstrated G × E interactions. Plasticity in tree growth and litter chemistry was significantly related to the variation in soil nutrient pools and processes across environments, reflecting tight plant–soil linkages. These data overall suggest that plant genetic variation can have differential affects on carbon storage and nitrogen cycling, with implications for understanding the role of genetic variation in plant–soil feedback as well as management plans for conservation and restoration of forest habitats with a changing climate.     

Within‐host interactions shape virulence‐related traits of trematode genotypes

Within‐host interactions between co‐infecting parasites can significantly influence the evolution of key parasite traits, such as virulence (pathogenicity of infection). The type of interaction is expected to predict the direction of selection, with antagonistic interactions favouring more virulent genotypes and synergistic interactions less virulent genotypes. Recently, it has been suggested that virulence can further be affected by the genetic identity of co‐infecting partners (G × G interactions), complicating predictions on disease dynamics. Here, we used a natural host–parasite system including a fish host and a trematode parasite to study the effects of G × G interactions on infection virulence. We exposed rainbow trout (Oncorhynchus mykiss) either to single genotypes or to mixtures of two genotypes of the eye fluke Diplostomum pseudospathaceum and estimated parasite infectivity (linearly related to pathogenicity of infection, measured as coverage of eye cataracts) and relative cataract coverage (controlled for infectivity). We found that both traits were associated with complex G × G interactions, including both increases and decreases from single infection to co‐infection, depending on the genotype combination. In particular, combinations where both genotypes had low average infectivity and relative cataract coverage in single infections benefited from co‐infection, while the pattern was opposite for genotypes with higher performance. Together, our results show that infection outcomes vary considerably between single and co‐infections and with the genetic identity of the co‐infecting parasites. This can result in variation in parasite fitness and consequently impact evolutionary dynamics of host–parasite interactions.     

Phenotypic integration plasticity in Daphnia magna: an integral facet of G × E interactions

Phenotypic integration can be defined as the network of multivariate relationships among behavioural, physiological and morphological traits that describe the organism. Phenotypic integration plasticity refers to the change in patterns of phenotypic integration across environments or ontogeny. Because studies of phenotypic plasticity have predominantly focussed on single traits, a G × E interaction is typically perceived as differences in the magnitude of trait expression across two or more environments. However, many plastic responses involve coordinated responses in multiple traits, raising the possibility that relative differences in trait expression in different environments are an important, but often overlooked, source of G × E interaction. Here, we use phenotypic change vectors to statistically compare the multivariate life‐history plasticity of six Daphnia magna clones collected from four disparate European populations. Differences in the magnitude of plastic responses were statistically distinguishable for two of the six clones studied. However, differences in phenotypic integration plasticity were statistically distinguishable for all six of the clones studied, suggesting that phenotypic integration plasticity is an important component of G × E interactions that may be missed unless appropriate multivariate analyses are used.     

Heritability of male attractiveness persists despite evidence for unreliable sexual signals in Drosophila simulans

Sexual signals can be used to attract mates, but to be honest indicators of signaller quality they need to convey information reliably. However, environmental variation and genotype‐by‐environment (G × E) interactions have the potential to compromise the reliability of sexual signals. Here, we test the reliability of cuticular hydrocarbons (CHCs) as signals of heritable aspects of male attractiveness in Drosophila simulans. We examined the heritability of male attractiveness and a measure of the difference between fathers' and sons' CHC profiles across dietary and temperature environments. Our results show that environmental heterogeneity disrupts the similarity of some components of father and son CHC profile. However, overall male attractiveness is heritable within and across environments, so that sire attractiveness is a good predictor of son attractiveness even with environmental heterogeneity. This suggests that although some male CHC signals are unreliable, attractive genotypes retain their attractiveness across environments on average.     

Genotype‐by‐genotype specificity remains robust to average temperature variation in an aphid/endosymbiont/parasitoid system

Genotype‐by‐genotype interactions demonstrate the existence of variation upon which selection acts in host–parasite systems at respective resistance and infection loci. These interactions can potentially be modified by environmental factors, which would entail that different genotypes are selected under different environmental conditions. In the current study, we checked for a G × G × E interaction in the context of average temperature and the genotypes of asexual lines of the endoparasitoid wasp Lysiphlebus fabarum and isolates of Hamiltonella defensa, a protective secondary endosymbiont of the wasp's host, the black bean aphid Aphis fabae. We exposed genetically identical aphids harbouring different isolates of H. defensa to three asexual lines of the parasitoid and measured parasitism success under three different temperatures (15, 22 and 29 °C). Although there was clear evidence for increased susceptibility to parasitoids at the highest average temperature and a strong G × G interaction between the host's symbionts and the parasitoids, no modifying effect of temperature, that is, no significant G × G × E interaction, was detected. This robustness of the observed specificity suggests that the relative fitness of different parasitoid genotypes on hosts protected by particular symbionts remains uncomplicated by spatial or temporal variation in temperature, which should facilitate biological control strategies.     

Polymorphisms in a desaturase 2 ortholog associate with cuticular hydrocarbon and male mating success variation in a natural population of Drosophila serrata

Elucidating the nature of genetic variation underlying both sexually selected traits and the fitness components of sexual selection is essential to understanding the broader consequences of sexual selection as an evolutionary process. To date, there have been relatively few attempts to connect the genetic variance in sexually selected traits with segregating DNA sequence polymorphisms. We set out to address this in a well‐characterized sexual selection system – the cuticular hydrocarbons (CHCs) of Drosophila serrata – using an indirect association study design that allowed simultaneous estimation of the genetic variance in CHCs, sexual fitness and single nucleotide polymorphism (SNP) effects in an outbred population. We cloned and sequenced an ortholog of the D. melanogaster desaturase 2 gene, previously shown to affect CHC biosynthesis in D. melanogaster, and associated 36 SNPs with minor allele frequencies > 0.02 with variance in CHCs and sexual fitness. Three SNPs had significant multivariate associations with CHC phenotype (q‐value < 0.05). At these loci, minor alleles had multivariate effects on CHCs that were weakly associated with the multivariate direction of sexual selection operating on these traits. Two of these SNPs had pleiotropic associations with male mating success, suggesting these variants may underlie responses to sexual selection due to this locus. There were 15 significant male mating success associations (q‐value < 0.1), and interestingly, we detected a nonrandom pattern in the relationship between allele frequency and direction of effect on male mating success. The minor‐frequency allele usually reduced male mating success, suggesting a positive association between male mating success and total fitness at this locus.     

Family‐level variation in early life‐cycle traits of kelp

Temperate kelp forests (Laminarians) are threatened by temperature stress due to ocean warming and photoinhibition due to increased light associated with canopy loss. However, the potential for evolutionary adaptation in kelp to rapid climate change is not well known. This study examined family‐level variation in physiological and photosynthetic traits in the early life‐cycle stages of the ecologically important Australasian kelp Ecklonia radiata and the response of E. radiata families to different temperature and light environments using a family × environment design. There was strong family‐level variation in traits relating to morphology (surface area measures, branch length, branch count) and photosynthetic performance (F_v/F_m) in both haploid (gametophyte) and diploid (sporophyte) stages of the life‐cycle. Additionally, the presence of family × environment interactions showed that offspring from different families respond differently to temperature and light in the branch length of male gametophytes and oogonia surface area of female gametophytes. Negative responses to high temperatures were stronger for females vs. males. Our findings suggest E. radiata may be able to respond adaptively to climate change but studies partitioning the narrow vs. broad sense components of heritable variation are needed to establish the evolutionary potential of E. radiata to adapt under climate change.     

The complex interplay between macronutrient intake,cuticular hydrocarbon expression and mating success in male decorated crickets

The condition dependence of male sexual traits plays a central role in sexual selection theory. Relatively little, however, is known about the condition dependence of chemical signals used in mate choice and their subsequent effects on male mating success. Furthermore, few studies have isolated the specific nutrients responsible for condition‐dependent variation in male sexual traits. Here, we used nutritional geometry to determine the effect of protein (P) and carbohydrate (C) intake on male cuticular hydrocarbon (CHC) expression and mating success in male decorated crickets (Gryllodes sigillatus). We show that both traits are maximized at a moderate‐to‐high intake of nutrients in a P:C ratio of 1 : 1.5. We also show that female precopulatory mate choice exerts a complex pattern of linear and quadratic sexual selection on this condition‐dependent variation in male CHC expression. Structural equation modelling revealed that although the effect of nutrient intake on mating success is mediated through condition‐dependent CHC expression, it is not exclusively so, suggesting that other traits must also play an important role. Collectively, our results suggest that the complex interplay between nutrient intake, CHC expression and mating success plays an important role in the operation of sexual selection in G. sigillatus.     

Premating isolation is determined by larval rearing substrates in cactophilic Drosophila mojavensis. IX. Host plant and population specific epicuticular hydrocarbon expression influences mate choice and sexual selection

Sexual signals in cactophilic Drosophila mojavensis include cuticular hydrocarbons (CHCs), contact pheromones that mediate female discrimination of males during courtship. CHCs, along with male courtship songs, cause premating isolation between diverged populations, and are influenced by genotype × environment interactions caused by different host cacti. CHC profiles of mated and unmated adult flies from a Baja California and a mainland Mexico population of D. mojavensis reared on two host cacti were assayed to test the hypothesis that male CHCs mediate within‐population female discrimination of males. In multiple choice courtship trials, mated and unmated males differed in CHC profiles, indicating that females prefer males with particular blends of CHCs. Mated and unmated females significantly differed in CHC profiles as well. Adults in the choice trials had CHC profiles that were significantly different from those in pair‐mated adults from no‐choice trials revealing an influence of sexual selection. Females preferred different male CHC blends in each population, but the influence of host cactus on CHC variation was significant only in the mainland population indicating population‐specific plasticity in CHCs. Different groups of CHCs mediated female choice‐based sexual selection in each population suggesting that geographical and ecological divergence has the potential to promote divergence in mate communication systems.     

Why does the magnitude of genotype‐by‐environment interaction vary?

Genotype‐by‐environment interaction (G × E), that is, genetic variation in phenotypic plasticity, is a central concept in ecology and evolutionary biology. G×E has wide‐ranging implications for trait development and for understanding how organisms will respond to environmental change. Although G × E has been extensively documented, its presence and magnitude vary dramatically across populations and traits. Despite this, we still know little about why G × E is so evident in some traits and populations, but minimal or absent in others. To encourage synthetic research in this area, we review diverse hypotheses for the underlying biological causes of variation in G × E. We extract common themes from these hypotheses to develop a more synthetic understanding of variation in G × E and suggest some important next steps.     

Breeding for robustness: investigating the genotype‐by‐environment interaction and micro‐environmental sensitivity of Genetically Improved Farmed Tilapia (Oreochromis niloticus)

Robustness has become a highly desirable breeding goal in the globalized agricultural market. Both genotype‐by‐environment interaction (G × E) and micro‐environmental sensitivity are important robustness components of aquaculture production, in which breeding stock is often disseminated to different environments. The objectives of this study were (i) to quantify the degree of G × E by assessing the growth performance of Genetically Improved Farmed Tilapia (GIFT) across three countries (Malaysia, India and China) and (ii) to quantify the genetic heterogeneity of environmental variance for body weight at harvest (BW) in GIFT as a measure of micro‐environmental sensitivity. Selection for BW was carried out for 13 generations in Malaysia. Subsets of 60 full‐sib families from Malaysia were sent to China and India after five and nine generations respectively. First, a multi‐trait animal model was used to analyse the BW in different countries as different traits. The results indicate a strong G × E. Second, a genetically structured environmental variance model, implemented using Bayesian inference, was used to analyse micro‐environmental sensitivity of BW in each country. The analysis revealed the presence of genetic heterogeneity of both BW and its environmental variance in all environments. The presence of genetic variation in residual variance of BW implies that the residual variance can be modified by selection. Incorporating both G × E and micro‐environmental sensitivity information may help in selecting robust genotypes with high performance across environments and resilience to environmental fluctuations.     

Genetic divergence and phenotypic plasticity contribute to variation in cuticular hydrocarbons in the seaweed fly Coelopa frigida

Cuticular hydrocarbons (CHCs) form the boundary between insects and their environments and often act as essential cues for species, mate, and kin recognition. This complex polygenic trait can be highly variable both among and within species, but the causes of this variation, especially the genetic basis, are largely unknown. In this study, we investigated phenotypic and genetic variation of CHCs in the seaweed fly, Coelopa frigida, and found that composition was affected by both genetic (sex and population) and environmental (larval diet) factors. We subsequently conducted behavioral trials that show CHCs are likely used as a sexual signal. We identified general shifts in CHC chemistry as well as individual compounds and found that the methylated compounds, mean chain length, proportion of alkenes, and normalized total CHCs differed between sexes and populations. We combined these data with whole genome resequencing data to examine the genetic underpinnings of these differences. We identified 11 genes related to CHC synthesis and found population‐level outlier SNPs in 5 that are concordant with phenotypic differences. Together these results reveal that the CHC composition of C. frigida is dynamic, strongly affected by the larval environment, and likely under natural and sexual selection.     

Evolution of sexually dimorphic pheromone profiles coincides with increased number of male‐specific chemosensory organs in Drosophila prolongata

Binary communication systems that involve sex‐specific signaling and sex‐specific signal perception play a key role in sexual selection and in the evolution of sexually dimorphic traits. The driving forces and genetic changes underlying such traits can be investigated in systems where sex‐specific signaling and perception have emerged recently and show evidence of potential coevolution. A promising model is found in Drosophila prolongata, which exhibits a species‐specific increase in the number of male chemosensory bristles. We show that this transition coincides with recent evolutionary changes in cuticular hydrocarbon (CHC) profiles. Long‐chain CHCs that are sexually monomorphic in the closest relatives of D. prolongata (D. rhopaloa, D. carrolli, D. kurseongensis, and D. fuyamai) are strongly male‐biased in this species. We also identify an intraspecific female‐limited polymorphism, where some females have male‐like CHC profiles. Both the origin of sexually dimorphic CHC profiles and the female‐limited polymorphism in D. prolongata involve changes in the relative amounts of three mono‐alkene homologs, 9‐tricosene, 9‐pentacosene, and 9‐heptacosene, all of which share a common biosynthetic origin and point to a potentially simple genetic change underlying these traits. Our results suggest that pheromone synthesis may have coevolved with chemosensory perception and open the way for reconstructing the origin of sexual dimorphism in this communication system.     

Between‐sex genetic covariance constrains the evolution of sexual dimorphism in Drosophila melanogaster

Males and females share much of their genome, and as a result, intralocus sexual conflict is generated when selection on a shared trait differs between the sexes. This conflict can be partially or entirely resolved via the evolution of sex‐specific genetic variation that allows each sex to approach, or possibly achieve, its optimum phenotype, thereby generating sexual dimorphism. However, shared genetic variation between the sexes can impose constraints on the independent expression of a shared trait in males and females, hindering the evolution of sexual dimorphism. Here, we examine genetic constraints on the evolution of sexual dimorphism in Drosophila melanogaster cuticular hydrocarbon (CHC) expression. We use the extended G matrix, which includes the between‐sex genetic covariances that constitute the B matrix, to compare genetic constraints on two sets of CHC traits that differ in the extent of their sexual dimorphism. We find significant genetic constraints on the evolution of further dimorphism in the least dimorphic traits, but no such constraints for the most dimorphic traits. We also show that the genetic constraints on the least dimorphic CHCs are asymmetrical between the sexes. Our results suggest that there is evidence both for resolved and ongoing sexual conflict in D. melanogaster CHC profiles.     

Investigating the production of sexual resting structures in a plant pathogen reveals unexpected self‐fertility and genotype‐by‐environment effects

The sexual stage of pathogens governs recombination patterns and often also provides means of surviving the off‐season. Despite its importance for evolutionary potential and between‐season epidemiology, sexual systems have not been carefully investigated for many important pathogens, and what generates variation in successful sexual reproduction of pathogens remains unexplored. We surveyed the sexually produced resting structures (chasmothecia) across 86 natural populations of fungal pathogen Podosphaera plantaginis (Ascomycota) naturally infecting Plantago lanceolata in the Åland archipelago, southwestern Finland. For this pathosystem, these resting structures are a key life‐history stage, as more than half of the local pathogen populations go extinct every winter. We uncovered substantial variation in the level of chasmothecia produced among populations, ranging from complete absence to presence on all infected leaves. We found that chasmothecia developed within clonal isolates (single‐strain cultures). Additionally, these clonal isolates all contained both MAT1‐1‐1 and MAT1‐2‐1 genes that characterize mating types in Ascomycetes. Hence, contrary to expectations, we conclude that this species is capable of haploid selfing. In controlled inoculations, we discovered that pathogen genotypes varied in their tendency to produce chasmothecia. Production of chasmothecia was also affected by ambient temperature (E) and by the interaction between temperature and pathogen genotype (G × E). These G, E and G × E effects found both at a European scale and within the Åland archipelago may partly explain the high variability observed among populations in chasmothecia levels. Consequently, they may be key drivers of the evolutionary potential and epidemiology of this highly dynamic pathosystem.     

Population genomics reveals a fine‐scale recombination landscape for genetic improvement of cotton

Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population‐scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population‐scaled recombination rate between Gossypium hirsutum and G. arboreum and sub‐genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome‐wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high‐recombination regions than were those of yield and early maturity traits. These results provide insights into the population‐scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.     

Genome‐wide association and genomic prediction for biomass yield in a genetically diverse Miscanthus sinensis germplasm panel phenotyped at five locations in Asia and North America

To improve the efficiency of breeding of Miscanthus for biomass yield, there is a need to develop genomics‐assisted selection for this long‐lived perennial crop by relating genotype to phenotype and breeding value across a broad range of environments. We present the first genome‐wide association (GWA) and genomic prediction study of Miscanthus that utilizes multilocation phenotypic data. A panel of 568 Miscanthus sinensis accessions was genotyped with 46,177 single nucleotide polymorphisms (SNPs) and evaluated at one subtropical and five temperate locations over 3 years for biomass yield and 14 yield‐component traits. GWA and genomic prediction were performed separately for different years of data in order to assess reproducibility. The analyses were also performed for individual field trial locations, as well as combined phenotypic data across groups of locations. GWA analyses identified 27 significant SNPs for yield, and a total of 504 associations across 298 unique SNPs across all traits, sites, and years. For yield, the greatest number of significant SNPs was identified by combining phenotypic data across all six locations. For some of the other yield‐component traits, greater numbers of significant SNPs were obtained from single site data, although the number of significant SNPs varied greatly from site to site. Candidate genes were identified. Accounting for population structure, genomic prediction accuracies for biomass yield ranged from 0.31 to 0.35 across five northern sites and from 0.13 to 0.18 for the subtropical location, depending on the estimation method. Genomic prediction accuracies of all traits were similar for single‐location and multilocation data, suggesting that genomic selection will be useful for breeding broadly adapted M. sinensis as well as M. sinensis optimized for specific climates. All of our data, including DNA sequences flanking each SNP, are publicly available. By facilitating genomic selection in M. sinensis and Miscanthus × giganteus, our results will accelerate the breeding of these species for biomass in diverse environments.