A novel LncRNA-based prognostic score reveals TP53-dependent subtype of lung adenocarcinoma with poor survival

The prognostic signatures play an essential role in the era of personalised therapy for cancer patients including lung adenocarcinoma (LUAD). Long noncoding RNA (LncRNA), a relatively novel class of RNA, has shown to play a crucial role in all the areas of cancer biology. Here, we developed and validated a robust LncRNA-based prognostic signature for LUAD patients using three different cohorts. In the discovery cohort, four LncRNAs were identified with 10% false discovery rate and a hazard ratio of >10 using univariate Cox regression analysis. A risk score, generated from the four LncRNAs’ expression, was found to be a significant predictor of survival in the discovery and validation cohort (p = 9.97 × 10 ⁻⁸ and 1.41 × 10 ⁻³, respectively). Further optimisation of four LncRNAs signature in the validation cohort, generated a three LncRNAs prognostic score (LPS), which was found to be an independent predictor of survival in both the cohorts ( p = 1.00 × 10 ⁻⁶ and 7.27 × 10 ⁻⁴, respectively). The LPS also significantly divided survival in clinically important subsets, including Stage I ( p = 9.00 × 10 ⁻⁴ and 4.40 × 10 ⁻², respectively), KRAS wild-type (WT), KRAS mutant ( p = 4.00 × 10 ⁻³ and 4.30 × 10 ⁻², respectively) and EGFR WT ( p = 2.00 × 10 ⁻⁴). In multivariate analysis LPS outperformed, eight previous prognosticators. Further, individual members of LPS showed a significant correlation with survival in microarray data sets. Mutation analysis showed that high-LPS patients have a higher mutation rate and inactivation of the TP53 pathway. In summary, we identified and validated a novel LncRNA signature LPS for LUAD.     

Systematic construction and validation of an immune prognostic model for lung adenocarcinoma

Lung adenocarcinoma (LUAD), the most common non‐small‐cell lung cancer, is characterized by a dense lymphocytic infiltrate, which indicates that the immune system plays an active role in the development and growth of this cancer. However, no investigations to date have proposed robust models for predicting survival outcome for patients with LUAD in terms of tumour immunology. A total of 761 LUAD patients were included in this study, in which the database of The Cancer Genome Atlas (TCGA) was utilized for discovery, and the Gene Expression Omnibus (GEO) database was utilized for validation. Bioinformatics analysis and R language tools were utilized to construct an immune prognostic model and annotate biological functions. Lung adenocarcinoma showed a weakened immune phenotype compared with adjacent normal tissues. Immune‐related gene sets were profiled, an immune prognostic model based on 2 immune genes (ANLN and F2) was developed with the TCGA database to distinguish cases as having a low or high risk of unfavourable prognosis, and the model was verified with the GEO database. The model was prognostically significant in stratified cohorts, including stage I‐II, stage III‐IV and epidermal growth factor receptor (EGFR) mutant subsets, and was considered to be an independent prognostic factor for LUAD. Furthermore, the low‐ and high‐risk groups showed marked differences in tumour‐infiltrating leucocytes, tumour mutation burden, aneuploidy and PD‐L1 expression. In conclusion, an immune prognostic model was proposed for LUAD that is capable of independently identifying patients at high risk for poor survival, suggesting a relationship between local immune status and prognosis.     

Identification of necroptosis‐related gene signature and characterization of tumour microenvironment infiltration in non‐small‐cell lung cancer

Necroptosis is a programmed necrosis in a caspase‐independent fashion. The role of necroptosis‐related genes (NRGs) in lung cancer remains unknow. Herein, we classified TCGA‐LUAD cohort into two necroptosis‐related subtypes (C1 and C2) by consensus clustering analysis. The result showed that subtype C1 had a favourable prognosis and higher infiltration levels of immune cells. Moreover, subtype C1 was more activated in immune‐associated pathways. Then, we established an NRG prognosis model (NRG score) composed of six NRGs (RIPK3, MLKL, TLR2, TLR4, TNFRSF1A, NDRG2) and divided the cohort into low‐ and high‐risk group. We found that the NRG score was associated with prognosis, tumour immune microenvironment and tumour mutation burden. We also constructed an accurate nomogram model to improve the clinical applicability of NRG score. The result indicated that NRG score may be an independent prognostic marker for lung cancer patients. Taken together, we established a prognosis model that may deepen the understanding of NRGs in lung cancer and provide a basis for developing more effective immunotherapy strategies.     

