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

LINC01116-dependent upregulation of RNA polymerase I transcription drives oncogenic phenotypes in lung adenocarcinoma

BACKGROUND

Hyperactive RNA Polymerase I (Pol I) transcription is canonical in cancer, associated with malignant proliferation, poor prognosis, epithelial-mesenchymal transition, and chemotherapy resistance. Despite its significance, the molecular mechanisms underlying Pol I hyperactivity remain unclear. This study aims to elucidate the role of long noncoding RNAs (lncRNAs) in regulating Pol I transcription in lung adenocarcinoma (LUAD).

METHODS

Bioinformatics analyses were applied to identify lncRNAs interacting with Pol I transcriptional machinery. Fluorescence in situ hybridization was employed to examine the nucleolar localization of candidate lncRNA in LUAD cells. RNA immunoprecipitation assay validated the interaction between candidate lncRNA and Pol I components. Chromatin isolation by RNA purification and Chromatin Immunoprecipitation (ChIP) were utilized to confirm the interactions of candidate lncRNA with Pol I transcriptional machinery and the rDNA core promoter. Functional analyses, including lncRNA knock-in and knockdown, inhibition of Pol I transcription, quantitative PCR, cell proliferation, clonogenicity, apoptosis, cell cycle, wound-healing, and invasion assays, were performed to determine the effect of candidate lncRNA on Pol I transcription and associated malignant phenotypes in LUAD cells. ChIP assays and luminometry were used to investigate the transcriptional regulation of the candidate lncRNA.

RESULTS

We demonstrate that oncogenic LINC01116 scaffolds essential Pol I transcription factors TAF1A and TAF1D, to the ribosomal DNA promoter, and upregulate Pol I transcription. Crucially, LINC01116-driven Pol I transcription activation is essential for its oncogenic activities. Inhibition of Pol I transcription abrogated LINC01116-induced oncogenic phenotypes, including increased proliferation, cell cycle progression, clonogenicity, reduced apoptosis, increased migration and invasion, and drug sensitivity. Conversely, LINC01116 knockdown reversed these effects. Additionally, we show that LINC01116 upregulation in LUAD is driven by the oncogene c-Myc, a known Pol I transcription activator, indicating a functional regulatory feedback loop within the c-Myc-LINC01116-Pol I transcription axis.

CONCLUSION

Collectively, our findings reveal, for the first time, that LINC01116 enhances Pol I transcription by scaffolding essential transcription factors to the ribosomal DNA promoter, thereby driving oncogenic activities in LUAD. We propose the c-Myc-LINC01116-Pol I axis as a critical oncogenic pathway and a potential therapeutic target for modulating Pol I transcription in LUAD.

The long non-coding RNA ROSALIND protects the mitochondrial translational machinery from oxidative damage

Upregulation of mitochondrial respiration coupled with high ROS-scavenging capacity is a characteristic shared by drug-tolerant cells in several cancers. As translational fidelity is essential for cell fitness, protection of the mitochondrial and cytosolic ribosomes from oxidative damage is pivotal. While mechanisms for recognising and repairing such damage exist in the cytoplasm, the corresponding process in the mitochondria remains unclear.By performing Ascorbate PEroXidase (APEX)-proximity ligation assay directed to the mitochondrial matrix followed by isolation and sequencing of RNA associated to mitochondrial proteins, we identified the nuclear-encoded lncRNA ROSALIND as an interacting partner of ribosomes. ROSALIND is upregulated in recurrent tumours and its expression can discriminate between responders and non-responders to immune checkpoint blockade in a melanoma cohort of patients. Featuring an unusually high G content, ROSALIND serves as a substrate for oxidation. Consequently, inhibiting ROSALIND leads to an increase in ROS and protein oxidation, resulting in severe mitochondrial respiration defects. This, in turn, impairs melanoma cell viability and increases immunogenicity both in vitro and ex vivo in preclinical humanised cancer models. These findings underscore the role of ROSALIND as a novel ROS buffering system, safeguarding mitochondrial translation from oxidative stress, and shed light on potential therapeutic strategies for overcoming cancer therapy resistance.

HAND2-AS1 plays a tumor-suppressive role in hepatoblastoma through the negative regulation of CDK1

Objective

Hepatoblastoma (HB) is the most commonly seen pediatric liver malignancy. The preliminary experiment of our research group found that cyclin dependent kinase 1 (CDK1) was upregulated in HB. By in silico analysis, long noncoding RNA (lncRNA) HAND2 antisense RNA 1 (HAND2-AS1) was determined as the research object. Herein, HAND2-AS1 expression in HB and its effect and mechanism on HB were extensively investigated.

Methods

CDK1-related lncRNAs were searched using the microarray data from the Gene Expression Omnibus (GEO) database and Gene Expression Profiling Interactive Analysis (GEPIA) online database. qRT-PCR, Western blot, and immunohistochemistry were performed to determine the mRNA expression and protein levels of target genes. MTT, flow cytometry and DAPI staining assays were conducted to measure proliferation activity, cell cycle progression, and apoptosis of HB cells. The interaction between lncRNA and protein was determined by RNA pull-down and FISH assays. Luciferase assay was applied to identify whether HAND2-AS1 stimulates the transcription of CDK1. CDK1 mRNA stability was detected through actinomycin D assay. Aycloheximide assay was used to detect the CDK1 protein stability.

Results

HAND2-AS1 was downregulated in HB tissues and cells. HAND2-AS1 overexpression impeded HB cells proliferation activity and cycle progression while inducing cell apoptosis of HB cells, while knockdown of HAND2-AS1 emerged the opposite effect. HAND2-AS1 negatively correlated with CDK1. HAND2-AS1 downregulated CDK1 expression by affecting the transcriptional activity, mRNA and protein stability of CDK1. Furthermore, HAND2-AS1 impeded HB cell proliferation and cycle progression while inducing cell apoptosis by downregulating CDK1.

Conclusion

Our research highlights that HAND2-AS1 can exert a tumor-suppressive effect on HB through the negative regulation of CDK1, and the HAND2-AS1/CDK1 is expected to be a diagnostic molecular marker and therapeutic target for HB in clinical practice.

Hierarchical lncRNA regulatory network in early-onset severe preeclampsia

BACKGROUND

Recent studies have shown that several long non-coding RNAs (lncRNAs) in the placenta are associated with preeclampsia (PE). However, the extent to which lncRNAs may contribute to the pathological progression of PE is unclear.

RESULTS

Here, we report a hierarchical regulatory network involved in early-onset severe PE (EOSPE). We have carried out transcriptome sequencing on the placentae from patients and normal subjects to identify the differentially expressed genes (DEGs), including some lncRNAs (DElncRNAs). We then constructed a high-quality hierarchical regulatory network of lncRNAs, transcription factors (TFs), and target DEGs, containing 1851 lncRNA-TF interactions and 6901 TF-promoter interactions. The lncRNA-to-target regulatory interactions were further validated by the triplex structures between the DElncRNAs and the promoters of the target DEGs. The DElncRNAs in the regulatory network were clustered into 3 clusters, one containing DElncRNAs correlated with the blood pressure, including FLNB-AS1 with targeting 27.89% (869/3116) DEGs in EOSPE. We further demonstrated that FLNB-AS1 could bind the transcription factor JUNB to regulate a series members of the HIF-1 signaling pathway in trophoblast cells.

CONCLUSIONS

Our results suggest that the differential expression of lncRNAs may perturb the lncRNA-TF-DEG hierarchical regulatory network, leading to the dysregulation of many genes involved in EOSPE. Our study provides a new strategy and a valuable resource for studying the mechanism underlying gene dysregulation in EOSPE patients.

Structure-based mechanism of riboregulation of the metabolic enzyme SHMT1

RNA can directly control protein activity in a process called riboregulation; only a few mechanisms of riboregulation have been described in detail, none of which have been characterized on structural grounds. Here, we present a comprehensive structural, functional, and phylogenetic analysis of riboregulation of cytosolic serine hydroxymethyltransferase (SHMT1), the enzyme interconverting serine and glycine in one-carbon metabolism. We have determined the cryoelectron microscopy (cryo-EM) structure of human SHMT1 in its free- and RNA-bound states, and we show that the RNA modulator competes with polyglutamylated folates and acts as an allosteric switch, selectively altering the enzyme's reactivity vs. serine. In addition, we identify the tetrameric assembly and a flap structural motif as key structural elements necessary for binding of RNA to eukaryotic SHMT1. The results presented here suggest that riboregulation may have played a role in evolution of eukaryotic SHMT1 and in compartmentalization of one-carbon metabolism. Our findings provide insights for RNA-based therapeutic strategies targeting this cancer-linked metabolic pathway.

Pyroptosis is involved in the immune microenvironment regulation of unexplained recurrent miscarriage

