Comparative regulomics supports pervasive selection on gene dosage following whole genome duplication

Evolutionary screening of precision oncology biomarkers and its applications in prognostic model construction

Biomarker screening is critical for precision oncology. However, one of the main challenges in precision oncology is that the screened biomarkers often fail to achieve the expected clinical effects and are rarely approved by regulatory authorities. Considering the close association between cancer pathogenesis and the evolutionary events of organisms, we first explored the evolutionary feature underlying clinically approved biomarkers, and two evolutionary features of approved biomarkers (Ohnologs and specific evolutionary stages of genes) were identified. Subsequently, we utilized evolutionary features for screening potential prognostic biomarkers in four common cancers: head and neck squamous cell carcinoma, liver hepatocellular carcinoma, lung adenocarcinoma, and lung squamous cell carcinoma. Finally, we constructed an evolution-strengthened prognostic model (ESPM) for cancers. These models can predict cancer patients' survival time across different cancer cohorts effectively and perform better than conventional models. In summary, our study highlights the application potentials of evolutionary information in precision oncology biomarker screening.

Draft genome sequence and tissue expression panel of Pacific saury (Cololabis saira)

Pacific saury (Cololabis saira) is an important fish in several countries. Notably, the catch of this fish has markedly decreased recently, which might be due to environmental changes, including feeding habitat changes. However, no clear correlation has been observed. Therefore, the physiological basis of Pacific saury in relation to possible environmental factors must be understood. We sequenced the genome of Pacific saury and extracted RNA from nine tissues (brain, eye, gill, anterior/posterior guts, kidney, liver, muscle, and ovary). In 1.09 Gb assembled genome sequences, a total of 26,775 protein-coding genes were predicted, of which 26,241 genes were similar to known genes in a public database. Transcriptome analysis revealed that 24,254 genes were expressed in at least one of the nine tissues, and 7,495 were highly expressed in specific tissues. Based on the similarity of the expression profiles to those of model organisms, the transcriptome obtained was validated to reflect the characteristics of each tissue. Thus, the present genomic and transcriptomic data serve as useful resources for molecular studies on Pacific saury. In particular, we emphasize that the gene expression data, which serve as the tissue expression panel of this species, can be employed in comparative transcriptomics on marine environmental responses.

Gene Function is a Driver of Activin Signaling Pathway Evolution Following Whole-Genome Duplication in Rainbow Trout (Oncorhynchus mykiss)

The genomes of plant and animal species are influenced by ancestral whole-genome duplication (WGD) events, which have profound impacts on the regulation and function of gene networks. To gain insight into the consequences of WGD events, we characterized the sequence conservation and expression patterns of ohnologs in the highly duplicated activin receptor signaling pathway in rainbow trout (RBT). The RBT activin receptor signaling pathway is defined by tissue-specific expression of inhibitors and ligands and broad expression of receptors and Co-Smad signaling molecules. Signaling pathway ligands exhibited shared expression, while inhibitors and Smad signaling molecules primarily express a single dominant ohnolog. Our findings suggest that gene function influences ohnolog evolution following duplication of the activin signaling pathway in RBT.

Identification and effective regulation of scarb1 gene involved in pigmentation change in autotetraploid Carassius auratus

The autotetraploid Carassius auratus (4nRR, 4 n=200, RRRR) is derived from whole-genome duplication of Carassius auratus red var. (RCC, 2 n=100, RR). In the current study, we demonstrated that chromatophores and pigment changes directly caused the coloration and variation of 4nRR skin (red in RCC, brownish-yellow in 4nRR). To further explore the molecular mechanisms underlying coloration formation and variation in 4nRR, we performed transcriptome profiling and molecular functional verification in RCC and 4nRR. Results revealed that scarb1, associated with carotenoid metabolism, underwent significant down-regulation in 4nRR. Efficient editing of this candidate pigment gene provided clear evidence of its significant role in RCC coloration. Subsequently, we identified four divergent scarb1 homeologs in 4nRR: two original scarb1 homeologs from RCC and two duplicated ones. Notably, three of these homeologs possessed two highly conserved alleles, exhibiting biased and allele-specific expression in the skin. Remarkably, after precise editing of both the original and duplicated scarb1 homeologs and/or alleles, 4nRR individuals, whether singly or multiply mutated, displayed a transition from brownish-yellow skin to a cyan-gray phenotype. Concurrently, the proportional areas of the cyan-gray regions displayed a gene-dose correlation. These findings illustrate the subfunctionalization of duplicated scarb1, with all scarb1 genes synergistically and equally contributing to the pigmentation of 4nRR. This is the first report concerning the functional differentiation of duplicated homeologs in an autopolyploid fish, substantially enriching our understanding of coloration formation and change within this group of organisms.

Models for the retention of duplicate genes and their biological underpinnings

Gene content in genomes changes through several different processes, with gene duplication being an important contributor to such changes. Gene duplication occurs over a range of scales from individual genes to whole genomes, and the dynamics of this process can be context dependent. Still, there are rules by which genes are retained or lost from genomes after duplication, and probabilistic modeling has enabled characterization of these rules, including their context-dependence. Here, we describe the biology and corresponding mathematical models that are used to understand duplicate gene retention and its contribution to the set of biochemical functions encoded in a genome.

Intergrative metabolomic and transcriptomic analyses reveal the potential regulatory mechanism of unique dihydroxy fatty acid biosynthesis in the seeds of an industrial oilseed crop Orychophragmus violaceus

BACKGROUND

