Compositional and Functional Analysis of the Microbiome in Tissue and Saliva of Oral Squamous Cell Carcinoma

Next-generation sequencing technologies generate very large numbers of short reads. Even with very deep genome coverage, short read lengths cause problems in de novo assemblies. The use of paired-end libraries with a fragment size shorter than twice the read length provides an opportunity to generate much longer reads by overlapping and merging read pairs before assembling a genome. We present FLASH, a fast computational tool to extend the length of short reads by overlapping paired-end reads from fragment libraries that are sufficiently short. We tested the correctness of the tool on one million simulated read pairs, and we then applied it as a pre-processor for genome assemblies of Illumina reads from the bacterium Staphylococcus aureus and human chromosome 14. FLASH correctly extended and merged reads >99% of the time on simulated reads with an error rate of <1%. With adequately set parameters, FLASH correctly merged reads over 90% of the time even when the reads contained up to 5% errors. When FLASH was used to extend reads prior to assembly, the resulting assemblies had substantially greater N50 lengths for both contigs and scaffolds. The FLASH system is implemented in C and is freely available as open-source code at http://www.cbcb.umd.edu/software/flash. t.magoc@gmail.com.

Background Most studies examining the commensal human oral microbiome are focused on disease or are limited in methodology. In order to diagnose and treat diseases at an early and reversible stage an in-depth definition of health is indispensible. The aim of this study therefore was to define the healthy oral microbiome using recent advances in sequencing technology (454 pyrosequencing). Results We sampled and sequenced microbiomes from several intraoral niches (dental surfaces, cheek, hard palate, tongue and saliva) in three healthy individuals. Within an individual oral cavity, we found over 3600 unique sequences, over 500 different OTUs or "species-level" phylotypes (sequences that clustered at 3% genetic difference) and 88 - 104 higher taxa (genus or more inclusive taxon). The predominant taxa belonged to Firmicutes (genus Streptococcus, family Veillonellaceae, genus Granulicatella), Proteobacteria (genus Neisseria, Haemophilus), Actinobacteria (genus Corynebacterium, Rothia, Actinomyces), Bacteroidetes (genus Prevotella, Capnocytophaga, Porphyromonas) and Fusobacteria (genus Fusobacterium). Each individual sample harboured on average 266 "species-level" phylotypes (SD 67; range 123 - 326) with cheek samples being the least diverse and the dental samples from approximal surfaces showing the highest diversity. Principal component analysis discriminated the profiles of the samples originating from shedding surfaces (mucosa of tongue, cheek and palate) from the samples that were obtained from solid surfaces (teeth). There was a large overlap in the higher taxa, "species-level" phylotypes and unique sequences among the three microbiomes: 84% of the higher taxa, 75% of the OTUs and 65% of the unique sequences were present in at least two of the three microbiomes. The three individuals shared 1660 of 6315 unique sequences. These 1660 sequences (the "core microbiome") contributed 66% of the reads. The overlapping OTUs contributed to 94% of the reads, while nearly all reads (99.8%) belonged to the shared higher taxa. Conclusions We obtained the first insight into the diversity and uniqueness of individual oral microbiomes at a resolution of next-generation sequencing. We showed that a major proportion of bacterial sequences of unrelated healthy individuals is identical, supporting the concept of a core microbiome at health.

Human papillomavirus (HPV), the causal agent of cervical cancer, appears to be involved in the etiology of cancer of the oral cavity and oropharynx. To investigate these associations, we conducted a multicenter case-control study of cancer of the oral cavity and oropharynx in nine countries. We recruited 1670 case patients (1415 with cancer of the oral cavity and 255 with cancer of the oropharynx) and 1732 control subjects and obtained an interview, oral exfoliated cells, and blood from all participants and fresh biopsy specimens from case patients. HPV DNA was detected by polymerase chain reaction (PCR). Antibodies against HPV16 L1, E6, and E7 proteins in plasma were detected with enzyme-linked immunosorbent assays. Multivariable models were used for case-control and case-case comparisons. HPV DNA was detected in biopsy specimens of 3.9% (95% confidence interval [CI] = 2.5% to 5.3%) of 766 cancers of the oral cavity with valid PCR results and 18.3% (95% CI = 12.0% to 24.7%) of 142 cancers of the oropharynx (oropharynx and tonsil combined) with valid PCR results. HPV DNA in cancer biopsy specimens was detected less frequently among tobacco smokers and paan chewers and more frequently among subjects who reported more than one sexual partner or who practiced oral sex. HPV16 DNA was found in 94.7% of HPV DNA-positive case patients. HPV DNA in exfoliated cells was not associated with cancer risk or with HPV DNA detection in biopsy specimens. Antibodies against HPV16 L1 were associated with risk for cancers of the oral cavity (odds ratio [OR] = 1.5, 95% CI = 1.1 to 2.1) and the oropharynx (OR = 3.5, 95% CI = 2.1 to 5.9). Antibodies against HPV16 E6 or E7 were also associated with risk for cancers of the oral cavity (OR = 2.9, 95% CI = 1.7 to 4.8) and the oropharynx (OR = 9.2, 95% CI = 4.8 to 17.7). HPV appears to play an etiologic role in many cancers of the oropharynx and possibly a small subgroup of cancers of the oral cavity. The most common HPV type in genital cancers (HPV16) was also the most common in these tumors. The mechanism of transmission of HPV to the oral cavity warrants further investigation.