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.

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.

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.

Genome-wide association study of café-au-lait macule number in neurofibromatosis type 1.

BACKGROUND: Neurofibromatosis type 1 (NF1) is a tumor-predisposition disorder that arises due to pathogenic variants in tumor suppressor NF1. NF1 has variable expressivity that may be due, at least in part, from heritable elements such as modifier genes; however, few genetic modifiers have been identified to date. METHODS: In this study, we performed a genome-wide association analysis of the number of café-au-lait macules (CALM) that are considered a tumor-like trait as a clinical phenotype modifying NF1. RESULTS: A borderline genome-wide significant association was identified in the discovery cohort (CALM1, N = 112) between CALM number and rs12190451 (and rs3799603, r2  = 1.0; p = 7.4 × 10-8 ) in the intronic region of RPS6KA2. Although, this association was not replicated in the second cohort (CALM2, N = 59) and a meta-analysis did not show significantly associated variants in this region, a significant corroboration score (0.72) was obtained for the RPS6KA2 signal in the discovery cohort (CALM1) using Complementary Pairs Stability Selection for Genome-Wide Association Studies (ComPaSS-GWAS) analysis, suggesting that the lack of replication may be due to heterogeneity of the cohorts rather than type I error. CONCLUSION: rs12190451 is located in a melanocyte-specific enhancer and may influence RPS6KA2 expression in melanocytes-warranting further functional studies.

Analysis of TP53 aflatoxin signature mutation in hepatocellular carcinomas from Guatemala: A cross-sectional study (2016-2017).

BACKGROUND AND AIMS: Guatemala has the highest incidence of hepatocellular carcinoma (HCC) in the Western hemisphere. The major risk factors in Guatemala are not well characterized, but the prevalence of hepatitis B virus (HBV) and hepatitis C virus (HCV) appears to be low, while the prevalence of aflatoxin (AFB1) exposure appears to be high. To examine whether AFB1 may contribute to the elevated incidence of HCC in Guatemala, this study examined the frequency of the AFB1-signature mutation in the TP53 gene (R249S) as well as other somatic mutations. In addition, we assessed whether the frequency of the TP53 mutation differed by sex. METHODS: Formalin-fixed, paraffin-embedded (FFPE) HCC tissues were obtained from three hospitals in Guatemala City between 2016 and 2017. In addition, tumor tissues preserved in RNAlater were also obtained. Sociodemographic and clinical information including HBV and HCV status were collected. Targeted sequencing of TP53 was performed in the FFPE samples, and a panel of 253 cancer-related genes was sequenced in the RNAlater samples. RESULTS: Ninety-one FFPE tissues were examined, from 52 men and 39 women. Median (IQR) age at diagnosis was 62 (51-70). Among those with known HBV and HCV status, two were HBV+ and three were HCV+. Overall, 47% of the HCCs had a TP53 mutation. The AFB1-signature R249S mutation was present in 24%. No overlap between the R249S mutation and HBV+ was observed in this cohort. Among 18 RNAlater samples examined, 44% had any TP53 mutation and 33% had the R249S mutation. Other somatic mutations were identified in known HCC driver genes. CONCLUSIONS: The presence of the TP53 R249S mutation in the samples studied suggests that AFB1 may contribute to the high incidence of HCC in Guatemala. The proportion of HBV+ tumors was low, suggesting that AFB1 may be associated with HCC in the absence of concomitant HBV infection. Further investigation of AFB1 and other risk factors for HCC in Guatemala is warranted.

Successful use of whole genome amplified DNA from multiple source types for high-density Illumina SNP microarrays.

BACKGROUND: The recommended genomic DNA input requirements for whole genome single nucleotide polymorphism microarrays can limit the scope of molecular epidemiological studies. We performed a large-scale evaluation of whole genome amplified DNA as input into high-density, whole-genome Illumina® Infinium® SNP microarray. RESULTS: Overall, 6622 DNA samples from 5970 individuals were obtained from three distinct biospecimen sources and genotyped using gDNA and/or wgaDNA inputs. When genotypes from the same individual were compared with standard, native gDNA input amount, we observed 99.94% mean concordance with wgaDNA input. CONCLUSIONS: Our results demonstrate that carefully conducted studies with wgaDNA inputs can yield high-quality genotyping results. These findings should enable investigators to consider expansion of ongoing studies using high-density SNP microarrays, currently challenged by small amounts of available DNA.

Effect of pre-analytic variables on the reproducibility of qPCR relative telomere length measurement.

