Divergence in gene regulation at young life history stages of whitefish (Coregonus sp.) and the emergence of genomic isolation

P-elements strengthen reproductive isolation within the Drosophila simulans species complex

Determining mechanisms that underlie reproductive isolation (RI) is key to understanding how species boundaries are maintained in nature. Transposable elements (TEs) are ubiquitous across eukaryotic genomes. However, the role of TEs in modulating the strength of RI between species is poorly understood. Several species of Drosophila have been found to harbor P-elements (PEs), yet only D. simulans is known to be currently polymorphic for their presence in wild populations. PEs can cause RI between PE-containing (P) and PE-lacking (M) lineages of the same species. However, it is unclear whether they also contribute to the magnitude of RI between species. Here, we use the simulans species complex to assess whether differences in PE status between D. simulans and its sister species, which do not harbor PEs, contribute to multiple barriers to gene flow between species. We show that crosses involving a P D. simulans father and an M mother from a sister species exhibit lower F1 female fecundity than crosses involving an M D. simulans father and an M sister-species mother. We also find that another TE, I-element, might play a minor role in determining the frequency of dysgenesis between species. Our results suggest that the presence of PEs in a species can strengthen isolation from its sister species, providing evidence that TEs can play a role in RI.

Differences among families in craniofacial shape at early life-stages of Arctic charr (Salvelinus alpinus)

BACKGROUND

Organismal fitness can be determined at early life-stages, but phenotypic variation at early life-stages is rarely considered in studies on evolutionary diversification. The trophic apparatus has been shown to contribute to sympatric resource-mediated divergence in several taxa. However, processes underlying diversification in trophic traits are poorly understood. Using phenotypically variable Icelandic Arctic charr (Salvelinus alpinus), we reared offspring from multiple families under standardized laboratory conditions and tested to what extent family (i.e. direct genetic and maternal effects) contributes to offspring morphology at hatching (H) and first feeding (FF). To understand the underlying mechanisms behind early life-stage variation in morphology, we examined how craniofacial shape varied according to family, offspring size, egg size and candidate gene expression.

RESULTS

