Search

Two-dimensional high-throughput data have become increasingly common in functional genomics studies, which raises new challenges in data analysis. Here, we introduce a new statistic called Zeta, initially developed to identify global splicing regulators from a two-dimensional RNAi screen, a high-throughput screen coupled with high-throughput functional readouts, and ZetaSuite, a software package to facilitate general application of the Zeta statistics. We compare our approach with existing methods using multiple benchmarked datasets and then demonstrate the broad utility of ZetaSuite in processing public data from large-scale cancer dependency screens and single-cell transcriptomics studies to elucidate novel biological insights.

Mammals have limited regenerative capacity, whereas some vertebrates, like fish and salamanders, are able to regenerate their organs efficiently. The regeneration in these species depends on cell dedifferentiation followed by proliferation. We generate a mouse model that enables the inducible expression of the four Yamanaka factors (Oct-3/4, Sox2, Klf4, and c-Myc, or 4F) specifically in hepatocytes. Transient in vivo 4F expression induces partial reprogramming of adult hepatocytes to a progenitor state and concomitantly increases cell proliferation. This is indicated by reduced expression of differentiated hepatic-lineage markers, an increase in markers of proliferation and chromatin modifiers, global changes in DNA accessibility, and an acquisition of liver stem and progenitor cell markers. Functionally, short-term expression of 4F enhances liver regenerative capacity through topoisomerase2-mediated partial reprogramming. Our results reveal that liver-specific 4F expression in vivo induces cellular plasticity and counteracts liver failure, suggesting that partial reprogramming may represent an avenue for enhancing tissue regeneration.

Glioblastoma (GBM) is the most common and lethal primary malignant brain tumor, containing GBM stem cells (GSCs) that contribute to therapeutic resistance and relapse. Exposing potential GSC vulnerabilities may provide therapeutic strategies against GBM. Here, we interrogated the role of adenosine-to-inosine (A-to-I) RNA editing mediated by adenosine deaminase acting on RNA 1 (ADAR1) in GSCs and found that both ADAR1 and global RNA editomes were elevated in GSCs compared with normal neural stem cells. ADAR1 inactivation or blocking of the upstream JAK/STAT pathway through TYK2 inhibition impaired GSC self-renewal and stemness. Downstream of ADAR1, RNA editing of the 3'-UTR of GM2A, a key ganglioside catabolism activator, proved to be critical, as interference with ganglioside catabolism and disruption of ADAR1 showed a similar functional impact on GSCs. These findings reveal that RNA editing links ganglioside catabolism to GSC self-renewal and stemness, exposing a potential vulnerability of GBM for therapeutic intervention.

The diverse functions of WASP, the deficiency of which causes Wiskott-Aldrich syndrome (WAS), remain poorly defined. We generated three isogenic WAS models using patient induced pluripotent stem cells and genome editing. These models recapitulated WAS phenotypes and revealed that WASP deficiency causes an upregulation of numerous RNA splicing factors and widespread altered splicing. Loss of WASP binding to splicing factor gene promoters frequently leads to aberrant epigenetic activation. WASP interacts with dozens of nuclear speckle constituents and constrains SRSF2 mobility. Using an optogenetic system, we showed that WASP forms phase-separated condensates that encompasses SRSF2, nascent RNA and active Pol II. The role of WASP in gene body condensates is corroborated by ChIPseq and RIPseq. Together our data reveal that WASP is a nexus regulator of RNA splicing that controls the transcription of splicing factors epigenetically and the dynamics of the splicing machinery through liquid-liquid phase separation.

Breast cancer is the most frequently occurred cancer type and the second cause of death in women worldwide. Alternative splicing (AS) is the process that generates more than one mRNA isoform from a single gene, and it plays a major role in expanding the human protein diversity. Aberrant AS contributes to breast cancer metastasis and resistance to chemotherapeutic interventions. Therefore, identifying cancer-specific isoforms is the prerequisite for therapeutic interventions intended to correct aberrantly expressed AS events. Here, we performed RNA-mediated oligonucleotide annealing, selection, and ligation coupled with next-generation sequencing (RASL-seq) in breast cancer cells, to identify global breast cancer-specific AS defects. By RT-PCR validation, we demonstrate the high accuracy of RASL-seq results. In addition, we analyzed identified AS events using the Cancer Genome Atlas (TCGA) database in a large number of non-pathological and breast tumor specimens and validated them in normal and breast cancer samples. Interestingly, aberrantly regulated AS cassette exons in cancer tissues do not encode for known functional domains but instead encode for amino acids constituting regions of intrinsically disordered protein portions characterized by high flexibility and prone to be subjected to post-translational modifications. Collectively, our results reveal novel AS errors occurring in human breast cancer, potentially affecting breast cancer-related biological processes.