Orychophragmus violaceus is a potentially important industrial oilseed crop due to the two 24-carbon dihydroxy fatty acids (diOH-FA) that was newly identified from its seed oil via a 'discontinuous elongation' process. Although many research efforts have focused on the diOH-FA biosynthesis mechanism and identified the potential co-expressed diacylglycerol acyltranferase (DGAT) gene associated with triacylglycerol (TAG)-polyestolides biosynthesis, the dynamics of metabolic changes during seed development of O. violaceus as well as its associated regulatory network changes are poorly understood.

RESULTS

In this study, by combining metabolome and transcriptome analysis, we identified that 1,003 metabolites and 22,479 genes were active across four stages of seed development, which were further divided into three main clusters based on the patterns of metabolite accumulation and/or gene expression. Among which, cluster2 was mostly related to diOH-FA biosynthesis pathway. We thus further constructed transcription factor (TF)-structural genes regulatory map for the genes associated with the flavonoids, fatty acids and diOH-FA biosynthesis pathway in this cluster. In particular, several TF families such as bHLH, B3, HD-ZIP, MYB were found to potentially regulate the metabolism associated with the diOH-FA pathway. Among which, multiple candidate TFs with promising potential for increasing the diOH-FA content were identified, and we further traced the evolutionary history of these key genes among species of Brassicaceae.

CONCLUSION

Taken together, our study provides new insight into the gene resources and potential relevant regulatory mechanisms of diOH-FA biosynthesis uniquely in seeds of O. violaceus, which will help to promote the downstream breeding efforts of this potential oilseed crop and advance the bio-lubricant industry.

Expectations of duplicate gene retention under the gene duplicability hypothesis

BACKGROUND

Gene duplication is an important process in evolution. What causes some genes to be retained after duplication and others to be lost is a process not well understood. The most prevalent theory is the gene duplicability hypothesis, that something about the function and number of interacting partners (number of subunits of protein complex, etc.), determines whether copies have more opportunity to be retained for long evolutionary periods. Some genes are also more susceptible to dosage balance effects following WGD events, making them more likely to be retained for longer periods of time. One would expect these processes that affect the retention of duplicate copies to affect the conditional probability ratio after consecutive whole genome duplication events. The probability that a gene will be retained after a second whole genome duplication event (WGD2), given that it was retained after the first whole genome duplication event (WGD1) versus the probability a gene will be retained after WGD2, given it was lost after WGD1 defines the probability ratio that is calculated.

RESULTS

Since duplicate gene retention is a time heterogeneous process, the time between the events (t1) and the time since the most recent event (t2) are relevant factors in calculating the expectation for observation in any genome. Here, we use a survival analysis framework to predict the probability ratio for genomes with different values of t1 and t2 under the gene duplicability hypothesis, that some genes are more susceptible to selectable functional shifts, some more susceptible to dosage compensation, and others only drifting. We also predict the probability ratio with different values of t1 and t2 under the mutational opportunity hypothesis, that probability of retention for certain genes changes in subsequent events depending upon how they were previously retained. These models are nested such that the mutational opportunity model encompasses the gene duplicability model with shifting duplicability over time. Here we present a formalization of the gene duplicability and mutational opportunity hypotheses to characterize evolutionary dynamics and explanatory power in a recently developed statistical framework.

CONCLUSIONS

This work presents expectations of the gene duplicability and mutational opportunity hypotheses over time under different sets of assumptions. This expectation will enable formal testing of processes leading to duplicate gene retention.

Epigenomic divergence correlates with sequence polymorphism in Arabidopsis paralogs

Processes affecting rates of sequence polymorphism are fundamental to the evolution of gene duplicates. The relationship between gene activity and sequence polymorphism can influence the likelihood that functionally redundant gene copies are co-maintained in stable evolutionary equilibria vs other outcomes such as neofunctionalization. Here, we investigate genic variation in epigenome-associated polymorphism rates in Arabidopsis thaliana and consider whether these affect the evolution of gene duplicates. We compared the frequency of sequence polymorphism and patterns of genetic differentiation between genes classified by exon methylation patterns: unmethylated (unM), gene-body methylated (gbM), and transposon-like methylated (teM) states, which reflect divergence in gene expression. We found that the frequency of polymorphism was higher in teM (transcriptionally repressed, tissue-specific) genes and lower in gbM (active, constitutively expressed) genes. Comparisons of gene duplicates were largely consistent with genome-wide patterns - gene copies that exhibit teM accumulate more variation, evolve faster, and are in chromatin states associated with reduced DNA repair. This relationship between expression, the epigenome, and polymorphism may lead to the breakdown of equilibrium states that would otherwise maintain genetic redundancies. Epigenome-mediated polymorphism rate variation may facilitate the evolution of novel gene functions in duplicate paralogs maintained over evolutionary time.

Allele surfing causes maladaptation in a Pacific salmon of conservation concern

How various factors, including demography, recombination or genome duplication, may impact the efficacy of natural selection and the burden of deleterious mutations, is a central question in evolutionary biology and genetics. In this study, we show that key evolutionary processes, including variations in i) effective population size (Ne) ii) recombination rates and iii) chromosome inheritance, have influenced the genetic load and efficacy of selection in Coho salmon (Oncorhynchus kisutch), a widely distributed salmonid species on the west coast of North America. Using whole genome resequencing data from 14 populations at different migratory distances from their southern glacial refugium, we found evidence supporting gene surfing, wherein reduced Ne at the postglacial recolonization front, leads to a decrease in the efficacy of selection and a surf of deleterious alleles in the northernmost populations. Furthermore, our results indicate that recombination rates play a prime role in shaping the load along the genome. Additionally, we identified variation in polyploidy as a contributing factor to within-genome variation of the load. Overall, our results align remarkably well with expectations under the nearly neutral theory of molecular evolution. We discuss the fundamental and applied implications of these findings for evolutionary and conservation genomics.

