The evolution and adaptive potential of transcriptional variation in sticklebacks--signatures of selection and widespread heritability

Genetic architecture and correlations between the gut microbiome and gut gene transcription in Chinook salmon (Oncorhynchus tshawytscha)

Population divergence through selection can drive local adaptation in natural populations which has implications for the effective restoration of declining and extirpated populations. However, adaptation to local environmental conditions is complicated when both the host and its associated microbiomes must respond via co-evolutionary change. Nevertheless, for adaptation to occur through selection, variation in both host and microbiome traits should include additive genetic effects. Here we focus on host immune function and quantify factors affecting variation in gut immune gene transcription and gut bacterial community composition in early life-stage Chinook salmon (Oncorhynchus tshawytscha). Specifically, we utilized a replicated factorial breeding design to determine the genetic architecture (sire, dam and sire-by-dam interaction) of gut immune gene transcription and microbiome composition. Furthermore, we explored correlations between host gut gene transcription and microbiota composition. Gene transcription was quantified using nanofluidic qPCR arrays (22 target genes) and microbiota composition using 16 S rRNA gene (V5-V6) amplicon sequencing. We discovered limited but significant genetic architecture in gut microbiota composition and transcriptional profiles. We also identified significant correlations between gut gene transcription and microbiota composition, highlighting potential mechanisms for functional interactions between the two. Overall, this study provides support for the co-evolution of host immune function and their gut microbiota in Chinook salmon, a species recognized as locally adapted. Thus, the inclusion of immune gene transcription profile and gut microbiome composition as factors in the development of conservation and commercial rearing practices may provide new and more effective approaches to captive rearing.

Genomic insights into selection for heterozygous alleles and woody traits in Populus tomentosa

Heterozygous alleles are widespread in outcrossing and clonally propagated woody plants. The variation in heterozygosity that underlies population adaptive evolution and phenotypic variation, however, remains largely unknown. Here, we describe a de novo chromosome-level genome assembly of Populus tomentosa, an economic and ecologically important native tree in northern China. By resequencing 302 natural accessions, we determined that the South subpopulation (Pop_S) encompasses the ancestral strains of P. tomentosa, while the Northwest subpopulation (Pop_NW) and Northeast subpopulation (Pop_NE) experienced different selection pressures during population evolution, resulting in significant population differentiation and a decrease in the extent of heterozygosity. Analysis of heterozygous selective sweep regions (HSSR) suggested that selection for lower heterozygosity contributed to the local adaptation of P. tomentosa by dwindling gene expression and genetic load in the Pop_NW and Pop_NE subpopulations. Genome-wide association studies (GWAS) revealed that 88 single nucleotide polymorphisms (SNPs) within 63 genes are associated with nine wood composition traits. Among them, the selection for the homozygous AA allele in PtoARF8 is associated with reductions in cellulose and hemicellulose contents by attenuating PtoARF8 expression, and the increase in lignin content is attributable to the selection for decreases in exon heterozygosity in PtoLOX3 during adaptive evolution of natural populations. This study provides novel insights into allelic variations in heterozygosity associated with adaptive evolution of P. tomentosa in response to the local environment and identifies a series of key genes for wood component traits, thereby facilitating genomic-based breeding of important traits in perennial woody plants.

Population transcriptogenomics highlights impaired metabolism and small population sizes in tree frogs living in the Chernobyl Exclusion Zone

BACKGROUND

Individual functional modifications shape the ability of wildlife populations to cope with anthropogenic environmental changes. But instead of adaptive response, human-altered environments can generate a succession of deleterious functional changes leading to the extinction of the population. To study how persistent anthropogenic changes impacted local species' population status, we characterised population structure, genetic diversity and individual response of gene expression in the tree frog Hyla orientalis along a gradient of radioactive contamination around the Chernobyl nuclear power plant.

RESULTS

We detected lower effective population size in populations most exposed to ionizing radiation in the Chernobyl Exclusion Zone that is not compensated by migrations from surrounding areas. We also highlighted a decreased body condition of frogs living in the most contaminated area, a distinctive transcriptomics signature and stop-gained mutations in genes involved in energy metabolism. While the association with dose will remain correlational until further experiments, a body of evidence suggests the direct or indirect involvement of radiation exposure in these changes.

CONCLUSIONS

Despite ongoing migration and lower total dose rates absorbed than at the time of the accident, our results demonstrate that Hyla orientalis specimens living in the Chernobyl Exclusion Zone are still undergoing deleterious changes, emphasizing the long-term impacts of the nuclear disaster.

The lack of genetic variation underlying thermal transcriptomic plasticity suggests limited adaptability of the Northern shrimp, Pandalus borealis

Introduction

Genetic variation underlies the populations’ potential to adapt to and persist in a changing environment, while phenotypic plasticity can play a key role in buffering the negative impacts of such change at the individual level.

Methods

We investigated the role of genetic variation in the thermal response of the northern shrimp Pandalus borealis, an ectotherm species distributed in the Arctic and North Atlantic Oceans. More specifically, we estimated the proportion transcriptomic responses explained by genetic variance of female shrimp from three origins after 30 days of exposure to three temperature treatments.

Results

We characterized the P. borealis transcriptome (170,377 transcripts, of which 27.48% were functionally annotated) and then detected a total of 1,607 and 907 differentially expressed transcripts between temperatures and origins, respectively. Shrimp from different origins displayed high but similar level of transcriptomic plasticity in response to elevated temperatures. Differences in transcript expression among origins were not correlated to population genetic differentiation or diversity but to environmental conditions at origin during sampling.

Discussion

The lack of genetic variation explaining thermal plasticity suggests limited adaptability in this species’ response to future environmental changes. These results together with higher mortality observed at the highest temperature indicate that the thermal niche of P. borealis will likely be restricted to higher latitudes in the future. This prediction concurs with current decreases in abundance observed at the southern edge of this species geographical distribution, as it is for other cold-adapted crustaceans.

Predicting recombination frequency from map distance

Map distance is one of the key measures in genetics and indicates the expected number of crossovers between two loci. Map distance is estimated from the observed recombination frequency using mapping functions, the most widely used of those, Haldane and Kosambi, being developed at the time when the number of markers was low and unobserved crossovers had a substantial effect on the recombination fractions. In contemporary high-density marker data, the probability of multiple crossovers between adjacent loci is negligible and different mapping functions yield the same result, that is, the recombination frequency between adjacent loci is equal to the map distance in Morgans. However, high-density linkage maps contain an interpretation problem: the map distance over a long interval is additive and its association with recombination frequency is not defined. Here, we demonstrate with high-density linkage maps from humans and stickleback fishes that the inverses of Haldane's and Kosambi's mapping functions systematically underpredict recombination frequencies from map distance. To remedy this, we formulate a piecewise function that yields more accurate predictions of recombination frequency from map distance. Our results demonstrate that the association between map distance and recombination frequency is context-dependent and without a universal solution.

Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