Aberrant alternative splicing (AS) is a hallmark of cancer and a potential target for novel anti-cancer therapeutics. Breast cancer-associated AS events are known to be linked to disease progression, metastasis, and survival of breast cancer patients. To identify altered AS programs occurring in metastatic breast cancer, we perform a global analysis of AS events by using RNA-mediated oligonucleotide annealing, selection, and ligation coupled with next-generation sequencing (RASL-seq). We demonstrate that, relative to low-metastatic, high-metastatic breast cancer cells show different AS choices in genes related to cancer progression. Supporting a global reshape of cancer-related splicing profiles in metastatic breast cancer we found an enrichment of RNA-binding motifs recognized by several splicing regulators, which have aberrant expression levels or activity during breast cancer progression, including SRSF1. Among SRSF1-regulated targets we found DCUN1D5, a gene for which skipping of exon 4 in its pre-mRNA introduces a premature termination codon (PTC), thus generating an unstable transcript degraded by nonsense-mediated mRNA decay (NMD). Significantly, distinct breast cancer subtypes show different DCUN1D5 isoform ratios with metastatic breast cancer expressing the highest level of the NMD-insensitive DCUN1D5 mRNA, thus showing high DCUN1D5 expression levels, which are ultimately associated with poor overall and relapse-free survival in breast cancer patients. Collectively, our results reveal global AS features of metastatic breast tumors, which open new possibilities for the treatment of these aggressive tumor types.

The paternal genome undergoes a massive exchange of histone with protamine for compaction into sperm during spermiogenesis. Upon fertilization, this process is potently reversed, which is essential for parental genome reprogramming and subsequent activation; however, it remains poorly understood how this fundamental process is initiated and regulated. Here, we report that the previously characterized splicing kinase SRPK1 initiates this life-beginning event by catalyzing site-specific phosphorylation of protamine, thereby triggering protamine-to-histone exchange in the fertilized oocyte. Interestingly, protamine undergoes a DNA-dependent phase transition to gel-like condensates and SRPK1-mediated phosphorylation likely helps open up such structures to enhance protamine dismissal by nucleoplasmin (NPM2) and enable the recruitment of HIRA for H3.3 deposition. Remarkably, genome-wide assay for transposase-accessible chromatin sequencing (ATAC-seq) analysis reveals that selective chromatin accessibility in both sperm and MII oocytes is largely erased in early pronuclei in a protamine phosphorylation-dependent manner, suggesting that SRPK1-catalyzed phosphorylation initiates a highly synchronized reorganization program in both parental genomes.

Heart performance relies on highly coordinated excitation-contraction (EC) coupling, and defects in this critical process may be exacerbated by additional genetic defects and/or environmental insults to cause eventual heart failure. Here we report a regulatory pathway consisting of the RNA binding protein RBFox2, a stress-induced microRNA miR-34a, and the essential EC coupler JPH2. In this pathway, initial cardiac defects diminish RBFox2 expression, which induces transcriptional repression of miR-34a, and elevated miR-34a targets Jph2 to impair EC coupling, which further manifests heart dysfunction, leading to progressive heart failure. The key contribution of miR-34a to this process is further established by administrating its mimic, which is sufficient to induce cardiac defects, and by using its antagomir to alleviate RBFox2 depletion-induced heart dysfunction. These findings elucidate a potential feed-forward mechanism to account for a critical transition to cardiac decompensation and suggest a potential therapeutic avenue against heart failure.

Mutations in several general pre-mRNA splicing factors have been linked to myelodysplastic syndromes (MDSs) and solid tumors. These mutations have generally been assumed to cause disease by the resultant splicing defects, but different mutations appear to induce distinct splicing defects, raising the possibility that an alternative common mechanism is involved. Here we report a chain of events triggered by multiple splicing factor mutations, especially high-risk alleles in SRSF2 and U2AF1, including elevated R-loops, replication stress, and activation of the ataxia telangiectasia and Rad3-related protein (ATR)-Chk1 pathway. We further demonstrate that enhanced R-loops, opposite to the expectation from gained RNA binding with mutant SRSF2, result from impaired transcription pause release because the mutant protein loses its ability to extract the RNA polymerase II (Pol II) C-terminal domain (CTD) kinase-the positive transcription elongation factor complex (P-TEFb)-from the 7SK complex. Enhanced R-loops are linked to compromised proliferation of bone-marrow-derived blood progenitors, which can be partially rescued by RNase H overexpression, suggesting a direct contribution of augmented R-loops to the MDS phenotype.