Unexplained recurrent miscarriage (URM) is a common pregnancy complication with few effective therapies. Moreover, little is known regarding the role of pyroptosis in the regulation of the URM immune microenvironment. To address this issue, gene expression profiles of publicly available placental datasets GSE22490 and GSE76862 were downloaded from the Gene Expression Omnibus database. Pyroptosis-related differentially expressed genes were identified and a total of 16 differentially expressed genes associated with pyroptosis were detected, among which 1 was upregulated and 15 were downregulated. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses indicated that the functionally enriched modules and pathways of these genes are closely related to immune and inflammatory responses. Four hub genes were identified: BTK, TLR8, NLRC4, and TNFSF13B. BTK, TLR8, and TNFSF13B were highly connected with immune cells, according to the correlation analysis of four hub genes and 20 different types of immune cells (p < 0.05). The four hub genes were used as research objects to construct the interaction networks. Chorionic villus tissue was used for quantitative real-time polymerase chain reaction and western blot to confirm the expression levels of hub genes, and the results showed that the expression of the four hub genes was significantly decreased in the chorionic villus tissue in the URM group. Collectively, the present study indicates that perhaps pyroptosis is essential to the diversity and complexity of the URM immune microenvironment, and provides a theoretical basis and research ideas for subsequent target gene verification and mechanism research.

Spatiotemporal constrained RNA-protein heterogeneous network for protein complex identification

The identification of protein complexes from protein interaction networks is crucial in the understanding of protein function, cellular processes and disease mechanisms. Existing methods commonly rely on the assumption that protein interaction networks are highly reliable, yet in reality, there is considerable noise in the data. In addition, these methods fail to account for the regulatory roles of biomolecules during the formation of protein complexes, which is crucial for understanding the generation of protein interactions. To this end, we propose a SpatioTemporal constrained RNA-protein heterogeneous network for Protein Complex Identification (STRPCI). STRPCI first constructs a multiplex heterogeneous protein information network to capture deep semantic information by extracting spatiotemporal interaction patterns. Then, it utilizes a dual-view aggregator to aggregate heterogeneous neighbor information from different layers. Finally, through contrastive learning, STRPCI collaboratively optimizes the protein embedding representations under different spatiotemporal interaction patterns. Based on the protein embedding similarity, STRPCI reweights the protein interaction network and identifies protein complexes with core-attachment strategy. By considering the spatiotemporal constraints and biomolecular regulatory factors of protein interactions, STRPCI measures the tightness of interactions, thus mitigating the impact of noisy data on complex identification. Evaluation results on four real PPI networks demonstrate the effectiveness and strong biological significance of STRPCI. The source code implementation of STRPCI is available from https://github.com/LI-jasm/STRPCI.

U3 snoRNA inter-regulates with DDX21 in the perichromosomal region to control mitosis

U3 snoRNA is essential for ribosome biogenesis during interphase. Upon mitotic onset, the nucleolus disassembles and U3 snoRNA relocates to the perichromosomal region (PR) to be considered as a chromosome passenger. Whether U3 controls mitosis remains unknown. Here, we demonstrate that U3 snoRNA is required for mitotic progression. We identified DDX21 as the predominant U3-binding protein during mitosis and confirmed that U3 snoRNA colocalizes with DDX21 in the PR. DDX21 knockdown induces mitotic catastrophe and similar mitotic defects caused by U3 snoRNA depletion. Interestingly, the uniform PR distribution of U3 snoRNA and DDX21 is interdependent. DDX21 functions in mitosis depending on its PR localization. Mechanistically, U3 snoRNA regulates DDX21 PR localization through maintaining its mobility. Moreover, Cy5-U3 snoRNA downsizes the fibrous condensates of His-DDX21 at proper molecular ratios in vitro. This work highlights the importance of the equilibrium between U3 snoRNA and DDX21 in PR formation and reveals the potential relationship between the PR assembly and mitotic regulation.

Prediction of protein-RNA interactions from single-cell transcriptomic data

Proteins are crucial in regulating every aspect of RNA life, yet understanding their interactions with coding and noncoding RNAs remains limited. Experimental studies are typically restricted to a small number of cell lines and a limited set of RNA-binding proteins (RBPs). Although computational methods based on physico-chemical principles can predict protein-RNA interactions accurately, they often lack the ability to consider cell-type-specific gene expression and the broader context of gene regulatory networks (GRNs). Here, we assess the performance of several GRN inference algorithms in predicting protein-RNA interactions from single-cell transcriptomic data, and propose a pipeline, called scRAPID (single-cell transcriptomic-based RnA Protein Interaction Detection), that integrates these methods with the catRAPID algorithm, which can identify direct physical interactions between RBPs and RNA molecules. Our approach demonstrates that RBP-RNA interactions can be predicted from single-cell transcriptomic data, with performances comparable or superior to those achieved for the well-established task of inferring transcription factor-target interactions. The incorporation of catRAPID significantly enhances the accuracy of identifying interactions, particularly with long noncoding RNAs, and enables the identification of hub RBPs and RNAs. Additionally, we show that interactions between RBPs can be detected based on their inferred RNA targets. The software is freely available at https://github.com/tartaglialabIIT/scRAPID.

Predicting nuclear G-quadruplex RNA-binding proteins with roles in transcription and phase separation

RNA-binding proteins are central for many biological processes and their characterization has demonstrated a broad range of functions as well as a wide spectrum of target structures. RNA G-quadruplexes are important regulatory elements occurring in both coding and non-coding transcripts, yet our knowledge of their structure-based interactions is at present limited. Here, using theoretical predictions and experimental approaches, we show that many chromatin-binding proteins bind to RNA G-quadruplexes, and we classify them based on their RNA G-quadruplex-binding potential. Combining experimental identification of nuclear RNA G-quadruplex-binding proteins with computational approaches, we build a prediction tool that assigns probability score for a nuclear protein to bind RNA G-quadruplexes. We show that predicted G-quadruplex RNA-binding proteins exhibit a high degree of protein disorder and hydrophilicity and suggest involvement in both transcription and phase-separation into membrane-less organelles. Finally, we present the G4-Folded/UNfolded Nuclear Interaction Explorer System (G4-FUNNIES) for estimating RNA G4-binding propensities at http://service.tartaglialab.com/new_submission/G4FUNNIES .

Endothelial downregulation of nuclear m6A reader YTHDC1 promotes pulmonary vascular remodeling in sugen hypoxia model of pulmonary hypertension

Background

Pulmonary hypertension (PH) is characterized with vascular remodeling, which is intiated by vascular endothelial dysfunction. N6-methyladenosine (m6A) modification mediates gene expression in many ways including mediating RNA degradation, splicing, nuclear export et al. m6A modification have been found to be associated with the development of PH. However, the role of m6A regulators in pulmonary artery endothelial cells (PAECs) dysfunction of PH is still under research.

Methods

The expression levels of m6A regulators in PAECs were analyzed with the single-cell sequencing Data(scRNA). Next, the target differentially expressed genes (DEGs) of m6A regulators in PAECs were functionally annotated. The analysis of cellular interactions included the examination of receptor-ligand pairs regulated by m6A regulators. Pseudo-time trajectory analyses and a ceRNA network involving lncRNAs, miRNAs, and mRNAs were conducted in PAECs. Furthermore, microarray data (GSE180169) for Sugen Hypoxia PH (SuHx PH) mouse models was screened for DEGs and m6A regulators in PAECs. Moreover, the expression of YTHDC1 in the lung samples of SuHx PH models was determined using immunofluorescence. In vitro, the mRNA expression of YTHDC1 in HPAECs under hypoxia conditions was detected. The effect of YTHDC1 recombinant protein on HPAEC proliferation was detected by Cell Counting Kit-8 (CCK8).

Results

Dysregulation of m6A regulators was observed in mouse PAECs. The m6A reader of YTHDC1 was decreased in PAECs in scRNA data and RNAseq data of isolated PAECs of SuHx PH models. Downregulation of YTHDC1 was caused by hypoxia in PAECs in vitro and similar results was observed in PAECs of SuHx PH mouse models. Next, YTHDC1 recombinant protein was found to inhibit HPAECs proliferation. The DEGs targeted by YTHDC1 were enriched in angiogenesis, endothelial cell migration, fluid shear stress, and stem cell maintenance. Analysis indicates that interactions among endothelial cells, smooth muscle cells, fibroblasts, and immune cells, mediated by specific YTHDC1 target genes (e.g., PTPRC-MRC1, ITBG2-ICAM1, COL4A1-CD44), contribute to PH development. Also, the YTHDC1 expression were consistent with Thioredoxin interacting protein (TXNIP). What's more, the predicted transcription factors showed that NFKB1, Foxd3 may be involved in the regulation of YTHDC1. Lastly, our data suggest that YTHDC1 may be involved in regulating PAECs dysfunction through lncRNA/miRNA/mRNA network.

Conclusion

For the first time, we analyzed changes in the expression and biological functions of m6A regulators in SuHx PH mouse models. We causatively linked YTHDC1 to PAECs dysfunction, providing novel insight into and opportunities to diagnose and treat PH.

G3BP1 Interact with JAK2 mRNA to Promote the Malignant Progression of Nasopharyngeal Carcinoma via Activating JAK2/STAT3 Signaling Pathway