Population‐specific responses of an insect herbivore to variation in host‐plant quality

Anthropogenic climate change poses a substantial challenge to many organisms, to which they need to respond to avoid fitness reductions. Investigating responses to environmental change is particularly interesting in herbivores, as they are potentially affected by indirect effects mediated via variation in host‐plant quality. We here use the herbivorous insect Pieris napi to investigate geographic variation in the response to variation in food quality. We performed a common garden experiment using replicated populations from Germany and Italy, and manipulated host quality by growing host plants at different temperature and water regimes. We found that feeding on plants grown at a higher temperature generally diminished the performance of P. napi, evidenced by a prolonged development time and reduced larval growth rate, body mass, fat content, and phenoloxidase activity. Genotype by environment interactions (G × E) were present in several performance traits, indicating that Italian populations (1) respond more strongly to variation in host‐plant quality and (2) are more sensitive to poor food quality than German ones. This may reflect a cost of the rapid lifestyle found in Italian populations. Consequently, German populations may be more resilient against environmental perturbations and may perhaps even benefit from warmer temperatures, while Italian populations will likely suffer from the concomitantly reduced host‐plant quality. Our study thus exemplifies how investigating G × E may help to better understand the vulnerability of populations to climate change.     

Quantitative restoration of immune defense in old animals determined by naive antigen‐specific CD8 T‐cell numbers

Older humans and animals often exhibit reduced immune responses to infection and vaccination, and this often directly correlates to the numbers and frequency of naive T (Tn) cells. We found such a correlation between reduced numbers of blood CD8+ Tn cells and severe clinical outcomes of West Nile virus (WNV) in both humans naturally exposed to, and mice experimentally infected with, WNV. To examine possible causality, we sought to increase the number of CD8 Tn cells by treating C57BL/6 mice with IL‐7 complexes (IL‐7C, anti‐IL‐7 mAb bound to IL‐7), shown previously to efficiently increase peripheral T‐cell numbers by homeostatic proliferation. T cells underwent robust expansion following IL‐7C administration to old mice increasing the number of total T cells (>fourfold) and NS4b:H‐2D^b‐restricted antigen‐specific CD8 T cells (twofold). This improved the numbers of NS4b‐specific CD8 T cells detected at the peak of the response against WNV, but not survival of WNV challenge. IL‐7C‐treated old animals also showed no improvement in WNV‐specific effector immunity (neutralizing antibody and in vivo T‐cell cytotoxicity). To test quantitative limits to which CD8 Tn cell restoration could improve protective immunity, we transferred graded doses of Ag‐specific precursors into old mice and showed that injection of 5400 (but not of 1800 or 600) adult naive WNV‐specific CD8 T cells significantly increased survival after WNV. These results set quantitative limits to the level of Tn reconstitution necessary to improve immune defense in older organisms and are discussed in light of targets of immune reconstitution.     

癌相关通路在肺腺癌中的共扰动机制

肺腺癌的发生涉及多个生物学功能通路的扰动,其遗传改变频繁地发生于MAPK信号、p53信号、Wnt信号、细胞周期和mTOR等通路的基因中。解析癌相关通路间的共扰动机制对我们理解癌机制以及寻找诊断标记具有重要意义。因此,本文基于肺腺癌突变谱数据,研究上述癌相关通路在肺腺癌中的共扰动机制。结果发现:在肺腺癌发生的过程中,MAPK信号、p53信号、Wnt信号、细胞周期和mTOR等通路同时被扰动。在不同的癌样本中,一对通路可能通过以下三种方式被共同扰动:(1)在两条通路中的不同基因间的共突变;(2)两条通路相互交叠基因的突变;(3)与两条通路同时具有频繁的互作关系的蛋白质的编码基因的突变。该结果提示,癌相关通路对在不同的样本中可能通过不同的方式被共扰动,这也可能是造成癌症异质性的重要原因之一。     