Expression dosage effects of a small number of genes after the artificial doubling of weeping forsythia

Whole genome doubling (WGD) plays a critical role in plant evolution, yet the mechanisms underlying the maintenance of overall equilibrium following an artificial doubling event, as well as its impact on phenotype and adaptability, remain unclear. By comparing the gene expression of naturally occurring weeping forsythia diploids and colchicine-induced autotetraploids under normal growth conditions and cold stress, we identified gene expression dosage responses resulting from ploidy change. Only a small proportion of effectively expressed genes showed dosage effect, and most genes did not exhibit significant expression differences. However, the genes that showed expression dosage effect were largely random. The autotetraploids had slower overall growth rates, possibly resulting from negative gene dosage effects on zeatin synthesis and multiple metabolic delays caused by other negative dosage genes. Our comparative analysis of cold response genes in diploids and autotetraploids revealed that genes related to "response to abscisic acid" and "cold acclimation" were key factors contributing to greater cold tolerance in the autotetraploids. In particular, gene expression related to "cold acclimation" might mitigate the effects of cold stress. Taken together, our findings suggested that overall gene expression equilibrium following WGD of weeping forsythia autotetraploids was achieved through the inactivation of the majority of duplicated genes. Our research provides new insights into the mechanisms regulating expression dosage balance following polyploidization events.

Phylogenetic modeling of enhancer shifts in African mole-rats reveals regulatory changes associated with tissue-specific traits

Changes in gene regulation are thought to underlie most phenotypic differences between species. For subterranean rodents such as the naked mole-rat, proposed phenotypic adaptations include hypoxia tolerance, metabolic changes, and cancer resistance. However, it is largely unknown what regulatory changes may associate with these phenotypic traits, and whether these are unique to the naked mole-rat, the mole-rat clade, or are also present in other mammals. Here, we investigate regulatory evolution in the heart and liver from two African mole-rat species and two rodent outgroups using genome-wide epigenomic profiling. First, we adapted and applied a phylogenetic modeling approach to quantitatively compare epigenomic signals at orthologous regulatory elements and identified thousands of promoter and enhancer regions with differential epigenomic activity in mole-rats. These elements associate with known mole-rat adaptations in metabolic and functional pathways and suggest candidate genetic loci that may underlie mole-rat innovations. Second, we evaluated ancestral and species-specific regulatory changes in the study phylogeny and report several candidate pathways experiencing stepwise remodeling during the evolution of mole-rats, such as the insulin and hypoxia response pathways. Third, we report nonorthologous regulatory elements overlap with lineage-specific repetitive elements and appear to modify metabolic pathways by rewiring of HNF4 and RAR/RXR transcription factor binding sites in mole-rats. These comparative analyses reveal how mole-rat regulatory evolution informs previously reported phenotypic adaptations. Moreover, the phylogenetic modeling framework we propose here improves upon the state of the art by addressing known limitations of inter-species comparisons of epigenomic profiles and has broad implications in the field of comparative functional genomics.

Transposon waves at the water-to-land transition

The major transitions in vertebrate evolution are associated with significant genomic reorganizations. In contrast to the evolutionary processes that occurred at the origin of vertebrates or prior to the radiation of teleost fishes, no whole-genome duplication events occurred during the water-to-land transition, and it remains an open question how did genome dynamics contribute to this prominent evolutionary event. Indeed, the recent sequencing of sarcopterygian and amphibian genomes has revealed that the extant lineages immediately preceding and succeeding this transition harbor an exceptional number of transposable elements and it is tempting to speculate that these sequences might have catalyzed the adaptations that enabled vertebrates to venture into land. Here, we review the genome dynamics associated with the major transitions in vertebrate evolution and discuss how the highly repetitive genomic landscapes revealed by recent efforts to characterize the genomes of amphibians and sarcopterygians argue for turbulent genome dynamics occurring before the water-to-land transition and possibly enabling it.

Genomic copy number variability at the genus, species and population levels impacts in situ ecological analyses of dinoflagellates and harmful algal blooms

The application of meta-barcoding, qPCR, and metagenomics to aquatic eukaryotic microbial communities requires knowledge of genomic copy number variability (CNV). CNV may be particularly relevant to functional genes, impacting dosage and expression, yet little is known of the scale and role of CNV in microbial eukaryotes. Here, we quantify CNV of rRNA and a gene involved in Paralytic Shellfish Toxin (PST) synthesis (sxtA4), in 51 strains of 4 Alexandrium (Dinophyceae) species. Genomes varied up to threefold within species and ~7-fold amongst species, with the largest (A. pacificum, 130 ± 1.3 pg cell-1 /~127 Gbp) in the largest size category of any eukaryote. Genomic copy numbers (GCN) of rRNA varied by 6 orders of magnitude amongst Alexandrium (102- 108 copies cell-1) and were significantly related to genome size. Within the population CNV of rRNA was 2 orders of magnitude (105 - 107 cell-1) in 15 isolates from one population, demonstrating that quantitative data based on rRNA genes needs considerable caution in interpretation, even if validated against locally isolated strains. Despite up to 30 years in laboratory culture, rRNA CNV and genome size variability were not correlated with time in culture. Cell volume was only weakly associated with rRNA GCN (20-22% variance explained across dinoflagellates, 4% in Gonyaulacales). GCN of sxtA4 varied from 0-102 copies cell-1, was significantly related to PSTs (ng cell-1), displaying a gene dosage effect modulating PST production. Our data indicate that in dinoflagellates, a major marine eukaryotic group, low-copy functional genes are more reliable and informative targets for quantification of ecological processes than unstable rRNA genes.

