Chromatin density and splicing destiny: on the cross-talk between chromatin structure and splicing

Exploring the influence of DNA methylation and single nucleotide polymorphisms of the Myostatin gene on growth traits in the hybrid grouper (Epinephelus fuscoguttatus (female) × Epinephelus polyphekadion (male))

Investigations into the correlation between growth characteristics and DNA methylation levels, along with genetic variations, can provide fundamental insights to enhance growth performance in groupers. The Myostatin (mstn) gene plays a vital role in regulating skeletal muscle development and growth. This study scrutinized the DNA methylation levels of the mstn gene across hybrid groupers (E. fuscoguttatus (♀) × E. polyphekadion (♂)) and their parental species, to evaluate its impact on growth attributes in grouper fish. The nucleotide sequence of the mstn gene was directly sequenced in the hybrid grouper, exhibiting different growth performance to identify the single nucleotide polymorphisms (SNPs) of the mstn gene and explore their correlation with growth characteristics. The findings revealed no significant differences in global DNA methylation levels within muscle tissue among the hybrid grouper and parents. However, significant differences in DNA methylation sites were discovered between the hybrid grouper and E. polyphekadion at sites 824 and 1521 (located at exon 2 and intron 2, respectively), and between E. fuscoguttatus and E. polyphekadion at site 1521. These variations could potentially influence the mRNA expression of the mstn gene. The study also identified that SNP g.1003 T > C in exon 2 of the mstn gene was significantly associated with various growth traits including body weight, total length, body length, head length, caudal peduncle height, and body height (p < 0.01). Specimens with the TT genotype at site 1003 demonstrated superior growth performance compared to those with the TC genotype. Furthermore, microstructural analyses of muscle tissue showed that the average area and diameter of muscle fibers in TT genotype individuals were significantly greater than those in TC genotype individuals. Therefore, this research provides robust evidence linking the DNA methylation level and polymorphisms of the mstn gene with growth traits, which could be beneficial for grouper breeding programs.

Predictive model of transcriptional elongation control identifies trans regulatory factors from chromatin signatures

Promoter-proximal Polymerase II (Pol II) pausing is a key rate-limiting step for gene expression. DNA and RNA-binding trans-acting factors regulating the extent of pausing have been identified. However, we lack a quantitative model of how interactions of these factors determine pausing, therefore the relative importance of implicated factors is unknown. Moreover, previously unknown regulators might exist. Here we address this gap with a machine learning model that accurately predicts the extent of promoter-proximal Pol II pausing from large-scale genome and transcriptome binding maps and gene annotation and sequence composition features. We demonstrate high accuracy and generalizability of the model by validation on an independent cell line which reveals the model's cell line agnostic character. Model interpretation in light of prior knowledge about molecular functions of regulatory factors confirms the interconnection of pausing with other RNA processing steps. Harnessing underlying feature contributions, we assess the relative importance of each factor, quantify their predictive effects and systematically identify previously unknown regulators of pausing. We additionally identify 16 previously unknown 7SK ncRNA interacting RNA-binding proteins predictive of pausing. Our work provides a framework to further our understanding of the regulation of the critical early steps in transcriptional elongation.

Elucidation of the Landscape of Alternatively Spliced Genes and Features in the Dorsal Striatum of Aggressive/Aggression-Deprived Mice in the Model of Chronic Social Conflicts

Both aggressive and aggression-deprived (AD) individuals represent pathological cases extensively studied in psychiatry and substance abuse disciplines. We employed the animal model of chronic social conflicts curated in our laboratory for over 30 years. In the study, we pursued the task of evaluation of the key events in the dorsal striatum transcriptomes of aggression-experienced mice and AD species, as compared with the controls, using RNA-seq profiling. We evaluated the alternative splicing-mediated transcriptome dynamics based on the RNA-seq data. We confined our attention to the exon skipping (ES) events as the major AS type for animals. We report the concurrent posttranscriptional and posttranslational regulation of the ES events observed in the phosphorylation cycles (in phosphoproteins and their targets) in the neuron-specific genes of the striatum. Strikingly, we found that major neurospecific splicing factors (Nova1, Ptbp1, 2, Mbnl1, 2, and Sam68) related to the alternative splicing regulation of cAMP genes (Darpp-32, Grin1, Ptpn5, Ppp3ca, Pde10a, Prkaca, Psd95, and Adora1) are upregulated specifically in aggressive individuals as compared with the controls and specifically AD animals, assuming intense switching between isoforms in the cAMP-mediated (de)phosphorylation signaling cascade. We found that the coding alternative splicing events were mostly attributed to synaptic plasticity and neural development-related proteins, while the nonsense-mediated decay-associated splicing events are mostly attributed to the mRNA processing of genes, including the spliceosome and splicing factors. In addition, considering the gene families, the transporter (Slc) gene family manifested most of the ES events. We found out that the major molecular systems employing AS for their plasticity are the ‘spliceosome’, ‘chromatin rearrangement complex’, ‘synapse’, and ‘neural development/axonogenesis’ GO categories. Finally, we state that approximately 35% of the exon skipping variants in gene coding regions manifest the noncoding variants subject to nonsense-mediated decay, employed as a homeostasis-mediated expression regulation layer and often associated with the corresponding gene expression alteration.

Adverse maternal environment affects hippocampal HTR2c variant expression and epigenetic characteristics in mouse offspring

BACKGROUND

An adverse maternal environment (AME) predisposes progeny towards cognitive impairment in humans and mice. Cognitive impairment associates with hippocampal dysfunction. An important regulator of hippocampal function is the hippocampal serotonergic system. Dysregulation of hippocampal serotonin receptor 2c (HTR2c) expression is linked with cognitive impairment. HTR2c contains multiple mRNA variants and isoforms that are epigenetically regulated including DNA methylation, histone modifications, and small nucleolar RNA MBII-52. We tested the hypotheses that AME increases HTR2c variant expression and alters epigenetic modifications along the HTR2c gene locus.

METHODS

We create an AME through maternal Western diet and prenatal environmental stress in the mouse. We analyzed hippocampal HTR2c and variants' expression, DNA methylation and histone modifications along the gene locus, and MBII-52 levels in postnatal day 21 offspring.

RESULTS

AME significantly increased the expressions of total HTR2c and full-length variants (V201 and V202) concurrently with an altered epigenetic profile along the HTR2c gene locus in male offspring hippocampi. Moreover, increased full-length variants' expression in AME males was in line with increased MBII-52 levels.

CONCLUSIONS

AME affects male offspring hippocampal expression of HTR2c and full-length variants via epigenetic mechanisms. Altered hippocampal HTR2c expression may contribute to cognitive impairment seen in adult males in this model.

IMPACT

The key message of our article is that an adverse maternal environment increases expression of total HTR2c mRNA and protein, alters proportions of HTR2c mRNA variants, and impacts HTR2c epigenetic modifications in male offspring hippocampi relative to controls. Our findings add to the literature by providing the first report of altered HTR2c mRNA variant expression in association with altered epigenetic modifications in the hippocampus of offspring mice exposed to an adverse maternal environment. Our findings suggest that an adverse maternal environment affects the expression of genes previously determined to regulate cognitive function through an epigenetic mechanism in a sex-specific manner.

PARP1′s Involvement in RNA Polymerase II Elongation: Pausing and Releasing Regulation through the Integrator and Super Elongation Complex

