ProteinHistorian: tools for the comparative analysis of eukaryote protein origin.

The evolutionary history of a protein reflects the functional history of its ancestors. Recent phylogenetic studies identified distinct evolutionary signatures that characterize proteins involved in cancer, Mendelian disease, and different ontogenic stages. Despite the potential to yield insight into the cellular functions and interactions of proteins, such comparative phylogenetic analyses are rarely performed, because they require custom algorithms. We developed ProteinHistorian to make tools for performing analyses of protein origins widely available. Given a list of proteins of interest, ProteinHistorian estimates the phylogenetic age of each protein, quantifies enrichment for proteins of specific ages, and compares variation in protein age with other protein attributes. ProteinHistorian allows flexibility in the definition of protein age by including several algorithms for estimating ages from different databases of evolutionary relationships. We illustrate the use of ProteinHistorian with three example analyses. First, we demonstrate that proteins with high expression in human, compared to chimpanzee and rhesus macaque, are significantly younger than those with human-specific low expression. Next, we show that human proteins with annotated regulatory functions are significantly younger than proteins with catalytic functions. Finally, we compare protein length and age in many eukaryotic species and, as expected from previous studies, find a positive, though often weak, correlation between protein age and length. ProteinHistorian is available through a web server with an intuitive interface and as a set of command line tools; this allows biologists and bioinformaticians alike to integrate these approaches into their analysis pipelines. ProteinHistorian's modular, extensible design facilitates the integration of new datasets and algorithms. The ProteinHistorian web server, source code, and pre-computed ages for 32 eukaryotic genomes are freely available under the GNU public license at http://lighthouse.ucsf.edu/ProteinHistorian/.

Transposable element expression in tumors is associated with immune infiltration and increased antigenicity.

Profound global loss of DNA methylation is a hallmark of many cancers. One potential consequence of this is the reactivation of transposable elements (TEs) which could stimulate the immune system via cell-intrinsic antiviral responses. Here, we develop REdiscoverTE, a computational method for quantifying genome-wide TE expression in RNA sequencing data. Using The Cancer Genome Atlas database, we observe increased expression of over 400 TE subfamilies, of which 262 appear to result from a proximal loss of DNA methylation. The most recurrent TEs are among the evolutionarily youngest in the genome, predominantly expressed from intergenic loci, and associated with antiviral or DNA damage responses. Treatment of glioblastoma cells with a demethylation agent results in both increased TE expression and de novo presentation of TE-derived peptides on MHC class I molecules. Therapeutic reactivation of tumor-specific TEs may synergize with immunotherapy by inducing inflammation and the display of potentially immunogenic neoantigens.

Informatics center for mouse genomics: the dissection of complex traits of the nervous system.

In recent years, there has been an explosion in the number of tools and techniques available to researchers interested in exploring the genetic basis of all aspects of central nervous system (CNS) development and function. Here, we exploit a powerful new reductionist approach to explore the genetic basis of the very significant structural and molecular differences between the brains of different strains of mice, called either complex trait or quantitative trait loci (QTL) analysis. Our specific focus has been to provide universal access over the web to tools for the genetic dissection of complex traits of the CNS--tools that allow researchers to map genes that modulate phenotypes at a variety of levels ranging from the molecular all the way to the anatomy of the entire brain. Our website, The Mouse Brain Library (MBL; http://mbl.org) is comprised of four interrelated components that are designed to support this goal: The Brain Library, iScope, Neurocartographer, and WebQTL. The centerpiece of the MBL is an image database of histologically prepared museum-quality slides representing nearly 2000 mice from over 120 strains--a library suitable for stereologic analysis of regional volume. The iScope provides fast access to the entire slide collection using streaming video technology, enabling neuroscientists to acquire high-magnification images of any CNS region for any of the mice in the MBL. Neurocartographer provides automatic segmentation of images from the MBL by warping precisely delineated boundaries from a 3D atlas of the mouse brain. Finally, WebQTL provides statistical and graphical analysis of linkage between phenotypes and genotypes.

RNA-seq Data: Challenges in and Recommendations for Experimental Design and Analysis.

RNA-seq is widely used to determine differential expression of genes or transcripts as well as identify novel transcripts, identify allele-specific expression, and precisely measure translation of transcripts. Thoughtful experimental design and choice of analysis tools are critical to ensure high-quality data and interpretable results. Important considerations for experimental design include number of replicates, whether to collect paired-end or single-end reads, sequence length, and sequencing depth. Common analysis steps in all RNA-seq experiments include quality control, read alignment, assigning reads to genes or transcripts, and estimating gene or transcript abundance. Our aims are two-fold: to make recommendations for common components of experimental design and assess tool capabilities for each of these steps. We also test tools designed to detect differential expression, since this is the most widespread application of RNA-seq. We hope that these analyses will help guide those who are new to RNA-seq and will generate discussion about remaining needs for tool improvement and development.

The synthetic genetic interaction network reveals small molecules that target specific pathways in Sacchromyces cerevisiae.

High-throughput elucidation of synthetic genetic interactions (SGIs) has contributed to a systems-level understanding of genetic robustness and fault-tolerance encoded in the genome. Pathway targets of various compounds have been predicted by comparing chemical-genetic synthetic interactions to a network of SGIs. We demonstrate that the SGI network can also be used in a powerful reverse pathway-to-drug approach for identifying compounds that target specific pathways of interest. Using the SGI network, the method identifies an indicator gene that may serve as a good candidate for screening a library of compounds. The indicator gene is selected so that compounds found to produce sensitivity in mutants deleted for the indicator gene are likely to abrogate the target pathway. We tested the utility of the SGI network for pathway-to-drug discovery using the DNA damage checkpoint as the target pathway. An analysis of the compendium of synthetic lethal interactions in yeast showed that superoxide dismutase 1 (SOD1) has significant SGI connectivity with a large subset of DNA damage checkpoint and repair (DDCR) genes in Saccharomyces cerevisiae, and minimal SGIs with non-DDCR genes. We screened a sod1Δ strain against three National Cancer Institute (NCI) compound libraries using a soft agar high-throughput halo assay. Fifteen compounds out of ∼3100 screened showed selective toxicity toward sod1Δ relative to the isogenic wild type (wt) strain. One of these, 1A08, caused a transient increase in growth in the presence of sublethal doses of DNA damaging agents, suggesting that 1A08 inhibits DDCR signaling in yeast. Genome-wide screening of 1A08 against the library of viable homozygous deletion mutants further supported DDCR as the relevant targeted pathway of 1A08. When assayed in human HCT-116 colorectal cancer cells, 1A08 caused DNA-damage resistant DNA synthesis and blocked the DNA-damage checkpoint selectively in S-phase.

GenomeMixer: a complex genetic cross simulator.

UNLABELLED: GenomeMixer is a cross-platform application that simulates meiotic recombination events for large and complex multigenerational genetic crosses among sexually reproducing diploid species and outputs simulated progeny to several standard mapping programs. AVAILABILITY: Documentation, C++ source, and binaries for Mac OS X and x86 Linux are freely available at http://www.nervenet.org/genome_mixer/. GenomeMixer can be compiled on any system with support for the Trolltech Qt toolkit, including Windows.