Genome-wide analysis of survival in early-stage non-small-cell lung cancer

Knockdown of STK39 inhibits lung cancer brain metastasis by suppressing the CPSF4/NFκB/COX2 pathway

PURPOSE

Lung cancer is the most common cancer worldwide, and approximately 30% of lung cancer patients will develop brain metastases. Serine/threonine kinase 39 (STK39) plays a significant role in various malignancies. However, the role and mechanism of STK39 in lung cancer brain metastasis have not been reported.

METHODS

The expression levels of STK39 in lung cancer cells were detected using quantitative reverse transcription PCR (RT-qPCR) and Western blotting. STK39 expression was knocked down in lung cancer cell lines PC9 and H1299 using RNA interference. Cell proliferation, apoptosis, cell cycle, migration, and invasion abilities were assessed using the CCK-8 assay, colony formation assay, flow cytometry, and Transwell chamber assay, respectively. Phosphoproteomics analysis was performed to identify phosphorylated target proteins of STK39 and associated signaling pathways. PC9 and H1299 cells with knocked-down STK39 were injected into nude mice via the common carotid artery to observe the formation of brain metastases. Finally, RT-qPCR and Western blotting were used to detect the expression of STK39, CPSF4/NFκB/COX2, and epithelial-mesenchymal transition (EMT) markers in lung cancer and brain metastasis tissues, and to analyze the correlation between STK39 expression and the size of metastatic tumors.

RESULTS

STK39 was highly expressed in lung cancer cell lines PC9 and H1299. Knockdown of STK39 inhibited proliferation, migration, and invasion of lung cancer cells, induced apoptosis, and caused cell cycle arrest. Phosphoproteomics and Phos-tag analyses showed that knockdown of STK39 significantly downregulated the expression of phosphorylated CPSF4 protein in PC9 and H1299 cells, along with significant downregulation of NFκB, COX2, and EMT markers. Knockdown of STK39 inhibited the formation of brain metastases by PC9 and H1299 cells in nude mice. Lung cancer brain metastasis tissues exhibited high expression of STK39, CPSF4, NFκB, and COX2, with their expression levels showing a significant positive correlation with the size of metastatic tumors.

CONCLUSION

STK39 is highly expressed in lung cancer brain metastasis tissues, and knockdown of STK39 significantly inhibits brain metastasis in experimental models, accompanied by the suppression of the CPSF4/NFκB/COX2 signaling pathway and EMT process. Therefore, STK39 may be a key factor promoting lung cancer brain metastasis and a potential therapeutic target.

Genetic Ancestry and Lung Cancer in Latin American Patients: A Crucial Step for Understanding a Diverse Population

Lung cancer is the second leading cause of cancer-related deaths in Latin America. While incidence and mortality rates are higher in other populations, the ``Hispanic paradox'' observed in US Hispanics reflects a lower mortality rate for mortality from non-small cell lung cancer (NSCLC) despite socioeconomic disparities, which may be related to epigenetic and cultural factors. Genetic studies have identified single nucleotide polymorphisms associated with ancestry as key contributors to lung cancer risk and outcomes, emphasizing the importance of genomic insights for early detection and personalized treatments. This narrative review explores the impact of genetic ancestry on lung cancer in Hispanic/Latino populations. We searched MEDLINE and Google Scholar for "((SNP) OR (germline) OR (variant)) AND (lung cancer) AND ((Hispanic) OR (Latin))," focusing on Latin American studies. We included articles published up to December 2024. Specific variation in genes such as XRCC1, CYP1A1, CYP1A2, SEMA3B, PADPRP, and mEPHX have been associated with increased lung cancer risk. Lung cancer incidence and prognosis vary significantly among Hispanics due to their diverse genetic ancestry. Understanding ancestry-specific genetic variations may help personalize treatment and improve outcomes for this population.

The role of EMG1 in lung adenocarcinoma progression: Implications for prognosis and immune cell infiltration

BACKGROUND AND AIMS

Lung adenocarcinoma (LUAD) is the most common and aggressive cancer with a high incidence. N1-specific pseudouridine methyltransferase (EMG1), a highly conserved nucleolus protein, plays an important role in the biological development of ribosomes. However, the role of EMG1 in the progression of LUAD is still unclear.

METHODS

The expression of EMG1 in LUAD cells, and LUAD tissues, and adjacent noncancerous tissues was quantified using real-time polymerase chain reaction (PCR) and western blotting. The roles of EMG1 in LUAD cell proliferation, migration, invasion and tumorigenicity were explored in vitro and in vivo. Western blot analysis to underlying molecular mechanism of EMG1 regulating the biological function of LUAD. EMG1 expression and its impact on tumor prognosis were analyzed using a range of databases including GEPIA, UALCAN, cBioPortal, LinkedOmics, and Kaplan-Meier Plotter.

RESULTS

EMG1 expression was elevated in LUAD patients compared to normal tissues, and EMG1 expression was strongly correlated with prognosis in LUAD patients. EMG1 expression correlated with age, gender, N stage, T stage, and pathologic stage. EMG1 expression was strongly positively correlated with MRPL51, PHB2, SNRPG, ATP5MD, and TPI1, and strongly negatively correlated with MACF1, DOCK9, RAPGEF2, SYNJ1, and KIDINS220, the major enrichment pathways for EMG1 and related genes include Cell cycle, DNA Replication and Pathways in cancer signaling pathways. EMG1 expression level was significantly increased in LUAD cell lines and tissues. Knockdown of EMG1 could inhibit LUAD cell proliferation, migration, invasion, and tumorigenicity. Besides, EMG1 overexpression could promote LUAD cell proliferation, migration, and invasion. High expression of EMG1 predicts poor prognosis in LUAD patients, and EMG1 may play an oncogenic role in the tumor microenvironment by participating in the infiltration of LUAD immune cells.

CONCLUSIONS

EMG1 regulated various functions in LUAD by directly mediating Akt/mTOR/p70s6k signaling pathways activation. The results suggest that EMG1 may be a novel biomarker for assessing prognosis and immune cell infiltration in LUAD.

Lung cancer in patients who have never smoked - an emerging disease

Lung cancer is the most common cause of cancer-related deaths globally. Although smoking-related lung cancers continue to account for the majority of diagnoses, smoking rates have been decreasing for several decades. Lung cancer in individuals who have never smoked (LCINS) is estimated to be the fifth most common cause of cancer-related deaths worldwide in 2023, preferentially occurring in women and Asian populations. As smoking rates continue to decline, understanding the aetiology and features of this disease, which necessitate unique diagnostic and treatment paradigms, will be imperative. New data have provided important insights into the molecular and genomic characteristics of LCINS, which are distinct from those of smoking-associated lung cancers and directly affect treatment decisions and outcomes. Herein, we review the emerging data regarding the aetiology and features of LCINS, particularly the genetic and environmental underpinnings of this disease as well as their implications for treatment. In addition, we outline the unique diagnostic and therapeutic paradigms of LCINS and discuss future directions in identifying individuals at high risk of this disease for potential screening efforts.

Genome-wide analysis study of gestational diabetes mellitus and related pathogenic factors in a Chinese Han population

BACKGROUND

Gestational diabetes mellitus (GDM) affects the metabolism of both the mother and fetus during and after pregnancy. Genetic factors are important in the pathogenesis of GDM, and associations vary by ethnicity. However, related studies about the relationship between the susceptibility genes and glucose traits remain limited in China. This study aimed to identify genes associated with GDM susceptibility in Chinese Han women and validate those findings using clinical data during pregnancy and postpartum period.

METHODS

A genome-wide association study (GWAS) of 398 Chinese Han women (199 each with and without GDM) was conducted and associations between single nucleotide polymorphisms (SNPs) and glucose metabolism were identified by searching public databases. Relationships between filtered differential SNPs and glucose metabolism were verified using clinical data during pregnancy. The GDM group were followed up postpartum to evaluate the progression of glucose metabolism.

RESULTS

We identified five novel SNPs with genome-wide significant associations with GDM: rs62069863 in TRPV3 gene and rs2232016 in PRMT6 gene were positive correlated with 1 h plasma glucose (1hPG) and 2 h plasma glucose (2hPG), rs1112718 in HHEX/EXOC6 gene and rs10460009 in LPIN2 gene were positive associated with fasting plasma glucose, 1hPG and 2hPG, rs927316 in GLIS3 gene was negative correlated with 2hPG. Of the 166 GDM women followed up postpartum, rs62069863 in TRPV3 gene was positively associated with fasting insulin, homoeostasis model assessment of insulin resistance.

CONCLUSIONS

The variants of rs62069863 in TRPV3 gene, rs2232016 in PRMT6 gene, rs1112718 in HHEX/EXOC6 gene, rs927316 in GLIS3 gene, and rs10460009 in LPIN2 gene were newly-identified susceptibility loci for GDM in the Chinese Han population. TRPV3 was associated with worse insulin resistance postpartum.

TRIAL REGISTRATION

This study was registered in the Chinese Clinical Trial Registry.

TRIAL REGISTRATION NUMBER

ChiCTR2100043762. Date of first registration: 28/02/2021.

To RNA-binding and beyond: Emerging facets of the role of Rbfox proteins in development and disease