RNA polymerase elongation along the gene body is tightly regulated to ensure proper transcription and alternative splicing events. Understanding the mechanism and factors critical in regulating the rate of RNA polymerase II elongation and processivity is clearly important. Recently we showed that PARP1, a well-known DNA repair protein, when bound to chromatin, regulates RNA polymerase II elongation. However, the mechanism by which it does so is not known. In the current study, we aimed to tease out how PARP1 regulates RNAPII elongation. We show, both in vivo and in vitro, that PARP1 binds directly to the Integrator subunit 3 (IntS3), a member of the elongation Integrator complex. The association between the two proteins is mediated via the C-terminal domain of PARP1 to the C-terminal domain of IntS3. Interestingly, the occupancy of IntS3 along two PARP1 target genes mimicked that of PARP1, suggesting a role in its recruitment/assembly of elongation factors. Indeed, the knockdown of PARP1 resulted in differential chromatin association and gene occupancy of IntS3 and other key elongation factors. Most of these PARP1-mediated effects were due to the physical presence of PARP1 rather than its PARylation activity. These studies argue that PARP1 controls the progressive RNAPII elongation complexes. In summary, we present a platform to begin to decipher PARP1′s role in recruiting/scaffolding elongation factors along the gene body regions during RNA polymerase II elongation and gene regulation.

Histone Marks-Dependent Effect on Alternative Splicing: New Perspectives for Targeted Splicing Modulation in Cancer?

Alternative splicing (AS) is a tightly regulated mechanism that generates the complex human proteome from a small number of genes. Cis-regulatory RNA motifs in exons and introns control AS, recruiting positive and negative trans-acting splicing regulators. At a higher level, chromatin affects splicing events. Growing evidence indicates that the popular histone code hypothesis can be extended to RNA-level processes, such as AS. In addition to nucleosome positioning, which can generate transcriptional barriers to shape the final splicing outcome, histone post-translational modifications can contribute to the detailed regulation of single exon inclusion/exclusion. A histone-based system can identify alternatively spliced chromatin stretches, affecting RNAPII elongation locally or recruiting splicing components via adaptor complexes. In tumor cells, several mechanisms trigger misregulated AS events and produce cancer-associated transcripts. On a genome-wide level, aberrant AS can be the consequence of dysfunctional epigenetic splicing code, including altered enrichment in histone post-translational modifications. This review describes the main findings related to the effect of histone modifications and variants on splicing outcome and how a dysfunctional epigenetic splicing code triggers aberrant AS in cancer. In addition, it highlights recent advances in programmable DNA-targeting technologies and their possible application for AS targeted epigenetic modulation.

The SWI/SNF subunit BRG1 affects alternative splicing by changing RNA binding factor interactions with nascent RNA

BRG1 and BRM are ATPase core subunits of the human SWI/SNF chromatin remodelling complexes mainly associated with transcriptional initiation. They also have a role in alternative splicing, which has been shown for BRM-containing SWI/SNF complexes at a few genes. Here, we have identified a subset of genes which harbour alternative exons that are affected by SWI/SNF ATPases by expressing the ATPases BRG1 and BRM in C33A cells, a BRG1- and BRM-deficient cell line, and analysed the effect on splicing by RNA sequencing. BRG1- and BRM-affected sub-sets of genes favouring both exon inclusion and exon skipping, with only a minor overlap between the ATPase. Some of the changes in alternative splicing induced by BRG1 and BRM expression did not require the ATPase activity. The BRG1-ATPase independent included exons displayed an exon signature of a high GC content. By investigating three genes with exons affected by the BRG-ATPase-deficient variant, we show that these exons accumulated phosphorylated RNA pol II CTD, both serine 2 and serine 5 phosphorylation, without an enrichment of the RNA polymerase II. The ATPases were recruited to the alternative exons, together with both core and signature subunits of SWI/SNF complexes, and promoted the binding of RNA binding factors to chromatin and RNA at the alternative exons. The interaction with the nascent RNP, however, did not reflect the association to chromatin. The hnRNPL, hnRNPU and SAM68 proteins associated with chromatin in cells expressing BRG1 and BRM wild type, but the binding of hnRNPU to the nascent RNP was excluded. This suggests that SWI/SNF can regulate alternative splicing by interacting with splicing-RNA binding factor and influence their binding to the nascent pre-mRNA particle.

Sperm Methylome Profiling Can Discern Fertility Levels in the Porcine Biomedical Model

