Negative feedback buffers effects of regulatory variants

Evolutionary trade-off and mutational bias could favor transcriptional over translational divergence within paralog pairs

How changes in the different steps of protein synthesis-transcription, translation and degradation-contribute to differences of protein abundance among genes is not fully understood. There is however accumulating evidence that transcriptional divergence might have a prominent role. Here, we show that yeast paralogous genes are more divergent in transcription than in translation. We explore two causal mechanisms for this predominance of transcriptional divergence: an evolutionary trade-off between the precision and economy of gene expression and a larger mutational target size for transcription. Performing simulations within a minimal model of post-duplication evolution, we find that both mechanisms are consistent with the observed divergence patterns. We also investigate how additional properties of the effects of mutations on gene expression, such as their asymmetry and correlation across levels of regulation, can shape the evolution of paralogs. Our results highlight the importance of fully characterizing the distributions of mutational effects on transcription and translation. They also show how general trade-offs in cellular processes and mutation bias can have far-reaching evolutionary impacts.

The missing link between genetic association and regulatory function

The genetic basis of most traits is highly polygenic and dominated by non-coding alleles. It is widely assumed that such alleles exert small regulatory effects on the expression of cis-linked genes. However, despite the availability of gene expression and epigenomic datasets, few variant-to-gene links have emerged. It is unclear whether these sparse results are due to limitations in available data and methods, or to deficiencies in the underlying assumed model. To better distinguish between these possibilities, we identified 220 gene-trait pairs in which protein-coding variants influence a complex trait or its Mendelian cognate. Despite the presence of expression quantitative trait loci near most GWAS associations, by applying a gene-based approach we found limited evidence that the baseline expression of trait-related genes explains GWAS associations, whether using colocalization methods (8% of genes implicated), transcription-wide association (2% of genes implicated), or a combination of regulatory annotations and distance (4% of genes implicated). These results contradict the hypothesis that most complex trait-associated variants coincide with homeostatic expression QTLs, suggesting that better models are needed. The field must confront this deficit and pursue this 'missing regulation.'

Antagonistic regulatory effects of a single cis-acting expression quantitative trait locus between transcription and translation of the MRPL43 gene

Background

Heterogeneity of expression quantitative trait locus (eQTL) effects have been shown across gene expression processes. Knowledge on how to produce the heterogeneity is quite limited. This study aims to examine fluctuations in differential gene expression by alleles of sequence variants across expression processes.

Results

Genome-wide eQTL analyses with transcriptome-wide gene expression data revealed 20 cis-acting eQTLs associated simultaneously with mRNA expression, ribosome occupancy, and protein abundance. A 97 kb-long eQTL signal for mitochondrial ribosomal protein L43 (MRPL43) covered the gene, showing a heterogeneous effect size on gene products across expression stages. One allele of the eQTL was associated with increased mRNA expression and ribosome occupancy but decreased protein abundance. We examined the heterogeneity and found that the eQTL can be attributed to the independent functions of three nucleotide variants, with a strong linkage. NC_000010.11:g.100987606G > T, upstream of MRPL43, may regulate the binding affinity of transcription factors. NC_000010.11:g.100986746C > G, 3 bp from an MRPL43 splice donor site, may alter the splice site. NC_000010.11:g.100978794A > G, in the isoform with a long 3′-UTR, may strengthen the binding affinity of the microRNA. Individuals with the TGG haplotype at these three variants had higher levels of mRNA expression and ribosome occupancy than individuals with the GCA haplotype but lower protein levels, producing the flipped effect throughout the expression process.

Conclusions

These findings suggest that multiple functional variants in a linkage exert their regulatory functions at different points in the gene expression process, producing a complexity of single eQTLs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01057-7.

Regulation of gene expression via translational buffering

Translation of an mRNA represents a critical step during the expression of protein-coding genes. As mechanisms governing post-transcriptional regulation of gene expression are progressively unveiled, it is becoming apparent that transcriptional programs are not fully reflected in the proteome. Herein, we highlight a previously underappreciated post-transcriptional mode of regulation of gene expression termed translational buffering. In principle, translational buffering opposes the impact of alterations in mRNA levels on the proteome. We further describe three types of translational buffering: compensation, which maintains protein levels e.g. across species or individuals; equilibration, which retains pathway stoichiometry; and offsetting, which acts as a reversible mechanism that maintains the levels of selected subsets of proteins constant despite genetic alteration and/or stress-induced changes in corresponding mRNA levels. While mechanisms underlying compensation and equilibration have been reviewed elsewhere, the principal focus of this review is on the less-well understood mechanism of translational offsetting. Finally, we discuss potential roles of translational buffering in homeostasis and disease.