Gene family expansion and functional diversification of chitinase and chitin synthase genes in Atlantic salmon (Salmo salar)

Chitin is one of the most abundant polysaccharides in nature, forming important structures in insects, crustaceans, and fungal cell walls. Vertebrates on the other hand are generally considered "nonchitinous" organisms, despite having highly conserved chitin metabolism-associated genes. Recent work has revealed that the largest group of vertebrates, the teleosts, have the potential to both synthesize and degrade endogenous chitin. Yet, little is known about the genes and proteins responsible for these dynamic processes. Here, we used comparative genomics, transcriptomics, and chromatin accessibility data to characterize the repertoire, evolution, and regulation of genes involved in chitin metabolism in teleosts, with a particular focus on Atlantic salmon. Reconstruction of gene family phylogenies provides evidence for an expansion of teleost and salmonid chitinase and chitin synthase genes after multiple whole-genome duplications. Analyses of multi-tissue gene expression data demonstrated a strong bias of gastrointestinal tract expression for chitin metabolism genes, but with different spatial and temporal tissue specificities. Finally, we integrated transcriptomes from a developmental time series of the gastrointestinal tract with chromatin accessibility data to identify putative transcription factors responsible for regulating chitin metabolism gene expression (CDX1 and CDX2) as well as tissue-specific divergence in the regulation of gene duplicates (FOXJ2). The findings presented here support the hypothesis that chitin metabolism genes in teleosts play a role in developing and maintaining a chitin-based barrier in the teleost gut and provide a basis for further investigations into the molecular basis of this barrier.

Dosage balance acts as a time-dependent selective barrier to subfunctionalization

BACKGROUND

Gene duplication is an important process for genome expansion, sometimes allowing for new gene functions to develop. Duplicate genes can be retained through multiple processes, either for intermediate periods of time through processes such as dosage balance, or over extended periods of time through processes such as subfunctionalization and neofunctionalization.

RESULTS

Here, we built upon an existing subfunctionalization Markov model by incorporating dosage balance to describe the interplay between subfunctionalization and dosage balance to explore selective pressures on duplicate copies. Our model incorporates dosage balance using a biophysical framework that penalizes the fitness of genetic states with stoichiometrically imbalanced proteins. These imbalanced states cause increased concentrations of exposed hydrophobic surface areas, which cause deleterious mis-interactions. We draw comparison between our Subfunctionalization + Dosage-Balance Model (Sub + Dos) and the previous Subfunctionalization-Only (Sub-Only) Model. This comparison includes how the retention probabilities change over time, dependent upon the effective population size and the selective cost associated with spurious interaction of dosage-imbalanced partners. We show comparison between Sub-Only and Sub + Dos models for both whole-genome duplication and small-scale duplication events.

CONCLUSION

These comparisons show that following whole-genome duplication, dosage balance serves as a time-dependent selective barrier to the subfunctionalization process, by causing an overall delay but ultimately leading to a larger portion of the genome retained through subfunctionalization. This higher percentage of the genome that is ultimately retained is caused by the alternative competing process, nonfunctionalization, being selectively blocked to a greater extent. In small-scale duplication, the reverse pattern is seen, where dosage balance drives faster rates of subfunctionalization, but ultimately leads to a smaller portion of the genome retained as duplicates. This faster rate of subfunctionalization is because the dosage balance of interacting gene products is negatively affected immediately after duplication and the loss of a duplicate restores the stoichiometric balance. Our findings provide support that the subfunctionalization of genes that are susceptible to dosage balance effects, such as proteins involved in complexes, is not a purely neutral process. With stronger selection against stoichiometrically imbalanced gene partners, the rates of subfunctionalization and nonfunctionalization slow; however, this ultimately leads to a greater proportion of subfunctionalized gene pairs.

Analyses of Genome Regulatory Evolution Following Whole-Genome Duplication Using the Phylogenetic EVE Model

Whole-genome duplications (WGDs) are important in shaping the evolution of complex genomes, including rewiring of genome regulation. To address key questions about how WGDs impact the evolution of genome regulation, we need to understand the relative importance of selection versus drift and temporal evolutionary dynamics. One promising class of statistical models that can help address such questions are phylogenetic Ornstein-Uhlenbeck (OU) models.Here we present a computational pipeline for the comparative phylogenetic analyses of genome regulation using an OU model. We have implemented this model in R and provide a step-by-step protocol for the use of this model, including example scripts and simulated test data. We provide the nonspecialist a brief overview of how this model works and how to perform tests for signatures of selection on genome regulation as well as power simulations to aid in experimental design and interpretation of results. We believe that these resources could help polyploidy research move forward in an era of rapidly increasing functional genomics data across the tree of life.

Molecular and cellular constraints on vesicle traffic evolution

