Effect of valproic acid on histone deacetylase expression in oral cancer (Review)

Oral squamous cell carcinoma (OSCC) is a frequent human malignancy that demonstrates a range of genetic and epigenetic alterations. Histone deacetylases (HDACs) are key epigenetic regulators of cell-cycle progression, differentiation and apoptosis and their dysregulation is implicated in cancer development. HDACs are promising targets for anticancer therapy through the utilisation of HDAC inhibitors (HDACis). OSCC cells have been shown to have low levels of histone acetylation, suggesting that HDACis may produce beneficial effects in patients with OSCC. Valproic acid (VPA) is a class I and IIa HDACi and, therefore, may be useful in anticancer therapy. VPA has been reported as a chemo-preventive epigenetic agent in individuals with high-risk oral dysplasia (OD) and thus associated with a reduced risk of HNSCC. It is hypothesised that HDAC inhibition by VPA triggers a change in the expression levels of different HDAC family gene-members. The present review summarises the current literature on HDAC expression changes in response to VPA in oral cancer patients and in vitro studies in an effort to better understand the potential epigenetic impact of VPA treatment. The present review outlined the need for exploring supportive evidence of the chemo-preventive role played by VPA-based epigenetic modification in treating oral pre-cancerous lesions and, thus, providing a novel tolerable chemotherapeutic strategy for patients with oral cancer.

KH-like Domains in PARP9/DTX3L and PARP14 Coordinate Protein-Protein Interactions to Promote Cancer Cell Survival

Certain members of the ADP-ribosyltransferase superfamily (ARTD or PARP enzymes) catalyse ADP-ribosylation in response to cellular stress, DNA damage and viral infection and are upregulated in various tumours. PARP9, its binding partner DTX3L and PARP14 protein levels are significantly correlated in head and neck squamous cell carcinoma (HNSCC) and other tumour types though a mechanism where PARP9/DTX3L regulates PARP14 post-transcriptionally. Depleting PARP9, DTX3L or PARP14 expression in HNSCC or HeLa cell lines decreases cell survival through a reduction of proliferation and an increase in apoptosis. A partial rescue of survival was achieved by expressing a PARP14 truncation containing a predicted eukaryotic type I KH domain. KH-like domains were also found in PARP9 and in DTX3L and contributed to protein-protein interactions between PARP9-DTX3L and PARP14-DTX3L. Homodimerization of DTX3L was also coordinated by a KH-like domain and was disrupted by site-specific mutation. Although, cell survival promoted by PARP14 did not require ADP-ribosyltransferase activity, interaction of DTX3L in vitro suppressed PARP14 auto-ADP-ribosylation and promoted trans-ADP-ribosylation of PARP9 and DTX3L. In summary, we characterised PARP9-DTX3L-PARP14 interactions important to pro-survival signalling in HNSCC cells, albeit in PARP14 catalytically independent fashion.

Non-coding RNA-directed therapeutics in lung cancer: Delivery technologies and clinical applications

Lung cancer is one of the most aggressive and deadliest health threats. There has been an increasing interest in non-coding RNA (ncRNA) recently, especially in the areas of carcinogenesis and tumour progression. However, ncRNA-directed therapies are still encountering obstacles on their way to the clinic. In the present article, we provide an overview on the potential of targeting ncRNA in the treatment of lung cancer. Then, we discuss the delivery challenges and recent approaches enabling the delivery of ncRNA-directed therapies to the lung cancer cells, where we illuminate some advanced technologies including chemically-modified oligonucleotides, nuclear targeting, and three-dimensional in vitro models. Furthermore, advanced non-viral delivery systems recruiting nanoparticles, biomimetic delivery systems, and extracellular vesicles are also highlighted. Lastly, the challenges limiting the clinical trials on the therapeutic targeting of ncRNAs in lung cancer and future directions to tackle them are explored.

Plasma protein biomarkers for early prediction of lung cancer

BACKGROUND

Individual plasma proteins have been identified as minimally invasive biomarkers for lung cancer diagnosis with potential utility in early detection. Plasma proteomes provide insight into contributing biological factors; we investigated their potential for future lung cancer prediction.

METHODS

The Olink® Explore-3072 platform quantitated 2941 proteins in 496 Liverpool Lung Project plasma samples, including 131 cases taken 1-10 years prior to diagnosis, 237 controls, and 90 subjects at multiple times. 1112 proteins significantly associated with haemolysis were excluded. Feature selection with bootstrapping identified differentially expressed proteins, subsequently modelled for lung cancer prediction and validated in UK Biobank data.

FINDINGS

For samples 1-3 years pre-diagnosis, 240 proteins were significantly different in cases; for 1-5 year samples, 117 of these and 150 further proteins were identified, mapping to significantly different pathways. Four machine learning algorithms gave median AUCs of 0.76-0.90 and 0.73-0.83 for the 1-3 year and 1-5 year proteins respectively. External validation gave AUCs of 0.75 (1-3 year) and 0.69 (1-5 year), with AUC 0.7 up to 12 years prior to diagnosis. The models were independent of age, smoking duration, cancer histology and the presence of COPD.

INTERPRETATION

The plasma proteome provides biomarkers which may be used to identify those at greatest risk of lung cancer. The proteins and the pathways are different when lung cancer is more imminent, indicating that both biomarkers of inherent risk and biomarkers associated with presence of early lung cancer may be identified.

FUNDING

Janssen Pharmaceuticals Research Collaboration Award; Roy Castle Lung Cancer Foundation.

Variable allelic expression of imprinted genes at the Peg13, Trappc9, Ago2 cluster in single neural cells

Genomic imprinting is an epigenetic process through which genes are expressed in a parent-of-origin specific manner resulting in mono-allelic or strongly biased expression of one allele. For some genes, imprinted expression may be tissue-specific and reliant on CTCF-influenced enhancer-promoter interactions. The Peg13 imprinting cluster is associated with neurodevelopmental disorders and comprises canonical imprinted genes, which are conserved between mouse and human, as well as brain-specific imprinted genes in mouse. The latter consist of Trappc9, Chrac1 and Ago2, which have a maternal allelic expression bias of ∼75% in brain. Findings of such allelic expression biases on the tissue level raise the question of how they are reflected in individual cells and whether there is variability and mosaicism in allelic expression between individual cells of the tissue. Here we show that Trappc9 and Ago2 are not imprinted in hippocampus-derived neural stem cells (neurospheres), while Peg13 retains its strong bias of paternal allele expression. Upon analysis of single neural stem cells and in vitro differentiated neurons, we find not uniform, but variable states of allelic expression, especially for Trappc9 and Ago2. These ranged from mono-allelic paternal to equal bi-allelic to mono-allelic maternal, including biased bi-allelic transcriptional states. Even Peg13 expression deviated from its expected paternal allele bias in a small number of cells. Although the cell populations consisted of a mosaic of cells with different allelic expression states, as a whole they reflected bulk tissue data. Furthermore, in an attempt to identify potential brain-specific regulatory elements across the Trappc9 locus, we demonstrate tissue-specific and general silencer activities, which might contribute to the regulation of its imprinted expression bias.

Gene-gene interaction of AhRwith and within the Wntcascade affects susceptibility to lung cancer

BACKGROUND

Aberrant Wnt signalling, regulating cell development and stemness, influences the development of many cancer types. The Aryl hydrocarbon receptor (AhR) mediates tumorigenesis of environmental pollutants. Complex interaction patterns of genes assigned to AhR/Wnt-signalling were recently associated with lung cancer susceptibility.

AIM

To assess the association and predictive ability of AhR/Wnt-genes with lung cancer in cases and controls of European descent.

METHODS

Odds ratios (OR) were estimated for genomic variants assigned to the Wnt agonist and the antagonistic genes DKK2, DKK3, DKK4, FRZB, SFRP4 and Axin2. Logistic regression models with variable selection were trained, validated and tested to predict lung cancer, at which other previously identified SNPs that have been robustly associated with lung cancer risk could also enter the model. Furthermore, decision trees were created to investigate variant × variant interaction. All analyses were performed for overall lung cancer and for subgroups.

RESULTS

No genome-wide significant association of AhR/Wnt-genes with overall lung cancer was observed, but within the subgroups of ever smokers (e.g., maker rs2722278 SFRP4; OR  = 1.20; 95% CI 1.13-1.27; p  = 5.6 × 10-10) and never smokers (e.g., maker rs1133683 Axin2; OR  = 1.27; 95% CI 1.19-1.35; p  = 1.0 × 10-12). Although predictability is poor, AhR/Wnt-variants are unexpectedly overrepresented in optimized prediction scores for overall lung cancer and for small cell lung cancer. Remarkably, the score for never-smokers contained solely two AhR/Wnt-variants. The optimal decision tree for never smokers consists of 7 AhR/Wnt-variants and only two lung cancer variants.