Ras-GTPase-activating protein (GAP)-binding protein 1 (G3BP1) is an RNA-binding protein implicated in various malignancies. However, its role in nasopharyngeal carcinoma (NPC) remains elusive. This study elucidates the potential regulation mechanisms of G3BP1 and its significance in NPC advancement. Through knockdown and overexpression approaches, we validate G3BP1's oncogenic role by promoting proliferation, migration, and invasion in vitro and in vivo. Moreover, G3BP1 emerges as a key regulator of the JAK2/STAT3 signaling pathway, augmenting JAK2 expression via mRNA binding. Notably, epigallocatechin gallate (EGCG), a green tea-derived antioxidant, counteracts G3BP1-mediated pathway activation. Clinical analysis reveals heightened G3BP1, JAK2, and p-STAT3 as powerful prognostic markers, with G3BP1's expression standing as an independent indicator of poorer outcomes for NPC patients. In conclusion, the study unveils the oncogenic prowess of G3BP1, its orchestration of the JAK2/STAT3 signaling pathway, and its pivotal role in NPC progression.

PABPN1 functions as a predictive biomarker in colorectal carcinoma

PURPOSE

PABPN1 acts as a modulator of poly(A) tail length and alternative polyadenylation. This research was aimed to explore the role of PABPN1 in colorectal cancer (CRC).

METHODS

Public databases were performed to analyze expression, location, roles of prognosis and tumor immunity and interaction with RNAs and proteins of PABPN1. To investigate PABPN1 expression in tissues, 78 CRC specimens were collected to conduct IHC, and 30 pairs of frozen CRC and corresponding adjacent normal tissues were used to conduct qRT-PCR and WB. In addition, in vitro experiments were then carried out to identify the role of PABPN1 in CRC.

RESULTS

Compared with normal tissues, PABPN1 expression was significant higher in CRC. Its high level predicted poor outcome of CRC. Th1 and Treg had significant negative relationships not only with PABPN1 expression, but also with six molecules interacting with PABPN1, including IFT172, KIAA0895L, RECQL4, WDR6, PABPC1 and NCBP1. In addition, PABPN1 had negative relationships with quite a few immune markers, such as CSF1R, IL-10, CCL2 and so on. In cellular experiments, silencing PABPN1 inhibited proliferation and promoted apoptosis in HCT-116 CRC cells.

CONCLUSION

In summary, PABPN1 might become a novel biomarker and correlate with tumor immunity in CRC.

Unveiling the role of PUS7-mediated pseudouridylation in host protein interactions specific for the SARS-CoV-2 RNA genome

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a positive single-stranded RNA virus, engages in complex interactions with host cell proteins throughout its life cycle. While these interactions enable the host to recognize and inhibit viral replication, they also facilitate essential viral processes such as transcription, translation, and replication. Many aspects of these virus-host interactions remain poorly understood. Here, we employed the catRAPID algorithm and utilized the RNA-protein interaction detection coupled with mass spectrometry technology to predict and validate the host proteins that specifically bind to the highly structured 5' and 3' terminal regions of the SARS-CoV-2 RNA. Among the interactions identified, we prioritized pseudouridine synthase PUS7, which binds to both ends of the viral RNA. Using nanopore direct RNA sequencing, we discovered that the viral RNA undergoes extensive post-transcriptional modifications. Modified consensus regions for PUS7 were identified at both terminal regions of the SARS-CoV-2 RNA, including one in the viral transcription regulatory sequence leader. Collectively, our findings offer insights into host protein interactions with the SARS-CoV-2 UTRs and highlight the likely significance of pseudouridine synthases and other post-transcriptional modifications in the viral life cycle. This new knowledge enhances our understanding of virus-host dynamics and could inform the development of targeted therapeutic strategies.

RNA sequestration driven by amyloid formation: the alpha synuclein case

Nucleic acids can act as potent modulators of protein aggregation, and RNA has the ability to either hinder or facilitate protein assembly, depending on the molecular context. In this study, we utilized a computational approach to characterize the physico-chemical properties of regions involved in amyloid aggregation. In various experimental datasets, we observed that while the core is hydrophobic and highly ordered, external regions, which are more disordered, display a distinct tendency to interact with nucleic acids. To validate our predictions, we performed aggregation assays with alpha-synuclein (aS140), a non-nucleic acid-binding amyloidogenic protein, and a mutant truncated at the acidic C-terminus (aS103), which is predicted to have a higher tendency to interact with RNA. For both aS140 and aS103, we observed an acceleration of aggregation upon RNA addition, with a significantly stronger effect for aS103. Due to favorable electrostatics, we noted an enhanced nucleic acid sequestration ability for the aggregated aS103, allowing it to entrap a larger amount of RNA compared to the aggregated wild-type counterpart. Overall, our research suggests that RNA sequestration might be a common phenomenon linked to protein aggregation, constituting a gain-of-function mechanism that warrants further investigation.

Nucleocytoplasmic mRNA redistribution accompanies RNA binding protein mislocalization in ALS motor neurons and is restored by VCP ATPase inhibition

Amyotrophic lateral sclerosis (ALS) is characterized by nucleocytoplasmic mislocalization of the RNA-binding protein (RBP) TDP-43. However, emerging evidence suggests more widespread mRNA and protein mislocalization. Here, we employed nucleocytoplasmic fractionation, RNA sequencing, and mass spectrometry to investigate the localization of mRNA and protein in induced pluripotent stem cell-derived motor neurons (iPSMNs) from ALS patients with TARDBP and VCP mutations. ALS mutant iPSMNs exhibited extensive nucleocytoplasmic mRNA redistribution, RBP mislocalization, and splicing alterations. Mislocalized proteins exhibited a greater affinity for redistributed transcripts, suggesting a link between RBP mislocalization and mRNA redistribution. Notably, treatment with ML240, a VCP ATPase inhibitor, partially restored mRNA and protein localization in ALS mutant iPSMNs. ML240 induced changes in the VCP interactome and lysosomal localization and reduced oxidative stress and DNA damage. These findings emphasize the link between RBP mislocalization and mRNA redistribution in ALS motor neurons and highlight the therapeutic potential of VCP inhibition.

Transactive response DNA-binding protein 43 is enriched at the centrosome in human cells

The centrosome, as the main microtubule organizing centre, plays key roles in cell polarity, genome stability and ciliogenesis. The recent identification of ribosomes, RNA-binding proteins and transcripts at the centrosome suggests local protein synthesis. In this context, we hypothesized that TDP-43, a highly conserved RNA binding protein involved in the pathophysiology of amyotrophic lateral sclerosis and frontotemporal lobar degeneration, could be enriched at this organelle. Using dedicated high magnification sub-diffraction microscopy on human cells, we discovered a novel localization of TDP-43 at the centrosome during all phases of the cell cycle. These results were confirmed on purified centrosomes by western blot and immunofluorescence microscopy. In addition, the co-localization of TDP-43 and pericentrin suggested a pericentriolar enrichment of the protein, leading us to hypothesize that TDP-43 might interact with local mRNAs and proteins. Supporting this hypothesis, we found four conserved centrosomal mRNAs and 16 centrosomal proteins identified as direct TDP-43 interactors. More strikingly, all the 16 proteins are implicated in the pathophysiology of TDP-43 proteinopathies, suggesting that TDP-43 dysfunction in this organelle contributes to neurodegeneration. This first description of TDP-43 centrosomal enrichment paves the way for a more comprehensive understanding of TDP-43 physiology and pathology.

Role of LncRNA MIR99AHG in breast cancer: Bioinformatic analysis and preliminary verification

Objective

This research was aimed to preliminarily explore the clinical roles and potential molecular mechanisms of MIR99AHG and its significant transcripts in breast cancer (BRCA).

Methods

Public databases were utilized to analyze the expression and prognostic roles of MIR99AHG and its transcripts. Relationships between MIR99AHG expression and immune cells infiltration were analyzed in Xiantao platform. In addition, co-expressed genes and interacting proteins of MIR99AHG were predicted. CancerSEA analyzed its relationship with functional states. Next, CNV status, DNA methylation, interacting transcription factors (TFs) and ceRNA network were analyzed to explore its possible mechanisms. Then, RNA ISH and FISH assays were used to detect its expression and location in BRCA tissues and cell lines, respectively. Finally, qRT-PCR was utilized to investigate MIR99AHG expression in cell lines.

Results

Compared with the corresponding normal tissues, MIR99AHG expression levels were lower in all BRCA subtypes, and luminal B's was the lowest one. And MIR99AHG expression was negatively related to the tumor stage. In addition, 4 transcripts (ENST00000619222.4, ENST00000418813.6, ENST00000602901.5 and ENST00000453910.5) of MIR99AHG showed significant differences in the expression. Databases also suggested that the high MIR99AHG expression levels indicated good prognosis, especially in patients without lymph node metastasis. Xiantao found that MIR99AHG was positively related to 17 immune cells and negatively linked with 2 immune cells. CancerSEA analysis showed no relationships between MIR99AHG and functional states. From GEPIA and BCIP databases, 19 co-expressed genes were highly related to these four significant transcripts of MIR99AHG. StarBase, RNAct and HDOCK showed that several tumor-associated proteins, including U2AF65, hnRNPC, AEBP2, CHIC1 and so on, might interact with MIR99AHG. Genetically, BRCA had a higher proportion of MIR99AHG CNV loss than CNV gain, and the high level of DNA methylation indicated a good prognosis. Furthermore, 19 TFs were predicted to combine with the promoter of MIR99AHG. Then, we screened out 10 miRNAs potentially interacting with the significant transcripts of MIR99AHG, and five were significantly increased in breast tumors compared to normal tissues, including miR-194-5p, miR-320 b and so on, which could combine 14 mRNAs. Through ISH and FISH assays, we verified that MIR99AHG was down-regulated in BRCA samples and cell lines in comparison to non-tumor tissues and mammary epithelial cell line (MCF10A), and MIR99AHG was located both in cytoplasm and nucleus. qRT-PCR assay also showed the lower expression of MIR99AHG in breast cancer cells than that in MCF10A.

