CAGI, the Critical Assessment of Genome Interpretation, establishes progress and prospects for computational genetic variant interpretation methods

BACKGROUND

The Critical Assessment of Genome Interpretation (CAGI) aims to advance the state-of-the-art for computational prediction of genetic variant impact, particularly where relevant to disease. The five complete editions of the CAGI community experiment comprised 50 challenges, in which participants made blind predictions of phenotypes from genetic data, and these were evaluated by independent assessors.

RESULTS

Performance was particularly strong for clinical pathogenic variants, including some difficult-to-diagnose cases, and extends to interpretation of cancer-related variants. Missense variant interpretation methods were able to estimate biochemical effects with increasing accuracy. Assessment of methods for regulatory variants and complex trait disease risk was less definitive and indicates performance potentially suitable for auxiliary use in the clinic.

CONCLUSIONS

Results show that while current methods are imperfect, they have major utility for research and clinical applications. Emerging methods and increasingly large, robust datasets for training and assessment promise further progress ahead.

Targeting IGF2BP3 enhances antileukemic effects of menin-MLL inhibition in MLL-AF4 leukemia

ABSTRACT

RNA-binding proteins (RBPs) are emerging as a novel class of therapeutic targets in cancer, including in leukemia, given their important role in posttranscriptional gene regulation, and have the unexplored potential to be combined with existing therapies. The RBP insulin-like growth factor 2 messenger RNA-binding protein 3 (IGF2BP3) has been found to be a critical regulator of MLL-AF4 leukemogenesis and represents a promising therapeutic target. Here, we study the combined effects of targeting IGF2BP3 and menin-MLL interaction in MLL-AF4-driven leukemia in vitro and in vivo, using genetic inhibition with CRISPR-Cas9-mediated deletion of Igf2bp3 and pharmacologic inhibition of the menin-MLL interaction with multiple commercially available inhibitors. Depletion of Igf2bp3 sensitized MLL-AF4 leukemia to the effects of menin-MLL inhibition on cell growth and leukemic initiating cells in vitro. Mechanistically, we found that both Igf2bp3 depletion and menin-MLL inhibition led to increased differentiation in vitro and in vivo, seen in functional readouts and by gene expression analyses. IGF2BP3 knockdown had a greater effect on increasing survival and attenuating disease than pharmacologic menin-MLL inhibition with small molecule MI-503 alone and showed enhanced antileukemic effects in combination. Our work shows that IGF2BP3 is an oncogenic amplifier of MLL-AF4-mediated leukemogenesis and a potent therapeutic target, providing a paradigm for targeting leukemia at both the transcriptional and posttranscriptional level.

An intronic RNA element modulates Factor VIII exon-16 splicing

Pathogenic variants in the human Factor VIII (F8) gene cause Hemophilia A (HA). Here, we investigated the impact of 97 HA-causing single-nucleotide variants on the splicing of 11 exons from F8. For the majority of F8 exons, splicing was insensitive to the presence of HA-causing variants. However, splicing of several exons, including exon-16, was impacted by variants predicted to alter exonic splicing regulatory sequences. Using exon-16 as a model, we investigated the structure-function relationship of HA-causing variants on splicing. Intriguingly, RNA chemical probing analyses revealed a three-way junction structure at the 3'-end of intron-15 (TWJ-3-15) capable of sequestering the polypyrimidine tract. We discovered antisense oligonucleotides (ASOs) targeting TWJ-3-15 partially rescue splicing-deficient exon-16 variants by increasing accessibility of the polypyrimidine tract. The apical stem loop region of TWJ-3-15 also contains two hnRNPA1-dependent intronic splicing silencers (ISSs). ASOs blocking these ISSs also partially rescued splicing. When used in combination, ASOs targeting both the ISSs and the region sequestering the polypyrimidine tract, fully rescue pre-mRNA splicing of multiple HA-linked variants of exon-16. Together, our data reveal a putative RNA structure that sensitizes F8 exon-16 to aberrant splicing.

Distinct oncogenic phenotypes in hematopoietic specific deletions of Trp53

Loss of function in the tumor suppressor gene TP53 is the most common alteration seen in human cancer. In mice, P53 deletion in all cells leads predominantly to the development of T-cell lymphomas, followed by B-cell lymphomas, sarcomas and teratomas. In order to dissect the role of P53 in the hematopoietic system, we generated and analyzed two different mouse models deficient for P53. A pan-hematopoietic P53 deletion mouse was created using Vav1-Cre based deletion; and a B-cell-specific deletion mouse was created using a CD19-Cre based deletion. The Vav1-P53CKO mice predominantly developed T-cell malignancies in younger mice, and myeloid malignancies in older mice. In T-cell malignancies, there was accelerated thymic cell maturation with overexpression of Notch1 and its downstream effectors. CD19-P53CKO mice developed marginal zone expansion in the spleen, followed by marginal zone lymphoma, some of which progressed to diffuse large B-cell lymphomas. Interestingly, marginal zone and diffuse large B-cell lymphomas had a unique gene expression signature characterized by activation of the PI3K pathway, compared with wild type marginal zone or follicular cells of the spleen. This study demonstrates lineage specific P53 deletion leading to distinct phenotypes secondary to unique gene expression programs set in motion.

Isoform-specific translational control is evolutionarily conserved in primates

Alternative splicing (AS) alters messenger RNA (mRNA) coding capacity, localization, stability, and translation. Here we use comparative transcriptomics to identify cis- acting elements coupling AS to translational control (AS-TC). We sequenced total cytosolic and polyribosome-associated mRNA from human, chimpanzee, and orangutan induced pluripotent stem cells (iPSCs), revealing thousands of transcripts with splicing differences between subcellular fractions. We found both conserved and species-specific polyribosome association patterns for orthologous splicing events. Intriguingly, alternative exons with similar polyribosome profiles between species have stronger sequence conservation than exons with lineage-specific ribosome association. These data suggest that sequence variation underlies differences in the polyribosome association. Accordingly, single nucleotide substitutions in luciferase reporters designed to model exons with divergent polyribosome profiles are sufficient to regulate translational efficiency. We used position specific weight matrices to interpret exons with species-specific polyribosome association profiles, finding that polymorphic sites frequently alter recognition motifs for trans- acting RNA binding proteins. Together, our results show that AS can regulate translation by remodeling the cis- regulatory landscape of mRNA isoforms.

