Integrative Analysis of Somatic Mutations in Non-coding Regions Altering RNA Secondary Structures in Cancer Genomes

The Secondary Structure of Potato Spindle Tuber Viroid Determines Its Infectivity in Nicotiana benthamiana

The function of RNAs is determined by their structure. However, studying the relationship between RNA structure and function often requires altering RNA sequences to modify the structures, which leads to the neglect of the importance of RNA sequences themselves. In our research, we utilized potato spindle tuber viroid (PSTVd), a circular-form non-coding infectious RNA, as a model with which to investigate the role of a specific rod-like structure in RNA function. By generating linear RNA transcripts with different start sites, we established 12 PSTVd forms with different secondary structures while maintaining the same sequence. The RNA secondary structures were predicted using the mfold tool and validated through native PAGE gel electrophoresis after in vitro RNA folding. Analysis using plant infection assays revealed that the formation of a correct rod-like structure is crucial for the successful infection of PSTVd. Interestingly, the inability of PSTVd forms with non-rod-like structures to infect plants could be partially compensated by increasing the amount of linear viroid RNA transcripts, suggesting the existence of additional RNA secondary structures, such as the correct rod-like structure, alongside the dominant structure in the RNA inoculum of these forms. Our study demonstrates the critical role of RNA secondary structures in determining the function of infectious RNAs.

An insight into the potential role of LINC00968 in luminal breast cancer: Case-control study and bioinformatics analysis

Background

Luminal A and B subtypes of breast cancer (BC) comprises up to 70% of all BC patients. LncRNAs can affect many biological and pathological processes, and dysregulation of them is related to human cancers. The potential role of lncRNA LINC00968 in luminal BC is still unclear.

Materials and methods

We analyzed the LINC00968 expression across 44 paired luminal BC tissues from the TCGA-BRCA RNA sequencing dataset. Besides, we used the GEPIA2 web server and GENEVESTIGATOR software, as well. Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) assay was performed to confirm the LINC00968 expression in 71 paired luminal BC tissues and two luminal A cell lines (MCF7 and T47D). Moreover, to better understanding the potential role of LINC00968 in luminal BC, computational data analyses including co-expression analysis, functional annotation analysis, and genetic alteration analysis have been done.

Results

The results of data analyses retrieved from BRCA dataset and databases revealed the significant downregulation of LINC00968 in luminal A and B BC. Also, the results of qRT-PCR in luminal BC tissues and cell lines confirmed the earlier data. LINC00968 expression was negatively associated with tumor stage and lymph node metastasis. Additionally, functional annotation analyses revealed that LINC00968 might be involved in vascular development and angiogenesis, extracellular matrix organization, and cell motility and migration. LINC00968 might play role in some cancer-related signaling pathways.

Conclusion

Our study found that downregulation of LINC00968 might promote tumorigenesis, invasion, and metastasis of luminal BC.

Genetic variants in the calcium signaling pathway participate in the pathogenesis of colorectal cancer through the tumor microenvironment

Background

Cancer risk is influenced by calcium signaling in intracellular and intercellular signaling pathways. However, the relationship between the calcium signaling pathway and colorectal cancer risk remains unknown. We aim to evaluate the role of genetic variants in calcium signaling pathway genes in colorectal cancer risk through the tumor microenvironment.

Methods

An analysis of genetic variants in the calcium signaling pathway was conducted using a case-control study that included 1150 colorectal cancer patients and 1342 non-cancer patients. Using the regression model, we assessed whether single-nucleotide polymorphisms (SNPs) increase the risk of colorectal cancer. We also performed a dual luciferase reporter gene assay using HCT116 cell lines and DLD1 cell lines to demonstrate the regulatory relationship between SNP and candidate risk gene. We evaluated the expression of candidate risk gene in different populations. In addition, we also evaluated candidate risk gene and 22 immune cells correlation studies.

Results

There was a significant association between the PDE1C rs12538364 T allele and colorectal cancer risk [odds ratio (OR) = 1.57, 95% confidence interval (CI) = 1.30 - 1.90, P = 3.07 × 10^-6, P _FDR = 0.004]. Mutation of intron region rs1538364 C to T locus reduces promoter activity of PDE1C in DLD1 and HCT116 cell lines (P < 0.05). We identified that PDE1C is significantly down-regulated in colorectal cancer, closely associated with 22 immune cells. Finally, we found that PDE1C could be the biomarker for individual immunotherapy of colorectal cancer.

Conclusion