The RNA-binding Fox-1 homologue (Rbfox) proteins represent an ancient family of splicing factors, conserved through evolution. All members share an RNA recognition motif (RRM), and a particular affinity for the GCAUG signature in target RNA molecules. The role of Rbfox, as a splice factor, deciding the tissue-specific inclusion/exclusion of an exon, depending on its binding position on the flanking introns, is well known. Rbfox often acts in concert with other splicing factors, and forms splicing regulatory networks. Apart from this canonical role, recent studies show that Rbfox can also function as a transcription co-factor, and affects mRNA stability and translation. The repertoire of Rbfox targets is vast, including genes involved in the development of tissue lineages, such as neurogenesis, myogenesis, and erythropoeiesis, and molecular processes, including cytoskeletal dynamics, and calcium handling. A second layer of complexity is added by the fact that Rbfox expression itself is regulated by multiple mechanisms, and, in vertebrates, exhibits tissue-specific expression. The optimum dosage of Rbfox is critical, and its misexpression is etiological to various disease conditions. In this review, we discuss the contextual roles played by Rbfox as a tissue-specific regulator for the expression of many important genes with diverse functions, through the lens of the emerging data which highlights its involvement in many human diseases. Furthermore, we explore the mechanistic details provided by studies in model organisms, with emphasis on the work with Drosophila. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > Splicing Regulation/Alternative Splicing.

Pan-cancer analysis of mRNA stability for decoding tumour post-transcriptional programs

Measuring mRNA decay in tumours is a prohibitive challenge, limiting our ability to map the post-transcriptional programs of cancer. Here, using a statistical framework to decouple transcriptional and post-transcriptional effects in RNA-seq data, we uncover the mRNA stability changes that accompany tumour development and progression. Analysis of 7760 samples across 18 cancer types suggests that mRNA stability changes are ~30% as frequent as transcriptional events, highlighting their widespread role in shaping the tumour transcriptome. Dysregulation of programs associated with >80 RNA-binding proteins (RBPs) and microRNAs (miRNAs) drive these changes, including multi-cancer inactivation of RBFOX and miR-29 families. Phenotypic activation or inhibition of RBFOX1 highlights its role in calcium signaling dysregulation, while modulation of miR-29 shows its impact on extracellular matrix organization and stemness genes. Overall, our study underlines the integral role of mRNA stability in shaping the cancer transcriptome, and provides a resource for systematic interrogation of cancer-associated stability pathways.

The role of mRNA stability in shaping the cancer transcriptome is revealed using a statistical analysis of transcriptomic data.

Structure-activity relationship studies of allosteric inhibitors of EYA2 tyrosine phosphatase

Human eyes absent (EYA) proteins possess Tyr phosphatase activity, which is critical for numerous cancer and metastasis promoting activities, making it an attractive target for cancer therapy. In this work, we demonstrate that the inhibitor-bound form of EYA2 does not favour binding to Mg2+ , which is indispensable for the Tyr phosphatase activity. We further describe characterization and optimization of this class of allosteric inhibitors. A series of analogues were synthesized to improve potency of the inhibitors and to elucidate structure-activity relationships. Two co-crystal structures confirm the binding modes of this class of inhibitors. Our medicinal chemical, structural, biochemical, and biophysical studies provide insight into the molecular interactions of EYA2 with these allosteric inhibitors. The compounds derived from this study are useful for exploring the function of the Tyr phosphatase activity of EYA2 in normal and cancerous cells and serve as reference compounds for screening or developing allosteric phosphatase inhibitors. Finally, the co-crystal structures reported in this study will aid in structure-based drug discovery against EYA2.

Associations between genetic loci, environment factors and mental disorders: a genome-wide survival analysis using the UK Biobank data

It is well-accepted that both environment and genetic factors contribute to the development of mental disorders (MD). However, few genetic studies used time-to-event data analysis to identify the susceptibility genetic variants associated with MD and explore the role of environment factors in these associations. In order to detect novel genetic loci associated with MD based on the time-to-event data and identify the role of environmental factors in them, this study recruited 376,806 participants from the UK Biobank cohort. The MD outcomes (including overall MD status, anxiety, depression and substance use disorders (SUD)) were defined based on in-patient hospital, self-reported and death registry data collected in the UK Biobank. SPACOX approach was used to identify the susceptibility loci for MD using the time-to-event data of the UK Biobank cohort. And then we estimated the associations between identified candidate loci, fourteen environment factors and MD through a phenome-wide association study and mediation analysis. SPACOX identified multiple candidate loci for overall MD status, depression and SUD, such as rs139813674 (P value = 8.39 × 10^-9, ZNF684) for overall MD status, rs7231178 (DCC, P value = 2.11 × 10^-9) for depression, and rs10228494 (FOXP2, P value = 6.58 × 10^-10) for SUD. Multiple environment factors could influence the associations between identified loci and MD, such as confide in others and felt hated. Our study identified novel candidate loci for MD, highlighting the strength of time-to-event data based genetic association studies. We also observed that multiple environment factors could influence the association between susceptibility loci and MD.

Concomitant genetic alterations are associated with plasma D-dimer level in patients with non-small cell lung cancer

Objective: D-dimer is correlated to the poor prognosis of non-small cell lung cancer. The study aimed to investigate the association between plasma D-dimer and concomitant mutations in non-small cell lung cancer. Methods: A total of 517 non-small cell lung cancer patients were recruited and tested for ALK, BRAF, EGFR, HER2/ERBB2, KRAS, MET, PIK3CA, RET and ROS1 mutation by next-generation sequencing. Multiple gene mutation information, clinical baseline data and laboratory test data were analyzed statistically. Results: All patients were divided into three groups: wild-type group, single-gene mutation group and concomitant mutation group. The analysis of D-dimer, uric acid, gender, family history, smoking history, histology and distant metastasis all showed significant differences in the three groups (p < 0.05). D-dimer was considered as a risk factor for concomitant mutations according to the unordered multiple logistic regression analysis. The receiver operating characteristic curve analysis indicated that D-dimer had an important predictive value for the occurrence of concomitant mutations (AUC: 0.94; sensitivity: 88.71%; specificity: 86.46). There was significantly shorter median progression-free survival in the concomitant mutation group compared with the single mutation group (7.70 months vs 14.00 months; p = 0.0133). Conclusion: Plasma D-dimer is significantly associated with concomitant mutations and may be regarded as a potent predictor of concomitant mutations for non-small cell lung cancer patients.

Inference for High-Dimensional Censored Quantile Regression

With the availability of high dimensional genetic biomarkers, it is of interest to identify heterogeneous effects of these predictors on patients' survival, along with proper statistical inference. Censored quantile regression has emerged as a powerful tool for detecting heterogeneous effects of covariates on survival outcomes. To our knowledge, there is little work available to draw inference on the effects of high dimensional predictors for censored quantile regression. This paper proposes a novel procedure to draw inference on all predictors within the framework of global censored quantile regression, which investigates covariate-response associations over an interval of quantile levels, instead of a few discrete values. The proposed estimator combines a sequence of low dimensional model estimates that are based on multi-sample splittings and variable selection. We show that, under some regularity conditions, the estimator is consistent and asymptotically follows a Gaussian process indexed by the quantile level. Simulation studies indicate that our procedure can properly quantify the uncertainty of the estimates in high dimensional settings. We apply our method to analyze the heterogeneous effects of SNPs residing in lung cancer pathways on patients' survival, using the Boston Lung Cancer Survivor Cohort, a cancer epidemiology study on the molecular mechanism of lung cancer.

Scalable and Robust Regression Methods for Phenome-Wide Association Analysis on Large-Scale Biobank Data

With the advances in genotyping technologies and electronic health records (EHRs), large biobanks have been great resources to identify novel genetic associations and gene-environment interactions on a genome-wide and even a phenome-wide scale. To date, several phenome-wide association studies (PheWAS) have been performed on biobank data, which provides comprehensive insights into many aspects of human genetics and biology. Although inspiring, PheWAS on large-scale biobank data encounters new challenges including computational burden, unbalanced phenotypic distribution, and genetic relationship. In this paper, we first discuss these new challenges and their potential impact on data analysis. Then, we summarize approaches that are scalable and robust in GWAS and PheWAS. This review can serve as a practical guide for geneticists, epidemiologists, and other medical researchers to identify genetic variations associated with health-related phenotypes in large-scale biobank data analysis. Meanwhile, it can also help statisticians to gain a comprehensive and up-to-date understanding of the current technical tool development.

Raddeanin A induced apoptosis of non-small cell lung cancer cells by promoting ROS-mediated STAT3 inactivation

PURPOSE

Non-small cell lung cancer (NSCLC) is a high-risk type of lung cancer. Raddeanin A exerts anti-tumor activity by regulating cell proliferation and apoptosis, but its role in NSCLC remains to be elucidated. This study was to investigate the effect of raddeanin A in NSCLC and its mechanism.

METHODS

The effect of raddeanin A (2, 4, 8, 10 μmol/L) on the viability, proliferation and apoptosis of A549 and H1299 cells was determined by cell counting kit-8, colony formation and flow cytometry assays, respectively. Next, western blot was performed to examine the protein expressions of cleaved caspase-3, Bax, phosphorylated signal transducer and activator of transcription 3 (p-STAT3) and STAT3. Subsequently, the intracellular reactive oxygen species (ROS) generation and mitochondrial membrane potential of NSCLC cells were detected by 2', 7'-dichlorofluorescein-diacetate (DCFH-DA) and JC-1 assay. Lastly, the effect of N-acetylcysteine (NAC) on the apoptosis, ROS generation, and STAT3 was evaluated by the above-mentioned assays again.

RESULTS

