Gene co-expression patterns in Atlantic salmon adipose tissue provide a molecular link among seasonal changes, energy balance and age at maturity

Sexual maturation in many fishes requires a major physiological change that involves a rapid transition between energy storage and usage. In Atlantic salmon, this transition for the initiation of maturation is tightly controlled by seasonality and requires a high-energy status. Lipid metabolism is at the heart of this transition since lipids are the main energy storing molecules. The balance between lipogenesis (lipid accumulation) and lipolysis (lipid use) determines energy status transitions. A genomic region containing a transcription co-factor of the Hippo pathway, vgll3, is the main determinant of maturation timing in Atlantic salmon. Interestingly, vgll3 acts as an inhibitor of adipogenesis in mice and its genotypes are potentially associated with seasonal heterochrony in lipid storage and usage in juvenile Atlantic salmon. Here, we explored changes in expression of more than 300 genes directly involved in the processes of adipogenesis, lipogenesis and lipolysis, as well as the Hippo pathway in the adipose tissue of immature and mature Atlantic salmon with distinct vgll3 genotypes. We found molecular evidence consistent with a scenario in which immature males with different vgll3 genotypes exhibit contrasting seasonal dynamics in their lipid profiles. We also identified components of the Hippo signalling pathway as potential major drivers of vgll3 genotype-specific differences in adipose tissue gene expression. This study demonstrates the importance of adipose gene expression patterns for directly linking environmental changes with energy balance and age at maturity through genetic factors bridging lipid metabolism, seasonality and sexual maturation.

The role of alternative splicing in adaptation and evolution

Regulation of gene expression plays a central role in adaptive divergence and evolution. Although the role of gene regulation in microevolutionary processes is gaining wide acceptance, most studies have only investigated the evolution of transcript levels, ignoring the potentially significant role of transcript structures. We argue that variation in alternative splicing plays an important and widely unexplored role in adaptation (e.g., by increasing transcriptome and/or proteome diversity, or buffering potentially deleterious genetic variation). New studies increasingly highlight the potential for independent evolution in alternative splicing and transcript level, providing alternative paths for selection to act upon. We propose that alternative splicing and transcript levels can provide contrasting, nonredundant mechanisms of equal importance for adaptive diversification of gene function and regulation.

Genetic Drift Dominates Genome-Wide Regulatory Evolution Following an Ancient Whole-Genome Duplication in Atlantic Salmon

Whole-genome duplications (WGD) have been considered as springboards that potentiate lineage diversification through increasing functional redundancy. Divergence in gene regulatory elements is a central mechanism for evolutionary diversification, yet the patterns and processes governing regulatory divergence following events that lead to massive functional redundancy, such as WGD, remain largely unknown. We studied the patterns of divergence and strength of natural selection on regulatory elements in the Atlantic salmon (Salmo salar) genome, which has undergone WGD 100-80 Ma. Using ChIPmentation, we first show that H3K27ac, a histone modification typical to enhancers and promoters, is associated with genic regions, tissue-specific transcription factor binding motifs, and with gene transcription levels in immature testes. Divergence in transcription between duplicated genes from WGD (ohnologs) correlated with difference in the number of proximal regulatory elements, but not with promoter elements, suggesting that functional divergence between ohnologs after WGD is mainly driven by enhancers. By comparing H3K27ac regions between duplicated genome blocks, we further show that a longer polyploid state post-WGD has constrained regulatory divergence. Patterns of genetic diversity across natural populations inferred from resequencing indicate that recent evolutionary pressures on H3K27ac regions are dominated by largely neutral evolution. In sum, our results suggest that post-WGD functional redundancy in regulatory elements continues to have an impact on the evolution of the salmon genome, promoting largely neutral evolution of regulatory elements despite their association with transcription levels. These results highlight a case where genome-wide regulatory evolution following an ancient WGD is dominated by genetic drift.

The Atlantic salmon whole blood transcriptome and how it relates to major locus maturation genotypes and other tissues

The Atlantic salmon (Salmo salar) is important to many ecosystems and local economies and has therefore become the focus of a broad range of research questions that have benefited from the availability of high-quality genomic resources. Albeit gene expression studies have been extensive for this species, the transcriptome information for Atlantic salmon whole blood has been lacking. A transcriptome of Atlantic salmon blood would be a valuable resource for future studies, especially those wishing to take non-lethal samples. Here, we report a whole blood transcriptome for Atlantic salmon constructed from twelve 8-month old salmon parr using RNA-seq. We identify transcriptomic proxies for the genotype at the major maturation timing locus vestigial-like 3 (vgll3). Differentially expressed genes between the early and late maturing genotypes showed overrepresented Gene Ontology (GO) terms with the strongest result linked to 13 ribosomal subunit genes. To assess how the whole blood gene expression profile relates to other tissues, we compare the blood transcriptome to the reference transcriptome of fourteen other tissue types using both a common PCA method and a novel method. The novel method compares transcriptomes when gene expression is visualised as a layer using thin-plate spline smoothers. Both methods found similar patterns with the blood transcriptome being quite unique compared to the transcription profiles of other tissues.

mRNA Extraction from Gill Tissue for RNA-sequencing