Detecting signatures of selection on gene expression

A substantial amount of phenotypic diversity results from changes in gene expression levels and patterns. Understanding how the transcriptome evolves is therefore a key priority in identifying mechanisms of adaptive change. However, in contrast to powerful models of sequence evolution, we lack a consensus model of gene expression evolution. Furthermore, recent work has shown that many of the comparative approaches used to study gene expression are subject to biases that can lead to false signatures of selection. Here we first outline the main approaches for describing expression evolution and their inherent biases. Next, we bridge the gap between the fields of phylogenetic comparative methods and transcriptomics to reinforce the main pitfalls of inferring selection on expression patterns and use simulation studies to show that shifts in tissue composition can heavily bias inferences of selection. We close by highlighting the multi-dimensional nature of transcriptional variation and identifying major unanswered questions in disentangling how selection acts on the transcriptome.

Comparison of Gene Expression Between Resistant and Susceptible Families Against VPAHPND and Identification of Biomarkers Used for Resistance Evaluation in Litopenaeus vannamei

Acute hepatopancreatic necrosis disease (AHPND) has caused a heavy loss to shrimp aquaculture since its outbreak. Vibrio parahaemolyticus (VP_AHPND) is regarded as one of the main pathogens that caused AHPND in the Pacific white shrimp Litopenaeus vannamei. In order to learn more about the mechanism of resistance to AHPND, the resistant and susceptible shrimp families were obtained through genetic breeding, and comparative transcriptome approach was used to analyze the gene expression patterns between resistant and susceptible families. A total of 95 families were subjected to VP_AHPND challenge test, and significant variations in the resistance of these families were observed. Three pairs of resistant and susceptible families were selected for transcriptome sequencing. A total of 489 differentially expressed genes (DEGs) that presented in at least two pairwise comparisons were screened, including 196 DEGs highly expressed in the susceptible families and 293 DEGs in the resistant families. Among these DEGs, 16 genes demonstrated significant difference in all three pairwise comparisons. Gene set enrichment analysis (GSEA) of all 27,331 expressed genes indicated that some energy metabolism processes were enriched in the resistant families, while signal transduction and immune system were enriched in the susceptible families. A total of 32 DEGs were further confirmed in the offspring of the detected families, among which 19 genes were successfully verified. The identified genes in this study will be useful for clarifying the genetic mechanism of shrimp resistance against Vibrio and will further provide molecular markers for evaluating the disease resistance of shrimp in the breeding program.

Gene expression responses to thermal shifts in the endangered lichen Lobaria pulmonaria

Anthropogenic climate change has led to unprecedented shifts in temperature across many ecosystems. In a context of rapid environmental changes, acclimation is an important process as it may influence the capacity of organisms to survive under novel thermal conditions. Mechanisms of acclimation could involve upregulation of stress response genes involved in protein folding, DNA damage repair and the regulation of signal transduction genes, along with a simultaneous downregulation of genes involved in growth or the cell cycle, in order to maintain cellular functions and equilibria. We transplanted Lobaria pulmonaria lichens originating from different forests to determine the relative effects of long-term acclimation and genetic factors on the variability in expression of mycobiont and photobiont genes. We found a strong response of the mycobiont and photobiont to high temperatures, regardless of sample origin. The green-algal photobiont had an overall lower response than the mycobiont. Gene expression of both symbionts was also influenced by acclimation to transplantation sites and by genetic factors. L. pulmonaria seems to have evolved powerful molecular pathways to deal with environmental fluctuations and stress and can acclimate to new habitats by transcriptomic convergence. Although L. pulmonaria has the molecular machinery to counteract short-term thermal stress, survival of lichens such as L. pulmonaria depends mostly on their long-term positive carbon balance, which can be compromised by higher temperatures and reduced precipitation, and both these outcomes have been predicted for Central Europe in connection with global climate change.

Local adaptation contributes to gene expression divergence in maize

Gene expression links genotypes to phenotypes, so identifying genes whose expression is shaped by selection will be important for understanding the traits and processes underlying local adaptation. However, detecting local adaptation for gene expression will require distinguishing between divergence due to selection and divergence due to genetic drift. Here, we adapt a QST−FST framework to detect local adaptation for transcriptome-wide gene expression levels in a population of diverse maize genotypes. We compare the number and types of selected genes across a wide range of maize populations and tissues, as well as selection on cold-response genes, drought-response genes, and coexpression clusters. We identify a number of genes whose expression levels are consistent with local adaptation and show that genes involved in stress response show enrichment for selection. Due to its history of intense selective breeding and domestication, maize evolution has long been of interest to researchers, and our study provides insight into the genes and processes important for in local adaptation of maize.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

Automated improvement of stickleback reference genome assemblies with Lep-Anchor software

We describe an integrative approach to improve contiguity and haploidy of a reference genome assembly and demonstrate its impact with practical examples. With two novel features of Lep-Anchor software and a combination of dense linkage maps, overlap detection and bridging long reads, we generated an improved assembly of the nine-spined stickleback (Pungitius pungitius) reference genome. We were able to remove a significant number of haplotypic contigs, detect more genetic variation and improve the contiguity of the genome, especially that of X chromosome. However, improved scaffolding cannot correct for mosaicism of erroneously assembled contigs, demonstrated by a de novo assembly of a 1.6-Mbp inversion. Qualitatively similar gains were obtained with the genome of three-spined stickleback (Gasterosteus aculeatus). Since the utility of genome-wide sequencing data in biological research depends heavily on the quality of the reference genome, the improved and fully automated approach described here should be helpful in refining reference genome assemblies.

Examining adaptive evolution of immune activity: opportunities provided by gastropods in the age of 'omics'

