Search

Esophageal squamous cell carcinoma (ESCC) shows remarkable variation in incidence that is not fully explained by known lifestyle and environmental risk factors. It has been speculated that an unknown exogenous exposure(s) could be responsible. Here we combine the fields of mutational signature analysis with cancer epidemiology to study 552 ESCC genomes from eight countries with varying incidence rates. Mutational profiles were similar across all countries studied. Associations between specific mutational signatures and ESCC risk factors were identified for tobacco, alcohol, opium and germline variants, with modest impacts on mutation burden. We find no evidence of a mutational signature indicative of an exogenous exposure capable of explaining differences in ESCC incidence. Apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like (APOBEC)-associated mutational signatures single-base substitution (SBS)2 and SBS13 were present in 88% and 91% of cases, respectively, and accounted for 25% of the mutation burden on average, indicating that APOBEC activation is a crucial step in ESCC tumor development.

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

While several high-penetrance melanoma risk genes are known, variation in these genes fail to explain melanoma susceptibility in a large proportion of high-risk families. As part of a melanoma family sequencing study, including 435 families from Mediterranean populations, we identified a novel NRAS variant (c.170A>C, p.D57A) in a melanoma-prone family. This variant is absent in exomes in gnomAD, ESP, UKBiobank, and the 1000 Genomes Project, as well as in 11 273 Mediterranean individuals and 109 melanoma-prone families from the US and Australia. This variant occurs in the GTP-binding pocket of NRAS. Differently from other RAS activating alterations, NRAS D57A expression is unable to activate MAPK-pathway both constitutively and after stimulation but enhances EGF-induced PI3K-pathway signaling in serum starved conditions in vitro. Consistent with in vitro data demonstrating that NRAS D57A does not enrich GTP binding, molecular modeling suggests that the D57A substitution would be expected to impair Mg2+ binding and decrease nucleotide-binding and GTPase activity of NRAS. While we cannot firmly establish NRAS c.170A>C (p.D57A) as a melanoma susceptibility variant, further investigation of NRAS as a familial melanoma gene is warranted.

Background & Aims Gallbladder cancer (GBC) is the most common type of biliary tract cancer, but the molecular mechanisms involved in gallbladder carcinogenesis remain poorly understood. In this study, we applied integrative genomics approaches to characterise GBC and explore molecular subtypes associated with patient survival. Methods We profiled the mutational landscape of GBC tumours (whole-exome sequencing on 92, targeted sequencing on 98, in total 190 patients). In a subset (n = 45), we interrogated the matched transcriptomes, DNA methylomes, and somatic copy number alterations. We explored molecular subtypes identified through clustering tumours by genes whose expression was associated with survival in 47 tumours and validated subtypes on 34 publicly available GBC cases. Results Exome analysis revealed TP53 was the most mutated gene. The overall mutation rate was low (median 0.82 Mut/Mb). APOBEC-mediated mutational signatures were more common in tumours with higher mutational burden. Aflatoxin-related signatures tended to be highly clonal (present in ≥50% of cancer cells). Transcriptome-wide survival association analysis revealed a 95-gene signature that stratified all GBC patients into 3 subtypes that suggested an association with overall survival post-resection. The 2 poor-survival subtypes were associated with adverse clinicopathologic features (advanced stage, pN1, pM1), immunosuppressive micro-environments (myeloid-derived suppressor cell accumulation, extensive desmoplasia, hypoxia) and T cell dysfunction, whereas the good-survival subtype showed the opposite features. Conclusion These data suggest that the tumour micro-environment and immune profiles could play an important role in gallbladder carcinogenesis and should be evaluated in future clinical studies, along with mutational profiles. Lay summary Gallbladder cancer is highly fatal, and its causes are poorly understood. We evaluated gallbladder tumours to see if there were differences between tumours in genetic information such as DNA and RNA. We found evidence of aflatoxin exposure in these tumours, and immune cells surrounding the tumours were associated with survival.

Background Blood type has been associated with the risk of gastric cancer, but few studies have examined the association with esophageal squamous cell carcinoma (ESCC). Methods We conducted a case-control study using genotyping data of Chinese individuals, including cases of 2022 ESCC, 1189 gastric cardia adenocarcinoma, 1161 gastric noncardia adenocarcinoma, and 2696 controls. Genetic blood type was imputed using three single nucleotide polymorphisms. We used logistic regression to examine the association between blood type and the risk of each cancer. Results Compared to blood type O, the risk of ESCC was significantly elevated for blood type B and AB, with the highest risk for type AB (OR, 95%CI: 1.34, 1.07–1.67). Analysis of genotype suggested that the association of ESCC was from carrying the B allele. Similarly, blood type was significantly associated with gastric noncardia adenocarcinoma (P P = 0.13). Conclusions This study provides novel insights into the association between blood type and the risk of ESCC and restricted previously observed association to only gastric noncardia cancer, providing important evidence to clarify the pattern of association and suggesting mechanisms of action.

Background: Skin cancer screening is routinely performed for members of melanoma-prone families, but longitudinal studies evaluating the efficacy of surveillance in this high-risk population are lacking. Methods: We evaluated thickness for first primary melanomas diagnosed in melanoma-prone families (≥2 individuals with melanoma) enrolled in NCT00040352 (NCI familial melanoma study) from 1976 through 2014; enrolled patients received routine skin cancer screening and education about skin self-exams. We used linear and ordinal logistic regression models adjusted for gender and age with a generalized estimating equations approach to report changes in thickness and tumor (T) stage over time, comparing outcomes for NCI cases diagnosed before (pre-study) versus after study participation (prospective) and for NCI cases versus nonfamilial cases [Surveillance, Epidemiology, and End Results (SEER) 9 registries]. Results: Tumor thickness was evaluated for 293 NCI (pre-study = 246; prospective = 47) patients. Compared with NCI pre-study cases, NCI prospective melanomas were thinner (0.6 vs. 1.1 mm; P < 0.001) and more likely to be T1 stage [39/47 (83%) vs. 98/246 (40%); P < 0.001]. Similar findings (P < 0.05) were observed for familial cases with and without germline CDKN2A and CDK4 mutations. Peters–Belson modeling suggested that calendar period effects of decreasing thickness in the general population (SEER 9) did not fully explain thickness trends in NCI families. Conclusions: Participation in a longitudinal surveillance program providing skin cancer screening and education about skin self-exams was associated with thinner melanomas for members of melanoma-prone families. Impact: The study findings support the clinical benefit of screening (physician and self) for this high-risk population.

