Comparing reference-based RNA-Seq mapping methods for non-human primate data.

BACKGROUND: The application of next-generation sequencing technology to gene expression quantification analysis, namely, RNA-Sequencing, has transformed the way in which gene expression studies are conducted and analyzed. These advances are of particular interest to researchers studying organisms with missing or incomplete genomes, as the need for knowledge of sequence information is overcome. De novo assembly methods have gained widespread acceptance in the RNA-Seq community for organisms with no true reference genome or transcriptome. While such methods have tremendous utility, computational cost is still a significant challenge for organisms with large and complex genomes. RESULTS: In this manuscript, we present a comparison of four reference-based mapping methods for non-human primate data. We utilize TopHat2 and GSNAP for mapping to the human genome, and Bowtie2 and Stampy for mapping to the human genome and transcriptome for a total of six mapping approaches. For each of these methods, we explore mapping rates and locations, number of detected genes, correlations between computed expression values, and the utility of the resulting data for differential expression analysis. CONCLUSIONS: We show that reference-based mapping methods indeed have utility in RNA-Seq analysis of mammalian data with no true reference, and the details of mapping methods should be carefully considered when doing so. Critical algorithm features include short seed sequences, the allowance of mismatches, and the allowance of gapped alignments in addition to splice junction gaps. Such features facilitate sensitive alignment of non-human primate RNA-Seq data to a human reference.

Development of a novel preclinical model of pneumococcal pneumonia in nonhuman primates.

Pneumococcal pneumonia is a leading cause of bacterial infection and death worldwide. Current diagnostic tests for detecting Streptococcus pneumoniae can be unreliable and can mislead clinical decision-making and treatment. To address this concern, we developed a preclinical model of pneumococcal pneumonia in nonhuman primates useful for identifying novel biomarkers, diagnostic tests, and therapies for human S. pneumoniae infection. Adult colony-bred baboons (n = 15) were infected with escalating doses of S. pneumoniae (Serotype 19A-7). We characterized the pathophysiological and serological profiles of healthy and infected animals over 7 days. Pneumonia was prospectively defined by the presence of three criteria: (1) change in white blood cell count, (2) isolation of S. pneumoniae from bronchoalveolar lavage fluid (BALF) or blood, and (3) concurrent signs/symptoms of infection. Animals given 10(9) CFU consistently met our definition and developed a phenotype of tachypnea, tachycardia, fever, hypoxemia, and radiographic lobar infiltrates at 48 hours. BALF and plasma cytokines, including granulocyte colony-stimulating factor, IL-6, IL-10, and IL-1ra, peaked at 24 to 48 hours. At necropsy, there was lobar consolidation with frequent pleural involvement. Lung histopathology showed alveolar edema and macrophage influx in areas of organizing pneumonia. Hierarchical clustering of peripheral blood RNA data at 48 hours correctly identified animals with and without pneumonia. Dose-dependent inoculation of baboons with S. pneumoniae produces a host response ranging from spontaneous clearance (10(6) CFU) to severe pneumonia (10(9) CFU). Selected BALF and plasma cytokine levels and RNA profiles were associated with severe pneumonia and may provide clinically useful parameters after validation.

A flexible statistical model for alignment of label-free proteomics data--incorporating ion mobility and product ion information.

BACKGROUND: The goal of many proteomics experiments is to determine the abundance of proteins in biological samples, and the variation thereof in various physiological conditions. High-throughput quantitative proteomics, specifically label-free LC-MS/MS, allows rapid measurement of thousands of proteins, enabling large-scale studies of various biological systems. Prior to analyzing these information-rich datasets, raw data must undergo several computational processing steps. We present a method to address one of the essential steps in proteomics data processing--the matching of peptide measurements across samples. RESULTS: We describe a novel method for label-free proteomics data alignment with the ability to incorporate previously unused aspects of the data, particularly ion mobility drift times and product ion information. We compare the results of our alignment method to PEPPeR and OpenMS, and compare alignment accuracy achieved by different versions of our method utilizing various data characteristics. Our method results in increased match recall rates and similar or improved mismatch rates compared to PEPPeR and OpenMS feature-based alignment. We also show that the inclusion of drift time and product ion information results in higher recall rates and more confident matches, without increases in error rates. CONCLUSIONS: Based on the results presented here, we argue that the incorporation of ion mobility drift time and product ion information are worthy pursuits. Alignment methods should be flexible enough to utilize all available data, particularly with recent advancements in experimental separation methods.

Transcription initiation patterns indicate divergent strategies for gene regulation at the chromatin level.

The application of deep sequencing to map 5' capped transcripts has confirmed the existence of at least two distinct promoter classes in metazoans: "focused" promoters with transcription start sites (TSSs) that occur in a narrowly defined genomic span and "dispersed" promoters with TSSs that are spread over a larger window. Previous studies have explored the presence of genomic features, such as CpG islands and sequence motifs, in these promoter classes, but virtually no studies have directly investigated the relationship with chromatin features. Here, we show that promoter classes are significantly differentiated by nucleosome organization and chromatin structure. Dispersed promoters display higher associations with well-positioned nucleosomes downstream of the TSS and a more clearly defined nucleosome free region upstream, while focused promoters have a less organized nucleosome structure, yet higher presence of RNA polymerase II. These differences extend to histone variants (H2A.Z) and marks (H3K4 methylation), as well as insulator binding (such as CTCF), independent of the expression levels of affected genes. Notably, differences are conserved across mammals and flies, and they provide for a clearer separation of promoter architectures than the presence and absence of CpG islands or the occurrence of stalled RNA polymerase. Computational models support the stronger contribution of chromatin features to the definition of dispersed promoters compared to focused start sites. Our results show that promoter classes defined from 5' capped transcripts not only reflect differences in the initiation process at the core promoter but also are indicative of divergent transcriptional programs established within gene-proximal nucleosome organization.

Gene by sex interaction for measures of obesity in the framingham heart study.

Obesity is an increasingly prevalent and severe health concern with a substantial heritable component and marked sex differences. We sought to determine if the effect of genetic variants also differed by sex by performing a genome-wide association study modeling the effect of genotype-by-sex interaction on obesity phenotypes. Genotype data from individuals in the Framingham Heart Study Offspring cohort were analyzed across five exams. Although no variants showed genome-wide significant gene-by-sex interaction in any individual exam, four polymorphisms displayed a consistent BMI association (P-values .00186 to .00010) across all five exams. These variants were clustered downstream of LYPLAL1, which encodes a lipase/esterase expressed in adipose tissue, a locus previously identified as having sex-specific effects on central obesity. Primary effects in males were in the opposite direction from females and were replicated in Framingham Generation 3. Our data support a sex-influenced association between genetic variation at the LYPLAL1 locus and obesity-related traits.