The construction and analysis of ceRNA network and patterns of immune infiltration in lung adenocarcinoma

Somatic Copy Number Alterations in Colorectal Cancer Lead to a Differentially Expressed ceRNA Network (ceRNet)

Colorectal cancer (CRC) represents the second deadliest malignancy worldwide. Around 75% of CRC patients exhibit high levels of chromosome instability that result in the accumulation of somatic copy number alterations. These alterations are associated with the amplification of oncogenes and deletion of tumor-ppressor genes and contribute to the tumoral phenotype in different malignancies. Even though this relationship is well known, much remains to be investigated regarding the effect of said alterations in long non-coding RNAs (lncRNAs) and, in turn, the impact these alterations have on the tumor phenotype. The present study aimed to evaluate the role of differentially expressed lncRNAs coded in regions with copy number alterations in colorectal cancer patient samples. We downloaded RNA-seq files of the Colorectal Adenocarcinoma Project from the The Cancer Genome Atlas (TCGA) repository (285 sequenced tumor tissues and 41 non-tumor tissues), evaluated differential expression, and mapped them over genome sequencing data with regions presenting copy number alterations. We obtained 78 differentially expressed (LFC > 1|< -1, padj < 0.05) lncRNAs, 410 miRNAs, and 5028 mRNAs and constructed a competing endogenous RNA (ceRNA) network, predicting significant lncRNA-miRNA-mRNA interactions. Said network consisted of 30 lncRNAs, 19 miRNAs, and 77 mRNAs. To understand the role that our ceRNA network played, we performed KEGG and GO analysis and found several oncogenic and anti-oncogenic processes enriched by the molecular players in our network. Finally, to evaluate the clinical relevance of the lncRNA expression, we performed survival analysis and found that C5orf64, HOTAIR, and RRN3P3 correlated with overall patient survival. Our results showed that lncRNAs coded in regions affected by SCNAs form a complex gene regulatory network in CCR.

Hypoxia-associated prognostic markers and competing endogenous RNA coexpression networks in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the most common form of non-small cell lung cancer (NSCLC). Hypoxia has been found in 50-60% of locally advanced solid tumors and is associated with poor prognosis in various tumors, including NSCLC. This study focused on hypoxia-associated molecular hallmarks in LUAD. Fifteen hypoxia-related genes were selected to define the hypoxia status of LUAD by ConsensusClusterPlus based on data from The Cancer Genome Atlas (TCGA). Then, we investigated the immune status under different hypoxia statuses. Subsequently, we constructed prognostic models based on hypoxia-related differentially expressed genes (DEGs), identified hypoxia-related microRNAs, lncRNAs and mRNAs, and built a network based on the competing endogenous RNA (ceRNA) theory. Two clusters (Cluster 1 and Cluster 2) were identified with different hypoxia statuses. Cluster 1 was defined as the hypoxia subgroup, in which all 15 hypoxia-associated genes were upregulated. The infiltration of CD4+ T cells and Tfh cells was lower, while the infiltration of regulatory T (Treg) cells, the expression of PD-1/PD-L1 and TMB scores were higher in Cluster 1, indicating an immunosuppressive status. Based on the DEGs, a risk signature containing 7 genes was established. Furthermore, three differentially expressed microRNAs (hsa-miR-9, hsa-miR-31, hsa-miR-196b) associated with prognosis under different hypoxia clusters and their related mRNAs and lncRNAs were identified, and a ceRNA network was built. This study showed that hypoxia was associated with poor prognosis in LUAD and explored the potential mechanism from the perspective of the gene signature and ceRNA theory.