Conclusion

These results indicate that MIR99AHG can be a favorable prognostic indicator for BRCA. ENST00000619222.4, ENST00000418813.6, ENST00000602901.5 and ENST00000453910.5 are significant transcripts and their down-regulation may play crucial roles in the progression of BRCA.

COVID-19 mRNA vaccines as hypothetical epigenetic players: Results from an in silico analysis, considerations and perspectives

OBJECTIVES

To investigate in silico the occurrence of epigenetic crosstalk by nucleotide sequence complementarity between the BNT162b2 mRNA vaccine and whole human genome, including coding and noncoding (nc)RNA genes. To correlate these results with those obtained with the original spike (S) gene of Severe Acute Respiratory Syndrome CoronaVirus-2 (SARS-CoV-2).

METHODS

The publicly available FASTA sequence of the BNT162b2 mRNA vaccine and the SARS-CoV-2 isolate Wuhan-Hu-1 S gene (NC_045512.2) were used separately as key input to the Ensembl.org library to evaluate base pair match to human GRCh38 genome. Human coding and noncoding genes harboring hits were assessed for functional activity and health effects using bioinformatics tools and GWAS databases.

RESULTS

The BLAT analysis against the human GRCh38 genome revealed a total of 37 hits for BNT162b2 mRNA and no hits for the SARS-CoV-2 S gene. More specifically, BNT162b2 mRNA matched 19 human genes whose protein products are variously involved in enzyme reactions, nucleotide or cation binding, signaling, and carrier functions. In BLASTN analysis of ncRNA genes, BNT162b2 mRNA and SARS-CoV-2 S gene matched 17 and 24 different human genomic regions, respectively. Overall, characterization of the matched noncoding sequences revealed stronger interference with epigenetic pathways for BNT162b2 mRNA compared with the original S gene.

CONCLUSION

This pivotal in silico analysis shows that SARS-CoV-2 S gene and the BNT162b2 mRNA vaccine exhibit Watson-Crick nucleotide complementarity with human coding or noncoding genes. Although they do not share the same complementarity pattern, both may disrupt epigenetic mechanisms in target cells, potentially leading to long-term complications.

ORC1 binds to cis-transcribed RNAs for efficient activation of replication origins

Cells must coordinate the activation of thousands of replication origins dispersed throughout their genome. Active transcription is known to favor the formation of mammalian origins, although the role that RNA plays in this process remains unclear. We show that the ORC1 subunit of the human Origin Recognition Complex interacts with RNAs transcribed from genes with origins in their transcription start sites (TSSs), displaying a positive correlation between RNA binding and origin activity. RNA depletion, or the use of ORC1 RNA-binding mutant, result in inefficient activation of proximal origins, linked to impaired ORC1 chromatin release. ORC1 RNA binding activity resides in its intrinsically disordered region, involved in intra- and inter-molecular interactions, regulation by phosphorylation, and phase-separation. We show that RNA binding favors ORC1 chromatin release, by regulating its phosphorylation and subsequent degradation. Our results unveil a non-coding function of RNA as a dynamic component of the chromatin, orchestrating the activation of replication origins.

The Potential Regulation of A-to-I RNA Editing on Genes in Parkinson's Disease

Parkinson's disease (PD) is characterized by dopaminergic neurodegeneration and an abnormal accumulation of α-synuclein aggregates. A number of genetic factors have been shown to increase the risk of PD. Exploring the underlying molecular mechanisms that mediate PD's transcriptomic diversity can help us understand neurodegenerative pathogenesis. In this study, we identified 9897 A-to-I RNA editing events associated with 6286 genes across 372 PD patients. Of them, 72 RNA editing events altered miRNA binding sites and this may directly affect miRNA regulations of their host genes. However, RNA editing effects on the miRNA regulation of genes are more complex. They can (1) abolish existing miRNA binding sites, which allows miRNAs to regulate other genes; (2) create new miRNA binding sites that may sequester miRNAs from regulating other genes; or (3) occur in the miRNA seed regions and change their targets. The first two processes are also referred to as miRNA competitive binding. In our study, we found 8 RNA editing events that may alter the expression of 1146 other genes via miRNA competition. We also found one RNA editing event that modified a miRNA seed region, which was predicted to disturb the regulation of four genes. Considering the PD-related functions of the affected genes, 25 A-to-I RNA editing biomarkers for PD are proposed, including the 3 editing events in the EIF2AK2, APOL6, and miR-4477b seed regions. These biomarkers may alter the miRNA regulation of 133 PD-related genes. All these analyses reveal the potential mechanisms and regulations of RNA editing in PD pathogenesis.

Effect of mechanical unloading on genome-wide DNA methylation profile of the failing human heart

Heart failure (HF) is characterized by global alterations in myocardial DNA methylation, yet little is known about the epigenetic regulation of the noncoding genome and potential reversibility of DNA methylation with left ventricular assist device (LVAD) therapy. Genome-wide mapping of myocardial DNA methylation in 36 patients with HF at LVAD implantation, 8 patients at LVAD explantation, and 7 nonfailing (NF) donors using a high-density bead array platform identified 2,079 differentially methylated positions (DMPs) in ischemic cardiomyopathy (ICM) and 261 DMPs in nonischemic cardiomyopathy (NICM). LVAD support resulted in normalization of 3.2% of HF-associated DMPs. Methylation-expression correlation analysis yielded several protein-coding genes that are hypomethylated and upregulated (HTRA1, FBXO16, EFCAB13, and AKAP13) or hypermethylated and downregulated (TBX3) in HF. A potentially novel cardiac-specific super-enhancer long noncoding RNA (lncRNA) (LINC00881) is hypermethylated and downregulated in human HF. LINC00881 is an upstream regulator of sarcomere and calcium channel gene expression including MYH6, CACNA1C, and RYR2. LINC00881 knockdown reduces peak calcium amplitude in the beating human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). These data suggest that HF-associated changes in myocardial DNA methylation within coding and noncoding genomes are minimally reversible with mechanical unloading. Epigenetic reprogramming strategies may be necessary to achieve sustained clinical recovery from heart failure.

Long noncoding RNA LEENE promotes angiogenesis and ischemic recovery in diabetes models

Impaired angiogenesis in diabetes is a key process contributing to ischemic diseases such as peripheral arterial disease. Epigenetic mechanisms, including those mediated by long noncoding RNAs (lncRNAs), are crucial links connecting diabetes and the related chronic tissue ischemia. Here we identify the lncRNA that enhances endothelial nitric oxide synthase (eNOS) expression (LEENE) as a regulator of angiogenesis and ischemic response. LEENE expression was decreased in diabetic conditions in cultured endothelial cells (ECs), mouse hind limb muscles, and human arteries. Inhibition of LEENE in human microvascular ECs reduced their angiogenic capacity with a dysregulated angiogenic gene program. Diabetic mice deficient in Leene demonstrated impaired angiogenesis and perfusion following hind limb ischemia. Importantly, overexpression of human LEENE rescued the impaired ischemic response in Leene-knockout mice at tissue functional and single-cell transcriptomic levels. Mechanistically, LEENE RNA promoted transcription of proangiogenic genes in ECs, such as KDR (encoding VEGFR2) and NOS3 (encoding eNOS), potentially by interacting with LEO1, a key component of the RNA polymerase II-associated factor complex and MYC, a crucial transcription factor for angiogenesis. Taken together, our findings demonstrate an essential role for LEENE in the regulation of angiogenesis and tissue perfusion. Functional enhancement of LEENE to restore angiogenesis for tissue repair and regeneration may represent a potential strategy to tackle ischemic vascular diseases.

In silico transcriptional analysis of asymptomatic and severe COVID-19 patients reveals the susceptibility of severe patients to other comorbidities and non-viral pathological conditions