Craniofacial shape (i.e. the Meckel's cartilage and hyoid arch) was more variable between families than within families both across and within developmental stages. Differences in craniofacial morphology between developmental stages correlated with offspring size, whilst within developmental stages only shape at FF correlated with offspring size, as well as female mean egg size. Larger offspring and offspring from females with larger eggs consistently had a wider hyoid arch and contracted Meckel's cartilage in comparison to smaller offspring.

CONCLUSIONS

This study provides evidence for family-level variation in early life-stage trophic morphology, indicating the potential for parental effects to facilitate resource polymorphism.

Many Options, Few Solutions: Over 60 My Snakes Converged on a Few Optimal Venom Formulations

Gene expression changes contribute to complex trait variations in both individuals and populations. However, the evolution of gene expression underlying complex traits over macroevolutionary timescales remains poorly understood. Snake venoms are proteinaceous cocktails where the expression of each toxin can be quantified and mapped to a distinct genomic locus and traced for millions of years. Using a phylogenetic generalized linear mixed model, we analyzed expression data of toxin genes from 52 snake species spanning the 3 venomous snake families and estimated phylogenetic covariance, which acts as a measure of evolutionary constraint. We find that evolution of toxin combinations is not constrained. However, although all combinations are in principle possible, the actual dimensionality of phylomorphic space is low, with envenomation strategies focused around only four major toxin families: metalloproteases, three-finger toxins, serine proteases, and phospholipases A2. Although most extant snakes prioritize either a single or a combination of major toxin families, they are repeatedly recruited and lost. We find that over macroevolutionary timescales, the venom phenotypes were not shaped by phylogenetic constraints, which include important microevolutionary constraints such as epistasis and pleiotropy, but more likely by ecological filtering that permits a small number of optimal solutions. As a result, phenotypic optima were repeatedly attained by distantly related species. These results indicate that venoms evolve by selection on biochemistry of prey envenomation, which permit diversity through parallelism, and impose strong limits, since only a few of the theoretically possible strategies seem to work well and are observed in extant snakes.

Mate choice and gene expression signatures associated with nutritional adaptation in the medfly (Ceratitis capitata)

Evolutionary responses to nutrition are key to understanding host shifts and the resulting potential for reproductive isolation. Experimental evolution has previously been used to describe the responses of the medfly (Ceratitis capitata) to larval diets with different nutritional properties. Within 30 generations this led to divergence in larval development time, egg to adult survival and adaptation in adult body size. Here we used mRNA-seq to identify differences in gene expression patterns in these same populations, using males from the 60^th generation of nutritional selection. We validated differential expression by using qRT-PCR and found that genes linked to metabolism, oxidative phosphorylation and proteolysis were significantly over-represented among the differentially expressed genes. The results provide the first genome-wide survey of the putative mechanisms underpinning evolved responses to nutritional adaptation. In addition, we tested the hypothesis that nutritional adaptation can alter mating patterns. We found evidence for assortative mating by diet at generation 60, but not 90. Hence, the pattern was variable across generations and there was no evidence overall for any isolating mating divergence between the lines. Overall, the results provide insight into the mechanisms underpinning dietary adaptation and extend our knowledge of which traits represent core responses to nutritional selection.

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.

Between semelparity and iteroparity: Empirical evidence for a continuum of modes of parity

The number of times an organism reproduces (i.e., its mode of parity) is a fundamental life-history character, and evolutionary and ecological models that compare the relative fitnesses of different modes of parity are common in life-history theory and theoretical biology. Despite the success of mathematical models designed to compare intrinsic rates of increase (i.e., density-independent growth rates) between annual-semelparous and perennial-iteroparous reproductive schedules, there is widespread evidence that variation in reproductive allocation among semelparous and iteroparous organisms alike is continuous. This study reviews the ecological and molecular evidence for the continuity and plasticity of modes of parity-that is, the idea that annual-semelparous and perennial-iteroparous life histories are better understood as endpoints along a continuum of possible strategies. I conclude that parity should be understood as a continuum of different modes of parity, which differ by the degree to which they disperse or concentrate reproductive effort in time. I further argue that there are three main implications of this conclusion: (1) that seasonality should not be conflated with parity; (2) that mathematical models purporting to explain the general evolution of semelparous life histories from iteroparous ones (or vice versa) should not assume that organisms can only display either an annual-semelparous life history or a perennial-iteroparous one; and (3) that evolutionary ecologists should base explanations of how different life-history strategies evolve on the physiological or molecular basis of traits underlying different modes of parity.

Adaptation and acclimation of traits associated with swimming capacity in Lake Whitefish (coregonus clupeaformis) ecotypes

Background

Improved performance in a given ecological niche can occur through local adaptation, phenotypic plasticity, or a combination of these mechanisms. Evaluating the relative importance of these two mechanisms is needed to better understand the cause of intra specific polymorphism. In this study, we reared populations of Lake Whitefish (Coregonus clupeaformis) representing the’normal’ (benthic form) and the ‘dwarf’ (derived limnetic form) ecotypes in two different conditions (control and swim-training) to test the relative importance of adaptation and acclimation in the differentiation of traits related to swimming capacity. The dwarf whitefish is a more active swimmer than the normal ecotype, and also has a higher capacity for aerobic energy production in the swimming musculature. We hypothesized that dwarf fish would show changes in morphological and physiological traits consistent with reductions in the energetic costs of swimming and maintenance metabolism.

Results

We found differences in traits predicted to decrease the costs of prolonged swimming and standard metabolic rate and allow for a more active lifestyle in dwarf whitefish. Dwarf whitefish evolved a more streamlined body shape, predicted to lead to a decreased drag, and a smaller brain, which may decrease their standard metabolic rate. Contrary to predictions, we also found evidence of acclimation in liver size and metabolic enzyme activities.

Conclusion

Results support the view that local adaptation has contributed to the genetically-based divergence of traits associated with swimming activity. Presence of post-zygotic barriers limiting gene flow between these ecotype pairs may have favoured repeated local adaptation to the limnetic niches.

Electronic supplementary material

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

Global transcriptional and miRNA insights into bases of heterosis in hybridization of Cyprinidae

Hybrid Megalobrama amblycephala × Culter alburnus represents a population newly formed by interspecific crossing between two different genera. Here we assessed the expression pattern of mRNA and small RNA in newly formed F1, F2 and their progenitors. Large amounts of nonadditively expressed protein-coding genes showed parental expression level dominance (ELD). Interestingly, the ELD pattern could inherit from F1 to F2, which guaranteed a stable appearance in progenies. The ELD-B genes were found to contribute to cell development, while the ELD-T genes were enriched in function of stress and adaptability. microRNAs (miRNA) also had similar expression patterns to genes. A high proportion of miRNAs showed nonadditive expression upon hybridization, and were found to target important genes with diverse roles potentially involved in stress adaption and development. Taken together, the gene and miRNA expression divergence contributes to heterosis in the newly formed hybrid, promising the successful existence of hybrid speciation.

Adaptation and acclimation of aerobic exercise physiology in Lake Whitefish ecotypes (Coregonus clupeaformis)

The physiological mechanisms underlying local adaptation in natural populations of animals, and whether the same mechanisms contribute to adaptation and acclimation, are largely unknown. Therefore, we tested for evolutionary divergence in aerobic exercise physiology in laboratory bred, size-matched crosses of ancestral, benthic, normal Lake Whitefish (Coregonus clupeaformis) and derived, limnetic, more actively swimming "dwarf" ecotypes. We acclimated fish to constant swimming (emulating limnetic foraging) and control conditions (emulating normal activity levels) to simultaneously study phenotypic plasticity. We found extensive divergence between ecotypes: dwarf fish generally had constitutively higher values of traits related to oxygen transport (ventricle size) and use by skeletal muscle (percent oxidative muscle, mitochondrial content), and also evolved differential plasticity of mitochondrial function (Complex I activity and flux through Complexes I-IV and IV). The effects of swim training were less pronounced than differences among ecotypes and the traits which had a significant training effect (ventricle protein content, ventricle malate dehydrogenase activity, and muscle Complex V activity) did not differ among ecotypes. Only one trait, ventricle mass, varied in a similar manner with acclimation and adaptation and followed a pattern consistent with genetic accommodation. Overall, the physiological and biochemical mechanisms underlying acclimation and adaptation to swimming activity in Lake Whitefish differ.

In absence of local adaptation, plasticity and spatially varying selection rule: a view from genomic reaction norms in a panmictic species (Anguilla rostrata)

Background

American eel (Anguilla rostrata) is one of the few species for which panmixia has been demonstrated at the scale of the entire species. As such, the development of long term local adaptation is impossible. However, both plasticity and spatially varying selection have been invoked in explaining how American eel may cope with an unusual broad scope of environmental conditions. Here, we address this question through transcriptomic analyses and genomic reaction norms of eels from two geographic origins reared in controlled environments.

Results

The null hypothesis of no difference in gene expression between eels from the two origins was rejected. Many unique transcripts and two out of seven gene clusters showed significant difference in expression, both at time of capture and after three months of common rearing. Differences in expression were observed at numerous genes representing many functional groups when comparing eels from a same origin reared under different salinity conditions. Plastic response to different rearing conditions varied among gene clusters with three clusters showing significant origin-environment interactions translating into differential genomic norms of reaction. Most genes and functional categories showing differences between origins were previously shown to be differentially expressed in a study comparing transcription profiles between adult European eels acclimated to different salinities.

Conclusions

These results emphasize that while plasticity in expression may be important, there is also a role for local genetic (and/or epigenetic) differences in explaining differences in gene expression between eels from different geographic origins. Such differences match those reported in genetically distinct populations in other fishes, both in terms of the proportion of genes that are differentially expressed and the diversity of biological functions involved. We thus propose that genetic differences between glass eels of different origins caused by spatially varying selection due to local environmental conditions translates into transcriptomic differences (including different genomic norms of reaction) which may in turn explain part of the phenotypic variance observed between different habitats colonized by eels.

Targeted sequence capture and resequencing implies a predominant role of regulatory regions in the divergence of a sympatric lake whitefish species pair (Coregonus clupeaformis)

Latest technological developments in evolutionary biology bring new challenges in documenting the intricate genetic architecture of species in the process of divergence. Sympatric populations of lake whitefish represent one of the key systems to investigate this issue. Despite the value of random genotype-by-sequencing methods and decreasing cost of sequencing technologies, it remains challenging to investigate variation in coding regions, especially in the case of recently duplicated genomes as in salmonids, as this greatly complicates whole genome resequencing. We thus designed a sequence capture array targeting 2773 annotated genes to document the nature and the extent of genomic divergence between sympatric dwarf and normal whitefish. Among the 2728 genes successfully captured, a total of 2182 coding and 10,415 noncoding putative single-nucleotide polymorphisms (SNPs) were identified after applying a first set of basic filters. A genome scan with a quality-refined selection of 2203 SNPs identified 267 outlier SNPs in 210 candidate genes located in genomic regions potentially involved in whitefish divergence and reproductive isolation. We found highly heterogeneous FST estimates among SNP loci. There was an overall low level of coding polymorphism, with a predominance of noncoding mutations among outliers. The heterogeneous patterns of divergence among loci confirm the porous nature of genomes during speciation with gene flow. Considering that few protein-coding mutations were identified as highly divergent, our results, along with previous transcriptomic studies, imply that changes in regulatory regions most likely had a greater role in the process of whitefish population divergence than protein-coding mutations. This study is the first to demonstrate the efficiency of large-scale targeted resequencing for a nonmodel species with such a large and unsequenced genome.

Gene coexpression networks reveal key drivers of phenotypic divergence in lake whitefish

A functional understanding of processes involved in adaptive divergence is one of the awaiting opportunities afforded by high-throughput transcriptomic technologies. Functional analysis of coexpressed genes has succeeded in the biomedical field in identifying key drivers of disease pathways. However, in ecology and evolutionary biology, functional interpretation of transcriptomic data is still limited. Here, we used Weighted Gene Co-Expression Network Analysis (WGCNA) to identify modules of coexpressed genes in muscle and brain tissue of a lake whitefish backcross progeny. Modules were connected to gradients of known adaptive traits involved in the ecological speciation process between benthic and limnetic ecotypes. Key drivers, that is, hub genes of functional modules related to reproduction, growth, and behavior were identified, and module preservation was assessed in natural populations. Using this approach, we identified modules of coexpressed genes involved in phenotypic divergence and their key drivers, and further identified a module part specifically rewired in the backcross progeny. Functional analysis of transcriptomic data can significantly contribute to the understanding of the mechanisms underlying ecological speciation. Our findings point to bone morphogenetic protein and calcium signaling as common pathways involved in coordinated evolution of trophic behavior, trophic morphology (gill rakers), and reproduction. Results also point to pathways implicating hemoglobins and constitutive stress response (HSP70) governing growth in lake whitefish.

The influence of parental effects on transcriptomic landscape during early development in brook charr (Salvelinus fontinalis, Mitchill)

Parental effects represent an important source of variation in offspring phenotypes. Depending on the specific mechanisms involved, parental effects may be caused to different degrees by either the maternal or the paternal parent, and these effects may in turn act at different stages of development. To detect parental effects acting on gene transcription regulation and length phenotype during ontogeny, the transcriptomic profiles of two reciprocal hybrids from Laval × Rupert and Laval × Domestic populations of brook charr were compared at hatching, yolk sac resorption and 15 weeks after exogenous feeding. Using a salmonid cDNA microarray, our results show that parental effects modulated gene expression among reciprocal hybrids only at the yolk sac resorption stage. In addition, Laval × Domestic and Laval × Rupert reciprocal hybrids differed in the magnitude of theses parental effects, with 199 and 630 differentially expressed transcripts, respectively. This corresponds to a maximum of 18.5% of the analyzed transcripts. These transcripts are functionally related to cell cycle, nucleic acid metabolism and intracellular protein traffic, which is consistent with observed differences associated with embryonic development and growth differences in other fish species. Our results thus illustrate how parental effects on patterns of gene transcription seem dependent on the genetic architecture of the parents. In addition, in absence of transcriptional differences, non-transcript deposits in the yolk sac could contribute to the observed length differences among the reciprocal hybrids before yolk sac resorption.

Stocking impacts the expression of candidate genes and physiological condition in introgressed brook charr (Salvelinus fontinalis) populations

Translocation of plants and animal populations between environments is one of the major forms of anthropogenic perturbation experienced by pristine populations, and consequently, human-mediated hybridization by stocking practices between wild and exogenous conspecifics is of increasing concern. In this study, we compared the expression of seven candidate genes involved in multifactorial traits and regulatory pathways for growth as a function of level of introgressive hybridization between wild and domestic brook charr to test the null hypothesis of no effect of introgression on wild fish. Our analyses revealed that the expression of two of the genes tested, cytochrome c oxidase VIIa and the growth hormone receptor isoform I, was positively correlated with the level of introgression. We also observed a positive relationship between the extent of introgression and physiological status quantified by the Fulton's condition index. The expression of other genes was influenced by other variables, including year of sampling (reflecting different thermal conditions), sampling method and lake of origin. This is the first demonstration in nature that introgression from stocked populations has an impact on the expression of genes playing a role in important biological functions that may be related with fitness in wild introgressed populations.

Exploring local immunological adaptation of two stickleback ecotypes by experimental infection and transcriptome-wide digital gene expression analysis

Understanding the extent of local adaptation in natural populations and the mechanisms that allow individuals to adapt to their native environment is a major avenue in molecular ecology research. Evidence for the frequent occurrence of diverging ecotypes in species that inhabit multiple ecological habitats is accumulating, but experimental approaches to understanding the biological pathways as well as the underlying genetic mechanisms are still rare. Parasites are invoked as one of the major selective forces driving evolution and are themselves dependent on the ecological conditions in a given habitat. Immunological adaptation to local parasite communities is therefore expected to be a key component of local adaptation in natural populations. Here, we use next-generation sequencing technology to compare the transcriptome-wide response of experimentally infected three-spined sticklebacks from a lake and a river population, which are known to evolve under selection by distinct parasite communities. By comparing overall gene expression levels as well as the activation of functional pathways in response to parasite exposure, we identified potential differences between the two stickleback populations at several levels. Our results suggest locally adapted patterns of gene regulation in response to parasite exposure, which may reflect different local optima in the trade-off between the benefits and the disadvantages of mounting an immune response because of quantitative differences of the local parasite communities.

Role of waterlogging-responsive genes in shaping interspecific differentiation between two sympatric oak species

Pedunculate (Quercus robur L.) and sessile oak (Quercus petreae Matt. Liebl.) are closely related species with a widely sympatric distribution in Europe. These two oak species are also known to display different ecological features, particularly related to their adaptation to soil waterlogging. Pedunculate oak grows in humid areas and can withstand high moisture content of the soil, whereas sessile oak requires drier soil with better drainage. The main goal of this study was to explore the role of gene expression contributing to differences in terms of waterlogging tolerance between these two species. We implemented a series of experiments aimed at evaluating whether differentially expressed genes between species are associated with their ecological preferences and underlie adaptive genetic divergence. Rooted cuttings of both species were grown in hydroponic conditions and subjected to gradual root hypoxia. White roots were sampled after 6, 12, 24 and 48 h. Real-time polymerase chain reaction (qPCR) was first used to monitor the expression of 10 known waterlogging-responsive genes, to identify discriminating sampling time points along the kinetics of hypoxia. Secondly, four subtractive suppressive hybridization libraries (sessile vs. pedunculate, pedunculate vs. sessile for early and late responses) were generated to isolate differentially expressed genes between species. A total of 2160 high-quality expressed sequence tags were obtained and annotated, and a subset of 45 genes were selected for qPCR analysis in a second independent factorial experimental design applying two stress durations per two species. Significant differences of gene expression between pedunculate and sessile oaks were detected, suggesting species-specific molecular strategies to respond to hypoxia. This study revealed that the ability of pedunculate oak to maintain glycolysis and fermentation under hypoxic conditions may help explain its tolerance to waterlogging.

Mapping gene expression in two Xenopus species: evolutionary constraints and developmental flexibility

Changes in gene expression are thought to be important for morphological evolution, though little is known about the nature or magnitude of the differences. Here, we examine Xenopus laevis and Xenopus tropicalis, two amphibians with very similar development, and ask how their transcriptomes compare. Despite separation for ~30-90 million years, there is strong conservation in gene expression in the vast majority of the expressed orthologs. Significant changes occur in the level of gene expression but changes in the timing of expression (heterochrony) were much less common. Differences in level were concentrated in the earliest embryonic stages. Changes in timing were prominently found in pathways that respond to selective features of the environment. We propose that different evolutionary rates across developmental stages may be explained by the stabilization of cell fate determination in the later stages.

Transcriptome profiling of embryonic development rate in rainbow trout advanced backcross introgression lines

In rainbow trout (Oncorhynchus mykiss) and other fishes, embryonic development rate is an ecologically and evolutionarily important trait that is closely associated with survival and physiological performance later in life. To identify genes differentially regulated in fast and slow-developing embryos of rainbow trout, we examined gene expression across developmental time points in rainbow trout embryos possessing alleles linked to a major quantitative trait loci (QTL) for fast versus slow embryonic development rate. Whole genome expression microarray analyses were conducted using embryos from a fourth generation backcross family, whereby each backcross generation involved the introgression of the fast-developing alleles for a major development rate QTL into a slow-developing clonal line of rainbow trout. Embryos were collected at 15, 19, and 28 days post-fertilization; sex and QTL genotype were determined using molecular markers, and cDNA from 48 embryos were used for microarray analysis. A total of 183 features were identified with significant differences between embryonic development rate genotypes. Genes associated with cell cycle growth, muscle contraction and protein synthesis were expressed significantly higher in embryos with the fast-developing allele (Clearwater) than those with the slow-developing allele (Oregon State University), which may associate with fast growth and early body mass construction in embryo development. Across time points, individuals with the fast-developing QTL allele appeared to have earlier onset of these developmental processes when compared to individuals with the slow development alleles, even as early as 15 days post-fertilization. Differentially expressed candidate genes chosen for linkage mapping were localized primarily to regions outside of the major embryonic development rate QTL, with the exception of a single gene (very low-density lipoprotein receptor precursor).

A phylogenetically based transcriptome age index mirrors ontogenetic divergence patterns

Parallels between phylogeny and ontogeny have been discussed for almost two centuries, and a number of theories have been proposed to explain such patterns. Especially elusive is the phylotypic stage, a phase during development where species within a phylum are particularly similar to each other. Although this has formerly been interpreted as a recapitulation of phylogeny, it is now thought to reflect an ontogenetic progression phase, where strong constraints on developmental regulation and gene interactions exist. Several studies have shown that genes expressed during this stage evolve at a slower rate, but it has so far not been possible to derive an unequivocal molecular signature associated with this stage. Here we use a combination of phylostratigraphy and stage-specific gene expression data to generate a cumulative index that reflects the evolutionary age of the transcriptome at given ontogenetic stages. Using zebrafish ontogeny and adult development as a model, we find that the phylotypic stage does indeed express the oldest transcriptome set and that younger sets are expressed during early and late development, thus faithfully mirroring the hourglass model of morphological divergence. Reproductively active animals show the youngest transcriptome, with major differences between males and females. Notably, ageing animals express increasingly older genes. Comparisons with similar data sets from flies and nematodes show that this pattern occurs across phyla. Our results indicate that an old transcriptome marks the phylotypic phase and that phylogenetic differences at other ontogenetic stages correlate with the expression of newly evolved genes.

Extent and scale of local adaptation in salmonid fishes: review and meta-analysis

What is the extent and scale of local adaptation (LA)? How quickly does LA arise? And what is its underlying molecular basis? Our review and meta-analysis on salmonid fishes estimates the frequency of LA to be ∼55-70%, with local populations having a 1.2 times average fitness advantage relative to foreign populations or to their performance in new environments. Salmonid LA is evident at a variety of spatial scales (for example, few km to>1000 km) and can manifest itself quickly (6-30 generations). As the geographic scale between populations increases, LA is generally more frequent and stronger. Yet the extent of LA in salmonids does not appear to differ from that in other assessed taxa. Moreover, the frequency with which foreign salmonid populations outperform local populations (∼23-35%) suggests that drift, gene flow and plasticity often limit or mediate LA. The relatively few studies based on candidate gene and genomewide analyses have identified footprints of selection at both small and large geographical scales, likely reflecting the specific functional properties of loci and the associated selection regimes (for example, local niche partitioning, pathogens, parasites, photoperiodicity and seasonal timing). The molecular basis of LA in salmonids is still largely unknown, but differential expression at the same few genes is implicated in the convergent evolution of certain phenotypes. Collectively, future research will benefit from an integration of classical and molecular approaches to understand: (i) species differences and how they originate, (ii) variation in adaptation across scales, life stages, population sizes and environmental gradients, and (iii) evolutionary responses to human activities.

Transcriptome-wide signature of hybrid breakdown associated with intrinsic reproductive isolation in lake whitefish species pairs (Coregonus spp. Salmonidae)

Genetic analyses of speciation have focused nearly exclusively on retrospective analyses of reproductive isolation between highly divergent species. Yet, a full understanding of the speciation process must encompass analysis of the consequences of genomic divergence in young lineages still capable of exchanging genes under natural conditions. The accumulation of conditionally neutral genetic variation may lead to the evolution of divergent gene networks. In a hybrid background, such mutations may no longer compensate one another, resulting in the appearance of selectively disadvantageous traits, including disruption of gene expression regulation. Here, we documented genome-wide patterns of gene expression divergence between young lineages of normal and dwarf lake whitefish and their backcross hybrids for which strong, yet incomplete post-zygotic isolation barriers exist. A significant proportion (33%) of backcross hybrids showed developmental abnormalities not seen in parental forms and eventually leading to death. Although the transcriptome of parental forms was nearly identical during embryonic development, suggesting a role for stabilizing selection, all hybrids displayed strongly divergent patterns of gene expression. By comparing healthy, surviving hybrids against moribund ones showing abnormal development, we observed that over 2000 genes were misregulated in these abnormal embryos. In particular, misregulation was significantly biased toward essential developmental genes, which were severely underexpressed. Furthermore, genes previously documented to be highly transgressive (exaggerated inter-individual variance) were almost invariably underexpressed in hybrids. Our results thus clearly showed a transcriptome-wide signature of hybrid breakdown in young, incipient species and demonstrated a persuasive link between misexpression of essential developmental genes and post-zygotic isolation.

The transcriptomics of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis spp., Salmonidae) divergence as revealed by next-generation sequencing