Telomeres, long nucleotide repeats and a protein complex at chromosome ends, shorten with each cell division and are susceptible to oxidative damage. Quantitative PCR (qPCR) is a widely-used technique to measure relative telomere length (RTL) in DNA samples but is challenging to optimize and significant lab-to-lab variability has been reported. In this study, we evaluated factors that may contribute to qPCR RTL measurement variability including DNA extraction methods, methods used for removing potential residual PCR inhibitors, sample storage conditions, and sample location in the PCR plate. Our results show that the DNA extraction and purification techniques, as well as sample storage conditions introduce significant variability in qPCR RTL results. We did not find significant differences in results based on sample location in the PCR plate or qPCR instrument used. These data suggest that lack of reproducibility in published association studies of RTL could be, in part, due to methodological inconsistencies. This study illustrates the importance of uniform sample handling, from DNA extraction through data generation and analysis, in using qPCR to determine RTL.

A comprehensive candidate gene approach identifies genetic variation associated with osteosarcoma.

BACKGROUND: Osteosarcoma (OS) is a bone malignancy which occurs primarily in adolescents. Since it occurs during a period of rapid growth, genes important in bone formation and growth are plausible modifiers of risk. Genes involved in DNA repair and ribosomal function may contribute to OS pathogenesis, because they maintain the integrity of critical cellular processes. We evaluated these hypotheses in an OS association study of genes from growth/hormone, bone formation, DNA repair, and ribosomal pathways. METHODS: We evaluated 4836 tag-SNPs across 255 candidate genes in 96 OS cases and 1426 controls. Logistic regression models were used to estimate the odds ratios (OR) and 95% confidence intervals (CI). RESULTS: Twelve SNPs in growth or DNA repair genes were significantly associated with OS after Bonferroni correction. Four SNPs in the DNA repair gene FANCM (ORs 1.9-2.0, P = 0.003-0.004) and 2 SNPs downstream of the growth hormone gene GH1 (OR 1.6, P = 0.002; OR 0.5, P = 0.0009) were significantly associated with OS. One SNP in the region of each of the following genes was significant: MDM2, MPG, FGF2, FGFR3, GNRH2, and IGF1. CONCLUSIONS: Our results suggest that several SNPs in biologically plausible pathways are associated with OS. Larger studies are required to confirm our findings.

Telomere length and variation in telomere biology genes in individuals with osteosarcoma.

Osteosarcoma, the most common primary bone tumor, occurs most frequently in adolescents. Chromosomal aneuploidy is common in osteosarcoma cells, suggesting underlying chromosomal instability. Telomeres, located at chromosome ends, are essential for genomic stability; several studies have suggested that germline telomere length (TL) is associated with cancer risk. We hypothesized that TL and/or common genetic variation in telomere biology genes may be associated with risk of osteosarcoma. We investigated TL in peripheral blood DNA and 713 single nucleotide polymorphisms (SNPs) from 39 telomere biology genes in 98 osteosarcoma cases and 69 orthopedic controls. For the genotyping component, we added 1363 controls from the Prostate, Lung, Colorectal, and Ovarian Cancer ScreeningTrial. Short TL was not associated with osteosarcoma risk overall (OR 1.39, P=0.67), although there was a statistically significant association in females (OR 4.35, 95% Cl 1.20-15.74, P=0.03). Genotype analyses identified seven SNPs in TERF1 significantly associated with osteosarcoma risk after Bonferroni correction by gene. These SNPs were highly linked and associated with a reduced risk of osteosarcoma (OR 0.48-0.53, P=0.0001-0.0006). We also investigated associations between TL and telomere gene SNPs in osteosarcoma cases and orthopedic controls. Several SNPs were associated with TL prior to Bonferroni correction; one SNP in NOLA2 and one in MEN1 were marginally non-significant after correction (P(adj)=0.057 and 0.066, respectively). This pilot-study suggests that females with short telomeres may be at increased risk of osteosarcoma, and that SNPs in TERF1 are inversely associated with osteosarcoma risk.

Genetic variation at chromosome 8q24 in osteosarcoma cases and controls.

Osteosarcoma is a primary bone malignancy that typically occurs during the pubertal growth spurt. Only a few small association studies have evaluated common germ line variation in individuals with osteosarcoma. The 8q24 chromosomal region contains several loci that are associated with risk of many different cancers. We conducted an association study of common single-nucleotide polymorphisms (SNPs) across 8q24 to explore the role this region may play in osteosarcoma risk. We genotyped 214 tag SNPs in 99 osteosarcoma cases and 1430 controls (65 controls from a hospital-based case-control study and 1365 controls from a population-based study). Additive, dominant and recessive genetic models were evaluated using unconditional logistic regression to estimate odds ratios (ORs) and 95% confidence intervals (CIs). Analyses of nine SNPs previously associated with cancer did not show strong statistically significant associations. Of the remaining 205 SNPs, 7 were statistically significant (P