CONCLUSIONS

The role of variants belonging to Wnt/AhR-pathways in lung cancer susceptibility may be underrated in main-effects association analysis. Complex interaction patterns in individuals of European descent have moderate predictive capacity for lung cancer or subgroups thereof, especially in never smokers.

Iam hiQ-a novel pair of accuracy indices for imputed genotypes

BACKGROUND

Imputation of untyped markers is a standard tool in genome-wide association studies to close the gap between directly genotyped and other known DNA variants. However, high accuracy with which genotypes are imputed is fundamental. Several accuracy measures have been proposed and some are implemented in imputation software, unfortunately diversely across platforms. In the present paper, we introduce Iam hiQ, an independent pair of accuracy measures that can be applied to dosage files, the output of all imputation software. Iam (imputation accuracy measure) quantifies the average amount of individual-specific versus population-specific genotype information in a linear manner. hiQ (heterogeneity in quantities of dosages) addresses the inter-individual heterogeneity between dosages of a marker across the sample at hand.

RESULTS

Applying both measures to a large case-control sample of the International Lung Cancer Consortium (ILCCO), comprising 27,065 individuals, we found meaningful thresholds for Iam and hiQ suitable to classify markers of poor accuracy. We demonstrate how Manhattan-like plots and moving averages of Iam and hiQ can be useful to identify regions enriched with less accurate imputed markers, whereas these regions would by missed when applying the accuracy measure info (implemented in IMPUTE2).

CONCLUSION

We recommend using Iam hiQ additional to other accuracy scores for variant filtering before stepping into the analysis of imputed GWAS data.

Genetic variants associated with mandibular osteoradionecrosis following radiotherapy for head and neck malignancy

BACKGROUND/AIM

Utilising radiotherapy in the management of head and neck cancer (HNC) often results in long term toxicities. Mandibular osteoradionecrosis (ORN) represents a late toxicity associated with significant morbidity. We aim to identify a panel of common genetic variants which can predict ORN to aid development of personalised radiotherapy protocols.

METHOD

Single nucleotide polymorphism (SNP) arrays were applied to DNA samples from patients who had prior HNC radiotherapy and minimum two years follow-up. A case cohort of mandibular ORN was compared to a control group of participants recruited to CRUK HOPON clinical trial. Relevant clinical parameters influencing ORN risk (e.g. smoking/alcohol) were collected. Significant associations from array data were internally validated using polymerase chain reaction (PCR) and pyrosequencing.

RESULTS

Following inclusion of 141 patients in the analysis (52 cases, 89 controls), a model predictive for ORN was developed; after controlling for alcohol consumption, smoking, and age, 4053 SNPs were identified as significant. This was reduced to a representative model of 18 SNPs achieving 92% accuracy. Following internal technical validation, a six SNP model (rs34798038, rs6011731, rs2348569, rs530752, rs7477958, rs1415848) was retained within multivariate regression analysis (ROC AUC 0.859). Of these, four SNPs (rs34798038 (A/G) (p 0.006), rs6011731 (C/T) (p 0.018), rs530752 (A/G) (p 0.046) and rs2348569 (G/G) (p 0.005)) were significantly associated with the absence of ORN.

CONCLUSION

This is the first genome wide association study in HNC using ORN as the endpoint and offers new insight into ORN pathogenesis. Subject to validation, these variants may guide patient selection for personalised radiotherapy strategies.

SAVER: sodium valproate for the epigenetic reprogramming of high-risk oral epithelial dysplasia-a phase II randomised control trial study protocol

Background

Sodium valproate (VPA) has been associated with a reduced risk of head and neck cancer development. The potential protective mechanism of action is believed to be via inhibition of histone deacetylase and subsequent epigenetic reprogramming. SAVER is a phase IIb open-label, randomised control trial of VPA as a chemopreventive agent in patients with high-risk oral epithelial dysplasia (OED). The aim of the trial is to gather preliminary evidence of the clinical and biological effects of VPA upon OED and assess the feasibility and acceptability of such a trial, with a view to inform a future definitive phase III study.

Methods

One hundred and ten patients with high-risk OED will be recruited from up to 10 secondary care sites in the UK and randomised into either VPA or observation only for 4 months. Women of childbearing potential will be excluded due to the teratogenic properties of VPA. Tissue and blood samples will be collected prior to randomisation and on the last day of the intervention/observation-only period (end of 4 months). Clinical measurement and additional safety bloods will be taken at multiple time points during the trial. The primary outcome will be a composite, surrogate endpoint of change in lesion size, change in grade of dysplasia and change in LOH profile at 8 key microsatellite regions. Feasibility outcomes will include recruitment targets, compliance with the study protocol and adverse effects. A qualitative sub-study will explore patient experience and perception of the trial.

Discussion

The current management options for patients with high-risk OED are limited and mostly include surgical resection and clinical surveillance. However, there remains little evidence whether surgery can effectively lead to a notable reduction in the risk of oral cancer development. Similarly, surveillance is associated with concerns regarding delayed diagnosis of OED progressing to malignancy. The SAVER trial provides an opportunity to investigate the effects of a repurposed, inexpensive and well-tolerated medication as a potential chemopreventive strategy for patients with high-risk OED. The clinical and biological findings of SAVER will inform the appropriateness, design and feasibility of a definitive phase III trial.

Trial registration

