Associations between genetic variants in mRNA splicing-related genes and risk of lung cancer: a pathway-based analysis from published GWASs

LSM2 is associated with a poor prognosis and promotes cell proliferation, migration, and invasion in skin cutaneous melanoma

BACKGROUND

Skin cutaneous melanoma (SKCM) is an extremely malignant tumor that is associated with a poor prognosis. LSM2 has been found to be related to different types of tumors; however, its role in SKCM is poorly defined. We aimed to determine the value of LSM2 as a prognostic biomarker for SKCM.

METHODS

The expression profile of LSM2 mRNA was compared between tumor and normal tissues in public databases, such as TCGA, GEO, and BioGPS. LSM2 protein expression was explored using immunohistochemistry (IHC) on a tissue microarray containing 44 SKCM tissues and 8 normal samples collected at our center. Kaplan-Meier analysis was performed to assess the prognostic value of LSM2 expression in patients with SKCM. SKCM cell lines with LSM2 knockdown were used to determine the effects of LSM2. Cell counting kit-8 (CCK8) and colony formation assays were conducted to assess cell proliferation, whereas wound healing and transwell assays were carried out to assess the migration and invasion abilities of SKCM cells.

RESULTS

LSM2 was more highly expressed at the mRNA and protein levels in SKCM than that in normal skin. Moreover, elevated expression of LSM2 was associated with shorter survival time and early recurrence in patients with SKCM. The in vitro results revealed that the silencing of LSM2 in SKCM cells significantly inhibited cell proliferation, migration, and invasion.

CONCLUSION

Overall, LSM2 contributes to malignant status and poor prognosis in patients with SKCM and may be identified as a novel prognostic biomarker and therapeutic target.

Identification of cancer risk assessment signature in patients with chronic obstructive pulmonary disease and exploration of the potential key genes

It is essential to assess the cancer risk for patients with chronic obstructive pulmonary disease (COPD). Comparing gene expression data from patients with lung cancer (a total of 506 samples) and those with cancer-adjacent normal lung tissues (a total of 370 samples), we generated a qualitative transcriptional signature consisting of 2046 gene pairs. The signature was verified in an evaluation dataset comprising 18 subjects with severe disease and 52 subjects with moderate disease (Wilcoxon rank-sum test; p = 7.33 × 10^-5). Similar results were obtained in other independent datasets. Among the gene pairs in the signature, 326 COPD stage-related gene pairs were identified based on Spearman's rank correlation tests and those gene pairs comprised 368 unique genes. Of these 368 genes, 16 genes were significantly dysregulated in COPD rat model data compared with control data. Some of these genes (Dhx16, Upf2, Notch3, Sec61a1, Dyrk2, and Hmmr) were altered when the COPD rat model was treated with traditional Chinese medicines (TCM), including Bufei Yishen formula, Bufei Jianpi formula, and Yiqi Zishen formula. Overall, the signature could predict the cancer incidence-risk of COPD and the identified key genes might provide guidance regarding both the treatment of COPD using TCM and the prevention of cancer in patients with COPD. KEY MESSAGESA cancer risk assessment signature was identified in patients with COPD.The signature is insensitive to batch effects and is well verified.COPD key genes identified in this study might play a crucial role in TCM treatment and cancer prevention.

Expression profile and prognostic values of LSM family in skin cutaneous melanoma

Background

The like-Smith (LSM) family is a group of RNA-binding proteins involved in RNA metabolism. However, their involvement in tumors, particularly skin cutaneous melanoma (SKCM), is not fully understood. In this study, we focused on the expression profiles and prognostic values of the LSM family in SKCM.

Methods

Raw data were downloaded from The Cancer Genome Atlas. The expression profile and prognostic value of LSM genes in SKCM were explored using the GEPIA, cBioPortal, and HPA databases. Protein–protein and gene–gene interaction analyses were performed using STRING and GeneMANIA. Enrichment and Cox regression analysis were conducted using R software. The TISIDB database was used to explore the relationship between LSMs and immunomodulators. Receiver operating characteristic curves and nomogram models were constructed to validate prognostic values.

Results

mRNA and protein expression levels of LSM2, LSM4, and LSM12 were significantly elevated in SKCM. The upregulated mRNA expression of LSM2 (p = 0.0013) and LSM4 (p = 0.0043) was significantly correlated with poor overall survival in patients with SKCM, whereas only LSM2 (p = 0.049) overexpression was markedly associated with worse disease-free survival. LSM2 overexpression was an independent risk factor (p = 0.013) and was confirmed to have a high prognostic value in SKCM using the receiver operating characteristic curve (AUC = 0.942) and nomogram models. All LSM genes were identified as genomic mutations, whereas alteration of LSM2 (p = 0.0153) significantly affected the overall survival in patients with SKCM. Significant correlations were observed between LSM family expression, immune cell infiltration, and immunomodulator. Furthermore, function and pathway enrichment analysis showed that the LSM family was mainly RNA binding proteins and involved in RNA splicing and degradation.

Conclusion

Expression profiles and prognostic values of LSM in SKCM were inconsistent. Among the LSM family, only LSM2 may serve as a potential poor prognosticator and immunotherapeutic target of SKCM.