Raddeanin A treatment had no obvious effect on 16HBE cells viability, but it inhibited viability and proliferation of A549 and H1299 cells, promoted the apoptosis, increased the protein expressions of cleaved caspase-3 and Bax, generated intracellular ROS, as well as decreased mitochondrial membrane potential and the expressions of p-STAT3 and STAT3 in A549 and H1299 cells. After cells treated with NAC, the effect of raddeanin A was reversed, as evidenced by the apoptosis and ROS generation were suppressed, and the expression of p-STAT3 was promoted.

CONCLUSION

Raddeanin A suppressed the proliferation and induced apoptosis of NSCLC cells via promoting the ROS-mediated STAT3 inactivation.

Potentially functional variants of ERAP1, PSMF1 and NCF2 in the MHC-I-related pathway predict non-small cell lung cancer survival

BACKGROUND

Cellular immunity against tumor cells is highly dependent on antigen presentation by major histocompatibility complex class I (MHC-I) molecules. However, few published studies have investigated associations between functional variants of MHC-I-related genes and clinical outcomes of lung cancer patients.

METHODS

We performed a two-phase Cox proportional hazards regression analysis by using two previously published genome-wide association studies to evaluate associations between genetic variants in the MHC-I-related gene set and the survival of non-small cell lung cancer (NSCLC) patients, followed by expression quantitative trait loci analysis.

RESULTS

Of the 7811 single-nucleotide polymorphisms (SNPs) in 89 genes of 1185 NSCLC patients in the discovery dataset of the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial, 24 SNPs remained statistically significant after validation in additional 984 NSCLC patients from the Harvard Lung Cancer Susceptibility Study. In a multivariate stepwise Cox model, three independent functional SNPs (ERAP1 rs469783 T > C, PSMF1 rs13040574 C > A and NCF2 rs36071574 G > A) remained significant with an adjusted hazards ratio (HR) of 0.83 [95% confidence interval (CI) = 0.77-0.89, P = 8.0 × 10^-7], 0.86 (0.80-0.93, P = 9.4 × 10^-5) and 1.31 (1.11-1.54, P = 0.001) for overall survival (OS), respectively. Further combined genotypes revealed a poor survival in a dose-response manner in association with the number of unfavorable genotypes (P_trend < 0.0001 and 0.0002 for OS and disease-specific survival, respectively). Also, ERAP1 rs469783C and PSMF1 rs13040574A alleles were associated with higher mRNA expression levels of their genes.

CONCLUSION

These potentially functional SNPs of the MHC-I-related genes may be biomarkers for NSCLC survival, possibly through modulating the expression of corresponding genes.

Genetic Determinants of Lung Cancer Prognosis in Never Smokers: A Pooled Analysis in the International Lung Cancer Consortium

BACKGROUND

Lung cancer remains the leading cause of cancer death worldwide, with 15% to 20% occurring in never smokers. To assess genetic determinants for prognosis among never smokers, we conducted a genome-wide investigation in the International Lung Cancer Consortium (ILCCO).

METHODS

Genomic and clinical data from 1,569 never-smoking patients with lung cancer of European ancestry from 10 ILCCO studies were included. HRs and 95% confidence intervals of overall survival were estimated. We assessed whether the associations were mediated through mRNA expression-based 1,553 normal lung tissues from the lung expression quantitative trait loci (eQTL) dataset and Genotype-Tissue Expression (GTEx). For cross-ethnicity generalization, we assessed the associations in a Japanese study (N = 887).

RESULTS

One locus at 13q22.2 was associated with lung adenocarcinoma survival at genome-wide level, with carriers of rs12875562-T allele exhibiting poor prognosis [HR = 1.71 (1.41-2.07), P = 3.60 × 10-8], and altered mRNA expression of LMO7DN in lung tissue (GTEx, P = 9.40 × 10-7; Lung eQTL dataset, P = 0.003). Furthermore, 2 of 11 independent loci that reached the suggestive significance level (P < 10-6) were significant eQTL affecting mRNA expression of nearby genes in lung tissues, including CAPZB at 1p36.13 and UBAC1 at 9q34.3. One locus encoding NWD2/KIAA1239 at 4p14 showed associations in both European [HR = 0.50 (0.38-0.66), P = 6.92 × 10-7] and Japanese populations [HR = 0.79 (0.67-0.94), P = 0.007].

CONCLUSIONS

Based on the largest genomic investigation on the lung cancer prognosis of never smokers to date, we observed that lung cancer prognosis is affected by inherited genetic variants.

IMPACT

We identified one locus near LMO7DN at genome-wide level and several potential prognostic genes with cis-effect on mRNA expression. Further functional genomics work is required to understand their role in tumor progression.

Potentially functional genetic variants in PLIN2, SULT2A1 and UGT1A9 genes of the ketone pathway and survival of nonsmall cell lung cancer

The ketone metabolism pathway is a principle procedure in physiological homeostasis and induces cancer cells to switch between glycolysis and oxidative phosphorylation for energy production. We conducted a two-phase analysis for associations between genetic variants in the ketone metabolism pathway genes and survival of nonsmall cell lung cancer (NSCLC) by analyzing genotyping data from two published genome-wide association studies (GWASs). In the discovery, we used a genotyping dataset from the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial in the multivariable Cox proportional hazards regression analysis. We used Bayesian false discovery probability (≤0.80) for multiple testing correction to evaluate associations between 25,819 (2,176 genotyped and 23,643 imputed) single-nucleotide polymorphisms (SNPs) in 162 genes and survival of 1,185 NSCLC patients. Subsequently, we validated the identified significant SNPs with an additional 984 NSCLC patients from the Harvard Lung Cancer Susceptibility GWAS study. Finally, we found that three independent and potentially functional SNPs in three different genes (i.e., PLIN2 rs7867814 G>A, SULT2A1 rs2547235 C>T and UGT1A9 rs2011404 C>T) were independently associated with risk of death from NSCLC, with a combined hazards ratio of 1.22 [95% confidence interval = 1.09-1.36 and p = 0.0003], 0.82 (0.74-0.91 and p = 0.0002) and 1.21 (1.10-1.33 and p = 0.0001), respectively. Additional expression quantitative trait loci analysis found that the survival-associated PLIN2 rs7867814 GA + AA genotypes, but not the genotypes of other two SNPs, were significantly associated with increased mRNA expression levels (p = 0.005). These results indicated that PLIN2 variants may be potential predictors of NSCLC survival through regulating the PLIN2 expression.

Set-based genetic association and interaction tests for survival outcomes based on weighted V statistics

With advancements in high-throughout technologies, studies have been conducted to investigate the role of massive genetic variants in human diseases. While set-based tests have been developed for binary and continuous disease outcomes, there are few computationally efficient set-based tests available for time-to-event outcomes. To facilitate the genetic association and interaction analyses of time-to-event outcomes, We develop a suite of multivariant tests based on weighted V statistics with or without considering potential genetic heterogeneity. In addition to the computation efficiency and nice asymptotic properties, all the new tests can deal with left truncation and competing risks in the survival data, and adjust for covariates. Simulation studies show that the new tests run faster, are more accurate in small samples, and account for confounding effect better than the existing multivariant survival tests. When the genetic effect is heterogeneous across individuals/subpopulations, the association test considering genetic heterogeneity is more powerful than the existing tests that do not account for genetic heterogeneity. Using the new methods, we perform a genome-wide association analysis of the genotype and age-to-Alzheimer's data from the Rush Memory and Aging Project and the Religious Orders Study. The analysis identifies two genes, APOE and APOC1, associated with age to Alzheimer's disease onset.

A Fast and Accurate Method for Genome-Wide Time-to-Event Data Analysis and Its Application to UK Biobank

With increasing biobanking efforts connecting electronic health records and national registries to germline genetics, the time-to-event data analysis has attracted increasing attention in the genetics studies of human diseases. In time-to-event data analysis, the Cox proportional hazards (PH) regression model is one of the most used approaches. However, existing methods and tools are not scalable when analyzing a large biobank with hundreds of thousands of samples and endpoints, and they are not accurate when testing low-frequency and rare variants. Here, we propose a scalable and accurate method, SPACox (a saddlepoint approximation implementation based on the Cox PH regression model), that is applicable for genome-wide scale time-to-event data analysis. SPACox requires fitting a Cox PH regression model only once across the genome-wide analysis and then uses a saddlepoint approximation (SPA) to calibrate the test statistics. Simulation studies show that SPACox is 76-252 times faster than other existing alternatives, such as gwasurvivr, 185-511 times faster than the standard Wald test, and more than 6,000 times faster than the Firth correction and can control type I error rates at the genome-wide significance level regardless of minor allele frequencies. Through the analysis of UK Biobank inpatient data of 282,871 white British European ancestry samples, we show that SPACox can efficiently analyze large sample sizes and accurately control type I error rates. We identified 611 loci associated with time-to-event phenotypes of 12 common diseases, of which 38 loci would be missed within a logistic regression framework with a binary phenotype defined as event occurrence status during the follow-up period.

Silencing of PRDX2 Inhibits the Proliferation and Invasion of Non-Small Cell Lung Cancer Cells