Gene expression divergence is one of the mechanisms thought to be involved in the emergence of incipient species. Next-generation sequencing has become an extremely valuable tool for the study of this process by allowing whole transcriptome sequencing, or RNA-Seq. We have conducted a 454 GS-FLX pyrosequencing experiment to refine our understanding of adaptive divergence between dwarf and normal lake whitefish species (Coregonus clupeaformis spp.). The objectives were to: (i) investigate transcriptomic divergence as measured by liver RNA-Seq; (ii) test the correlation between divergence in expression and sequence polymorphism; and (iii) investigate the extent of allelic imbalance. We also compared the results of RNA-seq with those of a previous microarray study performed on the same fish. Following de novo assembly, results showed that normal whitefish overexpressed more contigs associated with protein synthesis while dwarf fish overexpressed more contigs related to energy metabolism, immunity and DNA replication and repair. Moreover, 63 SNPs showed significant allelic imbalance, and this phenomenon prevailed in the recently diverged dwarf whitefish. Results also showed an absence of correlation between gene expression divergence as measured by RNA-Seq and either polymorphism rate or sequence divergence between normal and dwarf whitefish. This study reiterates an important role for gene expression divergence, and provides evidence for allele-specific expression divergence as well as evolutionary decoupling of regulatory and coding sequences in the adaptive divergence of normal and dwarf whitefish. It also demonstrates how next-generation sequencing can lead to a more comprehensive understanding of transcriptomic divergence in a young species pair.