Perspectives on Allele-Specific Expression

Diploidy has profound implications for population genetics and susceptibility to genetic diseases. Although two copies are present for most genes in the human genome, they are not necessarily both active or active at the same level in a given individual. Genomic imprinting, resulting in exclusive or biased expression in favor of the allele of paternal or maternal origin, is now believed to affect hundreds of human genes. A far greater number of genes display unequal expression of gene copies due to cis-acting genetic variants that perturb gene expression. The availability of data generated by RNA sequencing applied to large numbers of individuals and tissue types has generated unprecedented opportunities to assess the contribution of genetic variation to allelic imbalance in gene expression. Here we review the insights gained through the analysis of these data about the extent of the genetic contribution to allelic expression imbalance, the tools and statistical models for gene expression imbalance, and what the results obtained reveal about the contribution of genetic variants that alter gene expression to complex human diseases and phenotypes.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

Isoform-resolved correlation analysis between mRNA abundance regulation and protein level degradation

Profiling of biological relationships between different molecular layers dissects regulatory mechanisms that ultimately determine cellular function. To thoroughly assess the role of protein post‐translational turnover, we devised a strategy combining pulse stable isotope‐labeled amino acids in cells (pSILAC), data‐independent acquisition mass spectrometry (DIA‐MS), and a novel data analysis framework that resolves protein degradation rate on the level of mRNA alternative splicing isoforms and isoform groups. We demonstrated our approach by the genome‐wide correlation analysis between mRNA amounts and protein degradation across different strains of HeLa cells that harbor a high grade of gene dosage variation. The dataset revealed that specific biological processes, cellular organelles, spatial compartments of organelles, and individual protein isoforms of the same genes could have distinctive degradation rate. The protein degradation diversity thus dissects the corresponding buffering or concerting protein turnover control across cancer cell lines. The data further indicate that specific mRNA splicing events such as intron retention significantly impact the protein abundance levels. Our findings support the tight association between transcriptome variability and proteostasis and provide a methodological foundation for studying functional protein degradation.

Hybrid Incompatibilities and Transgressive Gene Expression Between Two Closely Related Subspecies of Drosophila

Genome-wide assays of expression between species and their hybrids have identified genes that become either over- or underexpressed relative to the parental species (i.e., transgressive). Transgressive expression in hybrids is of interest because it highlights possible changes in gene regulation linked to hybrid dysfunction. Previous studies in Drosophila that used long-diverged species pairs with complete or nearly complete isolation (i.e., full sterility and partial inviability of hybrids) and high-levels of genome misregulation have found correlations between expression and coding sequence divergence. The work highlighted the possible effects of directional selection driving sequence divergence and transgressive expression. Whether the same is true for taxa at early stages of divergence that have only achieved partial isolation remains untested. Here, we reanalyze previously published genome expression data and available genome sequence reads from a pair of partially isolated subspecies of Drosophila to compare expression and sequence divergence. We find a significant correlation in rates of expression and sequence evolution, but no support for directional selection driving transgressive expression in hybrids. We find that most transgressive genes in hybrids show no differential expression between parental subspecies and used SNP data to explore the role of stabilizing selection through compensatory mutations. We also examine possible misregulation through cascade effects that could be driven by interacting gene networks or co-option of off-target cis-regulatory elements.

Transcriptome and translatome co-evolution in mammals

Gene expression programs define shared and species-specific phenotypes, but their evolution remains largely uncharacterized beyond the transcriptome layer¹. Here we report an analysis of the co-evolution of translatomes and transcriptomes using ribosome-profling and matched RNA-sequencing data for three organs (brain, liver and testis) in fve mammals (human, macaque, mouse, opossum and platypus) and a bird (chicken). Our within-species analyses reveal that translational regulation is widespread in the diferent organs, in particular across the spermatogenic cell types of the testis. The between-species divergence in gene expression is around 20% lower at the translatome layer than at the transcriptome layer owing to extensive buffering between the expression layers, which especially preserved old, essential and housekeeping genes. Translational upregulation specifcally counterbalanced global dosage reductions during the evolution of sex chromosomes and the efects of meiotic sex-chromosome inactivation during spermatogenesis. Despite the overall prevalence of bufering, some genes evolved faster at the translatome layer—potentially indicating adaptive changes in expression; testis tissue shows the highest fraction of such genes. Further analyses incorporating mass spectrometry proteomics data establish that the co-evolution of transcriptomes and translatomes is refected at the proteome layer. Together, our work uncovers co-evolutionary patterns and associated selective forces across the expression layers, and provides a resource for understanding their interplay in mammalian organs.

