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 mitochondrial genome of the bioluminescent fish Malacosteusniger Ayres, 1848 (Stomiidae, Actinopterygii) is large and complex, and contains an inverted-repeat structure.

We determined the complete mitogenome sequence of the bioluminescent fish Malacosteusniger using long-read sequencing technologies. The 21,263 bp mitogenome features a complex structure with two copies of a 1198-bp inverted-repeat and a region of 2616-bp containing alternating copies of 16 and 26 bp repeat elements. Whole mitogenome phylogenies inferred from both nucleotide and amino-acid datasets place M.niger among Melanostomiinae. The need for additional complete mitogenome sequences from the subfamily Malacosteinae is discussed.

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.

NanoGBS: A Miniaturized Procedure for GBS Library Preparation.

High-throughput reduced-representation sequencing (RRS)-based genotyping methods, such as genotyping-by-sequencing (GBS), have provided attractive genotyping solutions in numerous species. Here, we present NanoGBS, a miniaturized and eco-friendly method for GBS library construction. Using acoustic droplet ejection (ADE) technology, NanoGBS libraries were constructed in tenfold smaller volumes compared to standard methods (StdGBS) and leading to a reduced use of plastics of up to 90%. A high-quality DNA library and SNP catalogue were obtained with extensive overlap (96%) in SNP loci and 100% agreement in genotype calls compared to the StdGBS dataset with a high level of accuracy (98.5%). A highly multiplexed pool of GBS libraries (768-plex) was sequenced on a single Ion Proton PI chip and yielded enough SNPs (~4K SNPs; 1.5 SNP per cM, on average) for many high-volume applications. Combining NanoGBS library preparation and increased multiplexing can dramatically reduce (72%) genotyping cost per sample. We believe that this approach will greatly facilitate the adoption of marker applications where extremely high throughputs are required and cost is still currently limiting.

The phenomics and expression quantitative trait locus mapping of brain transcriptomes regulating adaptive divergence in lake whitefish species pairs (Coregonus sp.).

We used microarrays and a previously established linkage map to localize the genetic determinants of brain gene expression for a backcross family of lake whitefish species pairs (Coregonus sp.). Our goals were to elucidate the genomic distribution and sex specificity of brain expression QTL (eQTL) and to determine the extent to which genes controlling transcriptional variation may underlie adaptive divergence in the recently evolved dwarf (limnetic) and normal (benthic) whitefish. We observed a sex bias in transcriptional genetic architecture, with more eQTL observed in males, as well as divergence in genome location of eQTL between the sexes. Hotspots of nonrandom aggregations of up to 32 eQTL in one location were observed. We identified candidate genes for species pair divergence involved with energetic metabolism, protein synthesis, and neural development on the basis of colocalization of eQTL for these genes with eight previously identified adaptive phenotypic QTL and four previously identified outlier loci from a genome scan in natural populations. Eighty-eight percent of eQTL-phenotypic QTL colocalization involved growth rate and condition factor QTL, two traits central to adaptive divergence between whitefish species pairs. Hotspots colocalized with phenotypic QTL in several cases, revealing possible locations where master regulatory genes, such as a zinc-finger protein in one case, control gene expression directly related to adaptive phenotypic divergence. We observed little evidence of colocalization of brain eQTL with behavioral QTL, which provides insight into the genes identified by behavioral QTL studies. These results extend to the transcriptome level previous work illustrating that selection has shaped recent parallel divergence between dwarf and normal lake whitefish species pairs and that metabolic, more than morphological, differences appear to play a key role in this divergence.

Transcriptional responses to environmental metal exposure in wild yellow perch (Perca flavescens) collected in lakes with differing environmental metal concentrations (Cd, Cu, Ni).

To investigate the mechanisms involved in metal stress in wild fish, yellow perch (Perca flavescens) were collected in eight lakes of the Rouyn-Noranda and Sudbury regions (Canada). Due to mining and smelting activities, these two regions indeed present a broad contamination gradient in metal concentrations (Cd, Cu, Zn and Ni; water, sediment and prey) and offer a unique research opportunity to investigate relationships between metal bioaccumulation and resulting deleterious effects in indigenous biota chronically exposed to metal mixtures. The expression level of genes encoding for proteins involved in metal detoxification (metallothioneins, mts), protein protection (heat shock protein-70, hsp-70), growth (insulin-like growth factor-1, igf-1), aerobic energy metabolism (cytochrome c oxydase, cco-1) and protection against oxidative stress (Cu/Zn superoxide dismutase, sod-1) were assessed in fish liver and muscle in association with protein and enzymatic assays for cytochrome c oxidase (CCO). Bioaccumulation of both Cd and Cu increased in response to higher ambient metal concentrations, but the two metals clearly have different modes of action. For Cd, changes in gene expression levels were more marked in the liver than in the dorsal muscle, whereas for Cu the opposite trend was observed. Hepatic Cd accumulation was linked to decreased cco-1 and sod-1 gene expression, whereas Cu accumulation was associated with a decrease in CCO enzymatic activity and an increase in total protein concentration and in cco-1, mts and hsp-70 gene expression levels. For Ni, no significant correlations were observed at the transcriptional level, but increasing hepatic Ni concentrations were significantly and positively correlated with protein concentrations and CCO activity. By coupling gene expression to biochemical and physiological endpoints, this work provides new insights into the mechanisms involved in metal stress and the adaptive response of fish chronically exposed to metal mixtures.

Genotyping-by-Sequencing on the Ion Torrent Platform in Barley.

The characterization of genetic polymorphism is a crucial step in both genetic studies and breeding programs. Genotyping-by-sequencing (GBS) constitutes one of the most attractive approaches for this purpose, especially in a genome as large as that of barley. The genome sequencing project undertaken by the International Barley Sequencing Consortium (IBSC) has produced a structured reference genome for the cultivar Morex [1] that can serve as an excellent resource for the analysis of GBS data. The genome assembly for this species [2] is thought to adequately capture the gene-rich portion of the genome (~80% of the entire genome). In this chapter, we describe the entire GBS process, from library preparation to the analysis of read data to produce a high-quality catalog of single nucleotide polymorphism (SNP) markers using the barley reference genome.

Differential timing of gene expression regulation between leptocephali of the two Anguilla eel species in the Sargasso Sea.

The unique life-history characteristics of North Atlantic catadromous eels have long intrigued evolutionary biologists, especially with respect to mechanisms that could explain their persistence as two ecologically very similar but reproductively and geographically distinct species. Differential developmental schedules during young larval stages have commonly been hypothesized to represent such a key mechanism. We performed a comparative analysis of gene expression by means of microarray experiments with American and European eel leptocephali collected in the Sargasso Sea in order to test the alternative hypotheses of (1) differential timing of gene expression regulation during early development versus (2) species-specific differences in expression of particular genes. Our results provide much stronger support for the former hypothesis since no gene showed consistent significant differences in expression levels between the two species. In contrast, 146 genes showed differential timings of expression between species, although the observed expression level differences between the species were generally small. Consequently, species-specific gene expression regulation seems to play a minor role in species differentiation. Overall, these results show that the basis of the early developmental divergence between the American and European eel is probably influenced by differences in the timing of gene expression regulation for genes involved in a large array of biological functions.

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.

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.