On the origin of species: insights from the ecological genomics of lake whitefish

In contrast to the large amount of ecological information supporting the role of natural selection as a main cause of population divergence and speciation, an understanding of the genomic basis underlying those processes is in its infancy. In this paper, we review the main findings of a long-term research programme that we have been conducting on the ecological genomics of sympatric forms of whitefish (Coregonus spp.) engaged in the process of speciation. We present this system as an example of how applying a combination of approaches under the conceptual framework of the theory of adaptive radiation has yielded substantial insight into evolutionary processes in a non-model species. We also discuss how the joint use of recent biotechnological developments will provide a powerful means to address issues raised by observations made to date. Namely, we present data illustrating the potential offered by combining next generation sequencing technologies with other genomic approaches to reveal the genomic bases of adaptive divergence and reproductive isolation. Given increasing access to these new genomic tools, we argue that non-model species studied in their ecological context such as whitefish will play an increasingly important role in generalizing knowledge of speciation.

Interspecific resources: a major tool for quantitative trait locus cloning and speciation research

Positional cloning of the quantitative trait locus (QTL) still encounters numerous difficulties, which explains why thousands of QTL have been mapped, while only a few have been identified at the molecular level. Here, we focus on a specific mapping tool that exists in plant and animal model species: interspecific recombinant congenic strains (IRCSs) or interspecific nearly isogenic lines (NILs). Such panels exhibit a much higher sequence diversity than intraspecific sets, thus enhancing the contrasts between phenotypes. In animals, it allows statistical significance to be reached even when using a limited number of individuals. Therefore, we argue that interspecific resources may constitute a major genetic tool for positional cloning and for understanding some bases of speciation mechanisms.

