Background Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap, we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80–100 million years ago. Results We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a gene dose reduction while avoiding accumulation of “toxic mutations”. Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one likely source of novel binding sites driving tissue-specific gains in expression. Conclusions Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.

The long-term evolutionary impacts of whole-genome duplication (WGD) are strongly influenced by the ensuing rediploidization process. Following autopolyploidization, rediploidization involves a transition from tetraploid to diploid meiotic pairing, allowing duplicated genes (ohnologs) to diverge genetically and functionally. Our understanding of autopolyploid rediploidization has been informed by a WGD event ancestral to salmonid fishes, where large genomic regions are characterized by temporally delayed rediploidization, allowing lineage-specific ohnolog sequence divergence in the major salmonid clades. Here, we investigate the long-term outcomes of autopolyploid rediploidization at genomewide resolution, exploiting a recent "explosion" of salmonid genome assemblies, including a new genome sequence for the huchen (Hucho hucho). We developed a genome alignment approach to capture duplicated regions across multiple species, allowing us to create 121,864 phylogenetic trees describing genome-wide ohnolog divergence across salmonid evolution. Using molecular clock analysis, we show that 61% of the ancestral salmonid genome experienced an initial "wave" of rediploidization in the late Cretaceous (85-106 Ma). This was followed by a period of relative genomic stasis lasting 17-39 My, where much of the genome remained tetraploid. A second rediploidization wave began in the early Eocene and proceeded alongside species diversification, generating predictable patterns of lineage-specific ohnolog divergence, scaling in complexity with the number of speciation events. Using gene set enrichment, gene expression, and codon-based selection analyses, we provide insights into potential functional outcomes of delayed rediploidization. This study enhances our understanding of delayed autopolyploid rediploidization and has broad implications for future studies of WGD events.

The grass subfamily Pooideae dominates the grass floras in cold temperate regions and has evolved complex physiological adaptations to cope with extreme environmental conditions like frost, winter, and seasonality. One such adaptation is cold acclimation, wherein plants increase their frost tolerance in response to gradually falling temperatures and shorter days in the autumn. However, understanding how complex traits like cold acclimation evolve remains a major challenge in evolutionary biology. Here, we investigated the evolution of cold acclimation in Pooideae and found that a phylogenetically diverse set of Pooideae species displayed cold acclimation capacity. However, comparing differential gene expression after cold treatment in transcriptomes of five phylogenetically diverse species revealed widespread species-specific responses of genes with conserved sequences. Furthermore, we studied the correlation between gene family size and number of cold-responsive genes as well as between selection pressure on coding sequences of genes and their cold responsiveness. We saw evidence of protein-coding and regulatory sequence evolution as well as the origin of novel genes and functions contributing toward evolution of a cold response in Pooideae. Our results reflect that selection pressure resulting from global cooling must have acted on already diverged lineages. Nevertheless, conservation of cold-induced gene expression of certain genes indicates that the Pooideae ancestor may have possessed some molecular machinery to mitigate cold stress. Evolution of adaptations to seasonally cold climates is regarded as particularly difficult. How Pooideae evolved to transition from tropical to temperate biomes sheds light on how complex traits evolve in the light of climate changes.

The long-term evolutionary impacts of whole genome duplication (WGD) are strongly influenced by the ensuing rediploidization process. Following autopolyploidization, rediploidization involves a transition from tetraploid to diploid meiotic pairing, allowing duplicated genes (ohnologues) to diverge genetically and functionally. Our understanding of autopolyploid rediploidization has been informed by a WGD event ancestral to salmonid fishes, where large genomic regions are characterized by temporally delayed rediploidization, allowing lineage-specific ohnologue sequence divergence in the major salmonid clades. Here, we investigate the long-term outcomes of autopolyploid rediploidization at genome-wide resolution, exploiting a recent 'explosion' of salmonid genome assemblies, including a new genome sequence for the huchen (Hucho hucho). We developed a genome alignment approach to capture duplicated regions across multiple species, allowing us to create 121,864 phylogenetic trees describing ohnologue divergence across salmonid evolution. Using molecular clock analysis, we show that 61% of the ancestral salmonid genome experienced an initial 'wave' of rediploidization in the late Cretaceous (85-106 Mya). This was followed by a period of relative genomic stasis lasting 17-39 My, where much of the genome remained in a tetraploid state. A second rediploidization wave began in the early Eocene and proceeded alongside species diversification, generating predictable patterns of lineage-specific ohnologue divergence, scaling in complexity with the number of speciation events. Finally, using gene set enrichment, gene expression, and codon-based selection analyses, we provide insights into potential functional outcomes of delayed rediploidization. Overall, this study enhances our understanding of delayed autopolyploid rediploidization and has broad implications for future studies of WGD events.

Background One objective of this study was to identify putative quantitative trait loci (QTL) that affect indicator phenotypes for growth, nitrogen, and carbon metabolism in muscle, liver, and adipose tissue, and for feed efficiency. Another objective was to perform an RNAseq analysis (184 fish from all families), to identify genes that are associated with carbon and nitrogen metabolism in the liver. The material consisted of a family experiment that was performed in freshwater and included 2281 individuals from 23 full-sib families. During the 12-day feed conversion test, families were randomly allocated to family tanks (50 fish per tank and 2 tanks per family) and fed a fishmeal-based diet labeled with the stable isotopes 15N and 13C at inclusion levels of 2 and 1%, respectively. Results Using a linear mixed-model algorithm, a QTL for pre-smolt growth was identified on chromosome 9 and a QTL for carbon metabolism in the liver was identified on chromosome 12 that was closely related to feed conversion ratio on a tank level. For the indicators of feed efficiency traits that were derived from the stable isotope ratios (15N and 13C) of muscle tissue and growth, no convincing QTL was detected, which suggests that these traits are polygenic. The transcriptomic analysis showed that high carbon and nitrogen metabolism was associated with individuals that convert protein from the feed more efficiently, primarily due to higher expression of the proteasome, lipid, and carbon metabolic pathways in liver. In addition, we identified seven transcription factors that were associated with carbon and nitrogen metabolism and located in the identified QTL regions. Conclusions Analyses revealed one QTL associated with pre-smolt growth and one QTL for carbon metabolism in the liver. Both of these traits are associated with feed efficiency. However, more accurate mapping of the putative QTL will require a more diverse family material. In this experiment, fish that have a high carbon and nitrogen metabolism in the liver converted protein from the feed more efficiently, potentially because of a higher expression of the proteasome, lipid, and carbon metabolic pathways in liver. Within the QTL regions, we detected seven transcription factors that were associated with carbon and nitrogen metabolism.