Search

Summary: 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.

Abstract 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.

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.

Abstract 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.

Abstract 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.

Abstract 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.

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

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.

Similarities in cancers can be studied to interrogate their etiology. Here, the authors use genome-wide association study summary statistics from six cancer types based on 296,215 cases and 301,319 controls of European ancestry, showing that solid tumours arising from different tissues share a degree of common germline genetic basis.

Abstract 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−9), general transcription factor IIH subunit H4 (GTF2H4) (p = 7.03 × 10−7), and Tenascin X (TNXB) (p = 6.30 × 10−6) 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).

Summary: 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.

Abstract 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.

Background: The high relapse and mortality rate of small-cell lung cancer (SCLC) fuels the need for epidemiologic study to aid in its prevention. Methods: We included 24 studies from the ILCCO collaboration. Random-effects panel logistic regression and cubic spline regression were used to estimate the effects of smoking behaviors on SCLC risk and explore their non-linearity. Further, we explored whether the risk of smoking on SCLC was mediated through COPD. Findings: Significant dose–response relationships of SCLC risk were observed for all quantitative smoking variables. Smoking pack-years were associated with a sharper increase of SCLC risk for pack-years ranged 0 to approximately 50. The former smokers with longer cessation showed a 43%quit_for_5–9 years to 89%quit_for_≥20 years declined SCLC risk vs. subjects who had quit smoking

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.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.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.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.

Detection of low copies of methylated DNA targets in clinical specimens is challenging. The quantitative Methylation-Specific PCR (qMSP) assays were designed to specifically amplify bisulphite-converted methylated DNA target sequences in the presence of an excess of unmethylated counterpart sequences. These qMSP assays are real-time PCR assays utilizing, sequence-specific primers and an intervening, also sequence specific, Taqman probe to cover an amplicon of approximately 100 bp in length. The use of Taqman probes bearing a minor groove binding (MGB) allow for the use of shorter probes and therefore facilitate design and significantly increases the analytical specificity of the reaction. In the context of the biomarker discovery program of the Liverpool Lung Project (LLP), ten gene promoters were selected. qMSP assays were developed, validated and used to screen 655 bronchial washings from patients with lung cancer and age/sex matched controls with non malignant lung disease (Nikolaidis et al., 2012).

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p ≤ 5 × 10⁻⁷). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1×10⁻⁸) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p =2 × 10⁻⁸) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5 × 10⁻⁸); rs1229984-ADH1B, p = 7 × 10⁻⁹; and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.