Fast transcriptional responses to domestication in the brook charr Salvelinus fontinalis

Domestication has been practiced for centuries yet directed toward relatively few terrestrial crops and animals. While phenotypic and quantitative genetic changes associated with domestication have been amply documented, little is known about the molecular changes underlying the phenotypic evolution during the process. Here, we have investigated the brook charr (Salvelinus fontinalis) responses to artificial selection by means of transcriptional analysis of approximately 32,000 cDNA features performed in both selected and control populations reared under identical environmental conditions during four generations. Our results indicate that selective breeding led to significant changes in the transcription of genes at the juvenile stage, where we observed 4.16% (156/3750) of differentially expressed genes between the two lines. No significant genes were revealed at the earlier life stage. Moreover, when comparing our results to those of previous studies on Atlantic salmon that compared lines that were selected for five to seven generations for similar traits (e.g., growth), genes with similar biological functions were found to be under selection in both studies. These observations indicate that (1) four generations of selection caused substantial changes in regulation of gene transcription between selected and control populations and (2) selective breeding for improving the same phenotypic traits (e.g., rapid growth) in brook charr and Atlantic salmon tended to select for the same changes in transcription profiles as the expression of a small and similar set of genes was affected by selection.

Linking genotypes to phenotypes and fitness: how mechanistic biology can inform molecular ecology

