Jordan Patterson

SOAPdenovo-Trans: De novo transcriptome assembly with short RNA-Seq reads

MOTIVATION Transcriptome sequencing has long been the favored method for quickly and inexpensively obtaining a large number of gene sequences from an organism with no reference genome. Owing to the rapid increase in throughputs and decrease in costs of next- generation sequencing, RNA-Seq in particular has become the method of choice. However, the very short reads (e.g. 2 × 90 bp paired ends) from next generation sequencing makes de novo assembly to recover complete or full-length transcript sequences an algorithmic challenge. RESULTS Here, we present SOAPdenovo-Trans, a de novo transcriptome assembler designed specifically for RNA-Seq. We evaluated its performance on transcriptome datasets from rice and mouse. Using as our benchmarks the known transcripts from these well-annotated genomes (sequenced a decade ago), we assessed how SOAPdenovo-Trans and two other popular transcriptome assemblers handled such practical issues as alternative splicing and variable expression levels. Our conclusion is that SOAPdenovo-Trans provides higher contiguity, lower redundancy and faster execution. AVAILABILITY AND IMPLEMENTATION Source code and user manual are available at http://sourceforge.net/projects/soapdenovotrans/.

Characterization of the Gut Microbiome Using 16S or Shotgun Metagenomics.

The advent of next generation sequencing (NGS) has enabled investigations of the gut microbiome with unprecedented resolution and throughput. This has stimulated the development of sophisticated bioinformatics tools to analyze the massive amounts of data generated. Researchers therefore need a clear understanding of the key concepts required for the design, execution and interpretation of NGS experiments on microbiomes. We conducted a literature review and used our own data to determine which approaches work best. The two main approaches for analyzing the microbiome, 16S ribosomal RNA (rRNA) gene amplicons and shotgun metagenomics, are illustrated with analyses of libraries designed to highlight their strengths and weaknesses. Several methods for taxonomic classification of bacterial sequences are discussed. We present simulations to assess the number of sequences that are required to perform reliable appraisals of bacterial community structure. To the extent that fluctuations in the diversity of gut bacterial populations correlate with health and disease, we emphasize various techniques for the analysis of bacterial communities within samples (α-diversity) and between samples (β-diversity). Finally, we demonstrate techniques to infer the metabolic capabilities of a bacteria community from these 16S and shotgun data.

Evaluating methods for isolating total RNA and predicting the success of sequencing phylogenetically diverse plant transcriptomes

Next-generation sequencing plays a central role in the characterization and quantification of transcriptomes. Although numerous metrics are purported to quantify the quality of RNA, there have been no large-scale empirical evaluations of the major determinants of sequencing success. We used a combination of existing and newly developed methods to isolate total RNA from 1115 samples from 695 plant species in 324 families, which represents >900 million years of phylogenetic diversity from green algae through flowering plants, including many plants of economic importance. We then sequenced 629 of these samples on Illumina GAIIx and HiSeq platforms and performed a large comparative analysis to identify predictors of RNA quality and the diversity of putative genes (scaffolds) expressed within samples. Tissue types (e.g., leaf vs. flower) varied in RNA quality, sequencing depth and the number of scaffolds. Tissue age also influenced RNA quality but not the number of scaffolds ≥1000 bp. Overall, 36% of the variation in the number of scaffolds was explained by metrics of RNA integrity (RIN score), RNA purity (OD 260/230), sequencing platform (GAIIx vs HiSeq) and the amount of total RNA used for sequencing. However, our results show that the most commonly used measures of RNA quality (e.g., RIN) are weak predictors of the number of scaffolds because Illumina sequencing is robust to variation in RNA quality. These results provide novel insight into the methods that are most important in isolating high quality RNA for sequencing and assembling plant transcriptomes. The methods and recommendations provided here could increase the efficiency and decrease the cost of RNA sequencing for individual labs and genome centers.

Effect of oral capsule– vs colonoscopy-delivered fecal microbiota transplantation on recurrent Clostridium difficile infection: A randomized clinical trial