Peroxiredoxin 2 (PRDX2), a member of the peroxiredoxin family of antioxidant enzymes, has been revealed to be an important player in cancer progression. However, the biological role of PRDX2 in the progression of non-small cell lung cancer (NSCLC) is poor reported. In the present study, the loss-of-function experiments were performed to investigate the specific role of PRDX2 in the growth and invasion of NSCLC. The results revealed that knockdown of PRDX2 by siRNA interference significantly suppressed the proliferation, migration, and invasion of A549 and H1299 cells, as well as diminished the activity of MMP9. Additionally, the decrease in PRDX2 expression significantly promoted apoptosis in NSCLC cells by downregulating expression of Bcl-2 and upregulating the expression of Bax, cleaved caspase 3 and cleaved caspase 9, but had no significant effect on the apoptosis of normal lung epithelial cells BEAS-2B. Moreover, PRDX2 inhibitor also inhibited the proliferation, migration, and invasion of A549 cells and promoted apoptosis. Further, our data demonstrated that silencing of PRDX2 markedly reduced the phosphorylation of Akt and mTOR and expression of downstream proteins Cyclin D1 and p70S6k. In conclusion, our findings indicate that PRDX2 exerts a prooncogenic role in the progression of NSCLC and might be a potential therapeutic target for NSCLC treatment.

Copy number variations associated with fetal congenital kidney malformations

Background

Congenital anomalies of the kidney and urinary tract (CAKUT) constitute 20–30% of all congenital malformations. Within the CAKUT phenotypic spectrum, renal hypodysplasia (RHD) is particularly severe. This study aimed to evaluate the applicability of single-nucleotide polymorphism (SNP) array test in prenatal diagnosis of RHD for improving prenatal genetic counseling and to search for evidence of a possible causative role of copy-number variations (CNVs) in RHD.

Results

We performed a systematic survey of CNV burden in 120 fetuses with RHD: 103 cases were isolated RHD and 17 were non-isolated RHD. Single-nucleotide polymorphism (SNP) array test was performed using the Affymetrix CytoScan HD platform. All annotated CNVs were validated by fluorescence in situ hybridization. We identified abnormal CNVs in 15 (12.5%) cases of RHD; of these CNVs, 11 were pathogenic and 4 were variants of uncertain significance. The detection rate of abnormal CNVs in non-isolated RHD was higher (29.4%, 5/17) than that in isolated RHD (9.7%, 10/103) (P = 0.060). Parents are more inclined to terminate the pregnancy if the fetuses have pathogenic results of the SNP-array test.

Conclusions

The variable phenotypes that abnormal CNVs may cause indicate the genetic counseling is needed for RHD cases.

MicroRNA‑299‑5p inhibits cell metastasis in breast cancer by directly targeting serine/threonine kinase 39

Numerous studies have demonstrated that microRNAs (miRNAs) play a key role in human carcinogenesis and metastasis. For example, miR‑299‑5p has previously been revealed to be dysregulated in several human cancers. However, the biological function of miR‑299‑5p in breast cancer remains unclear. The present study demonstrated that miR‑299‑5p was downregulated in breast cancer tissues and cell lines. The restoration of miR‑299‑5p expression suppressed cell migration and invasion, whereas inhibition of miR‑299‑5p promoted cell migration and invasion. In addition, in vivo studies demonstrated that miR‑299‑5p overexpression was able to inhibit tumour metastasis in nude mice. Mechanistically, through bioinformatics analysis and a dual‑luciferase assay, it was confirmed that miR‑299‑5p directly targets serine/threonine kinase 39 (STK39). Silencing STK39 inhibited cell metastasis and suppressed epithelial‑mesenchymal transition markers and matrix metalloproteinase expression, whereas restoration of STK39 expression was able to reverse miR‑299‑5p‑inhibited cell migration and invasion. Collectively, the results of the present study demonstrated that miR‑299‑5p supresses breast cancer cell migration and invasion by targeting STK39. These findings may provide novel insights into miR‑299‑5p and its potential diagnostic and therapeutic benefits in breast cancer.

Polymorphisms of rs4787050 and rs8045980 are associated with lung cancer risk in northeast Chinese female nonsmokers

Aim: We studied the association between two single-nucleotide polymorphisms (SNPs: rs4787050 and rs8045980) in RBFOX1 and lung cancer risk, and explored the interaction between the two SNPs and exposure to cooking oil fume on lung cancer risk in northeast Chinese female nonsmokers. Methods: Northeast Chinese female nonsmokers were enrolled into the study (people with lung cancer, 647; people without lung cancer, 675). All statistical analyses were performed using SPSS software. Results: The SNPs rs4787050 and rs8045980 showed a significant association with susceptibility to lung cancer. Moreover, cooking oil fume exposure was found to increase the risk of lung cancer. However, no gene-environment interactions were discovered. Conclusion: The present study revealed that rs4787050 and rs8045980 in RBFOX1 may be meaningful as a novel biomarker for lung cancer susceptibility.

Structural and Functional Analyses of an Allosteric EYA2 Phosphatase Inhibitor That Has On-Target Effects in Human Lung Cancer Cells

EYA proteins (EYA1-4) are critical developmental transcriptional cofactors that contain an EYA domain (ED) harboring Tyr phosphatase activity. EYA proteins are largely downregulated after embryogenesis but are reexpressed in cancers, and their Tyr phosphatase activity plays an important role in the DNA damage response and tumor progression. We previously identified a class of small-molecule allosteric inhibitors that specifically inhibit the Tyr phosphatase activity of EYA2. Herein, we determined the crystal structure of the EYA2 ED in complex with NCGC00249987 (a representative compound in this class), revealing that it binds to an induced pocket distant from the active site. NCGC00249987 binding leads to a conformational change of the active site that is unfavorable for Mg²⁺ binding, thereby inhibiting EYA2's Tyr phosphatase activity. We demonstrate, using genetic mutations, that migration, invadopodia formation, and invasion of lung adenocarcinoma cells are dependent on EYA2 Tyr phosphatase activity, whereas growth and survival are not. Further, we demonstrate that NCGC00249987 specifically targets migration, invadopodia formation, and invasion of lung cancer cells, but that it does not inhibit cell growth or survival. The compound has no effect on lung cancer cells carrying an EYA2 F290Y mutant that abolishes compound binding, indicating that NCGC00249987 is on target in lung cancer cells. These data suggest that the NCGC00249987 allosteric inhibitor can be used as a chemical probe to study the function of the EYA2 Tyr phosphatase activity in cells and may have the potential to be developed into an antimetastatic agent for cancers reliant on EYA2's Tyr phosphatase activity.

Tumor-associated antigens identified early in mouse mammary tumor development can be effective vaccine targets

Breast cancer vaccines composed of antigens identified by serological analysis of cDNA expression libraries (SEREX) induce antigen specific immune responses in patients but have had disappointing clinical benefits. While many attempts to modify the adjuvants and vaccine method have been tried, one issue not addressed was whether the SEREX tumor-associated antigens identified from late stages of disease were ideal targets. We questioned in the transgenic TgMMTV-neu mouse model whether the antigen repertoire is distinct between early and late stage breast cancer and whether the antigens identified via SEREX from transgenic mice with early or late stage tumors would elicit differential anti-tumor effects to address this question. Three early stage antigens, Pdhx, Stk39, and Otud6B, were identified from a SEREX screen of mice prior to development of palpable lesions. Formulated into a vaccine, each early antigen inhibited tumor growth (p < 0.0001). The antigens identified from mice with late stage tumors (Swap70, Gsn, and Arhgef2) were unable to inhibit tumor growth when used as vaccines (for example Gsn p = 0.26). Each of the three early stage antigens were essential for tumor survival in syngeneic mouse tumor cells and in human breast cancer cell lines across breast cancer subtypes. Silencing protein expression of the early antigens increased apoptosis (p < 0.0001 for all antigens in mouse and p < 0.05 for all antigens in human triple negative breast cancer) and decreased survival (p < 0.0001 for all antigens in mouse and human triple negative and HER2 positive breast cancer). Overexpression of the early stage antigens in women with breast cancer predicted worse prognosis (p = 0.03) while overexpression of late stage antigens did not impact prognosis (p = 0.09). These data suggest that antigens expressed earlier in breast tumor development and functionally relevant to breast tumor growth may be more effective targets for therapeutic breast cancer vaccines than antigens identified in later disease.

Predictive value of single-nucleotide polymorphism signature for recurrence in localised renal cell carcinoma: a retrospective analysis and multicentre validation study

BACKGROUND

Identification of high-risk localised renal cell carcinoma is key for the selection of patients for adjuvant treatment who are at truly higher risk of reccurrence. We developed a classifier based on single-nucleotide polymorphisms (SNPs) to improve the predictive accuracy for renal cell carcinoma recurrence and investigated whether intratumour heterogeneity affected the precision of the classifier.

METHODS

In this retrospective analysis and multicentre validation study, we used paraffin-embedded specimens from the training set of 227 patients from Sun Yat-sen University (Guangzhou, Guangdong, China) with localised clear cell renal cell carcinoma to examine 44 potential recurrence-associated SNPs, which were identified by exploratory bioinformatics analyses of a genome-wide association study from The Cancer Genome Atlas (TCGA) Kidney Renal Clear Cell Carcinoma (KIRC) dataset (n=114, 906 600 SNPs). We developed a six-SNP-based classifier by use of LASSO Cox regression, based on the association between SNP status and patients' recurrence-free survival. Intratumour heterogeneity was investigated from two other regions within the same tumours in the training set. The six-SNP-based classifier was validated in the internal testing set (n=226), the independent validation set (Chinese multicentre study; 428 patients treated between Jan 1, 2004 and Dec 31, 2012, at three hospitals in China), and TCGA set (441 retrospectively identified patients who underwent resection between 1998 and 2010 for localised clear cell renal cell carcinoma in the USA). The main outcome was recurrence-free survival; the secondary outcome was overall survival.