A thermostability perspective on enhancing physicochemical and cytological characteristics of octacalcium phosphate by doping iron and strontium

Investigation of thermostability will lead the groundbreaking of unraveling the mechanism of influence of ion-doping on the properties of calcium phosphates. In this work, octacalcium phosphate (OCP), a metastable precursor of biological apatite, was used as a stability model for doping ions (Fe³⁺ and Sr²⁺) with different ionic charges and radii. After treated under hot air at different temperatures (110–200 °C), the phase, morphology, structure, physicochemical properties, protein affinity, ions release, and cytological responses of the ion-doped OCPs were investigated comparatively. The results showed that the collapse of OCP crystals gradually occurred, accompanying with the dehydration of hydrated layers and the disintegration of plate-like crystals as the temperature increased. The collapsed crystals still retained the typical properties of OCP and the potential of conversion into hydroxyapatite. Compared to the undoped OCP, Fe-OCP, and Sr-OCP had lower and higher thermostability respectively, leading to different material surface properties and ions release. The adjusted thermostability of Fe-OCP and Sr-OCP significantly enhanced the adsorption of proteins (BSA and LSZ) and the cytological behavior (adhesion, spreading, proliferation, and osteogenic differentiation) of bone marrow mesenchymal stem cells to a varying extent under the synergistic effects of corresponding surface characteristics and early active ions release. This work paves the way for understanding the modification mechanism of calcium phosphates utilizing ion doping strategy and developing bioactive OCP-based materials for tissue repair.

•

OCP was used as a stability model for doping ions with different charges and radii.

•

Collapse of OCP crystals occurred with structural dehydration after heat treatment.

•

Fe and Sr doping altered the thermostability of OCP crystals in an opposite way.

•

The thermostable difference affected the surface properties and ion release of OCP.

•

Active surface and ion release of OCP synergistically mediated its biocompatibility.

Grp94 Regulates the Recruitment of Aneural AChR Clusters for the Assembly of Postsynaptic Specializations by Modulating ADF/Cofilin Activity and Turnover

Temperature is a physiological factor that affects neuronal growth and synaptic homeostasis at the invertebrate neuromuscular junctions (NMJs); however, whether temperature stress could also regulate the structure and function of the vertebrate NMJs remains unclear. In this study, we use Xenopus laevis primary cultures as a vertebrate model system for investigating the involvement of heat shock protein 90 (HSP90) family of stress proteins in NMJ development. First, cold temperature treatment or HSP90 inhibition attenuates the formation of aneural acetylcholine receptor (AChR) clusters, but increases their stability after they are formed, in cultured muscles. HSP90 inhibition specifically affects the stability of aneural AChR clusters and their associated intracellular scaffolding protein rapsyn, instead of causing a global change in cell metabolism and protein expression in Xenopus muscle cultures. Upon synaptogenic stimulation, a specific HSP90 family member, glucose-regulated protein 94 (Grp94), modulates the phosphorylation and dynamic turnover of actin depolymerizing factor (ADF)/cofilin at aneural AChR clusters, leading to the recruitment of AChR molecules from aneural clusters to the assembly of agrin-induced postsynaptic specializations. Finally, postsynaptic Grp94 knock-down significantly inhibits nerve-induced AChR clustering and postsynaptic activity in nerve-muscle co-cultures as demonstrated by live-cell imaging and electrophysiological recording, respectively. Collectively, this study suggests that temperature-dependent alteration in Grp94 expression and activity inhibits the assembly of postsynaptic specializations through modulating ADF/cofilin phosphorylation and activity at aneural AChR clusters, which prevents AChR molecules from being recruited to the postsynaptic sites via actin-dependent vesicular trafficking, at developing vertebrate NMJs.

mRNAs, proteins and the emerging principles of gene expression control