A combined Genotyping By Sequencing (GBS) and methylated DNA immunoprecipitation (MeDIP) protocol was used to identify—in parallel—genetic variation (Genomic-Wide Association Studies (GWAS) and epigenetic differences of Differentially Methylated Regions (DMR) in the genome of spermatozoa from the porcine animal model. Breeding boars with good semen quality (n = 11) and specific and well-documented differences in fertility (farrowing rate, FR) and prolificacy (litter size, LS) (n = 7) in artificial insemination programs, using combined FR and LS, were categorized as High Fertile (HF, n = 4) or Low Fertile (LF, n = 3), and boars with Unknown Fertility (UF, n = 4) were tested for eventual epigenetical similarity with those fertility-proven. We identified 165,944 Single Nucleotide Polymorphisms (SNPs) that explained 14–15% of variance among selection lines. Between HF and LF individuals (n = 7, 4 HF and 3 LF), we identified 169 SNPs with p ≤ 0.00015, which explained 58% of the variance. For the epigenetic analyses, we considered fertility and period of ejaculate collection (late-summer and mid-autumn). Approximately three times more DMRs were observed in HF than in LF boars across these periods. Interestingly, UF boars were clearly clustered with one of the other HF or LF groups. The highest differences in DMRs between HF and LF experimental groups across the pig genome were located in the chr 3, 9, 13, and 16, with most DMRs being hypermethylated in LF boars. In both HF and LF boars, DMRs were mostly hypermethylated in late-summer compared to mid-autumn. Three overlaps were detected between SNPs (p ≤ 0.0005, n = 1318) and CpG sites within DMRs. In conclusion, fertility levels in breeding males including FR and LS can be discerned using methylome analyses. The findings in this biomedical animal model ought to be applied besides sire selection for andrological diagnosis of idiopathic sub/infertility.

Putative Epigenetic Biomarkers of Stress in Red Blood Cells of Chickens Reared Across Different Biomes

Production animals are constantly subjected to early adverse environmental conditions that influence the adult phenotype and produce epigenetic effects. CpG dinucleotide methylation in red blood cells (RBC) could be a useful epigenetic biomarker to identify animals subjected to chronic stress in the production environment. Here we compared a reduced fraction of the RBC methylome of chickens exposed to social isolation to non-exposed. These experiments were performed in two different locations: Brazil and Sweden. The aim was to identify stress-associated DNA methylation profiles in RBC across these populations, in spite of the variable conditions to which birds are exposed in each facility and their different lineages. Birds were increasingly exposed to a social isolation treatment, combined with food and water deprivation, at random periods of the day from weeks 1-4 after hatching. We then collected the RBC DNA from individuals and compared a reduced fraction of their methylome between the experimental groups using two bioinformatic approaches to identify differentially methylated regions (DMRs): one using fixed-size windows and another that preselected differential peaks with MACS2. Three levels of significance were used (P ≤ 0.05, P ≤ 0.005, and P ≤ 0.0005) to identify DMRs between experimental groups, which were then used for different analyses. With both of the approaches more DMRs reached the defined significance thresholds in BR individuals compared to SW. However, more DMRs had higher fold change values in SW compared to BR individuals. Interestingly, ChrZ was enriched above expectancy for the presence of DMRs. Additionally, when analyzing the locations of these DMRs in relation to the transcription starting site (TSS), we found three peaks with high DMR presence: 10 kb upstream, the TSS itself, and 20-40 kb downstream. Interestingly, these peaks had DMRs with a high presence (>50%) of specific transcription factor binding sites. Three overlapping DMRs were found between the BR and SW population using the most relaxed p-value (P ≤ 0.05). With the most stringent p-value (P ≤ 0.0005), we found 7 and 4 DMRs between treatments in the BR and SW populations, respectively. This study is the first approximation to identify epigenetic biomarkers of long-term exposure to stress in different lineages of production animals.

PARPs and ADP-ribosylation in RNA biology: from RNA expression and processing to protein translation and proteostasis

In this review, Kim et al. discuss the importance of PARP family members and ADPRylation in gene regulation, mRNA processing, and protein abundance.

ADP-ribosylation (ADPRylation) is a posttranslational modification of proteins discovered nearly six decades ago, but many important questions remain regarding its molecular functions and biological roles, as well as the activity of the ADP-ribose (ADPR) transferase enzymes (PARP family members) that catalyze it. Growing evidence indicates that PARP-mediated ADPRylation events are key regulators of the protein biosynthetic pathway, leading from rDNA transcription and ribosome biogenesis to mRNA synthesis, processing, and translation. In this review we describe the role of PARP proteins and ADPRylation in all facets of this pathway. PARP-1 and its enzymatic activity are key regulators of rDNA transcription, which is a critical step in ribosome biogenesis. An emerging role of PARPs in alternative splicing of mRNAs, as well as direct ADPRylation of mRNAs, highlight the role of PARP members in RNA processing. Furthermore, PARP activity, stimulated by cellular stresses, such as viral infections and ER stress, leads to the regulation of mRNA stability and protein synthesis through posttranscriptional mechanisms. Dysregulation of PARP activity in these processes can promote disease states. Collectively, these results highlight the importance of PARP family members and ADPRylation in gene regulation, mRNA processing, and protein abundance. Future studies in these areas will yield new insights into the fundamental mechanisms and a broader utility for PARP-targeted therapeutic agents.

Regulation of alternative splicing by p300-mediated acetylation of splicing factors

Splicing of precursor mRNA (pre-mRNA) is an important regulatory step in gene expression. Recent evidence points to a regulatory role of chromatin-related proteins in alternative splicing regulation. Using an unbiased approach, we have identified the acetyltransferase p300 as a key chromatin-related regulator of alternative splicing. p300 promotes genome-wide exon inclusion in both a transcription-dependent and -independent manner. Using CD44 as a paradigm, we found that p300 regulates alternative splicing by modulating the binding of splicing factors to pre-mRNA. Using a tethering strategy, we found that binding of p300 to the CD44 promoter region promotes CD44v exon inclusion independently of RNAPII transcriptional elongation rate. Promoter-bound p300 regulates alternative splicing by acetylating splicing factors, leading to exclusion of hnRNP M from CD44 pre-mRNA and activation of Sam68. p300-mediated CD44 alternative splicing reduces cell motility and promotes epithelial features. Our findings reveal a chromatin-related mechanism of alternative splicing regulation and demonstrate its impact on cellular function.

RNA Polymerase II Phosphorylated on CTD Serine 5 Interacts with the Spliceosome during Co-transcriptional Splicing

The highly intronic nature of protein coding genes in mammals necessitates a co-transcriptional splicing mechanism as revealed by mNET-seq analysis. Immunoprecipitation of MNase-digested chromatin with antibodies against RNA polymerase II (Pol II) shows that active spliceosomes (both snRNA and proteins) are complexed to Pol II S5P CTD during elongation and co-transcriptional splicing. Notably, elongating Pol II-spliceosome complexes form strong interactions with nascent transcripts, resulting in footprints of approximately 60 nucleotides. Also, splicing intermediates formed by cleavage at the 5′ splice site are associated with nearly all spliced exons. These spliceosome-bound intermediates are frequently ligated to upstream exons, implying a sequential, constitutive, and U12-dependent splicing process. Finally, lack of detectable spliced products connected to the Pol II active site in human HeLa or murine lymphoid cells suggests that splicing does not occur immediately following 3′ splice site synthesis. Our results imply that most mammalian splicing requires exon definition for completion.

•

S5P CTD Pol II associates with the catalytic spliceosome

•

Elongating Pol II complexes protect about 60 newly synthesized nucleotides

•

Co-transcriptional splicing associated with dominant 5′ ss intermediate

•

U12-dependent introns are sequentially spliced in association with Pol II

A slow transcription rate causes embryonic lethality and perturbs kinetic coupling of neuronal genes

The rate of RNA polymerase II (RNAPII) elongation has an important role in the control of alternative splicing (AS); however, the in vivo consequences of an altered elongation rate are unknown. Here, we generated mouse embryonic stem cells (ESCs) knocked in for a slow elongating form of RNAPII We show that a reduced transcriptional elongation rate results in early embryonic lethality in mice. Focusing on neuronal differentiation as a model, we observed that slow elongation impairs development of the neural lineage from ESCs, which is accompanied by changes in AS and in gene expression along this pathway. In particular, we found a crucial role for RNAPII elongation rate in transcription and splicing of long neuronal genes involved in synapse signaling. The impact of the kinetic coupling of RNAPII elongation rate with AS is greater in ESC-differentiated neurons than in pluripotent cells. Our results demonstrate the requirement for an appropriate transcriptional elongation rate to ensure proper gene expression and to regulate AS during development.

Coupling of PARP1-mediated chromatin structural changes to transcriptional RNA polymerase II elongation and cotranscriptional splicing

Background

Recently, we showed that PARP1 is involved in cotranscriptional splicing, possibly by bridging chromatin to RNA and recruiting splicing factors. It also can influence alternative splicing decisions through the regulation of RNAPII elongation. In this study, we investigated the effect of PARP1-mediated chromatin changes on RNAPII movement, during transcription and alternative splicing.

Results

We show that RNAPII pauses at PARP1–chromatin structures within the gene body. Knockdown of PARP1 abolishes this RNAPII pausing, suggesting that PARP1 may regulate RNAPII elongation. Additionally, PARP1 alters nucleosome deposition and histone post-translational modifications at specific exon–intron boundaries, thereby affecting RNAPII movement. Lastly, genome-wide analyses confirmed that PARP1 influences changes in RNAPII elongation by either reducing or increasing the rate of RNAPII elongation depending on the chromatin context.

Conclusions

These studies suggest a context-specific effect of PARP1–chromatin binding on RNA polymerase movement and provide a platform to delineate PARP1’s role in RNA biogenesis and processing.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0261-1) contains supplementary material, which is available to authorized users.

The chromatin remodeler ZmCHB101 impacts alternative splicing contexts in response to osmotic stress