FINDINGS

Although intratumour heterogeneity was found in 48 (23%) of 206 cases in the internal testing set with complete SNP information, the predictive accuracy of the six-SNP-based classifier was similar in the three different regions of the training set (areas under the curve [AUC] at 5 years: 0·749 [95% CI 0·660-0·826] in region 1, 0·734 [0·651-0·814] in region 2, and 0·736 [0·649-0·824] in region 3). The six-SNP-based classifier precisely predicted recurrence-free survival of patients in three validation sets (hazard ratio [HR] 5·32 [95% CI 2·81-10·07] in the internal testing set, 5·39 [3·38-8·59] in the independent validation set, and 4·62 [2·48-8·61] in the TCGA set; all p<0·0001), independently of patient age or sex and tumour stage, grade, or necrosis. The classifier and the clinicopathological risk factors (tumour stage, grade, and necrosis) were combined to construct a nomogram, which had a predictive accuracy significantly higher than that of each variable alone (AUC at 5 years 0·811 [95% CI 0·756-0·861]).

INTERPRETATION

Our six-SNP-based classifier could be a practical and reliable predictor that can complement the existing staging system for prediction of localised renal cell carcinoma recurrence after surgery, which might enable physicians to make more informed treatment decisions about adjuvant therapy. Intratumour heterogeneity does not seem to hamper the accuracy of the six-SNP-based classifier as a reliable predictor of recurrence. The classifier has the potential to guide treatment decisions for patients at differing risks of recurrence.

FUNDING

National Key Research and Development Program of China, National Natural Science Foundation of China, Guangdong Provincial Science and Technology Foundation of China, and Guangzhou Science and Technology Foundation of China.

Down-regulation of FOXR2 inhibits non-small cell lung cancer cell proliferation and invasion through the Wnt/β-catenin signaling pathway

Forkhead box R2 (FOXR2), a new member of the FOX family, is an important player in a wide range of cellular processes such as proliferation, migration, differentiation and apoptosis. Recently, FOXR2 has been reported to be implicated in cancer development. However, the biological functions of FOXR2 in non-small cell lung cancer (NSCLC) remain unclear. In this study, we investigated the specific role of FOXR2 in NSCLC. The results showed that down-regulation of FOXR2 significantly inhibited NSCLC cell proliferation and invasion in vitro and suppressed NSCLC cell growth and metastasis in vivo. In addition, the decrease in FOXR2 expression markedly reduced the protein levels of β-catenin, cyclinD1 and c-Myc and hence inactivated the Wnt/β-catenin pathway in NSCLC cells. Taken together, we concluded that FOXR2 might be considered as a promising therapeutic target for NSCLC treatment.

A Decade of GWAS Results in Lung Cancer

Genome-wide association studies (GWAS) were successful to identify genetic factors robustly associated with lung cancer. This review aims to synthesize the literature in this field and accelerate the translation of GWAS discoveries into results that are closer to clinical applications. A chronologic presentation of published GWAS on lung cancer susceptibility, survival, and response to treatment is presented. The most important results are tabulated to provide a concise overview in one read. GWAS have reported 45 lung cancer susceptibility loci with varying strength of evidence and highlighted suspected causal genes at each locus. Some genetic risk loci have been refined to more homogeneous subgroups of lung cancer patients in terms of histologic subtypes, smoking status, gender, and ethnicity. Overall, these discoveries are an important step for future development of new therapeutic targets and biomarkers to personalize and improve the quality of care for patients. GWAS results are on the edge of offering new tools for targeted screening in high-risk individuals, but more research is needed if GWAS are to pay off the investment. Complementary genomic datasets and functional studies are needed to refine the underlying molecular mechanisms of lung cancer preliminarily revealed by GWAS and reach results that are medically actionable. Cancer Epidemiol Biomarkers Prev; 27(4); 363-79. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Genetic variants in cytokine signaling pathways and clinical outcomes in early-stage lung cancer patients

OBJECTIVES

The purpose of this study was to explore the role of single nucleotide polymorphisms (SNPs) in cytokine signaling genes and to compare them with clinical outcomes in surgical patients with non-small cell lung cancer (NSCLC).

METHODS

SNPs of the cytokine signaling pathway were analyzed using peripheral blood of 722 patients who underwent resection of stage I to III NSCLC between 1995 and 2009. Cox proportional hazard analyses were performed to identify SNPs associated with overall survival (OS) and risk of recurrence. Internal validation using bootstrap analysis selected SNPs for unfavorable genotype and survival tree analysis.

RESULTS

Seventeen and 9 SNPs were independently associated with OS and recurrence, respectively. Patients with ≥9 unfavorable genotypes experienced worse OS (median, 41 months) than patients with 7 to 8 (89 months) and ≤6 (153 months) after median follow-up of 71 months (P = 2.86 × 10^-23). Patients with ≤3 unfavorable genotypes had greater time to recurrence (median not reached) than those with 4 to 6 (114 months) and ≥7 (44 months; P = 1.3 × 10^-5). Survival tree analysis classified patients into 3 risk groups. Patients in the intermediate- (median OS, 82 months) and high-risk groups (43 months) had worse survival than the low-risk group (176 months; P = 5.51 × 10^-20). Median time to recurrence was worse in the intermediate- (114 months) and high-risk groups (58 months) than the low-risk group (median not reached; P = 2.52 × 10^-9).

CONCLUSIONS

Genetic variants in cytokine signaling pathways were associated with clinical outcomes in NSCLC patients treated with surgery individually and cumulatively. Further studies are necessary to elucidate our findings and translate them into the clinical setting.

Role of high expression levels of STK39 in the growth, migration and invasion of non-small cell type lung cancer cells

Non-small cell type lung cancer (NSCLC) is the most common malignancy and the leading cause of cancer related mortality. In this study, serine/threonine kinase 39 (STK39) was identified as an up-regulated gene in NSCLC tissues by next-generation RNA sequencing. Although STK39 gene polymorphisms may be prognostic of overall survival in patients with early stage NSCLC, the roles of STK39 in NSCLC cancer are poorly understood. In the current study, Genome Set Enrichment Analysis (GSEA) on the RNA-seq data of NSCLC specimens indicated that cancer-related process and pathways, including metastasis, cell cycle, apoptosis and p38 pathway, were significantly correlated with STK39 expression. STK39 expression was significantly increased in NSCLC cases and its protein expression was positively correlated with the poor tumor stage, large tumor size, advanced lymphnode metastasis and poor prognosis. Down-regulation of STK39 in NSCLC cells significantly decreased cell proliferation by blocking of cell cycle and inducing apoptosis. We also found that STK39 knockdown in NSCLC cells remarkably repressed cell migration and invasion. On the contrary, overexpression of STK39 in NSCLC cells had inverse effects on cell behaviors. Taken together, STK39 acts as a tumor oncogene in NSCLC and can be a potential biomarker of carcinogenesis.

Low expression of protocadherin7 (PCDH7) is a potential prognostic biomarker for primary non-muscle invasive bladder cancer

Bladder cancer is a heterogeneous disease with outcome difficult to predict, and novel predictive biomarkers are needed. PCDH7, a member of protocadherins family, functions as tumor suppressor in several human cancers. The human PCDH7 gene is localized in chromosome 4p15, which is often inactivated in human cancers, including bladder cancer. The aim of this study was to investigate the clinical significance of PCDH7 expression in non-muscle invasive bladder cancer (NMIBC). PCDH7 expression was examined using immunohistochemical staining in 199 primary NMIBC tissues and 25 normal bladder epithelial tissues. Then the relationship between PCDH7 expression and clinicopathologic features was evaluated. Kaplan-Meier survival analysis and Cox analysis was used to evaluate the correlation between PCDH7 expression and prognosis. PCDH7 expression in NMIBC tissues was significantly lower than that in normal bladder epithelial tissues (P < 0.001). Low PCDH7 expression correlated with advanced grade (P = 0.021) and larger tumor size (P = 0.044). Moreover, patients with low PCDH7 expression have shorter recurrence-free survival (P < 0.001), progression-free survival (P = 0.007) and overall survival (P = 0.011) than patients with high PCDH7 expression. Low PCDH7 expression is an independent predictor of recurrence-free survival (multivariate Cox analysis: P = 0.007), progression-free survival (multivariate Cox analysis: P = 0.014) and overall survival (multivariate Cox analysis: P = 0.004). The findings indicate that low PCDH7 expression is a potential prognostic biomarker for primary NMIBC.

Id-1 promotes migration and invasion of non-small cell lung cancer cells through activating NF-κB signaling pathway

Background

Numerous studies have shown that Id-1 (Inhibitor of differentiation 1) is upregulated in several cancers and associated with tumor malignant characters. However, the clinical significance and biological role of Id-1 in non-small cell lung cancer (NSCLC) remains unclear.

Methods

We used RT-PCR, Western blot and Immunohistochemistry to measure Id-1 expression in NSCLC tissues and matched adjacent noncancerous tissues. The expression pattern of Id-1 in NSCLC tissues was determined by scoring system of immunohistochemical analysis. The Kaplan-Meier method was used to calculate the survival curve, and log-rank test to determine statistical significance. The Id-1 gene was overexpressed or downreuglated with Lentiviral vectors in NSCLC cells. And, the migration ability of NSCLC cells was tested in a Transwell Boyden Chamber.