Parasites threaten all free-living organisms, including molluscs. Understanding the evolution of immune defence traits in natural host populations is crucial for predicting their long-term performance under continuous infection risk. Adaptive trait evolution requires that traits are subject to selection (i.e. contribute to organismal fitness) and that they are heritable. Despite broad interest in the evolutionary ecology of immune activity in animals, the understanding of selection on and evolutionary potential of immune defence traits is far from comprehensive. For instance, empirical observations are only rarely in line with theoretical predictions of immune activity being subject to stabilizing selection. This discrepancy may be because ecoimmunological studies can typically cover only a fraction of the complexity of an animal immune system. Similarly, molecular immunology/immunogenetics studies provide a mechanistic understanding of immunity, but neglect variation that arises from natural genetic differences among individuals and from environmental conditions. Here, we review the current literature on natural selection on and evolutionary potential of immune traits in animals, signal how merging ecological immunology and genomics will strengthen evolutionary ecological research on immunity, and indicate research opportunities for molluscan gastropods for which well-established ecological understanding and/or 'immune-omics' resources are already available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Transcriptome profiling of Lymnaea stagnalis (Gastropoda) for ecoimmunological research

BACKGROUND

Host immune function can contribute to numerous ecological/evolutionary processes. Ecoimmunological studies, however, typically use one/few phenotypic immune assays and thus do not consider the complexity of the immune system. Therefore, "omics" resources that allow quantifying immune activity across multiple pathways are needed for ecoimmunological models. We applied short-read based RNAseq (Illumina NextSeq 500, PE-81) to characterise transcriptome profiles of Lymnaea stagnalis (Gastropoda), a multipurpose model snail species. We used a genetically diverse snail stock and exposed individuals to immune elicitors (injury, bacterial/trematode pathogens) and changes in environmental conditions that can alter immune activity (temperature, food availability).

RESULTS

Immune defence factors identified in the de novo assembly covered elements broadly described in other gastropods. For instance, pathogen-recognition receptors (PRR) and lectins activate Toll-like receptor (TLR) pathway and cytokines that regulate cellular and humoral defences. Surprisingly, only modest diversity of antimicrobial peptides and fibrinogen related proteins were detected when compared with other taxa. Additionally, multiple defence factors that may contribute to the phenotypic immune assays used to quantify antibacterial activity and phenoloxidase (PO)/melanisation-type reaction in this species were found. Experimental treatments revealed factors from non-self recognition (lectins) and signalling (TLR pathway, cytokines) to effectors (e.g., antibacterial proteins, PO enzymes) whose transcription depended on immune stimuli and environmental conditions, as well as components of snail physiology/metabolism that may drive these effects. Interestingly, the transcription of many factors (e.g., PRR, lectins, cytokines, PO enzymes, antibacterial proteins) showed high among-individual variation.

CONCLUSIONS

Our results indicate several uniform aspects of gastropod immunity, but also apparent differences between L. stagnalis and some previously examined taxa. Interestingly, in addition to immune defence factors that responded to immune elicitors and changes in environmental conditions, many factors showed high among-individual variation across experimental snails. We propose that such factors are highly important to be included in future ecoimmunological studies because they may be the key determinants of differences in parasite resistance among individuals both within and between natural snail populations.

Mating-regulated atrial proteases control reinsemination rates in Anopheles gambiae females

Anopheles gambiae mosquitoes are the most important vectors of human malaria. The reproductive success of these mosquitoes relies on a single copulation event after which the majority of females become permanently refractory to further mating. This refractory behavior is at least partially mediated by the male-synthetized steroid hormone 20-hydroxyecdysone (20E), which is packaged together with other seminal secretions into a gelatinous mating plug and transferred to the female atrium during mating. In this study, we show that two 20E-regulated chymotrypsin-like serine proteases specifically expressed in the reproductive tract of An. gambiae females play an important role in modulating the female susceptibility to mating. Silencing these proteases by RNA interference impairs correct plug processing and slows down the release of the steroid hormone 20E from the mating plug. In turn, depleting one of these proteases, the Mating Regulated Atrial Protease 1 (MatRAP1), reduces female refractoriness to further copulation, so that a significant proportion of females mate again. Microscopy analysis reveals that MatRAP1 is localized on a previously undetected peritrophic matrix-like structure surrounding the mating plug. These data provide novel insight into the molecular mechanisms shaping the post-mating biology of these important malaria vectors.

Latitudinal variation in climate-associated genes imperils range edge populations

The ecological impacts of increasing global temperatures are evident in most ecosystems on Earth, but our understanding of how climatic variation influences natural selection and adaptive resilience across latitudes remains largely unknown. Latitudinal gradients allow testing general ecosystem-level theories relevant to climatic adaptation. We assessed differences in adaptive diversity of populations along a latitudinal region spanning highly variable temperate to subtropical climates. We generated and integrated information from environmental mapping, phenotypic variation and genome-wide data from across the geographical range of the rainbowfish Melanotaenia duboulayi, an emerging aquatic system for studies of climate change. We detected, after controlling for spatial population structure, strong interactions between genotypes and environment associated with variation in stream flow and temperature. Some of these hydroclimate-associated genes were found to interact within functional protein networks that contain genes of adaptive significance for projected future climates in rainbowfish. Hydroclimatic selection was also associated with variation in phenotypic traits, including traits known to affect fitness of rainbowfish exposed to different flow environments. Consistent with predictions from the "climatic variability hypothesis," populations exposed to extremes of important environmental variables showed stronger adaptive divergence and less variation in climate-associated genes compared to populations at the centre of the environmental gradient. Our findings suggest that populations that evolved at environmental range margins and at geographical range edges may be more vulnerable to changing climates, a finding with implications for predicting adaptive resilience and managing biodiversity under climate change.

Genetic Constraints, Transcriptome Plasticity, and the Evolutionary Response to Climate Change

In situ adaptation to climate change will be critical for the persistence of many ectotherm species due to their relative lack of dispersal capacity. Climate change is causing increases in both the mean and the variance of environmental temperature, each of which may act as agents of selection on different traits. Importantly, these traits may not be heritable or have the capacity to evolve independently from one another. When genetic constraints prevent the "baseline" values of thermal performance traits from evolving rapidly, phenotypic plasticity driven by gene expression might become critical. We review the literature for evidence that thermal performance traits in ectotherms are heritable and have genetic architectures that permit their unconstrained evolution. Next, we examine the relationship between gene expression and both the magnitude and duration of thermal stress. Finally, we identify genes that are likely to be important for adaptation to a changing climate and determine whether they show patterns consistent with thermal adaptation. Although few studies have measured narrow-sense heritabilities of thermal performance traits, current evidence suggests that the end points of thermal reaction norms (tolerance limits) are moderately heritable and have the potential to evolve rapidly. By contrast, performance at intermediate temperatures has substantially lower evolutionary potential. Moreover, evolution in many species appears to be constrained by genetic correlations such that populations can adapt to either increases in mean temperature or temperature variability, but not both. Finally, many species have the capacity for plastic expression of the transcriptome in response to temperature shifts, with the number of differentially expressed genes increasing with the magnitude, but not the duration, of thermal stress. We use these observations to develop a conceptual model that describes the likely trajectory of genome evolution in response to changes in environmental temperature. Our results indicate that extreme weather events, rather than gradual increases in mean temperature, are more likely to drive genetic and phenotypic change in wild ectotherms.

