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.

JAK2 V617F mutation and associated chromosomal alterations in primary and secondary myelofibrosis and post-HCT outcomes.

JAK2 V617F is the most common driver mutation in primary or secondary myelofibrosis for which allogeneic hematopoietic cell transplantation (HCT) is the only curative treatment. Knowledge of the prognostic utility of JAK2 alterations in the HCT setting is limited. We identified all patients with MF who received HCT between 2000 and 2016 and had a pre-HCT blood sample (N = 973) available at the Center of International Blood and Marrow Transplant Research biorepository. PacBio sequencing and single nucleotide polymorphism-array genotyping were used to identify JAK2V617F mutation and associated mosaic chromosomal alterations (mCAs), respectively. Cox proportional hazard models were used for HCT outcome analyses. Genomic testing was complete for 924 patients with MF (634 primary MF [PMF], 135 postpolycythemia vera [PPV-MF], and 155 postessential thrombocytopenia [PET-MF]). JAK2V617F affected 562 patients (57.6% of PMF, 97% of PPV-MF, and 42.6% of PET-MF). Almost all patients with mCAs involving the JAK2 region (97.9%) were JAK2V617-positive. In PMF, JAK2V617F mutation status, allele burden, or identified mCAs were not associated with disease progression/relapse, nonrelapse mortality (NRM), or overall survival. Almost all PPV-MF were JAK2V617F-positive (97%), with no association between HCT outcomes and mutation allele burden or mCAs. In PET-MF, JAK2V617F high mutation allele burden (≥60%) was associated with excess risk of NRM, restricted to transplants received in the era of JAK inhibitors (2013-2016; hazard ratio = 7.65; 95% confidence interval = 2.10-27.82; P = .002). However, allele burden was not associated with post-HCT disease progression/relapse or survival. Our findings support the concept that HCT can mitigate the known negative effect of JAK2V617F in patients with MF, particularly for PMF and PPV-MF.

Joint associations between established genetic susceptibility loci, pesticide exposures, and risk of prostate cancer.

More than 200 genetic variants have been independently associated with prostate cancer risk. Studies among farmers have also observed increased prostate cancer risk associated with exposure to specific organophosphate (fonofos, terbufos, malathion, dimethoate) and organochlorine (aldrin, chlordane) insecticides. We examined the joint associations between these pesticides, established prostate cancer loci, and prostate cancer risk among 1,162 cases (588 aggressive) and 2,206 frequency-matched controls nested in the Agricultural Health Study cohort. History of lifetime pesticide use was combined with a polygenic risk score (PRS) generated using 256 established prostate cancer risk variants. Logistic regression models estimated the joint associations of the pesticides, the PRS, and the 256 individual genetic variants with risk of total and aggressive prostate cancer. Likelihood ratio tests assessed multiplicative interaction. We observed interaction between ever use of fonofos and the PRS in relation to total and aggressive prostate cancer risk. Compared to the reference group (never use, PRS < median), men with ever use of fonofos and PRS > median had elevated risks of total (OR 1.35 [1.06-1.73], p-interaction = 0.03) and aggressive (OR 1.49 [1.09-2.04], p-interaction = 0.19) prostate cancer. There was also suggestion of interaction between pesticides and individual genetic variants occurring in regions associated with DNA damage response (CDH3, EMSY genes) and with variants related to altered androgen receptor-driven transcriptional programs critical for prostate cancer. Our study provides evidence that men with greater genetic susceptibility to prostate cancer may be at higher risk if they are also exposed to pesticides and suggests potential mechanisms by which pesticides may increase prostate cancer risk.

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.

Short-term changes in ultrasound tomography measures of breast density and treatment-associated endocrine symptoms after tamoxifen therapy.