Rescue of blood coagulation Factor VIII exon-16 mis-splicing by antisense oligonucleotides

The human Factor VIII ( F8 ) protein is essential for the blood coagulation cascade and specific F8 mutations cause the rare bleeding disorder Hemophilia A (HA). Here, we investigated the impact of HA-causing single-nucleotide mutations on F8 pre-mRNA splicing. We found that 14/97 (∼14.4%) coding sequence mutations tested in our study induced exon skipping. Splicing patterns of 4/11 (∼36.4%) F8 exons tested were especially sensitive to the presence of common disease-causing mutations. RNA-chemical probing analyses revealed a three-way junction structure at the 3' end of intron 15 (TWJ-3-15). TWJ-3-15 sequesters the polypyrimidine tract, a key determinant of 3' splice site strength. Using exon-16 of the F8 gene as a model, we designed specific antisense oligonucleotides (ASOs) that target TWJ-3-15 and identified three that promote the splicing of F8 exon-16. Interaction of TWJ-3-15 with ASOs increases accessibility of the polypyrimidine tract and inhibits the binding of hnRNPA1-dependent splicing silencing factors. Moreover, ASOs targeting TWJ-3-15 rescue diverse splicing-sensitive HA-causing mutations, most of which are distal to the 3' splice site being impacted. The TWJ-3-15 structure and its effect on mRNA splicing provide a model for HA etiology in patients harboring specific F8 mutations and provide a framework for precision RNA-based HA therapies.

A potential role for SARS-CoV-2 small viral RNAs in targeting host microRNAs and modulating gene expression

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus disease (COVID-19) in humans, with symptoms ranging from mild to severe, including fatality. The molecular mechanisms surrounding the effects of viral infection on the host RNA machinery remain poorly characterized. We used a comparative transcriptomics approach to investigate the effects of SARS-CoV-2 infection on the host mRNA and sRNA expression machinery in a human lung epithelial cell line (Calu-3) and an African green monkey kidney cell line (Vero-E6). Upon infection, we observed global changes in host gene expression and differential expression of dozens of host miRNAs, many with known links to viral infection and immune response. Additionally, we discovered an expanded landscape of more than a hundred SARS-CoV-2-derived small viral RNAs (svRNAs) predicted to interact with differentially expressed host mRNAs and miRNAs. svRNAs are derived from distinct regions of the viral genome and sequence signatures suggest they are produced by a non-canonical biogenesis pathway. 52 of the 67 svRNAs identified in Calu-3 cells are predicted to interact with differentially expressed miRNAs, with many svRNAs having multiple targets. Accordingly, we speculate that these svRNAs may play a role in SARS-CoV-2 propagation by modulating post-transcriptional gene regulation, and that methods for antagonizing them may have therapeutic value.

The RNA-binding protein IGF2BP3 is critical for MLL-AF4-mediated leukemogenesis

Despite recent advances in therapeutic approaches, patients with MLL-rearranged leukemia still have poor outcomes. Here, we find that the RNA-binding protein IGF2BP3, which is overexpressed in MLL-translocated leukemia, strongly amplifies MLL-Af4-mediated leukemogenesis. Deletion of Igf2bp3 significantly increases the survival of mice with MLL-Af4-driven leukemia and greatly attenuates disease, with a minimal impact on baseline hematopoiesis. At the cellular level, MLL-Af4 leukemia-initiating cells require Igf2bp3 for their function in leukemogenesis. At the molecular level, IGF2BP3 regulates a complex posttranscriptional operon governing leukemia cell survival and proliferation. IGF2BP3-targeted mRNA transcripts include important MLL-Af4-induced genes, such as those in the Hoxa locus, and the Ras signaling pathway. Targeting of transcripts by IGF2BP3 regulates both steady-state mRNA levels and, unexpectedly, pre-mRNA splicing. Together, our findings show that IGF2BP3 represents an attractive therapeutic target in this disease, providing important insights into mechanisms of posttranscriptional regulation in leukemia.

A method for campus-wide SARS-CoV-2 surveillance at a large public university

The systematic screening of asymptomatic and pre-symptomatic individuals is a powerful tool for controlling community transmission of infectious disease on college campuses. Faced with a paucity of testing in the beginning of the COVID-19 pandemic, many universities developed molecular diagnostic laboratories focused on SARS-CoV-2 diagnostic testing on campus and in their broader communities. We established the UC Santa Cruz Molecular Diagnostic Lab in early April 2020 and began testing clinical samples just five weeks later. Using a clinically-validated laboratory developed test (LDT) that avoided supply chain constraints, an automated sample pooling and processing workflow, and a custom laboratory information management system (LIMS), we expanded testing from a handful of clinical samples per day to thousands per day with the testing capacity to screen our entire campus population twice per week. In this report we describe the technical, logistical, and regulatory processes that enabled our pop-up lab to scale testing and reporting capacity to thousands of tests per day.

Treatment of horizontal silage bunker runoff using biochar amended vegetative filter strips

