Extensive cross-regulation of post-transcriptional regulatory networks in Drosophila

Novel roles for RNA binding proteins squid, hephaesteus, and Hrb27C in Drosophila oogenesis

BACKGROUND

Reproductive capacity in many organisms is maintained by germline stem cells (GSCs). A complex regulatory network influences stem cell fate, including intrinsic factors, local signals, and hormonal and nutritional cues. Posttranscriptional regulatory mechanisms ensure proper cell fate transitions, promoting germ cell differentiation to oocytes. As essential RNA binding proteins with constitutive functions in RNA metabolism, heterogeneous nuclear ribonucleoproteins (hnRNPs) have been implicated in GSC function and axis specification during oocyte development. HnRNPs support biogenesis, localization, maturation, and translation of nascent transcripts. Whether and individual hnRNPs specifically regulate GSC function has yet to be explored.

RESULTS

We demonstrate that hnRNPs are expressed in distinct patterns in the Drosophila germarium. We show that three hnRNPs, squid, hephaestus, and Hrb27C are cell-autonomously required in GSCs for their maintenance. Although these hnRNPs do not impact adhesion of GSCs to adjacent cap cells, squid and hephaestus (but not Hrb27C) are necessary for proper bone morphogenetic protein signaling in GSCs. Moreover, Hrb27C promotes proper GSC proliferation, whereas hephaestus promotes cyst division.

CONCLUSIONS

We find that hnRNPs are independently and intrinsically required in GSCs for their maintenance in adults. Our results support the model that hnRNPs play unique roles in stem cells essential for their self-renewal and proliferation.

Parallel evolution of a splicing program controlling neuronal excitability in flies and mammals

Alternative splicing increases neuronal transcriptomic complexity throughout animal phylogeny. To delve into the mechanisms controlling the assembly and evolution of this regulatory layer, we characterized the neuronal microexon program in Drosophila and compared it with that of mammals. In nonvertebrate bilaterians, this splicing program is restricted to neurons by the posttranscriptional processing of the enhancer of microexons (eMIC) domain in Srrm234. In Drosophila, this processing is dependent on regulation by Elav/Fne. eMIC deficiency or misexpression leads to widespread neurological alterations largely emerging from impaired neuronal activity, as revealed by a combination of neuronal imaging experiments and cell type-specific rescues. These defects are associated with the genome-wide skipping of short neural exons, which are strongly enriched in ion channels. We found no overlap of eMIC-regulated exons between flies and mice, illustrating how ancient posttranscriptional programs can evolve independently in different phyla to affect distinct cellular modules while maintaining cell-type specificity.

An "Omic" Overview of Fragile X Syndrome

Fragile X syndrome (FXS) is a neurodevelopmental disorder associated with a wide range of cognitive, behavioral and medical problems. It arises from the silencing of the fragile X mental retardation 1 (FMR1) gene and, consequently, in the absence of its encoded protein, FMRP (fragile X mental retardation protein). FMRP is a ubiquitously expressed and multifunctional RNA-binding protein, primarily considered as a translational regulator. Pre-clinical studies of the past two decades have therefore focused on this function to relate FMRP's absence to the molecular mechanisms underlying FXS physiopathology. Based on these data, successful pharmacological strategies were developed to rescue fragile X phenotype in animal models. Unfortunately, these results did not translate into humans as clinical trials using same therapeutic approaches did not reach the expected outcomes. These failures highlight the need to put into perspective the different functions of FMRP in order to get a more comprehensive understanding of FXS pathophysiology. This work presents a review of FMRP's involvement on noteworthy molecular mechanisms that may ultimately contribute to various biochemical alterations composing the fragile X phenotype.

RNAct: Protein-RNA interaction predictions for model organisms with supporting experimental data