Maize SWI3-type chromatin remodeler impacts alternative splicing contexts in response to osmotic stress by altering nucleosome density and affecting transcriptional elongation rate. Alternative splicing (AS) is commonly found in higher eukaryotes and is an important posttranscriptional regulatory mechanism to generate transcript diversity. AS has been widely accepted as playing essential roles in different biological processes including growth, development, signal transduction and responses to biotic and abiotic stresses in plants. However, whether and how chromatin remodeling complex functions in AS in plant under osmotic stress remains unknown. Here, we show that a maize SWI3D protein, ZmCHB101, impacts AS contexts in response to osmotic stress. Genome-wide analysis of mRNA contexts in response to osmotic stress using ZmCHB101-RNAi lines reveals that ZmCHB101 impacts alternative splicing contexts of a subset of osmotic stress-responsive genes. Intriguingly, ZmCHB101-mediated regulation of gene expression and AS is largely uncoupled, pointing to diverse molecular functions of ZmCHB101 in transcriptional and posttranscriptional regulation. We further found ZmCHB101 impacts the alternative splicing contexts by influencing alteration of chromatin and histone modification status as well as transcriptional elongation rates mediated by RNA polymerase II. Taken together, our findings suggest a novel insight of how plant chromatin remodeling complex impacts AS under osmotic stress .

Genome-wide mapping reveals conserved and diverged R-loop activities in the unusual genetic landscape of the African trypanosome genome

R-loops are stable RNA-DNA hybrids that have been implicated in transcription initiation and termination, as well as in telomere maintenance, chromatin formation, and genome replication and instability. RNA Polymerase (Pol) II transcription in the protozoan parasite Trypanosoma brucei is highly unusual: virtually all genes are co-transcribed from multigene transcription units, with mRNAs generated by linked trans-splicing and polyadenylation, and transcription initiation sites display no conserved promoter motifs. Here, we describe the genome-wide distribution of R-loops in wild type mammal-infective T. brucei and in mutants lacking RNase H1, revealing both conserved and diverged functions. Conserved localization was found at centromeres, rRNA genes and retrotransposon-associated genes. RNA Pol II transcription initiation sites also displayed R-loops, suggesting a broadly conserved role despite the lack of promoter conservation or transcription initiation regulation. However, the most abundant sites of R-loop enrichment were within the regions between coding sequences of the multigene transcription units, where the hybrids coincide with sites of polyadenylation and nucleosome-depletion. Thus, instead of functioning in transcription termination the most widespread localization of R-loops in T. brucei suggests a novel correlation with pre-mRNA processing. Finally, we find little evidence for correlation between R-loop localization and mapped sites of DNA replication initiation.

Epigenetic regulation of alternative splicing

Alternative splicing (AS) serves as an additional regulatory process for gene expression after transcription, and it generates distinct mRNA species, and even noncoding RNAs (ncRNAs), from one primary transcript. Generally, AS can be coupled with transcription and subjected to epigenetic regulation, such as DNA methylation and histone modifications. In addition, ncRNAs, especially long noncoding RNAs (lncRNAs), can be generated from AS and function as splicing factors ("interactors" or "hijackers") in AS. Recently, RNA modifications, such as the RNA N6-methyladenosine (m6A) modification, have been found to regulate AS. In this review, we summarize recent achievements related to the epigenetic regulation of AS.

The Krebs Cycle Connection: Reciprocal Influence Between Alternative Splicing Programs and Cell Metabolism

Alternative splicing is a pervasive mechanism that molds the transcriptome to meet cell and organism needs. However, how this layer of gene expression regulation is coordinated with other aspects of the cell metabolism is still largely undefined. Glucose is the main energy and carbon source of the cell. Not surprisingly, its metabolism is finely tuned to satisfy growth requirements and in response to nutrient availability. A number of studies have begun to unveil the connections between glucose metabolism and splicing programs. Alternative splicing modulates the ratio between M1 and M2 isoforms of pyruvate kinase in this way determining the choice between aerobic glycolysis and complete glucose oxidation in the Krebs cycle. Reciprocally, intermediates in the Krebs cycle may impact splicing programs at different levels by modulating the activity of 2-oxoglutarate-dependent oxidases. In this review we discuss the molecular mechanisms that coordinate alternative splicing programs with glucose metabolism, two aspects with profound implications in human diseases.

Human exonization through differential nucleosome occupancy

Nucleosomal modifications have been implicated in fundamental epigenetic regulation, whereas the roles of nucleosome binding in shaping changes through evolution remain to be addressed. Here we performed a comparative study to clarify the roles of nucleosome occupancy in exon origination. By profiling a high-resolution, cross-species mononucleosome landscape for mammalian tissues, we found nucleosome occupancy profiles are conserved across tissues and species. Further, through a phylogenetic approach, we found that the feature of differential nucleosome occupancy appears prior to the origination of new exons and, presumably, facilitates the origin of new exons by increasing the splice strength of the ancestral nonexonic regions through driving a local difference in GC content, which suggests the function of nucleosome binding in exonization.

Nucleosomal modifications have been implicated in fundamental epigenetic regulation, but the roles of nucleosome occupancy in shaping changes through evolution remain to be addressed. Here we present high-resolution nucleosome occupancy profiles for multiple tissues derived from human, macaque, tree shrew, mouse, and pig. Genome-wide comparison reveals conserved nucleosome occupancy profiles across both different species and tissue types. Notably, we found significantly higher levels of nucleosome occupancy in exons than in introns, a pattern correlated with the different exon–intron GC content. We then determined whether this biased occupancy may play roles in the origination of new exons through evolution, rather than being a downstream effect of exonization, through a comparative approach to sequentially trace the order of the exonization and biased nucleosome binding. By identifying recently evolved exons in human but not in macaque using matched RNA sequencing, we found that higher exonic nucleosome occupancy also existed in macaque regions orthologous to these exons. Presumably, such biased nucleosome occupancy facilitates the origination of new exons by increasing the splice strength of the ancestral nonexonic regions through driving a local difference in GC content. These data thus support a model that sites bound by nucleosomes are more likely to evolve into exons, which we term the “nucleosome-first” model.

From "Cellular" RNA to "Smart" RNA: Multiple Roles of RNA in Genome Stability and Beyond

Coding for proteins has been considered the main function of RNA since the "central dogma" of biology was proposed. The discovery of noncoding transcripts shed light on additional roles of RNA, ranging from the support of polypeptide synthesis, to the assembly of subnuclear structures, to gene expression modulation. Cellular RNA has therefore been recognized as a central player in often unanticipated biological processes, including genomic stability. This ever-expanding list of functions inspired us to think of RNA as a "smart" phone, which has replaced the older obsolete "cellular" phone. In this review, we summarize the last two decades of advances in research on the interface between RNA biology and genome stability. We start with an account of the emergence of noncoding RNA, and then we discuss the involvement of RNA in DNA damage signaling and repair, telomere maintenance, and genomic rearrangements. We continue with the depiction of single-molecule RNA detection techniques, and we conclude by illustrating the possibilities of RNA modulation in hopes of creating or improving new therapies. The widespread biological functions of RNA have made this molecule a reoccurring theme in basic and translational research, warranting it the transcendence from classically studied "cellular" RNA to "smart" RNA.

The histone variant H2A.Z promotes splicing of weak introns

In this study, Nissen et al. investigated the function of the highly conserved histone variant H2A.Z in pre-mRNA splicing using the intron-rich model yeast S. pombe. The findings suggest that H2A.Z occupancy promotes cotranscriptional splicing of suboptimal introns that may otherwise be discarded via proofreading ATPases.