In eukaryotic cells, the budding and fusion of intracellular transport vesicles is carefully orchestrated in space and time. Locally, a vesicle's source compartment, its cargo, and its destination compartment are controlled by dynamic multi-protein specificity modules. Globally, vesicle constituents must be recycled to ensure homeostasis of compartment compositions. The emergence of a novel vesicle pathway therefore requires new specificity modules as well as new recycling routes. Here, we review recent research on local (molecular) constraints on gene module duplication and global (cellular) constraints on intracellular recycling. By studying the evolution of vesicle traffic, we may discover general principles of how complex traits arise via multiple intermediate steps.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

WGDTree: a phylogenetic software tool to examine conditional probabilities of retention following whole genome duplication events

BACKGROUND

Multiple processes impact the probability of retention of individual genes following whole genome duplication (WGD) events. In analyzing two consecutive whole genome duplication events that occurred in the lineage leading to Atlantic salmon, a new phylogenetic statistical analysis was developed to examine the contingency of retention in one event based upon retention in a previous event. This analysis is intended to evaluate mechanisms of duplicate gene retention and to provide software to generate the test statistic for any genome with pairs of WGDs in its history.

RESULTS

Here a software package written in Python, 'WGDTree' for the analysis of duplicate gene retention following whole genome duplication events is presented. Using gene tree-species tree reconciliation to label gene duplicate nodes and differentiate between WGD and SSD duplicates, the tool calculates a statistic based upon the conditional probability of a gene duplicate being retained after a second whole genome duplication dependent upon the retention status after the first event. The package also contains methods for the simulation of gene trees with WGD events. After running simulations, the accuracy of the placement of events has been determined to be high. The conditional probability statistic has been calculated for Phalaenopsis equestris on a monocot species tree with a pair of consecutive WGD events on its lineage, showing the applicability of the method.

CONCLUSIONS

A new software tool has been created for the analysis of duplicate genes in examination of retention mechanisms. The software tool has been made available on the Python package index and the source code can be found on GitHub here: https://github.com/cnickh/wgdtree .

Sublinear scaling of the cellular proteome with ploidy

Ploidy changes are frequent in nature and contribute to evolution, functional specialization and tumorigenesis. Analysis of model organisms of different ploidies revealed that increased ploidy leads to an increase in cell and nuclear volume, reduced proliferation, metabolic changes, lower fitness, and increased genomic instability, but the underlying mechanisms remain poorly understood. To investigate how gene expression changes with cellular ploidy, we analyzed isogenic series of budding yeasts from 1N to 4N. We show that mRNA and protein abundance scales allometrically with ploidy, with tetraploid cells showing only threefold increase in protein abundance compared to haploids. This ploidy-dependent sublinear scaling occurs via decreased rRNA and ribosomal protein abundance and reduced translation. We demonstrate that the activity of Tor1 is reduced with increasing ploidy, which leads to diminished rRNA gene repression via a Tor1-Sch9-Tup1 signaling pathway. mTORC1 and S6K activity are also reduced in human tetraploid cells and the concomitant increase of the Tup1 homolog Tle1 downregulates the rDNA transcription. Our results suggest that the mTORC1-Sch9/S6K-Tup1/TLE1 pathway ensures proteome remodeling in response to increased ploidy.

The landscape of transcriptional and translational changes over 22 years of bacterial adaptation

Organisms can adapt to an environment by taking multiple mutational paths. This redundancy at the genetic level, where many mutations have similar phenotypic and fitness effects, can make untangling the molecular mechanisms of complex adaptations difficult. Here, we use the Escherichia coli long-term evolution experiment (LTEE) as a model to address this challenge. To understand how different genomic changes could lead to parallel fitness gains, we characterize the landscape of transcriptional and translational changes across 12 replicate populations evolving in parallel for 50,000 generations. By quantifying absolute changes in mRNA abundances, we show that not only do all evolved lines have more mRNAs but that this increase in mRNA abundance scales with cell size. We also find that despite few shared mutations at the genetic level, clones from replicate populations in the LTEE are remarkably similar in their gene expression patterns at both the transcriptional and translational levels. Furthermore, we show that the majority of the expression changes are due to changes at the transcriptional level with very few translational changes. Finally, we show how mutations in transcriptional regulators lead to consistent and parallel changes in the expression levels of downstream genes. These results deepen our understanding of the molecular mechanisms underlying complex adaptations and provide insights into the repeatability of evolution.

The HASTER lncRNA promoter is a cis-acting transcriptional stabilizer of HNF1A

The biological purpose of long non-coding RNAs (lncRNAs) is poorly understood. Haploinsufficient mutations in HNF1A homeobox A (HNF1A), encoding a homeodomain transcription factor, cause diabetes mellitus. Here, we examine HASTER, the promoter of an lncRNA antisense to HNF1A. Using mouse and human models, we show that HASTER maintains cell-specific physiological HNF1A concentrations through positive and negative feedback loops. Pancreatic β cells from Haster mutant mice consequently showed variegated HNF1A silencing or overexpression, resulting in hyperglycaemia. HASTER-dependent negative feedback was essential to prevent HNF1A binding to inappropriate genomic regions. We demonstrate that the HASTER promoter DNA, rather than the lncRNA, modulates HNF1A promoter-enhancer interactions in cis and thereby regulates HNF1A transcription. Our studies expose a cis-regulatory element that is unlike classic enhancers or silencers, it stabilizes the transcription of its target gene and ensures the fidelity of a cell-specific transcription factor program. They also show that disruption of a mammalian lncRNA promoter can cause diabetes mellitus.

Single cell transcriptomics of Atlantic salmon (Salmo salar L.) liver reveals cellular heterogeneity and immunological responses to challenge by Aeromonas salmonicida

