Unveiling circRNA-mediated ceRNA networks in ischemic stroke by integrative analysis of multi-source gene expression profiling.

Background: Circular RNAs (circRNAs) are emerging as potential therapeutic targets for ischemic stroke (IS) due to their regulatory roles in inflammation and apoptosis. This study aimed to develop a comprehensive and robust IS-specific competing endogenous RNA (ceRNA) network to facilitate the identification of novel diagnostic and therapeutic targets. Methods: We integrated expression data from 15 IS studies using the Rank-In algorithm to minimize batch effects. Differentially expressed circRNAs, miRNAs, and mRNAs were identified by comparing IS and control samples. Functional enrichment analysis of differentially expressed circRNA host genes revealed significantly enriched pathways and Gene Ontology (GO) terms relevant to IS pathogenesis. We further predicted miRNA-circRNA and mRNA-miRNA interactions, enabling the construction of a comprehensive ceRNA network to identify circRNA-related genes with diagnostic potential for IS. Results: Integrated analysis revealed 199 differentially expressed circRNAs, 103 miRNAs, and 1736 mRNAs in IS patients. Functional enrichment analysis implicated these molecules in relevant pathways like the neurotrophin signaling pathway and p53 signaling pathway. The constructed circRNA-miRNA-mRNA regulatory network provided insights into potential mechanisms underlying IS. Three circRNA-related genes (RGS2, CDK5R1, and NSF) displayed promising diagnostic potential for IS when combined. Conclusions: We successfully constructed a robust and informative IS-specific ceRNA network by integrating data from diverse sources. This network identified differentially expressed RNAs and revealed enriched pathways potentially involved in IS pathogenesis. Notably, our analysis identified CDK5R1, RGS2, and NSF as potential diagnostic biomarkers for IS. This study sheds light on a circRNA-mediated regulatory network with potential diagnostic and therapeutic implications for ischemic stroke.

PRMT5 orchestrates EGFR and AKT networks to activate NFκB and promote EMT.

Neuroblastoma remains a formidable challenge in pediatric oncology, representing 15% of cancer-related mortalities in children. Despite advancements in combinatorial and targeted treatments improving survival rates, nearly 50% of patients with high-risk neuroblastoma will ultimately succumb to their disease. Dysregulation of the epithelial-mesenchymal transition (EMT) is a key mechanism of tumor cell dissemination, resulting in metastasis and poor outcomes in many cancers. Our prior work identified PRMT5 as a key regulator of EMT via methylation of AKT at arginine 15, enhancing the expression of EMT-driving transcription factors and facilitating metastasis. Here, we identify that PRMT5 directly regulates the transcription of the epidermal growth factor receptor (EGFR). PRMT5, through independent modulation of the EGFR and AKT pathways, orchestrates the activation of NFκB, resulting in the upregulation of the pro-EMT transcription factors ZEB1, SNAIL, and TWIST1. Notably, EGFR and AKT form a compensatory feedback loop, reinforcing the expression of these EMT transcription factors. Small molecule inhibition of PRMT5 methyltransferase activity disrupts EGFR/AKT signaling, suppresses EMT transcription factor expression and ablates tumor growth in vivo . Our findings underscore the pivotal role of PRMT5 in the control of the EMT program in high-risk neuroblastoma.

TE-TSS: an integrated data resource of human and mouse transposable element (TE)-derived transcription start site (TSS).