Importance Fecal microbiota transplantation (FMT) is effective in preventing recurrent Clostridium difficile infection (RCDI). However, it is not known whether clinical efficacy differs by route of delivery. Objective To determine whether FMT by oral capsule is noninferior to colonoscopy delivery in efficacy. Design, Setting, and Participants Noninferiority, unblinded, randomized trial conducted in 3 academic centers in Alberta, Canada. A total of 116 adult patients with RCDI were enrolled between October 2014 and September 2016, with follow-up to December 2016. The noninferiority margin was 15%. Interventions Participants were randomly assigned to FMT by capsule or by colonoscopy at a 1:1 ratio. Main Outcomes and Measures The primary outcome was the proportion of patients without RCDI 12 weeks after FMT. Secondary outcomes included (1) serious and minor adverse events, (2) changes in quality of life by the 36-Item Short Form Survey on a scale of 0 (worst possible quality of life) to 100 (best quality of life), and (3) patient perception on a scale of 1 (not at all unpleasant) to 10 (extremely unpleasant) and satisfaction on a scale of 1 (best) to 10 (worst). Results Among 116 patients randomized (mean [SD] age, 58 [19] years; 79 women [68%]), 105 (91%) completed the trial, with 57 patients randomized to the capsule group and 59 to the colonoscopy group. In per-protocol analysis, prevention of RCDI after a single treatment was achieved in 96.2% in both the capsule group (51/53) and the colonoscopy group (50/52) (difference, 0%; 1-sided 95% CI, −6.1% to infinity; P < .001), meeting the criterion for noninferiority. One patient in each group died of underlying cardiopulmonary illness unrelated to FMT. Rates of minor adverse events were 5.4% for the capsule group vs 12.5% for the colonoscopy group. There was no significant between-group difference in improvement in quality of life. A significantly greater proportion of participants receiving capsules rated their experience as “not at all unpleasant” (66% vs 44%; difference, 22% [95% CI, 3%-40%]; P = .01). Conclusions and Relevance Among adults with RCDI, FMT via oral capsules was not inferior to delivery by colonoscopy for preventing recurrent infection over 12 weeks. Treatment with oral capsules may be an effective approach to treating RCDI. Trial Registration clinicaltrials.gov Identifier: NCT02254811

Metagenomic Analysis of Microbiome in Colon Tissue from Subjects with Inflammatory Bowel Diseases Reveals Interplay of Viruses and Bacteria

Abstract:Inflammatory bowel diseases (IBD), Crohns disease and ulcerative colitis, are poorly understood disorders affecting the intestinal tract. The current model for disease suggests that genetically susceptible patients develop intolerance to gut microflora, and chronic inflammation develops as a result of environmental insults. Although interest has mainly focused on studying genetic variants and gut bacterial flora, little is known about the potential of viral infection to contribute to disease. Accordingly, we conducted a metagenomic analysis to document the baseline virome in colonic biopsy samples from patients with IBD in order to assess the contribution of viral infection to IBD. Libraries were generated from colon RNA to create approximately 2 GB sequence data per library. Using a bioinformatic pipeline designed to detect viral sequences, more than 1000 viral reads were derived directly from tissue without any coculture or isolation procedure. Herein, we describe the complexity and abundance of viruses, bacteria/bacteriophage, and human endogenous retroviral sequences from 10 patients with IBD and 5 healthy subjects undergoing surveillance colonoscopy. Differences in gut microflora and the abundance of mammalian viruses and human endogenous retroviruses were readily detected in the metagenomic analyses. Specifically, patients with herpesviridae sequences in their colon demonstrated increased expression of human endogenous viral sequences and differences in the diversity of their microbiome. This study provides a promising metagenomic approach to describe the colonic microbiome that can be used to better understand virus–host and phage–bacteria interactions in IBD.

Identification of Hepatotropic Viruses from Plasma Using Deep Sequencing: A Next Generation Diagnostic Tool