The liver is a multitasking organ with essential functions for vertebrate health spanning metabolism and immunity. In contrast to mammals, our understanding of liver cellular heterogeneity and its role in regulating immunological status remains poorly defined in fishes. Addressing this knowledge gap, we generated a transcriptomic atlas of 47,432 nuclei isolated from the liver of Atlantic salmon (Salmo salar L.) contrasting control fish with those challenged with a pathogenic strain of Aeromonas salmonicida, a problematic bacterial pathogen in global aquaculture. We identified the major liver cell types and their sub-populations, revealing poor conservation of many hepatic cell marker genes utilized in mammals, while identifying novel heterogeneity within the hepatocyte, lymphoid, and myeloid lineages. This included polyploid hepatocytes, multiple T cell populations including γδ T cells, and candidate populations of monocytes/macrophages and dendritic cells. A dominant hepatocyte population radically remodeled its transcriptome following infection to activate the acute phase response and other defense functions, while repressing routine functions such as metabolism. These defense-specialized hepatocytes showed strong activation of genes controlling protein synthesis and secretion, presumably to support the release of acute phase proteins into circulation. The infection response further involved up-regulation of numerous genes in an immune-cell specific manner, reflecting functions in pathogen recognition and killing, antigen presentation, phagocytosis, regulation of inflammation, B cell differentiation and T cell activation. Overall, this study greatly enhances our understanding of the multifaceted role played by liver immune and non-immune cells in host defense and metabolic remodeling following infection and provides many novel cell-specific marker genes to empower future studies of this organ in fishes.

Pulmonate slug evolution is reflected in the de novo genome of Arion vulgaris Moquin-Tandon, 1855

Stylommatophoran pulmonate land slugs and snails successfully completed the water-to-land transition from an aquatic ancestor and flourished on land. Of the 30,000 estimated species, very few genomes have so far been published. Here, we assembled and characterized a chromosome-level genome of the “Spanish” slug, Arion vulgaris Moquin-Tandon, 1855, a notorious pest land slug in Europe. Using this reference genome, we conclude that a whole-genome duplication event occurred approximately 93–109 Mya at the base of Stylommatophora and might have promoted land invasion and adaptive radiation. Comparative genomic analyses reveal that genes related to the development of kidney, blood vessels, muscle, and nervous systems had expanded in the last common ancestor of land pulmonates, likely an evolutionary response to the terrestrial challenges of gravity and water loss. Analyses of A. vulgaris gene families and positively selected genes show the slug has evolved a stronger ability to counteract the greater threats of external damage, radiation, and water loss lacking a protective shell. Furthermore, a recent burst of long interspersed elements in the genome of A. vulgaris might affect gene regulation and contribute to rapid phenotype changes in A. vulgaris, which might be conducive to its rapid adaptation and invasiveness.

Molecular Regulation of Biosynthesis of Long Chain Polyunsaturated Fatty Acids in Atlantic Salmon

Salmon is a rich source of health-promoting omega-3 long chain polyunsaturated fatty acids (n-3 LC-PUFA), such as eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3). The LC-PUFA biosynthetic pathway in Atlantic salmon is one of the most studied compared to other teleosts. This has largely been due to the massive replacement of LC-PUFA-rich ingredients in aquafeeds with terrestrial plant oils devoid of these essential fatty acids (EFA) which ultimately pushed dietary content towards the minimal requirement of EFA. The practice would also reduce tissue content of n-3 LC-PUFA compromising the nutritional value of salmon to the human consumer. These necessitated detailed studies of endogenous biosynthetic capability as a contributor to these EFA. This review seeks to provide a comprehensive and concise overview of the current knowledge about the molecular genetics of PUFA biosynthesis in Atlantic salmon, highlighting the enzymology and nutritional regulation as well as transcriptional control networks. Furthermore, we discuss the impact of genome duplication on the complexity of salmon LC-PUFA pathway and highlight probable implications on endogenous biosynthetic capabilities. Finally, we have also compiled and made available a large RNAseq dataset from 316 salmon liver samples together with an R-script visualization resource to aid in explorative and hypothesis-driven research into salmon lipid metabolism.

Organellar transcripts dominate the cellular mRNA pool across plants of varying ploidy levels

Mitochondrial and plastid functions depend on coordinated expression of proteins encoded by genomic compartments that have radical differences in copy number of organellar and nuclear genomes. In polyploids, doubling of the nuclear genome may add challenges to maintaining balanced expression of proteins involved in cytonuclear interactions. Here, we use ribo-depleted RNA sequencing (RNA-seq) to analyze transcript abundance for nuclear and organellar genomes in leaf tissue from four different polyploid angiosperms and their close diploid relatives. We find that even though plastid genomes contain <1% of the number of genes in the nuclear genome, they generate the majority (69.9 to 82.3%) of messenger RNA (mRNA) transcripts in the cell. Mitochondrial genes are responsible for a much smaller percentage (1.3 to 3.7%) of the leaf mRNA pool but still produce much higher transcript abundances per gene compared to nuclear genome. Nuclear genes encoding proteins that functionally interact with mitochondrial or plastid gene products exhibit mRNA expression levels that are consistently more than 10-fold lower than their organellar counterparts, indicating an extreme cytonuclear imbalance at the RNA level despite the predominance of equimolar interactions at the protein level. Nevertheless, interacting nuclear and organellar genes show strongly correlated transcript abundances across functional categories, suggesting that the observed mRNA stoichiometric imbalance does not preclude coordination of cytonuclear expression. Finally, we show that nuclear genome doubling does not alter the cytonuclear expression ratios observed in diploid relatives in consistent or systematic ways, indicating that successful polyploid plants are able to compensate for cytonuclear perturbations associated with nuclear genome doubling.

