Landscape of Epstein-Barr virus gene expression and perturbations in cancer.

Epstein-Barr virus (EBV) is the causative agent for multiple neoplastic diseases of epithelial and lymphocytic origin1-3. The heterogeneity of the viral elements expressed and the mechanisms by which these coding and non-coding genes maintain cancer cell properties in vivo remain elusive4,5. Here we conducted a multi-modal transcriptomic analysis of EBV-associated neoplasms and identified that the ubiquitously expressed RPMS1 non-coding RNAs support cancer cell properties by disruption of the interferon response. Our map of EBV expression shows a variable, but pervasive expression of BNLF2 discerned from the overlapping LMP1 RNA in bulk sequencing data. Using long-read single-molecule sequencing, we identified three new viral elements within the RPMS1 gene. Furthermore, single-cell sequencing datasets allowed for the separation of cancer cells and healthy cells from the same tissue biopsy and the characterization of a microenvironment containing interferon gamma excreted by EBV-stimulated T-lymphocytes. In comparison with healthy epithelium, EBV-transformed cancer cells exhibited increased proliferation and inhibited immune response induced by the RPMS1-encoded microRNAs. Our atlas of EBV expression shows that the EBV-transformed cancer cells express high levels of non-coding RNAs originating from RPMS1 and that the oncogenic properties are maintained by RPMS1 microRNAs. Through bioinformatic disentanglement of single cells from cancer tissues we identified a positive feedback loop where EBV-activated immune cells stimulate cancer cells to proliferate, which in turn undergo viral reactivation and trigger an immune response.

Keratin 17 Impacts Global Gene Expression and Controls G2/M Cell Cycle Transition in Ionizing Radiation-Induced Skin Damage.

Keratin 17 (K17) is a cytoskeletal protein that is part of the intermediate filaments in epidermal keratinocytes. In K17-/- mice, ionizing radiation induced more severe hair follicle damage, whereas the epidermal inflammatory response was attenuated compared with that in wild-type mice. Both p53 and K17 have a major impact on global gene expression because over 70% of the differentially expressed genes in the skin of wild-type mice showed no expression change in p53-/- or K17-/- skin after ionizing radiation. K17 does not interfere with the dynamics of p53 activation; rather, global p53 binding in the genome is altered in K17-/- mice. The absence of K17 leads to aberrant cell cycle progression and mitotic catastrophe in epidermal keratinocytes, which is due to nuclear retention, thus reducing the degradation of B-Myb, a key regulator of the G2/M cell cycle transition. These results expand our understanding of the role of K17 in regulating global gene expression and ionizing radiation-induced skin damage.

Chemotherapy suppresses SHH gene expression via a specific enhancer.

Sonic hedgehog (SHH) signaling is a key regulator of embryonic development and tissue homeostasis that is involved in gastrointestinal (GI) cancer progression. Regulation of SHH gene expression is a paradigm of long-range enhancer function. Using the classical chemotherapy drug 5-fluorouracil (5FU) as an example, here we show that SHH gene expression is suppressed by chemotherapy. SHH is downstream of immediate early genes (IEGs), including Early growth response 1 (Egr1). A specific 139 kb upstream enhancer is responsible for its down-regulation. Knocking down EGR1 expression or blocking its binding to this enhancer renders SHH unresponsive to chemotherapy. We further demonstrate that down-regulation of SHH expression does not depend on 5FU's impact on nucleotide metabolism or DNA damage; rather, a sustained oxidative stress response mediates this rapid suppression. This enhancer is present in a wide range of tumors and normal tissues, thus providing a target for cancer chemotherapy and its adverse effects on normal tissues. We propose that SHH is a stress-responsive gene downstream of IEGs, and that traditional chemotherapy targets a specific enhancer to suppress its expression.

Development and Validation of Nine-RNA Binding Protein Signature Predicting Overall Survival for Kidney Renal Clear Cell Carcinoma.