Gene expression involves transcription, translation and the turnover of mRNAs and proteins. The degree to which protein abundances scale with mRNA levels and the implications in cases where this dependency breaks down remain an intensely debated topic. Here we review recent mRNA-protein correlation studies in the light of the quantitative parameters of the gene expression pathway, contextual confounders and buffering mechanisms. Although protein and mRNA levels typically show reasonable correlation, we describe how transcriptomics and proteomics provide useful non-redundant readouts. Integrating both types of data can reveal exciting biology and is an essential step in refining our understanding of the principles of gene expression control.

Unraveling cis and trans regulatory evolution during cotton domestication

Cis and trans regulatory divergence underlies phenotypic and evolutionary diversification. Relatively little is understood about the complexity of regulatory evolution accompanying crop domestication, particularly for polyploid plants. Here, we compare the fiber transcriptomes between wild and domesticated cotton (Gossypium hirsutum) and their reciprocal F1 hybrids, revealing genome-wide (~15%) and often compensatory cis and trans regulatory changes under divergence and domestication. The high level of trans evolution (54%-64%) observed is likely enabled by genomic redundancy following polyploidy. Our results reveal that regulatory variation is significantly associated with sequence evolution, inheritance of parental expression patterns, co-expression gene network properties, and genomic loci responsible for domestication traits. With respect to regulatory evolution, the two subgenomes of allotetraploid cotton are often uncoupled. Overall, our work underscores the complexity of regulatory evolution during fiber domestication and may facilitate new approaches for improving cotton and other polyploid plants.

Yeast genetic interaction screens in the age of CRISPR/Cas

The ease of performing both forward and reverse genetics in Saccharomyces cerevisiae, along with its stable haploid state and short generation times, has made this budding yeast the consummate model eukaryote for genetics. The major advantage of using budding yeast for reverse genetics is this organism's highly efficient homology-directed repair, allowing for precise genome editing simply by introducing DNA with homology to the chromosomal target. Although plasmid- and PCR-based genome editing tools are quite efficient, they depend on rare spontaneous DNA breaks near the target sequence. Consequently, they can generate only one genomic edit at a time, and the edit must be associated with a selectable marker. However, CRISPR/Cas technology is efficient enough to permit markerless and multiplexed edits in a single step. These features have made CRISPR/Cas popular for yeast strain engineering in synthetic biology and metabolic engineering applications, but it has not been widely employed for genetic screens. In this review, we critically examine different methods to generate multi-mutant strains in systematic genetic interaction screens and discuss the potential of CRISPR/Cas to supplement or improve on these methods.

Genomic evidence of Y chromosome microchimerism in the endometrium during endometriosis and in cases of infertility

BACKGROUND

Previous studies, which were primarily based on the fluorescent in-situ hybridisation (FISH) technique, revealed conflicting evidence regarding male foetal microchimerism in endometriosis. FISH is a relatively less sensitive technique, as it is performed on a small portion of the sample. Additionally, the probes used in the previous studies specifically detected centromeric and telomeric regions of Y chromosome, which are gene-sparse heterochromatised regions. In the present study, a panel of molecular biology tools such as qPCR, expression microarray, RNA-seq and qRT-PCR were employed to examine the Y chromosome microchimerism in the endometrium using secretory phase samples from fertile and infertile patients with severe (stage IV) ovarian endometriosis (OE) and without endometriosis.

METHODS

Microarray expression analysis followed by validation using RNA-seq and qRT-PCR experiments at the RNA levels and further validation at the DNA level by qPCR of target inserts for selected targets in eutopic endometrium samples obtained from control (CON) and stage IV ovarian endometriosis (OE), either from fertile (FCON and FOE; n = 30/each) or infertile (ICON and IOE; n = 30/each) women, were performed.

RESULTS

Six coding (AMELY, PCDH11, SRY, TGIF2LY, TSPY3, and USP9Y) and 10 non-coding (TTTY2, TTTY4C, TTTY5, TTTYY6, TTTY8, TTTY10, TTTY14, TTTY21, TTTY22, and TTTY23) genes exhibited a bimodal pattern of expression characterised by low expression in samples from fertile patients and high expression in samples from infertile patients. Seven coding MSY-linked genes (BAGE, CD24, EIF1AY, NLGN4Y, PRKY, VCY and ZFY) exhibited differential regulation in microarray analysis, and this change was validated by RNA-seq or qRT-PCR. DNA inserts for 7 genes in various samples were validated by qPCR. The prevalence and concentration of PCR-positive target inserts for BAGE, PRKY, TTTY9A and ZFY displayed higher values in the fertile, control (FCON) patients compared with the fertile, endometriosis patients (FOE).