Results

We found that Id-1 is generally expressed higher in NSCLC tissues compared with matched adjacent noncancerous tissues. We also found that high Id-1 expression in tumor tissues is significantly correlated with tumor progression and poor survival in NSCLC patients. Furthermore, our experimental data revealed that knockdown of Id-1 significantly suppressed the proliferation, migration and invasion of NSCLC cells, whereas ectopic expression of Id-1 promoted the malignant phenotype of NSCLC cells. Mechanistic study showed that NF-κB signaling pathway contributed to the effects of Id-1 in NSCLC cells. Moreover, blocking the NF-κB pathway significantly inhibited the tumor-promoting actions of Id-1 in NSCLC cells.

Conclusions

We identified a tumorigenic role of Id-1 in NSCLC and provided a novel therapeutic target for NSCLC patients.

EYA2 promotes lung cancer cell proliferation by downregulating the expression of PTEN

Lung cancer is the leading cause of cancer-related death worldwide. Despite advances have been made in diagnosis and therapeutic strategies, the prognosis of lung cancer is still very poor. Eyes absent transcriptional cofactor EYA2 has been shown to promote lung cancer cell growth, however, the underlying molecular mechanism is still not fully understood. In the present study, we found that EYA2 was up-regulated in lung cancer, and EYA2 led to increased cell proliferation by inhibiting Phosphatase and tensin homologue (PTEN) expression via modulation of miR-93. Additionally, survival analysis showed that lung cancer patients with higher EYA2 expression predicted a worse prognosis. Therefore, these findings demonstrate that EYA2 may play an important role in lung cancer occurrence and progression. Targeting EYA2 may provide a feasible approach in developing novel anticancer therapeutics.

Genomic Programming of Human Neonatal Dendritic Cells in Congenital Systemic and In Vitro Cytomegalovirus Infection Reveal Plastic and Robust Immune Pathway Biology Responses

Neonates and especially premature infants are highly susceptible to infection but still can have a remarkable resilience that is poorly understood. The view that neonates have an incomplete or deficient immune system is changing. Human neonatal studies are challenging, and elucidating host protective responses and underlying cognate pathway biology, in the context of viral infection in early life, remains to be fully explored. In both resource rich and poor settings, human cytomegalovirus (HCMV) is the most common cause of congenital infection. By using unbiased systems analyses of transcriptomic resources for HCMV neonatal infection, we find the systemic response of a preterm congenital HCMV infection, involves a focused IFN regulatory response associated with dendritic cells. Further analysis of transcriptional-programming of neonatal dendritic cells in response to HCMV infection in culture revealed an early dominant IFN-chemokine regulatory subnetworks, and at later times the plasticity of pathways implicated in cell-cycle control and lipid metabolism. Further, we identify previously unknown suppressed networks associated with infection, including a select group of GPCRs. Functional siRNA viral growth screen targeting 516-GPCRs and subsequent validation identified novel GPCR-dependent antiviral (ADORA1) and proviral (GPR146, RGS16, PTAFR, SCTR, GPR84, GPR85, NMUR2, FZ10, RDS, CCL17, and SORT1) roles. By contrast a gene family cluster of protocadherins is significantly differentially induced in neonatal cells, suggestive of possible immunomodulatory roles. Unexpectedly, programming responses of adult and neonatal dendritic cells, upon HCMV infection, demonstrated comparable quantitative and qualitative responses showing that functionally, neonatal dendritic cell are not overly compromised. However, a delay in responses of neonatal cells for IFN subnetworks in comparison with adult-derived cells are notable, suggestive of subtle plasticity differences. These findings support a set-point control mechanism rather than immaturity for explaining not only neonatal susceptibility but also resilience to infection. In summary, our findings show that neonatal HCMV infection leads to a highly plastic and functional robust programming of dendritic cells in vivo and in vitro. In comparison with adults, a minimal number of subtle quantitative and temporal differences may contribute to variability in host susceptibility and resilience, in a context dependent manner.

Protocadherin 7 inhibits cell migration and invasion through E-cadherin in gastric cancer

The protocadherin 7 is a member of the protocadherin family that expressed aberrantly in many types of human cancers. However, its expression, function, and underlying mechanisms are little known in gastric cancer. In this study, we detected protocadherin 7 expression in gastric cancer tissues and non-tumorous gastric mucosa tissues by real-time quantitative polymerase chain reaction and immunohistochemistry. The association of protocadherin 7 expression with the clinicopathological characteristics and the prognosis was subsequently analyzed. MTS ((3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium)) and transwell assays were performed to assess the effect of protocadherin 7 on proliferation, migration, and invasion in gastric cancer cell lines. Moreover, real-time quantitative polymerase chain reaction and western blot were used to detect the expression of epithelial-mesenchymal transition markers. Protocadherin 7 expression was decreased gradiently from normal tissue to gastric cancer, especially in gastric cancer tissue with lymph node metastasis. Low expression of protocadherin 7 was significantly associated with Lauren's classification ( p = 0.0005), lymph node metastases ( p = 0.0002), and tumor node metastasis stage ( p = 0.0221), as well as poor prognosis ( p < 0.05). Furthermore, down-regulation of protocadherin 7 in gastric cancer cell lines significantly increased their migration and invasion abilities (both p < 0.05), while it had no influence on the gastric cancer cell proliferation ( p > 0.05). Additionally, our results demonstrated that E-cadherin expression was down-regulated in gastric cancer cells with protocadherin 7 depletion. Our data indicated that protocadherin 7 may play important roles in the invasion and metastasis of gastric cancer, and protocadherin 7 could suppress cell migration and invasion through E-cadherin inhibition. Protocadherin 7 can serve as a novel biomarker for diagnostic and prognosis in patients with gastric cancer.

Scavenger Receptor A1 Prevents Metastasis of Non-Small Cell Lung Cancer via Suppression of Macrophage Serum Amyloid A1

Mechanisms of cross-talk between tumor cells and tumor-associated macrophages (TAM), which drive metastasis, are not fully understood. Scavenger receptor A1 (SR-A1) expressed primarily in macrophages has been associated with lung tumorigenesis. In this study, we used population genetics, transcriptomics, and functional analyses to uncover how SR-A1 is involved in lung cancer and its prognosis. SR-A1 genetic variants were investigated for possible association with survival of advanced stage NSCLC patients in the Harvard Lung Cancer Study cohort. Two SNPs (rs17484273, rs1484751) in SR-A1 were associated significantly with poor overall survival in this cohort. Data from The Cancer Genome Atlas showed considerable downregulation of SR-A1 in lung tumor tissues. The association of SR-A1 with prognosis was validated in animal models in the context of lung cancer metastasis. Macrophages derived from mice genetically deficient for SR-A1 exhibited accelerated metastasis in a model of lung cancer. On the other hand, tumor cell seeding, migration, and invasion, as well as macrophage accumulation in lung cancer tissue, were enhanced in SR-A1-deficient mice. SR-A1 deletion upregulated serum amyloid A1 (SAA1) in macrophages via MAPK/IκB/NFκB signaling. SAA1 promoted tumor cell invasion and macrophage migration in vitro and in vivo, but these effects were blocked by administration of an anti-SAA1 antibody. Overall, our findings show how SR-A1 suppresses lung cancer metastasis by downregulating SAA1 production in TAMs. Cancer Res; 77(7); 1586-98. ©2017 AACR.

Exome-Wide Association Study Identifies Low-Frequency Coding Variants in 2p23.2 and 7p11.2 Associated with Survival of Non-Small Cell Lung Cancer Patients

INTRODUCTION

A growing body of evidence has suggested that low-frequency or rare coding variants might have strong effects on the development and prognosis of cancer. Here, we aim to assess the role of low-frequency and rare coding variants in the survival of NSCLC in Chinese populations.

METHODS

We performed an exome-wide scan of 247,870 variants in 1008 patients with NSCLC and replicated the promising variants by using imputed genotype data of The Cancer Genome Atlas (TCGA) with a Cox regression model. Gene-based and pathway-based analysis were also performed for nonsynonymous or splice site variants. Additionally, analysis of gene expression data in the TCGA was used to increase the reliability of candidate loci and genes.

RESULTS

A low-frequency missense variant in chaperonin containing TCP1 subunit 6A gene (CCT6A) (rs33922584: adjusted hazard ratio [HR_adjusted] = 1.75, p = 6.06 × 10^-4) was significantly related to the survival of patients with NSCLC, which was further replicated by the TCGA samples (HR_adjusted = 4.19, p = 0.015). Interestingly, the G allele of rs33922584 was significantly associated with high expression of CCT6A (p = 0.019) that might induce the worse survival in the TCGA samples (HR_adjusted = 1.15, p = 0.047). Besides, rs117512489 in gene phospholipase B1 gene (PLB1) (HR = 2.02, p = 7.28 × 10^-4) was also associated with survival of the patients with NSCLC in our samples, but it was supported only by gene expression analysis in the TCGA (HR_adjusted = 1.15, p = 0.023). Gene-based and pathway-based analysis revealed a total of 32 genes, including CCT6A and 34 potential pathways might account for the survival of NSCLC, respectively.

CONCLUSION

These results provided more evidence for the important role of low-frequency or rare variants in the survival of patients with NSCLC.

PROTOCADHERIN 7 Acts through SET and PP2A to Potentiate MAPK Signaling by EGFR and KRAS during Lung Tumorigenesis