Multiple lines of evidence implicate chromatin in the regulation of premessenger RNA (pre-mRNA) splicing. However, the influence of chromatin factors on cotranscriptional splice site usage remains unclear. Here we investigated the function of the highly conserved histone variant H2A.Z in pre-mRNA splicing using the intron-rich model yeast Schizosaccharomyces pombe. Using epistatic miniarray profiles (EMAPs) to survey the genetic interaction landscape of the Swr1 nucleosome remodeling complex, which deposits H2A.Z, we uncovered evidence for functional interactions with components of the spliceosome. In support of these genetic connections, splicing-specific microarrays show that H2A.Z and the Swr1 ATPase are required during temperature stress for the efficient splicing of a subset of introns. Notably, affected introns are enriched for H2A.Z occupancy and more likely to contain nonconsensus splice sites. To test the significance of the latter correlation, we mutated the splice sites in an affected intron to consensus and found that this suppressed the requirement for H2A.Z in splicing of that intron. These data suggest that H2A.Z occupancy promotes cotranscriptional splicing of suboptimal introns that may otherwise be discarded via proofreading ATPases. Consistent with this model, we show that overexpression of splicing ATPase Prp16 suppresses both the growth and splicing defects seen in the absence of H2A.Z.

GT-rich promoters can drive RNA pol II transcription and deposition of H2A.Z in African trypanosomes

Genome‐wide transcription studies are revealing an increasing number of “dispersed promoters” that, unlike “focused promoters”, lack well‐conserved sequence motifs and tight regulation. Dispersed promoters are nevertheless marked by well‐defined chromatin structures, suggesting that specific sequence elements must exist in these unregulated promoters. Here, we have analyzed regions of transcription initiation in the eukaryotic parasite Trypanosoma brucei, in which RNA polymerase II transcription initiation occurs over broad regions without distinct promoter motifs and lacks regulation. Using a combination of site‐specific and genome‐wide assays, we identified GT‐rich promoters that can drive transcription and promote the targeted deposition of the histone variant H2A.Z in a genomic context‐dependent manner. In addition, upon mapping nucleosome occupancy at high resolution, we find nucleosome positioning to correlate with RNA pol II enrichment and gene expression, pointing to a role in RNA maturation. Nucleosome positioning may thus represent a previously unrecognized layer of gene regulation in trypanosomes. Our findings show that even highly dispersed, unregulated promoters contain specific DNA elements that are able to induce transcription and changes in chromatin structure.

RNA splicing in human disease and in the clinic

Defects at the level of the pre-mRNA splicing process represent a major cause of human disease. Approximately 15-50% of all human disease mutations have been shown to alter functioning of basic and auxiliary splicing elements. These elements are required to ensure proper processing of pre-mRNA splicing molecules, with their disruption leading to misprocessing of the pre-mRNA molecule and disease. The splicing process is a complex process, with much still to be uncovered before we are able to accurately predict whether a reported genomic sequence variant (GV) represents a splicing-associated disease mutation or a harmless polymorphism. Furthermore, even when a mutation is correctly identified as affecting the splicing process, there still remains the difficulty of providing an exact evaluation of the potential impact on disease onset, severity and duration. In this review, we provide a brief overview of splicing diagnostic methodologies, from in silico bioinformatics approaches to wet lab in vitro and in vivo systems to evaluate splicing efficiencies. In particular, we provide an overview of how the latest developments in high-throughput sequencing can be applied to the clinic, and are already changing clinical approaches.

Modulation the alternative splicing of GLA (IVS4+919G>A) in Fabry disease

While a base substitution in intron 4 of GLA (IVS4+919G>A) that causes aberrant alternative splicing resulting in Fabry disease has been reported, its molecular mechanism remains unclear. Here we reported that upon IVS4+919G>A transversion, H3K36me3 was enriched across the alternatively spliced region. PSIP1, an adapter of H3K36me3, together with Hsp70 and NONO were recruited and formed a complex with SF2/ASF and SRp20, which further promoted GLA splicing. Amiloride, a splicing regulator in cancer cells, could reverse aberrant histone modification patterns and disrupt the association of splicing complex with GLA. It could also reverse aberrant GLA splicing in a PP1-dependant manner. Our findings revealed the alternative splicing mechanism of GLA (IVS4+919G>A), and a potential treatment for this specific genetic type of Fabry disease by amiloride in the future.

Transposon-mediated Generation of Cellular and Mouse Models of Splicing Mutations to Assess the Efficacy of snRNA-based Therapeutics

Disease-causing splicing mutations can be rescued by variants of the U1 small nuclear RNA (U1snRNAs). However, the evaluation of the efficacy and safety of modified U1snRNAs as therapeutic tools is limited by the availability of cellular and animal models specific for a given mutation. Hence, we exploited the hyperactive Sleeping Beauty transposon system (SB100X) to integrate human factor IX (hFIX) minigenes into genomic DNA in vitro and in vivo. We generated stable HEK293 cell lines and C57BL/6 mice harboring splicing-competent hFIX minigenes either wild type (SChFIX-wt) or mutated (SChFIXex5-2C). In both models the SChFIXex5-2C variant, found in patients affected by Hemophilia B, displayed an aberrant splicing pattern characterized by exon 5 skipping. This allowed us to test, for the first time in a genomic DNA context, the efficacy of the snRNA U1-fix9, delivered with an adeno-associated virus (AAV) vector. With this approach, we showed rescue of the correct splicing pattern of hFIX mRNA, leading to hFIX protein expression. These data validate the SB100X as a versatile tool to quickly generate models of human genetic mutations, to study their effect in a stable DNA context and to assess mutation-targeted therapeutic strategies.

A Broad Set of Chromatin Factors Influences Splicing

Several studies propose an influence of chromatin on pre-mRNA splicing, but it is still unclear how widespread and how direct this phenomenon is. We find here that when assembled in vivo, the U2 snRNP co-purifies with a subset of chromatin-proteins, including histones and remodeling complexes like SWI/SNF. Yet, an unbiased RNAi screen revealed that the outcome of splicing is influenced by a much larger variety of chromatin factors not all associating with the spliceosome. The availability of this broad range of chromatin factors impacting splicing further unveiled their very context specific effect, resulting in either inclusion or skipping, depending on the exon under scrutiny. Finally, a direct assessment of the impact of chromatin on splicing using an in vitro co-transcriptional splicing assay with pre-mRNAs transcribed from a nucleosomal template, demonstrated that chromatin impacts nascent pre-mRNP in their competence for splicing. Altogether, our data show that numerous chromatin factors associated or not with the spliceosome can affect the outcome of splicing, possibly as a function of the local chromatin environment that by default interferes with the efficiency of splicing.

Splicing is an RNA editing step allowing to produce multiple transcripts from a single gene. The gene itself is organized in chromatin, associating DNA and multiple proteins. Some proteins regulating the compaction of the chromatin also affect RNA splicing. Yet, it was unclear whether these chromatin proteins were exceptions or whether chromatin very generally affected the outcome of splicing. Here, we show that a subset of chromatin proteins is physically in interaction with the enzyme responsible for RNA splicing. In addition, several chromatin proteins not found directly associated with the splicing machinery were also able to influence RNA splicing, suggesting that chromatin compaction very globally plays a role in splicing. This finding was confirmed using the first in vitro assay combining transcription and splicing in the context of chromatin; this assay showed that assembling DNA with chromatin proteins influences the efficiency of splicing.

Involvement of PARP1 in the regulation of alternative splicing