COVID-19 is a severe respiratory disease caused by SARS-CoV-2, a novel human coronavirus. Patients infected with SARS-CoV-2 exhibit heterogeneous symptoms that pose pragmatic hurdles for implementing appropriate therapy and management of the COVID-19 patients and their post-COVID complications. Thus, understanding the impact of infection severity at the molecular level in the host is vital to understand the host response and accordingly it's precise management. In the current study, we performed a comparative transcriptomics analysis of publicly available seven asymptomatic and eight severe COVID-19 patients. Exploratory data analysis employing Principal Component Analysis (PCA) showed the distinct clusters of asymptomatic and severe patients. Subsequently, the differential gene expression analysis using DESeq2 identified 1224 significantly upregulated genes (logFC≥ 1.5, p-adjusted value <0.05) and 268 significantly downregulated genes (logFC≤ −1.5, p-adjusted value <0.05) in severe samples in comparison to asymptomatic samples. Eventually, Gene Set Enrichment Analysis (GSEA) revealed the upregulation of anti-viral and anti-inflammatory pathways, secondary infections, Iron homeostasis, anemia, cardiac-related, etc.; while, downregulation of lipid metabolism, adaptive immune response, translation, recurrent respiratory infections, heme-biosynthetic pathways, etc. Conclusively, these findings provide insight into the enhanced susceptibility of severe COVID-19 patients to other health comorbidities including non-viral pathogenic infections, atherosclerosis, autoinflammatory diseases, anemia, male infertility, etc. owing to the activation of biological processes, pathways and molecular functions associated with them. We anticipate this study will facilitate the researchers in finding efficient therapeutic targets and eventually the clinicians in management of COVID-19 patients and post-COVID-19 effects in them.

circRNA_0067717 promotes paclitaxel resistance in nasopharyngeal carcinoma by acting as a scaffold for TRIM41 and p53

PURPOSE

Circular RNAs (circRNAs) play important roles in tumour progression. This study aimed to explore the mechanism of hsa_circ_0067717 (termed circRNA_0067717) promoting paclitaxel resistance in nasopharyngeal carcinoma (NPC).

METHODS

We assayed CNE-1 and HNE-2 parental cell lines and the corresponding paclitaxel-resistant NPC cell lines using circRNA microarrays. RNA pull-down assay, RNA immunoprecipitation, and RNA fluorescence in situ hybridization were used to identify the molecular mechanisms.

RESULTS

Here, we confirm that circRNA_0067717 is significantly upregulated in NPC paclitaxel-resistant cells and is associated with paclitaxel resistance in NPC. Mechanistically, circRNA_0067717 functions as a scaffold for TRIM41 protein (a ubiquitin E3 ligase) and p53 protein. In nasopharyngeal carcinoma paclitaxel-resistant cells, the highly expressed circRNA_0067717 can bind to more TRIM41 and p53 protein, promoting TRIM41-induced p53 ubiquitination and degradation, resulting in a decrease in p53 protein level. Moreover, the 1-176 nt area of circRNA_0067717 and the 301-425 nt region of circRNA_0067717 are the binding sites for p53 and TRIM41, respectively. The resistance of NPC cells to paclitaxel can be reduced by blocking these binding regions of circRNA_0067717.

CONCLUSION

We demonstrate that circRNA_0067717 acts as a scaffold for TRIM41 and p53, enhancing paclitaxel chemoresistance in NPC by promoting TRIM41-induced p53 degradation via ubiquitination.

Interaction between Long Noncoding RNAs and Syncytin-1/Syncytin-2 Genes and Transcripts: How Noncoding RNAs May Affect Pregnancy in Patients with Systemic Lupus Erythematosus

BACKGROUND

Patients with systemic lupus erythematosus (SLE) often suffer from obstetric complications not necessarily associated with the antiphospholipid syndrome. These events may potentially result from the reduced placental synthesis of the fusogenic proteins syncytin-1 and syncytin-2, observed in women with pregnancy-related disorders. SLE patients have an aberrant noncoding (nc)RNA signature that may in turn dysregulate the expression of syncytin-1 and syncytin-2 during placentation. The aim of this research is to computationally evaluate and characterize the interaction between syncytin-1 and syncytin-2 genes and human ncRNAs and to discuss the potential implications for SLE pregnancy adverse outcomes.

METHODS

The FASTA sequences of the syncytin-1 and syncytin-2 genes were used as inputs to the Ensembl.org library to find any alignments with human ncRNA genes and their transcripts, which were characterized for their tissue expression, regulatory activity on adjacent genes, biological pathways, and potential association with human disease.

RESULTS

BLASTN analysis revealed a total of 100 hits with human long ncRNAs (lncRNAs) for the syncytin-1 and syncytin-2 genes, with median alignment scores of 151 and 66.7, respectively. Only lncRNAs TP53TG1, TTTY14, and ENSG00000273328 were reported to be expressed in placental tissue. Dysregulated expression of lncRNAs TP53TG1, LINC01239, and LINC01320 found in this analysis has previously been described in SLE patients as well as in women with a high-risk pregnancy. In addition, some of the genes adjacent to lncRNAs aligned with syncytin-1 or syncytin-2 in a regulatory region might increase the risk of pregnancy complications or SLE.

CONCLUSIONS

This is the first computational study showing alignments between syncytin-1 and syncytin-2 genes and human lncRNAs. Whether this mechanism affects syncytiotrophoblast morphogenesis in SLE females is unknown and requires further investigation.

The PRALINE database: protein and Rna humAn singLe nucleotIde variaNts in condEnsates

SUMMARY

Biological condensates are membraneless organelles with different material properties. Proteins and RNAs are the main components, but most of their interactions are still unknown. Here, we introduce PRALINE, a database for the interrogation of proteins and RNAs contained in stress granules, processing bodies and other assemblies including droplets and amyloids. PRALINE provides information about the predicted and experimentally validated protein-protein, protein-RNA and RNA-RNA interactions. For proteins, it reports the liquid-liquid phase separation and liquid-solid phase separation propensities. For RNAs, it provides information on predicted secondary structure content. PRALINE shows detailed information on human single-nucleotide variants, their clinical significance and presence in protein and RNA binding sites, and how they can affect condensates' physical properties.

AVAILABILITY AND IMPLEMENTATION

PRALINE is freely accessible on the web at http://praline.tartaglialab.com.

fRNC: Uncovering the dynamic and condition-specific RBP-ncRNA circuits from multi-omics data

The RNA binding protein (RBP) and non-coding RNA (ncRNA) interacting networks are increasingly recognized as the main mechanism in gene regulation, and are tightly associated with cellular malfunction and disease. Here, we present fRNC, a systems biology tool to uncover the dynamic spectrum of RBP-ncRNA circuits (RNC) by integrating transcriptomics, interactomics and proteomics data. fRNC constructs the RBP-ncRNA network derived from CLIP-seq or PARE experiments. Given scoring on nodes and edges according to differential analysis of expression data, it finds an RNC containing global maximum significant RBPs and ncRNAs. Alternatively, it can also capture the locally maximum scoring RNC according to user-defined starting nodes with the greedy search. When compared with existing tools, fRNC can detect more accurate and robust sub-network with scalability. As shown in the cases of esophageal carcinoma, breast cancer and Alzheimer's disease, fRNC enables users to analyze the collective behaviors between RBP and the interacting ncRNAs, and reveal novel insights into the disease-associated processes. The fRNC R package is available at https://github.com/BioinformaticsSTU/fRNC.

LINC02870 facilitates SNAIL translation to promote hepatocellular carcinoma progression

Exploring the roles of long noncoding RNAs (lncRNAs) in tumorigenesis and metastasis could contribute to the recognition of novel diagnostic and therapeutic targets. LINC02870 is a novel lncRNA, whose role in tumors has not been reported. Herein, we focused on the function and mechanism of LINC02870 in human hepatocellular carcinoma (HCC). We first carried out a pan-cancer study of LINC02870 expression and its relationship to prognosis, and LINC02870 was determined to be a possible oncogene in HCC. Upregulated expressions of LINC02870 were also found in our HCC samples compared to the para-tumor samples. Moreover, overexpression of LINC02870 promoted the growth, migration, and invasion of HCC cells. Subsequently, binding proteins of LINC02870 were identified by a number of in silico analyses, including correlation analysis, signaling network analysis, and survival analysis. Intriguingly, the most promising binding protein of LINC02870 was predicted and confirmed to be eukaryotic translation initiation factor 4 gamma 1 (EIF4G1), an important component of the eukaryotic translation initiation factor 4F complex that initiates cap-dependent translation. Further investigation showed that LINC02870 increased the translation of SNAIL to induce malignant phenotypes in HCC cells. Additionally, HCC patients with higher expression levels of LINC02870 and EIF4G1 had shorter survival times than those with lower expression levels. Thus, our findings suggested that LINC02870 induced SNAIL translation and correlated with poor prognosis and tumor progression in HCC.

SPTBN1 abrogates renal clear cell carcinoma progression via glycolysis reprogramming in a GPT2-dependent manner

BACKGROUND

Renal clear cell carcinoma (ccRCC) is the most prevalent tumors worldwide. Discovering effective biomarkers is essential to monitor the prognosis and provide alternative clinical options. SPTBN1 is implicated in various cancerous processes. However, its role in ccRCC remains unelucidated. This study intends to explore the biological function and mechanism of SPTBN1 in ccRCC.

METHODS

Single-cell and bulk RNA-seq, tissue microarray, real-time quantitative PCR, and western blotting were applied to verify the expression and predictive value of SPTBN1 in ccRCC. Gain or loss of functional ccRCC cell line models were constructed, and in vitro and in vivo assays were performed to elucidate its tumorigenic phenotypes. Actinomycin D experiment, RNA immunoprecipitation (RIP), specific inhibitors, and rescue experiments were carried out to define the molecular mechanisms.

RESULTS