The trial is registered with the European Clinical Trials Database (Eudra-CT 2018-000197-30). (http://www.isrctn.com/ISRCTN12448611). The trial was prospectively registered on 24/04/2018.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13063-021-05373-8.

Cyclic AMP Response Element Modulator-α Suppresses PD-1 Expression and Promotes Effector CD4+ T Cells in Psoriasis

Effector CD4⁺ T lymphocytes contribute to inflammation and tissue damage in psoriasis, but the underlying molecular mechanisms remain poorly understood. The transcription factor CREMα controls effector T cell function in people with systemic autoimmune diseases. The inhibitory surface coreceptor PD-1 plays a key role in the control of effector T cell function and its therapeutic inhibition in patients with cancer can cause psoriasis. In this study, we show that CD4⁺ T cells from patients with psoriasis and psoriatic arthritis exhibit increased production of IL-17 but decreased expression of IL-2 and PD-1. In genetically modified mice and Jurkat T cells CREMα expression was linked to low PD-1 levels. We demonstrate that CREMα is recruited to the proximal promoter of PDCD1 in which it trans-represses gene expression and corecruits DNMT3a-mediating DNA methylation. As keratinocytes limit inflammation by PD-1 ligand expression and, in this study, reported reduced expression of PD-1 on CD4⁺ T cells is linked to low IL-2 and high IL-17A production, our studies reveal a molecular pathway in T cells from people with psoriasis that can deserve clinical exploitation.

TLR9 expression in chronic lymphocytic leukemia identifies a promigratory subpopulation and novel therapeutic target

Chronic lymphocytic leukemia (CLL) remains incurable despite B-cell receptor-targeted inhibitors revolutionizing treatment. This suggests that other signaling molecules are involved in disease escape mechanisms and resistance. Toll-like receptor 9 (TLR9) is a promising candidate that is activated by unmethylated cytosine guanine dinucleotide-DNA. Here, we show that plasma from patients with CLL contains significantly more unmethylated DNA than plasma from healthy control subjects (P < .0001) and that cell-free DNA levels correlate with the prognostic markers CD38, β2-microglobulin, and lymphocyte doubling time. Furthermore, elevated cell-free DNA was associated with shorter time to first treatment (hazard ratio, 4.0; P = .003). We also show that TLR9 expression was associated with in vitro CLL cell migration (P < .001), and intracellular endosomal TLR9 strongly correlated with aberrant surface expression (sTLR9; r = 0.9). In addition, lymph node-derived CLL cells exhibited increased sTLR9 (P = .016), and RNA-sequencing of paired sTLR9hi and sTLR9lo CLL cells revealed differential transcription of genes involved in TLR signaling, adhesion, motility, and inflammation in sTLR9hi cells. Mechanistically, a TLR9 agonist, ODN2006, promoted CLL cell migration (P < .001) that was mediated by p65 NF-κB and STAT3 transcription factor activation. Importantly, autologous plasma induced the same effects, which were reversed by a TLR9 antagonist. Furthermore, high TLR9 expression promoted engraftment and rapid disease progression in a NOD/Shi-scid/IL-2Rγnull mouse xenograft model. Finally, we showed that dual targeting of TLR9 and Bruton's tyrosine kinase (BTK) was strongly synergistic (median combination index, 0.2 at half maximal effective dose), which highlights the distinct role for TLR9 signaling in CLL and the potential for combined targeting of TLR9 and BTK as a more effective treatment strategy in this incurable disease.

A long intergenic non-coding RNA regulates nuclear localization of DNA methyl transferase-1

DNA methyl transferase-1 or DNMT1 maintains DNA methylation in the genome and is important for regulating gene expression in cells. Aberrant changes in DNMT1 activity and DNA methylation are commonly observed in cancers and many other diseases. Recently, a number of long intergenic non-protein-coding RNAs or lincRNAs have been shown to play a role in regulating DNMT1 activity. CCDC26 is a nuclear lincRNA that is frequently mutated in cancers and is a hotbed for disease-associated single nucleotide changes. However, the functional mechanism of CCDC26 is not understood. Here, we show that this lincRNA is concentrated on the nuclear periphery. Strikingly, in the absence of CCDC26 lincRNA, DNMT1 is mis-located in the cytoplasm, and the genomic DNA is significantly hypomethylated. This is accompanied by double-stranded DNA breaks and increased cell death. These results point to a previously unrecognized mechanism of lincRNA-mediated subcellular localization of DNMT1 and regulation of DNA methylation.

Rare deleterious germline variants and risk of lung cancer

Recent studies suggest that rare variants exhibit stronger effect sizes and might play a crucial role in the etiology of lung cancers (LC). Whole exome plus targeted sequencing of germline DNA was performed on 1045 LC cases and 885 controls in the discovery set. To unveil the inherited causal variants, we focused on rare and predicted deleterious variants and small indels enriched in cases or controls. Promising candidates were further validated in a series of 26,803 LCs and 555,107 controls. During discovery, we identified 25 rare deleterious variants associated with LC susceptibility, including 13 reported in ClinVar. Of the five validated candidates, we discovered two pathogenic variants in known LC susceptibility loci, ATM p.V2716A (Odds Ratio [OR] 19.55, 95%CI 5.04-75.6) and MPZL2 p.I24M frameshift deletion (OR 3.88, 95%CI 1.71-8.8); and three in novel LC susceptibility genes, POMC c.*28delT at 3' UTR (OR 4.33, 95%CI 2.03-9.24), STAU2 p.N364M frameshift deletion (OR 4.48, 95%CI 1.73-11.55), and MLNR p.Q334V frameshift deletion (OR 2.69, 95%CI 1.33-5.43). The potential cancer-promoting role of selected candidate genes and variants was further supported by endogenous DNA damage assays. Our analyses led to the identification of new rare deleterious variants with LC susceptibility. However, in-depth mechanistic studies are still needed to evaluate the pathogenic effects of these specific alleles.

Activating transcription factor-2 (ATF2) is a key determinant of resistance to endocrine treatment in an in vitro model of breast cancer

BACKGROUND

Activating transcription factor-2 (ATF2), a member of the leucine zipper family of DNA binding proteins, has been implicated as a tumour suppressor in breast cancer. However, its exact role in breast cancer endocrine resistance is still unclear. We have previously shown that silencing of ATF2 leads to a loss in the growth-inhibitory effects of tamoxifen in the oestrogen receptor (ER)-positive, tamoxifen-sensitive MCF7 cell line and highlighted that this multi-faceted transcription factor is key to the effects of tamoxifen in an endocrine sensitive model. In this work, we explored further the in vitro role of ATF2 in defining the resistance to endocrine treatment.

MATERIALS AND METHODS

We knocked down ATF2 in TAMR, LCC2 and LCC9 tamoxifen-resistant breast cancer cell lines as well as the parental tamoxifen sensitive MCF7 cell line and investigated the effects on growth, colony formation and cell migration. We also performed a microarray gene expression profiling (Illumina Human HT12_v4) to explore alterations in gene expression between MCF7 and TAMRs after ATF2 silencing and confirmed gene expression changes by quantitative RT-PCR.

RESULTS

By silencing ATF2, we observed a significant growth reduction of TAMR, LCC2 and LCC9 with no such effect observed with the parental MCF7 cells. ATF2 silencing was also associated with a significant inhibition of TAMR, LCC2 and LCC9 cell migration and colony formation. Interestingly, knockdown of ATF2 enhanced the levels of ER and ER-regulated genes, TFF1, GREB1, NCOA3 and PGR, in TAMR cells both at RNA and protein levels. Microarray gene expression identified a number of genes known to mediate tamoxifen resistance, to be differentially regulated by ATF2 in TAMR in relation to the parental MCF7 cells. Moreover, differential pathway analysis confirmed enhanced ER activity after ATF2 knockdown in TAMR cells.

CONCLUSION

These data demonstrate that ATF2 silencing may overcome endocrine resistance and highlights further the dual role of this transcription factor that can mediate endocrine sensitivity and resistance by modulating ER expression and activity.

Low protein intake during reproduction compromises the recovery of lactation-induced bone loss in female mouse dams without affecting skeletal muscles

Lactation-induced bone loss occurs due to high calcium requirements for fetal growth but skeletal recovery is normally achieved promptly postweaning. Dietary protein is vital for fetus and mother but the effects of protein undernutrition on the maternal skeleton and skeletal muscles are largely unknown. We used mouse dams fed with normal (N, 20%) or low (L, 8%) protein diet during gestation and lactation and maintained on the same diets (NN, LL) or switched from low to normal (LN) during a 28 d skeletal restoration period post lactation. Skeletal muscle morphology and neuromuscular junction integrity was not different between any of the groups. However, dams fed the low protein diet showed extensive bone loss by the end of lactation, followed by full skeletal recovery in NN dams, partial recovery in LN and poor bone recovery in LL dams. Primary osteoblasts from low protein diet fed mice showed decreased in vitro bone formation and decreased osteogenic marker gene expression; promoter methylation analysis by pyrosequencing showed no differences in Bmpr1a, Ptch1, Sirt1, Osx, and Igf1r osteoregulators, while miR-26a, -34a, and -125b expression was found altered in low protein fed mice. Therefore, normal protein diet is indispensable for maternal musculoskeletal health during the reproductive period.

Fibroblasts from Distinct Pancreatic Pathologies Exhibit Disease-Specific Properties

Although fibrotic stroma forms an integral component of pancreatic diseases, whether fibroblasts programmed by different types of pancreatic diseases are phenotypically distinct remains unknown. Here, we show that fibroblasts isolated from patients with pancreatic ductal adenocarcinoma (PDAC), chronic pancreatitis (CP), periampullary tumors, and adjacent normal (NA) tissue (N = 34) have distinct mRNA and miRNA profiles. Compared with NA fibroblasts, PDAC-associated fibroblasts were generally less sensitive to an antifibrotic stimulus (NPPB) and more responsive to positive regulators of activation such as TGFβ1 and WNT. Of the disease-associated fibroblasts examined, PDAC- and CP-derived fibroblasts shared greatest similarity, yet PDAC-associated fibroblasts expressed higher levels of tenascin C (TNC), a finding attributable to miR-137, a novel regulator of TNC. TNC protein and transcript levels were higher in PDAC tissue versus CP tissue and were associated with greater levels of stromal activation, and conditioned media from TNC-depleted PDAC-associated fibroblasts modestly increased both PDAC cell proliferation and PDAC cell migration, indicating that stromal TNC may have inhibitory effects on PDAC cells. Finally, circulating TNC levels were higher in patients with PDAC compared with CP. Our characterization of pancreatic fibroblast programming as disease-specific has consequences for therapeutic targeting and for the manner in which fibroblasts are used in research. SIGNIFICANCE: Primary fibroblasts derived from various types of pancreatic diseases possess and retain distinct molecular and functional characteristics in culture, providing a series of cellular models for treatment development and disease-specific research.

Overexpression of HURP mRNA in head and neck carcinoma and association with in vitro response to vinorelbine

HURP gene encodes the hepatoma upregulated protein (HURP), a microtubule associated protein regulating mitotic spindle dynamics, which promotes chromosomal congression and alignment during mitosis, with a potential role in tumorigenesis. In the present study, HURP mRNA expression was investigated by reverse transcription-quantitative PCR in oropharyngeal squamous cell carcinoma (OPSCC). Primary OPSCC tumors from 107 patients and 48 adjacent normal tissues, as well as 12 respiratory tract cancer cell lines (9 head and neck squamous cell carcinoma, 2 lung cancer and 1 normal bronchial) were utilised in the present study. mRNA expression levels of HURP were higher in malignant OPSCC tissues compared with in normal mucosa (P<1×10^-5) and significantly associated with sex and smoking status (P<0.0001). Vinorelbine in vitro toxicity at half-maximal inhibitory concentration (IC₅₀) was measured in the 11 cancer cell lines using an MTT assay. Sensitivity to vinorelbine was significantly correlated with HURP expression (r=0.636; P=0.035). The data indicated that HURP overexpression is frequent in OPSCC tissues and associated with smoking. The correlation between HURP mRNA expression and vinorelbine in vitro response suggests that HURP is a potential modulator of vinorelbine response; therefore, it should be explored for its possible predictive value for the efficiency of vinorelbine treatment in this type of cancer.

Integrated omics profiling reveals novel patterns of epigenetic programming in cancer-associated myofibroblasts

There is increasing evidence that stromal myofibroblasts play a key role in the tumour development however, the mechanisms by which they become reprogrammed to assist in cancer progression remain unclear. As cultured cancer-associated myofibroblasts (CAMs) retain an ability to enhance the proliferation and migration of cancer cells in vitro, it is possible that epigenetic reprogramming of CAMs within the tumour microenvironment may confer long-term pro-tumourigenic changes in gene expression. This study reports the first comparative multi-omics analysis of cancer-related changes in gene expression and DNA methylation in primary myofibroblasts derived from gastric and oesophageal tumours. In addition, we identify novel CAM-specific DNA methylation signatures, which are not observed in patient-matched adjacent tissue-derived myofibroblasts, or corresponding normal tissue-derived myofibroblasts. Analysis of correlated changes in DNA methylation and gene expression shows that different patterns of gene-specific DNA methylation have the potential to confer pro-tumourigenic changes in metabolism, cell signalling and differential responses to hypoxia. These molecular signatures provide new insights into potential mechanisms of stromal reprogramming in gastric and oesophageal cancer, while also providing a new resource to facilitate biomarker identification and future hypothesis-driven studies into mechanisms of stromal reprogramming and tumour progression in solid tumours.

Long non-coding RNA dysregulation is a frequent event in non-small cell lung carcinoma pathogenesis

Background

Long non-coding RNAs compose an important level of epigenetic regulation in normal physiology and disease. Despite the plethora of publications of lncRNAs in human cancer, the landscape is still unclear.

Methods

Microarray analysis in 44 NSCLC paired specimens was followed by qPCR-based validation in 29 (technical) and 38 (independent) tissue pairs. Cross-validation of the selected targets was achieved in 850 NSCLC tumours from TCGA datasets.

Results

Twelve targets were successfully validated by qPCR (upregulated: FEZF1-AS1, LINC01214, LINC00673, PCAT6, NUTM2A-AS1, LINC01929; downregulated: PCAT19, FENDRR, SVIL-AS1, LANCL1-AS1, ADAMTS9-AS2 and LINC00968). All of them were successfully cross validated in the TCGA datasets. Abnormal DNA methylation was observed in the promoters of FENDRR, FEZF1-AS1 and SVIL-AS1. FEZF1-AS1 and LINC01929 were associated with survival in the TCGA set.

Conclusions

Our study provides through multiple levels of internal and external validation, a comprehensive list of dysregulated lncRNAs in NSCLC. We therefore envisage this dataset to serve as an important source for the lung cancer research community assisting future investigations on the involvement of lncRNAs in the pathogenesis of the disease and providing novel biomarkers for diagnosis, prognosis and therapeutic stratification.

Publisher Correction: Shared heritability and functional enrichment across six solid cancers

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Lung Cancer Risk in Never-Smokers of European Descent is Associated With Genetic Variation in the 5p15.33 TERT-CLPTM1Ll Region

INTRODUCTION

Inherited susceptibility to lung cancer risk in never-smokers is poorly understood. The major reason for this gap in knowledge is that this disease is relatively uncommon (except in Asians), making it difficult to assemble an adequate study sample. In this study we conducted a genome-wide association study on the largest, to date, set of European-descent never-smokers with lung cancer.

METHODS

We conducted a two-phase (discovery and replication) genome-wide association study in never-smokers of European descent. We further augmented the sample by performing a meta-analysis with never-smokers from the recent OncoArray study, which resulted in a total of 3636 cases and 6295 controls. We also compare our findings with those in smokers with lung cancer.

RESULTS

We detected three genome-wide statistically significant single nucleotide polymorphisms rs31490 (odds ratio [OR]: 0.769, 95% confidence interval [CI]: 0.722-0.820; p value 5.31 × 10-16), rs380286 (OR: 0.770, 95% CI: 0.723-0.820; p value 4.32 × 10-16), and rs4975616 (OR: 0.778, 95% CI: 0.730-0.829; p value 1.04 × 10-14). All three mapped to Chromosome 5 CLPTM1L-TERT region, previously shown to be associated with lung cancer risk in smokers and in never-smoker Asian women, and risk of other cancers including breast, ovarian, colorectal, and prostate.

CONCLUSIONS

We found that genetic susceptibility to lung cancer in never-smokers is associated to genetic variants with pan-cancer risk effects. The comparison with smokers shows that top variants previously shown to be associated with lung cancer risk only confer risk in the presence of tobacco exposure, underscoring the importance of gene-environment interactions in the etiology of this disease.

Donor-Dependent and Other Nondefined Factors Have Greater Influence on the Hepatic Phenotype Than the Starting Cell Type in Induced Pluripotent Stem Cell Derived Hepatocyte-Like Cells

Drug‐induced liver injury is the greatest cause of post‐marketing drug withdrawal; therefore, substantial resources are directed toward triaging potentially dangerous new compounds at all stages of drug development. One of the major factors preventing effective screening of new compounds is the lack of a predictive in vitro model of hepatotoxicity. Primary human hepatocytes offer a metabolically relevant model for which the molecular initiating events of hepatotoxicity can be examined; however, these cells vary greatly between donors and dedifferentiate rapidly in culture. Induced pluripotent stem cell (iPSC)‐derived hepatocyte‐like cells (HLCs) offer a reproducible, physiologically relevant and genotypically normal model cell; however, current differentiation protocols produce HLCs with a relatively immature phenotype. During the reprogramming of somatic cells, the epigenome undergoes dramatic changes; however, this “resetting” is a gradual process, resulting in an altered differentiation propensity, skewed toward the lineage of origin, particularly in early passage cultures. We, therefore, performed a comparison of human hepatocyte‐ and dermal fibroblast‐derived iPSCs, assessing the impact of epigenetic memory at all stages of HLC differentiation. These results provide the first isogenic assessment of the starting cell type in human iPSC‐derived HLCs. Despite a trend toward improvement in hepatic phenotype in albumin secretion and gene expression, few significant differences in hepatic differentiation capacity were found between hepatocyte and fibroblast‐derived iPSCs. We conclude that the donor and inter‐clonal differences have a greater influence on the hepatocyte phenotypic maturity than the starting cell type. Therefore, it is not necessary to use human hepatocytes for generating iPSC‐derived HLCs. Stem Cells Translational Medicine2017;6:1321–1331

The expression signature of cancer-associated KRAB-ZNF factors identified in TCGA pan-cancer transcriptomic data

The KRAB-ZNF (Krüppel-associated box domain zinc finger) gene family is composed of a large number of highly homologous genes, gene isoforms, and pseudogenes. The proteins encoded by these genes, whose expression is often tissue-specific, act as epigenetic suppressors contributing to the addition of repressive chromatin marks and DNA methylation. Due to its high complexity, the KRAB-ZNF family has not been studied in sufficient detail, and the involvement of its members in carcinogenesis remains mostly unexplored. In this study, we aimed to provide a comprehensive description of cancer-associated KRAB-ZNFs using publicly available The Cancer Genome Atlas pan-cancer datasets. We analyzed 6727 tumor and normal tissue samples from 16 cancer types. Here, we showed that a small but distinctive cluster of 16 KRAB-ZNFs is commonly upregulated across multiple cancer cohorts in comparison to normal samples. We confirmed these observations in the independent panels of lung and breast cancer cell lines and tissues. This upregulation was also observed for most of the KRAB-ZNF splicing variants, whose expression is simultaneously upregulated in tumors compared to normal tissues. Finally, by analyzing the clinicopathological data for breast and lung cancers, we demonstrated that the expression of cancer-associated KRAB-ZNFs correlates with patient survival, tumor histology, and molecular subtyping. Altogether, our study allowed the identification and characterization of KRAB-ZNF factors that may have an essential function in cancer biology and thus potential to become novel oncologic biomarkers and treatment targets.

Shared heritability and functional enrichment across six solid cancers

Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (rg = 0.57, p = 4.6 × 10-8), breast and ovarian cancer (rg = 0.24, p = 7 × 10-5), breast and lung cancer (rg = 0.18, p =1.5 × 10-6) and breast and colorectal cancer (rg = 0.15, p = 1.1 × 10-4). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.

Whole-genome methylation profiling of the retinal pigment epithelium of individuals with age-related macular degeneration reveals differential methylation of the SKI, GTF2H4, and TNXB genes

Background

Age-related macular degeneration (AMD) is a degenerative disorder of the central retina and the foremost cause of blindness. The retinal pigment epithelium (RPE) is a primary site of disease pathogenesis. The genetic basis of AMD is relatively well understood; however, this knowledge is yet to yield a treatment for the most prevalent non-neovascular disease forms. Therefore, tissue-specific epigenetic mechanisms of gene regulation are of considerable interest in AMD. We aimed to identify differentially methylated genes associated with AMD in the RPE and differentiate local DNA methylation aberrations from global DNA methylation changes, as local DNA methylation changes may be more amenable to therapeutic manipulation.

Methods

Epigenome-wide association study and targeted gene expression profiling were carried out in RPE cells from eyes of human donors. We performed genome-wide DNA methylation profiling (Illumina 450k BeadChip array) on RPE cells from 44 human donor eyes (25 AMD and 19 normal controls). We validated the findings using bisulfite pyrosequencing in 55 RPE samples (30 AMD and 25 normal controls) including technical (n = 38) and independent replicate samples (n = 17). Long interspersed nucleotide element 1 (LINE-1) analysis was then applied to assess global DNA methylation changes in the RPE. RT-qPCR on independent donor RPE samples was performed to assess gene expression changes.

Results

Genome-wide DNA methylation profiling identified differential methylation of multiple loci including the SKI proto-oncogene (SKI) (p = 1.18 × 10⁻⁹), general transcription factor IIH subunit H4 (GTF2H4) (p = 7.03 × 10⁻⁷), and Tenascin X (TNXB) (p = 6.30 × 10⁻⁶) genes in AMD. Bisulfite pyrosequencing validated the differentially methylated locus cg18934822 in SKI, and cg22508626 within GTF2H4, and excluded global DNA methylation changes in the RPE in AMD. We further demonstrated the differential expression of SKI, GTF2H4, and TNXB in the RPE of independent AMD donors.

Conclusions

We report the largest genome-wide methylation analysis of RPE in AMD along with associated gene expression changes to date, for the first-time reaching genome-wide significance, and identified novel targets for functional and future therapeutic intervention studies. The novel differentially methylated genes SKI and GTF2H4 have not been previously associated with AMD, and regulate disease pathways implicated in AMD, including TGF beta signaling (SKI) and transcription-dependent DNA repair mechanisms (GTF2H4).

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0608-2) contains supplementary material, which is available to authorized users.

