Faculty Bibliography

Functional analysis of a clinical microbiome facilitates the elucidation of mechanisms by which microbiome perturbation can cause a phenotypic change in the patient. The direct approach for the analysis of the functional capacity of the microbiome is via shotgun metagenomics. An inexpensive method to estimate the functional capacity of a microbial community is through collecting 16S rRNA gene profiles then indirectly inferring the abundance of functional genes. This inference approach has been implemented in the PICRUSt and Tax4Fun software tools. However, those tools have important limitations since they rely on outdated functional databases and uncertain phylogenetic trees and require very specific data pre-processing protocols. Here we introduce Piphillin, a straightforward algorithm independent of any proposed phylogenetic tree, leveraging contemporary functional databases and not obliged to any singular data pre-processing protocol. When all three inference tools were evaluated against actual shotgun metagenomics, Piphillin was superior in predicting gene composition in human clinical samples compared to both PICRUSt and Tax4Fun (p<0.01 and p<0.001, respectively) and Piphillin's ability to predict disease associations with specific gene orthologs exhibited a 15% increase in balanced accuracy compared to PICRUSt. From laboratory animal samples, no performance advantage was observed for any one of the tools over the others and for environmental samples all produced unsatisfactory predictions. Our results demonstrate that functional inference using the direct method implemented in Piphillin is preferable for clinical biospecimens. Piphillin is publicly available for academic use at http://secondgenome.com/Piphillin.

Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures. Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice, very useful profiles can be obtained with ~0.1x genome coverage. We improve on previous methods by; first, implementing a combined correction for sequence mappability and GC content, and second, applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions were previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1,000 samples, most of which were obtained from the fixed tissue archives of over 25 institutions. We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches, and better copy number data than high-resolution microarrays at substantially lower cost.

beta-defensins are a family of important peptides of innate immunity, involved in host defense, immunomodulation, reproduction and pigmentation. Genes encoding beta-defensins show evidence of birth-and-death evolution, adaptation by amino acid sequence changes and extensive copy number variation (CNV) within humans and other species. The role of CNV in the adaptation of beta-defensins to new functions remains unclear, as does the adaptive role of copy number variation in general. Here, we fine-map CNV of a cluster of beta-defensins in humans and rhesus macaques. Remarkably, we found that the structure of the CNV is different between primates, with distinct mutational origins and CNV boundaries defined by retroviral long terminal repeat elements. While the human beta-defensin CNV region is 322 kb and encompasses several genes, including beta-defensins, a long non-coding RNA gene and testes-specific zinc-finger transcription factors, the orthologous region in the rhesus macaque shows CNV of a 20 kb region, containing only a single gene, the orthologue of the human beta-defensin-2 gene. Despite its independent origins, the range of gene copy numbers in the rhesus macaque is similar to humans. In addition, the rhesus macaque gene has been subject to divergent positive selection at the amino acid level following its initial duplication event between 3 and 9.5 million years ago, suggesting adaptation of this gene as the macaque successfully colonised novel environments outside Africa. Therefore, the molecular phenotype of beta-defensin-2 copy number variation has undergone convergent evolution, and this gene shows evidence of adaptation at the amino acid level in both primates.

BACKGROUND: Metastasis to the cervical (neck) lymph nodes is one of the most significant clinical factors responsible for death from oral squamous cell carcinoma (SCC). Therefore, the lymph nodes are frequently removed when the tumor is excised (neck dissection), even though the majority of patients will not benefit from the extra surgery. Two subtypes of oral SCC distinguished by the presence of tumor genomic aberrations +3q, -8p, +8q and/or +20 differ in risk for metastasis - high for the 3q8pq20 subtype, harboring one or more of the aberrations and low for the non-3q8pq20 subtype, lacking these alterations. A prior analysis of the literature suggested genes differentially methylated in the two subtypes. Therefore, the goal of this study was to further investigate the methylation status of candidate biomarkers of the non-3q8pq20 subtype, and evaluate their utility for identifying patients at low risk for metastasis. METHODS: Methylation status of genes in a cohort of 52 oral SCC patients with at least five year follow up was determined by pyrosequencing. Gene expression levels were determined by quantitative RT-PCR. Growth following re-expression of HOXA9 in cultured oral SCC cells was assessed by proliferation and colony formation assays. RESULTS: A pilot study evaluating methylation levels of HOXA9, MT1A and HOXA11 promoters in DNA from 12 tumors (six each of the 3q8pq20 and non-3q8pq20 subtypes) revealed that only HOXA9 was differentially methylated. Significant differences in methylation levels of HOXA9 were observed amongst the 52 oral SCCs with respect to genomic subtype and nodal status (p = 0.014, and p = 0.024, respectively, Wilcoxon rank sum test). High levels of HOXA9 methylation and low levels of expression in oral SCC cell lines were observed compared to HaCaT, a non-tumorigenic keratinocyte cell line. Re-expression of HOXA9 in the SCC4 oral cancer cell line resulted in diminished proliferation and colony formation. CONCLUSIONS: HOXA9 methylation is frequent in oral cancers and levels are higher in tumors with greater risk of metastasis. Expression of HOXA9 is low in cells with high levels of methylation and reduced expression appears to confer a growth advantage.

Individual bacteria and shifts in the composition of the microbiome have been associated with human diseases including cancer. To investigate changes in the microbiome associated with oral cancers, we profiled cancers and anatomically matched contralateral normal tissue from the same patient by sequencing 16S rDNA hypervariable region amplicons. In cancer samples from both a discovery and a subsequent confirmation cohort, abundance of Firmicutes (especially Streptococcus) and Actinobacteria (especially Rothia) was significantly decreased relative to contralateral normal samples from the same patient. Significant decreases in abundance of these phyla were observed for pre-cancers, but not when comparing samples from contralateral sites (tongue and floor of mouth) from healthy individuals. Weighted UniFrac principal coordinates analysis based on 12 taxa separated most cancers from other samples with greatest separation of node positive cases. These studies begin to develop a framework for exploiting the oral microbiome for monitoring oral cancer development, progression and recurrence.

In order to understand the nervous system, it is necessary to know the synaptic connections between the neurons, yet to date, only the wiring diagram of the adult hermaphrodite of the nematode Caenorhabditis elegans has been determined. Here, we present the wiring diagram of the posterior nervous system of the C. elegans adult male, reconstructed from serial electron micrograph sections. This region of the male nervous system contains the sexually dimorphic circuits for mating. The synaptic connections, both chemical and gap junctional, form a neural network with four striking features: multiple, parallel, short synaptic pathways directly connecting sensory neurons to end organs; recurrent and reciprocal connectivity among sensory neurons; modular substructure; and interneurons acting in feedforward loops. These features help to explain how the network robustly and rapidly selects and executes the steps of a behavioral program on the basis of the inputs from multiple sensory neurons.