Evolution and Plasticity of the Transcriptome Under Temperature Fluctuations in the Fungal Plant Pathogen Zymoseptoria tritici

Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.

Highly dynamic transcriptional reprogramming and shorter isoform shifts under acute stresses during biological invasions

Phenotypic plasticity has been increasingly recognized for its importance in adaptation to novel environments, and initial rapid plastic response to acute stresses usually serves as the stepping stone for future adaptation. Differential gene expression and alternative splicing have been proposed as two underlying mechanisms for rapid plastic response to environmental stresses. Here, we used an invasive model species, Ciona savignyi, to investigate the temporary plastic changes under temperature stresses on gene expression and alternative splicing. Our results revealed rapid and highly dynamic gene expression reprogramming and alternative splicing switch under acute stresses. Distinct transcriptional response profiles were triggered by two types of temperature stresses, showing resilience recovery and increasing divergence under heat and cold challenges, respectively. Interestingly, alternative exons were more inclined to be skipped under both heat and cold stresses, leading to shorter isoforms but with maintained Open Reading Frames (ORFs). Although similar response patterns were observed between differential gene expression and alternative splicing, low overlap between Differentially Expressed Genes (DEGs) and Differentially Alternative Spliced Genes (DASGs) suggests that distinct gene sets and associated functions should be involved in temperature challenges. Thus, alternative splicing should offer an additional layer of plastic response to environmental challenges. Finally, we identified key plastic genes involved in both gene expression regulation and alternative splicing. The results obtained here shed light on adaptation and accommodation mechanisms during biological invasions, particularly for acute environmental changes at early stages of biological invasions such as transport and introduction.

Additive and mostly adaptive plastic responses of gene expression to multiple stress in Tribolium castaneum

Gene expression is known to be highly responsive to the environment and important for adjustment of metabolism but there is also growing evidence that differences in gene regulation contribute to species divergence and differences among locally adapted populations. However, most studies so far investigated populations when divergence had already occurred. Selection acting on expression levels at the onset of adaptation to an environmental change has not been characterized. Understanding the mechanisms is further complicated by the fact that environmental change is often multivariate, meaning that organisms are exposed to multiple stressors simultaneously with potentially interactive effects. Here we use a novel approach by combining fitness and whole-transcriptome data in a large-scale experiment to investigate responses to drought, heat and their combination in Tribolium castaneum. We found that fitness was reduced by both stressors and their combined effect was almost additive. Expression data showed that stressor responses were acting independently and did not interfere physiologically. Since we measured expression and fitness within the same individuals, we were able to estimate selection on gene expression levels. We found that variation in fitness can be attributed to gene expression variation and that selection pressures were environment dependent and opposite between control and stress conditions. We could further show that plastic responses of expression were largely adaptive, i.e. in the direction that should increase fitness.

Genome-Wide Genotype-Expression Relationships Reveal Both Copy Number and Single Nucleotide Differentiation Contribute to Differential Gene Expression between Stickleback Ecotypes

Repeated and independent emergence of trait divergence that matches habitat differences is a sign of parallel evolution by natural selection. Yet, the molecular underpinnings that are targeted by adaptive evolution often remain elusive. We investigate this question by combining genome-wide analyses of copy number variants (CNVs), single nucleotide polymorphisms (SNPs), and gene expression across four pairs of lake and river populations of the three-spined stickleback (Gasterosteus aculeatus). We tested whether CNVs that span entire genes and SNPs occurring in putative cis-regulatory regions contribute to gene expression differences between sticklebacks from lake and river origins. We found 135 gene CNVs that showed a significant positive association between gene copy number and gene expression, suggesting that CNVs result in dosage effects that can fuel phenotypic variation and serve as substrates for habitat-specific selection. Copy number differentiation between lake and river sticklebacks also contributed to expression differences of two immune-related genes in immune tissues, cathepsin A and GIMAP7. In addition, we identified SNPs in cis-regulatory regions (eSNPs) associated with the expression of 1,865 genes, including one eSNP upstream of a carboxypeptidase gene where both the SNP alleles differentiated and the gene was differentially expressed between lake and river populations. Our study highlights two types of mutations as important sources of genetic variation involved in the evolution of gene expression and in potentially facilitating repeated adaptation to novel environments.

Modularity Facilitates Flexible Tuning of Plastic and Evolutionary Gene Expression Responses during Early Divergence

Gene expression changes have been recognized as important drivers of adaptation to changing environmental conditions. Little is known about the relative roles of plastic and evolutionary responses in complex gene expression networks during the early stages of divergence. Large gene expression data sets coupled with in silico methods for identifying coexpressed modules now enable systems genetics approaches also in nonmodel species for better understanding of gene expression responses during early divergence. Here, we combined gene coexpression analyses with population genetics to separate plastic and population (evolutionary) effects in expression networks using small salmonid populations as a model system. We show that plastic and population effects were highly variable among the six identified modules and that the plastic effects explained larger proportion of the total eigengene expression than population effects. A more detailed analysis of the population effects using a Q_ST - F_ST comparison across 16,622 annotated transcripts revealed that gene expression followed neutral expectations within modules and at the global level. Furthermore, two modules showed enrichment for genes coding for early developmental traits that have been previously identified as important phenotypic traits in thermal responses in the same model system indicating that coexpression analysis can capture expression patterns underlying ecologically important traits. We suggest that module-specific responses may facilitate the flexible tuning of expression levels to local thermal conditions. Overall, our study indicates that plasticity and neutral evolution are the main drivers of gene expression variance in the early stages of thermal adaptation in this system.

Convergent phenotypic evolution of the visual system via different molecular routes: How Neotropical cichlid fishes adapt to novel light environments

How predictable is evolution? This remains a fundamental but contested issue in evolutionary biology. When independent lineages colonize the same environment, we are presented with a natural experiment that allows us to ask if genetic and ecological differences promote species-specific evolutionary outcomes or whether species phenotypically evolve in a convergent manner in response to shared selection pressures. If so, are the molecular mechanisms underlying phenotypic convergence the same? In Nicaragua, seven species of cichlid fishes concurrently colonized two novel photic environments. Hence, their visual system represents a compelling model to address these questions, particularly since the adaptive value of phenotypic changes is well-understood. By analyzing retinal transcriptomes, we found that differential expression of genes responsible for color vision (cone opsins and cyp27c1) produced rapid and mostly convergent changes of predicted visual sensitivities. Notably, these changes occurred in the same direction in all species although there were differences in underlying gene expression patterns illustrating nonconvergence at the molecular level. Adaptive phenotypes evolved deterministically, even when species differ substantially in ecology and genetic variation. This provides strong evidence that phenotypic evolution of the visual system occurred in response to similar selective forces of the photic environment.