According to our findings, PDE1C may be a key factor contributing to colorectal cancer, thus improving individual immunotherapy for the disease. The potential mechanism by which polymorphisms in the calcium signaling pathway genes may participate in the pathogenesis of colorectal cancer through the tumor microenvironment.

RNA Secondary Structure Alteration Caused by Single Nucleotide Variants

A point mutation in coding RNA can cause not only an amino acid substitution but also a dynamic change of RNA secondary structure, leading to a dysfunctional RNA. Although in silico structure prediction has been used to detect structure-disrupting point mutations known as riboSNitches, exhaustive simulation of long RNAs is needed to detect a significant enrichment or depletion of riboSNitches in functional RNAs. Here, we have developed a novel algorithm Radiam (RNA secondary structure Analysis with Deletion, Insertion, And substitution Mutations) for a comprehensive riboSNitch analysis of long RNAs. Radiam is based on the ParasoR framework, which efficiently computes local RNA secondary structures for long RNAs. ParasoR can compute a variety of structure scores over globally consistent structures with maximal span constraints for the base pair distance. Using the reusable structure database made by ParasoR, Radiam performs an efficient recomputation of the secondary structures for mutated sequences. An exhaustive simulation of Radiam is expected to find reliable riboSNitch candidates on long RNAs by evaluating their statistical significance in terms of the change of local structure stability.

RNAsmc: A integrated tool for comparing RNA secondary structure and evaluating allosteric effects

RNA structure plays a crucial role in gene regulation, in RNA stability and the essential biological processes. RNA secondary structure (RSS) motifs are the basic building blocks for investigating the biological mechanisms of structure. Here, we present a strategy for structural motif-based dynamic alignment, namely, RNA secondary-structural motif-comparing (RNAsmc), to identify structural motifs and quantitatively evaluate their underlying molecular functions. RNAsmc also has strong robustness to sequence length, folding protocol and RNA structural profile by chemical probing. Notably, it is also applicable to quantify structural variation in special RNA editing events (SNVs or SNPs, fragment insertion or deletion, etc.). The findings indicate that RNAsmc can uncover the heterogeneity of RNA secondary structure and score for similarities among components, which provides an impetus to cluster RNA families and evaluate allosteric effects. We find that RNAsmc exhibits remarkable detection efficiency for experimentally-derived RiboSNitches. Finally, the pipeline was assembled into an R software package to serve as an automated toolkit to explore, align, and cluster RSS. It is freely available for download at https://CRAN.R-project.org/package=RNAsmc.

Expression profiles and functional prediction of long non-coding RNAs LINC01133, ZEB1-AS1 and ABHD11-AS1 in the luminal subtype of breast cancer

Background

Luminal breast cancer (BC) is the most frequent subtype accounting for more than 70% of BC. LncRNAs, a class of non-coding RNAs with more than 200 nucleotides, are involved in a variety of cellular processes and biological functions. Abberant expression is related to the development of various cancers, such as breast cancer. LINC01133, ZEB1-AS1, and ABHD11-AS1 were reported to be dysregulated in different cancers. However, their expression level in luminal BC remains poorly known. The aim of the present study was to evaluate the potential roles of these lncRNAs in BC, especially in luminal subtypes.

Methods

A comprehensive analysis was performed using the Lnc2Cancer database to identify novel cancer-associated lncRNA candidates. After conducting a literature review, three novel lncRNAs named LINC01133, ZEB1-AS1, and ABHD11-AS1 were chosen as target genes of the present study. Quantitative real‐time polymerase chain reaction (qRT-PCR) was used to evaluate the expression level of the mentioned lncRNAs in both luminal BC tissues and cell lines. Then, the correlation of the three mentioned lncRNAs expression with clinicopathological characteristics of the patients was studied. Moreover, several datasets were used to discover the potential roles and functions of LINC01133, ZEB1-AS1 and ABHD11-AS1 in luminal subtype of BC.

Results

According to the qRT-PCR assay, the expression levels of LINC01133 and ZEB1-AS1 were decreased in luminal BC tissues and cell lines. On the other hand, ABHD11-AS1 was upregulated in the above-mentioned samples. The expression levels of LINC01133, ZEB1-AS1, and ABHD11-AS1 were not associated with any of the clinical features. Also, the results obtained from the bioinformatics analyses were consistent with qRT-PCR data. Functional annotation of the co-expressed genes with the target lncRNAs, protein–protein interactions and significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways across luminal BC were also obtained using bioinformatics analysis.

Conclusions