Characterizing heterogeneity of non‐small cell lung tumour microenvironment to identify signature prognostic genes

Growing evidence has highlighted the immune response as an important feature of carcinogenesis and therapeutic efficacy in non‐small cell lung cancer (NSCLC). This study focused on the characterization of immune infiltration profiling in patients with NSCLC and its correlation with survival outcome. All TCGA samples were divided into three heterogeneous clusters based on immune cell profiles: cluster 1 (''low infiltration'' cluster), cluster 2 (''heterogeneous infiltration'' cluster) and cluster 3 (''high infiltration'' cluster). The immune cells were responsible for a significantly favourable prognosis for the ''high infiltration'' community. Cluster 1 had the lowest cytotoxic activity, tumour‐infiltrating lymphocytes and interferon‐gamma (IFN‐γ), as well as immune checkpoint molecules expressions. In addition, MHC‐I and immune co‐stimulator were also found to have lower cluster 1 expressions, indicating a possible immune escape mechanism. A total of 43 differentially expressed genes (DEGs) that overlapped among the groups were determined based on three clusters. Finally, based on a univariate Cox regression model, prognostic immune‐related genes were identified and combined to construct a risk score model able to predict overall survival (OS) rates in the validation datasets.     

miR‐126‐5p enhances radiosensitivity of lung adenocarcinoma cells by inhibiting EZH2 via the KLF2/BIRC axis

Radiotherapy is a common method for the treatment of lung adenocarcinoma, but it often fails due to the relative non‐susceptibility of lung adenocarcinoma cells to radiation. We aimed to discuss the related mechanisms by which miR‐126‐5p might mediate radiosensitivity of lung adenocarcinoma cells. The binding affinity between miR‐126‐5p and EZH2 and between KLF2 and BIRC5 was identified using multiple assays. A549 and H1650 cells treated with X‐ray were transfected with miR‐126‐5p mimic/inhibitor, oe‐EZH2, or si‐KLF2 to detect cell biological functions and radiosensitivity. Finally, lung adenocarcinoma nude mouse models were established. miR‐126‐5p and KLF2 were poorly expressed, while EZH2 and BIRC5 were upregulated in lung adenocarcinoma tissues and cells. miR‐126‐5p targeted EZH2 to promote the KLF2 expression so as to inhibit BIRC5 activation. Both in vitro and in vivo experiments verified that elevated miR‐126‐5p inhibited cell migration and promoted apoptosis to enhance the sensitivity of lung adenocarcinoma cells to radiotherapy via the EZH2/KLF2/BIRC5 axis. Collectively, miR‐126‐5p downregulated EZH2 to facilitate the sensitivity of lung adenocarcinoma cells to radiotherapy via KLF2/BIRC5.     

Genetic reduction of mTOR extends lifespan in a mouse model of Hutchinson‐Gilford Progeria syndrome

Hutchinson‐Gilford progeria syndrome (HGPS) is a rare accelerated aging disorder most notably characterized by cardiovascular disease and premature death from myocardial infarction or stroke. The majority of cases are caused by a de novo single nucleotide mutation in the LMNA gene that activates a cryptic splice donor site, resulting in production of a toxic form of lamin A with a 50 amino acid internal deletion, termed progerin. We previously reported the generation of a transgenic murine model of progeria carrying a human BAC harboring the common mutation, G608G, which in the single‐copy state develops features of HGPS that are limited to the vascular system. Here, we report the phenotype of mice bred to carry two copies of the BAC, which more completely recapitulate the phenotypic features of HGPS in skin, adipose, skeletal, and vascular tissues. We further show that genetic reduction of the mechanistic target of rapamycin (mTOR) significantly extends lifespan in these mice, providing a rationale for pharmacologic inhibition of the mTOR pathway in the treatment of HGPS.     

