Supplementary Online Content Hayes HB, Ahn J, Fan X, et al. Association of Oral Microbiome With Risk for Incident Head and Neck Squamous Cell Cancer [published online January 11, 2018]. JAMA Oncol. doi:10.1001/jamaoncol.2017.4777 eTable 1. Number of filtered sequence reads per sample eTable 2. Sequence reads per bacterial phylum, by cohort eTable 3a. Median counts and fold changes (FC) for the association1 between selected bacteria related to periodontal diseases and dental caries and risk of head and neck squamous cell cancer (HNSCC) in the 2 cohorts eTable 3b. Carriage rate and odds ratios (OR) for the association between selected bacteria related to periodontal disease and dental caries and risk of head and neck squamous cell cancer (HNSCC) in the 2 cohorts eTable 4. Median counts and fold changes (FC) for the association between selected taxa and risk of head and neck squamous cell cancer (HNSCC) according to smoking status in the 2 cohorts eTable 5. Median counts and fold changes (FC) for the association between selected taxa and risk of head and neck squamous cell cancer (HNSCC) according to time from oral wash sample collection to diagnosis in the 2 cohorts Table 6. Median counts and fold changes (FC) for the association between selected taxa and risk of head and neck squamous cell cancer (HNSCC) according to site of cancer (oral cavity, pharynx, and larynx) in the 2 cohorts eFigure. Overall oral microbiota structure according to subsequent occurrence of squamous-cell head and neck cancer (SCHNC) and SCHNC of the oral cavity, pharynx, and larynx in a case–control study nested within 2 cohorts This supplementary material has been provided by the authors to give readers additional information about their work. 1 1 Supplementary Table 1. Number of filtered sequence reads per sample Total Cases Controls (N=383) (N=129) (N=254) Cohort Mean ± SD Mean ± SD Mean ± SD CPS-II 9,145 ± 2,992 9,004 ± 2,767 9,217 ± 3,110 PLCO 7,744 ± 2,830 7,987 ± 3,460 7,620 ± 2,455 Total 8,373 ± 2,983 8,445 ± 3,196 8,337 ± 2,875 2 1. Poor-quality sequences were excluded using the default parameters of the QIIME script split_libraries.py 3 (minimum average quality score=25, minimum/maximum sequence length=200/1000 base pairs, no ambiguous 4 base calls, and no mismatches allowed in the primer sequence). 5 © 2017 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 6 Supplementary Table 2. Sequence reads per bacterial phylum, by cohort CPS-II Cohort PLCO Cohort All Subjects (N=172) (N=211) (N=383) Phylum Reads Percent (%) Reads Percent (%) Reads Percent (%) Actinobacteria 131,083 8.3 128,293 7.8 259,376 8.1 Proteobacteria 172,884 10.9 121,240 7.4 294,124 9.1 Bacteroidetes 182,466 11.6 167,913 10.3 350,379 10.9 Firmicutes 1,015,937 64.4 1,164,938 71.1 2,180,875 67.8 Fusobacteria 70,595 4.5 51,536 3.1 122,131 3.8 Others 5,388 0.3 5,001 0.3 10,389 0.3 Total 1,578,353 1,638,921 3,217,274 7 © 2017 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 8 Supplementary Table 3a. Median counts and fold changes (FC) for the association1 between selected bacteria related to 9 periodontal diseases and dental caries and risk of head and neck squamous cell cancer (HNSCC) in the 2 cohorts. 2 Transformed counts Median Maximum 3 4 Case Control FC 95% CI p q Cook’s D Porphyromonas;gingivalis 0.00 0.00 0.76 (0.42, 1.35) 0.3437 0.7767 8.83 Tannerella;forsythia 0.00 0.00 0.98 (0.59, 1.62) 0.9310 0.9796 5.60 Aggregatibacter;actinomycetemcomitans 0.00 0.00 0.97 (0.58, 1.63) 0.9099 0.9796 47.74 Streptococcus;mutans 5.49 7.82 0.66 (0.43, 1.02) 0.0635 0.4097 61.36 10 1. The association between taxonomic abundance and HNSCC was detected by DESeq function, adjusting for cohorts, matching factors [age at oral sample 11 collection (+/- 1 year), gender, and race (white, other)], smoking status (current, former, never), number of cigarettes per day for ever-smokers, alcohol 12 drinking status (never, ever), grams ethanol per day for ever-drinkers, and HPV-16 status. 13 2. Sequence read counts were normalized by dividing raw counts by DESeq size factors. 14 3. FDR-adjusted p value. FDR adjustment was conducted at each taxonomic level (i.e. class, genus) separately. 15 4. The maximum Cook’s Distance for DESeq model. 16 17 18 Supplementary Table 3b. Carriage1 rate and odds ratios (OR) for the association2 between selected bacteria related to 19 periodontal disease and dental caries and risk of head and neck squamous cell cancer (HNSCC) in the 2 cohorts Carriage, N (%) Cases Controls 2 2 2 OR 95% CI p-value (n=129) (n=254) Porphyromonas gingivalis 36 (27.9) 74 (29.1) 0.91 (0.51,1.61) 0.7392 Tannerella forsythia 49 (40.9) 104 (38.0) 0.91 (0.54, 1.54) 0.7199 Aggregatibacter actinomycetemcomitans 6 (4.7) 18 (7.1) 0.74 (0.24, 2.26) 0.5937 Streptococcus mutans 96 (74.4) 211 (83.1) 0.63 (0.33, 1.20) 0.1594 20 1. Selected taxa were coded as present/absent [1/0] for sequence reads from mouthwash samples. 21 2. Odds ratios, 95% confidence intervals, and p-values were calculated from logistic regression after controlling for matching factors [age at oral sample 22 collection (+/- 1 year), gender, and race (white, other)], smoking status (current, former, never), number of cigarettes per day for ever-smokers, alcohol 23 drinking status (never, ever), grams ethanol per day for ever-drinkers, and HPV-16 status. 24 25 © 2017 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 26 Supplementary Table 4. Median counts and fold changes (FC) for the association1 between selected taxa and risk of head 27 and neck squamous cell cancer (HNSCC) according to smoking status in the 2 cohorts Transformed 2 counts Taxon Median 3 Class;Order;Family;Genus;Species Case Control FC 95% CI p q Current smoker N=41 N=10 Actinobacteria;Corynebacteriales;Corynebacteriaceae;Corynebacterium 0.00 1.03 0.46 (0.11, 1.87) 0.2773 0.7571 Actinobacteria;Actinomycetales;Actinomycetaceae;Actinomyces;oral_taxon_170 0.00 0.00 2.84 (0.50, 15.98) 0.2368 0.7837 Betaproteobacteria 4.25 34.29 0.37 (0.15, 0.93) 0.0353 0.2434 Betaproteobacteria;Neisseriales 3.28 29.70 0.37 (0.16, 0.87) 0.0231 0.3935 Betaproteobacteria;Neisseriales;Neisseriaceae;Neisseria;sicca 0.00 4.12 0.14 (0.03, 0.71) 0.0174 0.2147 Betaproteobacteria;Neisseriales;Neisseriaceae;Kingella 1.10 3.62 0.58 (0.18, 1.89) 0.3663 0.7571 Bacteroidia;Bacteroidales;Prevotellaceae;Prevotella;nanceiensis 2.37 5.45 0.17 (0.04, 0.69) 0.0128 0.1787 Flavobacteriia;Flavobacteriales;Flavobacteriaceae;Capnocytophaga;leadbetteri 0.00 4.98 0.14 (0.04, 0.54) 0.0042 0.1505 Negativicutes Selenomonadales;Veillonellaceae;Selenomonas;sputigena 0.00 1.22 0.09 (0.02, 0.34) 0.0005 0.0296 Former smoker N=70 N=115 Actinobacteria;Corynebacteriales;Corynebacteriaceae;Corynebacterium 2.54 2.71 0.68 (0.44, 1.05) 0.0806 0.4238 Actinobacteria;Actinomycetales;Actinomycetaceae;Actinomyces;oral_taxon_170 0.00 0.65 1.67 (1.01, 2.76) 0.0444 0.4631 Betaproteobacteria 39.34 72.84 0.72 (0.55, 0.95) 0.0186 0.2046 Betaproteobacteria;Neisseriales 30.31 53.14 0.67 (0.48, 0.94) 0.0188 0.2765 Betaproteobacteria;Neisseriales;Neisseriaceae;Neisseria;sicca 5.22 10.80 0.56 (0.34, 0.93) 0.0242 0.3953 Betaproteobacteria;Neisseriales;Neisseriaceae;Kingella 4.71 5.48 0.54 (0.37, 0.81) 0.0026 0.1287 Bacteroidia;Bacteroidales;Prevotellaceae;Prevotella;nanceiensis 2.93 4.93 0.79 (0.49, 1.28) 0.3404 0.7986 Flavobacteriia;Flavobacteriales;Flavobacteriaceae;Capnocytophaga;leadbetteri 2.85 4.56 0.68 (0.44, 1.06) 0.0927 0.6192 Negativicutes Selenomonadales;Veillonellaceae;Selenomonas;sputigena 0.90 0.72 0.97 (0.61, 1.55) 0.8971 0.9709 Never smoker N=18 N=129 Actinobacteria;Corynebacteriales;Corynebacteriaceae;Corynebacterium 4.15 4.12 0.63 (0.35, 1.13) 0.1229 0.4902 Actinobacteria;Actinomycetales;Actinomycetaceae;Actinomyces;oral_taxon_170 1.97 1.95 0.96 (0.45, 2.04) 0.9075 0.9829 Betaproteobacteria 71.33 87.51 1.15 (0.83, 1.59) 0.3939 1.0000 Betaproteobacteria;Neisseriales 69.23 73.00 1.06 (0.68, 1.67) 0.7869 0.9598 Betaproteobacteria;Neisseriales;Neisseriaceae;Neisseria;sicca 18.51 15.52 0.97 (0.43, 2.18) 0.9452 0.9829 Betaproteobacteria;Neisseriales;Neisseriaceae;Kingella 5.50 5.81 1.07 (0.59, 1.93) 0.8164 0.9973 Bacteroidia;Bacteroidales;Prevotellaceae;Prevotella;nanceiensis 18.74 9.42 0.78 (0.37, 1.66) 0.5216 0.9702 Flavobacteriia;Flavobacteriales;Flavobacteriaceae;Capnocytophaga;leadbetteri 7.72 3.37 1.22 (0.61, 2.43) 0.5740 0.9702 Negativicutes Selenomonadales;Veillonellaceae;Selenomonas;sputigena 2.91 1.69 0.41 (0.19, 0.90) 0.0260 0.6494 28 Heterogeneity tests found no significant difference of fold changes between selected taxa-HNSCC associations according to smoking status by using Chi- 29 square test of Cochran’s Q statistic (p > 0.05). 30 1. The association between taxonomic abundance and HNSCC was detected by DESeq function, adjusting for cohorts, matching factors [age at oral sample 31 collection (+/- 1 year), gender, and race (white, other)], number of cigarettes per day for current- and former-smoker, alcohol drinking status (never, ever), 32 grams ethanol per day for ever-drinkers, and HPV-16 status.