Taken together, our findings disclosed the dysregulation of LINC01133, ZEB1-AS1, and ABHD11-AS1 in luminal BC. It was revealed that LINC01133 and ZEB1-AS1 expression was significantly downregulated in luminal BC tissues and cell lines, while ABHD11-AS1 was upregulated considerably in the mentioned tissues and cell lines. Also, bioinformatics and systems biology analyses have helped to identify the possible role of these lncRNAs in luminal BC. However, further analysis is needed to confirm the current findings.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03026-7.

Evaluation of common genetic variants in vitamin E-related pathway genes and colorectal cancer susceptibility

Vitamin E is effective for preventing the risk of cancer. However, few studies have elucidated the mechanism of vitamin E in cancer occurrence. Herein, we aimed to identify the genetic variants in vitamin E-related pathway genes associated with colorectal cancer risk. We applied logistic regression models to assess the association between single-nucleotide polymorphisms (SNPs) in vitamin E-related pathway genes and colorectal cancer risk in the Chinese and European population. The false discovery rate (FDR) method was used to correct multiple comparisons. The mRNA and protein expression analysis were evaluated in public database and in-house RNA-Seq data. SCARB1 rs73227586 was identified significantly increased risk of colorectal cancer in the Chinese population (odd ratio (OR) = 1.46, 95% confidence interval (CI) = 1.22-1.75, P = 2.99 × 10^-5). This finding was further validated in the European population (OR = 1.11, 95% CI = 1.02-1.20, P = 1.44 × 10^-2). Additionally, the mRNA and protein expression of SCARB1 were markedly up-regulated in colorectal tumor tissues. Moreover, rs73227586 T allele could increase the minimum free energy (MFE) and weaken binding ability to transcription factor ELL2. Our findings indicated that SCARB1 may play a carcinogenic role in colorectal cancer. Genetic variants in vitamin E-related pathway genes may concern to be predictors of colorectal cancer risk.

Cell-Free DNA Detection of Tumor Mutations in Heterogeneous, Localized Prostate Cancer Via Targeted, Multiregion Sequencing

Cell-free DNA (cfDNA) may allow for minimally invasive identification of biologically relevant genomic alterations and genetically distinct tumor subclones. Although existing biomarkers may detect localized prostate cancer, additional strategies interrogating genomic heterogeneity are necessary for identifying and monitoring aggressive disease. In this study, we aimed to evaluate whether circulating tumor DNA can detect genomic alterations present in multiple regions of localized prostate tumor tissue.

METHODS

Low-pass whole-genome and targeted sequencing with a machine-learning guided 2.5-Mb targeted panel were used to identify single nucleotide variants, small insertions and deletions (indels), and copy-number alterations in cfDNA. The majority of this study focuses on the subset of 21 patients with localized disease, although 45 total individuals were evaluated, including 15 healthy controls and nine men with metastatic castration-resistant prostate cancer. Plasma cfDNA was barcoded with duplex unique molecular identifiers. For localized cases, matched tumor tissue was collected from multiple regions (one to nine samples per patient) for comparison.

RESULTS

Somatic tumor variants present in heterogeneous tumor foci from patients with localized disease were detected in cfDNA, and cfDNA mutational burden was found to track with disease severity. Somatic tissue alterations were identified in cfDNA, including nonsynonymous variants in FOXA1, PTEN, MED12, and ATM. Detection of these overlapping variants was associated with seminal vesicle invasion (P = .019) and with the number of variants initially found in the matched tumor tissue samples (P = .0005).

CONCLUSION

Our findings demonstrate the potential of targeted cfDNA sequencing to detect somatic tissue alterations in heterogeneous, localized prostate cancer, especially in a setting where matched tumor tissue may be unavailable (ie, active surveillance or treatment monitoring).

Dynamic Molecular Epidemiology Reveals Lineage-Associated Single-Nucleotide Variants That Alter RNA Structure in Chikungunya Virus

Chikungunya virus (CHIKV) is an emerging Alphavirus which causes millions of human infections every year. Outbreaks have been reported in Africa and Asia since the early 1950s, from three CHIKV lineages: West African, East Central South African, and Asian Urban. As new outbreaks occurred in the Americas, individual strains from the known lineages have evolved, creating new monophyletic groups that generated novel geographic-based lineages. Building on a recently updated phylogeny of CHIKV, we report here the availability of an interactive CHIKV phylodynamics dataset, which is based on more than 900 publicly available CHIKV genomes. We provide an interactive view of CHIKV molecular epidemiology built on Nextstrain, a web-based visualization framework for real-time tracking of pathogen evolution. CHIKV molecular epidemiology reveals single nucleotide variants that change the stability and fold of locally stable RNA structures. We propose alternative RNA structure formation in different CHIKV lineages by predicting more than a dozen RNA elements that are subject to perturbation of the structure ensemble upon variation of a single nucleotide.