Horizontal silage bunkers produce leachate that contains contaminants that can be detrimental to the environment if released untreated. Vegetated filter strips are used to treat silage bunker runoff to prevent contamination of surface waters via infiltration, however increased infiltration poses risks to groundwater, particularly for nitrate (NO₃^-). Vegetated filter strip plots with a sandy loam soil, half of which are amended with biochar, were investigated to assess the treatment of silage bunker runoff over 20 application events. The subsurface effluent biological oxygen demand (BOD₅), chemical oxygen demand (COD), and total phosphorus (TP) were reduced on average by 40%, 46%, and 75%, respectively, and there was no statistical difference between treatments. The total nitrogen (TN) was reduced by 49 and 64% for control and biochar plots, respectively, which was significantly different between treatments. Biochar significantly reduced nitrate nitrogen (NO₃^--N) leaching by 40% compared to the control, however, the NO₃^--N concentration in leachate was still high ranging from 0.19 to 191.04 mg NO₃^--N L^-1 and 0.18-108.89 mg NO₃^--N L^-1 for control and biochar plots, respectively. A mass balance suggests the primary mechanism for a decrease in TN and NO₃^--N leaching from biochar amended plots was greater retention of NO₃^--N and organic N (ORG-N) within the soil/biochar matrix. The development of oxygenated functional groups and/or formation of organomineral layer on the biochar surface likely enhanced N retention.

Nitrate sorption to biochar following chemical oxidation

Biochar amendments can reduce nitrate (NO₃) leaching in agricultural soil. It has been hypothesized that functional groups on the biochar surface from oxidation can increase NO₃ sorption. This study evaluates the effect of chemical oxidation of biochar on NO₃ sorption characteristics. Eight biochars, made from wood and corn cobs, underwent sodium hypochlorite (NaClO) and hydrogen peroxide (H₂O₂) oxidation and then assessed for NO₃ sorption capacity using batch isotherm methods. The unoxidized and oxidized biochar produced at low temperatures (400 °C) had no significant NO₃ sorption. Oxidized biochars produced at higher temperatures (600 °C and 700 °C) had calculated maximum NO₃ sorption capacities (S_max) ranging from 0.50 to 3.97 mg NO₃-N g^-1. Biochar oxidations with 50 mmol NaClO g^-1 (N50) in combination with an acid wash (AW) had the largest estimated sorption capacities of 3.68, 3.97, and 1.46 mg NO₃-N g^-1 for CT_N50,AW, BW3_N50,AW, and CC3_N50,AW, respectively. Sorption capacity of wood-based biochars was higher than corn cob biochars due to increased oxidation as measured by total acid group content (TAGC). Wood biochar Smax values were correlated with ΔTAGC (R²= 0.86), with a slope of 1.2 μmol NO₃-N μmol TAGC^-1 suggesting that cationic bridging of NO₃ to oxidized sites is the primary mechanism for NO₃ sorption.

HNRNPA1 promotes recognition of splice site decoys by U2AF2 in vivo

Alternative pre-mRNA splicing plays a major role in expanding the transcript output of human genes. This process is regulated, in part, by the interplay of trans-acting RNA binding proteins (RBPs) with myriad cis-regulatory elements scattered throughout pre-mRNAs. These molecular recognition events are critical for defining the protein-coding sequences (exons) within pre-mRNAs and directing spliceosome assembly on noncoding regions (introns). One of the earliest events in this process is recognition of the 3′ splice site (3′ss) by U2 small nuclear RNA auxiliary factor 2 (U2AF2). Splicing regulators, such as the heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1), influence spliceosome assembly both in vitro and in vivo, but their mechanisms of action remain poorly described on a global scale. HNRNPA1 also promotes proofreading of 3′ss sequences though a direct interaction with the U2AF heterodimer. To determine how HNRNPA1 regulates U2AF–RNA interactions in vivo, we analyzed U2AF2 RNA binding specificity using individual-nucleotide resolution crosslinking immunoprecipitation (iCLIP) in control and HNRNPA1 overexpression cells. We observed changes in the distribution of U2AF2 crosslinking sites relative to the 3′ss of alternative cassette exons but not constitutive exons upon HNRNPA1 overexpression. A subset of these events shows a concomitant increase of U2AF2 crosslinking at distal intronic regions, suggesting a shift of U2AF2 to “decoy” binding sites. Of the many noncanonical U2AF2 binding sites, Alu-derived RNA sequences represented one of the most abundant classes of HNRNPA1-dependent decoys. We propose that one way HNRNPA1 regulates exon definition is to modulate the interaction of U2AF2 with decoy or bona fide 3′ss.

The 3' end of the story: deciphering combinatorial interactions that control mRNA fate

A new study investigates how microRNAs affect the binding of proteins to RNA.

Modulation of nonsense mediated decay by rapamycin

Rapamycin is a naturally occurring macrolide whose target is at the core of nutrient and stress regulation in a wide range of species. Despite well-established roles as an inhibitor of cap-dependent mRNA translation, relatively little is known about its effects on other modes of RNA processing. Here, we characterize the landscape of rapamycin-induced post-transcriptional gene regulation. Transcriptome analysis of rapamycin-treated cells reveals genome-wide changes in alternative mRNA splicing and pronounced changes in NMD-sensitive isoforms. We demonstrate that despite well-documented attenuation of cap-dependent mRNA translation, rapamycin can augment NMD of certain transcripts. Rapamycin-treatment significantly reduces the levels of both endogenous and exogenous Premature Termination Codon (PTC)-containing mRNA isoforms and its effects are dose-, UPF1- and 4EBP-dependent. The PTC-containing SRSF6 transcript exhibits a shorter half-life upon rapamycin-treatment as compared to the non-PTC isoform. Rapamycin-treatment also causes depletion of PTC-containing mRNA isoforms from polyribosomes, underscoring the functional relationship between translation and NMD. Enhanced NMD activity also correlates with an enrichment of the nuclear Cap Binding Complex (CBC) in rapamycin-treated cells. Our data demonstrate that rapamycin modulates global RNA homeostasis by NMD.

Improved genome assembly of American alligator genome reveals conserved architecture of estrogen signaling