MTHFD2 promotes tumorigenesis and metastasis in lung adenocarcinoma by regulating AKT/GSK‐3β/β‐catenin signalling

Recent studies have demonstrated that one‐carbon metabolism plays a significant role in cancer development. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2), a mitochondrial enzyme of one‐carbon metabolism, has been reported to be dysregulated in many cancers. However, the specific role and mechanism of MTHFD2 in lung adenocarcinoma (LUAD) still remains unclear. In this study, we evaluated the clinicopathological and prognostic values of MTHFD2 in LUAD patients. We conducted a series of functional experiments in vivo and in vitro to explore novel mechanism of MTHFD2 in LUAD. The results showed that MTHFD2 was significantly up‐regulated in LUAD tissues and predicted poor prognosis of LUAD patients. Knockdown of MTHFD2 dramatically inhibited cell proliferation and migration by blocking the cell cycle and inducing the epithelial‐mesenchymal transition (EMT). In addition, MTHFD2 knockdown suppressed LUAD growth and metastasis in cell‐derived xenografts. Mechanically, we found that MTHFD2 promoted LUAD cell growth and metastasis via AKT/GSK‐3β/β‐catenin signalling. Finally, we identified miR‐30a‐3p as a novel regulator of MTHFD2 in LUAD. Collectively, MTHFD2 plays an oncogenic role in LUAD progression and is a promising target for LUAD diagnosis and therapy.     

Characterization of novel nonacid glycosphingolipids as biomarkers of human gastric adenocarcinoma

Changes in glycosphingolipid structures have been shown to occur during the development of several types of human cancers, generating cancer-specific carbohydrate structures that could be used as biomarkers for diagnosis and therapeutic targeting. In this study, we characterized nonacid glycosphingolipids isolated from a human gastric adenocarcinoma by mass spectrometry, enzymatic hydrolysis, and by binding with a battery of carbohydrate-recognizing ligands. We show that the majority of the complex nonacid glycosphingolipids had type 2 (Galβ4GlcNAc) core chains (neolactotetraosylceramide, the Le^x, H type 2, x₂, and the P1 pentaosylceramides, and the Le^y, A type 2, and neolacto hexaosylceramides). We also found glycosphingolipids with type 1 (Galβ3GlcNAc) core (lactotetraosylceramide and the H type 1 pentaosylceramide) and globo (GalαGal) core chains (globotriaosylceramide and globotetraosylceramide). Interestingly, we characterized two complex glycosphingolipids as a P1 heptaosylceramide (Galα4Galβ4GlcNAcβ3Galβ4GlcNAcβ3Gal β4Glcβ1Cer) and a branched P1 decaosylceramide (Galα4Gal β4GlcNAcβ3(Galα4Galβ4GlcNAcβ6)Galβ4GlcNAcβ3Galβ4Glc β1Cer). These are novel glycosphingolipid structures and the first reported cases of complex glycosphingolipids larger than pentaosylceramide carrying the P1 trisaccharide. We propose that these P1 glycosphingolipids may represent potential biomarkers for the early diagnosis of gastric cancer.     

ANKRD49 promotes the invasion and metastasis of lung adenocarcinoma via a P38/ATF‐2 signalling pathway

Lung adenocarcinoma (LUAD) is the most challenging neoplasm to treat in clinical practice. Ankyrin repeat domain 49 protein (ANKRD49) is highly expressed in several carcinomas; however, its pattern of expression and role in LUAD are not known. Tissue microarrays, immunohistochemistry, χ ² test, Spearman correlation analysis, Kaplan–Meier, log‐rank test, and Cox''s proportional hazard model were used to analyse the clinical cases. The effect of ANKRD49 on the LUAD was investigated using CCK‐8, clonal formation, would healing, transwell assays, and nude mice experiment. Expressions of ANKRD49 and its associated downstream protein molecules were verified by real‐time PCR, Western blot, immunohistochemistry, and/or immunofluorescence analyses. ANKRD49 expression was highly elevated in LUAD. The survival rate and Cox''s modelling analysis indicated that there may be an independent prognostic indicator for LUAD patients. We also found that ANKRD49 promoted the invasion and migration in both in in vitro and in vivo assays, through upregulating matrix metalloproteinase (MMP)‐2 and MMP‐9 activities via the P38/ATF‐2 signalling pathway Our findings suggest that ANKRD49 is a latent biomarker for evaluating LUAD prognosis and promotes the metastasis of A549 cells via upregulation of MMP‐2 and MMP‐9 in a P38/ATF‐2 pathway‐dependent manner.     