Rare Variants in Known Susceptibility Loci and Their Contribution to Risk of Lung Cancer

BACKGROUND

Genome-wide association studies are widely used to map genomic regions contributing to lung cancer (LC) susceptibility, but they typically do not identify the precise disease-causing genes/variants. To unveil the inherited genetic variants that cause LC, we performed focused exome-sequencing analyses on genes located in 121 genome-wide association study-identified loci previously implicated in the risk of LC, chronic obstructive pulmonary disease, pulmonary function level, and smoking behavior.

METHODS

Germline DNA from 260 case patients with LC and 318 controls were sequenced by utilizing VCRome 2.1 exome capture. Filtering was based on enrichment of rare and potential deleterious variants in cases (risk alleles) or controls (protective alleles). Allelic association analyses of single-variant and gene-based burden tests of multiple variants were performed. Promising candidates were tested in two independent validation studies with a total of 1773 case patients and 1123 controls.

RESULTS

We identified 48 rare variants with deleterious effects in the discovery analysis and validated 12 of the 43 candidates that were covered in the validation platforms. The top validated candidates included one well-established truncating variant, namely, BRCA2, DNA repair associated gene (BRCA2) K3326X (OR = 2.36, 95% confidence interval [CI]: 1.38-3.99), and three newly identified variations, namely, lymphotoxin beta gene (LTB) p.Leu87Phe (OR = 7.52, 95% CI: 1.01-16.56), prolyl 3-hydroxylase 2 gene (P3H2) p.Gln185His (OR = 5.39, 95% CI: 0.75-15.43), and dishevelled associated activator of morphogenesis 2 gene (DAAM2) p.Asp762Gly (OR = 0.25, 95% CI: 0.10-0.79). Burden tests revealed strong associations between zinc finger protein 93 gene (ZNF93), DAAM2, bromodomain containing 9 gene (BRD9), and the gene LTB and LC susceptibility.