Transcriptome-wide high-throughput mapping of protein-RNA occupancy profiles using POP-seq

Interaction between proteins and RNA is critical for post-transcriptional regulatory processes. Existing high throughput methods based on crosslinking of the protein–RNA complexes and poly-A pull down are reported to contribute to biases and are not readily amenable for identifying interaction sites on non poly-A RNAs. We present Protein Occupancy Profile-Sequencing (POP-seq), a phase separation based method in three versions, one of which does not require crosslinking, thus providing unbiased protein occupancy profiles on whole cell transcriptome without the requirement of poly-A pulldown. Our study demonstrates that ~ 68% of the total POP-seq peaks exhibited an overlap with publicly available protein–RNA interaction profiles of 97 RNA binding proteins (RBPs) in K562 cells. We show that POP-seq variants consistently capture protein–RNA interaction sites across a broad range of genes including on transcripts encoding for transcription factors (TFs), RNA-Binding Proteins (RBPs) and long non-coding RNAs (lncRNAs). POP-seq identified peaks exhibited a significant enrichment (p value < 2.2e−16) for GWAS SNPs, phenotypic, clinically relevant germline as well as somatic variants reported in cancer genomes, suggesting the prevalence of uncharacterized genomic variation in protein occupied sites on RNA. We demonstrate that the abundance of POP-seq peaks increases with an increase in expression of lncRNAs, suggesting that highly expressed lncRNA are likely to act as sponges for RBPs, contributing to the rewiring of protein–RNA interaction network in cancer cells. Overall, our data supports POP-seq as a robust and cost-effective method that could be applied to primary tissues for mapping global protein occupancies.

Evaluation of the potential role of long non-coding RNA LINC00961 in luminal breast cancer: a case-control and systems biology study

Background

Luminal subtype is the most common subgroup of breast cancer (BC), accounting for more than 70% of this cancer. Long non-coding RNAs (lncRNAs) are a group of RNAs which play critical roles in diverse cellular processes. It is proved that dysregulation of them can contribute to the development of various cancers, including BC. LINC00961 was reported to be downregulated in several cancers, however, its expression level in BC remains largely unknown. The purpose of the present study was to investigate the possible role of LINC00961 in luminal A and B subtypes of BC.

Methods

To obtain novel lncRNAs associated with different cancers and differentially expressed lncRNAs (DElncRNAs) between BC tumor and normal tissues, Lnc2Cancer and GDC databases were used, respectively. After performing literature review, the expression level of the selected lncRNA (LINC00961) was evaluated in 79 luminal A and B BC specimens and adjacent non-cancerous tissues by Quantitative Reverse Transcription PCR (qRT-PCR). LINC00961 expression was also evaluated in two luminal A BC cell lines, compared to a normal breast cell line. The comparison of the differences between tumor and adjacent non-tumor samples was performed by paired sample t-test. Moreover, correlation analysis between LINC00961 expression and clinicopathological features was performed using the chi-square, fisher exact, and independent t-test. In order to investigate the possible roles of LINC00961 in luminal A and B BC, different bioinformatics analyses such as functional annotation of the LINC00961 co-expressed genes and protein–protein interaction (PPI) networks construction were also performed.

Results

LINC00961 was selected as a significant DElncRNA which had not been studied in BC. According to q-RT PCR assay, LINC00961 was downregulated in luminal BC tissues and cell lines. Its expression was correlated with smoking status and the age of menarche in luminal BC patients. Also, the results of the bioinformatics analysis were consistent with the data obtained from q-RT PCR assay. The final results indicated that LINC00961 might be involved in multiple cancer-associated pathways such as chemokine, Ras and PI3K–Akt signaling pathways, GPCR ligand binding, and signal transduction in luminal subtypes of BC. CDH5, GNG11, GNG8, SELL, S1PR1, CCL19, FYN, ACAN, CD3E, ACVRL1, CAV1, and PPARGC1A were identified as the top hub genes of the PPI networks across luminal subgroup.

Conclusion

Our findings suggested that LINC00961 was significantly downregulated in luminal A and B subtypes of BC. Moreover, bioinformatics analysis provided a basis for better identification of the potential role of LINC00961 in luminal subtype of BC.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.