Diana Tabbaa

Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons

Bacterial diversity among environmental samples is commonly assessed with PCR-amplified 16S rRNA gene (16S) sequences. Perceived diversity, however, can be influenced by sample preparation, primer selection, and formation of chimeric 16S amplification products. Chimeras are hybrid products between multiple parent sequences that can be falsely interpreted as novel organisms, thus inflating apparent diversity. We developed a new chimera detection tool called Chimera Slayer (CS). CS detects chimeras with greater sensitivity than previous methods, performs well on short sequences such as those produced by the 454 Life Sciences (Roche) Genome Sequencer, and can scale to large data sets. By benchmarking CS performance against sequences derived from a controlled DNA mixture of known organisms and a simulated chimera set, we provide insights into the factors that affect chimera formation such as sequence abundance, the extent of similarity between 16S genes, and PCR conditions. Chimeras were found to reproducibly form among independent amplifications and contributed to false perceptions of sample diversity and the false identification of novel taxa, with less-abundant species exhibiting chimera rates exceeding 70%. Shotgun metagenomic sequences of our mock community appear to be devoid of 16S chimeras, supporting a role for shotgun metagenomics in validating novel organisms discovered in targeted sequence surveys.

SNP Genotyping Using the Sequenom MassARRAY iPLEX Platform

The method for SNP genotyping described in this unit is based on the commercially available Sequenom MassARRAY platform. The assay consists of an initial locus‐specific PCR reaction, followed by single base extension using mass‐modified dideoxynucleotide terminators of an oligonucleotide primer which anneals immediately upstream of the polymorphic site of interest. Using MALDI‐TOF mass spectrometry, the distinct mass of the extended primer identifies the SNP allele. Curr. Protoc. Hum. Genet. 60:2.12.1‐2.12.18.

The African coelacanth genome provides insights into tetrapod evolution.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

Evaluation of 16s rDNA-based community profiling for human microbiome research

The Human Microbiome Project will establish a reference data set for analysis of the microbiome of healthy adults by surveying multiple body sites from 300 people and generating data from over 12,000 samples. To characterize these samples, the participating sequencing centers evaluated and adopted 16S rDNA community profiling protocols for ABI 3730 and 454 FLX Titanium sequencing. In the course of establishing protocols, we examined the performance and error characteristics of each technology, and the relationship of sequence error to the utility of 16S rDNA regions for classification- and OTU-based analysis of community structure. The data production protocols used for this work are those used by the participating centers to produce 16S rDNA sequence for the Human Microbiome Project. Thus, these results can be informative for interpreting the large body of clinical 16S rDNA data produced for this project.

Non-invasive mapping of the gastrointestinal microbiota identifies children with inflammatory bowel disease.

Background Pediatric inflammatory bowel disease (IBD) is challenging to diagnose because of the non-specificity of symptoms; an unequivocal diagnosis can only be made using colonoscopy, which clinicians are reluctant to recommend for children. Diagnosis of pediatric IBD is therefore frequently delayed, leading to inappropriate treatment plans and poor outcomes. We investigated the use of 16S rRNA sequencing of fecal samples and new analytical methods to assess differences in the microbiota of children with IBD and other gastrointestinal disorders. Methodology/Principal Findings We applied synthetic learning in microbial ecology (SLiME) analysis to 16S sequencing data obtained from i) published surveys of microbiota diversity in IBD and ii) fecal samples from 91 children and young adults who were treated in the gastroenterology program of Children’s Hospital (Boston, USA). The developed method accurately distinguished control samples from those of patients with IBD; the area under the receiver-operating-characteristic curve (AUC) value was 0.83 (corresponding to 80.3% sensitivity and 69.7% specificity at a set threshold). The accuracy was maintained among data sets collected by different sampling and sequencing methods. The method identified taxa associated with disease states and distinguished patients with Crohn’s disease from those with ulcerative colitis with reasonable accuracy. The findings were validated using samples from an additional group of 68 patients; the validation test identified patients with IBD with an AUC value of 0.84 (e.g. 92% sensitivity, 58.5% specificity). Conclusions/Significance Microbiome-based diagnostics can distinguish pediatric patients with IBD from patients with similar symptoms. Although this test can not replace endoscopy and histological examination as diagnostic tools, classification based on microbial diversity is an effective complementary technique for IBD detection in pediatric patients.

Paired-end sequencing of Fosmid libraries by Illumina

Eliminating the bacterial cloning step has been a major factor in the vastly improved efficiency of massively parallel sequencing approaches. However, this also has made it a technical challenge to produce the modern equivalent of the Fosmid- or BAC-end sequences that were crucial for assembling and analyzing complex genomes during the Sanger-based sequencing era. To close this technology gap, we developed Fosill, a method for converting Fosmids to Illumina-compatible jumping libraries. We constructed Fosmid libraries in vectors with Illumina primer sequences and specific nicking sites flanking the cloning site. Our family of pFosill vectors allows multiplex Fosmid cloning of end-tagged genomic fragments without physical size selection and is compatible with standard and multiplex paired-end Illumina sequencing. To excise the bulk of each cloned insert, we introduced two nicks in the vector, translated them into the inserts, and cleaved them. Recircularization of the vector via coligation of insert termini followed by inverse PCR generates a jumping library for paired-end sequencing with 101-base reads. The yield of unique Fosmid-sized jumps is sufficiently high, and the background of short, incorrectly spaced and chimeric artifacts sufficiently low, to enable applications such as mapping of structural variation and scaffolding of de novo assemblies. We demonstrate the power of Fosill to map genome rearrangements in a cancer cell line and identified three fusion genes that were corroborated by RNA-seq data. Our Fosill-powered assembly of the mouse genome has an N50 scaffold length of 17.0 Mb, rivaling the connectivity (16.9 Mb) of the Sanger-sequencing based draft assembly.

A cloning and expression system to probe T cell receptor specificity and assess functional avidity to neoantigens

Recent studies have highlighted the promise of targeting tumor neoantigens to generate potent antitumor immune responses and provide strong motivation for improving our understanding of antigen-T-cell receptor (TCR) interactions. Advances in single-cell sequencing technologies have opened the door for detailed investigation of the TCR repertoire, providing paired information from TCRα and TCRβ, which together determine specificity. However, a need remains for efficient methods to assess the specificity of discovered TCRs. We developed a streamlined approach for matching TCR sequences with cognate antigen through on-demand cloning and expression of TCRs and screening against candidate antigens. Here, we first demonstrate the systems capacity to identify viral-antigen-specific TCRs and compare the functional avidity of TCRs specific for a given antigen target. We then apply this system to identify neoantigen-specific TCR sequences from patients with melanoma treated with personalized neoantigen vaccines and characterize functional avidity of neoantigen-specific TCRs. Furthermore, we use a neoantigen-prediction pipeline to show that an insertion-deletion mutation in a putative chronic lymphocytic leukemia (CLL) driver gives rise to an immunogenic neoantigen mut-MGA, and use this approach to identify the mut-MGA-specific TCR sequence. This approach provides a means to identify and express TCRs, and then rapidly assess antigen specificity and functional avidity of a reconstructed TCR, which can be applied for monitoring antigen-specific T-cell responses, and potentially for guiding the design of effective T-cell-based immunotherapies.