Non-small cell lung cancer (NSCLC) is the leading cause of cancer-associated deaths worldwide. Given the efficacy of membrane proteins as therapeutic targets in human malignancies, we examined cell-surface receptors that may act as drivers of lung tumorigenesis. Here, we report that the PROTOCADHERIN PCDH7 is overexpressed frequently in NSCLC tumors where this event is associated with poor clinical outcome. PCDH7 overexpression synergized with EGFR and KRAS to induce MAPK signaling and tumorigenesis. Conversely, PCDH7 depletion suppressed ERK activation, sensitized cells to MEK inhibitors, and reduced tumor growth. PCDH7 potentiated ERK signaling by facilitating interaction of protein phosphatase PP2A with its potent inhibitor, the SET oncoprotein. By establishing an oncogenic role for PCDH7 in lung tumorigenesis, our results provide a rationale to develop novel PCDH7 targeting therapies that act at the cell surface of NSCLC cells to compromise their growth. Cancer Res; 77(1); 187-97. ©2016 AACR.

Prenatal sonographic diagnosis of urorectal septum malformation sequence and chromosomal microarray analysis: A case report and review of the literature

INTRODUCTION

Urorectal septum malformation sequence (URSMS) is a rare congenital abnormal syndrome that is caused by the incomplete division of the cloaca. Based on whether the cloaca membrane breaks down or not, the URSMS are classified as full and partial forms. The prenatal diagnosis of URSMS remains challenging because of poor recognition to this malformation and the relatively non-specific sonographic features. We report a prenatally sonographic diagnosed case of the partial URSMS, and review the literature to summarize the prenatal features.

CASE REPORT AND REVIEW

A 37-year old woman was referred at 24 weeks of gestation for fetal abdominal cyst. Detailed sonographic examination was done and revealed the vesicocolic fistula, distended colon, absence of perianal hypoechoic ring, pyelectasis, and small stomach bubble. The URSMS was suspected.Amniocentesis was done and karyotyping revealed 46,XY. Furthermore, chromosomal microarray analysis (CMA) was performed for the first time in URSMS and an alteration of 111.8Kb deletion was detected in 16p13.3 which was located inside the RBFOX1 gene. Parental studies showed that the deletion was inherited from the father who has nomal clinical phenotype.The woman elected to terminate the pregnancy at 25 weeks gestation and postmortem examination confirmed the diagnosis of partial URSMS.The published studies were reviewed and 28 cases of URSMS with conducted prenatal ultrasonography were collected in this report. The most common sonographic description, as suspicious signs of URSMS, were severe oligohydramnios or anhydramnios, urinary tract anomalies, fetal intra-abdominal cysts, and dilated bowel. Also, enterolithiasis and vesicocolic fistula were relatively infrequent but highly specific feature of URSMS.

CONCLUSIONS

URSMS is difficult to be diagnosed prenatally. However, it has characteristic features that can be detected by fetal ultrasonography, and a precise prenatal sonographic examination is crucial for diagnosing URSMS. Besides, more genomic profiling studies are needed to elucidate the causality.

Identification of hallmarks of lung adenocarcinoma prognosis using whole genome sequencing

In conjunction with clinical characteristics, prognostic biomarkers are essential for choosing optimal therapies to lower the mortality of lung adenocarcinoma. Whole genome sequencing (WGS) of 7 cancerous-noncancerous tissue pairs was performed to explore the comparative copy number variations (CNVs) associated with lung adenocarcinoma. The frequencies of top ranked CNVs were verified in an independent set of 114 patients and then the roles of target CNVs in disease prognosis were assessed in 313 patients. The WGS yielded 2604 CNVs. After frequency validation and biological function screening of top 10 CNVs, 9 mutant driver genes from 7 CNVs were further analyzed for an association with survival. Compared with the PBXIP1 amplified copy number, unamplified carriers had a 0.62-fold (95%CI = 0.43–0.91) decreased risk of death. Compared with an amplified TERT, those with an unamplified TERT had a 35% reduction (95% CI = 3%–56%) in risk of lung adenocarcinoma progression. Cases with both unamplified PBXIP1 and TERT had a median 34.32-month extension of overall survival and 34.55-month delay in disease progression when compared with both amplified CNVs. This study demonstrates that CNVs of TERT and PBXIP1 have the potential to translate into the clinic and be used to improve outcomes for patients with this fatal disease.

Pharmacogenomics in the treatment of lung cancer: an update

Significant advances have been made in the analysis of the human genome in the first decades of the 21st century and understanding of tumor biology has matured greatly. The identification of tumor-associated mutations and the pathways involved has led to the development of targeted anticancer therapies. However, the challenge now in using chemotherapy to treat nonsmall-cell lung cancer is to identify more molecular markers predictive of drug sensitivity and determine the optimal drug sequences in order to tailor treatment to each patient. This approach could permit selection of patients who could benefit most from a specific type of chemotherapy by matching their tumor and individual genetic profile. Nevertheless, this potential has been limited so far by reliance on the single biomarker approach, though this is now on the way to being overcome through whole genome studies.

Estimation of the Survival of Patients With Lung Squamous Cell Carcinoma Using Genomic Copy Number Aberrations

BACKGROUND

Estimation of the survival of patients with lung squamous cell carcinoma (SCC) on the basis of histopathology is inadequate. The aim of this study was to identify genomic regions with potential value for estimating the prognosis of these patients.

PATIENTS AND METHODS

Depending on their survival time, 100 patients with primary lung SCC were separated into high- or low-risk prognostic groups, and their copy number aberrations (CNAs) were analyzed using array-comparative genomic hybridization (array-CGH).

RESULTS

We identified 123 CNA regions that were significantly associated with survival. Among these regions, some have been reported previously (eg, amplifications of 8p12, 3q27.1, and loss of 9p21.3 and 13q34) but others have never been reported. For example, gains of 3q27.1, 5p13.2, and 5p13.3 were found to be associated with a favorable prognosis, but patients harboring gains of 11q23.3, 11q13.1, and 14q32.3, and deletions of 3p21.3 and 9p21.3 tended to have poor survival. Among the 123 CNA regions, 41 were further selected to construct a survival estimation model that could effectively separate SCC patients into high- or low-risk groups with an accuracy of 92%, sensitivity of 90%, and specificity of 94%. The results of the array-CGH were further validated in an independent cohort of 45 formalin-fixed, paraffin-embedded specimens using real-time polymerase chain reaction.

CONCLUSION

A number of CNA regions were found to be associated with the survival of SCC patients, and we were able to construct a model to estimate prognosis on the basis of these regions. Assessment of these CNAs could potentially assist in clinical decision-making regarding adjuvant therapy after surgery.

Focus on 16p13.3 Locus in Colon Cancer

Background

With one million new cases of colorectal cancer (CRC) diagnosed annually in the world, CRC is the third most commonly diagnosed cancer in the Western world. Patients with stage I-III CRC can be cured with surgery but are at risk for recurrence. Colorectal cancer is characterized by the presence of chromosomal deletions and gains. Large genomic profiling studies have however not been conducted in this disease. The number of a specific genetic aberration in a tumour sample could correlate with recurrence-free survival or overall survival, possibly leading to its use as biomarker for therapeutic decisions. At this point there are not sufficient markers for prediction of disease recurrence in colorectal cancer, which can be used in the clinic to discriminate between stage II patients who will benefit from adjuvant chemotherapy. For instance, the benefit of adjuvant chemotherapy has been most clearly demonstrated in stage III disease with an approximately 30 percent relative reduction in the risk of disease recurrence. The benefits of adjuvant chemotherapy in stage II disease are less certain, the risk for relapse is much smaller in the overall group and the specific patients at risk are hard to identify.

Materials and Methods

In this study, array-comparative genomic hybridization analysis (array-CGH) was applied to study high-resolution DNA copy number alterations in 93 colon carcinoma samples. These genomic data were combined with parameters like KRAS mutation status, microsatellite status and clinicopathological characteristics.

Results

Both large and small chromosomal losses and gains were identified in our sample cohort. Recurrent gains were found for chromosome 1q, 7, 8q, 13 and 20 and losses were mostly found for 1p, 4, 8p, 14, 15, 17p, 18, 21 and 22. Data analysis demonstrated that loss of chromosome 4 is linked to a worse prognosis in our patients series. Besides these alterations, two interesting small regions of overlap were identified, which could be associated with disease recurrence. Gain of the 16p13.3 locus (including the RNA binding protein, fox-1 homolog gene, RBFOX1) was linked with a worse recurrence-free survival in our patient cohort. On the other hand, loss of RBFOX1 was only found in patients without disease recurrence. Most interestingly, above mentioned characteristics were also found in stage II patients, for whom there is a high medical need for the identification of new prognostic biomarkers.

Conclusions

In conclusion, copy number variation of the 16p13.3 locus seems to be an important parameter for prediction of disease recurrence in colon cancer.

Tumour Necrosis Factor-α Gene Polymorphism Is Associated with Metastasis in Patients with Triple Negative Breast Cancer