Significance of Spliceosome-Related Genes in the Prediction of Prognosis and Treatment Strategies for Lung Adenocarcinoma

BACKGROUND

The leading cause of cancer-related fatalities globally is lung cancer; lung adenocarcinoma (LUAD) is the most common histological type in it. The spliceosome plays an important role in a majority of malignancies. However, it is yet unclear how spliceosome-related genes affect patients with LUAD in terms of treatment course and prognosis.

METHODS

Spliceosome-related genes were assessed from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database to obtain clinical information and gene expression in patients with LUAD. A spliceosome-related gene signature and prognostic model were constructed by using the least absolute shrinkage and selection operator (LASSO), time-dependent receiver operating characteristic (ROC), and nomogram. Immune infiltrate levels, mutation analysis, and pathway enrichment were predicted potential mechanisms of the signature by using single-sample gene set enrichment analysis (ssGSEA), Gene Set Cancer Analysis (GSCA) database, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO) database. Then, a protein-protein interaction (PPI) network and transcription factor- (TF-) hub gene and drug mining network were also established by Cytoscape software.

RESULTS

Firstly, we constructed a prognostic model for 11 spliceosome signature genes. Based on the prognostic risk score, we stratified patients with LUAD into high- and low-risk groups. The high- and low-risk groups were closely related to the OS, tumor immune infiltration level, immune checkpoint molecules, and tumor mutation burden (TMB) of LUAD patients. Based on PPI networks, we also predict relevant TF genes that may regulate signature prognostic genes. Finally, drugs including oxaliplatin, arsenic trioxide, cisplatin, and sunitinib were excavated for the treatment of the 11 spliceosome signature genes in LUAD patients.

CONCLUSION

In conclusion, this study is the first to explore the importance of spliceosome-related genes in the prognosis and treatment of LUAD. Through our study, we have innovatively provided potential prognosis genes and new therapeutic drug targets for the treatment of LUAD patients.

Genetic variants in splicing factor genes and susceptibility to bladder cancer

BACKGROUND

Genome-wide association studies have demonstrated that genetic variants are closely related to tumorigenesis and progression of cancer. However, the correlation between genetic variants in splicing factor genes and bladder cancer susceptibility remains unclear.

METHOD

A case-control study with 580 cases of bladder cancer and 1,101 controls was conducted to explore the association of single-nucleotide polymorphisms (SNPs) in splicing factors with bladder cancer susceptibility by logistic regression models, and multiple testing errors were justified by the false discovery rate (FDR) method. Next, we used the Cancer Genome Atlas (TCGA) database and the Gene Expression Omnibus (GEO) datasets to further analyze the differential expression of candidate genes.

RESULTS

We found that rs978416 G>A in RBFOX3 contributed to a reduced risk of bladder cancer [adjusted odds ratio (OR) = 0.72, 95% confidence internal (CI) = 0.62-0.84, P = 3.54 × 10^-5], especially in individuals who never smoked (P = 7.83 × 10^-5). Stratified analysis showed that the protective effect of rs978416 was more significant in the subgroup of low grade and non-muscle invasive bladder cancer. Furthermore, the RBFOX3 mRNA expression was decreased in bladder tumor tissues. However, the relatively high expression of RBFOX3 was related to a higher bladder cancer stage.

CONCLUSIONS

Our findings indicated that SNP rs978416 G>A in RBFOX3 may be related to bladder cancer predisposition in Chinese population and might serve as a novel biomarker for bladder cancer risk.

Potential Therapeutic and Prognostic Values of LSM Family Genes in Breast Cancer

In recent decades, breast cancer (BRCA) has become one of the most common diseases worldwide. Understanding crucial genes and their signaling pathways remain an enormous challenge in evaluating the prognosis and possible therapeutics. The "Like-Smith" (LSM) family is known as protein-coding genes, and its member play pivotal roles in the progression of several malignancies, although their roles in BRCA are less clear. To discover biological processes associated with LSM family genes in BRCA development, high-throughput techniques were applied to clarify expression levels of LSMs in The Cancer Genome Atlas (TCGA)-BRCA dataset, which was integrated with the cBioPortal database. Furthermore, we investigated prognostic values of LSM family genes in BCRA patients using the Kaplan-Meier database. Among genes of this family, LSM4 expression levels were highly associated with poor prognostic outcomes with a hazard ratio of 1.35 (95% confidence interval 1.21-1.51, p for trend = 3.4 × 10-7). MetaCore and GlueGo analyses were also conducted to examine transcript expression signatures of LSM family members and their coexpressed genes, together with their associated signaling pathways, such as "Cell cycle role of APC in cell cycle regulation" and "Immune response IL-15 signaling via MAPK and PI3K cascade" in BRCA. Results showed that LSM family members, specifically LSM4, were significantly correlated with oncogenesis in BRCA patients. In summary, our results suggested that LSM4 could be a prospective prognosticator of BRCA.

Differential alternative RNA splicing and transcription events between tumors from African American and White patients in The Cancer Genome Atlas

