Genomic profiles and clinical presentation of chordoma.

Chordoma is a rare bone cancer with variable clinical outcomes. Here, we recruited 184 sporadic chordoma patients from the US and Canada and collected their clinical and treatment data. The average age at diagnosis was 45.5 years (Range 5-78) and the chordoma site distribution was 49.2% clivus, 26.2% spinal, and 24.0% sacral. Most patients (97.5%) received surgery as the primary treatment, among whom 85.3% also received additional treatment. Except for the most prevalent cancers like prostate, lung, breast, and skin cancer, there was no discernible enrichment for any specific cancer type among patients or their family members. Among a subset of patients (N = 70) with tumor materials, we conducted omics analyses and obtained targeted panel sequencing and SNP array genotyping data for 51 and 49 patients, respectively. The most recurrent somatic driver mutations included PIK3CA (12%), followed by chromatin remodeling genes PBRM1 and SETD2. Amplification of the 6q27 region, containing the chordoma susceptibility gene TBXT, was detected in eight patients (16.3%). Clival patients appeared to be less likely to carry driver gene mutations, chromosome arm level deletion events (e.g., 5p, 5p, and 9p), or 6q27 amplification compared to sacral patients. After adjusting for age, sex, tumor site, and additional treatment, patients with somatic deletions of 14q (OR = 13.73, 95% CI 1.96-96.02, P = 0.008) and 18p (OR = 13.68, 95% CI 1.77-105.89, P = 0.012) were more likely to have persistent chordoma. The study highlights genomic heterogeneity in chordoma, potentially linked to location and clinical progression.

Frequency of pathogenic germline variants in pediatric medulloblastoma survivors.

Background: Medulloblastoma is the most common malignant brain tumor in children. Most cases are sporadic, but well characterized germline alterations in APC, ELP1, GPR161, PTCH1, SUFU, and TP53 predispose to medulloblastoma. However, knowledge about pathogenic/likely pathogenic (P/LP) variants that predispose to medulloblastoma vary based on genes evaluated, patient demographics, and pathogenicity definitions. Methods: Germline exome sequencing was conducted on 160 childhood survivors of medulloblastoma. Analyses focused on rare variants in 239 known cancer susceptibility genes (CSGs). P/LP variants were identified using ClinVar and InterVar. Variants of unknown significance in known medulloblastoma predisposing genes (APC, ELP1, GPR161, PTCH1, SUFU, TP53) were further classified for loss of function variants. We compared the frequency of P/LP variants in cases to that in 1,259 cancer-free adult controls. Results: Twenty cases (12.5%) had a P/LP variant in an autosomal dominant CSG versus 5% in controls (p=1.0 x10-3), and 10 (6.3%) of these were P/LP variants in a known medulloblastoma gene, significantly greater than 0.2% observed in controls (p=1.4x10-8). The CSGs with the most P/LP variants in cases, and significantly higher than controls, were ELP1 (p=3.0x10-4) and SUFU (p=1.4x10-3). Conclusion: Approximately one in eight pediatric medulloblastoma survivors had an autosomal dominant P/LP CSG variant. We confirm several known associated genes and identify novel genes that may be important in medulloblastoma.

Genomic characterization of cervical lymph node metastases in papillary thyroid carcinoma following the Chornobyl accident.

Childhood radioactive iodine exposure from the Chornobyl accident increased papillary thyroid carcinoma (PTC) risk. While cervical lymph node metastases (cLNM) are well-recognized in pediatric PTC, the PTC metastatic process and potential radiation association are poorly understood. Here, we analyze cLNM occurrence among 428 PTC with genomic landscape analyses and known drivers (131I-exposed = 349, unexposed = 79; mean age = 27.9 years). We show that cLNM are more frequent in PTC with fusion (55%) versus mutation (30%) drivers, although the proportion varies by specific driver gene (RET-fusion = 71%, BRAF-mutation = 38%, RAS-mutation = 5%). cLNM frequency is not associated with other characteristics, including radiation dose. cLNM molecular profiling (N = 47) demonstrates 100% driver concordance with matched primary PTCs and highly concordant mutational spectra. Transcriptome analysis reveals 17 differentially expressed genes, particularly in the HOXC cluster and BRINP3; the strongest differentially expressed microRNA also is near HOXC10. Our findings underscore the critical role of driver alterations and provide promising candidates for elucidating the biological underpinnings of PTC cLNM.