Tumour necrosis factor-α (TNF-α) is critical in the regulation of inflammation and tumour progression. TNF-α-308G > A is associated with constitutively elevated TNF-α expression. The purpose of this study was to assess the association between TNF-α-308G > A and breast cancer (BC) risk by subtype and the connection between genotypes and clinical features of BC. A total of 768 patients and 565 controls were enrolled in this study, and genotypes were detected using the TaqMan assay. No effect on susceptibility for any BC subtype was found for the TNF-α-308 polymorphism in our study or in the pooled meta-analysis. This polymorphism was shown to be associated with age at menarche in all BC and in progesterone receptor-negative BC. Interestingly, triple negative breast cancer (TNBC) patients with TNF-α-308A had an increased risk of distant tumour metastasis (OR = 3.80, 95% CI: 1.31-11.02, P = 0.009). Multi-regression analysis showed that TNF-α-308A was also a risk factor for distant tumour metastasis after adjustment for tumour size and lymph node metastasis status (OR = 6.26, 95% CI: 1.88-20.87, P = 0.003). These findings indicate that TNF-α might play a distinct role in the progression of TNBC, especially in distant tumour metastasis of TNBC.

Accurate computation of survival statistics in genome-wide studies

A key challenge in genomics is to identify genetic variants that distinguish patients with different survival time following diagnosis or treatment. While the log-rank test is widely used for this purpose, nearly all implementations of the log-rank test rely on an asymptotic approximation that is not appropriate in many genomics applications. This is because: the two populations determined by a genetic variant may have very different sizes; and the evaluation of many possible variants demands highly accurate computation of very small p-values. We demonstrate this problem for cancer genomics data where the standard log-rank test leads to many false positive associations between somatic mutations and survival time. We develop and analyze a novel algorithm, Exact Log-rank Test (ExaLT), that accurately computes the p-value of the log-rank statistic under an exact distribution that is appropriate for any size populations. We demonstrate the advantages of ExaLT on data from published cancer genomics studies, finding significant differences from the reported p-values. We analyze somatic mutations in six cancer types from The Cancer Genome Atlas (TCGA), finding mutations with known association to survival as well as several novel associations. In contrast, standard implementations of the log-rank test report dozens-hundreds of likely false positive associations as more significant than these known associations.

The SNP rs6500843 in 16p13.3 is associated with survival specifically among chemotherapy-treated breast cancer patients

We have utilized a two-stage study design to search for SNPs associated with the survival of breast cancer patients treated with adjuvant chemotherapy. Our initial GWS data set consisted of 805 Finnish breast cancer cases (360 treated with adjuvant chemotherapy). The top 39 SNPs from this stage were analyzed in three independent data sets: iCOGS (n=6720 chemotherapy-treated cases), SUCCESS-A (n=3596), and POSH (n=518). Two SNPs were successfully validated: rs6500843 (any chemotherapy; per-allele HR 1.16, 95% C.I. 1.08-1.26, p=0.0001, p(adjusted)=0.0091), and rs11155012 (anthracycline therapy; per-allele HR 1.21, 95% C.I. 1.08-1.35, p=0.0010, p(adjusted)=0.0270). The SNP rs6500843 was found to specifically interact with adjuvant chemotherapy, independently of standard prognostic markers (p(interaction)=0.0009), with the rs6500843-GG genotype corresponding to the highest hazard among chemotherapy-treated cases (HR 1.47, 95% C.I. 1.20-1.80). Upon trans-eQTL analysis of public microarray data, the rs6500843 locus was found to associate with the expression of a group of genes involved in cell cycle control, notably AURKA, the expression of which also exhibited differential prognostic value between chemotherapy-treated and untreated cases in our analysis of microarray data. Based on previously published information, we propose that the eQTL genes may be connected to the rs6500843 locus via a RBFOX1-FOXM1 -mediated regulatory pathway.

A genome wide association study on Newfoundland colorectal cancer patients' survival outcomes

Background

In this study we performed genome-wide association studies to identify candidate SNPs that may predict the risk of disease outcome in colorectal cancer.

Methods

Patient cohort consisted of 505 unrelated patients with Caucasian ancestry. Germline DNA samples were genotyped using the Illumina® human Omni-1quad SNP chip. Associations of SNPs with overall and disease free survivals were examined primarily for 431 patients with microsatellite instability-low (MSI-L) or stable (MSS) colorectal tumors using Cox proportional hazards method adjusting for clinical covariates. Bootstrap method was applied for internal validation of results. As exploratory analyses, association analyses for the colon (n = 334) and rectal (n = 171) cancer patients were also performed.

Results

As a result, there was no SNP that reached the genomewide significance levels (p < 5x10⁻⁸) in any of the analyses. A small number of genetic markers (n = 10) showed nominal associations (p <10⁻⁶) for MSS/MSI-L, colon, or rectal cancer patient groups. These markers were located in two non-coding RNA genes or intergenic regions and none were amino acid substituting polymorphisms. Bootstrap analysis for the MSS/MSI-L cohort data suggested the robustness of the observed nominal associations.

Conclusions

Likely due to small number of patients, our study did not identify an acceptable level of association of SNPs with outcome in MSS/MSI-L, colon, or rectal cancer patients. A number of SNPs with sub-optimal p-values were, however, identified; these loci may be promising and examined in other larger-sized patient cohorts.

Electronic supplementary material

The online version of this article (doi:10.1186/s40364-015-0031-6) contains supplementary material, which is available to authorized users.

Identification of novel vascular targets in lung cancer

BACKGROUND

Lung cancer remains the leading cause of cancer-related death, largely owing to the lack of effective treatments. A tumour vascular targeting strategy presents an attractive alternative; however, the molecular signature of the vasculature in lung cancer is poorly explored. This work aimed to identify novel tumour vascular targets in lung cancer.

METHODS

Enzymatic digestion of fresh tissue followed by endothelial capture with Ulex lectin-coated magnetic beads was used to isolate the endothelium from fresh tumour specimens of lung cancer patients. Endothelial isolates from the healthy and tumour lung tissue were subjected to whole human genome expression profiling using microarray technology.

RESULTS

Bioinformatics analysis identified tumour endothelial expression of angiogenic factors, matrix metalloproteases and cell-surface transmembrane proteins. Predicted novel tumour vascular targets were verified by RNA-seq, quantitative real-time PCR analysis and immunohistochemistry. Further detailed expression profiling of STEAP1 on 82 lung cancer patients confirmed STEAP1 as a novel target in the tumour vasculature. Functional analysis of STEAP1 using siRNA silencing implicates a role in endothelial cell migration and tube formation.

CONCLUSIONS

The identification of cell-surface tumour endothelial markers in lung is of interest in therapeutic antibody and vaccine development.

Interleukin polymorphisms associated with overall survival, disease-free survival, and recurrence in non-small cell lung cancer patients

Biomarkers based on germline DNA variations could have translational implications by identifying prognostic factors and sub-classifying patients to tailored, patient-specific treatment. To investigate the association between germline variations in interleukin (IL) genes and lung cancer outcomes, we genotyped 251 single nucleotide polymorphisms (SNPs) from 33 different IL genes in 651 non-small cell lung cancer (NSCLC) patients. Analyses were performed to investigate overall survival, disease-free survival, and recurrence. Our analyses revealed 24 different IL SNPs significantly associated with one or more of the lung cancer outcomes of interest. The GG genotype of IL16:rs7170924 was significantly associated with disease-free survival (HR = 0.65; 95% CI 0.50-0.83) and was the only SNP that produced a false discovery rate (FDR) of modest confidence that the association is unlikely to represent a false-positive result (FDR = 0.142). Classification and regression tree (CART) analyses were used to identify potential higher-order interactions. We restricted the CART analyses to the five SNPs that were significantly associated with multiple endpoints (IL1A:rs1800587, IL1B:rs1143634, IL8:s12506479, IL12A:rs662959, and IL13:rs1881457) and IL16:rs7170924 which had the lowest FDR. CART analyses did not yield a tree structure for overall survival; separate CART tree structures were identified for recurrence, based on three SNPs (IL13:rs1881457, IL1B:rs1143634, and IL12A:rs662959), and for disease-free survival, based on two SNPs (IL12A:rs662959 and IL16:rs7170924), which may suggest that these candidate IL SNPs have a specific impact on lung cancer progression and recurrence. These data suggest that germline variations in IL genes are associated with clinical outcomes in NSCLC patients.

A coding single-nucleotide polymorphism in lysine demethylase KDM4A associates with increased sensitivity to mTOR inhibitors

SNPs occur within chromatin-modulating factors; however, little is known about how these variants within the coding sequence affect cancer progression or treatment. Therefore, there is a need to establish their biochemical and/or molecular contribution, their use in subclassifying patients, and their impact on therapeutic response. In this report, we demonstrate that coding SNP-A482 within the lysine tridemethylase gene KDM4A/JMJD2A has different allelic frequencies across ethnic populations, associates with differential outcome in patients with non-small cell lung cancer (NSCLC), and promotes KDM4A protein turnover. Using an unbiased drug screen against 87 preclinical and clinical compounds, we demonstrate that homozygous SNP-A482 cells have increased mTOR inhibitor sensitivity. mTOR inhibitors significantly reduce SNP-A482 protein levels, which parallels the increased drug sensitivity observed with KDM4A depletion. Our data emphasize the importance of using variant status as candidate biomarkers and highlight the importance of studying SNPs in chromatin modifiers to achieve better targeted therapy.

SIGNIFICANCE

This report documents the first coding SNP within a lysine demethylase that associates with worse outcome in patients with NSCLC. We demonstrate that this coding SNP alters the protein turnover and associates with increased mTOR inhibitor sensitivity, which identifies a candidate biomarker for mTOR inhibitor therapy and a therapeutic target for combination therapy.