Analysis of natural female post-mating responses of Anopheles gambiae and Anopheles coluzzii unravels similarities and differences in their reproductive ecology

Anopheles gambiae and An. coluzzii, the two most important malaria vectors in sub-Saharan Africa, are recently radiated sibling species that are reproductively isolated even in areas of sympatry. In females from these species, sexual transfer of male accessory gland products, including the steroid hormone 20-hydroxyecdysone (20E), induces vast behavioral, physiological, and transcriptional changes that profoundly shape their post-mating ecology, and that may have contributed to the insurgence of post-mating, prezygotic reproductive barriers. As these barriers can be detected by studying transcriptional changes induced by mating, we set out to analyze the post-mating response of An. gambiae and An. coluzzii females captured in natural mating swarms in Burkina Faso. While the molecular pathways shaping short- and long-term mating-induced changes are largely conserved in females from the two species, we unravel significant inter-specific differences that suggest divergent regulation of key reproductive processes such as egg development, processing of seminal secretion, and mating behavior, that may have played a role in reproductive isolation. Interestingly, a number of these changes occur in genes previously shown to be regulated by the sexual transfer of 20E and may be due to divergent utilization of this steroid hormone in the two species.

Strong phenotypic plasticity limits potential for evolutionary responses to climate change

Phenotypic plasticity, the expression of multiple phenotypes from one genome, is a widespread adaptation to short-term environmental fluctuations, but whether it facilitates evolutionary adaptation to climate change remains contentious. Here, we investigate seasonal plasticity and adaptive potential in an Afrotropical butterfly expressing distinct phenotypes in dry and wet seasons. We assess the transcriptional architecture of plasticity in a full-factorial analysis of heritable and environmental effects across 72 individuals, and reveal pervasive gene expression differences between the seasonal phenotypes. Strikingly, intra-population genetic variation for plasticity is largely absent, consistent with specialisation to a particular environmental cue reliably predicting seasonal transitions. Under climate change, deteriorating accuracy of predictive cues will likely aggravate maladaptive phenotype-environment mismatches and increase selective pressures on reaction norms. However, the observed paucity of genetic variation for plasticity limits evolutionary responses, potentially weakening prospects for population persistence. Thus, seasonally plastic species may be especially vulnerable to climate change.

Environmental and genetic determinants of transcriptional plasticity in Chinook salmon

Variation in gene transcription is widely believed to be the mechanistic basis of phenotypically plastic traits; however, comparatively little is known about the inheritance patterns of transcriptional variation that would allow us to predict its response to selection. In addition, acclimation to different environmental conditions influences acute transcriptional responses to stress and it is unclear if these effects are heritable. To address these gaps in knowledge, we assayed levels of messenger RNA for 14 candidate genes at rest and in response to a 24-h confinement stress for 72 half-sib families of Chinook salmon reared in two different environments (hatchery and semi-natural stream channel). We observed extensive plasticity for mRNA levels of metabolic and stress response genes and demonstrated that mRNA level plasticity due to rearing environment affects mRNA level plasticity in response to stress. These effects have important implications for natural populations experiencing multiple stressors. We identified genotype-by-environment interactions for mRNA levels that were dominated by maternal effects; however, mRNA level response to challenge also exhibited a non-additive genetic basis. Our results indicate that while plasticity for mRNA levels can evolve, predicting the outcome of selection will be difficult. The inconsistency in genetic architecture among treatment groups suggests there is considerable cryptic genetic variation for gene expression.

Host Genotype and Microbiota Contribute Asymmetrically to Transcriptional Variation in the Threespine Stickleback Gut

Recent studies of interactions between hosts and their resident microbes have revealed important ecological and evolutionary consequences that emerge from these complex interspecies relationships, including diseases that occur when the interactions go awry. Given the preponderance of these interactions, we hypothesized that effects of the microbiota on gene expression in the developing gut—an important aspect of host biology—would be pervasive, and that these effects would be both comparable in magnitude to and contingent on effects of the host genetic background. To evaluate the effects of the microbiota, host genotype, and their interaction on gene expression in the gut of a genetically diverse, gnotobiotic host model, the threespine stickleback (Gasterosteus aculeatus), we compared RNA-seq data among 84 larval fish. Surprisingly, we found that stickleback population and family differences explained substantially more gene expression variation than the presence of microbes. Expression levels of 72 genes, however, were affected by our microbiota treatment. These genes, including many associated with innate immunity, comprise a tractable subset of host genetic factors for precise, systems-level study of host–microbe interactions in the future. Importantly, our data also suggest subtle signatures of a statistical interaction between host genotype and the microbiota on expression patterns of genetic pathways associated with innate immunity, coagulation and complement cascades, focal adhesion, cancer, and peroxisomes. These genotype-by-environment interactions may prove to be important leads to the understanding of host genetic mechanisms commonly at the root of sometimes complex molecular relationships between hosts and their resident microbes.

Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish

Understanding genomic signatures of divergent selection underlying long-term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated with lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated with inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated with inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat stress and immune responses. Outcomes of gene-temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene-environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.

The roles of plasticity and evolutionary change in shaping gene expression variation in natural populations of extremophile fish

The notorious plasticity of gene expression responses and the complexity of environmental gradients complicate the identification of adaptive differences in gene regulation among populations. We combined transcriptome analyses in nature with common-garden and exposure experiments to establish cause-effect relationships between the presence of a physiochemical stressor and expression differences, as well as to test how evolutionary change and plasticity interact to shape gene expression variation in natural systems. We studied two evolutionarily independent population pairs of an extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulphide (H₂ S)-rich springs and adjacent nontoxic habitats and assessed genomewide expression patterns of wild-caught and common-garden-raised individuals exposed to different concentrations of H₂ S. We found that 7.7% of genes that were differentially expressed between sulphidic and nonsulphidic ecotypes remained differentially expressed in the laboratory, indicating that sources of selection other than H₂ S-or plastic responses to other environmental factors-contribute substantially to gene expression patterns observed in the wild. Concordantly differentially expressed genes in the wild and the laboratory were primarily associated with H₂ S detoxification, sulphur processing and metabolic physiology. While shared, ancestral plasticity played a minor role in shaping gene expression variation observed in nature, we documented evidence for evolved population differences in the constitutive expression as well as the H₂ S inducibility of candidate genes. Mechanisms underlying gene expression variation also varied substantially across the two ecotype pairs. These results provide a springboard for studying evolutionary modifications of gene regulatory mechanisms that underlie expression variation in locally adapted populations.

Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops

The use of transcriptome data in the study of the population genetics of a species can capture faint signals of both genetic variation and expression variation and can provide a broad picture of a species' genomic response to environmental conditions. In this study, we characterized the genetic and expression diversity of Miscanthus lutarioriparius by comparing more than 16,225 transcripts obtained from 78 individuals, belonging to 10 populations distributed across the species' entire geographic range. We only observed a low level of nucleotide diversity (π = 0.000434) among the transcriptome data of these populations, which is consistent with highly conserved sequences of functional elements and protein-coding genes captured with this method. Tests of population divergence using the transcriptome data were consistent with previous microsatellite data but proved to be more sensitive, particularly if gene expression variation was considered as well. For example, the analysis of expression data showed that genes involved in photosynthetic processes and responses to temperature or reactive oxygen species stimuli were significantly enriched in certain populations. This differential gene expression was primarily observed among populations and not within populations. Interestingly, nucleotide diversity was significantly negatively correlated with expression diversity within populations, while this correlation was positive among populations. This suggests that genetic and expression variation play separate roles in adaptation and population persistence. Combining analyses of genetic and gene expression variation represents a promising approach for studying the population genetics of wild species and may uncover both adaptive and nonadaptive processes.

Challenges and advances for transcriptome assembly in non-model species

Analyses of high-throughput transcriptome sequences of non-model organisms are based on two main approaches: de novo assembly and genome-guided assembly using mapping to assign reads prior to assembly. Given the limits of mapping reads to a reference when it is highly divergent, as is frequently the case for non-model species, we evaluate whether using blastn would outperform mapping methods for read assignment in such situations (>15% divergence). We demonstrate its high performance by using simulated reads of lengths corresponding to those generated by the most common sequencing platforms, and over a realistic range of genetic divergence (0% to 30% divergence). Here we focus on gene identification and not on resolving the whole set of transcripts (i.e. the complete transcriptome). For simulated datasets, the transcriptome-guided assembly based on blastn recovers 94.8% of genes irrespective of read length at 0% divergence; however, assignment rate of reads is negatively correlated with both increasing divergence level and reducing read lengths. Nevertheless, we still observe 92.6% of recovered genes at 30% divergence irrespective of read length. This analysis also produces a categorization of genes relative to their assignment, and suggests guidelines for data processing prior to analyses of comparative transcriptomics and gene expression to minimize potential inferential bias associated with incorrect transcript assignment. We also compare the performances of de novo assembly alone vs in combination with a transcriptome-guided assembly based on blastn both via simulation and empirically, using data from a cyprinid fish species and from an oak species. For any simulated scenario, the transcriptome-guided assembly using blastn outperforms the de novo approach alone, including when the divergence level is beyond the reach of traditional mapping methods. Combining de novo assembly and a related reference transcriptome for read assignment also addresses the bias/error in contigs caused by the dependence on a related reference alone. Empirical data corroborate these findings when assembling transcriptomes from the two non-model organisms: Parachondrostoma toxostoma (fish) and Quercus pubescens (plant). For the fish species, out of the 31,944 genes known from D. rerio, the guided and de novo assemblies recover respectively 20,605 and 20,032 genes but the performance of the guided assembly approach is much higher for both the contiguity and completeness metrics. For the oak, out of the 29,971 genes known from Vitis vinifera, the transcriptome-guided and de novo assemblies display similar performance, but the new guided approach detects 16,326 genes where the de novo assembly only detects 9,385 genes.

Heritable gene expression differences between lake and stream stickleback include both parallel and antiparallel components

The repeated phenotypic patterns that characterize populations undergoing parallel evolution provide support for a deterministic role of adaptation by natural selection. Determining the level of parallelism also at the genetic level is thus central to our understanding of how natural selection works. Many studies have looked for repeated genomic patterns in natural populations, but work on gene expression is less common. The studies that have examined gene expression have found some support for parallelism, but those studies almost always used samples collected from the wild that potentially confounds the effects of plasticity with heritable differences. Here we use two independent pairs of lake and stream threespine stickleback (Gasterosteus aculeatus) raised in common garden conditions to assess both parallel and antiparallel (that is, similar versus different directions of lake–stream expression divergence in the two watersheds) heritable gene expression differences as measured by total RNA sequencing. We find that more genes than expected by chance show either parallel (22 genes, 0.18% of expressed genes) or antiparallel (24 genes, 0.20% of expressed genes) lake–stream expression differences. These results correspond well with previous genomic studies in stickleback ecotype pairs that found similar levels of parallelism. We suggest that parallelism might be similarly constrained at the genomic and transcriptomic levels.

Structure and stability of genetic variance-covariance matrices: A Bayesian sparse factor analysis of transcriptional variation in the three-spined stickleback

The genetic variance-covariance matrix (G) is a quantity of central importance in evolutionary biology due to its influence on the rate and direction of multivariate evolution. However, the predictive power of empirically estimated G-matrices is limited for two reasons. First, phenotypes are high-dimensional, whereas traditional statistical methods are tuned to estimate and analyse low-dimensional matrices. Second, the stability of G to environmental effects and over time remains poorly understood. Using Bayesian sparse factor analysis (BSFG) designed to estimate high-dimensional G-matrices, we analysed levels variation and covariation in 10,527 expressed genes in a large (n = 563) half-sib breeding design of three-spined sticklebacks subject to two temperature treatments. We found significant differences in the structure of G between the treatments: heritabilities and evolvabilities were higher in the warm than in the low-temperature treatment, suggesting more and faster opportunity to evolve in warm (stressful) conditions. Furthermore, comparison of G and its phenotypic equivalent P revealed the latter is a poor substitute of the former. Most strikingly, the results suggest that the expected impact of G on evolvability-as well as the similarity among G-matrices-may depend strongly on the number of traits included into analyses. In our results, the inclusion of only few traits in the analyses leads to underestimation in the differences between the G-matrices and their predicted impacts on evolution. While the results highlight the challenges involved in estimating G, they also illustrate that by enabling the estimation of large G-matrices, the BSFG method can improve predicted evolutionary responses to selection.

Transcriptional mechanisms underlying life-history responses to climate change in the three-spined stickleback