CONCLUSION

Our results extend the catalogue of regions associated with LC and highlight the importance of germline rare coding variants in LC susceptibility.

Detection and quantitation of copy number variation in the voltage-gated sodium channel gene of the mosquito Culex quinquefasciatus

Insecticide resistance is typically associated with alterations to the insecticidal target-site or with gene expression variation at loci involved in insecticide detoxification. In some species copy number variation (CNV) of target site loci (e.g. the Ace-1 target site of carbamate insecticides) or detoxification genes has been implicated in the resistance phenotype. We show that field-collected Ugandan Culex quinquefasciatus display CNV for the voltage-gated sodium channel gene (Vgsc), target-site of pyrethroid and organochlorine insecticides. In order to develop field-applicable diagnostics for Vgsc CN, and as a prelude to investigating the possible association of CN with insecticide resistance, three assays were compared for their accuracy in CN estimation in this species. The gold standard method is droplet digital PCR (ddPCR), however, the hardware is prohibitively expensive for widespread utility. Here, ddPCR was compared to quantitative PCR (qPCR) and pyrosequencing. Across all platforms, CNV was detected in ≈10% of mosquitoes, corresponding to three or four copies (per diploid genome). ddPCR and qPCR-Std-curve yielded similar predictions for Vgsc CN, indicating that the qPCR protocol developed here can be applied as a diagnostic assay, facilitating monitoring of Vgsc CN in wild populations and the elucidation of association between the Vgsc CN and insecticide resistance.

AURKA mRNA expression is an independent predictor of poor prognosis in patients with non-small cell lung cancer

Deregulation of mitotic spindle genes has been reported to contribute to the development and progression of malignant tumours. The aim of the present study was to explore the association between the expression profiles of Aurora kinases (AURKA, AURKB and AURKC), cytoskeleton-associated protein 5 (CKAP5), discs large-associated protein 5 (DLGAP5), kinesin-like protein 11 (KIF11), microtubule nucleation factor (TPX2), monopolar spindle 1 kinase (TTK), and β-tubulins (TUBB) and (TUBB3) genes and clinicopathological characteristics in human non-small cell lung carcinoma (NSCLC). Reverse transcription-quantitative polymerase chain reaction-based RNA gene expression profiles of 132 NSCLC and 44 adjacent wild-type tissues were generated, and Cox's proportional hazard regression was used to examine associations. With the exception of AURKC, all genes exhibited increased expression in NSCLC tissues. Of the 10 genes examined, only AURKA was significantly associated with prognosis in NSCLC. Multivariate Cox's regression analysis demonstrated that AURKA mRNA expression [hazard ratio (HR), 1.81; 95% confidence interval (CI), 1.16-2.84; P=0.009], age (HR, 1.03; 95% CI, 1.00-1.06; P=0.020), pathological tumour stage 2 (HR, 2.43; 95% CI, 1.16-5.10; P=0.019) and involvement of distal nodes (pathological node stage 2) (HR, 3.14; 95% CI, 1.24-7.99; P=0.016) were independent predictors of poor prognosis in patients with NSCLC. Poor prognosis of patients with increased AURKA expression suggests that those patients may benefit from surrogate therapy with AURKA inhibitors.

Aurora B expression modulates paclitaxel response in non-small cell lung cancer

Background:

Taxanes are mitotic poisons widely used in the treatment of non-small cell lung cancer (NSCLC), however, little is known about potential molecular modulators of response to these compounds. Aurora B (AURKB) is a critical regulator of the mitotic spindle assembly, previously shown overexpressed in NSCLC. Here we investigated the hypothesis that AURKB expression modulates the efficacy of taxanes in NSCLC cells.

Methods:

AURKB mRNA expression was determined by qPCR in 132 frozen NSCLC tissues and nine NSCLC cell lines. Aurora B expression was knocked down in cell lines using multiple shRNA constructs. Barasertib was used to specifically inhibit AURKB activity, determined by the level of H3S10 phosphorylation.

Results:

Frequent AURKB mRNA upregulation was observed in NSCLC tissues (P<0.0001), being more prominent in squamous carcinomas (P<0.0001). Aurora B expression in cell lines strongly correlated with sensitivity to both docetaxel (P=0.004) and paclitaxel (P=0.007). Aurora B knockdown derivatives consistently showed a dose-dependent association between low-AURKB expression and resistance to paclitaxel. Specific chemical inhibition of Aurora B activity also demonstrated a strong dose-dependent efficiency in triggering paclitaxel resistance.

Conclusions:

Aurora B activity is an important modulator of taxane response in NSCLC cells. This may lead to further insights into taxane sensitivity of NSCLC tumours.

Intraoperative Sentinel Lymph Node Evaluation: Implications of Cytokeratin 19 Expression for the Adoption of OSNA in Oral Squamous Cell Carcinoma

BACKGROUND

Intraoperative analysis of sentinel lymph nodes would enhance the care of early-stage oral squamous cell carcinoma (OSCC). We determined the frequency and extent of cytokeratin 19 (CK19) expression in OSCC primary tumours and surrounding tissues to explore the feasibility of a "clinic-ready" intraoperative diagnostic test (one step nucleic acid amplification-OSNA, sysmex).