Specialized chromatin structures such as nucleosomes with specific histone modifications decorate exons in eukaryotic genomes, suggesting a functional connection between chromatin organization and the regulation of pre-mRNA splicing. Through profiling the functional location of Poly (ADP) ribose polymerase, we observed that it is associated with the nucleosomes at exon/intron boundaries of specific genes, suggestive of a role for this enzyme in alternative splicing. Poly (ADP) ribose polymerase has previously been implicated in the PARylation of splicing factors as well as regulation of the histone modification H3K4me3, a mark critical for co-transcriptional splicing. In light of these studies, we hypothesized that interaction of the chromatin-modifying factor, Poly (ADP) ribose polymerase with nucleosomal structures at exon–intron boundaries, might regulate pre-mRNA splicing. Using genome-wide approaches validated by gene-specific assays, we show that depletion of PARP1 or inhibition of its PARylation activity results in changes in alternative splicing of a specific subset of genes. Furthermore, we observed that PARP1 bound to RNA, splicing factors and chromatin, suggesting that Poly (ADP) ribose polymerase serves as a gene regulatory hub to facilitate co-transcriptional splicing. These studies add another function to the multi-functional protein, Poly (ADP) ribose polymerase, and provide a platform for further investigation of this protein’s function in organizing chromatin during gene regulatory processes.

Regulation of alternative splicing through coupling with transcription and chromatin structure

Alternative precursor messenger RNA (pre-mRNA) splicing plays a pivotal role in the flow of genetic information from DNA to proteins by expanding the coding capacity of genomes. Regulation of alternative splicing is as important as regulation of transcription to determine cell- and tissue-specific features, normal cell functioning, and responses of eukaryotic cells to external cues. Its importance is confirmed by the evolutionary conservation and diversification of alternative splicing and the fact that its deregulation causes hereditary disease and cancer. This review discusses the multiple layers of cotranscriptional regulation of alternative splicing in which chromatin structure, DNA methylation, histone marks, and nucleosome positioning play a fundamental role in providing a dynamic scaffold for interactions between the splicing and transcription machineries. We focus on evidence for how the kinetics of RNA polymerase II (RNAPII) elongation and the recruitment of splicing factors and adaptor proteins to chromatin components act in coordination to regulate alternative splicing.

Trm9-Catalyzed tRNA Modifications Regulate Global Protein Expression by Codon-Biased Translation

Post-transcriptional modifications of transfer RNAs (tRNAs) have long been recognized to play crucial roles in regulating the rate and fidelity of translation. However, the extent to which they determine global protein production remains poorly understood. Here we use quantitative proteomics to show a direct link between wobble uridine 5-methoxycarbonylmethyl (mcm⁵) and 5-methoxy-carbonyl-methyl-2-thio (mcm⁵s²) modifications catalyzed by tRNA methyltransferase 9 (Trm9) in tRNA^Arg(UCU) and tRNA^Glu(UUC) and selective translation of proteins from genes enriched with their cognate codons. Controlling for bias in protein expression and alternations in mRNA expression, we find that loss of Trm9 selectively impairs expression of proteins from genes enriched with AGA and GAA codons under both normal and stress conditions. Moreover, we show that AGA and GAA codons occur with high frequency in clusters along the transcripts, which may play a role in modulating translation. Consistent with these results, proteins subject to enhanced ribosome pausing in yeast lacking mcm⁵U and mcm⁵s²U are more likely to be down-regulated and contain a larger number of AGA/GAA clusters. Together, these results suggest that Trm9-catalyzed tRNA modifications play a significant role in regulating protein expression within the cell.

Here we present evidence for a more complicated role for transfer RNAs (tRNAs) than as mere adapters that link the genetic code in messenger RNA (mRNA) to the amino acid sequence of a protein during translation. tRNAs have long been known to be modified with dozens of different chemical structures other than the 4 canonical ribonucleosides, though the role of these modifications in controlling translation is poorly understood. By quantifying the expression of thousands of proteins in the yeast S. cerevisiae, we identified a mechanistic link between modified ribonucleosides located at the wobble position of two tRNAs, tRNA^Arg(UCU) and tRNA^Glu(UUC), and the translation of proteins derived from genes enriched with codons read by these tRNAs: AGA and GAA. In cells lacking the enzyme that inserts these modifications, tRNA methyltransferase 9 (Trm9), we found a significant reduction in proteins from genes enriched with AGA and GAA codons and with runs of these codons. Also, mRNAs enriched with runs of AGA and GAA codons are subject to stalled translation on ribosomes in yeast lacking mcm⁵U and mcm⁵s²U. Together, these results reveal a distinct role for Trm9-catalyzed tRNA modifications in selectively regulating the expression of proteins enriched with AGA and GAA codons.

Global analysis of physical and functional RNA targets of hnRNP L reveals distinct sequence and epigenetic features of repressed and enhanced exons

HnRNP L is a ubiquitous splicing-regulatory protein that is critical for the development and function of mammalian T cells. Previous work has identified a few targets of hnRNP L-dependent alternative splicing in T cells and has described transcriptome-wide association of hnRNP L with RNA. However, a comprehensive analysis of the impact of hnRNP L on mRNA expression remains lacking. Here we use next-generation sequencing to identify transcriptome changes upon depletion of hnRNP L in a model T-cell line. We demonstrate that hnRNP L primarily regulates cassette-type alternative splicing, with minimal impact of hnRNP L depletion on transcript abundance, intron retention, or other modes of alternative splicing. Strikingly, we find that binding of hnRNP L within or flanking an exon largely correlates with exon repression by hnRNP L. In contrast, exons that are enhanced by hnRNP L generally lack proximal hnRNP L binding. Notably, these hnRNP L-enhanced exons share sequence and context features that correlate with poor nucleosome positioning, suggesting that hnRNP may enhance inclusion of a subset of exons via a cotranscriptional or epigenetic mechanism. Our data demonstrate that hnRNP L controls inclusion of a broad spectrum of alternative cassette exons in T cells and suggest both direct RNA regulation as well as indirect mechanisms sensitive to the epigenetic landscape.

Splicing Regulation: A Molecular Device to Enhance Cancer Cell Adaptation

Alternative splicing (AS) represents a major resource for eukaryotic cells to expand the coding potential of their genomes and to finely regulate gene expression in response to both intra- and extracellular cues. Cancer cells exploit the flexible nature of the mechanisms controlling AS in order to increase the functional diversity of their proteome. By altering the balance of splice isoforms encoded by human genes or by promoting the expression of aberrant oncogenic splice variants, cancer cells enhance their ability to adapt to the adverse growth conditions of the tumoral microenvironment. Herein, we will review the most relevant cancer-related splicing events and the underlying regulatory mechanisms allowing tumour cells to rapidly adapt to the harsh conditions they may face during the occurrence and development of cancer.

New insights into the expression and functions of the Kaposi's sarcoma-associated herpesvirus long noncoding PAN RNA

