Germline pathogenic variants in neuroblastoma patients are enriched in BARD1 and predict worse survival

BACKGROUND

Neuroblastoma is an embryonal cancer of the developing sympathetic nervous system. The genetic contribution of rare pathogenic or likely pathogenic germline variants in patients without a family history remains unclear.

METHODS

Germline DNA sequencing was performed on 786 neuroblastoma patients. The frequency of rare cancer predisposition gene pathogenic or likely pathogenic variants in patients was compared with 2 cancer-free control cohorts. Matched tumor DNA sequencing was evaluated for second hits, and germline DNA array data from 5585 neuroblastoma patients and 23 505 cancer-free control children were analyzed to identify rare germline copy number variants. Patients with germline pathogenic or likely pathogenic variants were compared with those without to test for association with clinical characteristics, tumor features, and survival.

RESULTS

We observed 116 pathogenic or likely pathogenic variants involving 13.9% (109 of 786) of neuroblastoma patients, representing a statistically significant excess burden compared with cancer-free participants (odds ratio [OR] = 1.60, 95% confidence interval [CI] = 1.27 to 2.00). BARD1 harbored the most statistically significant enrichment of pathogenic or likely pathogenic variants (OR = 32.30, 95% CI = 6.44 to 310.35). Rare germline copy number variants disrupting BARD1 were identified in patients but absent in cancer-free participants (OR = 29.47, 95% CI = 1.52 to 570.70). Patients harboring a germline pathogenic or likely pathogenic variant had a worse overall survival compared with those without (P = 8.6 x 10-3).

CONCLUSIONS

BARD1 is an important neuroblastoma predisposition gene harboring both common and rare germline pathogenic or likely pathogenic variations. The presence of any germline pathogenic or likely pathogenic variant in a cancer predisposition gene was independently predictive of worse overall survival. As centers move toward paired tumor-normal sequencing at diagnosis, efforts should be made to centralize data and provide an infrastructure to support cooperative longitudinal prospective studies of germline pathogenic variation.

Human leukocyte antigen-DQA1*04:01 and rs2040406 variants are associated with elevated risk of childhood Burkitt lymphoma

Burkitt lymphoma (BL) is responsible for many childhood cancers in sub-Saharan Africa, where it is linked to recurrent or chronic infection by Epstein-Barr virus or Plasmodium falciparum. However, whether human leukocyte antigen (HLA) polymorphisms, which regulate immune response, are associated with BL has not been well investigated, which limits our understanding of BL etiology. Here we investigate this association among 4,645 children aged 0-15 years, 800 with BL, enrolled in Uganda, Tanzania, Kenya, and Malawi. HLA alleles are imputed with accuracy >90% for HLA class I and 85-89% for class II alleles. BL risk is elevated with HLA-DQA1*04:01 (adjusted odds ratio [OR] = 1.61, 95% confidence interval [CI] = 1.32-1.97, P = 3.71 × 10-6), with rs2040406(G) in HLA-DQA1 region (OR = 1.43, 95% CI = 1.26-1.63, P = 4.62 × 10-8), and with amino acid Gln at position 53 versus other variants in HLA-DQA1 (OR = 1.36, P = 2.06 × 10-6). The associations with HLA-DQA1*04:01 (OR = 1.29, P = 0.03) and rs2040406(G) (OR = 1.68, P = 0.019) persist in mutually adjusted models. The higher risk rs2040406(G) variant for BL is associated with decreased HLA-DQB1 expression in eQTLs in EBV transformed lymphocytes. Our results support the role of HLA variation in the etiology of BL and suggest that a promising area of research might be understanding the link between HLA variation and EBV control.

Sickle cell allele HBB-rs334(T) is associated with decreased risk of childhood Burkitt lymphoma in East Africa