Individuals of African ancestry suffer disproportionally from higher incidence, aggressiveness, and mortality for particular cancers. This disparity likely results from an interplay among differences in multiple determinants of health, including differences in tumor biology. We used The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA aggregate expression datasets and identified differential alternative RNA splicing and transcription events (ARS/T) in cancers between self-identified African American (AA) and White (W) patients. We found that retained intron events were enriched among race-related ARS/T. In addition, on average, 12% of the most highly ranked race-related ARS/T overlapped between any two analyzed cancers. Moreover, the genes undergoing race-related ARS/T functioned in cancer-promoting pathways, and a number of race-related ARS/T were associated with patient survival. We built a web-application, CanSplice, to mine genomic datasets by self-identified race. The race-related targets have the potential to aid in the development of new biomarkers and therapeutics to mitigate cancer disparity.

Is Hydrogen Sulfide a Concern During Treatment of Lung Adenocarcinoma With Ammonium Tetrathiomolybdate?

Ammonium tetrathiomolybdate (ATTM) has been used in breast cancer therapy for copper chelation, as elevated copper promotes tumor growth. ATTM is also an identified H₂S donor and endogenous H₂S facilitates VitB₁₂-induced S-adenosylmethionine (SAM) generation, which have been confirmed in m⁶A methylation and lung cancer development. The m⁶A modification was recently shown to participate in lung adenocarcinoma (LUAD) progression. These conflicting analyses of ATTM's anticancer vs. H₂S's carcinogenesis suggest that H₂S should not be ignored during LUAD's treatment with ATTM. This study was aimed to explore ATTM's effects on LUAD cells and mechanisms associated with H₂S and m⁶A. It was found that treatment with ATTM inhibited cell growth at high concentrations, while enhanced cell growth at low concentrations in three LUAD cell lines (A549, HCC827, and PC9). However, another copper chelator triethylenetetramine, without H₂S releasing activity, was not found to induce cell growth. Low ATTM concentrations also elevated m⁶A content in A549 cells. Analysis of differentially expressed genes in TCGA cohort indicated that m⁶A writer METTL3 and reader YTHDF1 were upregulated while eraser FTO was downregulated in LUAD tissues, consistent with the findings of protein expression in patient tissues. ATTM treatment of A549 cells significantly increased METTL3/14 and YTHDF1 while decreased FTO expression. Furthermore, inhibition of m⁶A with shMETTL3 RNA significantly attenuated eukaryotic translation initiation factor (eIF) expressions in A549 cells. Correlation analysis indicated that small nuclear ribonucleic protein PRPF6 was positively expressed with YTHDF1 in LUAD tissues. Knockdown of YTHDF1 partially blocked both basal and ATTM-induced PRPF6 expression, as well as A549 cell growth. Lastly, ATTM treatment not only raised intracellular H₂S content but also upregulated H₂S-producing enzymes. Exogenous H₂S application mimicked ATTM's aforementioned effects, but the effects could be weakened by zinc-induced H₂S scavenging. Collectively, H₂S impedes ATTM-induced anticancer effects through YTHDF1-dependent PRPF6 m⁶A methylation in lung adenocarcinoma cells.

Associations of genetic variants in endocytic trafficking of epidermal growth factor receptor super pathway with risk of nonsyndromic cleft lip with or without cleft palate

BACKGROUND

The genetic etiology of nonsyndromic cleft lip with or without cleft palate (NSCL/P) has not been fully clarified to date. Epidermal growth factor receptor (EGFR) was reportedly involved in its biological establishment and regulation of cell migration during the embryonic stage.

METHODS

We selected a super pathway of endocytic trafficking of EGFR and investigated the associations of single-nucleotide polymorphisms (SNPs) in the super pathway with the risk of NSCL/P by analyzing our published genome-wide association study (GWAS) data from 504 NSCL/P individuals and 455 controls. After the false discovery rate (FDR) control, we conducted linkage disequilibrium (LD) analyses and conditional regression analyses to obtain independent lead SNPs. We performed LD analyses between the lead SNPs and the reported SNPs to find novel ones from our study. We annotated the lead SNPs and investigated their mapped genes in silico.

RESULTS

A total of 82 SNPs showed a statistical association with the risk of NSCL/P after FDR control. They contained three reported SNPs which were g.117068049G>A (rs7078160), g.117086783C>G (rs10886040), and g.117101266G>T (rs17095681). Four independent lead SNPs were obtained, including g.116979803 T>C (rs1905539) and g.117037960A>G (rs7902502) at 10q25.3, g.35720163G>C (rs75656820) at 17q12, and g.156864512G>A (rs1800877) at 1q23.1. Three of them were in low LD (r²  < 0.5) with the reported SNPs except g.117037960A>G (rs7902502), so these three were newly identified. Lead SNPs were mapped to three genes: SHTN1, AP2B1, and NTRK1. The three genes were relatively more highly expressed in the human craniofacial region and in the proximal maxillary location during the craniofacial development stage of the embryonic mouse.

CONCLUSION

Our results suggested that SHTN1, AP2B1, and NTRK1 might be associated with the development of NSCL/P.