BACKGROUND: CDKN2A, CDK4, and POT1 are well-established melanoma-susceptibility genes. OBJECTIVE: We evaluated melanoma histopathology for individuals with germline mutations of CDKN2A, CDK4, and POT1. METHODS: We assessed histopathology for melanomas diagnosed in melanoma-prone families (≥2 individuals with melanoma) from the United States, Italy, and Spain. Comparisons between mutation carriers and non-carriers (no mutation) were adjusted for age, sex, Breslow depth, and correlations among individuals within the same family. RESULTS: Histologic slides were evaluated for 290 melanomas (139 from 132 non-carriers, 122 from 68 CDKN2A carriers, 10 from 6 CDK4 carriers, 19 from 16 POT1 carriers). Superficial spreading was the predominant subtype for all groups. Spitzoid morphology (>25% of tumor) was observed in the majority of invasive melanomas from POT1 carriers (10 of 15 [67%], P

Background Height and body mass index (BMI) have both been positively associated with melanoma risk, although findings for BMI have been less consistent than height. It remains unclear, however, whether these associations reflect causality or are due to residual confounding by environmental and lifestyle risk factors. We re-evaluated these associations using a two-sample Mendelian randomization (MR) approach. Methods We identified single nucleotide polymorphisms (SNPs) for BMI and height from separate genome-wide association study (GWAS) meta-analyses. We obtained melanoma SNPs from the most recent melanoma GWAS meta-analysis comprising 12 874 cases and 23 203 controls. We used the inverse variance-weighted estimator to derive separate causal risk estimates across all SNP instruments for BMI and height. Results Based on the combined estimate derived from 730 SNPs for BMI, we found no evidence of an association between genetically predicted BMI and melanoma [odds ratio (OR) per one standard deviation (1 SD) (4.6 kg/m2) increase in BMI 1.00, 95% confidence interval (CI): 0.91–1.11]. In contrast, we observed a positive association between genetically-predicted height (derived from a pooled estimate of 3290 SNPs) and melanoma risk [OR 1.08, 95% CI: 1.02–1.13, per 1 SD (9.27 cm) increase in height]. Sensitivity analyses using two alternative MR methods yielded similar results. Conclusions These findings provide no evidence for a causal association between higher BMI and melanoma, but support the notion that height is causally associated with melanoma risk. Mechanisms through which height influences melanoma risk remain unclear, and it remains possible that the effect could be mediated through diverse pathways including growth factors and even socioeconomic status.

The incidence of esophageal squamous cell carcinoma (ESCC) is disproportionately high in the eastern corridor of Africa and parts of Asia. Emerging research has identified a potential association between poor oral health and ESCC. One proposed biological pathway linking poor oral health and ESCC involves the alteration of the microbiome. Thus, we performed an integrated analysis of four independent sequencing efforts of ESCC tumors from patients from high- and low-incidence regions of the world. Using whole genome sequencing (WGS) and RNA sequencing (RNAseq) of ESCC tumors and WGS of synchronous collections of saliva specimens from 61 patients in Tanzania, we identified a community of bacteria, including members of the genera Fusobacterium, Selenomonas, Prevotella, Streptococcus, Porphyromonas, Veillonella, and Campylobacter, present at high abundance in ESCC tumors. We then characterized the microbiome of 238 ESCC tumor specimens collected in two additional independent sequencing efforts consisting of patients from other high-ESCC incidence regions (Tanzania, Malawi, Kenya, Iran, China). This analysis revealed a similar tumor enrichment of the ESCC-associated bacterial community in these cancers. Because these genera are traditionally considered members of the oral microbiota, we explored if there is a relationship between the synchronous saliva and tumor microbiomes of ESCC patients in Tanzania. Comparative analyses revealed that paired saliva and tumor microbiomes are significantly similar with a specific enrichment of Fusobacterium and Prevotella in the tumor microbiome. Together, these data indicate that cancer-associated oral bacteria are associated with ESCC tumors at the time of diagnosis and support a model in which oral bacteria are present in high abundance in both saliva and tumors of ESCC patients. Longitudinal studies of the pre-diagnostic oral microbiome are needed to investigate whether these cross-sectional similarities reflect temporal associations.

Background Genetic and environmental factors are important determinants of nasopharyngeal carcinoma (NPC). NPC is associated with Epstein-Barr virus (EBV) infection. Studies have reported familial aggregation of NPC, but evidence has been mixed for elevated rates of cancers other than NPC. Methods The authors reassessed their previous evaluation of familial aggregation of cancer in 348 high-risk Taiwanese multiplex families with 2 or more NPC cases enrolled between 1980 and 2003. Participants were linked to the Taiwan National Cancer Registry and National Death Registry to identify cancers. Results In all, 2590 individuals contributed 37,959 person-years over an average of 15 years of follow-up; 314 incident cancers were identified. The authors computed multiple primary standardized incidence ratios (MP-SIRs) to evaluate the overall risk and the risk of infection-associated, EBV-associated, and individual cancers. The overall MP-SIR was 1.24 (95% confidence interval [CI], 1.10-1.38). The exclusion of excess NPC risk led to an overall MP-SIR of 1.11 (95% CI, 0.98-1.25). Similarly, the risk of cancers associated with infectious agents was driven by the excess in NPC, and its exclusion led to an MP-SIR of 1.22 (95% CI, 0.99-1.48) for infection-associated cancers and to an MP-SIR of 1.18 (95% CI, 0.72-1.82) for EBV-associated cancers. The authors observed a significant excess of second cancers among NPC cases (oral cancer, mouth cancer, tongue cancer, gum cancer, nasal cavity cancer, bone cancer, and non-Hodgkin lymphoma). Conclusions This reassessment of the largest NPC multiplex family study confirms the presence of NPC coaggregation within families in Taiwan but does not provide evidence for a broader familial syndrome involving NPC and other tumors. Among NPC cases, elevated rates of secondary cancers, mostly at the, head and neck and hematopoietic cancers suggest radiation treatment effects on subsequent cancer risk.

To assess the age at menarche and menopause of women with xeroderma pigmentosum, a DNA repair disease with premature aging, in a longitudinal natural history study.We conducted a natural history study that reviewed medical records for gynecologic and reproductive health of all female patients with xeroderma pigmentosum aged older than 9 years examined at the National Institutes of Health (NIH). We performed gynecologic and laboratory examinations on a subset of the patients. Women in a second subset, who could not be examined, were interviewed using a questionnaire. Women who were deceased or lost to follow-up formed a third subset.Sixty females with xeroderma pigmentosum aged older than 9 years (median 29 years, range 10-61 years) were evaluated at the NIH from 1971 to 2018. Of these 60, 31 had history, questionnaire, record review, and gynecologic evaluation; 14 had record review and questionnaire interview by telephone; and 15 had only NIH record review. Menarche in females with xeroderma pigmentosum occurred at a median age of 12.0 years (range 9-17 years), which was comparable with the U.S. general population. Among the 18 patients with menopause, the median age at menopause of 29.5 years (range 18-49.5 years) was more than 20 years younger than in the U.S. general population (52.9 years). Premature menopause (before age 40 years) occurred in 14 of the 45 (31%) women aged 18 years or older, and primary ovarian insufficiency was documented in nine of them. There were 32 live births among 21 of the women, five of whom subsequently developed premature menopause.Females with xeroderma pigmentosum in our study had a normal age at menarche and were fertile but had increased incidence of premature menopause. Premature menopause, a symptom of premature aging, should be considered for gynecologic and reproductive health as well as implicating DNA repair in maintaining normal ovarian function.ClinicalTrials.gov, NCT00001813.

Background: Genetic susceptibility is associated with nasopharyngeal carcinoma (NPC). We previously identified rare variants potentially involved in familial NPC and common variants significantly associated with sporadic NPC. Methods: We conducted targeted gene sequencing of 20 genes [16 identified from the study of multiplex families, three identified from a pooled analysis of NPC genome-wide association study (GWAS), and one identified from both studies] among 819 NPC cases and 938 controls from two case–control studies in Taiwan (independent from previous studies). A targeted, multiplex PCR primer panel was designed using the custom Ion AmpliSeq Designer v4.2 targeting the regions of the selected genes. Gene-based and single-variant tests were conducted. Results: We found that NPC was associated with combined common and rare variants in CDKN2A/2B (P = 1.3 × 10−4), BRD2 (P = 1.6 × 10−3), TNFRSF19 (P = 4.0 × 10−3), and CLPTM1L/TERT (P = 5.4 × 10−3). Such associations were likely driven by common variants within these genes, based on gene-based analyses evaluating common variants and rare variants separately (e.g., for common variants of CDKN2A/2B, P = 4.6 × 10−4; for rare variants, P = 0.04). We also observed a suggestive association with rare variants in HNRNPU (P = 3.8 × 10−3) for NPC risk. In addition, we validated four previously reported NPC risk–associated SNPs. Conclusions: Our findings confirm previously reported associated variants and suggest that some common variants in genes previously linked to familial NPC are associated with the development of sporadic NPC. Impact: NPC-associated genes, including CLPTM1L/TERT, BRD2, and HNRNPU, suggest a role for telomere length maintenance in NPC etiology.

Background People at high risk of developing melanoma are usually identified by pigmentary and naevus phenotypes. Objective We examined whether associations of these phenotypes with melanoma risk differed by ambient sun exposure or participant characteristics in two population‐based, case–control studies with comparable ancestry but different ambient sun exposure. Methods Data were analysed from 616 cases and 496 controls from the Australian Melanoma Family Study and 2012 cases and 504 controls from the Leeds (UK) case–control study. Questionnaire, interview and dermatological skin examination data were collected using the same measurement protocols. Relative risks were estimated as odds ratios using unconditional logistic regression, adjusted for potential confounders. Results Hair and skin colour were the strongest pigmentary phenotype risk factors. All associations of pigmentary phenotype with melanoma risk were similar across countries. The median number of clinically assessed naevi was approximately three times higher in Australia than Leeds, but the relative risks for melanoma associated with each additional common or dysplastic naevus were higher for Leeds than Australia, especially for naevi on the upper and lower limbs. Higher naevus counts on the head and neck were associated with a stronger relative risk for melanoma for women than men. The two countries had similar relative risks for melanoma based on self‐reported naevus density categories, but personal perceptions of naevus number differed by country. There was no consistent evidence of interactions between phenotypes on risk. Conclusions Classifying people at high risk of melanoma based on their number of naevi should ideally take into account their country of residence, type of counts (clinical or self‐reported), body site on which the naevus counts are measured and sex. The presence of naevi may be a stronger indicator of a genetic predisposition in the UK than in Australia based on less opportunity for sun exposure to influence naevus development.

Purpose More than half of the familial cutaneous melanomas have unknown genetic predisposition. This study aims at characterizing a novel melanoma susceptibility gene. Methods We performed exome and targeted sequencing in melanoma-prone families without any known melanoma susceptibility genes. We analyzed the expression of candidate gene DENND5A in melanoma samples in relation to pigmentation and UV signature. Functional studies were carried out using microscopic approaches and zebrafish model. Results We identified a novel DENND5A truncating variant that segregated with melanoma in a Swedish family and 2 additional rare DENND5A variants, 1 of which segregated with the disease in an American family. We found that DENND5A is significantly enriched in pigmented melanoma tissue. Our functional studies show that loss of DENND5A function leads to decrease in melanin content in vitro and pigmentation defects in vivo. Mechanistically, harboring the truncating variant or being suppressed leads to DENND5A losing its interaction with SNX1 and its ability to transport the SNX1-associated vesicles from melanosomes. Consequently, untethered SNX1-premelanosome protein and redundant tyrosinase are redirected to lysosomal degradation by default, causing decrease in melanin content. Conclusion Our findings provide evidence of a physiological role of DENND5A in the skin context and link its variants to melanoma susceptibility.

Esophageal squamous cell carcinoma (ESCC) shows a remarkable variation in incidence which is not fully explained by known lifestyle and environmental risk factors. It has been speculated that an unknown exogenous exposure(s) could be responsible. Here we combine the fields of mutational signature analysis with cancer epidemiology to study 552 ESCC genomes from eight countries with varying incidence rates. The mutational profiles of ESCC were similar across all countries studied. Associations between specific mutational signatures and ESCC risk factors were identified for tobacco, alcohol, opium and germline variants, with modest impacts on mutation burden. We find no evidence of a mutational signature indicative of an exogenous exposure capable of explaining the differences in ESCC incidence. APOBEC associated mutational signatures SBS2 and SBS13 were present in 88% and 91% of cases respectively and accounted for a quarter of the mutation burden on average, indicating that activation of APOBEC is a crucial step in ESCC tumor development.

Background: Chordoma is a rare bone cancer with an unknown etiology. TBXT is the only chordoma susceptibility gene identified to date, germline single nucleotide variants and copy number variants in TBXT have been associated with chordoma susceptibility in familial and sporadic chordoma. However, the genetic susceptibility of chordoma remains largely unknown. In this study, we investigated rare germline genetic variants in genes involved in TBXT/chordoma-related signaling pathways and other biological processes in chordoma patients from North America and China. Methods: We identified variants that were very rare in general population and internal control datasets and showed evidence for pathogenicity in 265 genes in a whole exome sequencing (WES) dataset of 138 chordoma patients of European ancestry and in a whole genome sequencing (WGS) dataset of 80 Chinese patients with skull base chordoma. Results: Rare and likely pathogenic variants were identified in 32 of 138 European ancestry patients (23%), including genes that are part of notochord development, PI3K/AKT/mTOR, Sonic Hedgehog, SWI/SNF complex and mesoderm development pathways. Rare pathogenic variants in COL2A1, EXT1, PDK1, LRP2, TBXT and TSC2, among others, were also observed in Chinese patients. Conclusion: We identified several rare loss-of-function and predicted deleterious missense variants in germline DNA from patients with chordoma, which may influence chordoma predisposition and reflect a complex susceptibility, warranting further investigation in large studies.

The genetics of pancreatic ductal adenocarcinoma (PDAC) is complex with patients reported to harbor germline pathogenic variants (PVs) in many different genes. PDAC patients with familial pancreatic cancer (FPC) are more likely to carry germline PVs but there is no consensus main gene involved in FPC. We performed a systematic review of publications from PubMed and Scopus reporting PVs in patients with FPC, sporadic pancreatic cancer (SPC) and unselected cohorts of PDAC patients undergoing genetic testing and calculated a cumulative prevalence of PVs for each gene evaluated across these three groups of patients. When available, variants in the selected publications were reclassified according to the American College of Medical Genetics and Genomics classification system and used for prevalence calculations if classified as pathogenic or likely pathogenic. We observed an increased prevalence of PVs in FPC compared to SPC or unselected PDAC patients for most of the 41 genes reported. The genes with the highest prevalence of carriers of PVs in FPC were ATM, BRCA2, and CDKN2A. BRCA2 and ATM showed the highest prevalence of PVs in both SPC and unselected PDAC cohorts. Several genes with the highest prevalence of PVs are involved in breast and ovarian cancer suggesting strong overlap with underlying genetics in these disorders but no single gene was predominant. More research is needed to further understand the risk of PDAC associated with these many diverse genes.

Genome-wide association studies have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP, and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). As ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWAS identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor gene (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. As AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus, and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression, and altering melanocyte growth phenotypes upon exposure.

While expression quantitative trait loci (eQTL) have been powerful in identifying susceptibility genes from genome-wide association studies (GWAS) findings, most trait-associated loci are not explained by eQTL alone. Alternative QTLs including DNA methylation QTL (meQTL) are emerging, but cell-type-specific meQTL using cells of disease origin has been lacking. Here we established an meQTL dataset using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization using meQTL, eQTL, and mRNA splice-junction QTL from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified susceptibility genes at 63% of melanoma GWAS loci. Among three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes are preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTL and melanoma meQTL identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTL identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIPseq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTL.

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

Background Pediatric cancers are the leading cause of death by disease in children despite improved survival rates overall. The contribution of germline genetic susceptibility to pediatric cancer survivors has not been extensively characterized. We assessed the frequency of pathogenic or likely pathogenic (P/LP) variants in 5451 long-term pediatric cancer survivors from the Childhood Cancer Survivor Study. Methods Exome sequencing was conducted on germline DNA from 5451 pediatric cancer survivors (cases who survived ≥5 years from diagnosis; n = 5105 European) and 597 European cancer-free adults (controls). Analyses focused on comparing the frequency of rare P/LP variants in 237 cancer-susceptibility genes and a subset of 60 autosomal dominant high-to-moderate penetrance genes, for both case-case and case-control comparisons. Results Of European cases, 4.1% harbored a P/LP variant in high-to-moderate penetrance autosomal dominant genes compared with 1.3% in controls (2-sided P = 3 × 10-4). The highest frequency of P/LP variants was in genes typically associated with adult onset rather than pediatric cancers, including BRCA1/2, FH, PALB2, PMS2, and CDKN2A. A statistically significant excess of P/LP variants, after correction for multiple tests, was detected in patients with central nervous system cancers (NF1, SUFU, TSC1, PTCH2), Wilms tumor (WT1, REST), non-Hodgkin lymphoma (PMS2), and soft tissue sarcomas (SDHB, DICER1, TP53, ERCC4, FGFR3) compared with other pediatric cancers. Conclusion In long-term pediatric cancer survivors, we identified P/LP variants in cancer-susceptibility genes not previously associated with pediatric cancer as well as confirmed known associations. Further characterization of variants in these genes in pediatric cancer will be important to provide optimal genetic counseling for patients and their families.

Background: Melanoma aetiology has been proposed to have two pathways, which are determined by naevi and type of sun exposure and related to the anatomical site where melanoma develops. Objectives: We examined associations with melanoma by anatomical site for a comprehensive set of risk factors including pigmentary and naevus phenotypes, ultraviolet radiation exposure and polygenic risk. Methods: We analysed harmonized data from 2617 people with incident first invasive melanoma and 975 healthy controls recruited through two population-based case–control studies in Australia and the UK. Questionnaire data were collected by interview using a single protocol, and pathway-specific polygenic risk scores were derived from DNA samples. We estimated adjusted odds ratios using unconditional logistic regression that compared melanoma cases at each anatomical site with all controls. Results: When cases were compared with control participants, there were stronger associations for many naevi vs. no naevi for melanomas on the trunk, and upper and lower limbs than on the head and neck (P-heterogeneity < 0·001). Very fair skin (vs. olive/brown skin) was more weakly related to melanoma on the trunk than to melanomas at other sites (P-heterogeneity = 0·04). There was no significant difference by anatomical site for polygenic risk. Increased weekday sun exposure was positively associated with melanoma on the head and neck but not on other sites. Conclusions: We found evidence of aetiological heterogeneity for melanoma, supporting the dual pathway hypothesis. These findings enhance understanding of risk factors for melanoma and can guide prevention and skin examination education and practices.

Importance Risk of sebaceous carcinoma (SC), a rare skin cancer associated with Muir-Torre syndrome, is elevated among solid organ transplant recipients (SOTRs). However, population studies evaluating this association and assessing survival for posttransplant cases are lacking, and further understanding of SC epidemiology in this immunosuppressed population could provide etiologic and clinical insights. Objective To assess SC incidence and patient survival after solid organ transplantation. Design, Setting, and Participants This cohort study, conducted from January 1, 1987, to December 31, 2017, used data from the Transplant Cancer Match Study, which links transplant and cancer registry data for 17 states and 1 metropolitan area in the United States. Altogether, these registries account for approximately 46% of all US transplants. Data on demographic and transplant characteristics as well as induction and initial maintenance immunosuppressive therapies were obtained from the transplant registry. Standardized incidence ratios (SIRs) comparing SC incidence among SOTRs to the general population were calculated. Incidence rate ratios (IRRs) comparing SC risk between SOTR subgroups were calculated using multivariate Poisson regression. Cox regression was used to compare overall survival between SC cases in SOTRs and other individuals. Main Outcomes and Measures Sebaceous carcinoma incidence and overall patient survival after transplantation compared with the general population. Results A total of 326 282 transplant procedures were performed for 301 075 patients (No. [%] age at transplant, 126 550 [38.8%] aged 0-44 years; 82 394 [25.3%] aged 45-54 years; 82 082 [25.5%] aged 55-64 years; 35 256 [10.8%] aged ≥65 years; 201 354 male patients [61.7%]; 202 557 White patients [62.1%]). A total of 102 SCs were diagnosed in 301 075 SOTRs, corresponding to a 25-fold increased incidence (SIR, 24.8; 95% CI, 20.2-30.1). Incidence was especially elevated among lung recipients (SIR, 47.7; 95% CI, 20.6-94.0) and after a posttransplant diagnosis of cutaneous squamous cell carcinoma (SIR, 104.0; 95% CI, 62.8-163.0). Among SOTRs, factors independently associated with SC risk included male sex (IRR, 2.46; 95% CI, 1.48-4.07; P

Sebaceous carcinoma is an aggressive skin cancer with a 5-year overall survival rate of 78% for localized/regional disease and 50% for metastatic disease. The incidence of this cancer has been increasing in the United States for several decades, but the underlying reasons for this increase are unclear. In this article, we review the epidemiology and genetics of sebaceous carcinoma, including recent population data and tumor genomic analyses that provide new insights into underlying tumor biology. We further discuss emerging evidence of a possible viral etiology for this cancer. Finally, we review the clinical implications of recent advances in sebaceous carcinoma research for screening, prevention, and treatment.

Background Esophageal squamous cell carcinoma (ESCC) is a leading cause of cancer death worldwide and in China. We know miRNAs influence gene expression in tumorigenesis, but it is unclear how miRNAs affect gene expression or influence survival at the genome-wide level in ESCC. Methods We performed miRNA and mRNA expression arrays in 113 ESCC cases with tumor/normal matched tissues to identify dysregulated miRNAs, to correlate miRNA and mRNA expressions, and to relate miRNA and mRNA expression changes to survival and clinical characteristics. Results Thirty-nine miRNAs were identified whose tumor/normal tissue expression ratios showed dysregulation (28 down- and 11 up-regulated by at least two-fold with P P P NF1, ASXL1, HSPA4, TGOLN2, BAIAP2, EZH2, CHAF1A, SUPT7L. Conclusions Our characterization and integrated analysis of genome-wide miRNA and gene expression in ESCC provides insights into the expression of miRNAs and their relation to regulation of RNA targets in ESCC tumorigenesis, and suggest opportunities for the future development of miRs and mRNAs as biomarkers for early detection, diagnosis, and prognosis in ESCC.

BACKGROUND: Epstein-Barr virus (EBV) positivity is associated with better gastric cancer prognosis and is found in a relatively fixed 9% of tumors worldwide. AIMS: We aimed to examine the EBV status of gastric adenocarcinomas in a very high-incidence population and to compare prevalence between cardia and non-cardia anatomic subsites. METHODS: We evaluated 1035 adult gastric adenocarcinoma cases presenting during 1997-2005 to the Shanxi Cancer Hospital in Taiyuan, Shanxi Province, China. EBV-encoded RNA was detected in alcohol-fixed paraffin-embedded tumor specimens by in situ hybridization. Associations were assessed in case-case comparisons using the chi-squared test for categorical variables and the Mann-Whitney U Test for continuous variables, with p-values < 0.05 considered statistically significant. Adjusted odds ratios were calculated using logistic regression, and mortality hazard ratios (HRs) were estimated by Cox proportional hazards regression. RESULTS: Sixty-four percent of the evaluated cancers were found in the cardia. Cardia tumor localization was associated with male sex, advanced primary tumor stage, better differentiated histology, and intestinal-type Lauren classification. Four percent of the non-cardia and only 0.9% of cardia cancers were EBV-positive. EBV-positivity was associated with better overall survival (adjusted HR 0.30, 95% CI 0.14-0.63). CONCLUSIONS: Our study highlights unusually low EBV prevalence in gastric adenocarcinoma among a high-incidence population, particularly for cardia cancers. These findings suggest a unique risk factor profile for the high incidence of gastric cancer in this population.

Chordoma is a rare bone tumor with an unknown etiology and high recurrence rate. Here we conduct whole genome sequencing of 80 skull-base chordomas and identify PBRM1, a SWI/SNF (SWItch/Sucrose Non-Fermentable) complex subunit gene, as a significantly mutated driver gene. Genomic alterations in PBRM1 (12.5%) and homozygous deletions of the CDKN2A/2B locus are the most prevalent events. The combination of PBRM1 alterations and the chromosome 22q deletion, which involves another SWI/SNF gene (SMARCB1), shows strong associations with poor chordoma-specific survival (Hazard ratio [HR] = 10.55, 95% confidence interval [CI] = 2.81-39.64, p = 0.001) and recurrence-free survival (HR = 4.30, 95% CI = 2.34-7.91, p = 2.77 × 10−6). Despite the low mutation rate, extensive somatic copy number alterations frequently occur, most of which are clonal and showed highly concordant profiles between paired primary and recurrence/metastasis samples, indicating their importance in chordoma initiation. In this work, our findings provide important biological and clinical insights into skull-base chordoma., Skull base chordomas are treated with surgery and chemotherapy but often recur due to incomplete resection, understanding the molecular underpinnings of the tumours may provide additional therapeutic strategies. Here, the authors carry out whole genome sequencing of 80 skull base chordoma tumours and identify the SWI/SNF component—PBRM1—as a frequently mutated gene.

Sebaceous carcinoma (SC) is an aggressive skin tumor. Although ultraviolet radiation (UVR) is an important risk factor for some skin cancer types, no population-level study has evaluated for an association between UVR and SC risk. Herein, we examined satellite-based ambient UVR in relation to SC incidence using Surveillance, Epidemiology, and End Results 18 cancer registry data (2000–2016). There were 3503 microscopically confirmed cases of SC diagnosed during the study period. For non-Hispanic whites, there was an association between increasing ambient UVR and SC risk (incidence rate ratio [per UVR quartile] = 1.15, 95% confidence interval = 1.11 to 1.19; two-sided P

Meta-analysis of 36,760 cases and 375,188 controls identifies 54 loci associated with susceptibility to cutaneous melanoma. Further analysis combining nevus count and hair color GWAS results provide insights into the genetic architecture of melanoma.Most genetic susceptibility to cutaneous melanoma remains to be discovered. Meta-analysis genome-wide association study (GWAS) of 36,760 cases of melanoma (67% newly genotyped) and 375,188 controls identified 54 significant (P < 5 x 10(-8)) loci with 68 independent single nucleotide polymorphisms. Analysis of risk estimates across geographical regions and host factors suggests the acral melanoma subtype is uniquely unrelated to pigmentation. Combining this meta-analysis with GWAS of nevus count and hair color, and transcriptome association approaches, uncovered 31 potential secondary loci for a total of 85 cutaneous melanoma susceptibility loci. These findings provide insights into cutaneous melanoma genetic architecture, reinforcing the importance of nevogenesis, pigmentation and telomere maintenance, together with identifying potential new pathways for cutaneous melanoma pathogenesis.

Background Nevoid basal cell carcinoma syndrome (NBCCS) is an autosomal dominant disorder with variable expression and nearly complete penetrance. PTCH1 is the major susceptibility locus and has no known hot spots or genotype–phenotype relationships. Methods We evaluated 18 NBCCS National Cancer Institute (NCI) families plus PTCH1 data on 333 NBCCS disease‐causing mutations (DM) reported in the Human Gene Mutation Database (HGMD). National Cancer Institute families underwent comprehensive genomic evaluation, and clinical data were extracted from NCI and HGMD cases. Genotype–phenotype relationships were analyzed using Fisher's exact tests focusing on mutation type and PTCH1 domains. Results PTCH1 pathogenic mutations were identified in 16 of 18 NCI families, including three previously mutation‐negative families. PTCH1 mutations were spread across the gene with no hot spot. After adjustment for multiple tests, a statistically significant genotype–phenotype association was observed for developmental delay and gross deletion–insertions (p = 9.0 × 10−6), and suggestive associations between falx cerebri calcification and all transmembrane domains (p = 0.002) and severe outcomes and gross deletion–insertions (p = 4.0 × 10−4). Conclusion Overall, 89% of our NCI families had a pathogenic PTCH1 mutation. The identification of PTCH1 mutations in previously mutation‐negative families underscores the importance of repeated testing when new technologies become available. Additional clinical information linked to mutation databases would enhance follow‐up and future studies of genotype–phenotype relationships.

Background In the United States, only approximately 0.4% of all melanomas are diagnosed in patients aged Methods For this non-population-based study, families were followed prospectively for up to 40 years. A total of 60 families with ≥ 3 patients with melanoma were included for analysis: 30 CDKN2A mutation-positive (CDKN2A+) and 30 CDKN2A mutation-negative (CDKN2A-) families. Age at the time of first melanoma and number of melanomas were obtained for each patient and summarized by family or sets (CDKN2A + vs CDKN2A-). For set comparisons and categorical variables (occurrence of melanoma in pediatric patients, number of melanomas, number of patients with single or multiple melanomas), the Pearson chi-square or Fisher exact test was used. Results Regardless of CDKN2A status, melanoma-prone families were found to have 6-fold to 28-fold higher percentages of patients with pediatric melanoma compared with the general population of patients with melanoma in the United States. Within CDKN2A + families, pediatric patients with melanoma were significantly more likely to have multiple melanomas compared with their relatives who were diagnosed at age >20 years (71% vs 38%, respectively; P = .004). CDKN2A + families had significantly higher percentages of pediatric patients with melanoma compared with CDKN2A- families (11.1% vs 2.5%; P = .004). Conclusions These observations have implications for the prevention of melanoma as well as clinical care for its early detection. Children in melanoma-prone families should have careful sun protection from an early age and skin surveillance to reduce their risk of melanoma.

Purpose: Carriers of CDKN2A mutations have high risks of melanoma and certain other cancers. In this study we examined the occurrence of tumors among CDKN2A wild type (wt) members of melanoma-prone families with CDKN2A mutations. Methods: Swedish and US melanoma-prone families with CDKN2A mutations were included. Data was collected on tumors diagnosed among family members. Among the CDKN2A mutated families, members with CDKN2A wt status who were diagnosed with melanoma were designated phenocopies. Results: Of patients with melanoma in the CDKN2A mutated families (n = 266), 7.1%, were seen among members with CDKN2A wt status (phenocopy rate). Among the CDKN2A wt family members of the CDKN2A mutated families (n = 256), 7.4% were diagnosed with melanoma. The prospective relative risk for melanomas was significantly higher among the CDKN2A wt subjects compared with population-based controls (7.4 (95% confidence interval 1.7–33.2)), while no elevated risks of nonmelanoma cancers were seen and their offspring did not have significantly elevated risks of melanoma or other cancers. Conclusion: Members of CDKN2A mutation carrying families who test negative for their family’s mutation have moderately increased risk for melanoma and should, in addition to being considered for continuing dermatologic surveillance, be encouraged to follow sun safety recommendations and practice skin self-exams. (Less)

OBJECTIVE To gain insight into the role of germline genetics in the development of chordoma, the authors evaluated data from 2 sets of patients with familial chordoma, those with and without a germline duplication of the T gene (T-dup+ vs T-dup−), which was previously identified as a susceptibility mechanism in some families. The authors then compared the patients with familial tumors to patients with sporadic chordoma in the US general population reported to the National Cancer Institute’s Surveillance, Epidemiology, and End Results (SEER) program through 2015. METHODS Evaluation of family members included review of personal and family medical history, physical and neurological examination, and pre- and postcontrast MRI of the skull base and spine. Sixteen patients from 6 white families with chordoma had a chordoma diagnosis at family referral. Screening MR images of 35 relatives revealed clival lesions in 6, 4 of which were excised and confirmed to be chordoma. Thus, data were available for 20 patients with histologically confirmed familial chordoma. There were 1759 patients with histologically confirmed chordoma in SEER whose race was known. RESULTS The median age at chordoma diagnosis differed across the groups: it was lowest in T-dup+ familial patients (26.8 years, range 5.3–68.4 years); intermediate in T-dup− patients (46.2 years, range 11.8–60.1 years); and highest in SEER patients (57 years, range 0–98 years). There was a marked preponderance of skull base tumors in patients with familial chordoma (93% in T-dup+ and 83% in T-dup−) versus 38% in the SEER program (37% in white, 53% in black, and 48.5% in Asian/Pacific Islander/American Indian/Alaska Native patients). Furthermore, 29% of white and 16%–17% of nonwhite SEER patients had mobile-spine chordoma, versus no patients in the familial group. Several T-dup+ familial chordoma patients had putative second/multiple primary chordomas. CONCLUSIONS The occurrence of young age at diagnosis, skull base presentation, or multiple primary chordomas should encourage careful review of family history for patients diagnosed with chordoma as well as screening of at-risk family members by MRI for early detection of chordoma. Furthermore, given genetic predisposition in some patients with familial chordoma, identification of a specific mutation in a family will permit surveillance to be limited to mutation carriers—and consideration should be given for imaging the entire neuraxis in any chordoma patient presenting at an early age or with a blood relative with chordoma. Finally, future studies should explore racial differences in age at diagnosis and presenting site in chordoma.

The application of next-generation sequencing (NGS) technologies in cancer research has accelerated the discovery of somatic mutations; however, progress in the identification of germline variation associated with cancer risk is less clear. We conducted a systematic literature review of cancer genetic susceptibility studies that used NGS technologies at an exome/genome-wide scale to obtain a fuller understanding of the research landscape to date and to inform future studies. The variability across studies on methodologies and reporting was considerable. Most studies sequenced few high-risk (mainly European) families, used a candidate analysis approach, and identified potential cancer-related germline variants or genes in a small fraction of the sequenced cancer cases. This review highlights the importance of establishing consensus on standards for the application and reporting of variants filtering strategies. It also describes the progress in the identification of cancer-related germline variation to date. These findings point to the untapped potential in conducting studies with appropriately sized and racially diverse families and populations, combining results across studies and expanding beyond a candidate analysis approach to advance the discovery of genetic variation that accounts for the unexplained cancer heritability.

Eleven high-evidence single-nucleotide polymorphisms (SNPs) at nine loci for gastric cancer (GC) risk were reported, but their associations with survival remain unknown. In this study, we examined associations between SNP and GC survival by anatomic location and histology among 1147 incident cases from the Shanxi Upper Gastrointestinal Genetics Project. We further examined whether SNPs were expression quantitative trait loci in normal and tumor gastric tissues, and whether tumor versus normal somatic mRNA differences in 126 cases were associated with survival. No SNPs were associated with GC survival overall. However, subtype-specific associations were observed for gastric cardia adenocarcinomas at MUC1/TRIM46/1q22 rs2070803 [HRAA versus GA+GG = 2.16; 95% confidence interval (CI) = 1.24-3.78; P = 0.0068] and LTA/TNF/6p21.33 rs1799724 (HRTT+CT versus CC = 1.30; 95% CI = 1.07-1.57; P = 0.0077), and for diffuse-type GC at PSCA/8q24.3 rs2294008 (HRTT versus CT+CC = 1.99; 95% CI = 1.33-2.97; P = 7.8E-04). Rs2294008T was a cis-expression quantitative trait loci for PSCA, upregulating mRNA in normal gastric (β = 0.60; P = 5.7E-21) and GC (β = 0.30; P = 0.0089) tissues. Cases in the highest quartile (the smallest downregulation of tumor PSCA) had shortest survival than cases with the most downregulated PSCA (median survival of 0.47 years in the highest quartile versus 3.73 years in the lowest quartile; hazard ratio = 9.70; 95% CI = 2.46-38.4; P = 0.0012). Less striking effects for mRNA levels were observed for MTX1 at 1q22 in gastric cardia adenocarcinoma and for JRK at 8q24.3 in diffuse GC. Our results suggest three high-evidence GC risk loci have prognostic importance in GC subtypes. Future studies in well-characterized independent populations are warranted to validate our findings and further investigate the clinical utility of these variants in predicting GC prognosis.

Based on our initial genome-wide association study (GWAS) on esophageal squamous cell carcinoma (ESCC) in Han Chinese, we conducted a follow-up study to examine the single nucleotide polymorphisms (SNPs) associated with family history (FH) of upper gastrointestinal cancer (UGI) cancer in cases with ESCC. We evaluated the association between SNPs and FH of UGI cancer among ESCC cases in a stage-1 case-only analysis of the National Cancer Institute (NCI, 541 cases with FH and 1399 without FH) and Henan GWAS (493 cases with FH and 869 without FH) data (discovery phase). The top SNPs (or their surrogates) from discovery were advanced to a stage-2 evaluation in additional Henan subjects (2801 cases with FH and 3136 without FH, replication phase). A total of 19 SNPs were associated with FH of UGI cancer in ESCC cases with P −5 in the stage-1 meta-analysis of NCI and Henan GWAS data. In stage-2, the association for rs79747906 (located at 18p11.31, P = 5.79 × 10−6 in discovery) was replicated (P = 0.006), with a pooled-OR of 1.59 (95%CI: 1.11-2.28). We identified potential genetic variants associated with FH of UGI cancer. Our findings may provide important insights into new low-penetrance susceptibility regions involved in the susceptibility of families with multiple UGI cancer cases.

Little is known about the link between gastric microbiota and the epidemiology of gastric cancer. In order to determine the epidemiologic and clinical relevance of gastric microbiota, we used 16 S ribosomal RNA gene sequencing analysis to characterize the composition and structure of the gastric microbial community of 80 paired samples (non-malignant and matched tumor tissues) from gastric cardia adenocarcinoma (GCA) patients in Shanxi, China. We also used PICRUSt to predict microbial functional profiles. Compared to patients without family history of upper gastrointestinal (UGI) cancer in the non-malignant gastric tissue microbiota, patients with family history of UGI cancer had higher Helicobacter pylori (Hp) relative abundance (median: 0.83 vs. 0.38, p = 0.01) and lower alpha diversity (median observed species: 51 vs. 85, p = 0.01). Patients with higher (vs. lower) tumor grade had higher Hp relative abundance (0.73 vs. 0.18, p = 0.03), lower alpha diversity (observed species, 66 vs. 89, p = 0.01), altered beta diversity (weighted UniFrac, p = 0.002) and significant alterations in relative abundance of five KEGG functional modules in non-malignant gastric tissue microbiota. Patients without metastases had higher relative abundance of Lactobacillales than patients with metastases (0.05 vs. 0.01, p = 0.04) in non-malignant gastric tissue microbiota. These associations were observed in non-malignant tissues but not in tumor tissues. In conclusion, this study showed a link of gastric microbiota to a major gastric cancer risk factor and clinical features in GCA patients from Shanxi, China. Studies with both healthy controls and gastric cardia and noncardia cancer cases across different populations are needed to further examine the association between gastric cancer and the microbiota.

Background: Although rare in the general population, highly penetrant germline mutations in CDKN2A are responsible for 5%-40% of melanoma cases reported in melanoma-prone families. We sought to determine whether MELPREDICT was generalizable to a global series of families with melanoma and whether performance improvements can be achieved.Methods: In total, 2116 familial melanoma cases were ascertained by the international GenoMEL Consortium. We recapitulated the MELPREDICT model within our data (GenoMELPREDICT) to assess performance improvements by adding phenotypic risk factors and history of pancreatic cancer. We report areas under the curve (AUC) with 95% confidence intervals (CIs) along with net reclassification indices (NRIs) as performance metrics.Results: MELPREDICT performed well (AUC 0.752, 95% CI 0.730-0.775), and GenoMELPREDICT performance was similar (AUC 0.748, 95% CI 0.726-0.771). Adding a reported history of pancreatic cancer yielded discriminatory improvement (P < .0001) in GenoMELPREDICT (AUC 0.772, 95% CI 0.750-0.793, NRI 0.40). Including phenotypic risk factors did not improve performance.Conclusion: The MELPREDICT model functioned well in a global data set of familial melanoma cases. Adding pancreatic cancer history improved model prediction. GenoMELPREDICT is a simple tool for predicting CDKN2A mutational status among melanoma patients from melanoma-prone families and can aid in directing these patients to receive genetic testing or cancer risk counseling.

Genetic susceptibility is likely involved in nasopharyngeal carcinoma (NPC), a cancer caused by Epstein-Barr virus (EBV) infection. Understanding of genetic factors involved in NPC and how they contribute to EBV-induced carcinogenesis is limited. We conducted whole-exome capture/sequencing among 251 individuals from 97 multiplex families from Taiwan (205 affected, 21 obligate carriers, and 25 unaffected) using SeqCap EZ Human Exome Library v3.0 and Illumina HiSeq. Aligned sequences were filtered to identify likely-to-be-functional deleterious variants that co-segregated with disease. Ingenuity Pathway analysis was performed. Circulating magnesium levels were measured in 13 individuals in 2 families with NIPAL1 mutations and in 197 sporadic NPC cases and 237 controls. We identified variants in 12 genes likely involved in cancer pathogenesis, viral infection or immune responses to infection. These included genes postulated to be involved in magnesium transport (NIPAL1), EBV cell entry (ITGB6), modulation of EBV infection (BCL2L12, NEDD4L), telomere biology (CLPTM1L, BRD2, HNRNPU), modulation of cAMP signaling (RAPGEF3), DNA repair (PRKDC, MLH1), and Notch signaling (NOTCH1, DLL3). Pathway based analysis demonstrated enrichment for Notch signaling genes (p-value = 0.0006). Evaluation of individuals within NIPAL1 families suggested lower serum magnesium in NPC compared to unaffected members. A significant reduction in serum magnesium levels was observed among sporadic NPC cases compared to controls (7.1% NPC/1.7% controls below normal range; OR = 4.5; 95% CI = 1.4,14) and is consistent with findings demonstrating a role for magnesium channeling in T-cell responses to EBV. We identified novel genes associated with NPC that point to new areas of inquiry to better understand genetic factors that determine the fate of viral infections and/or otherwise predisposes to NPC.

CDKN2A and CDK4 are well-established melanoma susceptibility genes, but their effect on tumor location and distribution is unknown. We used a case-case study design to assess for differences in tumor location between mutation carriers (CDKN2A = 141 patients, 348 melanomas; CDK4 = 15 patients, 54 melanomas) and noncarriers (104 patients, 157 melanomas) in US melanoma-prone families. Associations between groups were assessed with chi-square tests. Odds ratios (ORs) for tumor location were adjusted for diagnosis age, gender, and superficial spreading subtype. Models included random effects to account for within individual and family correlations. Compared with having a truncal melanoma, CDK4 (vs. noncarriers: lower extremities OR = 14.5, 95% confidence interval [CI] = 5.02-42.0, P0.001; upper extremities OR = 6.88, 95% CI = 2.37-19.9, P0.001; head and neck OR = 18.6, 95% CI = 4.04-85.2, P0.001) and CDKN2A (vs. noncarriers: lower extremities OR = 3.01, 95% CI = 1.56-5.82, P0.05; upper extremities OR = 1.91, 95% CI = 1.03-3.52, P 0.05; head and neck OR = 5.40, 95% CI = 2.10-13.9, P0.001) carriers had higher odds of developing melanoma at all other sites. Similar findings were observed for analyses stratified by gender, age, and first versus subsequent melanoma diagnoses. Further studies are needed to understand the biology underlying these genotype-associated patterns of tumor development, which could provide new insights into melanoma treatment and prevention.

Background Prior research has established that the prevalence of pathogenic/likely pathogenic (P/LP) variants across all of the American College of Medical Genetics (ACMG) Secondary Findings (SF) genes is approximately 0.8–5%. We investigated the prevalence of P/LP variants in the 24 ACMG SF v2.0 cancer genes in a family-based cancer research cohort (n = 1173) and in cancer-free ethnicity-matched controls (n = 982). Methods We used InterVar to classify variants and subsequently conducted a manual review to further examine variants of unknown significance (VUS). Results In the 24 genes on the ACMG SF v2.0 list associated with a cancer phenotype, we observed 8 P/LP unique variants (8 individuals; 0.8%) in controls and 11 P/LP unique variants (14 individuals; 1.2%) in cases, a non-significant difference. We reviewed 115 VUS. The median estimated per-variant review time required was 30 min; the first variant within a gene took significantly (p = 0.0009) longer to review (median = 60 min) compared with subsequent variants (median = 30 min). The concordance rate was 83.3% for the variants examined by two reviewers. Conclusion The 115 VUS required database and literature review, a time- and labor-intensive process hampered by the difficulty in interpreting conflicting P/LP determinations. By rigorously investigating the 24 ACMG SF v2.0 cancer genes, our work establishes a benchmark P/LP variant prevalence rate in a familial cancer cohort and controls. Electronic supplementary material The online version of this article (10.1186/s13073-018-0607-5) contains supplementary material, which is available to authorized users.

The risk of pancreatic cancer (PC) is increased in melanoma-prone families but the causal relationship between germline CDKN2A mutations and PC risk is uncertain, suggesting the existence of non-CDKN2A factors. One genetic possibility involves patients having mutations in multiple high-risk PC-related genes; however, no systematic examination has yet been conducted. We used next-generation sequencing data to examine 24 putative PC-related genes in 43 PC patients with and 23 PC patients without germline CDKN2A mutations and 1001 controls. For each gene and the four pathways in which they occurred, we tested whether PC patients (overall or CDKN2A+ and CDKN2A- cases separately) had an increased number of rare nonsynonymous variants. Overall, we identified 35 missense variants in PC patients, 14 in CDKN2A+ and 21 in CDKN2A- PC cases. We found nominally significant associations for mismatch repair genes (MLH1, MSH2, MSH6, PMS2) in all PC patients and for ATM, CPA1, and PMS2 in CDKN2A- PC patients. Further, nine CDKN2A+ and four CDKN2A- PC patients had rare potentially deleterious variants in multiple PC-related genes. Loss-of-function variants were only observed in CDKN2A- PC patients, with ATM having the most pathogenic variants. Also, ATM variants (n = 5) were only observed in CDKN2A- PC patients with a family history that included digestive system tumors. Our results suggest that a subset of PC patients may have increased risk because of germline mutations in multiple PC-related genes.

Hodgkin lymphoma shows strong familial aggregation but no major susceptibility genes have been identified to date. The goal of this study was to identify high-penetrance variants using whole exome sequencing in 17 Hodgkin lymphoma prone families with three or more affected cases or obligate carriers (69 individuals), followed by targeted sequencing in an additional 48 smaller HL families (80 individuals). Alignment and variant calling were performed using standard methods. Dominantly segregating, rare, coding or potentially functional variants were further prioritized based on predicted deleteriousness, conservation, and potential importance in lymphoid malignancy pathways. We selected 23 genes for targeted sequencing. Only the p.A1065T variant in KDR (kinase insert domain receptor) also known as VEGFR2 (vascular endothelial growth factor receptor 2) was replicated in two independent Hodgkin lymphoma families. KDR is a type III receptor tyrosine kinase, the main mediator of vascular endothelial growth factor induced proliferation, survival, and migration. Its activity is associated with several diseases including lymphoma. Functional experiments have shown that p.A1065T, located in the activation loop, can promote constitutive autophosphorylation on tyrosine in the absence of vascular endothelial growth factor and that the kinase activity was abrogated after exposure to kinase inhibitors. A few other promising mutations were identified but appear to be "private". In conclusion, in the largest sequenced cohort of Hodgkin lymphoma families to date, we identified a causal mutation in the KDR gene. While independent validation is needed, this mutation may increase downstream tumor cell proliferation activity and might be a candidate for targeted therapy.

Epstein–Barr virus (EBV) is an obligatory factor in the development of nasopharyngeal carcinoma (NPC), and anti-EBV IgA antibodies are elevated many years prior to the development of NPC. Nearly all adults are infected with EBV, but only a few develop cancer, suggesting that additional co-factors, including genetic susceptibility, must be required for the disease to manifest. Individuals were selected from the Taiwan Family Study, a cohort of 3389 individuals from NPC multiplex families. Primary analyses were conducted among 671 individuals from 69 pedigrees with the strongest family history of disease (>3 NPC-affected family members). The likelihood that a given family member carried a NPC susceptibility variant was estimated using Mendelian segregation rules, assuming a dominant mode of inheritance. We compared anti-EBV IgA antibody seropositivity between family members predicted to be carriers of NPC-linked genetic variants and those with a lower likelihood of carrying such variants. Obligate carriers of NPC susceptibility variants (100 % predicted probability of harbouring the genetic mutation) were nine-fold more likely to be anti-EBV IgA positive compared to family members predicted not to carry disease-causing variants (OR=9.2; P-trend

Most expression quantitative trait loci (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type specific regulatory landscape of human melanocytes, which give rise to melanoma but account for cis-eQTL SNPs prior to LD-pruning and 4,997 eGenes (FDRtrans-eQTLs with a pigmentation-associated SNP for four genes, likely through its cis-regulation of IRF4, encoding a transcription factor implicated in human pigmentation phenotypes. Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as melanoma-associated variants identified through genome-wide association studies (GWAS). Co-localization of melanoma GWAS variants and eQTLs from melanocyte and skin eQTL datasets identified candidate melanoma susceptibility genes for six known GWAS loci including unique genes identified by the melanocyte dataset. Further, a transcriptome-wide association study using published melanoma GWAS data uncovered four new loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings and present a robust database for genomic research of melanoma risk and melanocyte biology.

In a previous whole exome sequencing of patients from 41 families with Hodgkin lymphoma, we identified two families with distinct heterozygous rare coding variants in POT1 (D224N and Y36H), both in a highly conserved region of the gene. POT1 D224N mutant did not bind to a single-stranded telomere oligonucleotide in vitro suggesting the mutation perturbs POT1’s ability to bind to the telomeric G-rich overhang. Human HT1080 cells expressing POT1 D224N and lymphoblastoid cells carrying Y36H both showed increased telomere length and fragility in comparison to wild type cells. This strongly suggests that mutant POT1 causes chromosome instability and may play a role in lymphomagenesis in these families.

The International Agency for Research on Cancer (IARC) and the US National Cancer Institute (NCI) have initiated a series of cancer-focused seminars [Scelo G, Hofmann JN, Banks RE et al. International cancer seminars: a focus on kidney cancer. Ann Oncol 2016; 27(8): 1382-1385]. In this, the second seminar, IARC and NCI convened a workshop in order to examine the state of the current science on esophageal squamous cell carcinoma etiology, genetics, early detection, treatment, and palliation, was reviewed to identify the most critical open research questions. The results of these discussions were summarized by formulating a series of 'difficult questions', which should inform and prioritize future research efforts.

Known high-risk cutaneous malignant melanoma (CMM) genes account for melanoma risk in

© 2017 Yu, Torres, Hu, Medrano-Guzman, Herrera-Goepfert, Humphrys, Wang, Wang, Ding, Ravel, Taylor, Abnet and Goldstein. Helicobacter pylori (Hp) is the primary cause of gastric cancer but we know little of its relative abundance and other microbes in the stomach, especially at the time of gastric cancer diagnosis. Here we characterized the taxonomic and derived functional profiles of gastric microbiota in two different sets of gastric cancer patients, and compared them with microbial profiles in other body sites. Paired non-malignant and tumor tissues were sampled from 160 gastric cancer patients with 80 from China and 80 from Mexico. The 16S rRNA gene V3–V4 region was sequenced using MiSeq platform for taxonomic profiles. PICRUSt was used to predict functional profiles. Human Microbiome Project was used for comparison. We showed that Hp is the most abundant member of gastric microbiota in both Chinese and Mexican samples (51 and 24%, respectively), followed by oral-associated bacteria. Taxonomic (phylum-level) profiles of stomach microbiota resembled oral microbiota, especially when the Helicobacter reads were removed. The functional profiles of stomach microbiota, however, were distinct from those found in other body sites and had higher inter-subject dissimilarity. Gastric microbiota composition did not differ by Hp colonization status or stomach anatomic sites, but did differ between paired non-malignant and tumor tissues in either Chinese or Mexican samples. Our study showed that Hp is the dominant member of the non-malignant gastric tissue microbiota in many gastric cancer patients. Our results provide insights on the gastric microbiota composition and function in gastric cancer patients, which may have important clinical implications.

Mechanisms involved in wound healing play some role in carcinogenesis in multiple organs, likely by creating a chronic inflammatory milieu. This study sought to assess the role of genetic markers in selected inflammation-related genes involved in wound healing (interleukin (IL)-1a, IL-1b, IL-1 Receptor type I (IL-1Ra), IL-1 Receptor type II (IL-1Rb), tumour necrosis factor (TNF)-α, tumour necrosis factor receptor superfamily member (TNFRSF)1A, nuclear factor kappa beta (NF-kB)1, NF-kB2, inducible nitric oxide synthase (iNOS), cyclooxygenase (COX)-2, hypoxia induced factor (HIF)-1α, vascular endothelial growth factor (VEGF)A and P-53) in risk to oesophageal squamous cell carcinoma (OSCC).We genotyped 125 tag single nucleotide polymorphism (SNP)s in 410 cases and 377 age and sex matched disease-free individuals from Nutritional Intervention Trial (NIT) cohort, and 546 cases and 556 controls individually matched for age, sex and neighbourhood from Shanxi case-control study, both conducted in high-risk areas of north-central China (1985-2007). Cox proportional-hazard models and conditional logistic regression models were used for SNPs analyses for NIT and Shanxi, respectively. Fisher's inverse test statistics were used to obtain gene-level significance.Multiple SNPs were significantly associated with OSCC in both studies, however, none retained their significance after a conservative Bonferroni adjustment. Empiric p-values for tag SNPs in VEGFA in NIT were highly concentrated in the lower tail of the distribution, suggesting this gene may be influencing risk. Permutation tests confirmed the significance of this pattern. At the gene level, VEGFA yielded an empiric significance (P=0.027) in NIT. We also observed some evidence for interaction between environmental factors and some VEGFA tag SNPs.Our finding adds further evidence for a potential role for markers in the VEGFA gene in the development and progression of early precancerous lesions of oesophagus.