The accessibility of new genomic resources, high-throughput molecular technologies and analytical approaches such as genome scans have made finding genes contributing to fitness variation in natural populations an increasingly feasible task. Once candidate genes are identified, we argue that it is necessary to take a mechanistic approach and work up through the levels of biological organization to fully understand the impacts of genetic variation at these candidate genes. We demonstrate how this approach provides testable hypotheses about the causal links among levels of biological organization, and assists in designing relevant experiments to test the effects of genetic variation on phenotype, whole-organism performance capabilities and fitness. We review some of the research programs that have incorporated mechanistic approaches when examining naturally occurring genetic and phenotypic variation and use these examples to highlight the value of developing a comprehensive understanding of the relationship between genotype and fitness. We give suggestions to guide future research aimed at uncovering and understanding the genetic basis of adaptation and argue that further integration of mechanistic approaches will help molecular ecologists better understand the evolution of natural populations.

Major disruption of gene expression in hybrids between young sympatric anadromous and resident populations of brook charr (Salvelinus fontinalis Mitchill)

Genome-wide analyses of the transcriptome have suggested that male-biased genes are the first targets of genomic incompatibilities (g.i.) in inter-specific hybrids. However, those studies have almost invariably focused on Drosophila species that diverged at least 0.9 Ma, and with sterile male hybrids. Here, we use microarrays to analyse patterns of gene expression in very closely related (divergence <12,000 years), sympatric, but ecologically divergent anadromous and resident populations of brook charr (Salvelinus fontinalis) and their F(1) hybrids. Our results show a dramatic breakdown of gene expression patterns in hybrids compared with their parental relatives. Several disrupted genes are related to energetic metabolism, immune response, osmoregulation and protection against oxidative stress, and none has sex-biased functions. Besides, pure individuals show no expression differences at most of the genes disrupted in hybrids, which may suggest the operation of some form of stabilizing selection. Taken together, these results both confirm the idea that perturbations of regulatory networks represent a significant source of g.i. and support the suggestion that developmental pathways can diverge through time without any manifest change in the phenotypic outcome. While the role of other evolutionary forces (e.g. genetic drift) cannot be ruled out, this study suggests that ecological selective processes may provide the initial driving force behind disruption of gene expression in inter-specific hybrids.