Although breast density decline with tamoxifen therapy is associated with greater therapeutic benefit, limited data suggest that endocrine symptoms may also be associated with improved breast cancer outcomes. However, it is unknown whether endocrine symptoms are associated with reductions in breast density after tamoxifen initiation. We evaluated treatment-associated endocrine symptoms and breast density change among 74 women prescribed tamoxifen in a 12-month longitudinal study. Treatment-associated endocrine symptoms and sound speed measures of breast density, assessed via novel whole breast ultrasound tomography (m/s), were ascertained before tamoxifen (T0) and at 1-3 (T1), 4-6 (T2), and 12 months (T3) after initiation. CYP2D6 status was genotyped, and tamoxifen metabolites were measured at T3. Using multivariable linear regression, we estimated mean change in breast density by treatment-associated endocrine symptoms adjusting for age, race, menopausal status, body mass index, and baseline density. Significant breast density declines were observed in women with treatment-associated endocrine symptoms (mean change (95% confidence interval) at T1:-0.26 m/s (-2.17,1.65); T2:-2.12 m/s (-4.02,-0.22); T3:-3.73 m/s (-5.82,-1.63); p-trend = 0.004), but not among women without symptoms (p-trend = 0.18) (p-interaction = 0.02). Similar declines were observed with increasing symptom frequency (p-trends for no symptoms = 0.91; low/moderate symptoms = 0.03; high symptoms = 0.004). Density declines remained among women with detectable tamoxifen metabolites or intermediate/efficient CYP2D6 metabolizer status. Emergent/worsening endocrine symptoms are associated with significant, early declines in breast density after tamoxifen initiation. Further studies are needed to assess whether these observations predict clinical outcomes. If confirmed, endocrine symptoms may be a proxy for tamoxifen response and useful for patients and providers to encourage adherence.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding.

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

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.

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.

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.

Massively parallel reporter assays and variant scoring identified functional variants and target genes for melanoma loci and highlighted cell-type specificity.

The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRAs) by using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRAs, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (eQTLs or meQTLs, respectively) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription-factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRAs along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional.

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.

Next-generation sequencing errors due to genetic variation in WRAP53 encoding TCAB1 on chromosome 17.

Next-generation sequencing (NGS) is a valuable tool, but has limitations in sequencing through repetitive runs of single nucleotides (homopolymers). Pathogenic germline variants in WRAP53 encoding telomere Cajal body protein 1 (TCAB1) are a known cause of dyskeratosis congenita. We identified a significant NGS error in WRAP53, c.1562dup, p.Ala522Glyfs*8 (rs755116516 G>-/GG/GGG) that did not validate by Sanger sequencing. This error occurs because rs755116516 G>-/GG/GGG (Chr17:7,606,714) is polymorphic, and variants at this site challenge the ability of NGS to accurately call the correct number of nucleotides in a homopolymer run. This was further complicated by the fact that chr17:7,606,721 (rs769202794) is multiallelic G>A, C, T, and that chr17:7,606,722 is also multiallelic (rs7640C>A/G/T and rs373064567C>delC). In addition to the expert interpretation of potentially clinically actionable variants, it recommended that all variants in regions of the genome with homopolymers be validated by Sanger sequencing before clinical action.

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.

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.

The oral microbiome and breast cancer and nonmalignant breast disease, and its relationship with the fecal microbiome in the Ghana Breast Health Study.

The oral microbiome, like the fecal microbiome, may be related to breast cancer risk. Therefore, we investigated whether the oral microbiome was associated with breast cancer and nonmalignant breast disease, and its relationship with the fecal microbiome in a case-control study in Ghana. A total of 881 women were included (369 breast cancers, 93 nonmalignant cases and 419 population-based controls). The V4 region of the 16S rRNA gene was sequenced from oral and fecal samples. Alpha-diversity (observed amplicon sequence variants [ASVs], Shannon index and Faith's Phylogenetic Diversity) and beta-diversity (Bray-Curtis, Jaccard and weighted and unweighted UniFrac) metrics were computed. MiRKAT and logistic regression models were used to investigate the case-control associations. Oral sample alpha-diversity was inversely associated with breast cancer and nonmalignant breast disease with odds ratios (95% CIs) per every 10 observed ASVs of 0.86 (0.83-0.89) and 0.79 (0.73-0.85), respectively, compared to controls. Beta-diversity was also associated with breast cancer and nonmalignant breast disease compared to controls (P ≤ .001). The relative abundances of Porphyromonas and Fusobacterium were lower for breast cancer cases compared to controls. Alpha-diversity and presence/relative abundance of specific genera from the oral and fecal microbiome were strongly correlated among breast cancer cases, but weakly correlated among controls. Particularly, the relative abundance of oral Porphyromonas was strongly, inversely correlated with fecal Bacteroides among breast cancer cases (r = -.37, P ≤ .001). Many oral microbial metrics were strongly associated with breast cancer and nonmalignant breast disease, and strongly correlated with fecal microbiome among breast cancer cases, but not controls.

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.