Three problems in the genetics of speciation by selection

Speciation is the process by which barriers to gene flow evolve between populations. Although we now know that speciation is largely driven by natural selection, knowledge of the agents of selection and the genetic and genomic mechanisms that facilitate divergence is required for a satisfactory theory of speciation. In this essay, we highlight three advances/problems in our understanding of speciation that have arisen from studies of the genes and genomic regions that underlie the evolution of reproductive isolation. First, we describe how the identification of “speciation” genes makes it possible to identify the agents of selection causing the evolution of reproductive isolation, while also noting that the link between the genetics of phenotypic divergence and intrinsic postzygotic reproductive barriers remains tenuous. Second, we discuss the important role of recombination suppressors in facilitating speciation with gene flow, but point out that the means and timing by which reproductive barriers become associated with recombination cold spots remains uncertain. Third, we establish the importance of ancient genetic variation in speciation, although we argue that the focus of speciation studies on evolutionarily young groups may bias conclusions in favor of ancient variation relative to new mutations.

Molecular Evolution of the Activating Transcription Factors Shapes the Adaptive Cellular Responses to Oxidative Stress

Reactive oxygen species (ROS) play an essential part in physiology of individual cell. ROS can cause damage to various biomolecules, including DNA. The systems that have developed to harness the impacts of ROS are antique evolutionary adaptations that are intricately linked to almost every aspect of cellular function. This research reveals the idea that during evolution, rather than being largely conserved, the molecular pathways reacting to oxidative stress have intrinsic flexibility. The coding sequences of the ATF2, ATF3, ATF4, and ATF6 genes were aligned to examine selection pressure on the genes, which were shown to be very highly conserved among vertebrate species. A total of 33 branches were explicitly evaluated for their capacity to diversify selection. After accounting for multiple testing, significance was determined using the likelihood ratio test with a threshold of p ≤ 0.05. Positive selection signs in these genes were detected across vertebrate lineages. In the selected test branches of our phylogeny, the synonymous rate variation revealed evidence (LRT, p value = 0.011 ≤ 0.05) of gene-wide episodic diversifying selection. As a result, there is evidence that diversifying selection occurred at least once on at least one test branch. These findings indicate that the activities of ROS-responsive systems are also theoretically flexible and may be altered by environmental selection pressure. By determining where the genes encoding these processes are "targeted" during evolution, we may better understand the mechanism of adaptation to oxidative stress during evolution.

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.

Structural and Evolutionary Adaptations of Nei-Like DNA Glycosylases Proteins Involved in Base Excision Repair of Oxidative DNA Damage in Vertebrates

Oxidative stress is a type of stress that damages DNA and can occur from both endogenous and exogenous sources. Damage to DNA caused by oxidative stress can result in base modifications that promote replication errors and the formation of sites of base loss, which pose unique challenges to the preservation of genomic integrity. However, the adaptive evolution of the DNA repair mechanism is poorly understood in vertebrates. This research aimed to explore the evolutionary relationships, physicochemical characteristics, and comparative genomic analysis of the Nei-like glycosylase gene family involved in DNA base repair in the vertebrates. The genomic sequences of NEIL1, NEIL2, and NEIL3 genes were aligned to observe selection constraints in the genes, which were relatively low conserved across vertebrate species. The positive selection signals were identified in these genes across the vertebrate lineages. We identified that only about 2.7% of codons in these genes were subjected to positive selection. We also revealed that positive selection pressure was increased in the Fapy-DNA-glyco and H2TH domain, which are involved in the base excision repair of DNA that has been damaged by oxidative stress. Gene structure, motif, and conserved domain analysis indicated that the Nei-like glycosylase genes in mammals and avians are evolutionarily low conserved compared to other glycosylase genes in other “vertebrates” species. This study revealed that adaptive selection played a critical role in the evolution of Nei-like glycosylase in vertebrate species. Systematic comparative genome analyses will give key insights to elucidate the links between DNA repair and the development of lifespan in various organisms as more diverse vertebrate genome sequences become accessible.

Adaptive Variation in Homolog Number Within Transcript Families Promotes Expression Divergence in Reef-Building Coral

Gene expression, especially in multi-species experiments, is used to gain insight into the genetic basis of how organisms adapt and respond to changing environments. However, evolutionary processes which can influence gene expression patterns between species such as the presence of paralogs which arise from gene duplication events are rarely accounted for. Paralogous transcripts can alter the transcriptional output of a gene and thus exclusion of these transcripts can obscure important biological differences between species. To address this issue, we investigated how differences in transcript family size is associated with divergent gene expression patterns in five species of Caribbean reef-building corals. We demonstrate that transcript families that are rapidly evolving in terms of size have increased levels of expression divergence. Additionally, these rapidly evolving transcript families are enriched for multiple biological processes, with genes involved in the coral innate immune system demonstrating pronounced variation in homolog number between species. Overall, this investigation demonstrates the importance of incorporating paralogous transcripts when comparing gene expression across species by influencing both transcriptional output and the number of transcripts within biological processes. As this investigation was based on transcriptome assemblies, additional insights into the relationship between gene duplications and expression patterns will likely emergence once more genome assemblies are available for study.

Genome-Wide Reconstruction of Rediploidization Following Autopolyploidization across One Hundred Million Years of Salmonid 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 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 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.