The American alligator, Alligator mississippiensis, like all crocodilians, has temperature-dependent sex determination, in which the sex of an embryo is determined by the incubation temperature of the egg during a critical period of development. The lack of genetic differences between male and female alligators leaves open the question of how the genes responsible for sex determination and differentiation are regulated. Insight into this question comes from the fact that exposing an embryo incubated at male-producing temperature to estrogen causes it to develop ovaries. Because estrogen response elements are known to regulate genes over long distances, a contiguous genome assembly is crucial for predicting and understanding their impact. We present an improved assembly of the American alligator genome, scaffolded with in vitro proximity ligation (Chicago) data. We use this assembly to scaffold two other crocodilian genomes based on synteny. We perform RNA sequencing of tissues from American alligator embryos to find genes that are differentially expressed between embryos incubated at male- versus female-producing temperature. Finally, we use the improved contiguity of our assembly along with the current model of CTCF-mediated chromatin looping to predict regions of the genome likely to contain estrogen-responsive genes. We find that these regions are significantly enriched for genes with female-biased expression in developing gonads after the critical period during which sex is determined by incubation temperature. We thus conclude that estrogen signaling is a major driver of female-biased gene expression in the post-temperature sensitive period gonads.

RNA-binding protein IGF2BP3 targeting of oncogenic transcripts promotes hematopoietic progenitor proliferation

Posttranscriptional control of gene expression is important for defining both normal and pathological cellular phenotypes. In vitro, RNA-binding proteins (RBPs) have recently been shown to play important roles in posttranscriptional regulation; however, the contribution of RBPs to cell specification is not well understood. Here, we determined that the RBP insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3) is specifically overexpressed in mixed lineage leukemia-rearranged (MLL-rearranged) B-acute lymphoblastic leukemia (B-ALL), which constitutes a subtype of this malignancy associated with poor prognosis and high risk of relapse. IGF2BP3 was required for the survival of B-ALL cell lines, as knockdown led to decreased proliferation and increased apoptosis. Enforced expression of IGF2BP3 provided murine BM cells with a strong survival advantage, led to proliferation of hematopoietic stem and progenitor cells, and skewed hematopoietic development to the B cell/myeloid lineage. Cross-link immunoprecipitation and high throughput sequencing uncovered the IGF2BP3-regulated transcriptome, which includes oncogenes MYC and CDK6 as direct targets. IGF2BP3 regulated transcripts via targeting elements within 3' untranslated regions (3'UTR), and enforced IGF2BP3 expression in mice resulted in enhanced expression of Myc and Cdk6 in BM. Together, our data suggest that IGF2BP3-mediated targeting of oncogenic transcripts may represent a critical pathogenetic mechanism in MLL-rearranged B-ALL and support IGF2BP3 and its cognate RNA-binding partners as potential therapeutic targets in this disease.

The RNAissance family: SR proteins as multifaceted regulators of gene expression

Serine and arginine-rich (SR) proteins play multiple roles in the eukaryotic gene expression pathway. Initially described as constitutive and alternative splicing factors, now it is clear that SR proteins are key determinants of exon identity and function as molecular adaptors, linking the pre-messenger RNA (pre-mRNA) to the splicing machinery. In addition, now SR proteins are implicated in many aspects of mRNA and noncoding RNA (ncRNA) processing well beyond splicing. These unexpected roles, including RNA transcription, export, translation, and decay, may prove to be the rule rather than the exception. To simply define, this family of RNA-binding proteins as splicing factors belies the broader roles of SR proteins in post-transcriptional gene expression.

Exon identity crisis: disease-causing mutations that disrupt the splicing code

Cis-acting RNA elements control the accurate expression of human multi-exon protein coding genes. Single nucleotide variants altering the fidelity of this regulatory code and, consequently, pre-mRNA splicing are expected to contribute to the etiology of numerous human diseases.

MutPred Splice: machine learning-based prediction of exonic variants that disrupt splicing

We have developed a novel machine-learning approach, MutPred Splice, for the identification of coding region substitutions that disrupt pre-mRNA splicing. Applying MutPred Splice to human disease-causing exonic mutations suggests that 16% of mutations causing inherited disease and 10 to 14% of somatic mutations in cancer may disrupt pre-mRNA splicing. For inherited disease, the main mechanism responsible for the splicing defect is splice site loss, whereas for cancer the predominant mechanism of splicing disruption is predicted to be exon skipping via loss of exonic splicing enhancers or gain of exonic splicing silencer elements. MutPred Splice is available at http://mutdb.org/mutpredsplice.

Frac-seq reveals isoform-specific recruitment to polyribosomes

Pre-mRNA splicing is required for the accurate expression of virtually all human protein coding genes. However, splicing also plays important roles in coordinating subsequent steps of pre-mRNA processing such as polyadenylation and mRNA export. Here, we test the hypothesis that nuclear pre-mRNA processing influences the polyribosome association of alternative mRNA isoforms. By comparing isoform ratios in cytoplasmic and polyribosomal extracts, we determined that the alternative products of ∼30% (597/1954) of mRNA processing events are differentially partitioned between these subcellular fractions. Many of the events exhibiting isoform-specific polyribosome association are highly conserved across mammalian genomes, underscoring their possible biological importance. We find that differences in polyribosome association may be explained, at least in part by the observation that alternative splicing alters the cis-regulatory landscape of mRNAs isoforms. For example, inclusion or exclusion of upstream open reading frames (uORFs) in the 5′UTR as well as Alu-elements and microRNA target sites in the 3′UTR have a strong influence on polyribosome association of alternative mRNA isoforms. Taken together, our data demonstrate for the first time the potential link between alternative splicing and translational control of the resultant mRNA isoforms.

Site identification in high-throughput RNA-protein interaction data

MOTIVATION

Post-transcriptional and co-transcriptional regulation is a crucial link between genotype and phenotype. The central players are the RNA-binding proteins, and experimental technologies [such as cross-linking with immunoprecipitation- (CLIP-) and RIP-seq] for probing their activities have advanced rapidly over the course of the past decade. Statistically robust, flexible computational methods for binding site identification from high-throughput immunoprecipitation assays are largely lacking however.

