Search

Little is known about the association of diet with risk of bladder cancer. This might be due to the fact that the majority of studies have focused on single food items, rather than dietary patterns, which may better capture any influence of diet on bladder cancer risk. We aimed to investigate the association between a measure of Western dietary pattern and bladder cancer risk. Associations between adherence to a Western dietary pattern and risk of developing bladder cancer were assessed by pooling data from 13 prospective cohort studies in the “BLadder cancer Epidemiology and Nutritional Determinants” (BLEND) study and applying Cox regression analysis. Dietary data from 580 768 study participants, including 3401 incident cases, and 577 367 noncases were analyzed. A direct and significant association was observed between higher adherence to a Western dietary pattern and risk of bladder cancer (hazard ratio (HR) comparing highest with lowest tertile scores: 1.54, 95% confidence interval (CI): 1.37, 1.72; P‐trend = .001). This association was observed for men (HR comparing highest with lowest tertile scores: 1.72; 95% CI: 1.51, 1.96; P‐trend = .001), but not women (P‐het = .001). Results were consistent with HR above 1.00 after stratification on cancer subtypes (nonmuscle‐invasive and muscle‐invasive bladder cancer). We found evidence that adherence to a Western dietary pattern is associated with an increased risk of bladder cancer for men but not women., What's new? Does diet affect bladder‐cancer risk? Individual foods are rarely eaten in isolation, but little is known about the impact of overall dietary habits. In this large, prospective study, the authors found that greater adherence to a Western dietary pattern was associated with a significantly increased risk of bladder cancer in men. (Surprisingly, the same effect was not seen in women.) Further research is needed to identify the specific food types responsible and their mechanisms of bladder carcinogenesis. However, education to encourage changes in general dietary habits may provide a valuable public‐health benefit.