Gene family evolution underlies cell-type diversification in the hypothalamus of teleosts

Hundreds of cell types form the vertebrate brain but it is largely unknown how similar cellular repertoires are between or within species or how cell-type diversity evolves. To examine cell-type diversity across and within species, we performed single-cell RNA sequencing of ~130,000 hypothalamic cells from zebrafish (Danio rerio) and surface and cave morphs of Mexican tetra (Astyanax mexicanus). We found that over 75% of cell types were shared between zebrafish and Mexican tetra, which diverged from a common ancestor over 150 million years ago. Shared cell types displayed shifts in paralogue expression that were generated by subfunctionalization after genome duplication. Expression of terminal effector genes, such as neuropeptides, was more conserved than the expression of their associated transcriptional regulators. Species-specific cell types were enriched for the expression of species-specific genes and characterized by the neofunctionalization of expression patterns of members of recently expanded or contracted gene families. Comparisons between surface and cave morphs revealed differences in immune repertoires and transcriptional changes in neuropeptidergic cell types associated with genomic differences. The single-cell atlases presented here are a powerful resource to explore hypothalamic cell types and reveal how gene family evolution and shifts in paralogue expression contribute to cellular diversity.

Foxq2 determines blue cone identity in zebrafish

Most vertebrate lineages retain a tetrachromatic visual system, which is supported by a functional combination of spectrally distinct multiple cone photoreceptors, ultraviolet (UV), blue, green, and red cones. The blue cone identity is ensured by selective expression of blue (sws2) opsin, and the mechanism is poorly understood because sws2 gene has been lost in mammalian species such as mouse, whose visual system has been extensively studied. Here, we pursued loss-of-function studies on transcription factors expressed predominantly in zebrafish cone photoreceptors and identified Foxq2 as a blue cone–specific factor driving sws2 gene expression. Foxq2 has dual functions acting as an activator of sws2 transcription and as a suppressor of UV (sws1) opsin transcription in blue cones. A wide range of vertebrate species retain both foxq2 and sws2 genes. We propose that Foxq2-dependent sws2 expression is a prevalent regulatory mechanism that was acquired at the early stage of vertebrate evolution.

Partitioning of gene expression among zebrafish photoreceptor subtypes

Vertebrate photoreceptors are categorized into two broad classes, rods and cones, responsible for dim- and bright-light vision, respectively. While many molecular features that distinguish rods and cones are known, gene expression differences among cone subtypes remain poorly understood. Teleost fishes are renowned for the diversity of their photoreceptor systems. Here, we used single-cell RNA-seq to profile adult photoreceptors in zebrafish, a teleost. We found that in addition to the four canonical zebrafish cone types, there exist subpopulations of green and red cones (previously shown to be located in the ventral retina) that express red-shifted opsin paralogs (opn1mw4 or opn1lw1) as well as a unique combination of cone phototransduction genes. Furthermore, the expression of many paralogous phototransduction genes is partitioned among cone subtypes, analogous to the partitioning of the phototransduction paralogs between rods and cones seen across vertebrates. The partitioned cone-gene pairs arose via the teleost-specific whole-genome duplication or later clade-specific gene duplications. We also discovered that cone subtypes express distinct transcriptional regulators, including many factors not previously implicated in photoreceptor development or differentiation. Overall, our work suggests that partitioning of paralogous gene expression via the action of differentially expressed transcriptional regulators enables diversification of cone subtypes in teleosts.

A chromosome-anchored genome assembly for Lake Trout (Salvelinus namaycush)

Here, we present an annotated, chromosome‐anchored, genome assembly for Lake Trout (Salvelinus namaycush) – a highly diverse salmonid species of notable conservation concern and an excellent model for research on adaptation and speciation. We leveraged Pacific Biosciences long‐read sequencing, paired‐end Illumina sequencing, proximity ligation (Hi‐C) sequencing, and a previously published linkage map to produce a highly contiguous assembly composed of 7378 contigs (contig N50 = 1.8 Mb) assigned to 4120 scaffolds (scaffold N50 = 44.975 Mb). Long read sequencing data were generated using DNA from a female double haploid individual. 84.7% of the genome was assigned to 42 chromosome‐sized scaffolds and 93.2% of Benchmarking Universal Single Copy Orthologues were recovered, putting this assembly on par with the best currently available salmonid genomes. Estimates of genome size based on k‐mer frequency analysis were highly similar to the total size of the finished genome, suggesting that the entirety of the genome was recovered. A mitochondrial genome assembly was also produced. Self‐versus‐self synteny analysis allowed us to identify homeologs resulting from the salmonid specific autotetraploid event (Ss4R) as well as regions exhibiting delayed rediploidization. Alignment with three other salmonid genomes and the Northern Pike (Esox lucius) genome also allowed us to identify homologous chromosomes in related taxa. We also generated multiple resources useful for future genomic research on Lake Trout, including a repeat library and a sex‐averaged recombination map. A novel RNA sequencing data set for liver tissue was also generated in order to produce a publicly available set of annotations for 49,668 genes and pseudogenes. Potential applications of these resources to population genetics and the conservation of native populations are discussed.

Gene-dosage issues: a recurrent theme in whole genome duplication events

Two recent studies have addressed the long-term consequences of whole genome duplications (WGD). Specifically, they analyzed transcriptomes of the plant Arabidopsis thaliana and of four salmonids to assess the impact of WGD on gene expression. These studies point to commonalities in gene expression adjustments after polyploidization that we outline and discuss below.

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.