METHODS

Two cohorts were assembled: cohort 1, OSCC with stage and site that closely match cases suitable for sentinel lymph node biopsy (SLNB); cohort 2, HNSCC with sufficient fresh tumour tissue available for the OSNA assay (>50 mg). CK19 assays included qRT-PCR, RNA in situ hybridisation (ISH), and immunohistochemistry (IHC), as well as OSNA.

RESULTS

CK19 mRNA expression was detected with variable sensitivity, depending on method, in 60-80% of primary OSCC tumours, while protein expression was observed in only 50% of tumours. Discordance between different techniques indicated that OSNA was more sensitive than qRT-PCR or RNA-ISH, which in turn were more sensitive than IHC. OSNA results showed CK19 expression in 80% of primary cases, so if used for diagnosis of lymph node metastasis would lead to a false-negative result in 20% of patients with cervical lymph node metastases.

CONCLUSIONS

OSNA in its current form is not suitable for use in OSCC SLNB due to inadequate expression of the CK19 target in all case. However, the same assay technology would likely be very promising if applied using a more ubiquitous squamous epithelial target.

Length of Variable Numbers of Tandem Repeats in the Carboxyl Ester Lipase (CEL) Gene May Confer Susceptibility to Alcoholic Liver Cirrhosis but Not Alcoholic Chronic Pancreatitis

BACKGROUND

Carboxyl-ester lipase (CEL) contributes to fatty acid ethyl ester metabolism, which is implicated in alcoholic pancreatitis. The CEL gene harbours a variable number of tandem repeats (VNTR) region in exon 11. Variation in this VNTR has been linked to monogenic pancreatic disease, while conflicting results were reported for chronic pancreatitis (CP). Here, we aimed to investigate a potential association of CEL VNTR lengths with alcoholic CP.

METHODS

Overall, 395 alcoholic CP patients, 218 patients with alcoholic liver cirrhosis (ALC) serving as controls with a comparable amount of alcohol consumed, and 327 healthy controls from Germany and the United Kingdom (UK) were analysed by determination of fragment lengths by capillary electrophoresis. Allele frequencies and genotypes of different VNTR categories were compared between the groups.

RESULTS

Twelve repeats were overrepresented in UK ACP patients (P = 0.04) compared to controls, whereas twelve repeats were enriched in German ALC compared to alcoholic CP patients (P = 0.03). Frequencies of CEL VNTR lengths of 14 and 15 repeats differed between German ALC patients and healthy controls (P = 0.03 and 0.008, respectively). However, in the genotype and pooled analysis of VNTR lengths no statistical significant association was depicted. Additionally, the 16-16 genotype as well as 16 repeats were more frequent in UK ALC than in alcoholic CP patients (P = 0.034 and 0.02, respectively). In all other calculations, including pooled German and UK data, allele frequencies and genotype distributions did not differ significantly between patients and controls or between alcoholic CP and ALC.

CONCLUSIONS

We did not obtain evidence that CEL VNTR lengths are associated with alcoholic CP. However, our results suggest that CEL VNTR lengths might associate with ALC, a finding that needs to be clarified in larger cohorts.

UHRF1 regulation of the Keap1-Nrf2 pathway in pancreatic cancer contributes to oncogenesis

The cellular defence protein Nrf2 is a mediator of oncogenesis in pancreatic ductal adenocarcinoma (PDAC) and other cancers. However, the control of Nrf2 expression and activity in cancer is not fully understood. We previously reported the absence of Keap1, a pivotal regulator of Nrf2, in ∼70% of PDAC cases. Here we describe a novel mechanism whereby the epigenetic regulator UHRF1 suppresses Keap1 protein levels. UHRF1 expression was observed in 20% (5 of 25) of benign pancreatic ducts compared to 86% (114 of 132) of pancreatic tumours, and an inverse relationship between UHRF1 and Keap1 levels in PDAC tumours (n = 124) was apparent (p = 0.002). We also provide evidence that UHRF1-mediated regulation of the Nrf2 pathway contributes to the aggressive behaviour of PDAC. Depletion of UHRF1 from PDAC cells decreased growth and enhanced apoptosis and cell cycle arrest. UHRF1 depletion also led to reduced levels of Nrf2-regulated downstream proteins and was accompanied by heightened oxidative stress, in the form of lower glutathione levels and increased reactive oxygen species. Concomitant depletion of Keap1 and UHRF1 restored Nrf2 levels and reversed cell cycle arrest and the increase in reactive oxygen species. Mechanistically, depletion of UHRF1 reduced global and tumour suppressor promoter methylation in pancreatic cancer cell lines, and KEAP1 gene promoter methylation was reduced in one of three cell lines examined. Thus, methylation of the KEAP1 gene promoter may contribute to the suppression of Keap1 protein levels by UHRF1, although our data suggest that additional mechanisms need to be explored. Finally, we demonstrate that K-Ras drives UHRF1 expression, establishing a novel link between this oncogene and Nrf2-mediated cellular protection. Since UHRF1 over-expression occurs in other cancers, its ability to regulate the Keap1-Nrf2 pathway may be critically important to the malignant behaviour of these cancers.

Identification of lung cancer histology-specific variants applying Bayesian framework variant prioritization approaches within the TRICL and ILCCO consortia

Large-scale genome-wide association studies (GWAS) have likely uncovered all common variants at the GWAS significance level. Additional variants within the suggestive range (0.0001> P > 5×10(-8)) are, however, still of interest for identifying causal associations. This analysis aimed to apply novel variant prioritization approaches to identify additional lung cancer variants that may not reach the GWAS level. Effects were combined across studies with a total of 33456 controls and 6756 adenocarcinoma (AC; 13 studies), 5061 squamous cell carcinoma (SCC; 12 studies) and 2216 small cell lung cancer cases (9 studies). Based on prior information such as variant physical properties and functional significance, we applied stratified false discovery rates, hierarchical modeling and Bayesian false discovery probabilities for variant prioritization. We conducted a fine mapping analysis as validation of our methods by examining top-ranking novel variants in six independent populations with a total of 3128 cases and 2966 controls. Three novel loci in the suggestive range were identified based on our Bayesian framework analyses: KCNIP4 at 4p15.2 (rs6448050, P = 4.6×10(-7)) and MTMR2 at 11q21 (rs10501831, P = 3.1×10(-6)) with SCC, as well as GAREM at 18q12.1 (rs11662168, P = 3.4×10(-7)) with AC. Use of our prioritization methods validated two of the top three loci associated with SCC (P = 1.05×10(-4) for KCNIP4, represented by rs9799795) and AC (P = 2.16×10(-4) for GAREM, represented by rs3786309) in the independent fine mapping populations. This study highlights the utility of using prior functional data for sequence variants in prioritization analyses to search for robust signals in the suggestive range.

Strategies for validation and testing of DNA methylation biomarkers

DNA methylation is a stable covalent epigenetic modification of primarily CpG dinucleotides that has recently gained considerable attention for its use as a biomarker in different clinical settings, including disease diagnosis, prognosis and therapeutic response prediction. Although the advent of genome-wide DNA methylation profiling in primary disease tissue has provided a manifold resource for biomarker development, only a tiny fraction of DNA methylation-based assays have reached clinical testing. Here, we provide a critical overview of different analytical methods that are suitable for biomarker validation, including general study design considerations, which might help to streamline epigenetic marker development. Furthermore, we highlight some of the recent marker validation studies and established markers that are currently commercially available for assisting in clinical management of different cancers.

Abstract B20: Global hypomethylation of esophageal myofibroblasts is associated with accelerated tumor growth

Cancer associated myofibroblasts (CAMs) originating from either tissue-resident fibroblast precursors or from bone marrow derived mesenchymal stromal cells, contribute to the tumor microenvironment by the secretion of growth factors, cytokines, proteases and ECM proteins. Exaggerated CAM reactions have been associated with adverse outcomes in various cancers.

The role of altered DNA CpG-island methylation in promoting tumorigenesis has been widely studied in cancer cells; altered DNA methylation in cancer stromal cells including CAMs has also been described but its biological significance remains uncertain. Our present objective was to determine the significance of DNA methylation in influencing the properties of esophageal CAMs. In initial studies we demonstrated different DNA methylation patterns in seven low passage primary myofibroblast lines from esophageal adenocarcinomas (n=3) and squamous cell carcinoma (n=4) patients in comparison to matched myofibroblasts derived from adjacent, macroscopically tumor-free, tissue (ATMs). Further, we found changes in global DNA methylation assessed in LINE1 elements by pyrosequencing analysis between CAMs and matched ATMs. Cluster analysis of the myofibroblast methylome determined by Illumina 450k human DNA methylation array revealed distinct methylation patterns in CAMs and ATMs. In order to generate a model for studying the biological consequences of epigenetic modification of myofibroblasts we treated cultured primary low passage normal esophageal tissue myofibroblasts (NTMs) with the DNA methyl transferase inhibitor, 5'Aza-2'deoxycytidin (DAC). A three-day exposure of NTMs to 0.1μM DAC induced a transient 10% decrease of global DNA methylation assessed by LINE1 element analysis; DAC (0.1 μM) also decreased myofibroblast proliferation by 20 %, and this persisted up to 28 days after treatment, and it stably increased migration in transwell migration assays by 50%. In a nude mouse xenograft model using esophageal squamous carcinoma cells (OE21) inclusion of CAMs stimulated tumor growth more than the inclusion of ATMs; moreover, myofibroblasts pretreated with DAC increased tumor growth compared with xenografts of OE21cells and untreated myofibroblasts.