The diagnostic and prognostic value of pseudogene SIGLEC17P in lung adenocarcinoma and a preliminary functional study

Among malignant tumors, lung adenocarcinoma (LUAD) is the leading cause of death worldwide. This study explored the diagnostic, prognostic value, and preliminary functional verification of sialic acid binding Ig like lectin 17, pseudogene (SIGLEC17P) in LUAD. Prognostic lncRNAs for LUAD were identified by The Cancer Genome Atlas and quantitative real-time PCR (qRT-PCR) was used to detect the expression of SIGLEC17P in LUAD and paracarcinoma tissues. Subsequently, lentiviral vectors were used to overexpress SIGLEC17P in A549 and H1299 cells. The effects of SIGLEC17P overexpression on the proliferation, migration, and invasiveness of LUAD cells (A549 and H1299) were evaluated by Cell Counting Kit-8, wound healing, and transwell migration assays, respectively. Bioinformatics analyses were performed to reveal the potential pathways in which SIGLEC17P is involved in LUAD. qRT-PCR results revealed low SIGLEC17P expression in LUAD tissues and a significant association with the N stage, T stage, and tumor node metastasis stage. Furthermore, the receiver operating characteristic curve demonstrated a reliable diagnostic value. The proliferation, migration, and invasion of LUAD cells were inhibited by overexpression of SIGLEC17P. Bioinformatics analyses suggested that SIGLEC17P might exert antioncogenic effects in LUAD through the mir-20-3p/ADH1B or mir-4476-5p/DPYSL axis. In summary, our results revealed that SIGLEC17P acts as a prognostic biomarker, independent prognostic factor, and potential therapeutic target for patients with LUAD.     

Climate and habitat configuration limit range expansion and patterns of dispersal in a non‐native lizard

Invasive species are one of the main causes of biodiversity loss worldwide. As introduced, populations increase in abundance and geographical range, so does the potential for negative impacts on native communities. As such, there is a need to better understand the processes driving range expansion as species become established in recipient landscapes. Through an investigation into capacity for population growth and range expansion of introduced populations of a non‐native lizard (Podarcis muralis), we aimed to demonstrate how multi‐scale factors influence spatial spread, population growth, and invasion potential in introduced species. We collated location records of P. muralis presence in England, UK through data collected from field surveys and a citizen science campaign. We used these data as input for presence‐background models to predict areas of climate suitability at a national‐scale (5 km resolution), and fine‐scale habitat suitability at the local scale (2 m resolution). We then integrated local models into an individual‐based modeling platform to simulate population dynamics and forecast range expansion for 10 populations in heterogeneous landscapes. National‐scale models indicated climate suitability has restricted the species to the southern parts of the UK, primarily by a latitudinal cline in overwintering conditions. Patterns of population growth and range expansion were related to differences in local landscape configuration and heterogeneity. Growth curves suggest populations could be in the early stages of exponential growth. However, annual rates of range expansion are predicted to be low (5–16 m). We conclude that extensive nationwide range expansion through secondary introduction is likely to be restricted by currently unsuitable climate beyond southern regions of the UK. However, exponential growth of local populations in habitats providing transport pathways is likely to increase opportunities for regional expansion. The broad habitat niche of P. muralis, coupled with configuration of habitat patches in the landscape, allows populations to increase locally with minimal dispersal.     

Development of a prognostic model based on an immunogenomic landscape analysis of medulloblastoma

