Baheti, Saurabh

Abstract Motivation: RNA-seq has become the method of choice to quantify genes and exons, discover novel transcripts and detect fusion genes. However, reliable variant identification from RNA-seq data remains challenging because of the complexities of the transcriptome, the challenges of accurately mapping exon boundary spanning reads and the bias introduced during the sequencing library preparation. Method: We developed RVboost, a novel method specific for RNA variant prioritization. RVboost uses several attributes unique in the process of RNA library preparation, sequencing and RNA-seq data analyses. It uses a boosting method to train a model of ‘good quality’ variants using common variants from HapMap, and prioritizes and calls the RNA variants based on the trained model. We packaged RVboost in a comprehensive workflow, which integrates tools of variant calling, annotation and filtering. Results: RVboost consistently outperforms the variant quality score recalibration from the Genome Analysis Tool Kit and the RNA-seq variant-calling pipeline SNPiR in 12 RNA-seq samples using ground-truth variants from paired exome sequencing data. Several RNA-seq–specific attributes were identified as critical to differentiate true and false variants, including the distance of the variant positions to exon boundaries, and the percent of the reads supporting the variant in the first six base pairs. The latter identifies false variants introduced by the random hexamer priming during the library construction. Availability and implementation: The RVboost package is implemented to readily run in Mac or Linux environments. The software and user manual are available at http://bioinformaticstools.mayo.edu/research/rvboost/. Supplementary information: Supplementary data are available at Bioinformatics online.

Abstract MOTIVATION: Transcriptome sequencing (RNA-seq) is the most economical and information-rich next generation sequencing technology to study genomic medicine. Genomic features such as gene expression, transcript expression, fusion transcripts, long non-coding RNAs and expressed single nucleotide variants (eSNVs) can be obtained from RNA-seq data. There are several computational workflows to obtain gene counts, transcript counts and fusion gene candidates, but not a comprehensive system to identify expressed SNVs from RNA-sequencing data. METHODS: We have developed a novel computational system (The eSNV-Detect) with parallel alignment strategies to identify expressed SNVs accurately from tumor transcriptome sequencing data. Sequence-based features like reads supporting reference allele, alternate allele, ratio of alternate reads to total read depth, base qualities, mapping qualities, forward strand and reverse strand supporting reads etc. along with a variety of databases such as dbSNP, 1000 genome, 5400 exome, refGene, avsift, ljb_phylop from UCSC were used to call SNVs from RNA-sequencing data. For the expressed SNVs detected in tumors, we can also identify the amino acid change and classify the protein domains using our method. RESULTS: We have applied the eSNV-Detect method to call novel expressed SNVs in two tumors (one HER2 positive and one ER positive breast tumor). BWA and Bowtie were used to align the two breast tumors. Several filtering criteria's were applied to the tumors to remove any known (germ line) or low confidence SNVs. We randomly selected 83 and 31 eSNVs for the HER2 and ER positive tumors respectively and performed Sanger sequencing validation in tumor and matched normal tissues of the tumor. We were able to validate 81 out of 83 eSNVs in the HER2 tumor sample, and 29 out of 31 eSNVs in the ER positive tumor sample. Of the eSNVs that we validated 39/81 and 10/31 were localized to functional protein domains in HER2 and ER positive tumors respectively. We are also in the process of validating eSNVs detected from RNA-seq with the tumor and blood exome sequencing data we have for the same patients. Citation Format: Xiaojia Tang, Saurabh Baheti, Khader Shameer, Kevin J. Thompson, Jean-Pierre Kocher, Edith A. Perez, Aubrey E. Thompson, Krishna R. Kalari. The eSNV-Detect: a computational system to identify expressed single nucleotide variants from tumor transcriptomes. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3185. doi:10.1158/1538-7445.AM2013-3185

Abstract Motivation: Overexpression of HER2 (the product of the ERBB2 gene) occurs in about 15% of all breast tumors. We have undertaken to use next generation transcriptome sequencing technology to identify genomic features that are unique to HER2+ tumors. Interactome mapping was then used to integrate the genes associated with these features into a transcriptome landscape model, with a view towards identifying nodes of interaction that might be targeted in HER2+ tumors. Methods: We performed 50nt paired-end RNA-sequencing for 24 breast tumors: 8 each HER2+, ER+, triple negatives (TN). In addition to breast adenocarcinomas, we also sequenced 8 early passage non-transformed HMEC cell lines as normal controls. Reads from RNA sequencing were aligned to the genome and exon junctions using TopHat software. Gene counts were summarized and annotations were performed using our in-house programs. Tukey's test was used to obtain genes or transcripts that are specific to HER2 tumor group compared to other tumors/normal. A combination of bioinformatics softwares and algorithms were used to identify SNVs. Results: Only 13527 genes with median gene count greater than 16 reads in at least one of the 4 groups were considered for differential gene expression and splicing analysis. Some 696 genes were differentially expressed in HER2+ tumors compared to ER+, TN tumors and HMECs. We identified 272 alternately spliced transcripts for which the HER2+ tumors exhibited a mean transcript expression ratio significantly different from the means of other tumor groups. We also identified 4735 expressed single nucleotide variants that are uniquely associated with HER2+ tumors compared to other tumors/normal groups. Among these 3579/4735 sequence polymorphisms were not present in the 1000 genome germline sequence database or in the dbSNP132 database of naturally occurring germline polymorphisms. Integration of all the genes obtained from genomic feature analyses (differential expression, alternative splicing, single nucleotide variance) has been carried out to indentify biological processes that are specific to the HER2+ tumor subtype. Key nodes and pathways that are specific to HER2+ tumors will be evaluated for association with clinical outcome in a large series of patients who have received HER2-targeted therapy. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 4926. doi:1538-7445.AM2012-4926

Abstract Summary: Reduced representation bisulfite sequencing (RRBS) is a cost-effective approach for genome-wide methylation pattern profiling. Analyzing RRBS sequencing data is challenging and specialized alignment/mapping programs are needed. Although such programs have been developed, a comprehensive solution that provides researchers with good quality and analyzable data is still lacking. To address this need, we have developed a Streamlined Analysis and Annotation Pipeline for RRBS data (SAAP-RRBS) that integrates read quality assessment/clean-up, alignment, methylation data extraction, annotation, reporting and visualization. This package facilitates a rapid transition from sequencing reads to a fully annotated CpG methylation report to biological interpretation. Availability and implementation: SAAP-RRBS is freely available to non-commercial users at the web site http://ndc.mayo.edu/mayo/research/biostat/stand-alone-packages.cfm. Contact: baheti.saurabh@mayo.edu or sun.zhifu@mayo.edu Supplementary Information: Supplementary data are available at Bioinformatics online.