RESULTS

We introduce a method for site identification which provides four key advantages over previous methods: (i) it can be applied on all variations of CLIP and RIP-seq technologies, (ii) it accurately models the underlying read-count distributions, (iii) it allows external covariates, such as transcript abundance (which we demonstrate is highly correlated with read count) to inform the site identification process and (iv) it allows for direct comparison of site usage across cell types or conditions.

AVAILABILITY AND IMPLEMENTATION

We have implemented our method in a software tool called Piranha. Source code and binaries, licensed under the GNU General Public License (version 3) are freely available for download from http://smithlab.usc.edu.

CONTACT

andrewds@usc.edu

SUPPLEMENTARY INFORMATION

Supplementary data available at Bioinformatics online.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

Predicting sequence and structural specificities of RNA binding regions recognized by splicing factor SRSF1

BACKGROUND

RNA-binding proteins (RBPs) play diverse roles in eukaryotic RNA processing. Despite their pervasive functions in coding and noncoding RNA biogenesis and regulation, elucidating the sequence specificities that define protein-RNA interactions remains a major challenge. Recently, CLIP-seq (Cross-linking immunoprecipitation followed by high-throughput sequencing) has been successfully implemented to study the transcriptome-wide binding patterns of SRSF1, PTBP1, NOVA and fox2 proteins. These studies either adopted traditional methods like Multiple EM for Motif Elicitation (MEME) to discover the sequence consensus of RBP's binding sites or used Z-score statistics to search for the overrepresented nucleotides of a certain size. We argue that most of these methods are not well-suited for RNA motif identification, as they are unable to incorporate the RNA structural context of protein-RNA interactions, which may affect to binding specificity. Here, we describe a novel model-based approach--RNAMotifModeler to identify the consensus of protein-RNA binding regions by integrating sequence features and RNA secondary structures.

RESULTS

As an example, we implemented RNAMotifModeler on SRSF1 (SF2/ASF) CLIP-seq data. The sequence-structural consensus we identified is a purine-rich octamer 'AGAAGAAG' in a highly single-stranded RNA context. The unpaired probabilities, the probabilities of not forming pairs, are significantly higher than negative controls and the flanking sequence surrounding the binding site, indicating that SRSF1 proteins tend to bind on single-stranded RNA. Further statistical evaluations revealed that the second and fifth bases of SRSF1octamer motif have much stronger sequence specificities, but weaker single-strandedness, while the third, fourth, sixth and seventh bases are far more likely to be single-stranded, but have more degenerate sequence specificities. Therefore, we hypothesize that nucleotide specificity and secondary structure play complementary roles during binding site recognition by SRSF1.

CONCLUSION

In this study, we presented a computational model to predict the sequence consensus and optimal RNA secondary structure for protein-RNA binding regions. The successful implementation on SRSF1 CLIP-seq data demonstrates great potential to improve our understanding on the binding specificity of RNA binding proteins.

Prioritizing single-nucleotide variations that potentially regulate alternative splicing

Recent evidence suggests that many complex diseases are caused by genetic variations that play regulatory roles in controlling gene expression. Most genetic studies focus on nonsynonymous variations that can alter the amino acid composition of a protein and are therefore believed to have the highest impact on phenotype. Synonymous variations, however, can also play important roles in disease pathogenesis by regulating pre-mRNA processing and translational control. In this study, we systematically survey the effects of single-nucleotide variations (SNVs) on binding affinity of RNA-binding proteins (RBPs). Among the 10,113 synonymous SNVs identified in 697 individuals in the 1,000 Genomes Project and distributed by Genetic Analysis Workshop 17 (GAW17), we identified 182 variations located in alternatively spliced exons that can significantly change the binding affinity of nine RBPs whose binding preferences on 7-mer RNA sequences were previously reported. We found that the minor allele frequencies of these variations are similar to those of nonsynonymous SNVs, suggesting that they are in fact functional. We propose a workflow to identify phenotype-associated regulatory SNVs that might affect alternative splicing from exome-sequencing-derived genetic variations. Based on the affecting SNVs on the quantitative traits simulated in GAW17, we further identified two and four functional SNVs that are predicted to be involved in alternative splicing regulation in traits Q1 and Q2, respectively.

Loss of exon identity is a common mechanism of human inherited disease

It is widely accepted that at least 10% of all mutations causing human inherited disease disrupt splice-site consensus sequences. In contrast to splice-site mutations, the role of auxiliary cis-acting elements such as exonic splicing enhancers (ESE) and exonic splicing silencers (ESS) in human inherited disease is still poorly understood. Here we use a top-down approach to determine rates of loss or gain of known human exonic splicing regulatory (ESR) sequences associated with either disease-causing mutations or putatively neutral single nucleotide polymorphisms (SNPs). We observe significant enrichment toward loss of ESEs and gain of ESSs among inherited disease-causing variants relative to neutral polymorphisms, indicating that exon skipping may play a prominent role in aberrant gene regulation. Both computational and biochemical approaches underscore the relevance of exonic splicing enhancer loss and silencer gain in inherited disease. Additionally, we provide direct evidence that both SRp20 (SRSF3) and possibly PTB (PTBP1) are involved in the function of a splicing silencer that is created de novo by a total of 83 different inherited disease mutations in 67 different disease genes. Taken together, we find that ~25% (7154/27,681) of known mis-sense and nonsense disease-causing mutations alter functional splicing signals within exons, suggesting a much more widespread role for aberrant mRNA processing in causing human inherited disease than has hitherto been appreciated.

Genome-wide prediction of cis-acting RNA elements regulating tissue-specific pre-mRNA alternative splicing

Background