Medulloblastoma (MB) is one of the most common central nervous system tumors in children. At present, the vital role of immune abnormalities has been proved in tumorigenesis and progression. However, the immune mechanism in MB is still poorly understood. In the present study, 51 differentially expressed immune-related genes (DE-IRGs) and 226 survival associated immune-related genes (Sur-IRGs) were screened by an integrated analysis of multi-array. Moreover, the potential pathways were enriched by functional analysis, such as ‘cytokine–cytokine receptor interaction’, ‘Ras signaling pathway’, ‘PI3K-Akt signaling pathway’ and ‘pathways in cancer’. Furthermore, 10 core IRGs were identified from DE-IRGs and Sur-IRGs. And the potential regulatory mechanisms of core IRGs were also explored. Additionally, a new prognostic model, including 7 genes (HDGF, CSK, PNOC, S100A13, RORB, FPR1, and ICAM2) based on IRGs, was established by multivariable COX analysis. In summary, our study revealed the underlying immune mechanism of MB. Moreover, we developed a prognostic model associated with clinical characteristics and could reflect the infiltration of immune cells.     

Robust prognostic model based on immune infiltration‐related genes and clinical information in ovarian cancer

Immune infiltration of ovarian cancer (OV) is a critical factor in determining patient''s prognosis. Using data from TCGA and GTEx database combined with WGCNA and ESTIMATE methods, 46 genes related to OV occurrence and immune infiltration were identified. Lasso and multivariate Cox regression were applied to define a prognostic score (IGCI score) based on 3 immune genes and 3 types of clinical information. The IGCI score has been verified by K‐M curves, ROC curves and C‐index on test set. In test set, IGCI score (C‐index = 0.630) is significantly better than AJCC stage (C‐index = 0.541, p < 0.05) and CIN25 (C‐index = 0.571, p < 0.05). In addition, we identified key mutations to analyse prognosis of patients and the process related to immunity. Chi‐squared tests revealed that 6 mutations are significantly (p < 0.05) related to immune infiltration: BRCA1, ZNF462, VWF, RBAK, RB1 and ADGRV1. According to mutation survival analysis, we found 5 key mutations significantly related to patient prognosis (p < 0.05): CSMD3, FLG2, HMCN1, TOP2A and TRRAP. RB1 and CSMD3 mutations had small p‐value (p < 0.1) in both chi‐squared tests and survival analysis. The drug sensitivity analysis of key mutation showed when RB1 mutation occurs, the efficacy of six anti‐tumour drugs has changed significantly (p < 0.05).     

Simvastatin inhibits stem cell proliferation in human leiomyoma via TGF‐β3 and Wnt/β‐Catenin pathways

Uterine leiomyoma (UL) is the most common gynaecologic tumour, affecting an estimated 70 to 80% of women. Leiomyomas develop from the transformation of myometrial stem cells into leiomyoma stem (or tumour‐initiating) cells. These cells undergo self‐renewal and differentiation to mature cells, both are necessary for the maintenance of tumour stem cell niche and tumour growth, respectively. Wnt/β‐catenin and TGF‐β/SMAD pathways, both overactive in UL, promote stem cell self‐renewal, crosstalk between stem and mature cells, cellular proliferation, extracellular matrix (ECM) accumulation and drive overall UL growth. Recent evidence suggests that simvastatin, an antihyperlipidemic drug, may have anti‐leiomyoma properties. Herein, we investigated the effects of simvastatin on UL stem cells. We isolated leiomyoma stem cells by flow cytometry using DyeCycle Violet staining and Stro‐1/CD44 surface markers. We found that simvastatin inhibits proliferation and induces apoptosis in UL stem cells. In addition, it also suppressed the expression of the stemness markers Nanog, Oct4 and Sox2. Simvastatin significantly decreased the production of the key ECM proteins, collagen 1 and fibronectin. Finally, it inhibited genes and/or proteins expression of TGF‐β1, 2 and 3, SMAD2, SMAD4, Wnt4, β‐Catenin, LRP6, AXIN2 and Cyclin D1 in UL stem cells, all are key drivers of the TGF‐β3/SMAD2 and Wnt4/β‐Catenin pathways. Thus, we have identified a novel stem cell‐targeting anti‐leiomyoma simvastatin effect. Further studies are needed to replicate these findings in vivo.