SPTBN1 was down-regulated in ccRCC and knockdown of SPTBN1 displayed a remarkably oncogenic role both in vitro and in vivo; while overexpressing SPTBN1 reversed this effect. SPTBN1 mediated ccRCC progression via the pathway of glutamate pyruvate transaminase 2 (GPT2)-dependent glycolysis. The expression of GPT2 was significantly negatively correlated with that of SPTBN1. As an RNA binding protein SPTBN1, regulated the mRNA stability of GPT2.

CONCLUSION

Our research demonstrated that SPTBN1 is significantly down-regulated in ccRCC. SPTBN1 knockdown promotes ccRCC progression via activating GPT2-dependent glycolysis. SPTBN1 may serve as a therapeutic target for the treatment of ccRCC.

CRIT: Identifying RNA-binding protein regulator in circRNA life cycle via non-negative matrix factorization

Circular RNAs (circRNAs) are endogenous non-coding RNAs that regulate gene expression and participate in carcinogenesis. However, the RNA-binding proteins (RBPs) involved in circRNAs biogenesis and modulation remain largely unclear. We developed the circRNA regulator identification tool (CRIT), a non-negative matrix-factorization-based pipeline to identify regulating RBPs in cancers. CRIT uncovered 73 novel regulators across thousands of samples by effectively leveraging genomics data and functional annotations. We demonstrated that known RBPs involved in circRNA control are significantly enriched in these predictions. Analysis of circRNA-RBP interactions using two large cross-linking immunoprecipitation (CLIP) databases, we validated the consistency between CRIT prediction and the CLIP experiments. Furthermore, newly discovered RBPs are functionally connected with authentic circRNA regulators by various biological associations, such as physical interaction, similar binding motifs, common transcription factor modulation, and co-expression. When analyzing RNA sequencing (RNA-seq) datasets after short hairpin RNA (shRNA)/small interfering RNA (siRNA) knockdown, we found several novel RBPs that can affect global circRNA expression, which strengthens their role in the circRNA life cycle. The above evidence provided independent confirmation that CRIT is a useful tool to capture RBPs in circRNA processing. Finally, we show that authentic regulators are more likely the core splicing proteins and peripheral factors and usually harbor more alterations in the vast majority of cancers.

Identification of Novel Kinases of Tau Using Fluorescence Complementation Mass Spectrometry (FCMS)

Hyperphosphorylation of the microtubule-associated protein Tau is a major hallmark of Alzheimer's disease and other tauopathies. Understanding the protein kinases that phosphorylate Tau is critical for the development of new drugs that target Tau phosphorylation. At present, the repertoire of the Tau kinases remains incomplete, and methods to uncover novel upstream protein kinases are still limited. Here, we apply our newly developed proteomic strategy, fluorescence complementation mass spectrometry, to identify novel kinase candidates of Tau. By constructing Tau- and kinase-fluorescent fragment library, we detected 59 Tau-associated kinases, including 23 known kinases of Tau and 36 novel candidate kinases. In the validation phase using in vitro phosphorylation, among 15 candidate kinases we attempted to purify and test, four candidate kinases, OXSR1 (oxidative-stress responsive gene 1), DAPK2 (death-associated protein kinase 2), CSK (C-terminal SRC kinase), and ZAP70 (zeta chain of T-cell receptor-associated protein kinase 70), displayed the ability to phosphorylate Tau in time-course experiments. Furthermore, coexpression of these four kinases along with Tau increased the phosphorylation of Tau in human neuroglioma H4 cells. We demonstrate that fluorescence complementation mass spectrometry is a powerful proteomic strategy to systematically identify potential kinases that can phosphorylate Tau in cells. Our discovery of new candidate kinases of Tau can present new opportunities for developing Alzheimer's disease therapeutic strategies.

LncRNA PVT-1 promotes osteosarcoma cancer stem-like properties through direct interaction with TRIM28 and TSC2 ubiquitination

Osteosarcoma, the most common pediatric bone tumor, is an aggressive heterogeneous malignancy defined by complex chromosomal aberrations. Overall survival rates remain at ~70%, but patients with chemoresistant or metastatic disease have extremely poor outcomes of <30%. A subgroup of tumors harbor amplification of chromosome 8q24.2 and increased expression of the oncogenic long noncoding RNA (lncRNA) Plasmacytoma Variant Translocation-1 (PVT-1), which is associated with an extremely poor clinical prognosis. This study demonstrates that PVT-1 is critical for osteosarcoma tumor-initiation potential. Chromatin Hybridization by RNA Purification analysis identified Tripartite-Motif Containing Family 28 (TRIM28) as a novel PVT-1 binding partner. Mechanistically, co-immunoprecipitation studies showed the PVT-1/TRIM28 complex binds and increases SUMOylation of phosphatidylinositol 3-kinase catalytic subunit type 3 (Vps34), which leads to enhanced ubiquitination and degradation of tumor suppressor complex 2 (TSC2), thus contributing to increased self-renewal and stem cell phenotypes. Furthermore, we identified that osteosarcoma cells with increased PVT-1 have enhanced sensitivity to the SUMOylation inhibitor, TAK-981. Altogether, this study elucidated a role for PVT-1 in the enhancement of cancer stem-like behaviors, including migration and invasion, in osteosarcoma, and identified the novel PVT-1/TRIM28 axis signaling cascade as a potential therapeutic target for osteosarcoma treatment.

Identification of a Novel Long Non-coding RNA, lnc-ATMIN-4:2, and its Clinicopathological and Prognostic Significance in Advanced Gastric Cancer

BACKGROUND/AIM

Long non-coding RNAs (lncRNAs) are emerging as significant regulators of gene expression and a novel promising biomarker for cancer diagnosis and prognosis. This study identified a novel, differentially expressed lncRNA in advanced gastric cancer (AGC), Inc-ATMIN-4:2, and evaluated its clinicopathological and prognostic significance.

PATIENTS AND METHODS

Whole transcriptome sequencing was performed to identify differentially expressed lncRNAs in AGC tissue samples. We also analyzed lnc-ATMIN-4:2 expression in 317 patients with AGC using RNA in situ hybridization.

RESULTS

High (>30 dots) lnc-ATMIN-4:2 expression significantly correlated with younger age, poorly differentiated histology, diffuse type, deeper invasion depth, perineural invasion, lymph node metastasis, and higher stage group. In addition, high lnc-ATMIN-4:2 expression was significantly associated with worse overall survival in patients with AGC.

CONCLUSION

This study elucidated the significance of lncRNAs in AGC and indicated the value of lnc-ATMIN-4:2 expression as a predictive biomarker for the overall survival of patients with AGC.

RNA splicing regulators play critical roles in neurogenesis

Alternative RNA splicing increases transcript diversity in different cell types and under varying conditions. It is executed with the help of RNA splicing regulators (RSRs), which are operationally defined as RNA-binding proteins (RBPs) that regulate alternative splicing, but not directly catalyzing the chemical reactions of splicing. By systematically searching for RBPs and manually identifying those that regulate splicing, we curated 305 RSRs in the human genome. Surprisingly, most of the RSRs are involved in neurogenesis. Among these RSRs, we focus on nine families (PTBP, NOVA, RBFOX, ELAVL, CELF, DBHS, MSI, PCBP, and MBNL) that play essential roles in the neurogenic pathway. A better understanding of their functions will provide novel insights into the role of splicing in brain development, health, and disease. This comprehensive review serves as a stepping-stone to explore the diverse and complex set of RSRs as fundamental regulators of neural development. This article is categorized under: RNA-Based Catalysis &gt; RNA Catalysis in Splicing and Translation RNA Interactions with Proteins and Other Molecules &gt; Protein-RNA Interactions: Functional Implications RNA Processing &gt; Splicing Regulation/Alternative Splicing.

Transcriptomic profiling of Indian breast cancer patients revealed subtype-specific mRNA and lncRNA signatures

Breast cancer (BC) is one of the leading causes of cancer-associated death in women. Despite the progress in therapeutic regimen, resistance and recurrence of breast cancer have affected the overall survival of patients. The present signatures, such as PAM50 and Oncotype DX, do not segregate the Indian breast samples based on molecular subtypes. This study aims at finding signatures of long noncoding RNA (lncRNA) and mRNA in Indian breast cancer patients using RNA-seq. We have analyzed the survival based on the menopausal and hormone status of 380 Indian breast cancer patients, and of these, we have sequenced and analyzed matched tumor–normal transcriptome of 17 (pre- and postmenopausal) Indian breast cancer patients representing six different subtypes, namely, four patients in triple-positive, three patients in estrogen receptor–positive (ER+ve), three patients in estrogen and progesterone receptors–positive (ER+ve, PR+ve), two patients in human epidermal growth factor receptor (Her2+ve), three patients in triple-negative, and one patient in ER+ve and Her2+ve subtypes. We have identified a 25 mRNA–27 lncRNA gene set, which segregated the subtypes in our data. A pathway analysis of the differentially expressed genes revealed downregulated ECM interaction and upregulated immune regulation, cell cycle, DNA damage response and repair, and telomere elongation in premenopausal women. Postmenopausal women showed downregulated metabolism, innate immune system, upregulated translation, sumoylation, and AKT2 activation. A Kaplan–Meier survival analysis revealed that menopausal status, grade of the tumor, and hormonal status displayed statistically significant effects (p < 0.05) on the risk of mortality due to breast cancer. Her2+ve patients showed low overall survival. One of the unique lncRNA-mRNA pairs specific to the EP-subtype, SNHG12 and EPB41, showed interaction, which correlates with their expression level; SNHG12 is downregulated and EPB41 is upregulated in EP samples.