CONCLUSION

Several coding and non-coding MSY-linked genes displayed microchimerism as evidenced by the presence of their respective DNA inserts, along with their differential transcript expression, in the endometrium during endometriosis and in cases of infertility.

CRISPR/Cas9-induced monoallelic mutations in the cytosolic AGPase large subunit gene APL2 induce the ectopic expression of APL2 and the corresponding small subunit gene APS2b in rice leaves

The first committed step in the endosperm starch biosynthetic pathway is catalyzed by the cytosolic glucose-1-phosphate adenylyl transferase (AGPase) comprising large and small subunits encoded by the OsAPL2 and OsAPS2b genes, respectively. OsAPL2 is expressed solely in the endosperm so we hypothesized that mutating this gene would block starch biosynthesis in the endosperm without affecting the leaves. We used CRISPR/Cas9 to create two heterozygous mutants, one with a severely truncated and nonfunctional AGPase and the other with a C-terminal structural modification causing a partial loss of activity. Unexpectedly, we observed starch depletion in the leaves of both mutants and a corresponding increase in the level of soluble sugars. This reflected the unanticipated expression of both OsAPL2 and OsAPS2b in the leaves, generating a complete ectopic AGPase in the leaf cytosol, and a corresponding decrease in the expression of the plastidial small subunit OsAPS2a that was only partially complemented by an increase in the expression of OsAPS1. The new cytosolic AGPase was not sufficient to compensate for the loss of plastidial AGPase, most likely because there is no wider starch biosynthesis pathway in the leaf cytosol and because pathway intermediates are not shuttled between the two compartments.

Genetics of trans-regulatory variation in gene expression

Heritable variation in gene expression forms a crucial bridge between genomic variation and the biology of many traits. However, most expression quantitative trait loci (eQTLs) remain unidentified. We mapped eQTLs by transcriptome sequencing in 1012 yeast segregants. The resulting eQTLs accounted for over 70% of the heritability of mRNA levels, allowing comprehensive dissection of regulatory variation. Most genes had multiple eQTLs. Most expression variation arose from trans-acting eQTLs distant from their target genes. Nearly all trans-eQTLs clustered at 102 hotspot locations, some of which influenced the expression of thousands of genes. Fine-mapped hotspot regions were enriched for transcription factor genes. While most genes had a local eQTL, most of these had no detectable effects on the expression of other genes in trans. Hundreds of non-additive genetic interactions accounted for small fractions of expression variation. These results reveal the complexity of genetic influences on transcriptome variation in unprecedented depth and detail.

The Evolution of Gene Expression in cis and trans

There is abundant variation in gene expression between individuals, populations, and species. The evolution of gene regulation and expression within and between species is thought to frequently contribute to adaptation. Yet considerable evidence suggests that the primary evolutionary force acting on variation in gene expression is stabilizing selection. We review here the results of recent studies characterizing the evolution of gene expression occurring in cis (via linked polymorphisms) or in trans (through diffusible products of other genes) and their contribution to adaptation and response to the environment. We review the evidence for buffering of variation in gene expression at the level of both transcription and translation, and the possible mechanisms for this buffering. Lastly, we summarize unresolved questions about the evolution of gene regulation.

Post-translational buffering leads to convergent protein expression levels between primates

BACKGROUND

Differences in gene regulation between human and closely related species influence phenotypes that are distinctly human. While gene regulation is a multi-step process, the majority of research concerning divergence in gene regulation among primates has focused on transcription.

RESULTS

To gain a comprehensive view of gene regulation, we surveyed genome-wide ribosome occupancy, which reflects levels of protein translation, in lymphoblastoid cell lines derived from human, chimpanzee, and rhesus macaque. We further integrated messenger RNA and protein level measurements collected from matching cell lines. We find that, in addition to transcriptional regulation, the major factor determining protein level divergence between human and closely related species is post-translational buffering. Inter-species divergence in transcription is generally propagated to the level of protein translation. In contrast, gene expression divergence is often attenuated post-translationally, potentially mediated through post-translational modifications.

CONCLUSIONS

Results from our analysis indicate that post-translational buffering is a conserved mechanism that led to relaxation of selective constraint on transcript levels in humans.

Characterization of a low expression haplotype in canine glutathione S-transferase (GSTT1) and its prevalence in golden retrievers