Human genes undergo various patterns of pre-mRNA splicing across different tissues. Such variation is primarily regulated by trans-acting factors that bind on exonic and intronic cis-acting RNA elements (CAEs). Here we report a computational method to mechanistically identify cis-acting RNA elements that contribute to the tissue-specific alternative splicing pattern. This method is an extension of our previous model, SplicingModeler, which predicts the significant CAEs that contribute to the splicing differences between two tissues. In this study, we introduce tissue-specific functional levels estimation step, which allows evaluating regulatory functions of predicted CAEs that are involved in more than two tissues.

Results

Using a publicly available Affymetrix Genechip^® Human Exon Array dataset, our method identifies 652 cis-acting RNA elements (CAEs) across 11 human tissues. About one third of predicted CAEs can be mapped to the known RBP (RNA binding protein) binding sites or match with other predicted exonic splicing regulator databases. Interestingly, the vast majority of predicted CAEs are in intronic regulatory regions. A noticeable exception is that many exonic elements are found to regulate the alternative splicing between cerebellum and testes. Most identified elements are found to contribute to the alternative splicing between two tissues, while some are important in multiple tissues. This suggests that genome-wide alternative splicing patterns are regulated by a combination of tissue-specific cis-acting elements and "general elements" whose functional activities are important but differ across multiple tissues.

Conclusion

In this study, we present a model-based computational approach to identify potential cis-acting RNA elements by considering the exon splicing variation as the combinatorial effects of multiple cis-acting regulators. This methodology provides a novel evaluation on the functional levels of cis-acting RNA elements by estimating their tissue-specific functions on various tissues.

DNA Methyltransferase protein synthesis is reduced in CXXC finger protein 1-deficient embryonic stem cells

CXXC finger protein 1 (CFP1) binds to unmethylated CpG dinucleotides and is required for embryogenesis. CFP1 is also a component of the Setd1A and Setd1B histone H3K4 methyltransferase complexes. Murine embryonic stem (ES) cells lacking CFP1 fail to differentiate, and exhibit a 70% reduction in global genomic cytosine methylation and a 50% reduction in DNA methyltransferase (DNMT1) protein and activity. This study investigated the underlying mechanism for reduced DNMT1 expression in CFP1-deficient ES cells. DNMT1 transcript levels were significantly elevated in ES cells lacking CFP1, despite the observed reduction in DNMT1 protein levels. To address the posttranscriptional mechanisms by which CFP1 regulates DNMT1 protein activity, pulse/chase analyses were carried out, demonstrating a modest reduction in DNMT1 protein half-life in CFP1-deficient ES cells. Additionally, global protein synthesis was decreased in ES cells lacking CFP1, contributing to a reduction in the synthesis of DNMT1 protein. ES cells lacking CFP1 were found to contain elevated levels of phosphorylated eIF2alpha, and an accompanying reduction in translation initiation as revealed by a lower level of polyribosomes. These results reveal a novel role for CFP1 in the regulation of translation initiation, and indicate that loss of CFP1 function leads to decreased DNMT1 protein synthesis and half-life.

Splicing factor SFRS1 recognizes a functionally diverse landscape of RNA transcripts

Metazoan genes are encrypted with at least two superimposed codes: the genetic code to specify the primary structure of proteins and the splicing code to expand their proteomic output via alternative splicing. Here, we define the specificity of a central regulator of pre-mRNA splicing, the conserved, essential splicing factor SFRS1. Cross-linking immunoprecipitation and high-throughput sequencing (CLIP-seq) identified 23,632 binding sites for SFRS1 in the transcriptome of cultured human embryonic kidney cells. SFRS1 was found to engage many different classes of functionally distinct transcripts including mRNA, miRNA, snoRNAs, ncRNAs, and conserved intergenic transcripts of unknown function. The majority of these diverse transcripts share a purine-rich consensus motif corresponding to the canonical SFRS1 binding site. The consensus site was not only enriched in exons cross-linked to SFRS1 in vivo, but was also enriched in close proximity to splice sites. mRNAs encoding RNA processing factors were significantly overrepresented, suggesting that SFRS1 may broadly influence the post-transcriptional control of gene expression in vivo. Finally, a search for the SFRS1 consensus motif within the Human Gene Mutation Database identified 181 mutations in 82 different genes that disrupt predicted SFRS1 binding sites. This comprehensive analysis substantially expands the known roles of human SR proteins in the regulation of a diverse array of RNA transcripts.

Using RNase sequence specificity to refine the identification of RNA-protein binding regions

Massively parallel pyrosequencing is a high-throughput technology that can sequence hundreds of thousands of DNA/RNA fragments in a single experiment. Combining it with immunoprecipitation-based biochemical assays, such as cross-linking immunoprecipitation (CLIP), provides a genome-wide method to detect the sites at which proteins bind DNA or RNA. In a CLIP-pyrosequencing experiment, the resolutions of the detected protein binding regions are partially determined by the length of the detected RNA fragments (CLIP amplicons) after trimming by RNase digestion. The lengths of these fragments usually range from 50-70 nucleotides. Many genomic regions are marked by multiple RNA fragments. In this paper, we report an empirical approach to refine the localization of protein binding regions by using the distribution pattern of the detected RNA fragments and the sequence specificity of RNase digestion. We present two regions to which multiple amplicons map as examples to demonstrate this approach.

The splicing factor SF2/ASF regulates translation initiation by enhancing phosphorylation of 4E-BP1

The SR protein SF2/ASF has been initially characterized as a splicing factor but has also been shown to mediate postsplicing activities such as mRNA export and translation. Here we demonstrate that SF2/ASF promotes translation initiation of bound mRNAs and that this activity requires the presence of the cytoplasmic cap-binding protein eIF4E. SF2/ASF promotes translation initiation by suppressing the activity of 4E-BP, a competitive inhibitor of cap-dependent translation. This activity is mediated by interactions of SF2/ASF with both mTOR and the phosphatase PP2A, two key regulators of 4E-BP phosphorylation. These findings suggest the model whereby SF2/ASF functions as an adaptor protein to recruit the signaling molecules responsible for regulation of cap-dependent translation of specific mRNAs. Taken together, these data suggest a novel mechanism for the activation of translation initiation of a subset of mRNAs bound by the shuttling protein SF2/ASF.