Burkitt lymphoma (BL) is an aggressive B-cell lymphoma that significantly contributes to childhood cancer burden in sub-Saharan Africa. Plasmodium falciparum, which causes malaria, is geographically associated with BL, but the evidence remains insufficient for causal inference. Inference could be strengthened by demonstrating that mendelian genes known to protect against malaria-such as the sickle cell trait variant, HBB-rs334(T)-also protect against BL. We investigated this hypothesis among 800 BL cases and 3845 controls in four East African countries using genome-scan data to detect polymorphisms in 22 genes known to affect malaria risk. We fit generalized linear mixed models to estimate odds ratios (OR) and 95% confidence intervals (95% CI), controlling for age, sex, country, and ancestry. The ORs of the loci with BL and P. falciparum infection among controls were correlated (Spearman's ρ = 0.37, p = .039). HBB-rs334(T) was associated with lower P. falciparum infection risk among controls (OR = 0.752, 95% CI 0.628-0.9; p = .00189) and BL risk (OR = 0.687, 95% CI 0.533-0.885; p = .0037). ABO-rs8176703(T) was associated with decreased risk of BL (OR = 0.591, 95% CI 0.379-0.992; p = .00271), but not of P. falciparum infection. Our results increase support for the etiological correlation between P. falciparum and BL risk.

Factors associated with dementia attitudes in an adolescent cohort: Structural Equation Modelling

Objective

Dementia-related stigma is a key barrier to people living well with dementia, leading to social isolation and poor well-being. Adolescents represent an under-researched group that will become future carers and healthcare workers for the estimated 83.2 million people who will be living with dementia by 2030. Understanding the factors involved in dementia attitude formation in adolescents is useful for the development of evidence-based anti-stigma initiatives. This study aims to identify predictors of dementia attitudes in adolescents.

Methods

This is a cross-sectional study using secondary data analysis. 470 participants aged 12-15 years old from secondary schools in the Southeast of England, United Kingdom completed validated questionnaires relating to dementia attitudes (KIDS and Brief A-ADS) as well as demographic information. Multiple regressions were employed as well as a path analysis via a structural equation model to test for direct and mediatory effects.

Results

Multiple regression models revealed that being female, having higher levels of contact with dementia, and higher levels of empathy are positively associated with dementia attitudes in adolescents (p<0.05). Within the accepted structural equation model, empathy was a key mediator between contact and dementia attitudes.

Conclusion

This study highlights the pivotal role that contact with dementia can have in influencing dementia attitudes in adolescents with empathy serving as a mediator between contact and dementia attitudes. Interventions that use contact should consider how to stimulate empathetic responses to ultimately shape dementia attitudes.

Mosaic chromosomal alterations in peripheral blood leukocytes of children in sub-Saharan Africa

In high-income countries, mosaic chromosomal alterations in peripheral blood leukocytes are associated with an elevated risk of adverse health outcomes, including hematologic malignancies. We investigate mosaic chromosomal alterations in sub-Saharan Africa among 931 children with Burkitt lymphoma, an aggressive lymphoma commonly characterized by immunoglobulin-MYC chromosomal rearrangements, 3822 Burkitt lymphoma-free children, and 674 cancer-free men from Ghana. We find autosomal and X chromosome mosaic chromosomal alterations in 3.4% and 1.7% of Burkitt lymphoma-free children, and 8.4% and 3.7% of children with Burkitt lymphoma (P-values = 5.7×10-11 and 3.74×10-2, respectively). Autosomal mosaic chromosomal alterations are detected in 14.0% of Ghanaian men and increase with age. Mosaic chromosomal alterations in Burkitt lymphoma cases include gains on chromosomes 1q and 8, the latter spanning MYC, while mosaic chromosomal alterations in Burkitt lymphoma-free children include copy-neutral loss of heterozygosity on chromosomes 10, 14, and 16. Our results highlight mosaic chromosomal alterations in sub-Saharan African populations as a promising area of research.

Maker communities and the COVID-19 pandemic: a longitudinal analysis of Thingiverse's response to supply shortages

The COVID-19 pandemic profoundly affected various aspects of daily life, particularly the supply and demand of essential goods, resulting in critical shortages. This included personal protective equipment for medical professionals and the general public. To address these shortages, online 'maker communities' emerged, aiming to develop and locally manufacture critical products. While some organized efforts existed, the majority of initiatives originated from individuals and groups on platforms like Thingiverse. This paper presents a longitudinal analysis of Thingiverse, one of the largest maker community websites, to examine the pandemic's effects. Our findings reveal a surge in community output during the initial lockdown periods in major contributing nations (primarily those in the Western Hemisphere), followed by a subsequent decline. Additionally, throughout 2020, pandemic-related products dominated uploads and interactions during this period. Based on these observations, we propose recommendations to expedite the community's ability to support local, national and international responses to future disasters.

Genome-wide Association Study of Bladder Cancer Reveals New Biological and Translational Insights

BACKGROUND

Genomic regions identified by genome-wide association studies (GWAS) for bladder cancer risk provide new insights into etiology.

OBJECTIVE

To identify new susceptibility variants for bladder cancer in a meta-analysis of new and existing genome-wide genotype data.

DESIGN, SETTING, AND PARTICIPANTS

Data from 32 studies that includes 13,790 bladder cancer cases and 343,502 controls of European ancestry were used for meta-analysis.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSES

Log-additive associations of genetic variants were assessed using logistic regression models. A fixed-effects model was used for meta-analysis of the results. Stratified analyses were conducted to evaluate effect modification by sex and smoking status. A polygenic risk score (PRS) was generated on the basis of known and novel susceptibility variants and tested for interaction with smoking.

RESULTS AND LIMITATIONS

Multiple novel bladder cancer susceptibility loci (6p.22.3, 7q36.3, 8q21.13, 9p21.3, 10q22.1, 19q13.33) as well as improved signals in three known regions (4p16.3, 5p15.33, 11p15.5) were identified, bringing the number of independent markers at genome-wide significance (p < 5 × 10-8) to 24. The 4p16.3 (FGFR3/TACC3) locus was associated with a stronger risk for women than for men (p-interaction = 0.002). Bladder cancer risk was increased by interactions between smoking status and genetic variants at 8p22 (NAT2; multiplicative p value for interaction [pM-I] = 0.004), 8q21.13 (PAG1; pM-I = 0.01), and 9p21.3 (LOC107987026/MTAP/CDKN2A; pM-I = 0.02). The PRS based on the 24 independent GWAS markers (odds ratio per standard deviation increase 1.49, 95% confidence interval 1.44-1.53), which also showed comparable results in two prospective cohorts (UK Biobank, PLCO trial), revealed an approximately fourfold difference in the lifetime risk of bladder cancer according to the PRS (e.g., 1st vs 10th decile) for both smokers and nonsmokers.

CONCLUSIONS

We report novel loci associated with risk of bladder cancer that provide clues to its biological underpinnings. Using 24 independent markers, we constructed a PRS to stratify lifetime risk. The PRS combined with smoking history, and other established risk factors, has the potential to inform future screening efforts for bladder cancer.

PATIENT SUMMARY

We identified new genetic markers that provide biological insights into the genetic causes of bladder cancer. These genetic risk factors combined with lifestyle risk factors, such as smoking, may inform future preventive and screening strategies for bladder cancer.

Genome-wide association study of lung adenocarcinoma in East Asia and comparison with a European population

Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (Pinteraction = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications.

Abstract 3055: Associations of the mucosal microbiome and circulating bile acids with colorectal adenoma among average-risk women

Background: The gut microbiome (GM) and the metabolites it produces and regulates, such as bile acids (BAs), may act individually or interact to influence colorectal cancer development. We used data from a multi-center screening study to investigate the associations of the normal colon tissue microbiome and circulating BAs with colorectal adenoma.

Methods: We used data from the multi-center Colorectal Neoplasia Screening with Colonoscopy in Average-Risk Women Regional Navy/Army Medical Centers study (CONCeRN), comprising women undergoing complete colonoscopies. We individually matched women with adenoma (n=165 cases) to those without adenoma (n=311 controls). We extracted DNA using the Animal Tissue DNA Extraction Kit (AutoGene) and sequenced the V4 region of the 16S rRNA gene to characterize normal colonic mucosal bacteria, including alpha and beta diversity, taxonomic abundance, and co-abundance groups (CAGs). Fasting serum collected before colonoscopy underwent targeted quantitative liquid chromatography-tandem mass spectrometry to measure 13 primary and secondary BAs. We used multivariable conditional logistic regression, partial Pearson’s correlations, and linear regression to estimate interrelationships among the GM, BAs, and adenoma.

Results: Multiple genera were statistically significantly associated with adenomas. For example, for the centered log-ratio transformed relative abundance of Phascolarctobacterium, Parvimonas, and Porphymonas, the ORs (95% CIs) were 0.75 (0.60-0.95, P=0.02), 0.71 (0.55-0.92; P=0.01), and 0.67 (0.49-0.93; P=0.02), respectively. A CAG comprising phylum Firmicutes, Actinobacteria, Proteobacteria, and Fusobacteria was less abundant among cases than controls (log2folddiff=-0.80, P=0.002). Alpha and beta diversity associations were close to null and not statistically significant. However, alpha diversity metrics were inversely associated with primary BAs [e.g., for Shannon (R2=-0.13, P=0.01)] and positively associated with secondary BAs [e.g., for Shannon (R2=0.16, P=0.001)]. Though no circulating BAs were statistically significantly associated with adenomas, those in the highest relative to lowest tertile of taurine-conjugated (TC) BAs, primary BAs, and secondary BAs had a 44% (95% CI: 0.85-2.46), 24% (95% CI: 0.73-2.12), and 11% (95% CI: 0.61-2.01) higher odds of adenoma, respectively. We found that TCBAs, primary BAs, and secondary BAs were inversely associated with the abundance of Parvimonas [e.g., for TCBAs (R2=-0.11, P=0.02)] and Porphyromonas [e.g., for TCBAs (R2=-0.11, P=0.02)].

Conclusions: We found short-chain fatty acid-producing and oral-originating bacteria were inversely associated with adenomas. Our oral bacteria findings are contrary to findings from prior studies. Our findings suggest potential interrelationships among the GM and BAs in adenoma development.

Citation Format: Maria F. Gomez, Doratha A. Byrd, Stephanie R. Hogue, Jessica R. Burns, Nathan Smith, Joshua Sampson, Erikka Loftfield, Andrew Warner, Casey Dagnall, Kristine Jones, Belynda Hicks, Yunhu Wan, Youngchul Kim, Jin Xu, Rashmi Sinha, Emily Vogtmann. Associations of the mucosal microbiome and circulating bile acids with colorectal adenoma among average-risk women [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3055.

DNA methylation age in paired tumor and adjacent normal breast tissue in Chinese women with breast cancer

BACKGROUND

Few studies have examined epigenetic age acceleration (AA), the difference between DNA methylation (DNAm) predicted age and chronological age, in relation to somatic genomic features in paired cancer and normal tissue, with less work done in non-European populations. In this study, we aimed to examine DNAm age and its associations with breast cancer risk factors, subtypes, somatic genomic profiles including mutation and copy number alterations and other aging markers in breast tissue of Chinese breast cancer (BC) patients from Hong Kong.

METHODS

We performed genome-wide DNA methylation profiling of 196 tumor and 188 paired adjacent normal tissue collected from Chinese BC patients in Hong Kong (HKBC) using Illumina MethylationEPIC array. The DNAm age was calculated using Horvath's pan-tissue clock model. Somatic genomic features were based on data from RNA sequencing (RNASeq), whole-exome sequencing (WES), and whole-genome sequencing (WGS). Pearson's correlation (r), Kruskal-Wallis test, and regression models were used to estimate associations of DNAm AA with somatic features and breast cancer risk factors.

RESULTS

DNAm age showed a stronger correlation with chronological age in normal (Pearson r = 0.78, P < 2.2e-16) than in tumor tissue (Pearson r = 0.31, P = 7.8e-06). Although overall DNAm age or AA did not vary significantly by tissue within the same individual, luminal A tumors exhibited increased DNAm AA (P = 0.004) while HER2-enriched/basal-like tumors exhibited markedly lower DNAm AA (P = < .0001) compared with paired normal tissue. Consistent with the subtype association, tumor DNAm AA was positively correlated with ESR1 (Pearson r = 0.39, P = 6.3e-06) and PGR (Pearson r = 0.36, P = 2.4e-05) gene expression. In line with this, we found that increasing DNAm AA was associated with higher body mass index (P = 0.039) and earlier age at menarche (P = 0.035), factors that are related to cumulative exposure to estrogen. In contrast, variables indicating extensive genomic instability, such as TP53 somatic mutations, high tumor mutation/copy number alteration burden, and homologous repair deficiency were associated with lower DNAm AA.

CONCLUSIONS

Our findings provide additional insights into the complexity of breast tissue aging that is associated with the interaction of hormonal, genomic, and epigenetic mechanisms in an East Asian population.

Understanding the experience of time for dementia education programme on undergraduate radiography students

INTRODUCTION

Dementia is a global health priority, which requires the healthcare workforce to have the necessary attitudes and skills to deliver person-centred care to people with dementia. Radiographers have frequent contact with people with dementia, and undergraduate training is potentially an optimal time to deliver dementia education. Time for Dementia is an education programme in which undergraduate healthcare students visit a person with dementia and their carer over a two-year period to gain an in-depth understanding of the condition. The aim of this study was to understand undergraduate radiography students' experiences of undertaking the Time for Dementia (TFD) programme.

METHODS

Two focus groups were undertaken with 14 radiography students who were half-way through the TFD programme. Data was analysed using thematic analysis.

RESULTS

Three key themes were constructed from the analysis: A Holistic Learning Experience, Transferring Learning into Practice and Preparedness & Expectations. Participants discussed the value from learning directly from people with dementia and their carers, reporting an increase in their awareness and understanding of dementia as well as the impact of caring for somebody with the condition. Participants were able to identify learning to take into practice such as person-centred care, compassion, and patience. Challenges to learning were also identified.

CONCLUSIONS

This study suggests that a longitudinal, experiential education programme provides radiography students with the opportunity to develop a more holistic understanding of dementia and the impact it may have on the individual and their family members.

IMPLICATIONS FOR PRACTICE

Experiential dementia teaching is of value to radiography students, however preparation and learning support should fit with previous personal and teaching experience.

Prospective and Cross-sectional Associations of the Rectal Tissue Microbiome with Colorectal Adenoma Recurrence

BACKGROUND

The gut microbiome is plausibly associated with colorectal cancer risk; however, previous studies mostly investigated this association cross-sectionally. We investigated cross-sectional and prospective associations of the rectal tissue microbiome with adenoma recurrence in the Polyp Prevention Trial (PPT).

METHODS

PPT is a 4-year randomized clinical trial of the effect of a dietary intervention on adenoma recurrence among community members. We extracted DNA from rectal biopsies at baseline, end of year 1, and end of year 4 among 455 individuals and sequenced the V4 region of the 16S rRNA gene. At each timepoint, we investigated associations of alpha diversity, beta diversity, and presence and relative abundance of select taxa with adenoma recurrence using multivariable logistic regression.

RESULTS

Variation in beta diversity was primarily explained by subject and minimally by year of collection or time between biopsy and colonoscopy. Cross-sectionally, year 4 alpha diversity was strongly, inversely associated with adenoma prevalence [ORQ3 vs. Q1 Shannon index = 0.40 (95% confidence interval, CI: 0.21-0.76)]. Prospective alpha diversity associations (i.e., baseline/year 1 alpha diversity with adenoma recurrence 3-4 years later) were weak or null, as were cross-sectional and prospective beta diversity-adenoma associations. Bacteroides abundance was more strongly, positively associated with adenoma prevalence cross-sectionally than prospectively.

CONCLUSIONS

Rectal tissue microbiome profiles may be associated with prevalent adenomas, with little evidence supporting prospective associations.

IMPACT

Additional prospective studies, with serial fecal and tissue samples, to explore microbiome-colorectal cancer associations are needed. Eventually, it may be possible to use microbiome characteristics as intervenable risk factors or screening tools.

Gene Expression Profiling Identifies Two Chordoma Subtypes Associated with Distinct Molecular Mechanisms and Clinical Outcomes

PURPOSE

Chordoma is a rare bone tumor with a high recurrence rate and limited treatment options. The aim of this study was to identify molecular subtypes of chordoma that may improve clinical management.

EXPERIMENTAL DESIGN

We conducted RNA sequencing in 48 tumors from patients with Chinese skull-base chordoma and identified two major molecular subtypes. We then replicated the classification using a NanoString panel in 48 patients with chordoma from North America.

RESULTS

Tumors in one subtype were more likely to have somatic mutations and reduced expression in chromatin remodeling genes, such as PBRM1 and SETD2, whereas the other subtype was characterized by the upregulation of genes in epithelial-mesenchymal transition and Sonic Hedgehog pathways. IHC staining of top differentially expressed genes between the two subtypes in 312 patients with Chinese chordoma with long-term follow-up data showed that the expression of some markers such as PTCH1 was significantly associated with survival outcomes.

CONCLUSIONS

Our findings may improve the understanding of subtype-specific tumorigenesis of chordoma and inform clinical prognostication and targeted options.

Gene-Level Associations in Patients With and Without Pathogenic Germline Variants in CDKN2A and Pancreatic Cancer

PURPOSE

Pancreatic ductal adenocarcinoma (PDAC) is a component of familial melanoma due to germline pathogenic variants (GPVs) in CDKN2A. However, it is unclear what role this gene or other genes play in its etiology.

MATERIALS AND METHODS

We analyzed 189 cancer predisposition genes using parametric rare-variant association (RVA) tests and nonparametric permutation tests to identify gene-level associations in PDAC for patients with (CDKN2A+) and without (CDKN2A-) GPV. Exome sequencing was performed on 84 patients with PDAC, 47 CDKN2A+ and 37 CDKN2A-. After variant filtering, various RVA tests and permutation tests were run separately by CDKN2A status. Genes with the strongest nominal associations were evaluated in patients with PDAC from The Cancer Genome Atlas and the UK Biobank (UKB). A secondary analysis including only GPV from UKB was also performed.

RESULTS

In RVA tests, ERCC4 and RET showed the most compelling evidence as plausible PDAC candidate genes for CDKN2A+ patients. In contrast, the findings in CDKN2A- patients provided evidence for HMBS, EPCAM, and MRE11 as potential new candidate genes and confirmed ATM, BRCA2, and PALB2 as PDAC genes, consistent with findings in The Cancer Genome Atlas and the UKB. As expected, CDKN2A- patients were more likely to harbor GPVs from the 189 genes investigated. When including only GPVs from UKB, significant associations with PDAC were seen for ATM, BRCA2, and CDKN2A.

CONCLUSION

These results suggest that variants in other genes likely play a role in PDAC in all patients and that PDAC in CDKN2A+ patients has a distinct etiology from PDAC in CDKN2A- patients.

Rare germline variants in PALB2 and BRCA2 in familial and sporadic chordoma

Chordoma is a rare bone tumor with genetic risk factors largely unknown. We conducted a whole-exome sequencing (WES) analysis of germline DNA from 19 familial chordoma cases in five pedigrees and 137 sporadic chordoma patients and identified 17 rare germline variants in PALB2 and BRCA2, whose products play essential roles in homologous recombination (HR) and tumor suppression. One PALB2 variant showed disease cosegregation in a family with four affected people or obligate gene carrier. Chordoma cases had a significantly increased burden of rare variants in both genes when compared to population-based controls. Four of the six PALB2 variants identified from chordoma patients modestly affected HR function and three of the 11 BRCA2 variants caused loss of function in experimental assays. These results, together with previous reports of abnormal morphology and Brachyury expression of the notochord in Palb2 knockout mouse embryos and genomic signatures associated with HR defect and HR gene mutations in advanced chordomas, suggest that germline mutations in PALB2 and BRCA2 may increase chordoma susceptibility. Our data shed light on the etiology of chordoma and support the previous finding that PARP-1 inhibitors may be a potential therapy for some chordoma patients.

Integrated Analysis of Coexpression and Exome Sequencing to Prioritize Susceptibility Genes for Familial Cutaneous Melanoma

The application of whole-exome sequencing has led to the identification of high- and moderate-risk variants that contribute to cutaneous melanoma susceptibility. However, confirming disease-causing variants remains challenging. We applied a gene coexpression network analysis to prioritize the candidate genes identified from whole-exome sequencing of 34 melanoma-prone families, with at least three affected members sequenced per family (N = 119 cases). A coexpression network was constructed from genotype-tissue expression project, skin melanoma from the cancer genome atlas, and primary melanocyte cultures. We performed module-specific enrichment and focused on modules associated with pigmentation processes because they are the best-studied and most well-known risk factors for melanoma susceptibility. We found that pigmentation-associated modules across the four expression datasets examined were enriched for well-known melanoma susceptibility genes plus genes associated with pigmentation. We also used network properties to prioritize genes within pigmentation modules as candidate susceptibility genes. Integrating information from coexpression network analysis and variant prioritization, we identified 36 genes (such as DCT, TPCN2, TRPM1, ATP10A, and EPHA5) as potential melanoma risk genes in the families. Our approach also allowed us to link families with private gene mutations on the basis of gene coexpression patterns and thereby may provide an innovative perspective in gene identification in high-risk families.

Genetic testing in severe aplastic anemia is required for optimal hematopoietic cell transplant outcomes

Patients with severe aplastic anemia (SAA) can have an unrecognized inherited bone marrow failure syndrome (IBMFS) because of phenotypic heterogeneity. We curated germline genetic variants in 104 IBMFS-associated genes from exome sequencing performed on 732 patients who underwent hematopoietic cell transplant (HCT) between 1989 and 2015 for acquired SAA. Patients with pathogenic or likely pathogenic (P/LP) variants fitting known disease zygosity patterns were deemed unrecognized IBMFS. Carriers were defined as patients with a single P/LP variant in an autosomal recessive gene or females with an X-linked recessive P/LP variant. Cox proportional hazard models were used for survival analysis with follow-up until 2017. We identified 113 P/LP single-nucleotide variants or small insertions/deletions and 10 copy number variants across 42 genes in 121 patients. Ninety-one patients had 105 in silico predicted deleterious variants of uncertain significance (dVUS). Forty-eight patients (6.6%) had an unrecognized IBMFS (33% adults), and 73 (10%) were carriers. No survival difference between dVUS and acquired SAA was noted. Compared with acquired SAA (no P/LP variants), patients with unrecognized IBMFS, but not carriers, had worse survival after HCT (IBMFS hazard ratio [HR], 2.13; 95% confidence interval[CI], 1.40-3.24; P = .0004; carriers HR, 0.96; 95% CI, 0.62-1.50; P = .86). Results were similar in analyses restricted to patients receiving reduced-intensity conditioning (n = 448; HR IBMFS = 2.39; P = .01). The excess mortality risk in unrecognized IBMFS attributed to death from organ failure (HR = 4.88; P < .0001). Genetic testing should be part of the diagnostic evaluation for all patients with SAA to tailor therapeutic regimens. Carriers of a pathogenic variant in an IBMFS gene can follow HCT regimens for acquired SAA.

PLCOjs, a FAIR GWAS web SDK for the NCI Prostate, Lung, Colorectal and Ovarian Cancer Genetic Atlas project

MOTIVATION

The Division of Cancer Epidemiology and Genetics (DCEG) and the Division of Cancer Prevention (DCP) at the National Cancer Institute (NCI) have recently generated genome-wide association study (GWAS) data for multiple traits in the Prostate, Lung, Colorectal, and Ovarian (PLCO) Genomic Atlas project. The GWAS included 110 000 participants. The dissemination of the genetic association data through a data portal called GWAS Explorer, in a manner that addresses the modern expectations of FAIR reusability by data scientists and engineers, is the main motivation for the development of the open-source JavaScript software development kit (SDK) reported here.

RESULTS

The PLCO GWAS Explorer resource relies on a public stateless HTTP application programming interface (API) deployed as the sole backend service for both the landing page's web application and third-party analytical workflows. The core PLCOjs SDK is mapped to each of the API methods, and also to each of the reference graphic visualizations in the GWAS Explorer. A few additional visualization methods extend it. As is the norm with web SDKs, no download or installation is needed and modularization supports targeted code injection for web applications, reactive notebooks (Observable) and node-based web services.

AVAILABILITY AND IMPLEMENTATION

code at https://github.com/episphere/plco; project page at https://episphere.github.io/plco.

Comparison of somatic mutation landscapes in Chinese versus European breast cancer patients

Recent genomic studies suggest that Asian breast cancer (BC) may have distinct somatic features; however, most comparisons of BC genomic features across populations did not account for differences in age, subtype, and sequencing methods. In this study, we analyzed whole-exome sequencing (WES) data to characterize somatic copy number alterations (SCNAs) and mutation profiles in 98 Hong Kong BC (HKBC) patients and compared with those from The Cancer Genome Atlas of European ancestry (TCGA-EA, N = 686), which had similar distributions of age at diagnosis and PAM50 subtypes as in HKBC. We developed a two-sample Poisson model to compare driver gene selection pressure, which reflects the effect sizes of cancer driver genes, while accounting for differences in sample size, sequencing platforms, depths, and mutation calling methods. We found that somatic mutation and SCNA profiles were overall very similar between HKBC and TCGA-EA. The selection pressure for small insertions and deletions (indels) in GATA3 (false discovery rate (FDR) corrected p < 0.01) and single-nucleotide variants (SNVs) in TP53 (nominal p = 0.02, FDR corrected p = 0.28) was lower in HKBC than in TCGA-EA. Among the 13 signatures of single-base substitutions (SBS) that are common in BC, we found a suggestively higher contribution of SBS18 and a lower contribution of SBS1 in HKBC than in TCGA-EA, while the two APOBEC-induced signatures showed similar prevalence. Our results suggest that the genomic landscape of BC was largely very similar between HKBC and TCGA-EA, despite suggestive differences in some driver genes and mutational signatures that warrant future investigations in large and diverse Asian populations.

Germline ATM variants predispose to melanoma: a joint analysis across the GenoMEL and MelaNostrum consortia

PURPOSE

Ataxia-Telangiectasia Mutated (ATM) has been implicated in the risk of several cancers, but establishing a causal relationship is often challenging. Although ATM single-nucleotide polymorphisms have been linked to melanoma, few functional alleles have been identified. Therefore, ATM impact on melanoma predisposition is unclear.

METHODS

From 22 American, Australian, and European sites, we collected 2,104 familial, multiple primary (MPM), and sporadic melanoma cases who underwent ATM genotyping via panel, exome, or genome sequencing, and compared the allele frequency (AF) of selected ATM variants classified as loss-of-function (LOF) and variants of uncertain significance (VUS) between this cohort and the gnomAD non-Finnish European (NFE) data set.

RESULTS

LOF variants were more represented in our study cohort than in gnomAD NFE, both in all (AF = 0.005 and 0.002, OR = 2.6, 95% CI = 1.56-4.11, p < 0.01), and familial + MPM cases (AF = 0.0054 and 0.002, OR = 2.97, p < 0.01). Similarly, VUS were enriched in all (AF = 0.046 and 0.033, OR = 1.41, 95% CI = 1.6-5.09, p < 0.01) and familial + MPM cases (AF = 0.053 and 0.033, OR = 1.63, p < 0.01). In a case-control comparison of two centers that provided 1,446 controls, LOF and VUS were enriched in familial + MPM cases (p = 0.027, p = 0.018).

CONCLUSION

This study, describing the largest multicenter melanoma cohort investigated for ATM germline variants, supports the role of ATM as a melanoma predisposition gene, with LOF variants suggesting a moderate-risk.

Cigarette Smoking and Opium Use in Relation to the Oral Microbiota in Iran

Cigarettes and opium contain chemicals and particulate matter that may modify the oral microbiota. This study aimed to investigate the association between cigarette and opium use with the oral microbiota. A total of 558 participants were recruited from Iran between 2011 and 2015. Individuals were categorized as never cigarette nor opium users, ever cigarette-only smokers, ever opium-only users, and ever both cigarette and opium users. Participants provided saliva samples for 16S rRNA gene sequencing. Logistic regression, microbiome regression-based kernel association test (MiRKAT), and zero-inflated beta regression models were calculated. For every increase in 10 observed amplicon sequence variants (ASVs), the odds for being a cigarette-only smoker, opium-only user, and both user compared to never users decreased by 9% (odds ratio [OR] = 0.91; 95% confidence interval [95% CI] = 0.86 to 0.97), 13% (OR = 0.87; 95% CI = 0.75 to 1.01), and 12% (OR = 0.88; 95% CI = 0.80 to 0.96), respectively. The microbial communities differed by cigarette and opium use as indicated by MiRKAT models testing the three beta-diversity matrices (P < 0.05 for all). Three genera were less likely and one genus was more likely to be detected in cigarette-only smokers or opium-only users than in never users. The relative abundance of the phylum Actinobacteria (never, 14.78%; both, 21.20%) was higher and the phyla Bacteroidetes (never, 17.63%; both, 11.62%) and Proteobacteria (never, 9.06%; both, 3.70%) were lower in users of both cigarettes and opium, while the phylum Firmicutes (never, 54.29%; opium, 65.49%) was higher in opium-only users. Cigarette and opium use was associated with lower alpha-diversity, overall oral microbiota community composition, and both the presence and relative abundance of multiple taxa. IMPORTANCE Cigarette smoking and opium use are associated with periodontal disease caused by specific bacteria such as Porphyromonas gingivalis, which suggests a link between cigarette smoking and opium use and the oral microbiota. Alterations of the oral microbiota in cigarette smokers compared to nonsmokers have been reported, but this has not been studied across diverse populations. Additionally, the association of opium use with the oral microbiota has not been investigated to date. We conducted this study to investigate differences in the oral microbiota between ever users of cigarettes only, opium only, and both cigarettes and opium and never users of cigarettes and opium in Iran. Lower alpha-diversity, distinct overall oral microbial communities, and the presence and relative abundance of multiple taxa have been found for users of cigarettes and/or opium.

Genomic Classification and Clinical Outcome in Rhabdomyosarcoma: A Report From an International Consortium

PURPOSE

Rhabdomyosarcoma is the most common soft tissue sarcoma of childhood. Despite aggressive therapy, the 5-year survival rate for patients with metastatic or recurrent disease remains poor, and beyond PAX-FOXO1 fusion status, no genomic markers are available for risk stratification. We present an international consortium study designed to determine the incidence of driver mutations and their association with clinical outcome.

PATIENTS AND METHODS

Tumor samples collected from patients enrolled on Children's Oncology Group trials (1998-2017) and UK patients enrolled on malignant mesenchymal tumor and RMS2005 (1995-2016) trials were subjected to custom-capture sequencing. Mutations, indels, gene deletions, and amplifications were identified, and survival analysis was performed.

RESULTS

DNA from 641 patients was suitable for analyses. A median of one mutation was found per tumor. In FOXO1 fusion-negative cases, mutation of any RAS pathway member was found in > 50% of cases, and 21% had no putative driver mutation identified. BCOR (15%), NF1 (15%), and TP53 (13%) mutations were found at a higher incidence than previously reported and TP53 mutations were associated with worse outcomes in both fusion-negative and FOXO1 fusion-positive cases. Interestingly, mutations in RAS isoforms predominated in infants < 1 year (64% of cases). Mutation of MYOD1 was associated with histologic patterns beyond those previously described, older age, head and neck primary site, and a dismal survival. Finally, we provide a searchable companion database (ClinOmics), containing all genomic variants, and clinical annotation including survival data.

CONCLUSION

This is the largest genomic characterization of clinically annotated rhabdomyosarcoma tumors to date and provides prognostic genetic features that refine risk stratification and will be incorporated into prospective trials.

Tracing Lung Cancer Risk Factors Through Mutational Signatures in Never-Smokers

Epidemiologic studies often rely on questionnaire data, exposure measurement tools, and/or biomarkers to identify risk factors and the underlying carcinogenic processes. An emerging and promising complementary approach to investigate cancer etiology is the study of somatic "mutational signatures" that endogenous and exogenous processes imprint on the cellular genome. These signatures can be identified from a complex web of somatic mutations thanks to advances in DNA sequencing technology and analytical algorithms. This approach is at the core of the Sherlock-Lung study (2018-ongoing), a retrospective case-only study of over 2,000 lung cancers in never-smokers (LCINS), using different patterns of mutations observed within LCINS tumors to trace back possible exposures or endogenous processes. Whole genome and transcriptome sequencing, genome-wide methylation, microbiome, and other analyses are integrated with data from histological and radiological imaging, lifestyle, demographic characteristics, environmental and occupational exposures, and medical records to classify LCINS into subtypes that could reveal distinct risk factors. To date, we have received samples and data from 1,370 LCINS cases from 17 study sites worldwide and whole-genome sequencing has been completed on 1,257 samples. Here, we present the Sherlock-Lung study design and analytical strategy, also illustrating some empirical challenges and the potential for this approach in future epidemiologic studies.

Rare Germline Variants in Chordoma-Related Genes and Chordoma Susceptibility

Simple Summary

Chordoma is an extremely rare bone cancer that has not been fully characterized and few risk factors have been identified, highlighting the need for improving our understanding of the disease biology. Our study aims to identify chordoma susceptibility genes by investigating 265 genes involved in chordoma-related signaling pathways and other biological processes on germline DNA of 138 chordoma patients of European ancestry compared to internal control datasets and general population databases. Results were intersected with whole genome sequencing data from 80 skull-base chordoma patients of Chinese ancestry. Several rare loss-of-function and predicted deleterious missense variants were enriched in chordoma cases in both datasets, suggesting a complex model of pathways potentially involved in chordoma development and susceptibility, warranting further investigation in larger studies.

Abstract

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.

Abstract PO-055: Molecular characterization of papillary thyroid cancer in relation to ionizing radiation dose following the Chernobyl accident

The 1986 Chernobyl nuclear power plant accident increased papillary thyroid cancer (PTC) incidence in surrounding regions, particularly for 131I-exposed children. To investigate the contribution of environmental radiation to PTC characteristics and improve understanding of radiation-induced carcinogenesis, we analyzed genomic, transcriptomic, and epigenomic characteristics of 440 pathologically-confirmed fresh-frozen PTCs from Ukraine (359 with estimated childhood or in utero 131I exposure and 81 from unexposed children born after March 1987) and matched normal tissue (non-tumor thyroid tissue and/or blood). Mean age at PTC was 28.0 years (range: 10.0-45.6). Among 131I-exposed individuals, mean radiation dose was 250 mGy (range: 11.0-8,800). In multivariable models adjusted for age at PTC and sex, we observed radiation dose-dependent enrichment of fusion drivers (P=6.6 × 10−8), nearly all occurring in the MAPK pathway, as well as increases in small deletions (P=8.0 × 10−9) and simple/balanced structural variants (P=1.2 × 10−14). Further analyses demonstrated even stronger associations for those small deletions and simple/balanced structural variants that were clonal and bore hallmarks of non-homologous end-joining repair (deletions: P=4.9 × 10−31; simple/balanced structural variants: P=5.5 × 10−19). In contrast, radiation dose was not associated with subclonal small deletions (P=0.82) or subclonal simple/balanced structural variants (P=0.91). Additionally, radiation dose was not associated with TINS (locally templated insertions), which are characteristic of alt-end-joining repair (P=0.69). The effects of radiation on genomic alterations with more pronounced for those younger at exposure. Analyses generally were consistent with a linear radiation dose-response for all molecular characteristics except clonal small deletions. Analyses of transcriptomic and epigenomic features demonstrated strong associations with the PTC driver gene but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth following environmental radiation exposure.

Citation Format: Lindsay M. Morton, Danielle Karyadi, Chip Stewart, Tetiana Bogdanova, Eric Dawson, Mia Steinberg, Jieqiong Dai, Stephen Hartley, Sara Schonfeld, Joshua Sampson, Yosi Maruvka, Vidushi Kapoor, Dale Ramsden, Juan Carvajal-Garcia, Chuck Perou, Joel Parker, Marko Krznaric, Meredith Yeager, Joseph Boland, Amy Hutchinson, Belynda Hicks, Casey Dagnall, Julie Gastier-Foster, Jay Bowen, Olivia Lee, Mitchell Machiela, Elizabeth Cahoon, Alina Brenner, Kiyohiko Mabuchi, Vladimir Drozdovitch, Sergii Masiuk, Mykola Chepurny, Liudmyla Yu Zurnadzhy, Maureen Hatch, Amy Berrington de Gonzalez, Gerry Thomas, Mykola Tronko, Gad Getz, Stephen Chanock. Molecular characterization of papillary thyroid cancer in relation to ionizing radiation dose following the Chernobyl accident [abstract]. In: Proceedings of the AACR Virtual Special Conference on Radiation Science and Medicine; 2021 Mar 2-3. Philadelphia (PA): AACR; Clin Cancer Res 2021;27(8_Suppl):Abstract nr PO-055.

Whole genome sequencing of skull-base chordoma reveals genomic alterations associated with recurrence and chordoma-specific survival

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.

Using whole-exome sequencing and protein interaction networks to prioritize candidate genes for germline cutaneous melanoma susceptibility

Although next-generation sequencing has demonstrated great potential for novel gene discovery, confirming disease-causing genes after initial discovery remains challenging. Here, we applied a network analysis approach to prioritize candidate genes identified from whole-exome sequencing analysis of 98 cutaneous melanoma patients from 27 families. Using a network propagation method, we ranked candidate genes by their similarity to known disease genes in protein-protein interaction networks and identified gene clusters with functional connectivity. Using this approach, we identified several new candidate susceptibility genes that warrant future investigations such as NGLY1, IL1RN, FABP2, PRKDC, and PROSER2. The propagated network analysis also allowed us to link families that did not have common underlying genes but that carried variants in genes that interact on protein-protein interaction networks. In conclusion, our study provided an analysis perspective for gene prioritization in the context of genetic heterogeneity across families and prioritized top potential candidate susceptibility genes in our dataset.

Age-related DNA methylation in paired normal and tumour breast tissue in Chinese breast cancer patients

Age-related DNA methylation is a potential mechanism contributing to breast cancer development. Studies of primarily Caucasian women have identified many CpG sites of age-related methylation in non-diseased breast tissue possibly driving cancer development over time. There is a paucity of studies involving Asian women whose ages at breast cancer onset are usually younger than Caucasians. We identified the 181 most consistent age-related methylation events in non-diseased breast tissue across published studies. Age-related methylation events were measured in adjacent normal and breast tumour tissue in an exclusively Asian population at the previously identified age-related methylation sites. Age-related methylation was found in 118 probes in adjacent normal breast tissue. Methylation of 99% of these sites was increased with age and predominantly located on CpG islands in promoter regions. To ascertain biological relevance to breast cancer, we focused on the 37 sites with overall higher methylation in tumour compared to adjacent normal samples. Some sites positively related to age, including AQP5 and CORO6, inversely correlated with gene expression. Several others have known involvement in suppression of carcinogenesis including GPC5 and SST, suggesting that perturbation of epigenetic regulation at these sites due to ageing may contribute to the progression of carcinogenesis. This study highlights an age-related methylation landscape in non-tumour tissue, consistent not just across studies, but also across different populations. We present candidate age-related methylation sites warranting further investigation as potential epigenetic drivers of breast cancer. They may serve as potential targets of site-specific demethylation intervention strategies for the prevention of age-related breast cancer.

Differential susceptibility to allostatic load and educational inequalities in coronary heart disease

Background

Differential exposure to lifestyle factors may mediate the association between education and cardiovascular disease (CVD). However, differential susceptibility (the effect of exposure to the same “dose” of risk factors differs across groups) may also elevate CVD risk but the causal pathways remain unclear. Allostatic Load (AL) is a marker of cumulative biological burden resulting from mal-adaptation to chronic stressors. We aimed to examine the role of differential exposure and susceptibility to AL and other factors in coronary heart disease (CHD) educational gradients in Europe.

Methods

51,328 35-74-year-old participants originally free of CVD from 21 European cohorts in the BiomarCaRE consortium were identified and followed for a median of 10 years to their first CHD event. We defined an AL score as the sum of z-scores of 8 markers from the cardiovascular, metabolic, and inflammatory systems. To investigate the mediating role of AL (and smoking, alcohol and BMI) on educational differences in CHD incidence we applied marginal structural models and three-way decomposition on gender-specific additive hazards models.

Results

AL was a significant mediator of the association between educational status and CHD. The highest proportion mediated was observed in women, with 28% (95%CI 20% to 44%) attributable to differential exposure and 8% (95%CI 0% to 16%) to differential susceptibility. In men, AL mediated 16% of the increased CHD risk in the less educated, with 2% (95%CI 0%-6%) attributable to differential susceptibility. The effects of smoking, alcohol and BMI were relatively small for men and women, with a limited role of differential susceptibility.

Conclusions

While we found evidence of differential susceptibility to AL on CHD, effects were modest and the mediating effect of AL (and other lifestyle factors) was predominately via differential exposure. Controlling disproportionate exposure to AL may help reduce CHD morbidity among those with lower education.

Key messages

Educational inequalities in coronary heart disease in Europe were predominantly driven by differential exposure rather than susceptibility to allostatic load and other lifestyle factors. Controlling disproportionate exposure to AL may help reduce coronary heart disease morbidity among those with lower education.

Abstract 146: DNA methylation age of paired tumor-normal breast tissue in Chinese women with breast cancer

Background: Epigenetic age, captured by DNA methylation changes, is thought to be more informative with respect to disease risk and progression, than chronological age. A few studies have associated epigenetic age acceleration (EAA), the difference between DNA methylation age (DNAm age) and chronological age, with BC risk and subtypes. However, most of these studies were conducted in Western populations. In this study, we examined EAA in an Asian population, in which BC incidence rate is lower and age at BC onset is younger compared to most Western populations.

Methods: We performed genome-wide DNA methylation profiling of 97 tumor and 89 paired distant normal tissue samples collected from BC patients in Hong Kong (HK) using an Illumina MethylationEPIC array. Two independent datasets were used to compare results: The Cancer Genome Atlas (TCGA, n=525 tumor and n=88 adjacent normal) and healthy women from the Komen tissue bank (n=59). DNAm age was calculated using Horvath's model based on 353 CpGs. The significance of EAA was tested using a simple linear regression model. We also used a multivariate regression model to test for equality of slopes (EAA rates) among different BC subtypes or datasets, considering interaction terms between age and subtype/dataset.

Results: The average age at BC diagnosis was 58 years old (range: 33-81) and the distribution of the molecular subtype based on PAM50 was 43.0%, 29.1%, and 27.9% for luminal A, luminal B, and HER2-enriched/basal tumors, respectively. As expected, DNAm age showed a stronger correlation with chronological age in normal tissue (r=0.78, p&lt;.0001) than in tumor tissue (r=0.27, p=0.0075). The average EAAs in normal and tumor samples were 9.43 and 6.02 years, respectively. Among different BC subtypes, EAA in normal tissue did not vary by subtypes. However, in tumor samples, luminal patients showed positive EAA (average 10/13 years in luminal A/B, respectively), while HER2-enriched/basal patients showed a negative EAA (average -8 years), although the rate of EAA did not vary significantly by subtype. Analyses on TCGA data produced consistent results. When comparing the rate of DNAm acceleration in normal tissue of HK, TCGA, and Komen, HK patients had a significantly different rate (βHK=0.46) compared to TCGA (βTCGA=0.65, p=0.001) and Komen (βKomen=0.80, p&lt;.0001).

Conclusion: Consistent with previous studies, we found that EAA in tumor samples varied across tumor subtypes. We also found that HK BC patients' epigenetic age accelerated at a different rate compared to predominantly white TCGA and Komen women, suggesting a potential racial biological difference. Large studies in other Asian populations are warranted to confirm our findings, which may provide biological insight into racial heterogeneity of BC, especially with regard to age at onset.

Citation Format: Hela Koka, Bin Zhu, Shelly Lap Ah Tse, Difei Wang, Maeve Kiely, Jennifer Lyn Guida, Priscilla Lee, Feng Wang, Cherry Wu, Koon Ho Tsang, Wing-cheong Chan, Sze Hong Law, Eric Karlins, Bin Zhu, Amy Hutchinson, Belynda Hicks, Xiaohong R. Yang. DNA methylation age of paired tumor-normal breast tissue in Chinese women with breast cancer [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 146.

Genetic Consequences of the Transatlantic Slave Trade in the Americas

According to historical records of transatlantic slavery, traders forcibly deported an estimated 12.5 million people from ports along the Atlantic coastline of Africa between the 16th and 19th centuries, with global impacts reaching to the present day, more than a century and a half after slavery's abolition. Such records have fueled a broad understanding of the forced migration from Africa to the Americas yet remain underexplored in concert with genetic data. Here, we analyzed genotype array data from 50,281 research participants, which-combined with historical shipping documents-illustrate that the current genetic landscape of the Americas is largely concordant with expectations derived from documentation of slave voyages. For instance, genetic connections between people in slave trading regions of Africa and disembarkation regions of the Americas generally mirror the proportion of individuals forcibly moved between those regions. While some discordances can be explained by additional records of deportations within the Americas, other discordances yield insights into variable survival rates and timing of arrival of enslaved people from specific regions of Africa. Furthermore, the greater contribution of African women to the gene pool compared to African men varies across the Americas, consistent with literature documenting regional differences in slavery practices. This investigation of the transatlantic slave trade, which is broad in scope in terms of both datasets and analyses, establishes genetic links between individuals in the Americas and populations across Atlantic Africa, yielding a more comprehensive understanding of the African roots of peoples of the Americas.

SomaticCombiner: improving the performance of somatic variant calling based on evaluation tests and a consensus approach

It is challenging to identify somatic variants from high-throughput sequence reads due to tumor heterogeneity, sub-clonality, and sequencing artifacts. In this study, we evaluated the performance of eight primary somatic variant callers and multiple ensemble methods using both real and synthetic whole-genome sequencing, whole-exome sequencing, and deep targeted sequencing datasets with the NA12878 cell line. The test results showed that a simple consensus approach can significantly improve performance even with a limited number of callers and is more robust and stable than machine learning based ensemble approaches. To fully exploit the multi-callers, we also developed a software package, SomaticCombiner, that can combine multiple callers and integrates a new variant allelic frequency (VAF) adaptive majority voting approach, which can maintain sensitive detection for variants with low VAFs.

Genetic and epigenetic intratumor heterogeneity impacts prognosis of lung adenocarcinoma

Intratumor heterogeneity (ITH) of genomic alterations may impact prognosis of lung adenocarcinoma (LUAD). Here, we investigate ITH of somatic copy number alterations (SCNAs), DNA methylation, and point mutations in lung cancer driver genes in 292 tumor samples from 84 patients with LUAD. LUAD samples show substantial SCNA and methylation ITH, and clonal architecture analyses present congruent evolutionary trajectories for SCNAs and DNA methylation aberrations. Methylation ITH mapping to gene promoter areas or tumor suppressor genes is low. Moreover, ITH composed of genetic and epigenetic mechanisms altering the same cancer driver genes is shown in several tumors. To quantify ITH for valid statistical association analyses, we develope an average pairwise ITH index (APITH), which does not depend on the number of samples per tumor. Both APITH indexes for SCNAs and methylation aberrations show significant associations with poor prognosis. This study further establishes the important clinical implications of genetic and epigenetic ITH in LUAD.

Many tumors are known to be heterogeneous. Here, the authors examined multiple samples from 84 patients with lung adenocarcinoma and demonstrate that the intratumor heterogeneity of methylation and copy number associates with poor prognosis.

Population Frequency of Fanconi Pathway Gene Variants and Their Association with Survival After Hematopoietic Cell Transplantation for Severe Aplastic Anemia

Severe aplastic anemia (SAA) is most frequently immune-mediated; however, rare inherited bone marrow failure syndromes, such as Fanconi anemia (FA), may be causal and can present as aplastic anemia (AA). FA is primarily an autosomal recessive disorder caused by the presence of 2 pathogenic variants in a single FA/BRCA DNA repair pathway gene. Patients with SAA often undergo genetic testing during clinical evaluation that may identify single deleterious alleles in FA pathway genes. We quantified the rate of germline single deleterious alleles in 22 FA genes using both a general population database (3234 variants, 125,748 exomes) and in a cohort of patients with SAA undergoing hematopoietic cell transplantation (HCT) (21 variants in 730 patients). The variants were classified as deleterious using in silico tools (REVEL, MetaSVM, VEP) and database resources (ClinVar, LOVD-FA). We found similar rates of single deleterious alleles in FA genes in both groups (2.6% and 2.9%). The presence of a single deleterious variant in a gene for FA in SAA HCT recipients did not affect the overall survival after HCT (hazard ratio, 0.85; 95% CI, 0.37 to 1.95; P  = 0.71), or post-HCT cancer risk (P = 0.52). Our results demonstrate that the identification of a germline monoallelic deleterious variant in an FA gene in patients with idiopathic SAA does not influence the outcome of HCT. Our findings suggest that there is no need for special treatment considerations for patients with SAA and a single deleterious FA allele identified on genetic testing.

Genome-wide association meta-analyses combining multiple risk phenotypes provide insights into the genetic architecture of cutaneous melanoma susceptibility

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

Genome-wide Association Study Identifies HLA-DPB1 as a Significant Risk Factor for Severe Aplastic Anemia

Severe aplastic anemia (SAA) is a rare disorder characterized by hypoplastic bone marrow and progressive pancytopenia. The etiology of acquired SAA is not understood but is likely related to abnormal immune responses and environmental exposures. We conducted a genome-wide association study of individuals with SAA genetically matched to healthy controls in discovery (359 cases, 1,396 controls) and validation sets (175 cases, 1,059 controls). Combined analyses identified linked SNPs in distinct blocks within the major histocompatibility complex on 6p21. The top SNP encodes p.Met76Val in the P4 binding pocket of the HLA class II gene HLA-DPB1 (rs1042151A>G, odds ratio [OR] 1.75, 95% confidence interval [CI] 1.50-2.03, p = 1.94 × 10-13) and was associated with HLA-DP cell surface expression in healthy individuals (p = 2.04 × 10-6). Phylogenetic analyses indicate that Val76 is not monophyletic and likely occurs in conjunction with different HLA-DP binding groove conformations. Imputation of HLA-DPB1 alleles revealed increased risk of SAA associated with Val76-encoding alleles DPB1∗03:01, (OR 1.66, p = 1.52 × 10-7), DPB1∗10:01 (OR 2.12, p = 0.0003), and DPB1∗01:01 (OR 1.60, p = 0.0008). A second SNP near HLA-B, rs28367832G>A, reached genome-wide significance (OR 1.49, 95% CI 1.22-1.78, p = 7.27 × 10-9) in combined analyses; the association remained significant after excluding cases with clonal copy-neutral loss-of-heterozygosity affecting class I HLA genes (8.6% of cases and 0% of controls). SNPs in the HLA class II gene HLA-DPB1 and possibly class I (HLA-B) are associated with SAA. The replacement of Met76 to Val76 in certain HLA-DPB1 alleles might influence risk of SAA through mechanisms involving DP peptide binding specificity, expression, and/or other factors affecting DP function.

Oral microbial community composition is associated with pancreatic cancer: A case-control study in Iran

BACKGROUND

Oral microbiota may be related to pancreatic cancer risk because periodontal disease, a condition linked to multiple specific microbes, has been associated with increased risk of pancreatic cancer. We evaluated the association between oral microbiota and pancreatic cancer in Iran.

METHODS

A total of 273 pancreatic adenocarcinoma cases and 285 controls recruited from tertiary hospitals and a specialty clinic in Tehran, Iran provided saliva samples and filled out a questionnaire regarding demographics and lifestyle characteristics. DNA was extracted from saliva and the V4 region of the 16S rRNA gene was PCR amplified and sequenced on the MiSeq. The sequencing data were processed using the DADA2 plugin in QIIME 2 and taxonomy was assigned against the Human Oral Microbiome Database. Logistic regression and MiRKAT models were calculated with adjustment for potential confounders.

RESULTS

No association was observed for alpha diversity with an average of 91.11 (standard deviation [SD] 2.59) sequence variants for cases and 89.42 (SD 2.58) for controls. However, there was evidence for an association between beta diversity and case status. The association between the Bray-Curtis dissimilarity and pancreatic cancer was particularly strong with a MiRKAT P-value of .000142 and specific principal coordinate vectors had strong associations with cancer risk. Several specific taxa were also associated with case status after adjustment for multiple comparisons.

CONCLUSION

The overall microbial community appeared to differ between pancreatic cancer cases and controls. Whether these reflect differences evident before development of pancreatic cancer will need to be evaluated in prospective studies.

Immune gene expression profiling reveals heterogeneity in luminal breast tumors

Background

Heterogeneity of immune gene expression patterns of luminal breast cancer (BC), which is clinically heterogeneous and overall considered as low immunogenic, has not been well studied especially in non-European populations. Here, we aimed at characterizing the immune gene expression profile of luminal BC in an Asian population and associating it with patient characteristics and tumor genomic features.

Methods

We performed immune gene expression profiling of tumor and adjacent normal tissue in 92 luminal BC patients from Hong Kong using RNA-sequencing data and used unsupervised consensus clustering to stratify tumors. We then used luminal patients from The Cancer Genome Atlas (TCGA, N = 564) and a Korean breast cancer study (KBC, N = 112) as replication datasets.

Results

Based on the expression of 130 immune-related genes, luminal tumors were stratified into three distinct immune subtypes. Tumors in one subtype showed higher level of tumor-infiltrating lymphocytes (TILs), characterized by T cell gene activation, higher expression of immune checkpoint genes, higher nonsynonymous mutation burden, and higher APOBEC-signature mutations, compared with other luminal tumors. The high-TIL subtype was also associated with lower ESR1/ESR2 expression ratio and increasing body mass index. The comparison of the immune profile in tumor and matched normal tissue suggested a tumor-derived activation of specific immune responses, which was only seen in high-TIL patients. Tumors in a second subtype were characterized by increased expression of interferon-stimulated genes and enrichment for TP53 somatic mutations. The presence of three immune subtypes within luminal BC was replicated in TCGA and KBC, although the pattern was more similar in Asian populations. The germline APOBEC3B deletion polymorphism, which is prevalent in East Asian populations and was previously linked to immune activation, was not associated with immune subtypes in our study. This result does not support the hypothesis that the germline APOBEC3B deletion polymorphism is the driving force for immune activation in breast tumors in Asian populations.

Conclusion

Our findings suggest that immune gene expression and associated genomic features could be useful to further stratify luminal BC beyond the current luminal A/B classification and a subset of luminal BC patients may benefit from checkpoint immunotherapy, at least in Asian populations.

Leukocyte telomere length in patients with myotonic dystrophy type I: a pilot study

Myotonic dystrophy type I (DM1) is an autosomal dominant disease of which clinical manifestations resemble premature aging. We evaluated the contribution of telomere length in pathogenesis in 361 DM1 patients (12 with serial measurements) and 223 unaffected relative controls using qPCR assay. While no differences in baseline leukocyte relative telomere length (RTL) was noted, the data suggested an accelerated RTL attrition in DM1 (discovery cohort: T/S change/year = −0.013 in DM1 vs. −0.005 in controls, P = 0.04); similar trend was noted in validation cohort. Further investigations are needed to examine the role of TL in the pathophysiology of DM1.

Survey on Scientific Shared Resource Rigor and Reproducibility

Shared scientific resources, also known as core facilities, support a significant portion of the research conducted at biomolecular research institutions. The Association of Biomolecular Resource Facilities (ABRF) established the Committee on Core Rigor and Reproducibility (CCoRRe) to further its mission of integrating advanced technologies, education, and communication in the operations of shared scientific resources in support of reproducible research. In order to first assess the needs of the scientific shared resource community, the CCoRRe solicited feedback from ABRF members via a survey. The purpose of the survey was to gain information on how U.S. National Institutes of Health (NIH) initiatives on advancing scientific rigor and reproducibility influenced current services and new technology development. In addition, the survey aimed to identify the challenges and opportunities related to implementation of new reporting requirements and to identify new practices and resources needed to ensure rigorous research. The results revealed a surprising unfamiliarity with the NIH guidelines. Many of the perceived challenges to the effective implementation of best practices (i.e., those designed to ensure rigor and reproducibility) were similarly noted as a challenge to effective provision of support services in a core setting. Further, most cores routinely use best practices and offer services that support rigor and reproducibility. These services include access to well-maintained instrumentation and training on experimental design and data analysis as well as data management. Feedback from this survey will enable the ABRF to build better educational resources and share critical best-practice guidelines. These resources will become important tools to the core community and the researchers they serve to impact rigor and transparency across the range of science and technology.

Pre-transplant short telomeres are associated with high mortality risk after unrelated donor haematopoietic cell transplant for severe aplastic anaemia

Telomeres are essential for chromosomal stability and markers of biological age. We evaluated the effect of pre-transplant short (<10th percentile-for-age) or very short (<5th or <1st percentile-for-age) leucocyte telomere length on survival after unrelated donor haematopoietic cell transplantation (HCT) for acquired severe aplastic anaemia (SAA). Patient pre-transplant blood samples and clinical data were available at the Center for International Blood and Marrow Transplant Research. We used quantitative real time polymerase chain reaction to measure relative telomere length (RTL) in 490 SAA patients who received HCT between 1990 and 2013 (median age = 20 years). One hundred and twelve patients (22·86%) had pre-HCT RTL <10th percentile-for-age, with the majority below the 5th percentile (N = 80, 71·43%). RTL <10th percentile-for-age was associated with a higher risk of post-HCT mortality (hazard ratio [HR] = 1·78, 95% confidence interval [CI]=1·18-2·69, P = 0·006) compared with RTL ≥50th percentile; no survival differences were noted in longer RTL categories (P > 0·10). Time-dependent effects for post-HCT mortality were only observed in relation to very short RTL; HR comparing RTL <5th versus ≥5th percentile = 1·38, P = 0·15 for the first 12 months after HCT, and HR = 3·91, P < 0·0001, thereafter, P-heterogeneity = 0·008; the corresponding HRs for RTL <1st versus ≥1st percentile = 1·29, P = 0·41, and HR = 5·18, P < 0·0001, P-heterogeneity = 0·005. The study suggests a potential role for telomere length in risk stratification of SAA patients in regard to their HCT survival.

Tuberculosis infection and lung adenocarcinoma: Mendelian randomization and pathway analysis of genome-wide association study data from never-smoking Asian women

We investigated whether genetic susceptibility to tuberculosis (TB) influences lung adenocarcinoma development among never-smokers using TB genome-wide association study (GWAS) results within the Female Lung Cancer Consortium in Asia. Pathway analysis with the adaptive rank truncated product method was used to assess the association between a TB-related gene-set and lung adenocarcinoma using GWAS data from 5512 lung adenocarcinoma cases and 6277 controls. The gene-set consisted of 31 genes containing known/suggestive associations with genetic variants from previous TB-GWAS. Subsequently, we followed-up with Mendelian Randomization to evaluate the association between TB and lung adenocarcinoma using three genome-wide significant variants from previous TB-GWAS in East Asians. The TB-related gene-set was associated with lung adenocarcinoma (p = 0.016). Additionally, the Mendelian Randomization showed an association between TB and lung adenocarcinoma (OR = 1.31, 95% CI: 1.03, 1.66, p = 0.027). Our findings support TB as a causal risk factor for lung cancer development among never-smoking Asian women.

Abstract 1638: Whole-exome sequencing and protein interaction networks to prioritize candidate genes for cutaneous melanoma susceptibility

Background: Known cutaneous malignant melanoma (CMM) genes account for melanoma risk in less than 40% of melanoma-prone families, suggesting the existence of additional risk genes or other modifiers. Whole exome sequencing (WES) of high-risk families usually results in the identification of a large number of potential disease-causing genes. However, prioritizing these genes and finding the true disease-causing ones have been challenging. We hypothesize that physically interacting proteins might have a similar effect on cell function and therefore, might produce a shared phenotype or increase the risk when mutated. To test this hypothesis, we applied a network approach integrating WES and protein-protein interaction (PPI) data to detect and prioritize rare variants in melanoma-prone families.

Methods: We conducted WES analysis on germline DNA of 212 patients from 45 CMM families (2-5 affected people in a family) without known mutations. Variants that were rare (&lt;1% in public and our in-house control datasets), loss-of-function or missense predicted to be deleterious, and showed complete disease co-segregation (present in all affected members within a family) were included in the network analysis. For this analysis, we only focused on genes that interact with known CMM genes (gene list in Goldstein AM, et al. Hum Mol Genet, 2017) which we used as seed proteins. We then applied a degree-aware algorithm (DADA) to rank the set of candidate genes from WES with respect to the set of known CMM genes using human PPI networks. The goal of DADA is to prioritize a candidate set of genes based on their association level with genes known to be related to the disease. Network proximity calculations of the candidates with respect to the seed set were based on Random Walk with Restarts method. We used high quality and scored human proteins interactions collected from InWeb_IM. In the primary analysis, the network contains 625,641 interactions, aggregated from 8 source databases and spanning 87% of reviewed human proteins. We also conducted additional analyses with CMM PPI subnetworks generated using high confidence interaction data from different sources to compare results.

Results: The WES analysis described above identified 546 genes that harbored rare variants, likely affect protein function and co-segregated in affected members in CMM families. After applying DADA across all networks, the rankings of these genes were fairly consistent. Top genes included both known CMM genes (such as MC1R, CDKN2B, KITLG, and TINF2) and genes that have not yet been associated with CMM susceptibility (such as FOXK1, AR, LEF1, RGL4, PABPC1, CHEK2, CTBP2, and BLM). These results highlight the importance of known CMM genes and their networks in CMM susceptibility, and demonstrate the potential value of network approaches in gene prioritization. Further evaluation of novel genes is in progress.

Citation Format: Sally Yepes, Margaret Tucker, Hela Koka, Kristine Jones, Aurelie Vogt, Laurie Burdette, Wen Luo, Bin Zhu, Meredith Yeager, Belynda Hicks, Neal D. Freedman, Stephen J. Chanock, Alisa M. Goldstein, Xiaohong R. Yang. Whole-exome sequencing and protein interaction networks to prioritize candidate genes for cutaneous melanoma susceptibility [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1638.

Abstract 4159: Gonadal mosaicism in a family with TP53-associated Li-Fraumeni syndrome

Li-Fraumeni syndrome (LFS) is an autosomal dominant cancer predisposition syndrome with exceptionally high lifetime cancer risks, caused primarily by germline pathogenic variants in TP53. Sarcomas, pre-menopausal breast cancer, brain tumors and adrenocortical carcinomas are the ‘core’ LFS cancers among a multitude of other cancers that can occur. Individuals with LFS are also at risk of developing multiple primary cancers over their lifetime. Clinical criteria have been established to denote LFS and ‘Li-Fraumeni like’ (LFL) families to facilitate genetic testing for TP53. De novo pathogenic germline variants in TP53 are estimated to cause between 7-25% of LFS. Gonadal mosaicism is an alternate explanation for presumed de novo variants and is yet to be described in LFS.

We present a proband who developed WHO grade II astrocytoma at age 18 years, who had a deceased brother with osteosarcoma diagnosed at age 13 years, and a healthy older sibling. The only other cancer history was leukemia at age 62 years in the paternal grandmother, who also had polycythemia vera. The family meets Eele’s criteria for LFL. The proband’s germline genetic testing revealed a heterozygous TP53 c.743G&gt;A p.R248Q pathogenic variant. Subsequent clinical genetic testing of the parents failed to reveal the variant, and the proband was presumed to have a de novo variant. Given the characteristic LFS-associated cancers in the brothers, we investigated further and identified the same TP53 pathogenic variant in the brother’s osteosarcoma tissue and pathology-confirmed benign liver tissue. TP53 sequencing of the brother’s tumor did not show loss of heterozygosity. Parental relationships to the proband and brother were confirmed through STR-based identifiler profiling. Ultra-deep sequencing (&gt;50,000x) of DNA derived from the parents’ blood, saliva and father’s benign colon polyp did not identify the variant, thus ruling out somatic mosaicism in the parents. We were unable to test gonadal tissues from the parents. In the setting of two siblings with the same pathogenic variant in TP53 which is absent in both parents, this family presents a case of assumed gonadal mosaicism resulting from a post-zygotic event.

The possibility of gonadal mosaicism presents complexity in genetic counselling of families with inherited disorders, and a conundrum in syndromes such as TP53-associated LFS where the presumed de novo rate is high and cancer risks are significantly elevated. The presented family suggests that further investigation of individuals with presumed de novo variation, especially in the setting of suspicious family cancer history, is required to accurately identify family members at-risk.

Citation Format: Payal P. Khincha, Kristine Jones, Kedest Teshome, Belynda Hicks, Phuong L. Mai, Sharon A. Savage. Gonadal mosaicism in a family with TP53-associated Li-Fraumeni syndrome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4159.

Subsequent neoplasm risk associated with rare variants in DNA repair and clinical radiation sensitivity syndrome genes: A report from the Childhood Cancer Survivor Study

10028

Background: Radiotherapy for childhood cancer is associated with strikingly elevated risk for developing subsequent neoplasms (SNs). Whether mutations in DNA repair and radiation sensitivity genes modulate SN risks is largely unknown. Methods: We conducted whole-exome sequencing in 5105 long-term childhood cancer survivors of European descent (mean follow-up = 32.7 years). SnpEff and ClinVar identified potentially damaging rare variants in 476 DNA repair or radiation sensitivity genes. Conditional logistic regression assessed SN risk associated with these variants aggregated by gene or pathway (N = 155 with ≥5 carriers). Controls were matched on sex, childhood cancer type and diagnosis age, radiation dose to the SN site, and survival. Exact p-values were calculated by permutation. Analyses used all survivors or subsets stratified on radiation dose. Results: A total of 1108 (21.7%) survivors developed at least one SN type known to be related to ionizing radiation exposure (e.g., breast cancer, basal cell carcinoma, meningioma, thyroid cancer, sarcoma). Radiation-related SN risk was associated with homologous recombination (HR) gene variants for SN sites that received > 0- < 10 Gy (20.9% cases, 11.0% controls; odds ratio [OR] = 2.20, 95% confidence interval [CI] 1.52-3.18; P = 1.41x10-4), most notably for FANCM (3.1% cases, 0.5% controls; OR = 9.91, 95%CI 3.73-26.4; P = 2.85x10-4). For radiation-related SNs at sites with higher doses (≥10 Gy), associations were not observed for the HR pathway (14.4% cases, 12.4% controls, P = 0.17) but were observed for two individual genes implicated in double-strand DNA break repair, EXO1 (1.8% cases, 0.4% controls; OR = 6.50, 95%CI 3.49-12.1; P = 7.43x10-4) and NEIL3 (0% cases, 1.0% controls; P = 3.23x10-4). Conclusions: In this discovery study, we identified dose-specific novel associations between SN risk after radiotherapy for childhood cancer and potentially damaging rare variants in genes involved in double-strand DNA break repair, particularly HR. If replicated, these results could impact long-term screening of childhood cancer survivors and risk-benefit assessments of treatment approaches.

SeqPlus sequencing methodology enables robust whole-genome sequencing, true variant detection, and novel genomic insights from archival esophageal carcinoma FFPE samples

e13016

Background: Whole-genome sequencing (WGS) of formalin-fixed, paraffin-embedded (FFPE) samples could enable novel insights from archival sample collections, yet robust FFPE WGS is challenged by fragmented DNA, uneven genomic coverage & sequencing artifacts attributed to FFPE fixation. We report our proprietary extraction & library preparation methodology (SeqPlus) with high quality, uniform WGS sequencing performance comparable to that from fresh-frozen samples. Methods: We analyzed 20 paired esophageal carcinoma (EC) samples i.e., primary tumors & matched germline samples to assess SeqPlus performance on 10-15-year-old FFPE tissues, measure variant concordance between WGS and a high-depth sequencing panel (269 genes, 400x coverage) & identify novel genomic features. Results: At a targeted 70x WGS tumor sequencing depth, 93% of the genome was covered by ³ 20 reads, 99% of bases had 10x coverage & average duplicate reads were 31%. We noted similar transition/transversion ratios & mutational spectra as from fresh-frozen EC specimens, suggesting that extraction & library preparation contributes to prior FFPE artifacts. Concordance of tumor-specific SNVs & indels derived from WGS & targeted panel was high at 86%. All 76 targeted panel-detected variants above the WGS limit of detection (mutant allele frequency [MAF] > 10%) were detected by WGS, 2 variants (2 tumors) were detected only by WGS, and 12 variants at MAF ≤ 6% (9 tumors) were only detected by the targeted panel. Tumor WGS yielded SNV, indels & CNV findings beyond variants detected by targeted sequencing. WGS enabled detection of 10.4 putative cancer variants per tumor compared to 12 variants per patient from frozen specimens and a median of 7 (up to 16) cancer-associated variants in genes outside the targeted panel. WGS copy number analysis revealed CCND1, EGFR, TP63, and SOX2amplification, CDKN2A/B deletion and additional unrecognized genomic aberrations. Conclusions: Our study reinforces the utility of high-quality, uniform WGS sequencing of archival FFPE cancer samples with SeqPlus and unlocks the potential for massive-scale retrospective genomic analysis of archived pathology samples with associated clinical & outcomes data.

Blood DNA methylation and breast cancer risk: a meta-analysis of four prospective cohort studies

BACKGROUND

Environmental and genetic factors play an important role in the etiology of breast cancer. Several small blood-based DNA methylation studies have reported risk associations with methylation at individual CpGs and average methylation levels; however, these findings require validation in larger prospective cohort studies. To investigate the role of blood DNA methylation on breast cancer risk, we conducted a meta-analysis of four prospective cohort studies, including a total of 1663 incident cases and 1885 controls, the largest study of blood DNA methylation and breast cancer risk to date.

METHODS

We assessed associations with methylation at 365,145 CpGs present in the HumanMethylation450 (HM450K) Beadchip, after excluding CpGs that did not pass quality controls in all studies. Each of the four cohorts estimated odds ratios (ORs) and 95% confidence intervals (CI) for the association between each individual CpG and breast cancer risk. In addition, each study assessed the association between average methylation measures and breast cancer risk, adjusted and unadjusted for cell-type composition. Study-specific ORs were combined using fixed-effect meta-analysis with inverse variance weights. Stratified analyses were conducted by age at diagnosis (< 50, ≥ 50), estrogen receptor (ER) status (+/-), and time since blood collection (< 5, 5-10, > 10 years). The false discovery rate (q value) was used to account for multiple testing.

RESULTS

The average age at blood draw ranged from 52.2 to 62.2 years across the four cohorts. Median follow-up time ranged from 6.6 to 8.4 years. The methylation measured at individual CpGs was not associated with breast cancer risk (q value > 0.59). In addition, higher average methylation level was not associated with risk of breast cancer (OR = 0.94, 95% CI = 0.85, 1.05; P = 0.26; P for study heterogeneity = 0.86). We found no evidence of modification of this association by age at diagnosis (P = 0.17), ER status (P = 0.88), time since blood collection (P = 0.98), or CpG location (P = 0.98).

CONCLUSIONS

Our data indicate that DNA methylation measured in the blood prior to breast cancer diagnosis in predominantly postmenopausal women is unlikely to be associated with substantial breast cancer risk on the HM450K array. Larger studies or with greater methylation coverage are needed to determine if associations exist between blood DNA methylation and breast cancer risk.

Germline mutations in Protection of Telomeres 1 in two families with Hodgkin lymphoma

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.

Prevalence of pathogenic/likely pathogenic variants in the 24 cancer genes of the ACMG Secondary Findings v2.0 list in a large cancer cohort and ethnicity-matched controls

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.