Adaptation is thought to proceed in part through spatial and temporal changes in gene expression. Fish species such as the threespine stickleback are powerful vertebrate models to study the genetic architecture of adaptive changes in gene expression since divergent adaptation to different environments is common, they are abundant and easy to study in the wild and lab, and have well-established genetic and genomic resources. Fish gills, due to their respiratory and osmoregulatory roles, show many physiological adaptations to local water chemistry, including differences in gene expression. However, obtaining high-quality RNA using popular column-based extraction methods can be challenging from small tissue samples high in cartilage and bone such as fish gills. Here, we describe a bead-based mRNA extraction and transcriptome RNA-seq protocol that does not use purification columns. The protocol can be readily scaled according to sample size for the purposes of diverse gene expression experiments using animal or plant tissue.

Predominance of cis-regulatory changes in parallel expression divergence of sticklebacks

Regulation of gene expression is thought to play a major role in adaptation, but the relative importance of cis- and trans- regulatory mechanisms in the early stages of adaptive divergence is unclear. Using RNAseq of threespine stickleback fish gill tissue from four independent marine-freshwater ecotype pairs and their F1 hybrids, we show that cis-acting (allele-specific) regulation consistently predominates gene expression divergence. Genes showing parallel marine-freshwater expression divergence are found near to adaptive genomic regions, show signatures of natural selection around their transcription start sites and are enriched for cis-regulatory control. For genes with parallel increased expression among freshwater fish, the quantitative degree of cis- and trans-regulation is also highly correlated across populations, suggesting a shared genetic basis. Compared to other forms of regulation, cis-regulation tends to show greater additivity and stability across different genetic and environmental contexts, making it a fertile substrate for the early stages of adaptive evolution.

Dissection of expression-quantitative trait locus and allele specificity using a haploid/diploid plant system - insights into compensatory evolution of transcriptional regulation within populations

Regulation of gene expression plays a central role in translating genotypic variation into phenotypic variation. Dissection of the genetic basis of expression variation is key to understanding how expression regulation evolves. Such analyses remain challenging in contexts where organisms are outbreeding, highly heterozygous and long-lived such as in the case of conifer trees. We developed an RNA sequencing (RNA-seq)-based approach for both expression-quantitative trait locus (eQTL) mapping and the detection of cis-acting (allele-specific) vs trans-acting (non-allele-specific) eQTLs. This method can be potentially applied to many conifers. We used haploid and diploid meiotic seed tissues of a single self-fertilized white spruce (Picea glauca) individual to dissect eQTLs according to linkage and allele specificity. The genetic architecture of local eQTLs linked to the expressed genes was particularly complex, consisting of cis-acting, trans-acting and, surprisingly, compensatory cis-trans effects. These compensatory effects influence expression in opposite directions and are neutral when combined in homozygotes. Nearly half of local eQTLs were under compensation, indicating that close linkage between compensatory cis-trans factors is common in spruce. Compensated genes were overrepresented in developmental and cell organization functions. Our haploid-diploid eQTL analysis in spruce revealed that compensatory cis-trans eQTLs segregate within populations and evolve in close genetic linkage.

Conifer genomics and adaptation: at the crossroads of genetic diversity and genome function

Conifers have been understudied at the genomic level despite their worldwide ecological and economic importance but the situation is rapidly changing with the development of next generation sequencing (NGS) technologies. With NGS, genomics research has simultaneously gained in speed, magnitude and scope. In just a few years, genomes of 20-24 gigabases have been sequenced for several conifers, with several others expected in the near future. Biological insights have resulted from recent sequencing initiatives as well as genetic mapping, gene expression profiling and gene discovery research over nearly two decades. We review the knowledge arising from conifer genomics research emphasizing genome evolution and the genomic basis of adaptation, and outline emerging questions and knowledge gaps. We discuss future directions in three areas with potential inputs from NGS technologies: the evolutionary impacts of adaptation in conifers based on the adaptation-by-speciation model; the contributions of genetic variability of gene expression in adaptation; and the development of a broader understanding of genetic diversity and its impacts on genome function. These research directions promise to sustain research aimed at addressing the emerging challenges of adaptation that face conifer trees.

Transcriptome profiling in conifers and the PiceaGenExpress database show patterns of diversification within gene families and interspecific conservation in vascular gene expression

Background

Conifers have very large genomes (13 to 30 Gigabases) that are mostly uncharacterized although extensive cDNA resources have recently become available. This report presents a global overview of transcriptome variation in a conifer tree and documents conservation and diversity of gene expression patterns among major vegetative tissues.

Results

An oligonucleotide microarray was developed from Picea glauca and P. sitchensis cDNA datasets. It represents 23,853 unique genes and was shown to be suitable for transcriptome profiling in several species. A comparison of secondary xylem and phelloderm tissues showed that preferential expression in these vascular tissues was highly conserved among Picea spp. RNA-Sequencing strongly confirmed tissue preferential expression and provided a robust validation of the microarray design. A small database of transcription profiles called PiceaGenExpress was developed from over 150 hybridizations spanning eight major tissue types. In total, transcripts were detected for 92% of the genes on the microarray, in at least one tissue. Non-annotated genes were predominantly expressed at low levels in fewer tissues than genes of known or predicted function. Diversity of expression within gene families may be rapidly assessed from PiceaGenExpress. In conifer trees, dehydrins and late embryogenesis abundant (LEA) osmotic regulation proteins occur in large gene families compared to angiosperms. Strong contrasts and low diversity was observed in the dehydrin family, while diverse patterns suggested a greater degree of diversification among LEAs.

Conclusion

Together, the oligonucleotide microarray and the PiceaGenExpress database represent the first resource of this kind for gymnosperm plants. The spruce transcriptome analysis reported here is expected to accelerate genetic studies in the large and important group comprised of conifer trees.