R-loop, a three-stranded RNA/DNA structure, has been linked to induced genome instability and regulated gene expression. To enable precision analysis of R-loops in vivo, we develop an RNase-H-based approach; this reveals predominant R-loop formation near gene promoters with strong G/C skew and propensity to form G-quadruplex in non-template DNA, corroborating with all biochemically established properties of R-loops. Transcription perturbation experiments further indicate that R-loop induction correlates to transcriptional pausing. Interestingly, we note that most mapped R-loops are each linked to a nearby free RNA end; by using a ribozyme to co-transcriptionally cleave nascent RNA, we demonstrate that such a free RNA end coupled with a G/C-skewed sequence is necessary and sufficient to induce R-loop. These findings provide a topological solution for RNA invasion into duplex DNA and suggest an order for R-loop initiation and elongation in an opposite direction to that previously proposed.

MicroRNA (miRNA) biogenesis is known to be modulated by a variety of RNA-binding proteins (RBPs), but in most cases, individual RBPs appear to influence the processing of a small subset of target miRNAs. Here, we report that the RNA-binding NONO-PSF heterodimer binds a large number of expressed pri-miRNAs in HeLa cells to globally enhance pri-miRNA processing by the Drosha-DGCR8 Microprocessor. NONO and PSF are key components of paraspeckles organized by the long noncoding RNA (lncRNA) NEAT1. We further demonstrate that NEAT1 also has a profound effect on global pri-miRNA processing. Mechanistic dissection reveals that NEAT1 broadly interacts with the NONO-PSF heterodimer as well as many other RBPs and that multiple RNA segments in NEAT1, including a 'pseudo pri-miRNA' near its 3' end, help attract the Microprocessor. These findings suggest a 'bird nest' model in which an lncRNA orchestrates efficient processing of potentially an entire class of small noncoding RNAs in the nucleus.

Trimethylation of histone H3K36 is a chromatin mark associated with active gene expression, which has been implicated in coupling transcription with mRNA splicing and DNA damage response. SETD2 is a major H3K36 trimethyltransferase, which has been implicated as a tumor suppressor in mammals. Here, we report the regulation of SETD2 protein stability by the proteasome system, and the identification of SPOP, a key subunit of the CUL3 ubiquitin E3 ligase complex, as a SETD2-interacting protein. We demonstrate that SPOP is critically involved in SETD2 stability control and that the SPOP/CUL3 complex is responsible for SETD2 polyubiquitination both in vivo and in vitro ChIP-Seq analysis and biochemical experiments demonstrate that modulation of SPOP expression confers differential H3K36me3 on SETD2 target genes, and induce H3K36me3-coupled alternative splicing events. Together, these findings establish a functional connection between oncogenic SPOP and tumor suppressive SETD2 in the dynamic regulation of gene expression on chromatin.

The utility of the mitochondrial genomes (mitogenomes) in analyzing the evolutionary history of animals has been proven. Five deep-sea corals (Bathypathes sp.1, Bathypathes sp.2, Schizopathidae 1, Trissopathes sp., and Leiopathes sp.) were collected in the South China Sea (SCS). Initially, the structures and collinearity of the five deep-sea coral mitogenomes were analyzed. The gene arrangements in the five deep-sea coral mitogenomes were similar to those in the order Antipatharia, which evidenced their conservation throughout evolutionary history. Additionally, to elucidate the slow evolutionary rates in Hexacorallia mitogenomes, we conducted comprehensive analyses, including examining phylogenetic relationships, performing average nucleotide identity (ANI) analysis, and assessing GC-skew dissimilarity combining five deep-sea coral mitogenomes and 522 reference Hexacorallia mitogenomes. Phylogenetic analysis using 13 conserved proteins revealed that species clustered together at the order level, and they exhibited interspersed distributions at the family level. The ANI results revealed that species had significant similarities (identity > 85%) within the same order, while species from different orders showed notable differences (identity < 80%). The investigation of the Hexacorallia mitogenomes also highlighted that the GC-skew dissimilarity was highly significant at the order level, but not as pronounced at the family level. These results might be attributed to the slow evolution rate of Hexacorallia mitogenomes and provide evidence of mitogenomic diversity. Furthermore, divergence time analysis revealed older divergence times assessed via mitogenomes compared with nuclear data, shedding light on significant evolutionary events shaping distinct orders within Hexacorallia corals. Those findings provide new insights into understanding the slow evolutionary rates of deep-sea corals in all lineages of Hexacorallia using their mitogenomes. [ABSTRACT FROM AUTHOR]