We conducted an unbiased metagenomics survey using plasma from patients with chronic hepatitis B, chronic hepatitis C, autoimmune hepatitis (AIH), non-alcoholic steatohepatitis (NASH), and patients without liver disease (control). RNA and DNA libraries were sequenced from plasma filtrates enriched in viral particles to catalog virus populations. Hepatitis viruses were readily detected at high coverage in patients with chronic viral hepatitis B and C, but only a limited number of sequences resembling other viruses were found. The exception was a library from a patient diagnosed with hepatitis C virus (HCV) infection that contained multiple sequences matching GB virus C (GBV-C). Abundant GBV-C reads were also found in plasma from patients with AIH, whereas Torque teno virus (TTV) was found at high frequency in samples from patients with AIH and NASH. After taxonomic classification of sequences by BLASTn, a substantial fraction in each library, ranging from 35% to 76%, remained unclassified. These unknown sequences were assembled into scaffolds along with virus, phage and endogenous retrovirus sequences and then analyzed by BLASTx against the non-redundant protein database. Nearly the full genome of a heretofore-unknown circovirus was assembled and many scaffolds that encoded proteins with similarity to plant, insect and mammalian viruses. The presence of this novel circovirus was confirmed by PCR. BLASTx also identified many polypeptides resembling nucleo-cytoplasmic large DNA viruses (NCLDV) proteins. We re-evaluated these alignments with a profile hidden Markov method, HHblits, and observed inconsistencies in the target proteins reported by the different algorithms. This suggests that sequence alignments are insufficient to identify NCLDV proteins, especially when these alignments are only to small portions of the target protein. Nevertheless, we have now established a reliable protocol for the identification of viruses in plasma that can also be adapted to other patient samples such as urine, bile, saliva and other body fluids.

Modeling the Yew Tree Tubulin and a Comparison of its Interaction with Paclitaxel to Human Tubulin

ABSTRACTPurposeTo explore possible ways in which yew tree tubulin is naturally resistant to paclitaxel. While the yew produces a potent cytotoxin, paclitaxel, it is immune to paclitaxel’s cytotoxic action.MethodsTubulin sequence data for plant species were obtained from Alberta 1000 Plants Initiative. Sequences were assembled with Trinity de novo assembly program and tubulin identified. Homology modeling using MODELLER software was done to generate structures for yew tubulin. Molecular dynamics simulations and molecular mechanics Poisson–Boltzmann calculations were performed with the Amber package to determine binding affinity of paclitaxel to yew tubulin. ClustalW2 program and PHYLIP package were used to perform phylogenetic analysis on plant tubulin sequences.ResultsWe specifically analyzed several important regions in tubulin structure: the high-affinity paclitaxel binding site, as well as the intermediate binding site and microtubule nanopores. Our analysis indicates that the high-affinity binding site contains several substitutions compared to human tubulin, all of which reduce the binding energy of paclitaxel.ConclusionsThe yew has achieved a significant reduction of paclitaxel’s affinity for its tubulin by utilizing several specific residue changes in the binding pocket for paclitaxel.

Fecal Microbial Transplant After Ileocolic Resection Reduces Ileitis but Restores Colitis in IL-10-/- Mice.

Background:Ileocolic resection (ICR) is frequently performed for Crohns disease; however, disease commonly recurs early in the neoterminal ileum. The aim of this study was to use the IL-10−/− mouse to determine the effects of ICR on gut microbiome and immune function and if postoperative fecal microbial transplant (FMT) would improve disease. Methods:ICR was performed in 129S1/SvlmJ IL10−/− mice followed by FMT using stool from wild-type mice. Sham-transplant mice received their own stool. Stool samples were collected on day 0, day 13 (after ICR), and day 27 (after FMT) for whole metagenome shot-gun sequencing. Mucosal-associated bacteria were quantified with quantitative PCR and visualized by fluorescent in situ hybridization. Tissue cytokines were measured with multiplex arrays and mononuclear phagocyte populations by flow cytometry. Results:Surgery induced microbial functional and taxonomic shifts, decreased diversity, and depleted Bacteroidia and Clostridia. ICR mice had reduced colitis but worse ileitis with bacterial overgrowth, increased translocation, and reduction in tissue macrophages. FMT prevented ileitis but restored colitis and allowed for a bloom of &ggr;-proteobacteria. In the colon, ICR and sham transplant were associated with recruitment of tolerogenic dendritic cells, whereas FMT shifted these immune cell subsets to control profiles along with increasing cytokine levels. Conclusions:This study suggests that surgical-induced immune dysfunction and microbial dysbiosis with impaired clearance may be the underlying cause of the early ulcerations found in the ileum of patients with Crohns disease after ICR. FMT has an immunostimulatory effect on the postoperative intestine, which was beneficial in preventing ileitis, but detrimental in restoring colonic injury after surgery.