Protein-RNA interactions are implicated in a number of physiological roles as well as diseases, with molecular mechanisms ranging from defects in RNA splicing, localization and translation to the formation of aggregates. Currently, ∼1400 human proteins have experimental evidence of RNA-binding activity. However, only ∼250 of these proteins currently have experimental data on their target RNAs from various sequencing-based methods such as eCLIP. To bridge this gap, we used an established, computationally expensive protein-RNA interaction prediction method, catRAPID, to populate a large database, RNAct. RNAct allows easy lookup of known and predicted interactions and enables global views of the human, mouse and yeast protein-RNA interactomes, expanding them in a genome-wide manner far beyond experimental data (http://rnact.crg.eu).

Neuronal upregulation of Prospero protein is driven by alternative mRNA polyadenylation and Syncrip-mediated mRNA stabilisation

During Drosophila and vertebrate brain development, the conserved transcription factor Prospero/Prox1 is an important regulator of the transition between proliferation and differentiation. Prospero level is low in neural stem cells and their immediate progeny, but is upregulated in larval neurons and it is unknown how this process is controlled. Here, we use single molecule fluorescent in situ hybridisation to show that larval neurons selectively transcribe a long prospero mRNA isoform containing a 15 kb 3' untranslated region, which is bound in the brain by the conserved RNA-binding protein Syncrip/hnRNPQ. Syncrip binding increases the stability of the long prospero mRNA isoform, which allows an upregulation of Prospero protein production. Adult flies selectively lacking the long prospero isoform show abnormal behaviour that could result from impaired locomotor or neurological activity. Our findings highlight a regulatory strategy involving alternative polyadenylation followed by differential post-transcriptional regulation.This article has an associated First Person interview with the first author of the paper.

UPF1 regulates the malignant biological behaviors of glioblastoma cells via enhancing the stability of Linc-00313

There is growing evidence that the long non-coding RNAs(lncRNAs) play an important role in the biological behaviors of glioblastoma cells. In this study, we elucidated the function and possible effect and molecular mechanisms of lncRNA-Linc-00313 on the biological behaviors of glioblastoma cells as well as UPF1 function as a RNA-binding protein to enhance its stability. Here, we used qRT-PCR and western blot to measure the expression, cell Transfection to disrupt the expression of genes, cell viability analysis, quantization of apoptosis, cell migration, and invasion assays, Reporter vectors construction and luciferase assays to investigate the malignant biological behaviors of cells, human lncRNA microarrays, RNA-Immunoprecipitation, dual-luciferase gene reporter assay, half-life assay and chromatin immunoprecipitation to verify the binding sites, tumor xenograft implantation for in vivo experiment, SPSS 18.0 statistical software for data statistics. UPF1 and Linc-00313 were both upregulated in glioma tissues and cells. Knockdown of UPF1 or Linc-00313 significantly inhibited malignant biological behaviors of glioma cells by regulating miR-342-3p and miR-485-5p, which are downregulated and functioned as tumor suppressors in glioma. Furthermore, Linc-00313 could acted as a competing endogenous RNA(ceRNA) to regulate the expression of Zic4 by binding to miR-342-3p and miR-485-5p. Interestingly, Zic4 could bind to the promoters of UPF1 and Linc-00313 respectively and upregulate the expression of them. These results indicated that a positive-feedback loop was formed in the regulation of the biological behaviors of glioma cells. The study is the first to prove that the UPF1-Linc-00313-miR-342-3p/miR-485-5p-Zic4-SHCBP1 pathway forms a positive-feedback loop and regulates the biological behaviors of U87 and U251 cells, which might provide a new therapeutic target for glioma.

Cluster analysis of replicated alternative polyadenylation data using canonical correlation analysis

BACKGROUND

Alternative polyadenylation (APA) has emerged as a pervasive mechanism that contributes to the transcriptome complexity and dynamics of gene regulation. The current tsunami of whole genome poly(A) site data from various conditions generated by 3' end sequencing provides a valuable data source for the study of APA-related gene expression. Cluster analysis is a powerful technique for investigating the association structure among genes, however, conventional gene clustering methods are not suitable for APA-related data as they fail to consider the information of poly(A) sites (e.g., location, abundance, number, etc.) within each gene or measure the association among poly(A) sites between two genes.

RESULTS

Here we proposed a computational framework, named PASCCA, for clustering genes from replicated or unreplicated poly(A) site data using canonical correlation analysis (CCA). PASCCA incorporates multiple layers of gene expression data from both the poly(A) site level and gene level and takes into account the number of replicates and the variability within each experimental group. Moreover, PASCCA characterizes poly(A) sites in various ways including the abundance and relative usage, which can exploit the advantages of 3' end deep sequencing in quantifying APA sites. Using both real and synthetic poly(A) site data sets, the cluster analysis demonstrates that PASCCA outperforms other widely-used distance measures under five performance metrics including connectivity, the Dunn index, average distance, average distance between means, and the biological homogeneity index. We also used PASCCA to infer APA-specific gene modules from recently published poly(A) site data of rice and discovered some distinct functional gene modules. We have made PASCCA an easy-to-use R package for APA-related gene expression analyses, including the characterization of poly(A) sites, quantification of association between genes, and clustering of genes.

CONCLUSIONS

By providing a better treatment of the noise inherent in repeated measurements and taking into account multiple layers of poly(A) site data, PASCCA could be a general tool for clustering and analyzing APA-specific gene expression data. PASCCA could be used to elucidate the dynamic interplay of genes and their APA sites among various biological conditions from emerging 3' end sequencing data to address the complex biological phenomenon.

Striking circadian neuron diversity and cycling of Drosophila alternative splicing

Although alternative pre-mRNA splicing (AS) significantly diversifies the neuronal proteome, the extent of AS is still unknown due in part to the large number of diverse cell types in the brain. To address this complexity issue, we used an annotation-free computational method to analyze and compare the AS profiles between small specific groups of Drosophila circadian neurons. The method, the Junction Usage Model (JUM), allows the comprehensive profiling of both known and novel AS events from specific RNA-seq libraries. The results show that many diverse and novel pre-mRNA isoforms are preferentially expressed in one class of clock neuron and also absent from the more standard Drosophila head RNA preparation. These AS events are enriched in potassium channels important for neuronal firing, and there are also cycling isoforms with no detectable underlying transcriptional oscillations. The results suggest massive AS regulation in the brain that is also likely important for circadian regulation.

The Fragile X Protein and Genome Function

The fragile X syndrome (FXS) arises from loss of expression or function of the FMR1 gene and is one of the most common monogenic forms of intellectual disability and autism. During the past two decades of FXS research, the fragile X mental retardation protein (FMRP) has been primarily characterized as a cytoplasmic RNA binding protein that facilitates transport of select RNA substrates through neural projections and regulation of translation within synaptic compartments, with the protein products of such mRNAs then modulating cognitive functions. However, the presence of a small fraction of FMRP in the nucleus has long been recognized. Accordingly, recent studies have uncovered several mechanisms or pathways by which FMRP influences nuclear gene expression and genome function. Some of these pathways appear to be independent of the classical role for FMRP as a regulator of translation and point to novel functions, including the possibility that FMRP directly participates in the DNA damage response and in the maintenance of genome stability. In this review, we highlight these advances and discuss how these new findings could contribute to our understanding of FMRP in brain development and function, the neural pathology of fragile X syndrome, and perhaps impact of future therapeutic considerations.

Trans-acting translational regulatory RNA binding proteins

The canonical molecular machinery required for global mRNA translation and its control has been well defined, with distinct sets of proteins involved in the processes of translation initiation, elongation and termination. Additionally, noncanonical, trans-acting regulatory RNA-binding proteins (RBPs) are necessary to provide mRNA-specific translation, and these interact with 5' and 3' untranslated regions and coding regions of mRNA to regulate ribosome recruitment and transit. Recently it has also been demonstrated that trans-acting ribosomal proteins direct the translation of specific mRNAs. Importantly, it has been shown that subsets of RBPs often work in concert, forming distinct regulatory complexes upon different cellular perturbation, creating an RBP combinatorial code, which through the translation of specific subsets of mRNAs, dictate cell fate. With the development of new methodologies, a plethora of novel RNA binding proteins have recently been identified, although the function of many of these proteins within mRNA translation is unknown. In this review we will discuss these methodologies and their shortcomings when applied to the study of translation, which need to be addressed to enable a better understanding of trans-acting translational regulatory proteins. Moreover, we discuss the protein domains that are responsible for RNA binding as well as the RNA motifs to which they bind, and the role of trans-acting ribosomal proteins in directing the translation of specific mRNAs. This article is categorized under: RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes Translation > Translation Regulation Translation > Translation Mechanisms.

Iterative random forests to discover predictive and stable high-order interactions

We developed a predictive, stable, and interpretable tool: the iterative random forest algorithm (iRF). iRF discovers high-order interactions among biomolecules with the same order of computational cost as random forests. We demonstrate the efficacy of iRF by finding known and promising interactions among biomolecules, of up to fifth and sixth order, in two data examples in transcriptional regulation and alternative splicing.

Genomics has revolutionized biology, enabling the interrogation of whole transcriptomes, genome-wide binding sites for proteins, and many other molecular processes. However, individual genomic assays measure elements that interact in vivo as components of larger molecular machines. Understanding how these high-order interactions drive gene expression presents a substantial statistical challenge. Building on random forests (RFs) and random intersection trees (RITs) and through extensive, biologically inspired simulations, we developed the iterative random forest algorithm (iRF). iRF trains a feature-weighted ensemble of decision trees to detect stable, high-order interactions with the same order of computational cost as the RF. We demonstrate the utility of iRF for high-order interaction discovery in two prediction problems: enhancer activity in the early Drosophila embryo and alternative splicing of primary transcripts in human-derived cell lines. In Drosophila, among the 20 pairwise transcription factor interactions iRF identifies as stable (returned in more than half of bootstrap replicates), 80% have been previously reported as physical interactions. Moreover, third-order interactions, e.g., between Zelda (Zld), Giant (Gt), and Twist (Twi), suggest high-order relationships that are candidates for follow-up experiments. In human-derived cells, iRF rediscovered a central role of H3K36me3 in chromatin-mediated splicing regulation and identified interesting fifth- and sixth-order interactions, indicative of multivalent nucleosomes with specific roles in splicing regulation. By decoupling the order of interactions from the computational cost of identification, iRF opens additional avenues of inquiry into the molecular mechanisms underlying genome biology.

Comprehensive analysis of nucleocytoplasmic dynamics of mRNA in Drosophila cells

Eukaryotic mRNAs undergo a cycle of transcription, nuclear export, and degradation. A major challenge is to obtain a global, quantitative view of these processes. Here we measured the genome-wide nucleocytoplasmic dynamics of mRNA in Drosophila cells by metabolic labeling in combination with cellular fractionation. By mathematical modeling of these data we determined rates of transcription, export and cytoplasmic decay for 5420 genes. We characterized these kinetic rates and investigated links with mRNA features, RNA-binding proteins (RBPs) and chromatin states. We found prominent correlations between mRNA decay rate and transcript size, while nuclear export rates are linked to the size of the 3'UTR. Transcription, export and decay rates are each associated with distinct spectra of RBPs. Specific classes of genes, such as those encoding cytoplasmic ribosomal proteins, exhibit characteristic combinations of rate constants, suggesting modular control. Binding of splicing factors is associated with faster rates of export, and our data suggest coordinated regulation of nuclear export of specific functional classes of genes. Finally, correlations between rate constants suggest global coordination between the three processes. Our approach provides insights into the genome-wide nucleocytoplasmic kinetics of mRNA and should be generally applicable to other cell systems.

All mRNAs start from production in the nucleus, undergo exportation through nuclear pores and finally are degraded in the cytoplasm. A comprehensive characterization of the kinetic rates of all mRNAs is an important prerequisite for a global understanding of the regulation of the transcriptome and the cell. By conducting a time-series experiment and building a mathematical model, we trace the dynamics of mRNAs from the nucleus to the cytoplasm and determine the rates at each kinetic step at transcriptome-wide level. This information allows us to associate mRNA kinetic rates with a wealth of biological features and made some intriguing discoveries. We show mRNA decay is positively linked to transcript length while mRNA export is negatively linked to the length of the 3' UTR. We show binding of splicing factors is associated with faster rates of mRNA export. We provide evidence for global coordination between nuclear export an decay of mRNA. We show genes sharing specific functions tend to have similar nucleoplasmic kinetics, in which ribosomal proteins possessing special kinetic features exclusively stand out. Altogether, our integrated approach to quantitatively determine the rates of kinetic steps on a gene-by-gene basis provides a blueprint to obtain the global understanding of RNA regulation.

Highly accessible AU-rich regions in 3' untranslated regions are hotspots for binding of regulatory factors

Post-transcriptional regulation is regarded as one of the major processes involved in the regulation of gene expression. It is mainly performed by RNA binding proteins and microRNAs, which target RNAs and typically affect their stability. Recent efforts from the scientific community have aimed at understanding post-transcriptional regulation at a global scale by using high-throughput sequencing techniques such as cross-linking and immunoprecipitation (CLIP), which facilitates identification of binding sites of these regulatory factors. However, the diversity in the experimental procedures and bioinformatics analyses has hindered the integration of multiple datasets and thus limited the development of an integrated view of post-transcriptional regulation. In this work, we have performed a comprehensive analysis of 107 CLIP datasets from 49 different RBPs in HEK293 cells to shed light on the complex interactions that govern post-transcriptional regulation. By developing a more stringent CLIP analysis pipeline we have discovered the existence of conserved regulatory AU-rich regions in the 3’UTRs where miRNAs and RBPs that regulate several processes such as polyadenylation or mRNA stability bind. Analogous to promoters, many factors have binding sites overlapping or in close proximity in these hotspots and hence the regulation of the mRNA may depend on their relative concentrations. This hypothesis is supported by RBP knockdown experiments that alter the relative concentration of RBPs in the cell. Upon AGO2 knockdown (KD), transcripts containing “free” target sites show increased expression levels compared to those containing target sites in hotspots, which suggests that target sites within hotspots are less available for miRNAs to bind. Interestingly, these hotspots appear enriched in genes with regulatory functions such as DNA binding and RNA binding. Taken together, our results suggest that hotspots are functional regulatory elements that define an extra layer of regulation of post-transcriptional regulatory networks.

All the cells in a given organism contain the same genome, yet their phenotype can be very diverse. The vast majority of this diversity arises from the differences in the expression of genes and proteins in them. One of the main mechanisms involved in controlling the protein and mRNA repertoire in cells is post-transcriptional regulation. The recent development of high-throughput sequencing techniques gives us now an unprecedented opportunity to investigate how post-transcriptional regulation works and which are the elements involved in defining the final set of mRNAs and proteins inside cells. In this work, we have performed a comprehensive computational analysis of several post-transcriptional regulators in a commonly used human cell line in order to understand which factors are involved in post-transcriptional regulation and how they coordinate their function. The results of our analysis show that this process is orchestrated around small regions in the mRNAs where many regulators bind and may compete with each other to regulate the mRNAs. The investigation and characterization of these regions gives us insight into the underlying combinatorial control that causes gene expression to differ across cell types and in diseases.

Bombyx mori P-element Somatic Inhibitor (BmPSI) Is a Key Auxiliary Factor for Silkworm Male Sex Determination

Manipulation of sex determination pathways in insects provides the basis for a wide spectrum of strategies to benefit agriculture and public health. Furthermore, insects display a remarkable diversity in the genetic pathways that lead to sex differentiation. The silkworm, Bombyx mori, has been cultivated by humans as a beneficial insect for over two millennia, and more recently as a model system for studying lepidopteran genetics and development. Previous studies have identified the B. mori Fem piRNA as the primary female determining factor and BmMasc as its downstream target, while the genetic scenario for male sex determination was still unclear. In the current study, we exploite the transgenic CRISPR/Cas9 system to generate a comprehensive set of knockout mutations in genes BmSxl, Bmtra2, BmImp, BmImp^M, BmPSI and BmMasc, to investigate their roles in silkworm sex determination. Absence of Bmtra2 results in the complete depletion of Bmdsx transcripts, which is the conserved downstream factor in the sex determination pathway, and induces embryonic lethality. Loss of BmImp or BmImp^M function does not affect the sexual differentiation. Mutations in BmPSI and BmMasc genes affect the splicing of Bmdsx and the female reproductive apparatus appeared in the male external genital. Intriguingly, we identify that BmPSI regulates expression of BmMasc, BmImp^M and Bmdsx, supporting the conclusion that it acts as a key auxiliary factor in silkworm male sex determination.

The sex determination system extremely diverse among organisms including insects in which even each order occupy a different manner of sex determination. The silkworm, Bombyx mori, is a lepidopteran model insect with economic importance. The mechanism of the silkworm sex determination has been in mystery for a long time until a Fem piRNA was identified as the primary female sex determinator recently. However, genetic and phenotypic proofs are urgently needed to fully exploit the mechanism, especially of the male sex determination. In the current study, we provided comprehensively genetic evidences by generating CRISPR/Cas9-mediated knockout mutations for those genes BmSxl, Bmtra2, BmImp, BmImp^M, BmPSI and BmMasc, which were considered to be involved in insect sex determination. The results showed that mutations of BmSxl, BmImp and BmImp^M had no physiological and morphological effects on the sexual development while Bmtra2 depletion caused Bmdsx splicing disappeared and induced embryonic lethality. Importantly, the BmPSI regulates expression of BmMasc, BmImp^M and Bmdsx, supporting the conclusion that it acts as a key auxiliary factor to regulate the male sex determination in the silkworm.

Identification of Genes Mediating Drosophila Follicle Cell Progenitor Differentiation by Screening for Modifiers of GAL4::UAS Variegation

The Drosophila melanogaster ovarian follicle cell lineage provides a powerful system for investigating how epigenetic changes contribute to differentiation. Downstream from an epithelial stem cell, follicle progenitors undergo nine mitotic cell cycles before transitioning to the endocycle and initiating differentiation. During their proliferative phase, follicle progenitors experience Lsd1-dependent changes in epigenetic stability that can be monitored using GAL4::UAS variegation. Eventually, follicle progenitors acquire competence to respond to Delta, a Notch ligand present in the environment, which signals them to cease division and initiate differentiation. The time required to acquire competence determines the duration of mitotic cycling and hence the final number of follicle cells. We carried out a screen for dominant modifiers of variegation spanning nearly 70% of Drosophila euchromatin to identify new genes influencing follicle progenitor epigenetic maturation. The eight genes found include chromatin modifiers, but also cell cycle regulators and transcription factors. Five of the modifier genes accelerate the acquisition of progenitor competence and reduce follicle cell number, however, the other three genes affect follicle cell number in an unexpected manner.

A high-throughput pipeline for the production of synthetic antibodies for analysis of ribonucleoprotein complexes

Post-transcriptional regulation of mRNAs plays an essential role in the control of gene expression. mRNAs are regulated in ribonucleoprotein (RNP) complexes by RNA-binding proteins (RBPs) along with associated protein and noncoding RNA (ncRNA) cofactors. A global understanding of post-transcriptional control in any cell type requires identification of the components of all of its RNP complexes. We have previously shown that these complexes can be purified by immunoprecipitation using anti-RBP synthetic antibodies produced by phage display. To develop the large number of synthetic antibodies required for a global analysis of RNP complex composition, we have established a pipeline that combines (i) a computationally aided strategy for design of antigens located outside of annotated domains, (ii) high-throughput antigen expression and purification in Escherichia coli, and (iii) high-throughput antibody selection and screening. Using this pipeline, we have produced 279 antibodies against 61 different protein components of Drosophila melanogaster RNPs. Together with those produced in our low-throughput efforts, we have a panel of 311 antibodies for 67 RNP complex proteins. Tests of a subset of our antibodies demonstrated that 89% immunoprecipitate their endogenous target from embryo lysate. This panel of antibodies will serve as a resource for global studies of RNP complexes in Drosophila. Furthermore, our high-throughput pipeline permits efficient production of synthetic antibodies against any large set of proteins.

Regulation of alternative splicing in Drosophila by 56 RNA binding proteins

Alternative splicing is regulated by RNA binding proteins (RBPs) that recognize pre-mRNA sequence elements and activate or repress adjacent exons. Here, we used RNA interference and RNA-seq to identify splicing events regulated by 56 Drosophila proteins, some previously unknown to regulate splicing. Nearly all proteins affected alternative first exons, suggesting that RBPs play important roles in first exon choice. Half of the splicing events were regulated by multiple proteins, demonstrating extensive combinatorial regulation. We observed that SR and hnRNP proteins tend to act coordinately with each other, not antagonistically. We also identified a cross-regulatory network where splicing regulators affected the splicing of pre-mRNAs encoding other splicing regulators. This large-scale study substantially enhances our understanding of recent models of splicing regulation and provides a resource of thousands of exons that are regulated by 56 diverse RBPs.