DNA methylation is a key epigenetic mechanism orchestrating gene expression networks in many biological processes. Nonetheless, studying the role of specific gene methylation events in fish faces challenges. In this study, we validate the regulation of DNA methylation on empty spiracles homeobox 2 (emx2) expression with decitabine treatment in Chinese tongue sole testis cells. We used the emx2 gene as the target gene and developed a new DNA methylation editing system by fusing dnmt3a with catalytic dead Cas9 (dCas9) and demonstrated its ability for sequence-specific DNA methylation editing. Results revealed that utilizing dCas9-dnmt3a to target emx2 promoter region led to increased DNA methylation levels and decreased emx2 expression in Chinese tongue sole testis cells. More importantly, the DNA methylation editing significantly suppressed the expression of MYC proto-oncogene, bHLH transcription factor (myc), one target gene of emx2. Furthermore, we assessed the off-target effects of dCas9-dnmt3a and confirmed no significant impact on the predicted off-target gene expression. Taken together, we developed the first DNA methylation editing system in marine species and demonstrated its effective editing ability in Chinese tongue sole cells. This provides a new strategy for both epigenetic research and molecular breeding of marine species. [ABSTRACT FROM AUTHOR]

Proteasome 26S Subunit, Non-ATPase 9 (psmd9) plays an important role in the balance of protamine and the stability of the nucleolar structure during spermatogenesis. In this study, we cloned the psmd9 of Cynoglossus semilaevis and analyzed its expression pattern. psmd9 was identified on the Z chromosome of C. semilaevis, which is considered an interesting candidate gene for spermatogenesis. qRT-PCR and FISH experiments showed that the psmd9 gene was significantly highly expressed in the testes. It is worth noting that the expression level of psmd9 in male fish testes is significantly higher than that in pseudomales. In order to further explore the role of psmd9 in spermatogenesis, a male testicular cell line was used as the experimental material. The results of the psmd9-RNAi and overexpression experiments showed that psmd9 had a synergistic effect with spermatogenesis-related genes dnd1, cfap69, dnah3 and dnajb13, but had an antagonistic effect with ccne2. Our findings offer a scientific foundation for comprehending the role of psmd9 in the spermatogenesis regulatory network of C. semilaevis. [ABSTRACT FROM AUTHOR]

Myelodysplastic syndromes (MDS) are heterogeneous myeloid disorders with prevalent mutations in several splicing factors, but the splicing programs linked to specific mutations or MDS in general remain to be systematically defined. We applied RASL-seq, a sensitive and cost-effective platform, to interrogate 5502 annotated splicing events in 169 samples from MDS patients or healthy individuals. We found that splicing signatures associated with normal hematopoietic lineages are largely related to cell signaling and differentiation programs, whereas MDS-linked signatures are primarily involved in cell cycle control and DNA damage responses. Despite the shared roles of affected splicing factors in the 3 splice site definition, mutations in U2AF1, SRSF2, and SF3B1 affect divergent splicing programs, and interestingly, the affected genes fall into converging cancer-related pathways. A risk score derived from 11 splicing events appears to be independently associated with an MDS prognosis and AML transformation, suggesting potential clinical relevance of altered splicing patterns in MDS.

The U2AF heterodimer has been well studied for its role in defining functional 3' splice sites in pre-mRNA splicing, but many fundamental questions still remain unaddressed regarding the function of U2AF in mammalian genomes. Through genome-wide analysis of U2AF-RNA interactions, we report that U2AF has the capacity to directly define ~88% of functional 3' splice sites in the human genome, but numerous U2AF binding events also occur in intronic locations. Mechanistic dissection reveals that upstream intronic binding events interfere with the immediate downstream 3' splice site associated either with the alternative exon, to cause exon skipping, or with the competing constitutive exon, to induce exon inclusion. We further demonstrate partial functional impairment with leukemia-associated mutations in U2AF35, but not U2AF65, in regulated splicing. These findings reveal the genomic function and regulatory mechanism of U2AF in both normal and disease states.