BACKGROUND: Recurrent mutations in the promoter of the telomerase reverse transcriptase (TERT) gene (C228T and C250T) detected in tumours and cells shed into urine of urothelial cancer (UC) patients are putative biomarkers for UC detection and monitoring. However, the possibility of detecting these mutations in cell-free circulating DNA (cfDNA) in blood and urine, or DNA from urinary exfoliated cells (cellDNA) with a single-gene sensitive assay has never been tested in a case-control setting. METHODS: We developed a single-plex assay (UroMuTERT) for the detection of low-abundance TERT promoter mutations. We tested 93 primary and recurrent UC cases and 94 controls recruited in France (blood, urine samples and tumours for the cases), and 50 primary UC cases and 50 controls recruited in Portugal (urinary exfoliated cell samples). We compared our assay with urine cytology. FINDINGS: In the French series, C228T or C250T were detected in urinary cfDNA or cellDNA in 81 cases (87·1%; 95% CI 78·6–93·2), and five controls (Specificity 94·7%; 95%CI 88·0–98·3), with 98·6% (95% CI 92·5–99·96) concordance in matched tumours. Detection rate in plasma cfDNA among cases was 7·1%. The UroMuTERT sensitivity was (i) highest for urinary cfDNA and cellDNA combined, (ii) consistent across primary and recurrent cases, tumour stages and grades, (iii) higher for low-risk non-muscle invasive UC (86·1%) than urine cytology (23·0%) (P

Background Although a potential inverse association between vegetable intake and bladder cancer risk has been reported, epidemiological evidence is inconsistent. This research aimed to elucidate the association between vegetable intake and bladder cancer risk by conducting a pooled analysis of data from prospective cohort studies. Methods Vegetable intake in relation to bladder cancer risk was examined by pooling individual-level data from 13 cohort studies, comprising 3203 cases among a total of 555,685 participants. Pooled multivariate hazard ratios (HRs), with corresponding 95% confidence intervals (CIs), were estimated using Cox proportional hazards regression models stratified by cohort for intakes of total vegetable, vegetable subtypes (i.e. non-starchy, starchy, green leafy and cruciferous vegetables) and individual vegetable types. In addition, a diet diversity score was used to assess the association of the varied types of vegetable intake on bladder cancer risk. Results The association between vegetable intake and bladder cancer risk differed by sex (P-interaction = 0.011) and smoking status (P-interaction = 0.038); therefore, analyses were stratified by sex and smoking status. With adjustment of age, sex, smoking, energy intake, ethnicity and other potential dietary factors, we found that higher intake of total and non-starchy vegetables were inversely associated with the risk of bladder cancer among women (comparing the highest with lowest intake tertile: HR = 0.79, 95% CI = 0.64–0.98, P = 0.037 for trend, HR per 1 SD increment = 0.89, 95% CI = 0.81–0.99; HR = 0.78, 95% CI = 0.63–0.97, P = 0.034 for trend, HR per 1 SD increment = 0.88, 95% CI = 0.79–0.98, respectively). However, no evidence of association was observed among men, and the intake of vegetable was not found to be associated with bladder cancer when stratified by smoking status. Moreover, we found no evidence of association for diet diversity with bladder cancer risk. Conclusion Higher intakes of total and non-starchy vegetable are associated with reduced risk of bladder cancer for women. Further studies are needed to clarify whether these results reflect causal processes and potential underlying mechanisms.

In the last few years, we published two special issues devoted to highlighting important scientific results in the field of bladder cancer research and clinical implications [...]

Epigenetic mechanisms such as aberrant DNA methylation (DNAme) are known to drive esophageal squamous cell carcinoma (ESCC), yet they remain poorly understood. Here, we studied tumor-specific DNAme in ESCC cases from nine high-incidence countries of Africa, Asia, and South America. Infinium MethylationEPIC array was performed on 108 tumors and 51 normal tissues adjacent to the tumors (NAT) in the discovery phase, and targeted pyrosequencing was performed on 132 tumors and 36 NAT in the replication phase. Top genes for replication were prioritized by weighting methylation results using RNA-sequencing data from The Cancer Genome Atlas and GTEx and validated by qPCR. Methylome analysis comparing tumor and NAT identified 6,796 differentially methylated positions (DMP) and 866 differential methylated regions (DMR), with a 30% methylation (Δβ) difference. The majority of identified DMPs and DMRs were hypermethylated in tumors, particularly in promoters and gene-body regions of genes involved in transcription activation. The top three prioritized genes for replication, PAX9, SIM2, and THSD4, had similar methylation differences in the discovery and replication sets. These genes were exclusively expressed in normal esophageal tissues in GTEx and downregulated in tumors. The specificity and sensitivity of these DNAme events in discriminating tumors from NAT were assessed. Our study identified novel, robust, and crucial tumor-specific DNAme events in ESCC tumors across several high-incidence populations of the world. Methylome changes identified in this study may serve as potential targets for biomarker discovery and warrant further functional characterization. Significance: This largest genome-wide DNA methylation study on ESCC from high-incidence populations of the world identifies functionally relevant and robust DNAme events that could serve as potential tumor-specific markers.

Background & Aims Although HBV is a major cause of death in Africa, its genetic variability has been poorly documented. This study aimed to address whether HBV genotype and surface gene variants are associated with HBV-related liver disease in The Gambia. Methods We conducted a case-control study nested in the Prevention of Liver Fibrosis and Cancer in Africa programme. Consecutive treatment-naive patients with chronic HBV infection and detectable viral load were recruited: 211 controls with no significant liver disease and 91 cases (56 cirrhosis and 35 HCC cases). HBV genotypes and surface gene variants were determined by Sanger sequencing or next-generation sequencing (NGS) in serum DNA. Aflatoxin B1 (AFB1)-specific codon 249 TP53 mutation was determined by NGS in circulating cell-free plasma DNA. Results In phylogenetic analysis, 85% of individuals carried HBV genotype E, 14% genotype A, and 1% A/E recombinant viruses. Surface gene variants were more frequently observed in cases (43% and 57% in cirrhosis and HCC cases, respectively) than controls (25%; p 2,000 IU/ml (OR 22.7 [8.0–64.9]), HBsAg levels, Graphical abstract, Highlights • The rarer HBV genotype A (vs. E) is associated with a higher risk of HCC (odds ratio 5.2). • A genotype-independent deletion hotspot was found in the preS2 domain (preS2Δ38–55). • PreS2Δ38–55 was significantly associated with HCC and is an independent risk factor.

Human papillomaviruses are currently classified in different genera. Mucosal HPVs belonging to the alpha genus have been clearly associated with carcinogenesis of the mucosal epithelium at different sites. Beta HPV types have been classified as cutaneous. Although findings indicate that some beta HPVs from species 1 and 2 play a role, together with UV irradiation, in skin cancer, very little is known about the transforming properties of most of the beta HPVs. This report shows the transforming activity of E6 and E7 from beta-3 HPV types. Moreover, it highlights that beta-3 HPVs share some biological properties more extensively with mucosal high-risk HPV16 than with beta-2 HPV38. This report provides new paradigms for a better understanding of the biology of the different HPV types and their possible association with lesions at mucosal and/or cutaneous epithelia., The beta human papillomaviruses (HPVs) are subdivided into 5 species (beta-1 to beta-5), and they were first identified in the skin. However, the beta-3 species appears to be more highly represented in the mucosal epithelia than in the skin. Functional studies have also highlighted that beta-3 HPV49 shares some functional similarities with mucosal high-risk (HR) HPV16. Here, we describe the characterization of the in vitro transforming properties of the entire beta-3 species, which includes three additional HPV types: HPV75, HPV76, and HPV115. HPV49, HPV75, and HPV76 E6 and E7 (E6/E7), but not HPV115 E6 and E7, efficiently inactivate the p53 and pRb pathways and immortalize or extend the life span of human foreskin keratinocytes (HFKs). As observed for HR HPV16, cell cycle deregulation mediated by beta-3 HPV E6/E7 expression leads to p16INK4a accumulation, whereas no p16INK4a was detected in beta-2 HPV38 E6/E7 HFKs. As shown for HPV49 E6, HPV75 and HPV76 E6s degrade p53 by an E6AP/proteasome-mediated mechanism. Comparative analysis of cellular gene expression patterns of HFKs containing E6 and E7 from HR HPV16, beta-3 HPV types, and beta-2 HPV38 further highlights the functional similarities of HR HPV16 and beta-3 HPV49, HPV75, and HPV76. The expression profiles of these four HPV HFKs show some similarities and diverge substantially from those of beta-3 HPV115 E6/E7 and beta-2 HPV38 E6/E7 HFKs. In summary, our data show that beta-3 HPV types share some mechanisms with HR HPV types and pave the way for additional studies aiming to evaluate their potential role in human pathologies. IMPORTANCE Human papillomaviruses are currently classified in different genera. Mucosal HPVs belonging to the alpha genus have been clearly associated with carcinogenesis of the mucosal epithelium at different sites. Beta HPV types have been classified as cutaneous. Although findings indicate that some beta HPVs from species 1 and 2 play a role, together with UV irradiation, in skin cancer, very little is known about the transforming properties of most of the beta HPVs. This report shows the transforming activity of E6 and E7 from beta-3 HPV types. Moreover, it highlights that beta-3 HPVs share some biological properties more extensively with mucosal high-risk HPV16 than with beta-2 HPV38. This report provides new paradigms for a better understanding of the biology of the different HPV types and their possible association with lesions at mucosal and/or cutaneous epithelia.

This review summarizes state-of-the-art knowledge in early-generation and novel urine biomarkers targeting the telomerase pathway for the detection and follow-up of bladder cancer (BC). The limitations of the assays detecting telomerase reactivation are discussed and the potential of transcription-activating mutations in the promoter of the TERT gene detected in the urine as promising simple non-invasive BC biomarkers is highlighted. Studies have shown good sensitivity and specificity of the urinary TERT promoter mutations in case-control studies and, more recently, in a pilot prospective cohort study, where the marker was detected up to 10 years prior to clinical diagnosis. However, large prospective cohort studies and intervention studies are required to fully validate their robustness and assess their clinical utility. Furthermore, it may be interesting to evaluate whether the clinical performance of urinary TERT promoter mutations could increase when combined with other simple urinary biomarkers. Finally, different approaches for assessment of TERT promoter mutations in urine samples are presented together with technical challenges, thus highlighting the need of careful technological validation and standardization of laboratory methods prior to translation into clinical practice.

Background: The DNA released into the bloodstream by malignant tumours· called circulating tumour DNA (ctDNA), is often a small fraction of total cell-free DNA shed predominantly by hematopoietic cells and is therefore challenging to detect. Understanding the biological properties of ctDNA is key to the investigation of its clinical relevance as a non-invasive marker for cancer detection and monitoring. Methods: We selected 40 plasma DNA samples of pancreatic cancer cases previously reported to carry a KRAS mutation at the ‘hotspot’ codon 12 and re-screened the cell-free DNA using a 4-size amplicons strategy (57 bp, 79 bp, 167 bp and 218 bp) combined with ultra-deep sequencing in order to investigate whether amplicon lengths could impact on the capacity of detection of ctDNA, which in turn could provide inference of ctDNA and non-malignant cell-free DNA size distribution. Findings: Higher KRAS amplicon size (167 bp and 218 bp) was associated with lower detectable cell-free DNA mutant allelic fractions (p

Previous studies have revealed a robust association between exposure to asbestos and human lung cancer. Accumulating evidence has highlighted the role of epigenome deregulation in the mechanism of carcinogen-induced malignancies. We examined the impact of asbestos on DNA methylation. Our genome-wide studies (using Illumina HumanMethylation450K BeadChip) of lung cancer tissue and paired normal lung from 28 asbestos-exposed or non-exposed patients, mostly smokers, revealed distinctive DNA methylation changes. We identified a number of differentially methylated regions (DMR) and differentially variable, differentially methylated CpGs (DVMC), with individual CpGs further validated by pyrosequencing in an independent series of 91 non-small cell lung cancer and paired normal lung. We discovered and validated BEND4, ZSCAN31 and GPR135 as significantly hypermethylated in lung cancer. DMRs in genes such as RARB (FDR 1.1 × 10-19 , mean change in beta [Δ] -0.09), GPR135 (FDR 1.87 × 10-8 , mean Δ -0.09) and TPO (FDR 8.58 × 10-5 , mean Δ -0.11), and DVMCs in NPTN, NRG2, GLT25D2 and TRPC3 (all with p

We report a novel CASP9 germline mutation that may increase susceptibility to the development of brain tumors. We identified this mutation in a family in which three brain tumors had developed within three generations, including two anaplastic astrocytomas occurring in cousins. The cousins were diagnosed at similar ages (29 and 31 years), and their tumors showed similar histological features. Genetic analysis revealed somatic IDH1 and TP53 mutations in both tumors. However, no germline TP53 mutations were detected, despite the fact that this family fulfills the criteria of Li-Fraumeni-like syndrome. Whole exome sequencing revealed a germline stop-gain mutation (R65X) in the CASP9 gene, which encodes caspase-9, a key molecule for the p53-dependent mitochondrial death pathway. This mutation was also detected in DNA extracted from blood samples from the two siblings who were each a parent of one of the affected cousins. Caspase-9 immunohistochemistry showed the absence of caspase-9 immunoreactivity in the anaplastic astrocytomas and normal brain tissues of the cousins. These observations suggest that CASP9 germline mutations may have played a role at least in part to the susceptibility of development of gliomas in this Li-Fraumeni-like family lacking a TP53 germline mutation.

Merkel cell polyomavirus (MCPyV) is the first human polyomavirus etiologically associated with Merkel cell carcinoma (MCC), a rare and aggressive form of skin cancer. Similar to other polyomaviruses, MCPyV encodes early T antigen genes, viral oncogenes required for MCC tumor growth. To identify the unique oncogenic properties of MCPyV, we analyzed the gene expression profiles in human spontaneously immortalized keratinocytes (NIKs) expressing the early genes from six distinct human polyomaviruses (PyVs), including MCPyV. A comparison of the gene expression profiles revealed 28 genes specifically deregulated by MCPyV. In particular, the MCPyV early gene downregulated the expression of the tumor suppressor gene N-myc downstream-regulated gene 1 (NDRG1) in MCPyV gene-expressing NIKs and hTERT-MCPyV gene-expressing human keratinocytes (HK) compared to their expression in the controls. In MCPyV-positive MCC cells, the expression of NDRG1 was downregulated by the MCPyV early gene, as T antigen knockdown rescued the level of NDRG1. In addition, NDRG1 overexpression in hTERT-MCPyV gene-expressing HK or MCC cells resulted in a decrease in the number of cells in S phase and cell proliferation inhibition. Moreover, a decrease in wound healing capacity in hTERT-MCPyV gene-expressing HK was observed. Further analysis revealed that NDRG1 exerts its biological effect in Merkel cell lines by regulating the expression of the cyclin-dependent kinase 2 (CDK2) and cyclin D1 proteins. Overall, NDRG1 plays an important role in MCPyV-induced cellular proliferation. IMPORTANCE Merkel cell carcinoma was first described in 1972 as a neuroendocrine tumor of skin, most cases of which were reported in 2008 to be caused by a PyV named Merkel cell polyomavirus (MCPyV), the first PyV linked to human cancer. Thereafter, numerous studies have been conducted to understand the etiology of this virus-induced carcinogenesis. However, it is still a new field, and much work is needed to understand the molecular pathogenesis of MCC. In the current work, we sought to identify the host genes specifically deregulated by MCPyV, as opposed to other PyVs, in order to better understand the relevance of the genes analyzed on the biological impact and progression of the disease. These findings open newer avenues for targeted drug therapies, thereby providing hope for the management of patients suffering from this highly aggressive cancer.

Renal cell carcinoma (RCC) has an undisputed genetic component and a stable 2:1 male to female sex ratio in its incidence across populations, suggesting possible sexual dimorphism in its genetic susceptibility. We conducted the first sex-specific genome-wide association analysis of RCC for men (3227 cases, 4916 controls) and women (1992 cases, 3095 controls) of European ancestry from two RCC genome-wide scans and replicated the top findings using an additional series of men (2261 cases, 5852 controls) and women (1399 cases, 1575 controls) from two independent cohorts of European origin. Our study confirmed sex-specific associations for two known RCC risk loci at 14q24.2 (DPF3) and 2p21(EPAS1). We also identified two additional suggestive male-specific loci at 6q24.3 (SAMD5, male odds ratio (OR(male)) = 0.83 [95% CI = 0.78-0.89], P(male) = 1.71 × 10(−8) compared with female odds ratio (OR(female)) = 0.98 [95% CI = 0.90–1.07], P(female) = 0.68) and 12q23.3 (intergenic, OR(male) = 0.75 [95% CI = 0.68-0.83], P(male) = 1.59 × 10(−8) compared with OR(female) = 0.93 [95% CI = 0.82–1.06], P(female) = 0.21) that attained genome-wide significance in the joint meta-analysis. Herein, we provide evidence of sex-specific associations in RCC genetic susceptibility and advocate the necessity of larger genetic and genomic studies to unravel the endogenous causes of sex bias in sexually dimorphic traits and diseases like RCC.

The utility of KRAS mutations in plasma circulating cell-free DNA (cfDNA) samples as non-invasive biomarkers for the detection of pancreatic cancer has never been evaluated in a large case-control series. We applied a KRAS amplicon-based deep sequencing strategy combined with analytical pipeline specifically designed for the detection of low-abundance mutations to screen plasma samples of 437 pancreatic cancer cases, 141 chronic pancreatitis subjects, and 394 healthy controls. We detected mutations in 21.1% (N=92) of cases, of whom 82 (89.1%) carried at least one mutation at hotspot codons 12, 13 or 61, with mutant allelic fractions from 0.08% to 79%. Advanced stages were associated with an increased proportion of detection, with KRAS cfDNA mutations detected in 10.3%, 17,5% and 33.3% of cases with local, regional and systemic stages, respectively. We also detected KRAS cfDNA mutations in 3.7% (N=14) of healthy controls and in 4.3% (N=6) of subjects with chronic pancreatitis, but at significantly lower allelic fractions than in cases. Combining cfDNA KRAS mutations and CA19-9 plasma levels on a limited set of case-control samples did not improve the overall performance of the biomarkers as compared to CA19-9 alone. Whether the limited sensitivity and specificity observed in our series of KRAS mutations in plasma cfDNA as biomarkers for pancreatic cancer detection are attributable to methodological limitations or to the biology of cfDNA should be further assessed in large case-control series.

Circulating tumor DNA (ctDNA) is emerging as a key potential biomarker for post-diagnosis surveillance but it may also play a crucial role in the detection of pre-clinical cancer. Small-cell lung cancer (SCLC) is an excellent candidate for early detection given there are no successful therapeutic options for late-stage disease, and it displays almost universal inactivation of TP53. We assessed the presence of TP53 mutations in the cell-free DNA (cfDNA) extracted from the plasma of 51 SCLC cases and 123 non-cancer controls. We identified mutations using a pipeline specifically designed to accurately detect variants at very low fractions. We detected TP53 mutations in the cfDNA of 49% SCLC patients and 11.4% of non-cancer controls. When stratifying the 51 initial SCLC cases by stage, TP53 mutations were detected in the cfDNA of 35.7% early-stage and 54.1% late-stage SCLC patients. The results in the controls were further replicated in 10.8% of an independent series of 102 non-cancer controls. The detection of TP53 mutations in 11% of the 225 non-cancer controls suggests that somatic mutations in cfDNA among individuals without any cancer diagnosis is a common occurrence, and poses serious challenges for the development of ctDNA screening tests., Highlights • CtDNA may play a crucial role in the detection of pre-clinical cancer. • TP53 mutations are detectable in the cfDNA of SCLC patients with early-stage tumors. • Detection of TP53 mutations in non-cancer controls poses serious challenges for the development of ctDNA screening tests. Cell-free DNA (cfDNA) has potential for monitoring response to treatment and relapse, but also for early detection. This is the first study reporting the detection of circulating-tumor DNA (ctDNA) in cases diagnosed with small-cell lung cancer (SCLC). Our results show that TP53 mutations are detectable in the cfDNA of SCLC patients including those with early-stage tumors. Importantly, we also provide evidence that cancer-like TP53 mutations are present in non-cancer controls, which poses serious challenges for the development of ctDNA screening tests.

Background: Detecting pre-clinical bladder cancer (BC) using urinary biomarkers may provide a valuable opportunity for screening and management. Telomerase reverse transcriptase (TERT) promoter mutations detectable in urine have emerged as promising BC biomarkers. Methods: We performed a nested case-control study within the population-based prospective Golestan Cohort Study (50,045 participants, followed up to 14 years) and assessed TERT promoter mutations in baseline urine samples from 38 asymptomatic individuals who subsequently developed primary BC and 152 matched controls using a Next-Generation Sequencing-based single-plex assay (UroMuTERT) and droplet digital PCR assays. Findings: Results were obtained for 30 cases and 101 controls. TERT promoter mutations were detected in 14 pre-clinical cases (sensitivity 46·67%) and none of the controls (specificity 100·00%). At an estimated BC cumulative incidence of 0·09% in the cohort, the positive and negative predictive values were 100·00% and 99·95% respectively. The mutant allelic fractions decreased with the time interval from urine collection until BC diagnosis (p = 0·033) but the mutations were detectable up to 10 years prior to clinical diagnosis. Interpretation: Our results provide the first evidence from a population-based prospective cohort study of the potential of urinary TERT promoter mutations as promising non-invasive biomarkers for early detection of BC. Further studies should validate this finding and assess their clinical utility in other longitudinal cohorts. Funding: French Cancer League, World Cancer Research Fund International, Cancer Research UK, Tehran University of Medical Sciences, the International Agency for Research on Cancer, and the U.S. National Cancer Institute. Keywords: Early detection, TERT promoter mutations, Urinary biomarker, Bladder cancer, Prospective cohort

BackgroundIn Latin America (LA), there is a high incidence rate of breast cancer (BC) in premenopausal women, and the genomic features of these BC remain unknown. Here, we aim to characterize the molecular features of BC in young LA women within the framework of the PRECAMA study, a multicenter population-based case-control study on breast cancer in premenopausal women.MethodsPathological tumor tissues were collected from incident cases from four LA countries. Immunohistochemistry (IHC) was performed centrally for ER, PR, HER2, Ki67, EGFR, CK5/6 and p53 protein markers. Targeted deep sequencing was done on genomic DNA extracted from formalin-fixed paraffin-embedded (FFPE) tumour tissues and their paired blood samples to screen for somatic mutations in eight genes frequently mutated in BC. A subset of samples was analyzed by exome sequencing to identify somatic mutational signatures.ResultsThe majority of cases were positive for ER or PR (168/233; 72%) and there were 21% triple negative (TN) cases, mainly of basal type. Most tumors were positive for Ki67 (189/233; 81%). In 126 sequenced cases,TP53andPIK3CAwere the most frequently mutated genes (32.5% and 21.4%, respectively), followed byAKT1(9.5%).TP53mutations were more frequent in HER2-enriched and TN IHC subtypes, whilePIK3CA/AKT1mutations were more frequent in ER positive tumors, as expected. Interestingly, a higher proportion of G:C>T:A mutations was observed inTP53gene in PRECAMA cases compared to TCGA and METABRIC breast cancer series (27% vs 14%). Exome-wide mutational patterns in 10 TN cases revealed alterations in signaling transduction pathways and major contributions of mutational signatures caused by altered DNA repair pathways.ConclusionsThis pilot results on PRECAMA tumors gives a preview of the molecular features of premenopausal BC in LA. Although, the overall mutation burden was as expected from data in other populations, mutational patterns observed inTP53and exome-wide suggested possible differences in mutagenic processes giving rise to these tumors compared to other populations. Further omics analyses of a larger number of cases in the near future will allow investigating relationships between these molecular features and risk factors.

We investigated how somatic changes in HNSCC interact with environmental and host risk factors and whether they influence the risk of HNSCC occurrence and outcome. 180-paired samples diagnosed as HNSCC in two high incidence regions of Europe and South America underwent targeted sequencing (14 genes) and evaluation of copy number alterations (SCNAs). TP53, PIK3CA, NOTCH1, TP63 and CDKN2A were the most frequently mutated genes. Cases were characterized by a low copy number burden with recurrent focal amplification in 11q13.3 and deletion in 15q22. Cases with low SCNAs showed an improved overall survival. We found significant correlations with decreased overall survival between focal amplified regions 4p16, 10q22 and 22q11, and losses in 12p12, 15q14 and 15q22. The mutational landscape in our cases showed an association to both environmental exposures and clinical characteristics. We confirmed that somatic copy number alterations are an important predictor of HNSCC overall survival.

// Pierre-Benoit Ancey 1 , Barbara Testoni 2 , Marion Gruffaz 2 , Marie-Pierre Cros 1 , Geoffroy Durand 3 , Florence Le Calvez-Kelm 3 , David Durantel 2 , Zdenko Herceg 1 , Hector Hernandez-Vargas 1 1 Epigenetics Group, International Agency for Research on Cancer (IARC), Lyon, France 2 INSERM U1052, Molecular Physiopathology and New Treatments of Viral Hepatitis, Centre de Recherche en Cancerologie (CRCL), Lyon, France 3 Genetic Cancer Susceptibility Group, International Agency for Research on Cancer (IARC), Lyon, France Correspondence to: Hector Hernandez-Vargas, e-mail: vargash@iarc.fr Keywords: HBV, hepatitis, methylome, HM450, transcriptome Received: August 17, 2015 Accepted: October 14, 2015 Published: November 02, 2015 ABSTRACT Viral infections are able to modify the host’s cellular programs, with DNA methylation being a biological intermediate in this process. The extent to which viral infections deregulate gene expression and DNA methylation is not fully understood. In the case of Hepatitis B virus (HBV), there is evidence for an interaction between viral proteins and the host DNA methylation machinery. We studied the ability of HBV to modify the host transcriptome and methylome, using naturally infected primary human hepatocytes to better mimic the clinical setting. Gene expression was especially sensitive to culture conditions, independently of HBV infection. However, we identified non-random changes in gene expression and DNA methylation occurring specifically upon HBV infection. There was little correlation between expression and methylation changes, with transcriptome being a more sensitive marker of time-dependent changes induced by HBV. In contrast, a set of differentially methylated sites appeared early and were stable across the time course experiment. Finally, HBV-induced DNA methylation changes were defined by a specific chromatin context characterized by CpG-poor regions outside of gene promoters. These data support the ability of HBV to modulate host cell expression and methylation programs. In addition, it may serve as a reference for studies addressing the genome-wide consequences of HBV infection in human hepatocytes.

// Sandra Perdomo 1, 2, * , Patrice H. Avogbe 1, * , Matthieu Foll 1 , Behnoush Abedi-Ardekani 1 , Violeta Lescher Facciolla 3 , Devasena Anantharaman 1,9 , Priscilia Chopard 1 , Florence Le Calvez-Kelm 1 , Marta Vilensky 4 , Jerry Polesel 5 , Ivana Holcatova 6 , Lorenzo Simonato 7 , Cristina Canova 7 , Pagona Lagiou 8 , James D. McKay 1 and Paul Brennan 1 1 International Agency for Research on Cancer (IARC), Lyon 69372, France 2 Institute of Nutrition, Genetics and Metabolism Research, Faculty of Medicine, Universidad El Bosque, Bogota 110121, Colombia 3 Departamento de Saude Coletiva, Faculdade de Ciencias Medicas da Santa Casa de Sao Paulo, Sao Paulo CEP01221-020, Brazil 4 Instituto Angel Roffo, Buenos Aires C1417DTB, Argentina 5 CRO Aviano National Cancer Institute, Aviano 33081, Italy 6 Charles University, 1st Faculty of Medicine, Prague 116 36, Czech Republic 7 Laboratory of Public Health and Population Studies, University of Padova, Padova 35122, Italy 8 Department of Hygiene, Epidemiology and Medical Statistics, School of Medicine, National and Kapodistrian University of Athens, Athens 11527, Greece 9 Current/present address: Rajiv Gandhi Centre for Biotechnology, Trivandrum 695010, India * These authors have contributed equally to this work Correspondence to: Paul Brennan, email: BrennanP@iarc.fr Keywords: ctDNA, head and neck, cancer, mutation, detection Abbreviations: cfDNA: cell free DNA; ctDNA: circulating tumor DNA; HNSCC: head and neck squamous cell carcinoma; AF: allelic frequency Received: April 27, 2017 Accepted: July 12, 2017 Published: August 07, 2017 ABSTRACT The use of non-invasive biomarkers such as circulating tumor DNA (ctDNA) in head and neck tumors may be of relevance in early diagnosis and eventually improved outcome. We evaluated two different approaches from two case series in Europe and South America including (i) targeted screening of ctDNA mutations, and (ii) detection of TP53 mutations in plasma and oral rinses without previous knowledge of mutational status in tumor samples. Targeted sequencing in 5 genes identified ctDNA mutations in plasma among 42% of HNSCC cases, 67% of who were early stage cases. No association was found between ctDNA mutation detection and overall survival. Sequencing of the entire coding region of the TP53 gene resulted in identification of TP53 mutations in 76% of tumor cases. However, concordance of mutation detection was low between tumor, oral rinses (11%) and plasma (2,7%) samples. Identification of 5 pathogenic TP53 mutations in oral rinses from 3 non-cancer controls gives additional evidence of mutation occurrence in individuals without a diagnosed cancer and presents an additional challenge for the development of ctDNA diagnostic assays.

Accepted manuscript version, licensed CC BY-NC-ND 4.0. Published version available in European Journal of Cancer, 75, 299-307. Aim of the study: A vast majority of human malignancies are associated with ageing, and age is a strong predictor of cancer risk. Recently, DNA methylation-based marker of ageing, known as ‘epigenetic clock’, has been linked with cancer risk factors. This study aimed to evaluate whether the epigenetic clock is associated with breast cancer risk susceptibility and to identify potential epigenetics-based biomarkers for risk stratification. Methods: Here, we profiled DNA methylation changes in a nested caseecontrol study embedded in the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort (n Z 960) using the Illumina HumanMethylation 450K BeadChip arrays and used the Horvath age estimation method to calculate epigenetic age for these samples. Intrinsic epigenetic age acceleration (IEAA) was estimated as the residuals by regressing epigenetic age on chronological age. Results: We observed an association between IEAA and breast cancer risk (OR, 1.04; 95% CI, 1.007e1.076, P Z 0.016). One unit increase in IEAA was associated with a 4% increased odds of developing breast cancer (OR, 1.04; 95% CI, 1.007e1.076). Stratified analysis based on menopausal status revealed that IEAA was associated with development of postmenopausal breast cancers (OR, 1.07; 95% CI, 1.020e1.11, P Z 0.003). In addition, methylome-wide analyses revealed that a higher mean DNA methylation at cytosine-phosphate-guanine (CpG) islands was associated with increased risk of breast cancer development (OR per 1 SD Z 1.20; 95 %CI: 1.03e1.40, P Z 0.02) whereas mean methylation levels at non-island CpGs were indistinguishable between cancer cases and controls. Conclusion: Epigenetic age acceleration and CpG island methylation have a weak, but statistically significant, association with breast cancer susceptibility.

To examine the diversity of somatic alterations and clonal evolution according to aggressiveness of disease, nineteen tumor-blood pairs of 'formerly bronchiolo-alveolar carcinoma (BAC)' which had been reclassified into preinvasive lesion (adenocarcinoma in situ; AIS), focal invasive lesion (minimally invasive adenocarcinoma; MIA), and invasive lesion (lepidic predominant adenocarcinoma; LPA and non-lepidic predominant adenocarcinoma; non-LPA) according to IASLC/ATS/ERS 2011 classification were explored by whole exome sequencing. Several distinct somatic alterations were observed compare to the lung adenocarcinoma study from the Cancer Genome Atlas (TCGA). There were higher numbers of tumors with significant APOBEC mutation fold enrichment (73% vs. 58% TCGA). The frequency of KRAS mutations was lower in our study (5% vs. 32% TCGA), while a higher number of mutations of RNA-splicing genes, RBM10 and U2AF1, were found (37% vs. 11% TCGA). We found neither mutational pattern nor somatic copy number alterations that were specific to AIS/MIA. We demonstrated that clonal cell fraction was the only distinctive feature that discriminated LPA/non-LPA from AIS/MIA. The broad range of clonal frequency signified a more branched clonal evolution at the time of diagnosis. Assessment of tumor clonal cell fraction might provide critical information for individualized therapy as a prognostic factor, however this needs further study.

In Brazil, breast cancer is a public health care problem due to its high incidence and mortality rates. In this study, we investigated the prevalence of hereditary breast cancer syndromes (HBCS) in a population-based cohort in Brazils southernmost capital, Porto Alegre. All participants answered a questionnaire about family history (FH) of breast, ovarian and colorectal cancer and those with a positive FH were invited for genetic cancer risk assessment (GCRA). If pedigree analysis was suggestive of HBCS, genetic testing of the BRCA1, BRCA2, TP53, and CHEK2 genes was offered. Of 902 women submitted to GCRA, 214 had pedigrees suggestive of HBCS. Fifty of them underwent genetic testing: 18 and 40 for BRCA1/BRCA2 and TP53 mutation screening, respectively, and 7 for CHEK2 1100delC testing. A deleterious BRCA2 mutation was identified in one of the HBOC probands and the CHEK2 1100delC mutation occurred in one of the HBCC families. No deleterious germline alterations were identified in BRCA1 or TP53. Although strict inclusion criteria and a comprehensive testing approach were used, the suspected genetic risk in these families remains unexplained. Further studies in a larger cohort are necessary to better understand the genetic component of hereditary breast cancer in Southern Brazil.

International audience; Growing evidence supports the existence of a subpopulation of cancer cells with stem cell characteristics within breast tumors. In spite of its potential clinical implications, an understanding of the mechanisms responsible for retaining the stem cell characteristics in these cells is still lacking. Here, we used the mammosphere model combined with DNA methylation bead arrays and quantitative gene expression to characterize the epigenetic mechanisms involved in the regulation of developmental pathways in putative breast cancer stem cells. Our results revealed that MCF7-derived mammospheres exhibit distinct CpG promoter methylation profiles in a specific set of genes, including those involved in Jak-STAT signaling pathway. Hypomethylation of several gene components of the Jak-STAT pathway was correlated with an increased expression in mammospheres relative to parental cells. Remarkably, cell sorting of the cells with a putative cancer stem cell phenotype (CD44+/CD24 low) suggests a constitutive activation of Jak-STAT pathway in these cells. These results show that Jak-STAT activation may represent a characteristic of putative breast cancer stem cells. In addition, they favor the concept that the expression of cancer stem-like pathways and the establishment and maintenance of defining properties of cancer stem cells are orchestrated by epigenetic mechanisms.

Schwannomas are benign nerve sheath tumors composed of well-differentiated Schwann cells. Other than frequent NF2 (neurofibromatosis type 2) mutations (50%-60%), their molecular pathogenesis is not fully understood. LATS1 and LATS2 are downstream molecules of NF2 and are negative regulators of the yes-associated protein (YAP) oncogene in the Hippo signaling pathway. We assessed mutations of the NF2, LATS1, and LATS2 genes, promoter methylation of LATS1 and LATS2, and expression of YAP and phosphorylated YAP in 82 cases of sporadic schwannomas. Targeted sequencing using the Ion Torrent Proton instrument revealed NF2 mutations in 45 cases (55%), LATS1 mutations in 2 cases (2%), and LATS2 mutations in 1 case (1%) of schwannoma. Methylation-specific polymerase chain reaction showed promoter methylation of LATS1 and LATS2 in 14 cases (17%) and 25 cases (30%), respectively. Overall, 62 cases (76%) had at least 1 alteration in the NF2, LATS1, and/or LATS2 genes. Immunohistochemistry revealed nuclear YAP expression in 18 of 42 cases of schwannoma (43%) and reduced cytoplasmic phosphorylated YAP expression in 15 of 49 cases of schwannoma (31%), all of which had at least 1 alteration in the NF2, LATS1, and/or LATS2 genes. These results suggest that an abnormal Hippo signaling pathway is involved in the pathogenesis of most sporadic schwannomas.

Background\ud Exposure to environmental toxins during embryonic development may lead to epigenetic changes that influence disease risk in later life. Aflatoxin is a contaminant of staple foods in sub-Saharan Africa, is a known human liver carcinogen and has been associated with stunting in infants.\ud Methods\ud We have measured aflatoxin exposure in 115 pregnant women in The Gambia and examined the DNA methylation status of white blood cells from their infants at 2-8 months old (mean 3.6 ± 0.9). Aflatoxin exposure in women was assessed using an ELISA method to measure aflatoxin albumin (AF-alb) adducts in plasma taken at 1-16 weeks of pregnancy. Genome-wide DNA methylation of infant white blood cells was measured using the Illumina Infinium HumanMethylation450beadchip.\ud Results\ud AF-alb levels ranged from 3.9 to 458.4 pg/mg albumin. We found that aflatoxin exposure in the mothers was associated to DNA methylation in their infants for 71 CpG sites (false discovery rate < 0.05), with an average effect size of 1.7% change in methylation. Aflatoxin-associated differential methylation was observed in growth factor genes such as FGF12 and IGF1, and immune-related genes such as CCL28, TLR2 and TGFBI. Moreover, one aflatoxin-associated methylation region (corresponding to the miR-4520b locus) was identified.\ud Conclusions\ud This study shows that maternal exposure to aflatoxin during the early stages of pregnancy is associated with differential DNA methylation patterns of infants, including in genes related to growth and immune function. This reinforces the need for interventions to reduce aflatoxin exposure, especially during critical periods of fetal and infant development.

Background: Dietary exposure to cytotoxic and carcinogenic aristolochic acid (AA) causes severe nephropathy typically associated with urologic cancers. Monitoring of AA exposure uses biomarkers such as aristolactam-DNA adducts, detected by mass spectrometry in the kidney cortex, or the somatic A>T transversion pattern characteristic of exposure to AA, as revealed by previous DNA-sequencing studies using fresh-frozen tumors. Methods: Here, we report a low-coverage whole-exome sequencing method (LC-WES) optimized for multisample detection of the AA mutational signature, and demonstrate its utility in 17 formalin-fixed paraffin-embedded urothelial tumors obtained from 15 patients with endemic nephropathy, an environmental form of AA nephropathy. Results: LC-WES identified the AA signature, alongside signatures of age and APOBEC enzyme activity, in 15 samples sequenced at the average per-base coverage of approximately 10×. Analysis at 3 to 9× coverage revealed the signature in 91% of the positive samples. The exome-wide distribution of the predominant A>T transversions exhibited a stochastic pattern, whereas 83 cancer driver genes were enriched for recurrent nonsynonymous A>T mutations. In two patients, pairs of tumors from different parts of the urinary tract, including the bladder, harbored overlapping mutation patterns, suggesting tumor dissemination via cell seeding. Conclusions: LC-WES analysis of archived tumor tissues is a reliable method applicable to investigations of both the exposure to AA and its biologic effects in human carcinomas. Impact: By detecting cancers associated with AA exposure in high-risk populations, LC-WES can support future molecular epidemiology studies and provide evidence-base for relevant preventive measures. Cancer Epidemiol Biomarkers Prev; 24(12); 1873–81. ©2015 AACR.

Nat. Genet. 46, 736–741 (2014); published online 1 June 2014; corrected after print 23 January 2017 In the version of this article initially published, the name of author Florence Le Calvez-Kelm appeared incorrectly as Florence LeCalvez-Kelm. The error has been corrected in the HTML and PDF versionsof the article.

Approximately half of the familial aggregation of breast cancer remains unexplained. A multiple-case breast cancer family exome-sequencing study identified three likely pathogenic mutations in RINT1 (NM_021930.4) not present in public sequencing databases: RINT1 c.343C>T (p.Q115X), c.1132_1134del (p.M378del), and c.1207G>T (p.D403Y). On the basis of this finding, a population-based case–control mutation-screening study was conducted that identified 29 carriers of rare (minor allele frequency < 0.5%), likely pathogenic variants: 23 in 1,313 early-onset breast cancer cases and six in 1,123 frequency-matched controls [OR, 3.24; 95% confidence interval (CI), 1.29–8.17; P = 0.013]. RINT1 mutation screening of probands from 798 multiple-case breast cancer families identified four additional carriers of rare genetic variants. Analysis of the incidence of first primary cancers in families of women carrying RINT1 mutations estimated that carriers were at increased risk of Lynch syndrome–spectrum cancers [standardized incidence ratio (SIR), 3.35; 95% CI, 1.7–6.0; P = 0.005], particularly for relatives diagnosed with cancer under the age of 60 years (SIR, 10.9; 95% CI, 4.7–21; P = 0.0003). Significance: The work described in this study adds RINT1 to the growing list of genes in which rare sequence variants are associated with intermediate levels of breast cancer risk. Given that RINT1 is also associated with a spectrum of cancers with mismatch repair defects, these findings have clinical applications and raise interesting biological questions. Cancer Discov; 4(7); 804–15. ©2014 AACR. See related commentary by Ngeow and Eng, p. 762 This article is highlighted in the In This Issue feature, p. 745

Purpose: Tumor somatic mutation analysis is part of the standard management of metastatic lung cancer. However, physicians often have to deal with small biopsies and consequently with challenging mutation testing. Circulating free DNA (cfDNA) is a promising tool for accessing the tumor genome as a liquid biopsy. Here, we evaluated next-generation sequencing (NGS) on cfDNA samples obtained from a consecutive series of patients for the screening of a range of clinically relevant mutations. Experimental Design: A total of 107 plasma samples were collected from the BioCAST/IFCT-1002 lung cancer study (never-smokers cohort). Matched tumor DNA (tDNA) was obtained for 68 cases. Multiplex PCR-based assays were designed to target specific coding regions in EGFR, KRAS, BRAF, ERBB2, and PI3KCA genes, and amplicon sequencing was performed at deep coverage on the cfDNA/tDNA pairs using the NGS IonTorrent Personal Genome Machine Platform. Results: CfDNA concentration in plasma was significantly associated with both stage and number of metastatic sites. In tDNA, 50 mutations (36 EGFR, 5 ERBB2, 4 KRAS, 3 BRAF, and 2 PIK3CA) were identified, of which 26 were detected in cfDNA. Sensitivity of the test was 58% (95% confidence interval, 43%–71%) and the estimated specificity was 87% (62%–96%). Conclusion: These data demonstrate the feasibility and potential utility of mutation screening in cfDNA using IonTorrent NGS for the detection of a range of tumor biomarkers in patients with metastatic lung cancer. Clin Cancer Res; 20(17); 4613–24. ©2014 AACR.

Background Distinct subpopulations of neoplastic cells within tumors, including hepatocellular carcinoma (HCC), display pronounced ability to initiate new tumors and induce metastasis. Recent evidence suggests that signals from transforming growth factor beta (TGF-β) may increase the survival of these so called tumor initiating cells leading to poor HCC prognosis. However, how TGF-β establishes and modifies the key features of these cell subpopulations is not fully understood. Results In the present report we describe the differential DNA methylome of CD133-negative and CD133-expressing liver cancer cells. Next, we show that TGF-β is able to increase the proportion of CD133+ cells in liver cancer cell lines in a way that is stable and persistent across cell division. This process is associated with stable genome-wide changes in DNA methylation that persist through cell division. Differential methylation in response to TGF-β is under-represented at promoter CpG islands and enriched at gene bodies, including a locus in the body of the de novo DNA methyl-transferase DNMT3B gene. Moreover, phenotypic changes induced by TGF-β, including the induction of CD133, are impaired by siRNA silencing of de novo DNA methyl-transferases. Conclusions Our study reveals a self-perpetuating crosstalk between TGF-β signaling and the DNA methylation machinery, which can be relevant in the establishment of cellular phenotypes. This is the first indication of the ability of TGF-β to induce genome-wide changes in DNA methylation, resulting in a stable change in the proportion of liver cancer cell subpopulations. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-435) contains supplementary material, which is available to authorized users.

Exome DNA sequencing of blood samples from a Li-Fraumeni family with a TP53 germline mutation (codon 236 deletion) and multiple nervous system tumors revealed additional germline mutations. Missense mutations in the MSH4 DNA repair gene (c.2480T>A; p.I827N) were detected in three patients with gliomas (two anaplastic astrocytomas, two glioblastomas). Two family members without a TP53 germline mutation who developed peripheral schwannomas also carried the MSH4 germline mutation, and in addition, a germline mutation of the LATS1 gene (c.286C>T; p.R96W). LATS1 is a downstream mediator of the NF2 , but has not previously been found to be related to schwannomas. We therefore screened the entire coding sequence of the LATS1 gene in 65 sporadic schwannomas, 12 neurofibroma/schwannoma hybrid tumors, and 4 cases of schwannomatosis. We only found a single base deletion at codon 827 (exon 5) in a spinal schwannoma, leading to a stop at codon 835 (c.2480delG; p.*R827Kfs*8). Mutational loss of LATS1 function may thus play a role in some inherited schwannomas, but only exceptionally in sporadic schwannomas. This is the first study reporting a germline MSH4 mutation. Since it was present in all patients, it may have contributed to the subsequent acquisition of TP53 and LATS1 germline mutations.

Background A subset of breast cancer cells displays increased ability to self-renew and reproduce breast cancer heterogeneity. The characterization of these so-called putative breast tumor-initiating cells (BT-ICs) may open the road for novel therapeutic strategies. As microRNAs (miRNAs) control developmental programs in stem cells, BT-ICs may also rely on specific miRNA profiles for their sustained activity. To explore the notion that miRNAs may have a role in sustaining BT-ICs, we performed a comprehensive profiling of miRNA expression in a model of putative BT-ICs enriched by non-attachment growth conditions. Results We found breast cancer cells grown under non-attachment conditions display a unique pattern of miRNA expression, highlighted by a marked low expression of miR-30 family members relative to parental cells. We further show that miR-30a regulates non-attachment growth. A target screening revealed that miR-30 family redundantly modulates the expression of apoptosis and proliferation-related genes. At least one of these targets, the anti-apoptotic protein AVEN, was able to partially revert the effect of miR-30a overexpression. Finally, overexpression of miR-30a in vivo was associated with reduced breast tumor progression. Conclusions miR30-family regulates the growth of breast cancer cells in non-attachment conditions. This is the first analysis of target prediction in a whole family of microRNAs potentially involved in survival of putative BT-ICs.

The quality of tissue samples and extracted mRNA is a major source of variability in tumor transcriptome analysis using genome-wide expression microarrays. During and immediately after surgical tumor resection, tissues are exposed to metabolic, biochemical and physical stresses characterized as "warm ischemia". Current practice advocates cryopreservation of biosamples within 30 minutes of resection, but this recommendation has not been systematically validated by measurements of mRNA decay over time. Using Illumina HumanHT-12 v3 Expression BeadChips, providing a genome-wide coverage of over 24,000 genes, we have analyzed gene expression variation in samples of 3 hepatocellular carcinomas (HCC) and 3 lung carcinomas (LC) cryopreserved at times up to 2 hours after resection. RNA Integrity Numbers (RIN) revealed no significant deterioration of mRNA up to 2 hours after resection. Genome-wide transcriptome analysis detected non-significant gene expression variations of -3.5%/hr (95% CI: -7.0%/hr to 0.1%/hr; p = 0.054). In LC, no consistent gene expression pattern was detected in relation with warm ischemia. In HCC, a signature of 6 up-regulated genes (CYP2E1, IGLL1, CABYR, CLDN2, NQO1, SCL13A5) and 6 down-regulated genes (MT1G, MT1H, MT1E, MT1F, HABP2, SPINK1) was identified (FDR

Background: Pancreatic ductal adenocarcinoma (PDAC) is characterized by hotspot mutations in the KRAS gene (codons 12, 13 or 61) in 85-90% of cases. Codon 12 KRAS mutations have been detected in pancreatic juice, blood and stool samples from pancreatic cancer cases and represent promising biomarkers for early detection. However, the proportion of tumor-derived KRAS mutations in cell-free DNA fragments (cfDNA) has shown large variations, probably because of the heterogeneity in biosamples and assays tested. Deep sequencing technologies (NGS) allow the identification of low-abundance somatic variants, but have not previously been applied to the detection of KRAS hotspot mutations in cfDNA of PDAC cases. Moreover, variant calling methods have rarely been tested against cancer-free individuals so the proportion of false positives is unknown. We investigated whether deep sequencing of KRAS mutations at codons 12, 13 and 61 in plasma samples could represent a robust assay to distinguish pancreatic cancer from chronic pancreatitis and healthy controls. Methods: We developed an Ion Torrent-based NGS KRAS assay (partial exons 2 and 3, totalling 208bp) to screen cfDNA from plasma samples of 461 PDAC cases, 154 individuals with chronic pancreatitis and 421 healthy controls. cfDNA extraction (>4ng) and sequencing (>1000x coverage on average, and absence of systematic high sequencing error rates on the 208bp) performed well on 431 (93%) PDAC cases, 138 (90%) chronic pancreatitis, and 388 (95%) controls. We fit a robust negative-binomial regression to estimate the distribution of the sequencing errors at each DNA bp position and identified outlying samples, which were considered as KRAS positive when q-value Results: Sequencing of the serial dilutions of KRAS p.G12V mutated DNA indicated a detection threshold at a minor allele frequency of 0.2%. KRAS mutations in cfDNA were detected in 83 (19.3%) PDAC cases (73 on codon 12; 8 on codon 61; 1 on codon 13; and 1 multiple codons, i.e., similar in proportions as reported in tumor tissue from the International Cancer Genomic Consortium); 3 (2.2%) chronic pancreatitis (all on codon 12); and 8 (2.1%) healthy controls (4 on codon 12 and 4 on codon 61). Stage was significantly associated with the proportion of detected mutations in cancer cases (chi-squared p = 0.0005): the proportions of cases with detectable KRAS mutations in plasma were 7.9%, 14.9%, and 31.1% for local, regional, and advanced stages, respectively. Conclusions: The NGS-based KRAS mutation screening is a sensitive approach to detect low allelic fraction in plasma cfDNA, although its utility for early detection is still limited. However, it has the capacity to identify specific KRAS mutations, which could be useful in a panel of other non-invasive biomarkers. Citation Format: Florence Le Calvez-Kelm, Matthieu Foll, Magdalena B. Wozniak, Geoffroy Durand, Priscilia Chopard, Maroulio Pertesi, Tiffany Delhomme, Ivana Holcatova, Lenka Foretova, Vladimir Janout, Eleonora Fabianova, Maxime P. Vallée, Paul Brennan, James D. McKay, Graham Byrnes, Ghislaine Scélo. NGS-based detection of KRAS hotspot mutations in plasma cell-free DNA of pancreatic cancer cases. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3137.

Background The US National Lung Cancer Screening Trial (NLST) demonstrated in 2011 that screening with computed tomography (CT) scans could reduce lung cancer mortality by 20%, but with important financial costs and high number of false positives. The identification of novel biomarkers is a need to obtain the maximum benefit from CT screening. Given its economical and minimally invasive nature, screening for somatic mutations in circulating tumor DNA (ctDNA) using next-generation sequencing may complement existing screening tools. However, for application in early detection, variant detection must also be done agnostically, i.e. without prior knowledge from tumour tissue of the mutations expected and unfortunately, most currently available variant callers are not adapted for this task. Methods We performed multiplex PCR on circulating free DNA (cfDNA) extracted from the plasma of 35 lung squamous cell carcinoma (SCC) and 64 small-cell carcinoma (SCLC) patients. We additionally included 133 hospital controls to evaluate the specificity of ctDNA. We applied (>10,000X) Ion torrent targeted sequencing on the full-coding region of TP53 since this gene is known to be mutated in more than 70% and 90% of SCC and SCLC, respectively. Each amplification, library preparation, and sequencing was performed in duplicate to control for amplification and sequencing errors. Detecting mutations on ctDNA raises important statistical and bioinformatics challenges as it represents only a small fraction of cfDNA. We therefore developed and applied a method based on the idea that a data-derived model of sequencing errors has the potential to improve our ability to detect low-allelic fraction (AF) somatic variants. Results We detected TP53 non-synonymous coding mutations with AFs between 0.04% and 85% (median 1.7%) in 8 (23%) SCC patients, 28 (44%) SCLC patients, and 8 (6%) controls. We estimated odds ratios of 4.6 (p = 0.006) for SCC and 12.0 (p = 6.7×10-10) for SCLC. Observations in controls are surprising, but in this instance there was no information regarding a subsequent cancer diagnosis. Conclusion We show that it is possible to detect ctDNA in the cfDNA of lung cancer patients. Since only patients with early stage (I-IIA) SCC tumours were included, these results support the potential utility of the approach for early detection. Nevertheless, if such mutations are found prior to diagnosis has not been explored in a prospective study design with pre-diagnostic plasma samples and individuals without a cancer diagnosis through a follow-up period. Citation Format: Patrice H. Avogbe, Tiffany Delhomme, Noémie Leblay, Florence Le Calvez-Kelm, Priscilia Chopard, Valérie Gaborieau, Ghislaine Scelo, Behnoush Abedi-Ardekani, David Zaridze, Anush Mukeria, Graham Byrnes, Paul Brennan, Lynnette Fernandez-Cuesta, Matthieu Foll, James D. McKay. NGS-based screening for TP53 mutations in circulating cell-free DNA: A first step towards early detection of lung cancers. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3156.

Approximately half of the familial aggregation of breast cancer remains unexplained. This proportion is less for early-onset disease where familial aggregation is greater, suggesting that other susceptibility genes remain to be discovered. The majority of known breast cancer susceptibility genes are involved in the DNA double-strand break repair pathway. ABRAXAS is involved in this pathway and mutations in this gene impair BRCA1 recruitment to DNA damage foci and increase cell sensitivity to ionizing radiation. Moreover, a recurrent germline mutation was reported in Finnish high-risk breast cancer families. To determine if ABRAXAS could be a breast cancer susceptibility gene in other populations, we conducted a population-based case-control mutation screening study of the coding exons and exon/intron boundaries of ABRAXAS in the Breast Cancer Family Registry. In addition to the common variant p.Asp373Asn, sixteen distinct rare variants were identified. Although no significant difference in allele frequencies between cases and controls was observed for the identified variants, two variants, p.Gly39Val and p.Thr141Ile, were shown to diminish phosphorylation of gamma-H2AX in MCF7 human breast adenocarcinoma cells, an important biomarker of DNA double-strand breaks. Overall, likely damaging or neutral variants were evenly represented among cases and controls suggesting that rare variants in ABRAXAS may explain only a small proportion of hereditary breast cancer.

Hepatocellular carcinoma (HCC) may develop according to two major pathways, one involving HBV infection and TP53 mutation and the other characterized by HCV infection and CTNNB1 mutation. We have investigated HBV/HCV infections and TP53/CTNNB1 mutations in 26 HCC patients from Thailand. HBV DNA (genotype B or C) was detected in 19 (73%) of the cases, including 5 occult infections and 3 coinfections with HCV. TP53 and CTNNB1 mutations were not mutually exclusive, and most of TP53 mutations were R249S, suggesting a significant impact of aflatoxin-induced mutagenesis in HCC development.

BackgroundHepatocellular carcinoma (HCC) is characterized by late detection and fast progression, and it is believed that epigenetic disruption may be the cause of its molecular and clinicopathological heterogeneity. A better understanding of the global deregulation of methylation states and how they correlate with disease progression will aid in the design of strategies for earlier detection and better therapeutic decisions.Methods and findingsWe characterized the changes in promoter methylation in a series of 30 HCC tumors and their respective surrounding tissue and identified methylation signatures associated with major risk factors and clinical correlates. A wide panel of cancer-related gene promoters was analyzed using Illumina bead array technology, and CpG sites were then selected according to their ability to classify clinicopathological parameters. An independent series of HCC tumors and matched surrounding tissue was used for validation of the signatures. We were able to develop and validate a signature of methylation in HCC. This signature distinguished HCC from surrounding tissue and from other tumor types, and was independent of risk factors. However, aberrant methylation of an independent subset of promoters was associated with tumor progression and etiological risk factors (HBV or HCV infection and alcohol consumption). Interestingly, distinct methylation of an independent panel of gene promoters was strongly correlated with survival after cancer therapy.ConclusionOur study shows that HCC tumors exhibit specific DNA methylation signatures associated with major risk factors and tumor progression stage, with potential clinical applications in diagnosis and prognosis.

Hepatocellular carcinoma (HCC) is the second most common cause of death from cancer worldwide, estimated to be responsible for nearly 746,000 deaths in 2012 representing 9.1% of the total cancer-related deaths. Incidence rates of HCC continue to increase in North America and Europe with a slight decrease in Asia. It is currently accepted that epigenetic alterations are critical in cancer development in humans, including HCC. Among the different layers of epigenetic control, long non coding RNAs (lncRNAs) have been steadily gaining prominence in recent years. LncRNAs have emerged as critical regulators of the gene expression and growing evidence points to their tissues and tumour specificity and that their deregulation may represent early event in tumour development. Several recent reports have implicated lncRNAs in HCC development, such as HULC or H19, although the biological and molecular mechanisms by which lncRNAs contributes to liver cancer development remain poorly understood.. The purpose of this study is to perform a global profiling of lncRNAs in HCC using Next Generation Sequencing in order to improve the understanding of lncRNAs’ role in HCC and to provide new biologically sound candidate biomarkers for the identification of patients at higher risk of developping HCC. We have performed RNA-Sequencing (RNA-Seq) of tumor tissues and paired-matched surrounding non-tumor (cirrhotic) tissues from 10 patients and normal livers from 3 patients operated for liver metastases. Based on our first results, we observed a global decrease in the number of lncRNAs expressed in cirrhotic and tumor tissues compared with normal livers. We identified 41 lncRNAs differentially expressed in HCC compared to matched non-tumour tissues, most of whom have not been reported previously. We identified also more than 700 coding genes differentially expressed in HCC compared to their non-tumor counterparts. We compared our results on the coding trascriptome with RNA-Seq data from 44 cases of HCC (paired-matched tumor and non tumor tissues) available in The Cancer Genome Atlas (TCGA) dataset and found that more than 70% of the genes identified in our study overlap with differentially expressed genes identified in the TCGA dataset. Our results suggest that lncRNAs may be involved in the homeostasis of the hepatic tissue and that the loss of their expression is associated with tissue regeneration and transformation. Citation Format: Davide Degli Esposti, Catherine Voegele, Vincent Cahais, Nathalie Forey, Florence Le Calvez Kelm, James McKay, Hector Hernandez-Vargas, Zdenko Herceg. Deciphering the role of long noncoding RNAs in the progression of hepatocellular carcinoma using next-generation sequencing. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 171. doi:10.1158/1538-7445.AM2015-171

Introduction and objectives. Epigenetic modifications such as changes in DNA methylation have been shown as response to many toxic exposures. However, not much is known about the role of epigenetics in human lung carcinogenesis associated with asbestos exposure. Here, we explored the epigenome-wide DNA methylation (DNA methylome) patterns in lung tumors focused on asbestos-exposed and non-exposed patients who were mostly smokers. We sought to investigate whether methylation markers in lung cancer associate with exposure (asbestos, tobacco smoke) and clinico-pathological features. Study design, subjects, and methods. The DNA methylome was examined by Illumina Infinium 450K HumanMethylation technology covering more than 450 000 CpG sites. Validation experiments were done by pyrosequencing. DNA methylation was studied in lung tumors and normal lung as matched pairs obtained from 28 patients from lung surgery in Helsinki during years 1988-1996. Detailed data on the patients’ work and smoking histories, pulmonary fiber counts (14 asbestos-exposed and 14 non-exposed), and clinical data were available. As controls, normal lung specimens from patients with non-malignant lung disease were studied for genome-wide DNA methylation. Results. DNA methylation profiles differed in relation to sample type, tumor histology and exposure status. We identified more than 10.000 methylation variable positions (MVPs) that differed (adjusted P < 0.05) between lung tumors and matched normal lung. Differential methylation in 14 MVPs were significantly associated with asbestos exposure in lung tumor tissue whereas more than 100 MVPs were suggested as smoking-related. Differential methylation of genes such as those involved in cell stress and inflammation, cellular energy and lipid metabolism, and protein degradation processes were observed in lung tumors. Conclusions. Our analysis suggests that methylation profiles were globally affected in lung tumors in comparison to matched normal lung. We also found associations to asbestos exposure and tobacco smoking, suggesting a role for these exposures in modification of DNA methylome profiles. Citation Format: Eeva M. Kettunen, Hector Hernandez-Vargas, Marie-Pierre Cros, Geoffroy Durand, Florence Le Calvez-Kelm, Henrik Wolff, Kaisa Salmenkivi, Sisko L. Anttila, Zdenko Herceg, Kirsti Husgafvel-Pursiainen. Epigenome-wide DNA methylation profiling of lung tumors from asbestos-exposed patients. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1060. doi:10.1158/1538-7445.AM2015-1060

Triple negative breast cancer (TNBC), defined by the lack of expression of the estrogen receptor, progesterone receptor and human epidermal receptor 2, is an aggressive form of breast cancer that is more prevalent in certain populations, in particular in low- and middle-income regions. The detailed molecular features of TNBC in these regions remain unexplored as samples are mostly accessible as formalin-fixed paraffin embedded (FFPE) archived tissues, a challenging material for advanced genomic and transcriptomic studies. Using dedicated reagents and analysis pipelines, we performed whole exome sequencing and miRNA and mRNA profiling of 12 FFPE tumor tissues collected from pathological archives in Mexico. Sequencing analyses of the tumor tissues and their blood pairs identified TP53 and RB1 genes as the most frequently mutated genes, with a somatic mutation load of 1.7 mutations/exome Mb on average. Transcriptional analyses revealed an overexpression of growth-promoting signals (EGFR, PDGFR, VEGF, PIK3CA, FOXM1), a repression of cell cycle control pathways (TP53, RB1), a deregulation of DNA-repair pathways, and alterations in epigenetic modifiers through miRNA:mRNA network de-regulation. The molecular programs identified were typical of those described in basal-like tumors in other populations. This work demonstrates the feasibility of using archived clinical samples for advanced integrated genomics analyses. It thus opens up opportunities for investigating molecular features of tumors from regions where only FFPE tissues are available, allowing retrospective studies on the search for treatment strategies or on the exploration of the geographic diversity of breast cancer.

Dietary intake of Aristolochia plants containing aristolochic acid (AA, IARC Group 1 carcinogen), leads to renal disease and to recurrent urinary tract transitional cell carcinomas (TCC). Mutation screens identified frequent A>T transversions, alterations underrepresented in other cancer types, as the signature of AA exposure. AA-associated TCC is a preventable cancer yet millions are estimated at risk due to unregulated medicinal use of Aristolochia. Here we present an ultra-low coverage modification of whole-exome sequencing (ULC-WES) as a novel molecular epidemiology method that reliably detects the AA mutation signature in small amounts of DNA from FFPE tumor tissues and can identify new exposed populations. DNA was isolated from paraffin sections of 18 TCC tumors from endemic nephropathy (EN) patients from regions in Croatia and Bosnia linked to exposure to AA. Four non-EN upper tract TCCs were used as controls. 250 ng of DNA was used for sequencing libraries and exome enrichment (Nimblegen, Roche). Libraries were sequenced on Illumina HiSeq2500 at 10-12x coverage. After data processing, GATK was used to call single nucleotide variants (SNV). To assess platform independence, 7 samples were run on SOLID 5500XL (Life), using proprietary chemistry and pipeline for data processing and variant calling. Normal genetic variants were filtered out after SNV annotation. ULC-WES at 10x identified the A>T predominance in 17 of 18 EN tumors (94%), and at T pattern in 14 (78%) samples. In contrast, the A>T presence was T mutations in the genome (R^2=0.83 for mutation count correlated with the chromosome length) but also identified recurrent non-synonymous A>T mutations in tumor suppressors and oncogenes including TP53, ARID1A, BRCA1/2, AXIN1, GATA3 and NRAS, DNMT1, FLT3, CSF1R and PIK3CA. We next detected recurrent A>T mutations in a number of histone and DNA modifiers suggesting epigenetic deregulation of genome activity in the AA-associated TCCs. In 2 patients, pairs of distinct TCCs from the same surgery had markedly overlapping mutation patterns (624 SNVs shared between same-side pelvic and ureteric tumors in one case, and 191 shared between renal pelvic and bladder tumors in the other). Thus, ULC-WES provided evidence supporting the concept of tumor cell seeding along the urinary tract as a mechanism of TCC spread and recurrence. We conclude that ULC-WES is a sensitive and specific method to identify AA exposure in tumor tissue and is directly applicable to identification of new populations at risk. The method can ultimately facilitate evidence-based preventive measures in the affected regions. Acknowledgments: Genetic Services Platform, IARC; Genome Technology Center, NYU Langone Medical Center; NIH P01 ES04068, NIH Fogarty International Center 1R03TW007042, Croatian Ministry of Science. Citation Format: Xavier Castells, Sandra Karanovic, Magali Olivier, Maude Ardin, Florence Le Calvez-Kelm, Catherine Voegele, James McKay, Damir Dittrich, Karla Medak Tomic, Kathleen Dickman, Arthur P. Grollman, Bojan Jelakovic, Jiri Zavadil. Ultra-low coverage exome sequencing of FFPE tumor specimens identifies exposure to carcinogenic aristolochic acid. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 305. doi:10.1158/1538-7445.AM2014-305

Aberrant DNA methylation (DNAm) in blood has been used as an epigenetic marker of cancer risk although, the significance of DNAm changes in surrogate tissues as a marker of predicting breast cancer risk remains elusive. This may be potentially accounted by the fact that there are very few study cohorts to understand the significance of DNA methyation prospectively. We examined the potential of DNAm in peripheral blood as a marker of breast cancer risk in a nested case-control cohort of European Prospective Investigation into Cancer and Nutrition (EPIC) study. The study cohort involved 451 female incident breast cancer cases and matched controls (matched for center, age, date of blood collection and menopausal status). Majority of breast cancer cases were localized (47%) and receptor positive (83.1%) while a smaller proportion of samples were metastatic (14%) and 20% of the cases were triple negative. Methylation levels were measured across 485,777 CpG sites using the Illumina HumanMethylation 450K BeadChip array. The cross-reactive probes and polymorphic CpGs were filtered during the analysis. Methylation values for all the probes thus obtained were compared between cases and matched controls using logistic regression model with adjustments for potential confounders. We found differential methylation across 185 probes (P Citation Format: Srikant Ambatipudi, Veronique Chajes, Florence Le Calvez-Kelm, Cyrille Cuenin, Geoffroy Durand, Isabelle Romieu, Zdenko Herceg. Epigenome-wide association study to predict breast cancer risk in European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 386. doi:10.1158/1538-7445.AM2014-386

Background: The quality of tissue samples and extracted mRNA is a major source of variability in tumor transcriptome analysis using genome-wide expression microarrays. During and immediately after surgical tumor resection, tissues are exposed to metabolic, biochemical and physical stresses characterized as “warm ischemia”. Current practice advocates cryopreservation of biosamples within 30 minutes of resection, but this recommendation has not been systematically validated by measurements of mRNA decay over time. Methods: Using Illumina HumanHT-12 v3 Expression BeadChips, providing a genome-wide coverage of over 24,000 genes, we have analysed gene expression variation in samples of 3 hepatocellular carcinomas (HCC) and 3 lung carcinomas (LC) cryopreserved at times up to 2 hours after resection. Results: RNA Integrity Numbers (RIN) revealed no significant deterioration of mRNA up to 2 hours after resection. Genome-wide transcriptome analysis detected non-significant gene expression variations of -2.4%/hr (95% CI: -5.5%/hr to 0.6%/hr; p=0.12). In LC, no consistent gene expression pattern was detected in relation with warm ischemia. In HCC, a signature of 6 up-regulated genes (CYP2E1, IGLL1, CABYR, CLDN2, NQO1, SCL13A5) and 6 down-regulated (MT1G, MT1H, MT1E, MT1F, HABP2, SPINK1) was identified (FDR Conclusions: Overall, our observations support current recommendation of time to cryopreservation of up to 30 minutes and emphasize the need for identifying tissue-specific genes deregulated following resection to avoid misinterpreting expression changes induced by warm ischemia as pathologically significant changes. Citation Format: Florence Le Calvez-Kelm, Elodie Caboux, Maria Paciencia, Geoffroy Durand, Nivonirina Robinot, Magdalena Wozniak, Francoise Galateau-Salle, Graham Byrnes, Pierre Hainaut. Impact of delay to cryopreservation on RNA integrity and genome-wide expression profiles in resected tumor samples. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2000. doi:10.1158/1538-7445.AM2013-2000

Gene expression microarray and next generation sequencing efforts on conventional, clear cell renal cell carcinoma (ccRCC) have been mostly performed in North American and Western European populations, while the highest incidence rates are found in Central/Eastern Europe. We conducted whole-genome expression profiling on 101 pairs of ccRCC tumours and adjacent non-tumour renal tissue from Czech patients recruited within the “K2 Study”, using the Illumina HumanHT-12 v4 Expression BeadChips to explore the molecular variations underlying the biological and clinical heterogeneity of this cancer. Differential expression analysis identified 1650 significant probes (fold change ≥2 and false discovery rate

Renal cell carcinoma (RCC) is a poorly documented cancer, although time trends and geographical variations in incidence are intriguing. In addition to the frequent VHL mutation, recent studies have identified novel somatic alterations involved in RCC pathogenesis. Microarray and next generation sequencing efforts have been mostly performed on North American and Western European populations, while the highest incidence rates are found in Central/Eastern Europe. Our project has focused on the delineation of gene expression patterns specific to the Czech Republic population, where the incidence rates are the highest in the world. As a feasibility study, we performed whole genome expression profiling on 12 pairs of fresh frozen conventional clear cell RCC (ccRCC) tumors and adjacent non-tumor renal tissue using the Illumina HumanHT-12 v4 Expression BeadChips. To identify genes differentially expressed in ccRCC tumor versus normal tissue samples, the Significance Analysis of Microarrays (SAM) was used. Gene set enrichment analysis (GSEA) software was used for analysis of differentially expressed BioCarta categories and pathways associated with ccRCC pathogenesis. Following initial analyses a signature composed of 994 downregulated and 872 upregulated genes was found in ccRCC compared with normal tissue (FDR Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 5103. doi:1538-7445.AM2012-5103

Introduction Both protein-truncating variants and some missense substitutions in CHEK2 confer increased risk of breast cancer. However, no large-scale study has used full open reading frame mutation screening to assess the contribution of rare missense substitutions in CHEK2 to breast cancer risk. This absence has been due in part to a lack of validated statistical methods for summarizing risk attributable to large numbers of individually rare missense substitutions. Methods Previously, we adapted an in silico assessment of missense substitutions used for analysis of unclassified missense substitutions in BRCA1 and BRCA2 to the problem of assessing candidate genes using rare missense substitution data observed in case-control mutation-screening studies. The method involves stratifying rare missense substitutions observed in cases and/or controls into a series of grades ordered a priori from least to most likely to be evolutionarily deleterious, followed by a logistic regression test for trends to compare the frequency distributions of the graded missense substitutions in cases versus controls. Here we used this approach to analyze CHEK2 mutation-screening data from a population-based series of 1,303 female breast cancer patients and 1,109 unaffected female controls. Results We found evidence of risk associated with rare, evolutionarily unlikely CHEK2 missense substitutions. Additional findings were that (1) the risk estimate for the most severe grade of CHEK2 missense substitutions (denoted C65) is approximately equivalent to that of CHEK2 protein-truncating variants; (2) the population attributable fraction and the familial relative risk explained by the pool of rare missense substitutions were similar to those explained by the pool of protein-truncating variants; and (3) post hoc power calculations implied that scaling up case-control mutation screening to examine entire biochemical pathways would require roughly 2,000 cases and controls to achieve acceptable statistical power. Conclusions This study shows that CHEK2 harbors many rare sequence variants that confer increased risk of breast cancer and that a substantial proportion of these are missense substitutions. The study validates our analytic approach to rare missense substitutions and provides a method to combine data from protein-truncating variants and rare missense substitutions into a one degree of freedom per gene test.

The susceptibility gene for ataxia telangiectasia, ATM, is also an intermediate-risk breast-cancer-susceptibility gene. However, the spectrum and frequency distribution of ATM mutations that confer increased risk of breast cancer have been controversial. To assess the contribution of rare variants in this gene to risk of breast cancer, we pooled data from seven published ATM case-control mutation-screening studies, including a total of 1544 breast cancer cases and 1224 controls, with data from our own mutation screening of an additional 987 breast cancer cases and 1021 controls. Using an in silico missense-substitution analysis that provides a ranking of missense substitutions from evolutionarily most likely to least likely, we carried out analyses of protein-truncating variants, splice-junction variants, and rare missense variants. We found marginal evidence that the combination of ATM protein-truncating and splice-junction variants contribute to breast cancer risk. There was stronger evidence that a subset of rare, evolutionarily unlikely missense substitutions confer increased risk. On the basis of subset analyses, we hypothesize that rare missense substitutions falling in and around the FAT, kinase, and FATC domains of the protein may be disproportionately responsible for that risk and that a subset of these may confer higher risk than do protein-truncating variants. We conclude that a comparison between the graded distributions of missense substitutions in cases versus controls can complement analyses of truncating variants and help identify susceptibility genes and that this approach will aid interpretation of the data emerging from new sequencing technologies.