Gene expression divergence and hybrid misexpression between lake whitefish species pairs (Coregonus spp. Salmonidae)

Genomewide analyses of the transcriptome have confirmed that gene misexpression may underlie reproductive isolation mechanisms in interspecific hybrids. Here, using a 16,006 features cDNA microarray, we compared and contrasted gene expression divergence at two ontogenetic stages in incipient species of normal and dwarf whitefish (Coregonus clupeaformis) with that of first generation (normal x dwarf) and second-generation hybrid crosses (backcross: [normal x dwarf] x normal]. Our goal was to identify the main mode of action responsible for gene transcription and to discover key genes misexpressed in hybrids. Very few transcripts (five of 4,950 expressed) differed in mean expression level between parentals and hybrids at the embryonic stage, in contrast to 16-week-old juvenile fish for which 617 out of 5,359 transcripts differed significantly. We also found evidence for more misexpression in backcross hybrids whereby nonadditivity explained a larger fraction of hybrid inheritance patterns in backcross (54%) compared with F1-hybrids (9%). Gene expression in hybrids was more variable than in pure crosses and transgressive patterns of expression were ubiquitous in hybrids. In backcross embryos in particular, the expression of three key developmental genes involved in protein folding and mRNA translation was severely disrupted. Accordingly, gene misexpression in hybrids adds to other factors previously identified as contributing to the reproductive isolation of incipient species of lake whitefish.