Phenotypic plasticity, the ability of an organism to express different phenotypes depending on the environment, provides an important mechanism by which an animal population can persist under rapid climate change. We experimentally tested both life‐history and transcriptional responses of an ecological model species, the three‐spined stickleback, to warm acclimation at the southern edge of its European range. We explored cross‐environment genetic correlations of key life‐history traits in male sticklebacks exposed to long‐term temperature changes to examine whether the plasticity pattern was variable among genotypes by using a character‐state approach. We also studied gene expression plasticity by analysing both whole‐transcriptome and candidate gene expression in brain and liver. Male sticklebacks that developed under warmer conditions during winter were smaller in size and invested less in nuptial coloration at the beginning of the breeding season, showing similar responses across different genotypes. The lack of genetic variation in life‐history responses may limit any future evolution of the thermal reaction norm in the study population. After long‐term exposure to increased winter temperatures, genes responsible for several metabolic and oxidation–reduction processes were upregulated, and some hormone genes involved in growth and reproduction were downregulated in the brain. In the liver, there was no significantly represented gene ontology by the differentially expressed genes. Since a higher temperature leads to a higher resting metabolic rate, living in warmer environments may incur higher energetic costs for ectotherms to maintain cellular homoeostasis, resulting in negative consequences for life‐history traits. The expression of genes related to metabolism, cellular homoeostasis and regulatory signalling may underlie temperature‐induced changes in life history.

Large-scale SNP screenings identify markers linked with GCRV resistant traits through transcriptomes of individuals and cell lines in Ctenopharyngodon idella

Grass carp (Ctenopharyngodon idella) is an important economic species in freshwater aquaculture and its industry has been confined due to variety degeneration and frequent diseases. Marker-assisted selection is a feasible method for selective breeding of new varieties. Transcriptome data have greatly facilitated high-throughput single nucleotide polymorphism (SNP) marker discovery and phenotype association study. In this study, we gained a total of 25,981 and 5,775 high quality SNPs in two transcriptomes from individuals and cell lines, respectively. Comparative transcriptome analysis identified 413 and 832 grass carp reovirus (GCRV)-resistant-association SNPs as well as 1,381 and 1,606 GCRV-susceptible-association SNPs in individuals and cell lines, respectively. Integrated analysis indicated 22 genes with single SNP share common resistant/susceptible traits in two transcriptomes. Furthermore, we infected grass carp with GCRV, genotyping and association analyses were performed, and 9 in 22 SNPs were confirmed by PCR-RFLP. Meanwhile, mRNA expression profiles of 6 genes containing confirmed SNPs were examined by qRT-PCR. The results demonstrated that mRNA expressions were significant differences in resistant/susceptible individuals and cell lines. The present study develops an important strategy for high throughput screening of phenotype association genetic markers and the results will serve in grass carp breeding for GCRV resistance.

Loggerhead sea turtle embryos (Caretta caretta) regulate expression of stress response and developmental genes when exposed to a biologically realistic heat stress

Oviparous reptile embryos are expected to breach their critical thermal maxima if temperatures reach those predicted under current climate change models due to the lack of the maternal buffering processes and parental care. Heat-shock proteins (HSPs) are integral in the molecular response to thermal stress, and their expression is heritable, but the roles of other candidate families such as the heat-shock factors (HSFs) have not been determined in reptiles. Here, we subject embryonic sea turtles (Caretta caretta) to a biologically realistic thermal stress and employ de novo transcriptomic profiling of brain tissue to investigate the underlying molecular response. From a reference transcriptome of 302 293 transcripts, 179 were identified as differentially expressed between treatments. As anticipated, genes enriched in the heat-shock treatment were primarily associated with the Hsp families, or were genes whose products play similar protein editing and chaperone functions (e.g. bag3, MYOC and serpinh1). Unexpectedly, genes encoding the HSFs were not significantly upregulated under thermal stress, indicating their presence in unstressed cells in an inactive state. Genes that were downregulated under thermal stress were less well functionally defined but were associated with stress response, development and cellular organization, suggesting that developmental processes may be compromised at realistically high temperatures. These results confirm that genes from the Hsp families play vital roles in the thermal tolerance of developing reptile embryos and, in addition with a number of other genes, should be targets for evaluating the capacity of oviparous reptiles to respond adaptively to the effects of climate change.

Gene co-expression network connectivity is an important determinant of selective constraint

While several studies have investigated general properties of the genetic architecture of natural variation in gene expression, few of these have considered natural, outbreeding populations. In parallel, systems biology has established that a general feature of biological networks is that they are scale-free, rendering them buffered against random mutations. To date, few studies have attempted to examine the relationship between the selective processes acting to maintain natural variation of gene expression and the associated co-expression network structure. Here we utilised RNA-Sequencing to assay gene expression in winter buds undergoing bud flush in a natural population of Populus tremula, an outbreeding forest tree species. We performed expression Quantitative Trait Locus (eQTL) mapping and identified 164,290 significant eQTLs associating 6,241 unique genes (eGenes) with 147,419 unique SNPs (eSNPs). We found approximately four times as many local as distant eQTLs, with local eQTLs having significantly higher effect sizes. eQTLs were primarily located in regulatory regions of genes (UTRs or flanking regions), regardless of whether they were local or distant. We used the gene expression data to infer a co-expression network and investigated the relationship between network topology, the genetic architecture of gene expression and signatures of selection. Within the co-expression network, eGenes were underrepresented in network module cores (hubs) and overrepresented in the periphery of the network, with a negative correlation between eQTL effect size and network connectivity. We additionally found that module core genes have experienced stronger selective constraint on coding and non-coding sequence, with connectivity associated with signatures of selection. Our integrated genetics and genomics results suggest that purifying selection is the primary mechanism underlying the genetic architecture of natural variation in gene expression assayed in flushing leaf buds of P. tremula and that connectivity within the co-expression network is linked to the strength of purifying selection.

Different contributions of local- and distant-regulatory changes to transcriptome divergence between stickleback ecotypes

Differential gene expression can play an important role in phenotypic evolution and divergent adaptation. Although differential gene expression can be caused by both local- and distant-regulatory changes, we know little about their relative contribution to transcriptome evolution in natural populations. Here, we conducted expression quantitative trait loci (eQTL) analysis to investigate the genetic architecture underlying transcriptome divergence between marine and stream ecotypes of threespine sticklebacks (Gasterosteus aculeatus). We identified both local and distant eQTLs, some of which constitute hotspots, regions with a disproportionate number of significant eQTLs relative to the genomic background. The majority of local eQTLs including those in the hotspots caused expression changes consistent with the direction of transcriptomic divergence between ecotypes. Genome scan analysis showed that many local eQTLs overlapped with genomic regions of high differentiation. In contrast, nearly half of the distant eQTLs including those in the hotspots caused opposite expression changes, and few overlapped with regions of high differentiation, indicating that distant eQTLs may act as a constraint of transcriptome evolution. Finally, a comparison between two salinity conditions revealed that nearly half of eQTL hotspots were environment specific, suggesting that analysis of genetic architecture in multiple conditions is essential for predicting response to selection.

