Motivation: The Biological Reference Repository (BioR) is a toolkit for annotating variants. BioR stores public and user-specific annotation sources in indexed JSON-encoded flat files (catalogs). The BioR toolkit provides the functionality to combine and retrieve annotation from these catalogs via the command-line interface. Several catalogs from commonly used annotation sources and instructions for creating user-specific catalogs are provided. Commands from the toolkit can be combined with other UNIX commands for advanced annotation processing. We also provide instructions for the development of custom annotation pipelines. Availability and implementation: The package is implemented in Java and makes use of external tools written in Java and Perl. The toolkit can be executed on Mac OS X 10.5 and above or any Linux distribution. The BioR application, quickstart, and user guide documents and many biological examples are available at http://bioinformaticstools.mayo.edu . Contact:Kocher.JeanPierre@mayo.edu Supplementary information: Supplementary data are available at Bioinformatics online.
Chromatin immunoprecipitation (ChIP) followed by next-generation sequencing (ChIP-Seq) has been widely used to identify genomic loci of transcription factor (TF) binding and histone modifications. ChIP-Seq data analysis involves multiple steps from read mapping and peak calling to data integration and interpretation. It remains challenging and time-consuming to process large amounts of ChIP-Seq data derived from different antibodies or experimental designs using the same approach. To address this challenge, there is a need for a comprehensive analysis pipeline with flexible settings to accelerate the utilization of this powerful technology in epigenetics research. We have developed a highly integrative pipeline, termed HiChIP for systematic analysis of ChIP-Seq data. HiChIP incorporates several open source software packages selected based on internal assessments and published comparisons. It also includes a set of tools developed in-house. This workflow enables the analysis of both paired-end and single-end ChIP-Seq reads, with or without replicates for the characterization and annotation of both punctate and diffuse binding sites. The main functionality of HiChIP includes: (a) read quality checking; (b) read mapping and filtering; (c) peak calling and peak consistency analysis; and (d) result visualization. In addition, this pipeline contains modules for generating binding profiles over selected genomic features, de novo motif finding from transcription factor (TF) binding sites and functional annotation of peak associated genes. HiChIP is a comprehensive analysis pipeline that can be configured to analyze ChIP-Seq data derived from varying antibodies and experiment designs. Using public ChIP-Seq data we demonstrate that HiChIP is a fast and reliable pipeline for processing large amounts of ChIP-Seq data.
Structural variation (SV) represents a significant, yet poorly understood contribution to an individual’s genetic makeup. Advanced next-generation sequencing technologies are widely used to discover such variations, but there is no single detection tool that is considered a community standard. In an attempt to fulfil this need, we developed an algorithm, SoftSearch, for discovering structural variant breakpoints in Illumina paired-end next-generation sequencing data. SoftSearch combines multiple strategies for detecting SV including split-read, discordant read-pair, and unmated pairs. Co-localized split-reads and discordant read pairs are used to refine the breakpoints. We developed and validated SoftSearch using real and synthetic datasets. SoftSearch’s key features are 1) not requiring secondary (or exhaustive primary) alignment, 2) portability into established sequencing workflows, and 3) is applicable to any DNA-sequencing experiment (e.g. whole genome, exome, custom capture, etc.). SoftSearch identifies breakpoints from a small number of soft-clipped bases from split reads and a few discordant read-pairs which on their own would not be sufficient to make an SV call. We show that SoftSearch can identify more true SVs by combining multiple sequence features. SoftSearch was able to call clinically relevant SVs in the BRCA2 gene not reported by other tools while offering significantly improved overall performance.
11 March 2016 in Frontiers in Behavioral Neuroscience
The molecular mechanisms underlying the neuronal signaling changes in alcohol addiction and withdrawal are complex and multifaceted. The cortico-striatal circuit is highly implicated in these processes, and the striatum plays a significant role not only in the early stages of addiction, but in the developed-addictive state as well, including withdrawal symptoms. Transcriptional analysis is a useful method for determining changes in gene expression, however, the results do not always accurately correlate with protein levels. In this study, we employ label-free proteomic analysis to determine changes in protein expression within the striatum during chronic ethanol use and early withdrawal. The striatum, composed primarily of medium spiny GABAergic neurons, glutamatergic and dopaminergic nerve terminals and astrocytes, is relatively homogeneous for proteomic analysis. We were able to analyze more than 5,000 proteins from both the dorsal (caudate and putamen) and ventral (nucleus accumbens) striatum and identified significant changes following chronic intermittent ethanol exposure and acute (8 h) withdrawal compared to ethanol naïve and ethanol exposure groups respectively. Our results showed significant changes in proteins involved in glutamate and opioid peptide signaling, and also uncovered novel pathways including mitochondrial function and lipid/cholesterol metabolism, as revealed by changes in electron transport chain proteins and RXR activation pathways. These results will be useful in the development of novel treatments for alcohol withdrawal and thereby aid in recovery from alcohol use disorder.
CD36 expression in both immune and non-immune cells is known to be directly involved in cancer metastasis. Extracellular vesicles (EVs) secreted by malignant melanocytes play a vital role in developing tumor-promoting microenvironments, but it is unclear whether this is mediated through CD36. To understand the role of CD36 in melanoma, we first analyzed the SKCM dataset for clinical prognosis, evaluated the percentage of CD36 in lymphatic fluid-derived EVs (LEVs), and tested whether melanoma-derived EVs increase CD36 expression and induce M2-macrophage-like characteristics. Furthermore, we performed a multiplex immunofluorescence (MxIF) imaging analysis to evaluate the CD36 expression and its colocalization with various other cells in the lymph node (LN) of patients and control subjects. Our findings show that cutaneous melanoma patients have a worse clinical prognosis with high CD36 levels, and a higher percentage of CD36 in total LEVs were found at baseline in melanoma patients compared to control. We also found that monocytic and endothelial cells treated with melanoma EVs expressed more CD36 than untreated cells. Furthermore, melanoma-derived EVs can regulate immunosuppressive macrophage-like characteristics by upregulating CD36. The spatial imaging data show that cells in tumor-involved sentinel LNs exhibit a higher probability of CD36 expression than cells from control LNs, but this was not statistically significant. Conclusively, our findings demonstrated that CD36 plays a vital role in controlling the immunosuppressive microenvironment in the LN, which can promote the formation of a protumorigenic niche.