Translational control of C-terminal Src kinase (Csk) expression by PRL3 phosphatase

Phosphatase of regenerating liver 3 (PRL3) is up-regulated in cancer metastases. However, little is known of PRL3-mediated cellular signaling pathways. We previously reported that elevated PRL3 expression increases Src kinase activity, which likely contributes to the increased tumorigenesis and metastasis potential of PRL3. PRL3-induced Src activation is proposed to be indirect through down-regulation of Csk, a negative regulator of Src. Given the importance of PRL3 in tumor metastasis and the role of Csk in controlling Src activity, we addressed the mechanism by which PRL3 mediates Csk down-regulation. PRL3 is shown to exert a negative effect on Csk protein synthesis, rather than regulation of Csk mRNA levels or protein turnover. Interestingly, the preferential decrease in Csk protein synthesis is a consequence of increased eIF2 phosphorylation resulting from PRL3 expression. Reduced Csk synthesis also occurs in response to cellular stress that induces eIF2 phosphorylation, indicating that this regulatory mechanism may occur in response to a wider spectrum of cellular conditions known to direct translational control. Thus, we have uncovered a previously uncharacterized role for PRL3 in the gene-specific translational control of Csk expression.

A rapid and efficient protocol to purify biologically active recombinant proteins from mammalian cells

Here, we describe a simple and efficient method for the expression and purification of active recombinant proteins in mammalian cells. This method uses the expression of T7 epitope-tagged proteins in transiently transfected 293T cells grown in monolayer, followed by anti-T7-agarose affinity chromatography. This procedure yields approximately between 75 and 100 microg of biologically active protein/150 cm(2) flask that can be used for biochemical studies. We have tested this protocol for the expression of the prototype SR protein, SF2/ASF, which is a member of the SR protein family with a role in constitutive and alternative splicing. We show that SF2/ASF purified using this protocol is able to complement an S100 HeLa extract, demonstrating that is biologically active. Moreover, expression of a novel SR-related protein that it is required for the second step of pre-mRNA splicing also rendered an active protein. In summary, we present a protocol based on transient transfection of mammalian cells that results in easy purification of significant amounts of biologically active proteins.

Concerted regulation of nuclear and cytoplasmic activities of SR proteins by AKT

Serine/arginine-rich (SR) proteins are important regulators of mRNA splicing. Several postsplicing activities have been described for a subset of shuttling SR proteins, including regulation of mRNA export and translation. Using the fibronectin gene to study the links between signal-transduction pathways and SR protein activity, we show that growth factors not only modify the alternative splicing pattern of the fibronectin gene but also alter translation of reporter messenger RNAs in an SR protein-dependent fashion, providing two coregulated levels of isoform-specific amplification. These effects are inhibited by specific small interfering RNAs against SR proteins and are mediated by the AKT kinase, which elicits opposite effects to those evoked by overexpressing SR protein kinases Clk and SRPK. These results show how SR protein activity is modified in response to extracellular stimulation, leading to a concerted regulation of splicing and translation.

Reversible phosphorylation differentially affects nuclear and cytoplasmic functions of splicing factor 2/alternative splicing factor

The Ser/Arg-rich (SR) proteins constitute a family of highly conserved nuclear phosphoproteins that are involved in many steps of mRNA metabolism. Previously, we demonstrated that shuttling SR proteins can associate with translating ribosomes and enhance translation of reporter mRNAs both in vivo and in vitro. Here, we show that endogenous, cytoplasmic splicing factor 2/alternative splicing factor (SF2/ASF) associated with the translation machinery is hypophosphorylated, suggesting that the phosphorylation state of the Arg-Ser-rich (RS) domain may influence the role of SF2/ASF in cytoplasmic RNA processing. In agreement, we show that mutations mimicking a hypophosphorylated RS domain strongly increased SF2/ASF binding to cytoplasmic mRNA and its activity in translation. We also demonstrate that, whereas the RS domain is not required for the function of SF2/ASF in mRNA translation in vivo or in vitro, its second RNA recognition motif (RRM)2 plays a critical role in this process. Taken together, these data suggest that RS-domain phosphorylation may influence the association of SF2/ASF with mRNA, whereas RRM2 may play an important role in mediating protein-protein interactions during translation. These data are consistent with a model whereby reversible protein phosphorylation differentially regulates the subcellular localization and activity of shuttling SR proteins.

Multiple roles of arginine/serine-rich splicing factors in RNA processing

SR proteins (serine- and arginine-rich proteins) are an evolutionarily conserved family consisting of essential pre-mRNA splicing factors. Since their discovery and initial characterization, roles of SR proteins in pre-mRNA splicing and in subsequent steps of post-transcriptional gene expression have expanded significantly. The current hypotheses suggest that SR proteins are multifunctional adaptor molecules that may couple distinct steps of RNA metabolism. In the present study, we will provide an overview of the roles of SR proteins in different steps of post-transcriptional gene expression.

A novel role for shuttling SR proteins in mRNA translation

The Ser-Arg-rich (SR) proteins comprise a large family of nuclear phosphoproteins that are required for constitutive and alternative splicing. A subset of SR proteins shuttles continuously between the nucleus and the cytoplasm, suggesting that the role of shuttling SR proteins in gene expression may not be limited to nuclear pre-mRNA splicing, but may also include unknown cytoplasmic functions. Here, we show that shuttling SR proteins, in particular SF2/ASF, associate with translating ribosomes and stimulate translation when tethered to a reporter mRNA in Xenopus oocytes. Moreover, SF2/ASF enhances translation of reporter mRNAs in HeLa cells, and this activity is dependent on its ability to shuttle from the nucleus to the cytoplasm and is increased by the presence of an exonic-splicing enhancer. Furthermore, SF2/ASF can stimulate translation in vitro using a HeLa cell-free translation system. Thus, the association of SR proteins with translating ribosomes, as well as the stimulation of translation both in vivo and in vitro, strongly suggest a role for shuttling SR proteins in translation. We propose that shuttling SR proteins play multiple roles in the posttranscriptional expression of eukaryotic genes and illustrate how they may couple splicing and translation.