Regulatory Architecture of Gene Expression Variation in the Threespine Stickleback Gasterosteus aculeatus

Much adaptive evolutionary change is underlain by mutational variation in regions of the genome that regulate gene expression rather than in the coding regions of the genes themselves. An understanding of the role of gene expression variation in facilitating local adaptation will be aided by an understanding of underlying regulatory networks. Here, we characterize the genetic architecture of gene expression variation in the threespine stickleback (Gasterosteus aculeatus), an important model in the study of adaptive evolution. We collected transcriptomic and genomic data from 60 half-sib families using an expression microarray and genotyping-by-sequencing, and located expression quantitative trait loci (eQTL) underlying the variation in gene expression in liver tissue using an interval mapping approach. We identified eQTL for several thousand expression traits. Expression was influenced by polymorphism in both cis- and trans-regulatory regions. Trans-eQTL clustered into hotspots. We did not identify master transcriptional regulators in hotspot locations: rather, the presence of hotspots may be driven by complex interactions between multiple transcription factors. One observed hotspot colocated with a QTL recently found to underlie salinity tolerance in the threespine stickleback. However, most other observed hotspots did not colocate with regions of the genome known to be involved in adaptive divergence between marine and freshwater habitats.

The heritability of chimpanzee and human brain asymmetry

Human brains are markedly asymmetric in structure and lateralized in function, which suggests a relationship between these two properties. The brains of other closely related primates, such as chimpanzees, show similar patterns of asymmetry, but to a lesser degree, indicating an increase in anatomical and functional asymmetry during hominin evolution. We analysed the heritability of cerebral asymmetry in chimpanzees and humans using classic morphometrics, geometric morphometrics, and quantitative genetic techniques. In our analyses, we separated directional asymmetry and fluctuating asymmetry (FA), which is indicative of environmental influences during development. We show that directional patterns of asymmetry, those that are consistently present in most individuals in a population, do not have significant heritability when measured through simple linear metrics, but they have marginally significant heritability in humans when assessed through three-dimensional configurations of landmarks that reflect variation in the size, position, and orientation of different cortical regions with respect to each other. Furthermore, genetic correlations between left and right hemispheres are substantially lower in humans than in chimpanzees, which points to a relatively stronger environmental influence on left-right differences in humans. We also show that the level of FA has significant heritability in both species in some regions of the cerebral cortex. This suggests that brain responsiveness to environmental influences, which may reflect neural plasticity, has genetic bases in both species. These results have implications for the evolvability of brain asymmetry and plasticity among humans and our close relatives.

Epistasis and Pleiotropy Affect the Modularity of the Genotype-Phenotype Map of Cross-Resistance in HIV-1

The genotype–phenotype (GP) map is a central concept in evolutionary biology as it describes the mapping of molecular genetic variation onto phenotypic trait variation. Our understanding of that mapping remains partial, especially when trying to link functional clustering of pleiotropic gene effects with patterns of phenotypic trait co-variation. Only on rare occasions have studies been able to fully explore that link and tend to show poor correspondence between modular structures within the GP map and among phenotypes. By dissecting the structure of the GP map of the replicative capacity of HIV-1 in 15 drug environments, we provide a detailed view of that mapping from mutational pleiotropic variation to phenotypic co-variation, including epistatic effects of a set of amino-acid substitutions in the reverse transcriptase and protease genes. We show that epistasis increases the pleiotropic degree of single mutations and provides modularity to the GP map of drug resistance in HIV-1. Moreover, modules of epistatic pleiotropic effects within the GP map match the phenotypic modules of correlated replicative capacity among drug classes. Epistasis thus increases the evolvability of cross-resistance in HIV by providing more drug- and class-specific pleiotropic profiles to the main effects of the mutations. We discuss the implications for the evolution of cross-resistance in HIV.

Transcriptome sequencing reveals population differentiation in gene expression linked to functional traits and environmental gradients in the South African shrub Protea repens

Understanding the environmental and genetic mechanisms underlying locally adaptive trait variation across the ranges of species is a major focus of evolutionary biology. Combining transcriptome sequencing with common garden experiments on populations spanning geographical and environmental gradients holds promise for identifying such mechanisms. The South African shrub Protea repens displays diverse phenotypes in the wild along drought and temperature gradients. We grew plants from seeds collected at 19 populations spanning this species' range, and sequenced the transcriptomes of these plants to reveal gene pathways associated with adaptive trait variation. We related expression in co-expressed gene networks to trait phenotypes measured in the common garden and to source population climate. We found that expression in gene networks correlated with source-population environment and with plant traits. In particular, the activity of gene networks enriched for growth related pathways correlated strongly with source site minimum winter temperature and with leaf size, stem diameter and height in the garden. Other gene networks with enrichments for photosynthesis related genes showed associations with precipitation. Our results strongly suggest that this species displays population-level differences in gene expression that have been shaped by source population site climate, and that are reflected in trait variation along environmental gradients.

The adaptive potential of subtropical rainbowfish in the face of climate change: heritability and heritable plasticity for the expression of candidate genes

Whilst adaptation and phenotypic plasticity might buffer species against habitat degradation associated with global climate change, few studies making such claims also possess the necessary and sufficient data to support them. Doing so requires demonstration of heritable variation in traits affecting fitness under new environmental conditions. We address this issue using an emerging aquatic system to study adaptation to climate change, the crimson-spotted rainbowfish (Melanotaenia duboulayi), a freshwater species from a region of eastern Australia projected to be affected by marked temperature increases. Captive born M. duboulayi of known pedigree were used to assess the long-term effects of contemporary and 2070-projected summer temperatures on the expression of genes previously identified in a climate change transcriptomics (RNA-Seq) experiment. Nearly all genes responded to increasing temperature. Significant additive genetic variance explained a moderate proportion of transcriptional variation for all genes. Most genes also showed broad-sense genetic variation in transcriptional plasticity. Additionally, molecular pathways of candidate genes co-occur with genes inferred to be under climate-mediated selection in wild M. duboulayi populations. Together, these results indicate the presence of existing variation in important physiological traits, and the potential for adaptive responses to a changing thermal environment.