Glutathione S-transferase-theta (GSTT1) is a carcinogen detoxification enzyme, and low activity variants are associated with lymphoma in humans. We recently found a variant in the 3' untranslated region (UTR) of canine GSTT1, *101_102insT, which was predicted to change miRNA binding and was found in 5 of 17 golden retriever (GR) dogs with lymphoma but none of 14 healthy GRs. The aim of this study was to determine whether this variant led to decreased GSTT1 expression and was a discernible risk factor for lymphoma within the GR breed. On resequencing, *101_102insT appeared to be in complete linkage disequilibrium with 3 additional 3'UTR variants, leading to the inferred haplotype *3T>C; *101_102insT; *190C>A; *203T>C. In canine livers that were heterozygous for this variant haplotype, GSTT1 protein expression was significantly lower compared to the reference haplotype (densitometry .40 vs .64, P = .022), and GSTT1 transcript levels by qPCR were also significantly lower (fold difference .52, P = .012), without evidence of substantial allelic expression imbalance. The variant haplotype led to >50% decrease in expression in vitro (.31 ± .07 vs .64 ± .19; P = .019). We found no significant difference in minor allele frequencies between 71 GR dogs with lymphoma (MAF .162) and 33 healthy age-matched controls (MAF .136, P = .69). Our results indicate that the variant GSTT1 3'UTR haplotype containing *101_102insT reduces gene expression, which could lead to impaired carcinogen detoxification, but was not a detectable risk factor for lymphoma in GR dogs.

On the Dependency of Cellular Protein Levels on mRNA Abundance

The question of how genomic information is expressed to determine phenotypes is of central importance for basic and translational life science research and has been studied by transcriptomic and proteomic profiling. Here, we review the relationship between protein and mRNA levels under various scenarios, such as steady state, long-term state changes, and short-term adaptation, demonstrating the complexity of gene expression regulation, especially during dynamic transitions. The spatial and temporal variations of mRNAs, as well as the local availability of resources for protein biosynthesis, strongly influence the relationship between protein levels and their coding transcripts. We further discuss the buffering of mRNA fluctuations at the level of protein concentrations. We conclude that transcript levels by themselves are not sufficient to predict protein levels in many scenarios and to thus explain genotype-phenotype relationships and that high-quality data quantifying different levels of gene expression are indispensable for the complete understanding of biological processes.

Gene Regulation and Speciation

Understanding the genetic architecture of speciation is a major goal in evolutionary biology. Hybrid dysfunction is thought to arise most commonly through negative interactions between alleles at two or more loci. Divergence between interacting regulatory elements that affect gene expression (i.e., regulatory divergence) may be a common route for these negative interactions to arise. We review here how regulatory divergence between species can result in hybrid dysfunction, including recent theoretical support for this model. We then discuss the empirical evidence for regulatory divergence between species and evaluate evidence for misregulation as a source of hybrid dysfunction. Finally, we review unresolved questions in gene regulation as it pertains to speciation and point to areas that could benefit from future research.

Natural gene expression variation studies in yeast

The rise of sequence information across different yeast species and strains is driving an increasing number of studies in the emerging field of genomics to associate polymorphic variants, mRNA abundance and phenotypic differences between individuals. Here, we gathered evidence from recent studies covering several layers that define the genotype-phenotype gap, such as mRNA abundance, allele-specific expression and translation efficiency to demonstrate how genetic variants co-evolve and define an individual's genome. Moreover, we exposed several antecedents where inter- and intra-specific studies led to opposite conclusions, probably owing to genetic divergence. Future studies in this area will benefit from the access to a massive array of well-annotated genomes and new sequencing technologies, which will allow the fine breakdown of the complex layers that delineate the genotype-phenotype map. Copyright © 2016 John Wiley & Sons, Ltd.

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.

Buffering of Genetic Regulatory Networks in Drosophila melanogaster