Transposable elements (TEs) are abundant in the genome and serve as crucial regulatory elements. Some TEs function as epigenetically regulated promoters, and these TE-derived transcription start sites (TSSs) play a crucial role in regulating genes associated with specific functions, such as cancer and embryogenesis. However, the lack of an accessible database that systematically gathers TE-derived TSS data is a current research gap. To address this, we established TE-TSS, an integrated data resource of human and mouse TE-derived TSSs (http://xozhanglab.com/TETSS). TE-TSS has compiled 2681 RNA sequencing datasets, spanning various tissues, cell lines and developmental stages. From these, we identified 5768 human TE-derived TSSs and 2797 mouse TE-derived TSSs, with 47% and 38% being experimentally validated, respectively. TE-TSS enables comprehensive exploration of TSS usage in diverse samples, providing insights into tissue-specific gene expression patterns and transcriptional regulatory elements. Furthermore, TE-TSS compares TE-derived TSS regions across 15 mammalian species, enhancing our understanding of their evolutionary and functional aspects. The establishment of TE-TSS facilitates further investigations into the roles of TEs in shaping the transcriptomic landscape and offers valuable resources for comprehending their involvement in diverse biological processes.

Allele-specific binding (ASB) analyzer for annotation of allele-specific binding SNPs.

BACKGROUND: Allele-specific binding (ASB) events occur when transcription factors (TFs) bind more favorably to one of the two parental alleles at heterozygous single nucleotide polymorphisms (SNPs). Evidence suggests that ASB events could reveal the impact of sequence variations on TF binding and may have implications for the risk of diseases. RESULTS: Here we present ASB-analyzer, a software platform that enables the users to quickly and efficiently input raw sequencing data to generate individual reports containing the cytogenetic map of ASB SNPs and their associated phenotypes. This interactive tool thereby combines ASB SNP identification, biological annotation, motif analysis, phenotype associations and report summary in one pipeline. With this pipeline, we identified 3772 ASB SNPs from thirty GM12878 ChIP-seq datasets and demonstrated that the ASB SNPs were more likely to be enriched at important sites in TF-binding domains. CONCLUSIONS: ASB-analyzer is a user-friendly tool that enables the detection, characterization and visualization of ASB SNPs. It is implemented in Python, R and bash shell and packaged in the Conda environment. It is available as an open-source tool on GitHub at https://github.com/Liying1996/ASBanalyzer .

BCL2 is a major regulator of haploidy maintenance in murine embryonic stem cells.

Mammalian haploid cells are important resources for forward genetic screening and are important in genetic medicine and drug development. However, the self-diploidization of murine haploid embryonic stem cells (haESCs) during daily culture or differentiation jeopardizes their use in genetic approaches. Here, we show that overexpression (OE) of an antiapoptosis gene, BCL2, in haESCs robustly ensures their haploidy maintenance in various situations, even under strict differentiation in vivo (embryonic 10.5 chimeric fetus or 21-day teratoma). Haploid cell lines of many lineages, including epiblasts, trophectodermal lineages, and neuroectodermal lineages, can be easily derived by the differentiation of BCL2-OE haESCs in vitro. Transcriptome analysis revealed that BCL2-OE activates another regulatory gene, Has2, which is also sufficient for haploidy maintenance. Together, our findings provide an effective and secure strategy to reduce diploidization during differentiation, which will contribute to the generation of haploid cell lines of the desired lineage and related genetic screening.

Epigenetic and chromosomal features drive transposon insertion in Drosophila melanogaster.

Transposons are mobile genetic elements prevalent in the genomes of most species. The distribution of transposons within a genome reflects the actions of two opposing processes: initial insertion site selection, and selective pressure from the host. By analyzing whole-genome sequencing data from transposon-activated Drosophila melanogaster, we identified 43 316 de novo and 237 germline insertions from four long-terminal-repeat (LTR) transposons, one LINE transposon (I-element), and one DNA transposon (P-element). We found that all transposon types favored insertion into promoters de novo, but otherwise displayed distinct insertion patterns. De novo and germline P-element insertions preferred replication origins, often landing in a narrow region around transcription start sites and in regions of high chromatin accessibility. De novo LTR transposon insertions preferred regions with high H3K36me3, promoters and exons of active genes; within genes, LTR insertion frequency correlated with gene expression. De novo I-element insertion density increased with distance from the centromere. Germline I-element and LTR transposon insertions were depleted in promoters and exons, suggesting strong selective pressure to remove transposons from functional elements. Transposon movement is associated with genome evolution and disease; therefore, our results can improve our understanding of genome and disease biology.

Arih2 regulates Hedgehog signaling through smoothened ubiquitylation and ER-associated degradation.

During Hedgehog signaling, the ciliary levels of Ptch1 and Smo are regulated by the pathway. At the basal state, Ptch1 localizes to cilia and prevents the ciliary accumulation and activation of Smo. Upon binding a Hedgehog ligand, Ptch1 exits cilia, relieving inhibition of Smo. Smo then concentrates in cilia, becomes activated and activates downstream signaling. Loss of the ubiquitin E3 ligase Arih2 elevates basal Hedgehog signaling, elevates the cellular level of Smo and increases basal levels of ciliary Smo. Mice express two isoforms of Arih2 with Arih2α found primarily in the nucleus and Arih2ß found on the cytoplasmic face of the endoplasmic reticulum (ER). Re-expression of ER-localized Arih2ß but not nuclear-localized Arih2α rescues the Arih2 mutant phenotypes. When Arih2 is defective, protein aggregates accumulate in the ER and the unfolded protein response is activated. Arih2ß appears to regulate the ER-associated degradation (ERAD) of Smo preventing excess and potentially misfolded Smo from reaching the cilium and interfering with pathway regulation.

PRMT5 activates AKT via methylation to promote tumor metastasis.

Protein arginine methyltransferase 5 (PRMT5) is the primary methyltransferase generating symmetric-dimethyl-arginine marks on histone and non-histone proteins. PRMT5 dysregulation is implicated in multiple oncogenic processes. Here, we report that PRMT5-mediated methylation of protein kinase B (AKT) is required for its subsequent phosphorylation at Thr308 and Ser473. Moreover, pharmacologic or genetic inhibition of PRMT5 abolishes AKT1 arginine 15 methylation, thereby preventing AKT1 translocation to the plasma membrane and subsequent recruitment of its upstream activating kinases PDK1 and mTOR2. We show that PRMT5/AKT signaling controls the expression of the epithelial-mesenchymal-transition transcription factors ZEB1, SNAIL, and TWIST1. PRMT5 inhibition significantly attenuates primary tumor growth and broadly blocks metastasis in multiple organs in xenograft tumor models of high-risk neuroblastoma. Collectively, our results suggest that PRMT5 inhibition augments anti-AKT or other downstream targeted therapeutics in high-risk metastatic cancers.

Deletion and replacement of long genomic sequences using prime editing.

Genomic insertions, duplications and insertion/deletions (indels), which account for ~14% of human pathogenic mutations, cannot be accurately or efficiently corrected by current gene-editing methods, especially those that involve larger alterations (>100 base pairs (bp)). Here, we optimize prime editing (PE) tools for creating precise genomic deletions and direct the replacement of a genomic fragment ranging from ~1 kilobases (kb) to ~10 kb with a desired sequence (up to 60 bp) in the absence of an exogenous DNA template. By conjugating Cas9 nuclease to reverse transcriptase (PE-Cas9) and combining it with two PE guide RNAs (pegRNAs) targeting complementary DNA strands, we achieve precise and specific deletion and repair of target sequences via using this PE-Cas9-based deletion and repair (PEDAR) method. PEDAR outperformed other genome-editing methods in a reporter system and at endogenous loci, efficiently creating large and precise genomic alterations. In a mouse model of tyrosinemia, PEDAR removed a 1.38-kb pathogenic insertion within the Fah gene and precisely repaired the deletion junction to restore FAH expression in liver.

Integration of high-resolution promoter profiling assays reveals novel, cell type-specific transcription start sites across 115 human cell and tissue types.

Accurate transcription start site (TSS) annotations are essential for understanding transcriptional regulation and its role in human disease. Gene collections such as GENCODE contain annotations for tens of thousands of TSSs, but not all of these annotations are experimentally validated nor do they contain information on cell type-specific usage. Therefore, we sought to generate a collection of experimentally validated TSSs by integrating RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression (RAMPAGE) data from 115 cell and tissue types, which resulted in a collection of approximately 50 thousand representative RAMPAGE peaks. These peaks are primarily proximal to GENCODE-annotated TSSs and are concordant with other transcription assays. Because RAMPAGE uses paired-end reads, we were then able to connect peaks to transcripts by analyzing the genomic positions of the 3' ends of read mates. Using this paired-end information, we classified the vast majority (37 thousand) of our RAMPAGE peaks as verified TSSs, updating TSS annotations for 20% of GENCODE genes. We also found that these updated TSS annotations are supported by epigenomic and other transcriptomic data sets. To show the utility of this RAMPAGE rPeak collection, we intersected it with the NHGRI/EBI genome-wide association study (GWAS) catalog and identified new candidate GWAS genes. Overall, our work shows the importance of integrating experimental data to further refine TSS annotations and provides a valuable resource for the biological community.

5'-Modifications improve potency and efficacy of DNA donors for precision genome editing.

Nuclease-directed genome editing is a powerful tool for investigating physiology and has great promise as a therapeutic approach to correct mutations that cause disease. In its most precise form, genome editing can use cellular homology-directed repair (HDR) pathways to insert information from an exogenously supplied DNA-repair template (donor) directly into a targeted genomic location. Unfortunately, particularly for long insertions, toxicity and delivery considerations associated with repair template DNA can limit HDR efficacy. Here, we explore chemical modifications to both double-stranded and single-stranded DNA-repair templates. We describe 5'-terminal modifications, including in its simplest form the incorporation of triethylene glycol (TEG) moieties, that consistently increase the frequency of precision editing in the germlines of three animal models (Caenorhabditis elegans, zebrafish, mice) and in cultured human cells.

An organ-on-a-chip model for pre-clinical drug evaluation in progressive non-genetic cardiomyopathy.

Angiotensin II (Ang II) presents a critical mediator in various pathological conditions such as non-genetic cardiomyopathy. Osmotic pump infusion in rodents is a commonly used approach to model cardiomyopathy associated with Ang II. However, profound differences in electrophysiology and pharmacokinetics between rodent and human cardiomyocytes may limit predictability of animal-based experiments. This study investigates the application of an Organ-on-a-chip (OOC) system in modeling Ang II-induced progressive cardiomyopathy. The disease model is constructed to recapitulate myocardial response to Ang II in a temporal manner. The long-term tissue cultivation and non-invasive functional readouts enable monitoring of both acute and chronic cardiac responses to Ang II stimulation. Along with mapping of cytokine secretion and proteomic profiles, this model presents an opportunity to quantitatively measure the dynamic pathological changes that could not be otherwise identified in animals. Further, we present this model as a testbed to evaluate compounds that target Ang II-induced cardiac remodeling. Through assessing the effects of losartan, relaxin, and saracatinib, the drug screening data implicated multifaceted cardioprotective effects of relaxin in restoring contractile function and reducing fibrotic remodeling. Overall, this study provides a controllable platform where cardiac activities can be explicitly observed and tested over the pathological process. The facile and high-content screening can facilitate the evaluation of potential drug candidates in the pre-clinical stage.

Genetic and epigenetic features of promoters with ubiquitous chromatin accessibility support ubiquitous transcription of cell-essential genes.

Gene expression is controlled by regulatory elements within accessible chromatin. Although most regulatory elements are cell type-specific, a subset is accessible in nearly all the 517 human and 94 mouse cell and tissue types assayed by the ENCODE consortium. We systematically analyzed 9000 human and 8000 mouse ubiquitously-accessible candidate cis-regulatory elements (cCREs) with promoter-like signatures (PLSs) from ENCODE, which we denote ubi-PLSs. These are more CpG-rich than non-ubi-PLSs and correspond to genes with ubiquitously high transcription, including a majority of cell-essential genes. ubi-PLSs are enriched with motifs of ubiquitously-expressed transcription factors and preferentially bound by transcriptional cofactors regulating ubiquitously-expressed genes. They are highly conserved between human and mouse at the synteny level but exhibit frequent turnover of motif sites; accordingly, ubi-PLSs show increased variation at their centers compared with flanking regions among the ∼186 thousand human genomes sequenced by the TOPMed project. Finally, ubi-PLSs are enriched in genes implicated in Mendelian diseases, especially diseases broadly impacting most cell types, such as deficiencies in mitochondrial functions. Thus, a set of roughly 9000 mammalian promoters are actively maintained in an accessible state across cell types by a distinct set of transcription factors and cofactors to ensure the transcriptional programs of cell-essential genes.

Investigating the Potential Roles of SINEs in the Human Genome.

Short interspersed nuclear elements (SINEs) are nonautonomous retrotransposons that occupy approximately 13% of the human genome. They are transcribed by RNA polymerase III and can be retrotranscribed and inserted back into the genome with the help of other autonomous retroelements. Because they are preferentially located close to or within gene-rich regions, they can regulate gene expression by various mechanisms that act at both the DNA and the RNA levels. In this review, we summarize recent findings on the involvement of SINEs in different types of gene regulation and discuss the potential regulatory functions of SINEs that are in close proximity to genes, Pol III-transcribed SINE RNAs, and embedded SINE sequences within Pol II-transcribed genes in the human genome. These discoveries illustrate how the human genome has exapted some SINEs into functional regulatory elements.

Optimized RNA-targeting CRISPR/Cas13d technology outperforms shRNA in identifying functional circRNAs.

Short hairpin RNAs (shRNAs) are used to deplete circRNAs by targeting back-splicing junction (BSJ) sites. However, frequent discrepancies exist between shRNA-mediated circRNA knockdown and the corresponding biological effect, querying their robustness. By leveraging CRISPR/Cas13d tool and optimizing the strategy for designing single-guide RNAs against circRNA BSJ sites, we markedly enhance specificity of circRNA silencing. This specificity is validated in parallel screenings by shRNA and CRISPR/Cas13d libraries. Using a CRISPR/Cas13d screening library targeting > 2500 human hepatocellular carcinoma-related circRNAs, we subsequently identify a subset of sorafenib-resistant circRNAs. Thus, CRISPR/Cas13d represents an effective approach for high-throughput study of functional circRNAs.

YAP1 Withdrawal in Hepatoblastoma Drives Therapeutic Differentiation of Tumor Cells to Functional Hepatocyte-Like Cells.

BACKGROUND AND AIMS: Despite surgical and chemotherapeutic advances, the 5-year survival rate for stage IV hepatoblastoma (HB), the predominant pediatric liver tumor, remains at 27%. Yes-associated protein 1 (YAP1) and ß-catenin co-activation occurs in 80% of children's HB; however, a lack of conditional genetic models precludes tumor maintenance exploration. Thus, the need for a targeted therapy remains unmet. Given the predominance of YAP1 and ß-catenin activation in HB, we sought to evaluate YAP1 as a therapeutic target in HB. APPROACH AND RESULTS: We engineered the conditional HB murine model using hydrodynamic injection to deliver transposon plasmids encoding inducible YAP1S127A , constitutive ß-cateninDelN90 , and a luciferase reporter to murine liver. Tumor regression was evaluated using bioluminescent imaging, tumor landscape characterized using RNA and ATAC sequencing, and DNA footprinting. Here we show that YAP1S127A withdrawal mediates more than 90% tumor regression with survival for 230+ days in mice. YAP1S127A withdrawal promotes apoptosis in a subset of tumor cells, and in remaining cells induces a cell fate switch that drives therapeutic differentiation of HB tumors into Ki-67-negative hepatocyte-like HB cells ("HbHeps") with hepatocyte-like morphology and mature hepatocyte gene expression. YAP1S127A withdrawal drives the formation of hbHeps by modulating liver differentiation transcription factor occupancy. Indeed, tumor-derived hbHeps, consistent with their reprogrammed transcriptional landscape, regain partial hepatocyte function and rescue liver damage in mice. CONCLUSIONS: YAP1S127A withdrawal, without silencing oncogenic ß-catenin, significantly regresses hepatoblastoma, providing in vivo data to support YAP1 as a therapeutic target for HB. YAP1S127A withdrawal alone sufficiently drives long-term regression in HB, as it promotes cell death in a subset of tumor cells and modulates transcription factor occupancy to reverse the fate of residual tumor cells to mimic functional hepatocytes.

Expanded encyclopaedias of DNA elements in the human and mouse genomes.

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

Chemical modifications of adenine base editor mRNA and guide RNA expand its application scope.

CRISPR-Cas9-associated base editing is a promising tool to correct pathogenic single nucleotide mutations in research or therapeutic settings. Efficient base editing requires cellular exposure to levels of base editors that can be difficult to attain in hard-to-transfect cells or in vivo. Here we engineer a chemically modified mRNA-encoded adenine base editor that mediates robust editing at various cellular genomic sites together with moderately modified guide RNA, and show its therapeutic potential in correcting pathogenic single nucleotide mutations in cell and animal models of diseases. The optimized chemical modifications of adenine base editor mRNA and guide RNA expand the applicability of CRISPR-associated gene editing tools in vitro and in vivo.

Comprehensive identification of alternative back-splicing in human tissue transcriptomes.

Circular RNAs (circRNAs) are covalently closed RNAs derived from back-splicing of genes across eukaryotes. Through alternative back-splicing (ABS), a single gene produces multiple circRNAs sharing the same back-splice site. Although many ABS events have recently been discovered, to what extent ABS involves in circRNA biogenesis and how it is regulated in different human tissues still remain elusive. Here, we reported an in-depth analysis of ABS events in 90 human tissue transcriptomes. We observed that ABS occurred for about 84% circRNAs. Interestingly, alternative 5' back-splicing occurs more prevalently than alternative 3' back-splicing, and both of them are tissue-specific, especially enriched in brain tissues. In addition, the patterns of ABS events in different brain regions are similar to each other and are more complex than the patterns in non-brain tissues. Finally, the intron length and abundance of Alu elements positively correlated with ABS event complexity, and the predominant circRNAs had longer flanking introns and more Alu elements than other circRNAs in the same ABS event. Together, our results represent a resource for circRNA research-we expanded the repertoire of ABS events of circRNAs in human tissue transcriptomes and provided insights into the complexity of circRNA biogenesis, expression, and regulation.