The Kaposi's sarcoma-associated herpesvirus (KSHV) is a clinically relevant pathogen associated with several human diseases that primarily affect immunocompromised individuals. KSHV encodes a noncoding polyadenylated nuclear (PAN) RNA that is essential for viral propagation and viral gene expression. PAN RNA is the most abundant viral transcript produced during lytic replication. The accumulation of PAN RNA depends on high levels of transcription driven by the Rta protein, a KSHV transcription factor necessary and sufficient for latent-to-lytic phase transition. In addition, KSHV uses several posttranscriptional mechanisms to stabilize PAN RNA. A cis-acting element, called the ENE, prevents PAN RNA decay by forming a triple helix with its poly(A) tail. The viral ORF57 and the cellular PABPC1 proteins further contribute to PAN RNA stability during lytic phase. PAN RNA functions are only beginning to be uncovered, but PAN RNA has been proposed to control gene expression by several different mechanisms. PAN RNA associates with the KSHV genome and may regulate gene expression by recruiting chromatin-modifying factors. Moreover, PAN RNA binds the viral latency-associated nuclear antigen (LANA) protein and decreases its repressive activity by sequestering it from the viral genome. Surprisingly, PAN RNA was found to associate with translating ribosomes, so this noncoding RNA may be additionally used to produce viral peptides. In this review, I highlight the mechanisms of PAN RNA accumulation and describe recent insights into potential functions of PAN RNA.

Role of six single nucleotide polymorphisms, risk factors in coronary disease, in OLR1 alternative splicing

The OLR1 gene encodes the oxidized low-density lipoprotein receptor (LOX-1), which is responsible for the cellular uptake of oxidized LDL (Ox-LDL), foam cell formation in atheroma plaques and atherosclerotic plaque rupture. Alternative splicing (AS) of OLR1 exon 5 generates two protein isoforms with antagonistic functions in Ox-LDL uptake. Previous work identified six single nucleotide polymorphisms (SNPs) in linkage disequilibrium that influence the inclusion levels of OLR1 exon 5 and correlate with the risk of cardiovascular disease. Here we use minigenes to recapitulate the effects of two allelic series (Low- and High-Risk) on OLR1 AS and identify one SNP in intron 4 (rs3736234) as the main contributor to the differences in exon 5 inclusion, while the other SNPs in the allelic series attenuate the drastic effects of this key SNP. Bioinformatic, proteomic, mutational and functional high-throughput analyses allowed us to define regulatory sequence motifs and identify SR protein family members (SRSF1, SRSF2) and HMGA1 as factors involved in the regulation of OLR1 AS. Our results suggest that antagonism between SRSF1 and SRSF2/HMGA1, and differential recognition of their regulatory motifs depending on the identity of the rs3736234 polymorphism, influence OLR1 exon 5 inclusion and the efficiency of Ox-LDL uptake, with potential implications for atherosclerosis and coronary disease.

Acetylation of intragenic histones on HPV16 correlates with enhanced HPV16 gene expression

We report that many histone modifications are unevenly distributed over the HPV16 genome in cervical cancer cells as well as in HPV16-immortalized keratinocytes. For example, H3K36me3 and H3K9Ac that are common in highly expressed cellular genes and over exons, were more common in the early than in the late region of the HPV16 genome. In contrast, H3K9me3, H4K20me3, H2BK5me1 and H4K16Ac were more frequent in the HPV16 late region. Furthermore, a region encompassing the HPV16 early polyadenylation signal pAE displayed high levels of histone H3 acetylation. Histone deacetylase (HDAC) inhibitors caused a 2- to 8-fold induction of HPV16 early and late mRNAs in cervical cancer cells and in immortalized keratinocytes, while at the same time increasing the levels of acetylated histones in the cells and on the HPV16 genome specifically. We concluded that increased histone acetylation on the HPV16 genome correlates with increased HPV16 gene expression.

Widespread exon skipping triggers degradation by nuclear RNA surveillance in fission yeast

Exon skipping is considered a principal mechanism by which eukaryotic cells expand their transcriptome and proteome repertoires, creating different splice variants with distinct cellular functions. Here we analyze RNA-seq data from 116 transcriptomes in fission yeast (Schizosaccharomyces pombe), covering multiple physiological conditions as well as transcriptional and RNA processing mutants. We applied brute-force algorithms to detect all possible exon-skipping events, which were widespread but rare compared to normal splicing events. Exon-skipping events increased in cells deficient for the nuclear exosome or the 5′-3′ exonuclease Dhp1, and also at late stages of meiotic differentiation when nuclear-exosome transcripts decreased. The pervasive exon-skipping transcripts were stochastic, did not increase in specific physiological conditions, and were mostly present at less than one copy per cell, even in the absence of nuclear RNA surveillance and during late meiosis. These exon-skipping transcripts are therefore unlikely to be functional and may reflect splicing errors that are actively removed by nuclear RNA surveillance. The average splicing rate by exon skipping was ∼0.24% in wild type and ∼1.75% in nuclear exonuclease mutants. We also detected approximately 250 circular RNAs derived from single or multiple exons. These circular RNAs were rare and stochastic, although a few became stabilized during quiescence and in splicing mutants. Using an exhaustive search algorithm, we also uncovered thousands of previously unknown splice sites, indicating pervasive splicing; yet most of these splicing variants were cryptic and increased in nuclear degradation mutants. This study highlights widespread but low frequency alternative or aberrant splicing events that are targeted by nuclear RNA surveillance.

DNA methylation and expression of KCNQ3 in bipolar disorder

OBJECTIVES

Accumulating evidence implicates the potassium voltage-gated channel, KQT-like subfamily, member 2 and 3 (KCNQ2 and KCNQ3) genes in the etiology of bipolar disorder (BPD). Reduced KCNQ2 or KCNQ3 gene expression might lead to a loss of inhibitory M-current and an increase in neuronal hyperexcitability in disease. The goal of the present study was to evaluate epigenetic and gene expression associations of the KCNQ2 and KCNQ3 genes with BPD.

METHODS

DNA methylation and gene expression levels of alternative transcripts of KCNQ2 and KCNQ3 capable of binding the ankyrin G (ANK3) gene were evaluated using bisulfite pyrosequencing and the quantitative real-time polymerase chain reaction in the postmortem prefrontal cortex of subjects with BPD and matched controls from the McLean Hospital. Replication analyses of DNA methylation findings were performed using prefrontal cortical DNA obtained from the Stanley Medical Research Institute.

RESULTS

Significantly lower expression was observed in KCNQ3, but not KCNQ2. DNA methylation analysis of CpGs within an alternative exonic region of KCNQ3 exon 11 demonstrated significantly lower methylation in BPD, and correlated significantly with KCNQ3 mRNA levels. Lower KCNQ3 exon 11 DNA methylation was observed in the Stanley Medical Research Institute replication cohort, although only after correcting for mood stabilizer status. Mood stabilizer treatment in rats resulted in a slight DNA methylation increase at the syntenic KCNQ3 exon 11 region, which subsequent analyses suggested could be the result of alterations in neuronal proportion.

CONCLUSION

The results of the present study suggest that epigenetic alterations in the KCNQ3 gene may be important in the etiopathogenesis of BPD and highlight the importance of controlling for medication and cellular composition-induced heterogeneity in psychiatric studies of the brain.

Adaptive thermal control of stem gravitropism through alternative RNA splicing in Arabidopsis