Regulatory variation in gene expression can be described by cis- and trans-genetic components. Here we used RNA-seq data from a population panel of Drosophila melanogaster test crosses to compare allelic imbalance (AI) in female head tissue between mated and virgin flies, an environmental change known to affect transcription. Indeed, 3048 exons (1610 genes) are differentially expressed in this study. A Bayesian model for AI, with an intersection test, controls type I error. There are ∼200 genes with AI exclusively in mated or virgin flies, indicating an environmental component of expression regulation. On average 34% of genes within a cross and 54% of all genes show evidence for genetic regulation of transcription. Nearly all differentially regulated genes are affected in cis, with an average of 63% of expression variation explained by the cis-effects. Trans-effects explain 8% of the variance in AI on average and the interaction between cis and trans explains an average of 11% of the total variance in AI. In both environments cis- and trans-effects are compensatory in their overall effect, with a negative association between cis- and trans-effects in 85% of the exons examined. We hypothesize that the gene expression level perturbed by cis-regulatory mutations is compensated through trans-regulatory mechanisms, e.g., trans and cis by trans-factors buffering cis-mutations. In addition, when AI is detected in both environments, cis-mated, cis-virgin, and trans-mated-trans-virgin estimates are highly concordant with 99% of all exons positively correlated with a median correlation of 0.83 for cis and 0.95 for trans We conclude that the gene regulatory networks (GRNs) are robust and that trans-buffering explains robustness.

Splicing-Mediated Autoregulation Modulates Rpl22p Expression in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, splicing is critical for expression of ribosomal protein genes (RPGs), which are among the most highly expressed genes and are tightly regulated according to growth and environmental conditions. However, knowledge of the precise mechanisms by which RPG pre-mRNA splicing is regulated on a gene-by-gene basis is lacking. Here we show that Rpl22p has an extraribosomal role in the inhibition of splicing of the RPL22B pre-mRNA transcript. A stem loop secondary structure within the intron is necessary for pre-mRNA binding by Rpl22p in vivo and splicing inhibition in vivo and in vitro and can rescue splicing inhibition in vitro when added in trans to splicing reactions. Splicing inhibition by Rpl22p may be partly attributed to the reduction of co-transcriptional U1 snRNP recruitment to the pre-mRNA at the RPL22B locus. We further demonstrate that the inhibition of RPL22B pre-mRNA splicing contributes to the down-regulation of mature transcript during specific stress conditions, and provide evidence hinting at a regulatory role for this mechanism in conditions of suppressed ribosome biogenesis. These results demonstrate an autoregulatory mechanism that fine-tunes the expression of the Rpl22 protein and by extension Rpl22p paralog composition according to the cellular demands for ribosome biogenesis.

Ribosomes are responsible for protein production in all living cells, serving as the grounds for the translation of genetic information from RNA to protein. Given the vital role of the ribosome in protein assembly, ribosome components are highly expressed and are subject to tight regulation. Some ribosomal proteins are also known to engage in extra-ribosomal activities. In our study, we demonstrate that the ribosomal protein Rpl22p is able to regulate its own expression by inhibiting the processing of its own RNA transcript, leading to degradation of the RNA. We also show that this self-imposed regulation plays a role in limiting RPL22 transcript levels in specific stress conditions. We suggest that this mechanism may impact the composition of ribosomes by influencing the availability of the Rpl22p paralogs.

Drought responsive gene expression regulatory divergence between upland and lowland ecotypes of a perennial C4 grass

Climatic adaptation is an example of a genotype-by-environment interaction (G×E) of fitness. Selection upon gene expression regulatory variation can contribute to adaptive phenotypic diversity; however, surprisingly few studies have examined how genome-wide patterns of gene expression G×E are manifested in response to environmental stress and other selective agents that cause climatic adaptation. Here, we characterize drought-responsive expression divergence between upland (drought-adapted) and lowland (mesic) ecotypes of the perennial C4 grass,Panicum hallii, in natural field conditions. Overall, we find that cis-regulatory elements contributed to gene expression divergence across 47% of genes, 7.2% of which exhibit drought-responsive G×E. While less well-represented, we observe 1294 genes (7.8%) with transeffects.Trans-by-environment interactions are weaker and much less common than cis G×E, occurring in only 0.7% oft rans-regulated genes. Finally, gene expression heterosis is highly enriched in expression phenotypes with significant G×E. As such, modes of inheritance that drive heterosis, such as dominance or overdominance, may be common among G×E genes. Interestingly, motifs specific to drought-responsive transcription factors are highly enriched in the promoters of genes exhibiting G×E and transregulation, indicating that expression G×E and heterosis may result from the evolution of transcription factors or their binding sites.P. hallii serves as the genomic model for its close relative and emerging biofuel crop, switchgrass (Panicum virgatum). Accordingly, the results here not only aid in the discovery of the genetic mechanisms that underlie local adaptation but also provide a foundation to improve switchgrass yield under water-limited conditions.