Our data indicate (a) modification of DNA methylation in esophageal cancer myofibroblasts, and (b) that DNA demethylation induced by DAC stimulates migration of CAMs as well as their capacity to promote xenograft growth in vivo. Understanding the mechanisms and consequences of epigenetic changes in esophageal stromal cells may be useful in developing future anticancer strategies.

Citation Format: Olivier T. Giger, Benjamin Brown, Triantafillos Liloglou, Puthen V. Jithesh, J Graham Dockray, Andrea Varro. Global hypomethylation of esophageal myofibroblasts is associated with accelerated tumor growth.. [abstract]. In: Abstracts: AACR Special Conference on Cellular Heterogeneity in the Tumor Microenvironment; 2014 Feb 26-Mar 1; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(1 Suppl):Abstract nr B20. doi:10.1158/1538-7445.CHTME14-B20

Squamous cell carcinoma of the head and neck outside the oropharynx is rarely human papillomavirus related

OBJECTIVES/HYPOTHESIS

The incidence of human papillomavirus (HPV)-driven disease beyond the oropharynx varies greatly in the reported literature.

STUDY DESIGN

Case series.

METHODS

Two hundred twenty-one samples were strictly classified to the subsites of oral cavity, larynx, or hypopharynx at the time of primary surgery. Formalin-fixed paraffin-embedded samples were subjected to a validated, tiered, diagnostic algorithm of p16 immunohistochemistry, high-risk HPV in situ hybridization, and quantitative polymerase chain reaction for HPV E6 DNA. An additional 60 oropharyngeal cases acted as an internal biological control.

RESULTS

An incidence of 4% of HPV-driven cases was observed across the subsites outside the oropharynx compared to 70% of tumors confined within it.

CONCLUSIONS

This is the first reporting of a broad range of nonoropharyngeal HPV rates using this validated diagnostic algorithm. It remains unclear whether patients with HPV-driven disease originating outside the oropharynx enjoy the same survival advantage apparent in those patients with oropharyngeal squamous cell carcinomas.

LEVEL OF EVIDENCE

4 Laryngoscope, 124:2739-2744, 2014.

SERPINE1 and SMA expression at the invasive front predict extracapsular spread and survival in oral squamous cell carcinoma

BACKGROUND

Extracapsular spread (ECS) in cervical lymph nodes is the single-most prognostic clinical variable in oral squamous cell carcinoma (OSCC), but diagnosis is possible only after histopathological examination. A promising biomarker in the primary tumour, alpha smooth muscle actin (SMA) has been shown to be highly prognostic, however, validated biomarkers to predict ECS prior to primary treatment are not yet available.

METHODS

In 102 OSCC cases, conventional imaging was compared with pTNM staging. SERPINE1, identified from expression microarray of primary tumours as a potential biomarker for ECS, was validated through mRNA expression, and by immunohistochemistry (IHC) on a tissue microarray from the same cohort. Similarly, expression of SMA was also compared with its association with ECS and survival. Expression was analysed separately in the tumour centre and advancing front; and prognostic capability determined using Kaplan-Meier survival analysis.

RESULTS

Immunohistochemistry indicated that both SERPINE1 and SMA expression at the tumour-advancing front were significantly associated with ECS (P<0.001). ECS was associated with expression of either or both proteins in all cases. SMA+/SERPINE1+ expression in combination was highly significantly associated with poor survival (P<0.001). MRI showed poor sensitivity for detection of nodal metastasis (56%) and ECS (7%). Both separately, and in combination, SERPINE1 and SMA were superior to MRI for the detection of ECS (sensitivity: SERPINE1: 95%; SMA: 82%; combination: 81%).

CONCLUSION

A combination of SMA and SERPINE1 IHC offer potential as prognostic biomarkers in OSCC. Our findings suggest that biomarkers at the invasive front are likely to be necessary in prediction of ECS or in therapeutic stratification.

Serum mannose-binding lectin concentration, but not genotype, is associated with Clostridium difficile infection recurrence: a prospective cohort study

Low mannose-binding lectin concentration, but not genotype, was associated with disease recurrence in a large prospective cohort of patients with Clostridium difficile infection.

Background. Mannose-binding lectin (MBL) plays a key role in the activation of the lectin-complement pathway of innate immunity, and its deficiency has been linked with several acute infections. However, its role in predisposing to, or modulating disease severity in, Clostridium difficile infection (CDI) has not been investigated.

Methods. We prospectively recruited 308 CDI case patients and 145 control patients with antibiotic-associated diarrhea (AAD). CDI outcome measures were disease severity, duration of symptoms, 30-day mortality, and 90-day recurrence. Serum concentrations of MBL were determined using a commercial enzyme-linked immunosorbent assay transferred to an electrochemiluminescence–based platform. MBL2 polymorphisms were typed using a combination of pyrosequencing and TaqMan genotyping assays.

Results. The frequency of the MBL2 genetic variants was similar to that reported in other white populations. MBL serum concentrations in CDI and AAD subjects were determined by MBL2 exonic variants B, C, and D and the haplotypes (LYPB, LYQC, and HYPD). There was no difference in either MBL concentrations or genotypes between cases and controls. MBL concentration, but not genotype, was a determinant of CDI recurrence (odds ratios, 3.18 [95% confidence interval {CI}, 1.40–7.24] and 2.61 [95% CI, 1.35–5.04] at the <50 ng/mL and <100 ng/mL cutoff points, respectively; P < .001). However, neither MBL concentration nor MBL2 genotype was linked with the other CDI outcomes.

Conclusions. Serum MBL concentration did not differentiate between CDI cases and AAD controls, but among CDI cases, MBL concentration, but not genotype, was associated with CDI recurrence, indicating that MBL acts as a modulator of disease, rather than a predisposing factor.

Exposure to secondhand tobacco smoke and lung cancer by histological type: a pooled analysis of the International Lung Cancer Consortium (ILCCO)

While the association between exposure to secondhand smoke and lung cancer risk is well established, few studies with sufficient power have examined the association by histological type. In this study, we evaluated the secondhand smoke-lung cancer relationship by histological type based on pooled data from 18 case-control studies in the International Lung Cancer Consortium (ILCCO), including 2,504 cases and 7,276 control who were never smokers and 10,184 cases and 7,176 controls who were ever smokers. We used multivariable logistic regression, adjusting for age, sex, race/ethnicity, smoking status, pack-years of smoking, and study. Among never smokers, the odds ratios (OR) comparing those ever exposed to secondhand smoke with those never exposed were 1.31 (95% CI: 1.17-1.45) for all histological types combined, 1.26 (95% CI: 1.10-1.44) for adenocarcinoma, 1.41 (95% CI: 0.99-1.99) for squamous cell carcinoma, 1.48 (95% CI: 0.89-2.45) for large cell lung cancer, and 3.09 (95% CI: 1.62-5.89) for small cell lung cancer. The estimated association with secondhand smoke exposure was greater for small cell lung cancer than for nonsmall cell lung cancers (OR=2.11, 95% CI: 1.11-4.04). This analysis is the largest to date investigating the relation between exposure to secondhand smoke and lung cancer. Our study provides more precise estimates of the impact of secondhand smoke on the major histological types of lung cancer, indicates the association with secondhand smoke is stronger for small cell lung cancer than for the other histological types, and suggests the importance of intervention against exposure to secondhand smoke in lung cancer prevention.

Quantitative measurement of Human Papillomavirus type 16 L1/L2 DNA methylation correlates with cervical disease grade

BACKGROUND

Persistent infection with Human Papillomavirus (HPV) type 16 causes the majority of cervical cancers. Genital HPV infection is very common, but neoplastic progression is uncommon. There is an urgent need to identify biomarkers associated with cervical neoplasia that can be used to triage women who test positive for HPV.

OBJECTIVES

To assess the ability of quantitative measurement of HPV16 DNA methylation to separate samples of different cytology grades and cervical cancers, and determine which of the assessed regions of the HPV genome and individual CpGs are most informative.

STUDY DESIGN