Cumulative studies have shown that RNA binding proteins (RBPs) play an important role in numerous malignant tumors and are related to the occurrence and progression of tumors. However, the role of RBPs in kidney renal clear cell carcinoma (KIRC) is not fully understood. In this study, we first downloaded gene expression data and corresponding clinical information of KIRC from the Cancer Genome Atlas (TCGA) database, International Cancer Genome Consortium (ICGC), and Gene Expression Omnibus (GEO) database, respectively. A total of 137 differentially expressed RBPs (DERBPs) were then identified between normal and tumor tissue, including 38 downregulated and 99 upregulated RBPs. Nine RBPs (EIF4A1, RPL36A, EXOSC5, RPL28, RPL13, RPS19, RPS2, EEF1A2, and OASL) were served as prognostic genes and exploited to construct a prognostic model through survival analysis. Kaplan-Meier curves analysis showed that the low-risk group had a better survival outcome when compared with the high-risk group. The area under the curve (AUC) value of the prognostic model was 0.713 in the TCGA data set (training data set), 0.706 in the ICGC data set, and 0.687 in the GSE29609 data set, respectively, confirming a good prognostic model. The prognostic model was also identified as an independent prognostic factor for KIRC survival by performing cox regression analysis. In addition, we also built a nomogram relying on age and the prognostic model and internal validation in the TCGA data set. The clinical benefit of the prognostic model was revealed by decision curve analysis (DCA). Gene set enrichment analysis revealed several crucial pathways (ERBB signaling pathway, pathways in cancer, MTOR signaling pathway, WNT signaling pathway, and TGF BETA signaling pathway) that may explain the underlying mechanisms of KIRC. Furthermore, potential drugs for KIRC treatment were predicted by the Connectivity Map (Cmap) database based on DERBPs, including several important drugs, such as depudecin and vorinostat, that could reverse KIRC gene expression, which may provide reference for the treatment of KIRC. In summary, we developed and validated a robust nine-RBP signature for KIRC prognosis prediction. A nomogram with risk score and age can be applied to promote the individualized prediction of overall survival in patients with KIRC. Moreover, the two drugs depudecin and vorinostat may contribute to KIRC treatment.

Identification of 12 immune-related lncRNAs and molecular subtypes for the clear cell renal cell carcinoma based on RNA sequencing data.

Clear cell renal cell carcinoma (ccRCC) is the most common type of renal cell carcinoma (RCC). Despite the existing extensive research, the molecular and pathogenic mechanisms of ccRCC are elusive. We aimed to identify the immune-related lncRNA signature and molecular subtypes associated with ccRCC. By integrating 4 microarray datasets from Gene Expression Omnibus database, we identified 49 immune-related genes. The corresponding immune-related lncRNAs were further identified in the TCGA dataset. 12-lncRNAs prognostic and independent signature was identified through survival analysis and survival difference between risk groups was further identified based on the risk score. Besides, we identified 3 molecular subtypes and survival analysis result showed that cluster 2 has a better survival outcome. Further, ssGSEA enrichment analysis for the immune-associated gene sets revealed that cluster 1 corresponded to a high immune infiltration level. While cluster 2 and cluster 3 corresponded to low and medium immune infiltration level, respectively. In addition, we validated the 12-lncRNA prognostic signature and molecular subtypes in an external validation dataset from the ICGC database. In summary, we identified a 12-lncRNA prognostic signature which may provide new insights into the molecular mechanisms of ccRCC and the molecular subtypes provided a theoretical basis for personalized treatment by clinicians.

Coupling of apical-basal polarity and planar cell polarity to interpret the Wnt signaling gradient in feather development.

Sensing a global directional cue to orient cell growth is crucial in tissue morphogenesis. An anterior-posterior gradient of Wnt signaling controls the helical growth of feather branches (barbs), and thus the formation of bilateral feathers. However, it remains unclear how the keratinocytes sense this gradient and orient barb growth. Here, we show that in chicken, owing to feather branching, the global Wnt gradient is subdivided into periodic barbs. Within each barb, the anterior barbule plate cells tilt before the posterior cells. The core planar cell polarity gene Prickle1 is involved, as knockdown of its expression resulted in no cell shape change and no barb tilting. Furthermore, perturbation of the Wnt gradient leads to diffusive Prickle1 expression and loss of barb orientation. Finally, the asymmetric distribution of Wnt6/Fzd10 is coordinated by the apical-basal polarity of the barbule plate keratinocytes, which is in turn regulated by the Par3/aPKC machinery. Our data elucidate a new mechanism through which the global Wnt signaling gradient is interpreted locally to construct complex spatial forms.