Oxidative Stress and Its Modulation by Ladostigil Alter the Expression of Abundant Long Non-Coding RNAs in SH-SY5Y Cells

Neurodegenerative disorders, brain injury, and the decline in cognitive function with aging are accompanied by a reduced capacity of cells in the brain to cope with oxidative stress and inflammation. In this study, we focused on the response to oxidative stress in SH-SY5Y, a human neuroblastoma cell line. We monitored the viability of the cells in the presence of oxidative stress. Such stress was induced by hydrogen peroxide or by Sin1 (3-morpholinosydnonimine) that generates reactive oxygen and nitrogen species (ROS and RNS). Both stressors caused significant cell death. Our results from the RNA-seq experiments show that SH-SY5Y cells treated with Sin1 for 24 h resulted in 94 differently expressed long non-coding RNAs (lncRNAs), including many abundant ones. Among the abundant lncRNAs that were upregulated by exposing the cells to Sin1 were those implicated in redox homeostasis, energy metabolism, and neurodegenerative diseases (e.g., MALAT1, MIAT, GABPB1-AS1, NEAT1, MIAT, GABPB1-AS1, and HAND2-AS1). Another group of abundant lncRNAs that were significantly altered under oxidative stress included cancer-related SNHG family members. We tested the impact of ladostigil, a bifunctional reagent with antioxidant and anti-inflammatory properties, on the lncRNA expression levels. Ladostigil was previously shown to enhance learning and memory in the brains of elderly rats. In SH-SY5Y cells, several lncRNAs involved in transcription regulation and the chromatin structure were significantly induced by ladostigil. We anticipate that these poorly studied lncRNAs may act as enhancers (eRNA), regulating transcription and splicing, and in competition for miRNA binding (ceRNA). We found that the induction of abundant lncRNAs, such as MALAT1, NEAT-1, MIAT, and SHNG12, by the Sin1 oxidative stress paradigm specifies only the undifferentiated cell state. We conclude that a global alteration in the lncRNA profiles upon stress in SH-SY5Y may shift cell homeostasis and is an attractive in vitro system to characterize drugs that impact the redox state of the cells and their viability.

Circulating U13 Small Nucleolar RNA as a Potential Biomarker in Huntington's Disease: A Pilot Study

Plasma small RNAs have been recently explored as biomarkers in Huntington’s disease (HD). We performed an exploratory study on nine HD patients, eight healthy subjects (HS), and five psychiatric patients (PP; to control for iatrogenic confounder effects) through an Affymetrix-Gene-Chip-miRNA-Array. We validated the results in an independent population of 23 HD, 15 pre-HD, 24 PP, 28 Alzheimer’s disease (AD) patients (to control the disease-specificity) and 22 HS through real-time PCR. The microarray results showed higher levels of U13 small nucleolar RNA (SNORD13) in HD patients than controls (fold change 1.54, p = 0.003 HD vs. HS, and 1.44, p = 0.0026 HD vs. PP). In the validation population, a significant increase emerged with respect to both pre-HD and the control groups (p < 0.0001). SNORD13 correlated with the status of the mutant huntingtin carrier (r = 0.73; p < 0.001) and the disease duration (r = 0.59; p = 0.003). The receiver operating characteristic (ROC) curve analysis showed the high accuracy of SNORD13 in discriminating HD patients from other groups (AUC = 0.963). An interactome and pathway analysis on SNORD13 revealed enrichments for factors relevant to HD pathogenesis. We report the unprecedented finding of a potential disease-specific role of SNORD13 in HD. It seems to peripherally report a ‘tipping point’ in the pathogenic cascade at the neuronal level.

Computational tools to study RNA-protein complexes

RNA is the key player in many cellular processes such as signal transduction, replication, transport, cell division, transcription, and translation. These diverse functions are accomplished through interactions of RNA with proteins. However, protein–RNA interactions are still poorly derstood in contrast to protein–protein and protein–DNA interactions. This knowledge gap can be attributed to the limited availability of protein-RNA structures along with the experimental difficulties in studying these complexes. Recent progress in computational resources has expanded the number of tools available for studying protein-RNA interactions at various molecular levels. These include tools for predicting interacting residues from primary sequences, modelling of protein-RNA complexes, predicting hotspots in these complexes and insights into derstanding in the dynamics of their interactions. Each of these tools has its strengths and limitations, which makes it significant to select an optimal approach for the question of interest. Here we present a mini review of computational tools to study different aspects of protein-RNA interactions, with focus on overall application, development of the field and the future perspectives.

Expansion and functional analysis of the SR-related protein family across the domains of life

Serine/arginine-rich (SR) proteins comprise a family of proteins that is predominantly found in eukaryotes and plays a prominent role in RNA splicing. A characteristic feature of SR proteins is the presence of an S/R-rich low-complexity domain (RS domain), often in conjunction with spatially distinct RNA recognition motifs (RRMs). To date, 52 human proteins have been classified as SR or SR-related proteins. Here, using an unbiased series of composition criteria together with enrichment for known RNA binding activity, we identified >100 putative SR-related proteins in the human proteome. This method recovers known SR and SR-related proteins with high sensitivity (∼94%), yet identifies a number of additional proteins with many of the hallmark features of true SR-related proteins. Newly identified SR-related proteins display slightly different amino acid compositions yet similar levels of post-translational modification, suggesting that these new SR-related candidates are regulated in vivo and functionally important. Furthermore, candidate SR-related proteins with known RNA-binding activity (but not currently recognized as SR-related proteins) are nevertheless strongly associated with a variety of functions related to mRNA splicing and nuclear speckles. Finally, we applied our SR search method to all available reference proteomes, and provide maps of RS domains and Pfam annotations for all putative SR-related proteins as a resource. Together, these results expand the set of SR-related proteins in humans, and identify the most common functions associated with SR-related proteins across all domains of life.

Comprehensive analysis of the expression of N6-methyladenosine RNA methylation regulators in pulmonary artery hypertension

Background: Pulmonary arterial hypertension (PAH) is a progressive disease characterized by pulmonary vascular remodeling. The development of PAH involves N6-methyladenosine (m6A) modification. However, the functional role of m6A regulators in PAH and the underlying regulatory mechanisms remain unknown so far.

Methods: Microarray data (GSE149713) for monocrotaline induced PAH (MCT-PAH) rat models were downloaded and screened for differentially expressed genes (DEGs) and m6A regulators. Next, we screened for differentially expressed m6A regulators in endothelial cells (ECs), smooth muscle cells (SMCs), fibroblasts, interstitial macrophages, NK cells, B cells, T cells, regulatory T cells (Tregs) using scRNA sequencing data. The target DEGs of m6A regulators in ECs, SMCs, fibroblasts, and Tregs were functionally annotated using the Gene Ontology (GO) functional analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. In addition, the cellular interaction analysis was performed to reveal the receptor—ligand pairs regulated by m6A regulators. Pseudo-time trajectory analyses were performed and a ceRNA network of lncRNAs-miRNAs-mRNAs was constructed in SMCs. Furthermore, the RNA transcriptome sequencing data for the SMCs isolated from idiopathic PAH (IPAH) patients (GSE144274) were validated for differentially expressed m6A regulators. Moreover, the HNRNPA2B1 levels in the lung samples from PAH patients and MCT-PAH were determined using immunohistochemistry.

Results: The m6A regulators were observed to be dysregulated in PAH. HNRNPA2B1expression level was increased in the PASMCs of scRNAs and IPAH patients. The target DEGs of HNRNPA2B1 were enriched in the regulation of muscle cell differentiation and vasculature development in PASMCs. The HNRNPA2B1 expression levels determined were consistent with the proliferation-related and collagen synthesis-related gene COL4A1. Moreover, the predicted transcription factors (TFs) foxd2/3 and NFκB could be involved in the regulation of HNRNPA2B1. HNRNPA2B1 might be regulating SMCs proliferation and phenotypic transition via rno-miR-330–3p/TGFβR3 and rno-miR-125a-3p/slc39a1. In addition, HNRNPA2B1 was observed to be highly expressed in the lung samples from MCT-PAH rat models and patients with PAH.

Conclusion: In summary, the present study identified certain key functional m6A regulators that are involved in pulmonary vascular remodeling. The investigation of m6A patterns might be promising and provide biomarkers for diagnosis and treatment of PAH in the future.

A high-throughput approach to predict A-to-I effects on RNA structure indicates a change of double-stranded content in non-coding RNAs