The Wilms' tumour protein (WT1) shuttles between nucleus and cytoplasm and is present in functional polysomes

Mutations of the Wilms' tumour-1 (WT1) gene in humans can lead to childhood kidney cancer, life-threatening glomerular nephropathy and gonadal dysgenesis. The WT1 protein is normally expressed in the developing genitourinary tract, heart, spleen and adrenal glands and is crucial for their development, however it's function at the molecular level is yet to be fully understood. The protein is predominantly nuclear and there is evidence that the two different isoforms of WT1 (-KTS and +KTS) are involved in two different steps of gene expression control: transcription and RNA processing. In this study we report a novel property of WT1, namely that it shuttles between the nucleus and cytoplasm. Moreover, western blot analysis showed that between 10 and 50% of total cellular WT1 can be detected in the cytoplasm depending on the cell type. A significant proportion of cytoplasmic WT1 is in association with ribonucleoprotein particles (RNPs), which strengthens the idea of its involvement in RNA metabolism. Furthermore, we report that WT1 is associated with actively translating polysomes, extending even further the potential roles of WT1 and opening the possibility that it is involved in the regulation of translation. Interestingly, despite the functional differences between two of the WT1 isoforms (+/-KTS) within the nucleus, both isoforms share the shuttling property and are found in translating polysomes.

Regulation of SR protein localization during development

Ser-Arg-rich (SR) proteins play numerous roles in spliceosome assembly and the regulation of splice-site selection. Whereas considerable attention has focused on the mechanistic details of SR protein activities, little is known concerning how these splicing regulators are controlled by the cell. Here we examined the subcellular localization of precursor mRNA splicing factors during early development of the nematode Ascaris lumbricoides. In the early embryo, before major zygotic gene activation, most SR proteins, along with RNA polymerase II, are localized in the cytoplasm. As development proceeds, we observe a significant decrease in the cytoplasmic levels of these factors and a concomitant increase in nuclear localization. In contrast, trimethylguanosine-capped small nuclear ribonucleoproteins are predominantly localized in the nucleus throughout this period. We previously showed that the phosphorylation state and activity of SR proteins are regulated during A. lumbricoides embryogenesis. These changes correlate with the onset of precursor mRNA splicing and zygotic transcription. Thus, a coordinate change in the subcellular localization of SR proteins and RNA polymerase II occurs at the transition from reliance on maternally deposited factors to embryonic expression. We propose that before zygotic gene activation, SR proteins and RNA polymerase II are stockpiled in the cytoplasm of early embryos, awaiting signals that lead to their activation.

SR proteins are required for nematode trans-splicing in vitro

SR (ser/arg) proteins have been shown to play roles in numerous aspects of pre-mRNA splicing, including modulation of alternative splicing, commitment of substrates to the splicing pathway, and splice site communication. The last of these, splice site communication, is particularly relevant to trans-splicing in which the 5' and 3' exons originate on separate molecules. The participation of SR proteins in naturally occurring, spliced leader RNA-dependent transsplicing has not been examined. Here, we have isolated SR proteins from an organism that performs both trans- and cis-splicing, the nematode Ascaris lumbricoides. To examine their activity in in vitro splicing reactions, we have also developed and characterized an SR protein-depleted whole-cell extract. When tested in this extract, the nematode SR proteins are required for both trans- and cis-splicing. In addition, the state of phosphorylation of the nematode SR proteins is critical to their activity in vitro. Interestingly, mammalian (HeLa) and A. lumbricoides SR proteins exhibit equivalent activities in cis-splicing, while the nematode SR proteins are much more active in trans-splicing. Thus, it appears that SR proteins purified from an organism that naturally trans-splices its pre-mRNAs promote this reaction to a greater extent than do their mammalian counterparts.

Developmental regulation of SR protein phosphorylation and activity

Serine/arginine-rich splicing factors (SR proteins) are substrates for serine phosphorylation that can regulate SR protein function. We have observed gross changes in SR protein phosphorylation during early development coincident with major zygotic gene activation in the nematode Ascaris lumbricoides. These differences correlate with large-scale changes in SR protein activity in promoting both trans- and cis-splicing. Importantly, inactive early stage extracts can be made splicing competent on addition of later stage SR proteins. These data suggest that changes in SR protein phosphorylation have a role in the activation of pre-mRNA splicing during early development.

A new approach to tuberculosis notification

A tuberculosis notification system was initiated which used several sources--the doctor, the diagnostic services of bacteriology, histopathology and radiology, and the pharmacy. The results of six months' notifications have shown that the pharmacy was by far the most complete source (74% of all new cases); doctors notified only 21% of all new cases. The combination of both doctors and the pharmacy was the most complete pair, notifying 85% of all new cases. It is suggested that consideration should be given to the use of pharmacies as a source of notification both in other settings and for other diseases.

Tolerance of debility in elderly dependants by supporters at home: its significance for hospital practice

Some 12% of all geriatric admissions to University College Hospital and Whittington Hospital are for patients whose relatives or friends can no longer cope with them at home. The person principally involved with home support was interviewed in 50 such cases. The causes of inability to cope were identified on a quantitative and qualitative basis. The supporters were asked to assess which of the problems identified would have to be alleviated to restore a tolerable situation at home; 46 (92%) were able to do so. Identification of the "alleviation factors" forms a therapeutic and prognostic guideline in this type of admission and may have far-reaching social and economic implications.