Hemophagocytic Lymphohistiocytosis Gene Variants in Severe Aplastic Anemia and Their Impact on Hematopoietic Cell Transplantation Outcomes.

Germline genetic testing for patients with severe aplastic anemia (SAA) is recommended to guide treatment, including the use of immunosuppressive therapy and/or adjustment of hematopoietic cell transplantation (HCT) modalities. Hemophagocytic lymphohistiocytosis (HLH) is a life-threatening hyperinflammatory condition often associated with cytopenias with autosomal recessive (AR) or X-linked recessive (XLR) inheritance. HLH is part of the SAA differential diagnosis, and genetic testing may identify variants in HLH genes in patients with SAA. The impact of pathogenic/likely pathogenic (P/LP) variants in HLH genes on HCT outcomes in SAA is unclear. In this study, we aimed to determine the frequency of HLH gene variants in a large cohort of patients with acquired SAA and to evaluate their association(s) with HCT outcomes. The Transplant Outcomes in Aplastic Anemia project, a collaboration between the National Cancer Institute and the Center for International Blood and Marrow Transplant Research, collected genomic and clinical data from 824 patients who underwent HCT for SAA between 1989 and 2015. We excluded 140 patients with inherited bone marrow failure syndromes and used exome sequencing data from the remaining 684 patients with acquired SAA to identify P/LP variants in 14 HLH-associated genes (11 AR, 3 XLR) curated using American College of Medical Genetics and Genomics/Association of Molecular Pathology (ACMG/AMP) criteria. Deleterious variants of uncertain significance (del-VUS) were defined as those not meeting the ACMG/AMP P/LP criteria but with damaging predictions in ≥3 of 5 meta-predictors (BayesDel, REVEL, CADD, MetaSVM, and/or EIGEN). The Kaplan-Meier estimator was used to calculate the probability of overall survival (OS) after HCT, and the cumulative incidence calculator was used for other HCT outcomes, accounting for relevant competing risks. There were 46 HLH variants in 49 of the 684 patients (7.2%). Seventeen variants in 19 patients (2.8%) were P/LP; 8 of these were loss-of-function variants. Among the 19 patients with P/LP HLH variants, 16 (84%) had monoallelic variants in genes with AR inheritance, and 3 had variants in XLR genes. PRF1 was the most frequently affected gene (in 8 of the 19 patients). We found no statistically significant differences in transplantation-related factors between patients with and those without P/LP HLH variants. The 5-year survival probability was 89% (95% confidence interval [CI], 72% to 99%) in patients with P/LP HLH variants and 70% (95% CI, 53% to 85%) in those with del-VUS HLH variants, compared to 66% (95% CI, 62% to 70%) in those without variants (P = .16, log-rank test). The median time to neutrophil engraftment was 16 days for patients with P/LP HLH variants and 18 days in those with del-VUS HLH variants or without variants combined (P = .01, Gray's test). No statistically significant associations between P/LP HLH variants and the risk of acute or chronic graft-versus-host disease were noted. In this large cohort of patients with acquired SAA, we found that 2.8% of patients harbored a P/LP variant in an HLH gene. No negative effects of HLH gene variants on post-HCT survival were noted. The small number of patients with P/LP HLH variants limits the study's ability to provide conclusive evidence; nonetheless, our data suggest that there is no need for special transplantation considerations for patients with SAA carrying P/LP variants.

Multi-ancestry genome-wide association study of kidney cancer identifies 63 susceptibility regions.

Here, in a multi-ancestry genome-wide association study meta-analysis of kidney cancer (29,020 cases and 835,670 controls), we identified 63 susceptibility regions (50 novel) containing 108 independent risk loci. In analyses stratified by subtype, 52 regions (78 loci) were associated with clear cell renal cell carcinoma (RCC) and 6 regions (7 loci) with papillary RCC. Notably, we report a variant common in African ancestry individuals ( rs7629500 ) in the 3' untranslated region of VHL, nearly tripling clear cell RCC risk (odds ratio 2.72, 95% confidence interval 2.23-3.30). In cis-expression quantitative trait locus analyses, 48 variants from 34 regions point toward 83 candidate genes. Enrichment of hypoxia-inducible factor-binding sites underscores the importance of hypoxia-related mechanisms in kidney cancer. Our results advance understanding of the genetic architecture of kidney cancer, provide clues for functional investigation and enable generation of a validated polygenic risk score with an estimated area under the curve of 0.65 (0.74 including risk factors) among European ancestry individuals.

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.

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.

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.

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.

Associations of Circulating Estrogens and Estrogen Metabolites with Fecal and Oral Microbiome in Postmenopausal Women in the Ghana Breast Health Study.

The human fecal and oral microbiome may play a role in the etiology of breast cancer through modulation of endogenous estrogen metabolism. This study aimed to investigate associations of circulating estrogens and estrogen metabolites with the fecal and oral microbiome in postmenopausal African women. A total of 117 women with fecal (N = 110) and oral (N = 114) microbiome data measured by 16S rRNA gene sequencing, and estrogens and estrogen metabolites data measured by liquid chromatography tandem mass spectrometry were included. The outcomes were measures of the microbiome and the independent variables were the estrogens and estrogen metabolites. Estrogens and estrogen metabolites were associated with the fecal microbial Shannon index (global P < 0.01). In particular, higher levels of estrone (ß = 0.36, P = 0.03), 2-hydroxyestradiol (ß = 0.30, P = 0.02), 4-methoxyestrone (ß = 0.51, P = 0.01), and estriol (ß = 0.36, P = 0.04) were associated with higher levels of the Shannon index, while 16alpha-hydroxyestrone (ß = -0.57, P < 0.01) was inversely associated with the Shannon index as indicated by linear regression. Conjugated 2-methoxyestrone was associated with oral microbial unweighted UniFrac as indicated by MiRKAT (P < 0.01) and PERMANOVA, where conjugated 2-methoxyestrone explained 2.67% of the oral microbial variability, but no other estrogens or estrogen metabolites were associated with any other beta diversity measures. The presence and abundance of multiple fecal and oral genera, such as fecal genera from families Lachnospiraceae and Ruminococcaceae, were associated with several estrogens and estrogen metabolites as indicated by zero-inflated negative binomial regression. Overall, we found several associations of specific estrogens and estrogen metabolites and the fecal and oral microbiome. IMPORTANCE Several epidemiologic studies have found associations of urinary estrogens and estrogen metabolites with the fecal microbiome. However, urinary estrogen concentrations are not strongly correlated with serum estrogens, a known risk factor for breast cancer. To better understand whether the human fecal and oral microbiome were associated with breast cancer risk via the regulation of estrogen metabolism, we conducted this study to investigate the associations of circulating estrogens and estrogen metabolites with the fecal and oral microbiome in postmenopausal African women. We found several associations of parent estrogens and several estrogen metabolites with the microbial communities, and multiple individual associations of estrogens and estrogen metabolites with the presence and abundance of multiple fecal and oral genera, such as fecal genera from families Lachnospiraceae and Ruminococcaceae, which have estrogen metabolizing properties. Future large, longitudinal studies to investigate the dynamic changes of the fecal and oral microbiome and estrogen relationship are needed.

Evaluation of alcohol-free mouthwash for studies of the oral microbiome.

Oral bacteria play important roles in human health and disease. Oral samples collected using ethanol-containing mouthwash are widely used for oral microbiome studies. However, ethanol is flammable and not ideal for transportation/storage in large quantities, and some individuals may avoid ethanol due to the burning sensation or due to various personal, medical, religious, and/or cultural factors. Here, we compared ethanol-free and ethanol-containing mouthwashes using multiple microbiome metrics and assessed the stability of the mouthwash samples stored up to 10 days before processing. Forty volunteers provided oral wash samples collected using ethanol-free and ethanol-containing mouthwashes. From each sample, one aliquot was immediately frozen, one was stored at 4°C for 5 days and frozen, while the third aliquot was stored for 5 days at 4°C and 5 days at ambient temperature to mimic shipping delays and then frozen. DNA was extracted, the 16S rRNA gene V4 region was amplified and sequenced, and bioinformatic processing was performed using QIIME 2. Microbiome metrics measured in the two mouthwash types were very similar, with intraclass correlation coefficients (ICCs) for alpha and beta diversity metrics greater than 0.85. Relative abundances of some taxa were significantly different, but ICCs of the top four most abundant phyla and genera were high (> 0.75) for the comparability of the mouthwashes. Stability during delayed processing was also high for both mouthwashes based on alpha and beta diversity measures and relative abundances of the top four phyla and genera (ICCs ≥ 0.90). These results demonstrate ethanol-free mouthwash performs similarly to ethanol-containing mouthwash for microbial analyses, and both mouthwashes are stable for at least 10 days without freezing prior to laboratory processing. Ethanol-free mouthwash is suitable for collecting and shipping oral wash samples, and these results have important implications for planning future epidemiologic studies of the oral microbiome.

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.

Oral microbiome and risk of incident head and neck cancer: A nested case-control study.

OBJECTIVES: This nested case-control study in the NIH-AARP Diet and Health Study was carried out to prospectively investigate the relationship of oral microbiome with head and neck cancer (HNC). MATERIALS AND METHODS: 56 incident HNC cases were identified, and 112 controls were incidence-density matched to cases. DNA extracted from pre-diagnostic oral wash samples was whole-genome shotgun metagenomic sequenced to measure the overall oral microbiome. ITS2 gene qPCR was used to measure the presence of fungi. ITS2 gene sequencing was performed on ITS2 gene qPCR positive samples. We computed taxonomic and functional alpha-diversity and beta-diversity metrics. The presence and relative abundance of groups of red-complex (e.g., Porphyromonas gingivalis) and/or orange-complex (e.g., Fusobacterium nucleatum) periodontal pathogens were compared between cases and controls using conditional logistic regression models and MiRKAT. RESULTS: Participants with higher taxonomic microbial alpha-diversity had a non-statistically significant decreased risk of HNC. No case-control differences were found for beta diversity by MiRKAT model (all p > 0.05). A greater relative abundance of red-complex periodontal pathogens (OR = 0.51, 95 % CI = 0.26-1.00), orange-complex (OR = 0.38, 95 % CI = 0.18-0.83), and both complexes' pathogens (OR = 0.32, 95 % CI = 0.14-0.75), were associated with reduced risk of HNC. The presence of oral fungi was also strongly associated with reduced risk of HNC compared with controls (OR = 0.39, 95 % CI = 0.17-0.92). CONCLUSION: Greater taxonomic alpha-diversity, the presence of oral fungi, and the presence or relative abundance of multiple microbial species, including the red- and orange-complex periodontal pathogens, were associated with reduced risk of HNC. Future studies with larger sample sizes are needed to evaluate these associations.

GWAS Explorer: an open-source tool to explore, visualize, and access GWAS summary statistics in the PLCO Atlas.

The Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial is a prospective cohort study of nearly 155,000 U.S. volunteers aged 55-74 at enrollment in 1993-2001. We developed the PLCO Atlas Project, a large resource for multi-trait genome-wide association studies (GWAS), by genotyping participants with available DNA and genomic consent. Genotyping on high-density arrays and imputation was performed, and GWAS were conducted using a custom semi-automated pipeline. Association summary statistics were generated from a total of 110,562 participants of European, African and Asian ancestry. Application programming interfaces (APIs) and open-source software development kits (SKDs) enable exploring, visualizing and open data access through the PLCO Atlas GWAS Explorer website, promoting Findable, Accessible, Interoperable, and Re-usable (FAIR) principles. Currently the GWAS Explorer hosts association data for 90 traits and >78,000,000 genomic markers, focusing on cancer and cancer-related phenotypes. New traits will be posted as association data becomes available. The PLCO Atlas is a FAIR resource of high-quality genetic and phenotypic data with many potential reuse opportunities for cancer research and genetic epidemiology.

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.

Impact of Transcript (p16/p14ARF) Alteration on Cancer Risk in CDKN2A Germline Pathogenic Variant Carriers.

BACKGROUND: Few studies have evaluated the relationship between CDKN2A germline pathogenic variants (GPV), transcript (p16/p14ARF) alteration, and cancer risk. METHODS: Standardized incidence ratios (SIRs) comparing cancer risk with the general population were calculated for 385 CDKN2A GPV carriers from 2 large cohorts (259 United States and 126 Swedish individuals) using Poisson regression; statistical significance was defined as P less than .002 (Bonferroni correction). Cumulative incidence is reported for melanoma and nonmelanoma cancer. RESULTS: Incidence was increased for melanoma (SIR = 159.8, 95% confidence interval [CI] = 132.1 to 193.2), pancreatic cancer (SIR = 24.1, 95% CI = 14.7 to 39.4), head and neck squamous cell carcinoma (SIR = 16.2, 95% CI = 9.5 to 27.6), and lung cancer (SIR = 5.6, 95% CI = 3.4 to 9.1) in GPV carriers. Similar associations were observed with p16 alteration. Combined p16 and p14ARF alteration was associated with increased incidence of esophageal cancer (SIR = 16.7, 95% CI = 5.7 to 48.9) and malignant peripheral nerve sheath tumor (SIR = 113.0, 95% CI = 16.4 to 780.9), although cancer events were limited (n < 5 for each malignancy). Cumulative incidence at age 70 years for melanoma and nonmelanoma cancer was 68.3% (95% CI = 68.0% to 68.6%) and 35.2% (95% CI = 34.9% to 35.6%), respectively. A total 89% of smoking-related cancers (lung, head and neck squamous cell carcinoma, pancreatic, esophageal) occurred in ever smokers. CONCLUSION: These findings highlight the impact of p16 and p14ARF alteration on cancer risk. Smoking was an important risk factor for smoking-related cancers in our study.

Germline-somatic JAK2 interactions are associated with clonal expansion in myelofibrosis.

Myelofibrosis is a rare myeloproliferative neoplasm (MPN) with high risk for progression to acute myeloid leukemia. Our integrated genomic analysis of up to 933 myelofibrosis cases identifies 6 germline susceptibility loci, 4 of which overlap with previously identified MPN loci. Virtual karyotyping identifies high frequencies of mosaic chromosomal alterations (mCAs), with enrichment at myelofibrosis GWAS susceptibility loci and recurrently somatically mutated MPN genes (e.g., JAK2). We replicate prior MPN associations showing germline variation at the 9p24.1 risk haplotype confers elevated risk of acquiring JAK2V617F mutations, demonstrating with long-read sequencing that this relationship occurs in cis. We also describe recurrent 9p24.1 large mCAs that selectively retained JAK2V617F mutations. Germline variation associated with longer telomeres is associated with increased myelofibrosis risk. Myelofibrosis cases with high-frequency JAK2 mCAs have marked reductions in measured telomere length - suggesting a relationship between telomere biology and myelofibrosis clonal expansion. Our results advance understanding of the germline-somatic interaction at JAK2 and implicate mCAs involving JAK2 as strong promoters of clonal expansion of those mutated clones.

The Oral Microbiome and Lung Cancer Risk: An Analysis of 3 Prospective Cohort Studies.

BACKGROUND: Previous studies suggested associations between the oral microbiome and lung cancer, but studies were predominantly cross-sectional and underpowered. METHODS: Using a case-cohort design, 1306 incident lung cancer cases were identified in the Agricultural Health Study; National Institutes of Health-AARP Diet and Health Study; and Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Referent subcohorts were randomly selected by strata of age, sex, and smoking history. DNA was extracted from oral wash specimens using the DSP DNA Virus Pathogen kit, the 16S rRNA gene V4 region was amplified and sequenced, and bioinformatics were conducted using QIIME 2. Hazard ratios and 95% confidence intervals were calculated using weighted Cox proportional hazards models. RESULTS: Higher alpha diversity was associated with lower lung cancer risk (Shannon index hazard ratio = 0.90, 95% confidence interval = 0.84 to 0.96). Specific principal component vectors of the microbial communities were also statistically significantly associated with lung cancer risk. After multiple testing adjustment, greater relative abundance of 3 genera and presence of 1 genus were associated with greater lung cancer risk, whereas presence of 3 genera were associated with lower risk. For example, every SD increase in Streptococcus abundance was associated with 1.14 times the risk of lung cancer (95% confidence interval = 1.06 to 1.22). Associations were strongest among squamous cell carcinoma cases and former smokers. CONCLUSIONS: Multiple oral microbial measures were prospectively associated with lung cancer risk in 3 US cohort studies, with associations varying by smoking history and histologic subtype. The oral microbiome may offer new opportunities for lung cancer prevention.