RNA molecules undergo a number of chemical modifications whose effects can alter their structure and molecular interactions. Previous studies have shown that RNA editing can impact the formation of ribonucleoprotein complexes and influence the assembly of membrane-less organelles such as stress-granules. For instance, N6-methyladenosine (m6A) enhances SG formation and N1-methyladenosine (m1A) prevents their transition to solid-like aggregates. Yet, very little is known about adenosine to inosine (A-to-I) modification that is very abundant in human cells and not only impacts mRNAs but also non-coding RNAs. Here, we built the CROSSalive predictor of A-to-I effects on RNA structure based on high-throughput in-cell experiments. Our method shows an accuracy of 90% in predicting the single and double-stranded content of transcripts and identifies a general enrichment of double-stranded regions caused by A-to-I in long intergenic non-coding RNAs (lincRNAs). For the individual cases of NEAT1, NORAD and XIST, we investigated the relationship between A-to-I editing and interactions with RNA-binding proteins using available CLIP data and catRAPID predictions. We found that A-to-I editing is linked to alteration of interaction sites with proteins involved in phase-separation, which suggests that RNP assembly can be influenced by A-to-I. CROSSalive is available at http://service.tartaglialab.com/new_submission/crossalive. This article is protected by copyright. All rights reserved.

lncRNA CYTOR promotes aberrant glycolysis and mitochondrial respiration via HNRNPC-mediated ZEB1 stabilization in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC), the most common malignancy of the oral and maxillofacial region, severely affects human health. However, current treatments for OSCC commonly show only a ~60% 5-year survival rate of patients with distant metastases, indicating an urgent need for targeted treatments for patients with advanced metastases. Here, we report a survival-related long non-coding RNA, CYTOR, which is highly expressed in the lesions of oral cancer patients. We found that CYTOR can promote both migration and invasion in oral cancer cells as well as the epithelial-mesenchymal transition (EMT). RNA-sequencing of CYTOR-knockdown oral cancer cells revealed that CYTOR can regulate mitochondrial respiration and RNA splicing. Mechanistically, we found that nuclear-localized CYTOR interacts with HNRNPC, resulting in stabilization of ZEB1 mRNAs by inhibiting the nondegradative ubiquitination of HNRNPC. By synthesizing CYTOR-targeting small interfering RNAs (siRNAs) encapsulated in Nanoscale Metal Organic Frameworks (NMOFs), we demonstrate the targeted suppression of CYTOR to inhibit invasion and metastasis of oral cancer cells in a nude mouse model. Cumulatively, this study reveals the potential role of the CYTOR-HNRNPC-ZEB1 axis in regulating mitochondrial metabolism and glycolysis of oral cancer cells, and illustrates the effective use of lncRNA targeting in anti-metastatic cancer therapies.

Probing TDP-43 condensation using an in silico designed aptamer

Aptamers are artificial oligonucleotides binding to specific molecular targets. They have a promising role in therapeutics and diagnostics but are often difficult to design. Here, we exploited the catRAPID algorithm to generate aptamers targeting TAR DNA-binding protein 43 (TDP-43), whose aggregation is associated with Amyotrophic Lateral Sclerosis. On the pathway to forming insoluble inclusions, TDP-43 adopts a heterogeneous population of assemblies, many smaller than the diffraction-limit of light. We demonstrated that our aptamers bind TDP-43 and used the tightest interactor, Apt-1, as a probe to visualize TDP-43 condensates with super-resolution microscopy. At a resolution of 10 nanometers, we tracked TDP-43 oligomers undetectable by standard approaches. In cells, Apt-1 interacts with both diffuse and condensed forms of TDP-43, indicating that Apt-1 can be exploited to follow TDP-43 phase transition. The de novo generation of aptamers and their use for microscopy opens a new page to study protein condensation.

RPpocket: An RNA–Protein Intuitive Database with RNA Pocket Topology Resources

RNA–protein complexes regulate a variety of biological functions. Thus, it is essential to explore and visualize RNA–protein structural interaction features, especially pocket interactions. In this work, we develop an easy-to-use bioinformatics resource: RPpocket. This database provides RNA–protein complex interactions based on sequence, secondary structure, and pocket topology analysis. We extracted 793 pockets from 74 non-redundant RNA–protein structures. Then, we calculated the binding- and non-binding pocket topological properties and analyzed the binding mechanism of the RNA–protein complex. The results showed that the binding pockets were more extended than the non-binding pockets. We also found that long-range forces were the main interaction for RNA–protein recognition, while short-range forces strengthened and optimized the binding. RPpocket could facilitate RNA–protein engineering for biological or medical applications.

Mixed model-based eQTL analysis reveals lncRNAs associated with regulation of genes involved in sex determination and spermatogenesis: The key to understanding human gender imbalance

BACKGROUND

An imbalance in the prenatal sex ratio in humans may be due to several factors affecting sperm physiology, including genetic features. In this study, we conducted a transcriptome-wide analysis of expression quantitative trait loci (eQTLs) to identify target genes associated with previously described QTLs associated with gender imbalance.

METHODS

A mixed model explaining polygenic effects by genomic covariance among individuals was used to identify the eQTLs using gene expression and genotype data from 462 European/African individuals.

RESULTS

Eight eGenes were associated with four QTLs (P < 4.00 × 10^-5), with strong associations found (P < 4.00 × 10^-8) between rs2485007 and eGenes ANKRD26P3 (P = 3.40 × 10^-9) and LINC00421 (P = 1.35 × 10^-9). ANKRD26P3 and LINC00421 are both lncRNAs associated with the control of testis-dominant genes PELP1, TAF15, NANOG, TEX14, TCF3, ZNF433, ZNF555, TEX37, FATE1, TCP11, and CYLC2 and Y-linked genes SRY and ZFY, as well as several genes with roles in spermatogenesis (ODF1, SPATC1, SPATA3, SPATA31E1, SPERT, SPATA16, MOSPD1, SPATA24, and SPO11) and sex determination (SOX family genes).

CONCLUSIONS

The above eGenes contribute directly or indirectly to gene regulation for sex determination and spermatogenesis, thereby serving as important functional clues for gender-biased selection.

Sequence complementarity between human noncoding RNAs and SARS-CoV-2 genes: What are the implications for human health?

Objectives

To investigate in silico the presence of nucleotide sequence complementarity between the RNA genome of Severe Acute Respiratory Syndrome CoronaVirus-2 (SARS-CoV-2) and human non-coding (nc)RNA genes.

Methods

The FASTA sequence (NC_045512.2) of each of the 11 SARS-CoV-2 isolate Wuhan-Hu-1 genes was retrieved from NCBI.nlm.nih.gov/gene and the Ensembl.org library interrogated for any base-pair match with human ncRNA genes. SARS-CoV-2 gene-matched human ncRNAs were screened for functional activity using bioinformatic analysis. Finally, associations between identified ncRNAs and human diseases were searched in GWAS databases.

Results

A total of 252 matches were found between the nucleotide sequence of SARS-CoV-2 genes and human ncRNAs. With the exception of two small nuclear RNAs, all of them were long non-coding (lnc)RNAs expressed mainly in testis and central nervous system under physiological conditions. The percentage of alignment ranged from 91.30% to 100% with a mean nucleotide alignment length of 17.5 ± 2.4. Thirty-three (13.09%) of them contained predicted R-loop forming sequences, but none of these intersected the complementary sequences of SARS-CoV-2. However, in 31 cases matches fell on ncRNA regulatory sites, whose adjacent coding genes are mostly involved in cancer, immunological and neurological pathways. Similarly, several polymorphic variants of detected non-coding genes have been associated with neuropsychiatric and proliferative disorders.

Conclusion

This pivotal in silico study shows that SARS-CoV-2 genes have Watson-Crick nucleotide complementarity to human ncRNA sequences, potentially disrupting ncRNA epigenetic control of target genes. It remains to be elucidated whether this could result in the development of human disease in the long term.

Identification of long non-coding RNAs and RNA binding proteins in breast cancer subtypes

Breast cancer is a heterogeneous disease classified into four main subtypes with different clinical outcomes, such as patient survival, prognosis, and relapse. Current genetic tests for the differential diagnosis of BC subtypes showed a poor reproducibility. Therefore, an early and correct diagnosis of molecular subtypes is one of the challenges in the clinic. In the present study, we identified differentially expressed genes, long non-coding RNAs and RNA binding proteins for each BC subtype from a public dataset applying bioinformatics algorithms. In addition, we investigated their interactions and we proposed interacting biomarkers as potential signature specific for each BC subtype. We found a network of only 2 RBPs (RBM20 and PCDH20) and 2 genes (HOXB3 and RASSF7) for luminal A, a network of 21 RBPs and 53 genes for luminal B, a HER2-specific network of 14 RBPs and 30 genes, and a network of 54 RBPs and 302 genes for basal BC. We validated the signature considering their expression levels on an independent dataset evaluating their ability to classify the different molecular subtypes with a machine learning approach. Overall, we achieved good performances of classification with an accuracy >0.80. In addition, we found some interesting novel prognostic biomarkers such as RASSF7 for luminal A, DCTPP1 for luminal B, DHRS11, KLC3, NAGS, and TMEM98 for HER2, and ABHD14A and ADSSL1 for basal. The findings could provide preliminary evidence to identify putative new prognostic biomarkers and therapeutic targets for individual breast cancer subtypes.