DNA methylation was quantified by pyrosequencing of bisulphite converted DNA from liquid based cytology samples from 17 women with normal cytology and 20 women with severe dyskaryosis, and from fixed tissue from 24 women with cervical cancer. Methylation was assessed in the HPV Long Control Region (LCR), E2 and L1/L2 regions.

RESULTS

In cervical cancers, increased HPV DNA methylation was present in all regions. Increased methylation was also observed in severely dyskaryotic relative to normal samples, but only in the E2 and L1/L2 regions. The ability of methylation based classifiers to separate the three classes of material was assessed by ROC curve analyses. The best separation between normal and dyskaryotic samples was achieved by assessment of the L1/L2 CpGs at nucleotide positions 5600 and 5609 (AUC=0.900, 95% CI: 0.793-1).

CONCLUSIONS

This study demonstrates the potential of quantification of HPV DNA methylation as a biomarker of cervical neoplasia. An algorithm considering methylation at specific L1/L2 CpGs appeared the most promising model.

A microRNA-based prediction algorithm for diagnosis of non-small lung cell carcinoma in minimal biopsy material

BACKGROUND

Diagnosis is jeopardised when limited biopsy material is available or histological quality compromised. Here we developed and validated a prediction algorithm based on microRNA (miRNA) expression that can assist clinical diagnosis of lung cancer in minimal biopsy material to improve clinical management.

METHODS

Discovery utilised Taqman Low Density Arrays (754 miRNAs) in 20 non-small cell lung cancer (NSCLC) tumour/normal pairs. In an independent set of 40 NSCLC patients, 28 miRNA targets were validated using qRT-PCR. A prediction algorithm based on eight miRNA targets was validated blindly in a third independent set of 47 NSCLC patients. The panel was also tested in formalin-fixed paraffin-embedded (FFPE) specimens from 20 NSCLC patients. The genomic methylation status of highly deregulated miRNAs was investigated by pyrosequencing.

RESULTS

In the final, frozen validation set the panel had very high sensitivity (97.5%), specificity (96.3%) and ROC-AUC (0.99, P=10(-15)). The panel provided 100% sensitivity and 95% specificity in FFPE tissue (ROC-AUC=0.97 (P=10(-6))). DNA methylation abnormalities contribute little to the deregulation of the miRNAs tested.

CONCLUSION

The developed prediction algorithm is a valuable potential biomarker for assisting lung cancer diagnosis in minimal biopsy material. A prospective validation is required to measure the enhancement of diagnostic accuracy of our current clinical practice.

A novel in silico reverse-transcriptomics-based identification and blood-based validation of a panel of sub-type specific biomarkers in lung cancer

Lung cancer accounts for the highest number of cancer-related deaths worldwide. Early diagnosis significantly increases the disease-free survival rate and a large amount of effort has been expended in screening trials and the development of early molecular diagnostics. However, a gold standard diagnostic strategy is not yet available. Here, based on miRNA expression profile in lung cancer and using a novel in silico reverse-transcriptomics approach, followed by analysis of the interactome; we have identified potential transcription factor (TF) markers that would facilitate diagnosis of subtype specific lung cancer. A subset of seven TF markers has been used in a microarray screen and was then validated by blood-based qPCR using stage-II and IV non-small cell lung carcinomas (NSCLC). Our results suggest that overexpression of HMGA1, E2F6, IRF1, and TFDP1 and downregulation or no expression of SUV39H1, RBL1, and HNRPD in blood is suitable for diagnosis of lung adenocarcinoma and squamous cell carcinoma sub-types of NSCLC. Here, E2F6 was, for the first time, found to be upregulated in NSCLC blood samples. The miRNA-TF-miRNA interaction based molecular mechanisms of these seven markers in NSCLC revealed that HMGA1 and TFDP1 play vital roles in lung cancer tumorigenesis. The strategy developed in this work is applicable to any other cancer or disease and can assist in the identification of potential biomarkers.

A prognostic DNA methylation signature for stage I non-small-cell lung cancer

PURPOSE

Non-small-cell lung cancer (NSCLC) is a tumor in which only small improvements in clinical outcome have been achieved. The issue is critical for stage I patients for whom there are no available biomarkers that indicate which high-risk patients should receive adjuvant chemotherapy. We aimed to find DNA methylation markers that could be helpful in this regard.

PATIENTS AND METHODS

A DNA methylation microarray that analyzes 450,000 CpG sites was used to study tumoral DNA obtained from 444 patients with NSCLC that included 237 stage I tumors. The prognostic DNA methylation markers were validated by a single-methylation pyrosequencing assay in an independent cohort of 143 patients with stage I NSCLC.

RESULTS

Unsupervised clustering of the 10,000 most variable DNA methylation sites in the discovery cohort identified patients with high-risk stage I NSCLC who had shorter relapse-free survival (RFS; hazard ratio [HR], 2.35; 95% CI, 1.29 to 4.28; P = .004). The study in the validation cohort of the significant methylated sites from the discovery cohort found that hypermethylation of five genes was significantly associated with shorter RFS in stage I NSCLC: HIST1H4F, PCDHGB6, NPBWR1, ALX1, and HOXA9. A signature based on the number of hypermethylated events distinguished patients with high- and low-risk stage I NSCLC (HR, 3.24; 95% CI, 1.61 to 6.54; P = .001).

CONCLUSION

The DNA methylation signature of NSCLC affects the outcome of stage I patients, and it can be practically determined by user-friendly polymerase chain reaction assays. The analysis of the best DNA methylation biomarkers improved prognostic accuracy beyond standard staging.

Lung cancer epigenetics: emerging biomarkers

Lung cancer is the leading cause of cancer-related deaths worldwide, and the 5-year survival rate is still very poor due to the scarcity of effective tools for early detection. The discovery of highly sensitive and specific biomarkers highlighting pathological changes early enough to allow clinical intervention is therefore of great importance. In the last decade, epigenetics and particularly research on DNA methylation have provided important information towards a better understanding of lung cancer pathogenesis. Novel and promising molecular biomarkers for diagnosis and prognosis of lung cancer are continuously emerging in this area, requiring further evaluation. This process includes extensive validation in prospective clinical trials before they can be routinely used in a clinical setting. This review summarizes the evidence on epigenetic biomarkers for lung cancer, focusing on DNA methylation.

Cytoglobin has bimodal: tumour suppressor and oncogene functions in lung cancer cell lines

Cytoglobin (CYGB) is frequently downregulated in many types of human malignancies, and its exogenous overexpression reduces proliferation of cancer cells. Despite its implied tumour suppressor (TSG) functions, its exact role in carcinogenesis remains unclear as CYGB upregulation is also associated with tumour hypoxia and aggressiveness. In this study, we explore the TSG role of CYGB, its influence on the phenotype of cancerous cells under stress conditions and the clinical significance of CYGB expression and promoter methylation in non-small cell lung cancer (NSCLC). DNA methylation-dependent expression silencing of CYGB is demonstrated in both clinical samples and cell lines. CYGB promoter was more frequently methylated in lung adenocarcinomas (P = 1.4 × 10(-4)). Demethylation by 5'-azadeoxycytidine partially restored CYGB expression in cell lines. Interestingly, trichostatin A triggered upregulation of CYGB expression in cancer cell lines and downregulation in non-tumourigenic ones. CYGB mRNA expression in NSCLC surgical specimens correlated with that of HIF1α and VEGFa (P < 1 × 10(-4)). Overexpression of CYGB in cancer cell lines reduced cell migration, invasion and anchorage-independent growth. Moreover, CYGB impaired cell proliferation, but only in the lung adenocarcinoma cell line (H358). Upon hydrogen peroxide treatment, CYGB protected cell viability, migratory potential and anchorage independence by attenuating oxidative injury. In hypoxia, CYGB overexpression decreased cell viability, augmented migration and anchorage independence in a cell-type-specific manner. In conclusion, CYGB revealed TSG properties in normoxia but promoted tumourigenic potential of the cells exposed to stress, suggesting a bimodal function in lung tumourigenesis, depending on cell type and microenvironmental conditions.

DNA methylation biomarkers in biological fluids for early detection of respiratory tract cancer

Cancers of the respiratory tract (lung and head and neck) share common aetiologies, risk factors and molecular characteristics. Epigenetic reprogramming is one of the hallmarks of cancer and DNA methylation is currently the best-studied form. There are a number of characteristics of DNA methylation, which seem advantageous in biomarker development. Early detection is still an unmet clinical care need, which guarantees to significantly reduce the mortality of patients with respiratory cancers. The application of such biomarkers in biological fluids being sampled in everyday clinical practice is a long-term demand. In this review we summarise the current literature on DNA methylation detection in bronchial washings, sputum, saliva, plasma and serum and discuss the potential of their clinical implementation. We also discuss important aspects of biomarker development and validation pointing to the appropriate route for a biomarker to reach clinical standards.