mRNA was extracted from biopsy samples and used for RNA-Seq library preparation following instructions in the Illumina
TruSeq RNA Library Prep Kit v2. Sequencing was run on an Illumina
High Seq-2000 in the Mayo Clinic Sequencing Core with 101bp paired end reads. Gene expression counts were obtained using the MAP-RSeq v.2.0.0 workflow (Kalari et al., 2014 (
link)) which is part of the Mayo Bioinformatics Core pipeline. MAP-RSeq consists of alignment with TopHat 2.0.12 (Kim et al., 2013 (
link)) against the human hg19 genome build and gene counts with the Subread package 1.4.4 (Liao et al., 2019 (
link)). Gene annotation files were obtained from Ensemble version 75. Gene counts were normalized using RPKM (Reads Per Kilobase per million Mapped reads). Differential expression analysis was performed using edgeR 2.6.2 (Robinson et al., 2010 (
link)). Pathway enrichment analyses were performed using R package RITAN (Rapid Integration of Term Annotation and Network resources,
https://www.bioconductor.org/packages/release/bioc/html/RITAN.html).
Mars R.A., Yang Y., Ward T., Houtti M., Priya S., Lekatz H.R., Tang X., Sun Z., Kalari K.R., Korem T., Bhattarai Y., Zheng T., Bar N., Frost G., Johnson A.J., van Treuren W., Han S., Ordog T., Grover M., Sonnenburg J., D’Amato M., Camilleri M., Elinav E., Segal E., Blekhman R., Farrugia G., Swann J.R., Knights D, & Kashyap P.C. (2020). Longitudinal Multi-omics Reveals Subset-Specific Mechanisms Underlying Irritable Bowel Syndrome. Cell, 182(6), 1460-1473.e17.