Novel mutations involving βI-, βIIA-, or βIVB-tubulin isotypes with functional resemblance to βIII-tubulin in breast cancer

Tubulin is the target for very widely used anti-tumor drugs, including Vinca alkaloids, taxanes, and epothilones, which are an important component of chemotherapy in breast cancer and other malignancies. Paclitaxel and other tubulin-targeting drugs bind to the β subunit of tubulin, which is a heterodimer of α and β subunits. β-Tubulin exists in the form of multiple isotypes, which are differentially expressed in normal and neoplastic cells and differ in their ability to bind to drugs. Among them, the βIII isotype is overexpressed in many aggressive and metastatic cancers and may serve as a prognostic marker in certain types of cancer. The underpinning mechanisms accounting for the overexpression of this isotype in cancer cells are unclear. To better understand the role of β-tubulin isotypes in cancer, we analyzed over 1000 clones from 90 breast cancer patients, sequencing their β-tubulin isotypes, in search of novel mutations. We have elucidated two putative emerging molecular subgroups of invasive breast cancer, each of which involve mutations in the βI-, βIIA-, or βIVB isotypes of tubulin that increase their structural, and possibly functional, resemblance to the βIII isotype. A unifying feature of the first of the two subgroups is the mutation of the highly reactive C239 residue of βI- or βIVB-tubulin to L239, R239, Y239, or P239, culminating in probable conversion of these isotypes from ROS-sensitive to ROS-resistant species. In the second subgroup, βI, βIIA, and βIVB have up to seven mutations to the corresponding residues in βIII-tubulin. Given that βIII-tubulin has emerged as a pro-survival factor, overexpression of this isotype may confer survival advantages to certain cancer cell types. In this mini-review, we bring attention to a novel mechanism by which cancer cells may undergo adaptive mutational changes involving alternate β-tubulin isotypes to make them acquire some of the pro-survival properties of βIII-tubulin. These “hybrid” tubulins, combining the sequences and/or properties of two wild-type tubulins (βIII and either βI, βIIA, or βIVB), are novel isotypes expressed solely in cancer cells and may contribute to the molecular understanding and stratification of invasive breast cancer and provide novel molecular targets for rational drug development.

Cerebrospinal Fluid in a Small Cohort of Patients with Multiple Sclerosis Was Generally Free of Microbial DNA

Multiple sclerosis (MS) is a common cause of non-traumatic neurologic disability with high incidence in many developed countries. Although the etiology of the disease remains elusive, it is thought to entail genetic and environmental causes, and microbial pathogens have also been envisioned as contributors to the phenotype. We conducted a metagenomic survey in cerebrospinal fluid (CSF) from 28 MS patients and 15 patients suffering other type of neurological conditions. We detected bacterial reads in eight out of the 15 non-MS patients and in a single MS patient, at an abundance >1% of total classified reads. Two patients were of special interest: one non-MS patient harbored ~73% bacterial reads, while an MS patient had ~83% bacterial reads. In the former case, Veillonella parvula, a bacterium occasionally found associated with meningitis was the predominant species, whilst Kocuria flava, apparently an environmental bacterium, predominated in the latter case. Thirty-four out of 43 samples contained <1% bacterial reads, which we regard as cross- or environmental contamination. A few viral reads corresponding to Epstein-Barr virus, cytomegalovirus, and parvovirus were also identified. Our results suggest that CSF of MS patients is often (but not always) free of microbial DNA.