Gravitropism is an important growth movement in response to gravity in virtually all higher plants: the roots showing positive gravitropism and the shoots showing negative gravitropism. The gravitropic orientation of plant organs is also influenced by environmental factors, such as light and temperature. It is known that a zinc finger (ZF)-containing transcription factor SHOOT GRAVITROPISM 5/INDETERMINATE DOMAIN 15 (SGR5/IDD15) mediates the early events of gravitropic responses occurring in inflorescence stems. We have recently found that SGR5 gene undergoes alternative splicing to produce 2 protein variants, the full-size SGR5α transcription factor and the truncated SGR5β form lacking functional ZF motifs. The SGR5β form inhibits SGR5α function possibly by forming nonfunctional heterodimers that are excluded from DNA binding. Notably, SGR5 alternative splicing is accelerated at high temperatures, resulting in a high-level accumulation of SGR5β proteins. Accordingly, transgenic plants overexpressing SGR5β exhibit a reduction in the negative gravitropism of inflorescence stems, as observed in the SGR5-defective mutant. It is proposed that the thermos-responsive alternative splicing of SGR5 gene provides an adaptation strategy by which plants protect the shoots from aerial heat frequently occurring in natural habitats.

Genetic variation and RNA binding proteins: tools and techniques to detect functional polymorphisms

At its most fundamental level the goal of genetics is to connect genotype to phenotype. This question is asked at a basic level evaluating the role of genes and pathways in genetic model organism. Increasingly, this question is being asked in the clinic. Genomes of individuals and populations are being sequenced and compared. The challenge often comes at the stage of analysis. The variant positions are analyzed with the hope of understanding human disease. However after a genome or exome has been sequenced, the researcher is often deluged with hundreds of potentially relevant variations. Traditionally, amino-acid changing mutations were considered the tractable class of disease-causing mutations; however, mutations that disrupt noncoding elements are the subject of growing interest. These noncoding changes are a major avenue of disease (e.g., one in three hereditary disease alleles are predicted to affect splicing). Here, we review some current practices of medical genetics, the basic theory behind biochemical binding and functional assays, and then explore technical advances in how variations that alter RNA protein recognition events are detected and studied. These advances are advances in scale-high-throughput implementations of traditional biochemical assays that are feasible to perform in any molecular biology laboratory. This chapter utilizes a case study approach to illustrate some methods for analyzing polymorphisms. The first characterizes a functional intronic SNP that deletes a high affinity PTB site using traditional low-throughput biochemical and functional assays. From here we demonstrate the utility of high-throughput splicing and spliceosome assembly assays for screening large sets of SNPs and disease alleles for allelic differences in gene expression. Finally we perform three pilot drug screens with small molecules (G418, tetracycline, and valproic acid) that illustrate how compounds that rescue specific instances of differential pre-mRNA processing can be discovered.

Reverse transcriptase and intron number evolution

BACKGROUND

Introns are universal in eukaryotic genomes and play important roles in transcriptional regulation, mRNA export to the cytoplasm, nonsense-mediated decay as both a regulatory and a splicing quality control mechanism, R-loop avoidance, alternative splicing, chromatin structure, and evolution by exon-shuffling.

METHODS

Sixteen complete fungal genomes were used 13 of which were sequenced and annotated by JGI. Ustilago maydis, Cryptococcus neoformans, and Coprinus cinereus (also named Coprinopsis cinerea) were from the Broad Institute. Gene models from JGI-annotated genomes were taken from the GeneCatalog track that contained the best representative gene models. Varying fractions of the GeneCatalog were manually curated by external users. For clarity, we used the JGI unique database identifier.

RESULTS

The last common ancestor of eukaryotes (LECA) has an estimated 6.4 coding exons per gene (EPG) and evolved into the diverse eukaryotic life forms, which is recapitulated by the development of a stem cell. We found a parallel between the simulated reverse transcriptase (RT)-mediated intron loss and the comparative analysis of 16 fungal genomes that spanned a wide range of intron density. Although footprints of RT (RTF) were dynamic, relative intron location (RIL) to the 5'-end of mRNA faithfully traced RT-mediated intron loss and revealed 7.7 EPG for LECA. The mode of exon length distribution was conserved in simulated intron loss, which was exemplified by the shared mode of 75 nt between fungal and Chlamydomonas genomes. The dominant ancient exon length was corroborated by the average exon length of the most intron-rich genes in fungal genomes and consistent with ancient protein modules being ~25 aa. Combined with the conservation of a protein length of 400 aa, the earliest ancestor of eukaryotes could have 16 EPG. During earlier evolution, Ascomycota's ancestor had significantly more 3'-biased RT-mediated intron loss that was followed by dramatic RTF loss. There was a down trend of EPG from more conserved to less conserved genes. Moreover, species-specific genes have higher exon-densities, shorter exons, and longer introns when compared to genes conserved at the phylum level. However, intron length in species-specific genes became shorter than that of genes conserved in all species after genomes experiencing drastic intron loss. The estimated EPG from the most frequent exon length is more than double that from the RIL method.

CONCLUSIONS

This implies significant intron loss during the very early period of eukaryotic evolution. De novo gene-birth contributes to shorter exons, longer introns, and higher exon-density in species-specific genes relative to conserved genes.

Insulators recruit histone methyltransferase dMes4 to regulate chromatin of flanking genes

Chromosomal domains in Drosophila are marked by the insulator-binding proteins (IBPs) dCTCF/Beaf32 and cofactors that participate in regulating long-range interactions. Chromosomal borders are further enriched in specific histone modifications, yet the role of histone modifiers and nucleosome dynamics in this context remains largely unknown. Here, we show that IBP depletion impairs nucleosome dynamics specifically at the promoters and coding sequence of genes flanked by IBP binding sites. Biochemical purification identifies the H3K36 histone methyltransferase NSD/dMes-4 as a novel IBP cofactor, which specifically co-regulates the chromatin accessibility of hundreds of genes flanked by dCTCF/Beaf32. NSD/dMes-4 presets chromatin before the recruitment of transcriptional activators including DREF that triggers Set2/Hypb-dependent H3K36 trimethylation, nucleosome positioning, and RNA splicing. Our results unveil a model for how IBPs regulate nucleosome dynamics and gene expression through NSD/dMes-4, which may regulate H3K27me3 spreading. Our data uncover how IBPs dynamically regulate chromatin organization depending on distinct cofactors.

Genome-Wide Analysis of Heat-Sensitive Alternative Splicing in Physcomitrella patens

Plant growth and development are constantly influenced by temperature fluctuations. To respond to temperature changes, different levels of gene regulation are modulated in the cell. Alternative splicing (AS) is a widespread mechanism increasing transcriptome complexity and proteome diversity. Although genome-wide studies have revealed complex AS patterns in plants, whether AS impacts the stress defense of plants is not known. We used heat shock (HS) treatments at nondamaging temperature and messenger RNA sequencing to obtain HS transcriptomes in the moss Physcomitrella patens. Data analysis identified a significant number of novel AS events in the moss protonema. Nearly 50% of genes are alternatively spliced. Intron retention (IR) is markedly repressed under elevated temperature but alternative donor/acceptor site and exon skipping are mainly induced, indicating differential regulation of AS in response to heat stress. Transcripts undergoing heat-sensitive IR are mostly involved in specific functions, which suggests that plants regulate AS with transcript specificity under elevated temperature. An exonic GAG-repeat motif in these IR regions may function as a regulatory cis-element in heat-mediated AS regulation. A conserved AS pattern for HS transcription factors in P. patens and Arabidopsis (Arabidopsis thaliana) reveals that heat regulation for AS evolved early during land colonization of green plants. Our results support that AS of specific genes, including key HS regulators, is fine-tuned under elevated temperature to modulate gene regulation and reorganize metabolic processes.