Search

BackgroundSingle-cell RNA-seq (scRNA-seq) profiling has revealed remarkable variation in transcription, suggesting that expression of many genes at the single-cell level is intrinsically stochastic and noisy. Yet, on the cell population level, a subset of genes traditionally referred to as housekeeping genes (HKGs) are found to be stably expressed in different cell and tissue types. It is therefore critical to question whether stably expressed genes (SEGs) can be identified on the single-cell level, and if so, how can their expression stability be assessed? We have previously proposed a computational framework for ranking expression stability of genes in single cells for scRNA-seq data normalization and integration. In this study, we perform detailed evaluation and characterization of SEGs derived from this framework.ResultsHere, we show that gene expression stability indices derived from the early human and mouse development scRNA-seq datasets and the "Mouse Atlas" dataset are reproducible and conserved across species. We demonstrate that SEGs identified from single cells based on their stability indices are considerably more stable than HKGs defined previously from cell populations across diverse biological systems. Our analyses indicate that SEGs are inherently more stable at the single-cell level and their characteristics reminiscent of HKGs, suggesting their potential role in sustaining essential functions in individual cells.ConclusionsSEGs identified in this study have immediate utility both for understanding variation and stability of single-cell transcriptomes and for practical applications such as scRNA-seq data normalization. Our framework for calculating gene stability index, "scSEGIndex," is incorporated into the scMerge Bioconductor R package (https://sydneybiox.github.io/scMerge/reference/scSEGIndex.html) and can be used for identifying genes with stable expression in scRNA-seq datasets.

Background The differentiation and maturation trajectories of fetal liver stem/progenitor cells (LSPCs) are not fully understood at single-cell resolution, and a priori knowledge of limited biomarkers could restrict trajectory tracking. Results We employed marker-free single-cell RNA-Seq to characterize comprehensive transcriptional profiles of 507 cells randomly selected from seven stages between embryonic day 11.5 and postnatal day 2.5 during mouse liver development, and also 52 Epcam-positive cholangiocytes from postnatal day 3.25 mouse livers. LSPCs in developing mouse livers were identified via marker-free transcriptomic profiling. Single-cell resolution dynamic developmental trajectories of LSPCs exhibited contiguous but discrete genetic control through transcription factors and signaling pathways. The gene expression profiles of cholangiocytes were more close to that of embryonic day 11.5 rather than other later staged LSPCs, cuing the fate decision stage of LSPCs. Our marker-free approach also allows systematic assessment and prediction of isolation biomarkers for LSPCs. Conclusions Our data provide not only a valuable resource but also novel insights into the fate decision and transcriptional control of self-renewal, differentiation and maturation of LSPCs. Electronic supplementary material The online version of this article (10.1186/s12864-017-4342-x) contains supplementary material, which is available to authorized users.

Learning from positive and unlabeled data frequently occurs in applications where only a subset of positive instances is available while the rest of the data are unlabeled. In such scenarios, often the goal is to create a discriminant model that can accurately classify both positive and negative data by modelling from labeled and unlabeled instances. In this study, we propose an adaptive sampling (AdaSampling) approach that utilises prediction probabilities from a model to iteratively update the training data. Starting with equal prior probabilities for all unlabeled data, our method "wraps" around a predictive model to iteratively update these probabilities to distinguish positive and negative instances in unlabeled data. Subsequently, one or more robust negative set(s) can be drawn from unlabeled data, according to the likelihood of each instance being negative, to train a single classification model or ensemble of models.

Tandem mass spectrometry is one of the most popular techniques for quantitation of proteomes. There exists a large variety of options in each stage of data preprocessing that impact the bias and variance of the summarized protein-level values. Using a newly released data set satisfying a replicated Latin squares design, a diverse set of performance metrics has been developed and implemented in a web-based application, Quantitative Performance Evaluator for Proteomics (QPEP). QPEP has the flexibility to allow users to apply their own method to preprocess this data set and share the results, allowing direct and straightforward comparison of new methodologies. Application of these new metrics to three case studies highlights that (i) the summarization of peptides to proteins is robust to the choice of peptide summary used, (ii) the differences between iTRAQ labels are stronger than the differences between experimental runs, and (iii) the commercial software ProteinPilot performs equivalently well at between-sample normalization to more complicated methods developed by academics. Importantly, finding (ii) underscores the benefits of using the principles of randomization and blocking to avoid the experimental measurements being confounded by technical factors. Data are available via ProteomeXchange with identifier PXD003608.

Protein post-translational modifications (PTMs) are crucial for signal transduction in cells. In order to understand key cell signaling events, identification of functionally important PTMs, which are more likely to be evolutionarily conserved, is necessary. In recent times, high-throughput mass spectrometry (MS) has made quantitative datasets in diverse species readily available, which has led to a growing need for tools to facilitate cross-species comparison of PTM data. Cross-species comparison of PTM sites is difficult since they often lie in structurally disordered protein domains. Current tools that address this can only map known PTMs between species based on previously annotated orthologous phosphosites and do not enable cross-species mapping of newly identified modification sites. Here, we describe an automated web-based tool, PhosphOrtholog, that accurately maps annotated and novel orthologous PTM sites from high-throughput MS-based experimental data obtained from different species without relying on existing PTM databases. Identification of conserved PTMs across species from large-scale experimental data increases our knowledgebase of evolutionarily conserved and functional PTM sites that influence most biological processes. In this Chapter, we illustrate with examples how to use PhosphOrtholog to map novel PTM sites from cross-species MS-based phosphoproteomics data.

Protein post-translational modifications (PTMs) are crucial for signal transduction in cells. In order to understand key cell signaling events, identification of functionally important PTMs, which are more likely to be evolutionarily conserved, is necessary. In recent times, high-throughput mass spectrometry (MS) has made quantitative datasets in diverse species readily available, which has led to a growing need for tools to facilitate cross-species comparison of PTM data. Cross-species comparison of PTM sites is difficult since they often lie in structurally disordered protein domains. Current tools that address this can only map known PTMs between species based on previously annotated orthologous phosphosites and do not enable cross-species mapping of newly identified modification sites. Here, we describe an automated web-based tool, PhosphOrtholog, that accurately maps annotated and novel orthologous PTM sites from high-throughput MS-based experimental data obtained from different species without relying on existing PTM databases. Identification of conserved PTMs across species from large-scale experimental data increases our knowledgebase of evolutionarily conserved and functional PTM sites that influence most biological processes. In this Chapter, we illustrate with examples how to use PhosphOrtholog to map novel PTM sites from cross-species MS-based phosphoproteomics data.

Heritable trait variation within a population of organisms is largely governed by DNA variations that impact gene transcription and protein function. Identifying genetic variants that affect complex functional traits is a primary aim of population genetics studies, especially in the context of human disease and agricultural production traits. The identification of alleles directly altering mRNA expression and thereby biological function is challenging due to difficulty in isolating direct effects of cis-acting genetic variations from indirect trans-acting genetic effects. Allele specific gene expression or allelic imbalance in gene expression (AI) occurring at heterozygous loci provides an opportunity to identify genes directly impacted by cis-acting genetic variants as indirect trans-acting effects equally impact the expression of both alleles. However, the identification of genes showing AI in the context of the expression of all genes remains a challenge due to a variety of technical and statistical issues. The current study focuses on the discovery of genes showing AI using single nucleotide polymorphisms as allelic reporters. By developing a computational and statistical process that addressed multiple analytical challenges, we ranked 5,809 genes for evidence of AI using RNA-Seq data derived from brown adipose tissue samples from a cohort of late gestation fetal lambs and then identified a conservative subgroup of 1,293 genes. Thus, AI was extensive, representing approximately 25% of the tested genes. Genes associated with AI were enriched for multiple Gene Ontology (GO) terms relating to lipid metabolism, mitochondrial function and the extracellular matrix. These functions suggest that cis-acting genetic variations causing AI in the population are preferentially impacting genes involved in energy homeostasis and tissue remodelling. These functions may contribute to production traits likely to be under genetic selection in the population. Refereed/Peer-reviewed

The increasing role played by liquid chromatography-mass spectrometry (LC-MS)-based proteomics in biological discovery has led to a growing need for quality control (QC) on the LC-MS systems. While numerous quality control tools have been developed to track the performance of LC-MS systems based on a pre-defined set of performance factors (e.g., mass error, retention time), the precise influence and contribution of the performance factors and their generalization property to different biological samples are not as well characterized. Here, a web-based application (QCMAP) is developed for interactive diagnosis and prediction of the performance of LC-MS systems across different biological sample types. Leveraging on a standardized HeLa cell sample run as QC within a multi-user facility, predictive models are trained on a panel of commonly used performance factors to pinpoint the precise conditions to a (un)satisfactory performance in three LC-MS systems. It is demonstrated that the learned model can be applied to predict LC-MS system performance for brain samples generated from an independent study. By compiling these predictive models into our web-application, QCMAP allows users to benchmark the performance of their LC-MS systems using their own samples and identify key factors for instrument optimization. QCMAP is freely available from: http://shiny.maths.usyd.edu.au/QCMAP/.

Background The molecular mechanisms underlying the development of the unusual echinoderm pentameral body plan and their likeness to mechanisms underlying the development of the bilateral plans of other deuterostomes are of interest in tracing body plan evolution. In this first study of the spatial expression of genes associated with Nodal and BMP2/4 signalling during the transition to pentamery in sea urchins, we investigate Heliocidaris erythrogramma, a species that provides access to the developing adult rudiment within days of fertilization. Results BMP2/4, and the putative downstream genes, Six1/2, Eya, Tbx2/3 and Msx were expressed in the earliest morphological manifestation of pentamery during development, the five hydrocoele lobes. The formation of the vestibular ectoderm, the specialized region overlying the left coelom that forms adult ectoderm, involved the expression of putative Nodal target genes Chordin, Gsc and BMP2/4 and putative BMP2/4 target genes Dlx, Msx and Tbx. The expression of Nodal, Lefty and Pitx2 in the right ectoderm, and Pitx2 in the right coelom, was as previously observed in other sea urchins. Conclusion That genes associated with Nodal and BMP2/4 signalling are expressed in the hydrocoele lobes, indicates that they have a role in the developmental transition to pentamery, contributing to our understanding of how the most unusual body plan in the Bilateria may have evolved. We suggest that the Nodal and BMP2/4 signalling cascades might have been duplicated or split during the evolution to pentamery. Electronic supplementary material The online version of this article (doi:10.1186/s12861-017-0145-1) contains supplementary material, which is available to authorized users.

A consistent difference in average expression level, often referred to as differential expression (DE), has long been used to identify genes useful for classification. However, recent cancer studies have shown that when transcription factors or epigenetic signals become deregulated, a change in expression variability (DV) of target genes is frequently observed. This suggests that assessing the importance of genes by either differential expression or variability alone potentially misses sets of important biomarkers that could lead to improved predictions and treatments. Here, we describe a new approach for assessing the importance of genes based on differential distribution (DD), which combines information from differential expression and differential variability into a unified metric. We show that feature ranking and selection stability based on DD can perform two to three times better than DE or DV alone, and that DD yields equivalent error rates to DE and DV. Finally, assessing genes via differential distribution produces a complementary set of selected genes to DE and DV, potentially opening up new categories of biomarkers.

Mass spectrometry (MS)-based quantitative phosphoproteomics has become a key approach for proteome-wide profiling of phosphorylation in tissues and cells. Traditional experimental design often compares a single treatment with a control, whereas increasingly more experiments are designed to compare multiple treatments with respect to a control. To this end, the development of bioinformatic tools that can integrate multiple treatments and visualise kinases and substrates under combinatorial perturbations is vital for dissecting concordant and/or independent effects of each treatment. Here, we propose a hypothesis driven kinase perturbation analysis (KinasePA) to annotate and visualise kinases and their substrates that are perturbed by various combinatorial effects of treatments in phosphoproteomics experiments. We demonstrate the utility of KinasePA through its application to two large-scale phosphoproteomics datasets and show its effectiveness in dissecting kinases and substrates within signalling pathways driven by unique combinations of cellular stimuli and inhibitors. We implemented and incorporated KinasePA as part of the "directPA" R package available from the comprehensive R archive network (CRAN). Furthermore, KinasePA also has an interactive web interface that can be readily applied to annotate user provided phosphoproteomics data (http://kinasepa.pengyiyang.org).

BackgroundData made available through large cancer consortia like The Cancer Genome Atlas make for a rich source of information to be studied across and between cancers. In recent years, network approaches have been applied to such data in uncovering the complex interrelationships between mutational and expression profiles, but lack direct testing for expression changes via mutation. In this pan-cancer study we analyze mutation and gene expression information in an integrative manner by considering the networks generated by testing for differences in expression in direct association with specific mutations. We relate our findings among the 19 cancers examined to identify commonalities and differences as well as their characteristics. ResultsUsing somatic mutation and gene expression information across 19 cancers, we generated mutation-expression networks per cancer. On evaluation we found that our generated networks were significantly enriched for known cancer-related genes, such as skin cutaneous melanoma (p

A key step in the analysis of mass spectrometry (MS)-based proteomics data is the inference of proteins from identified peptide sequences. Here we describe Re-Fraction, a novel machine learning algorithm that enhances deterministic protein identification. Re-Fraction utilizes several protein physical properties to assign proteins to expected protein fractions that comprise large-scale MS-based proteomics data. This information is then used to appropriately assign peptides to specific proteins. This approach is sensitive, highly specific, and computationally efficient. We provide algorithms and source code for the current version of Re-Fraction, which accepts output tables from the MaxQuant environment. Nevertheless, the principles behind Re-Fraction can be applied to other protein identification pipelines where data are generated from samples fractionated at the protein level. We demonstrate the utility of this approach through reanalysis of data from a previously published study and generate lists of proteins deterministically identified by Re-Fraction that were previously only identified as members of a protein group. We find that this approach is particularly useful in resolving protein groups composed of splice variants and homologues, which are frequently expressed in a cell- or tissue-specific manner and may have important biological consequences.

Motivation: Mass spectrometry (MS)-based proteomics is one of the most commonly used research techniques for identifying and characterizing proteins in biological and medical research. The identification of a protein is the critical first step in elucidating its biological function. Successful protein identification depends on various interrelated factors, including effective analysis of MS data generated in a proteomic experiment. This analysis comprises several stages, often combined in a pipeline or workflow. The first component of the analysis is known as spectra pre-processing. In this component, the raw data generated by the mass spectrometer is processed to eliminate noise and identify the mass-to-charge ratio (m/z) and intensity for the peaks in the spectrum corresponding to the presence of certain peptides or peptide fragments. Since all downstream analyses depend on the pre-processed data, effective pre-processing is critical to protein identification and characterization. There is a critical need for more robust pre-processing algorithms that perform well on tandem mass spectra under a variety of different conditions and can be easily integrated into sophisticated data analysis pipelines for practical wet-lab applications. Result: We have developed a new pre-processing algorithm. Based on wavelet theory, our method uses a dynamic peak model to identify peaks. It is designed to be easily integrated into a complete proteomic analysis workflow. We compared the method with other available algorithms using a reference library of raw MS and tandem MS spectra with known protein composition information. Our pre-processing algorithm results in the identification of significantly more peptides and proteins in the downstream analysis for a given false discovery rate. Availability: Software available at: http://www.maths.usyd.edu.au/u/penghao/index.html Contact: penghao.wang@sydney.edu.au

Human embryonic stem cell (hESC) lines are derived from the inner cell mass (ICM) of preimplantation human blastocysts obtained on days 5-6 following fertilization. Based on their derivation, they were once thought to be the equivalent of the ICM. Recently, however, studies in mice reported the derivation of mouse embryonic stem cell lines from the epiblast; these epiblast lines bear significant resemblance to human embryonic stem cell lines in terms of culture, differentiation potential and gene expression. In this study, we compared gene expression in human ICM cells isolated from the blastocyst and embryonic stem cells. We demonstrate that expression profiles of ICM clusters from single embryos and hESC populations were highly reproducible. Moreover, comparison of global gene expression between individual ICM clusters and human embryonic stem cells indicated that these two cell types are significantly different in regards to gene expression, with fewer than one half of all genes expressed in both cell types. Genes of the isolated human inner cell mass that are upregulated and downregulated are involved in numerous cellular pathways and processes; a subset of these genes may impart unique characteristics to hESCs such as proliferative and self-renewal properties.

Background Most biological processes are influenced by protein post-translational modifications (PTMs). Identifying novel PTM sites in different organisms, including humans and model organisms, has expedited our understanding of key signal transduction mechanisms. However, with increasing availability of deep, quantitative datasets in diverse species, there is a growing need for tools to facilitate cross-species comparison of PTM data. This is particularly important because functionally important modification sites are more likely to be evolutionarily conserved; yet cross-species comparison of PTMs is difficult since they often lie in structurally disordered protein domains. Current tools that address this can only map known PTMs between species based on known orthologous phosphosites, and do not enable the cross-species mapping of newly identified modification sites. Here, we addressed this by developing a web-based software tool, PhosphOrtholog (www.phosphortholog.com) that accurately maps protein modification sites between different species. This facilitates the comparison of datasets derived from multiple species, and should be a valuable tool for the proteomics community. Results Here we describe PhosphOrtholog, a web-based application for mapping known and novel orthologous PTM sites from experimental data obtained from different species. PhosphOrtholog is the only generic and automated tool that enables cross-species comparison of large-scale PTM datasets without relying on existing PTM databases. This is achieved through pairwise sequence alignment of orthologous protein residues. To demonstrate its utility we apply it to two sets of human and rat muscle phosphoproteomes generated following insulin and exercise stimulation, respectively, and one publicly available mouse phosphoproteome following cellular stress revealing high mapping and coverage efficiency. Although coverage statistics are dataset dependent, PhosphOrtholog increased the number of cross-species mapped sites in all our example data sets by more than double when compared to those recovered using existing resources such as PhosphoSitePlus. Conclusions PhosphOrtholog is the first tool that enables mapping of thousands of novel and known protein phosphorylation sites across species, accessible through an easy-to-use web interface. Identification of conserved PTMs across species from large-scale experimental data increases our knowledgebase of functional PTM sites. Moreover, PhosphOrtholog is generic being applicable to other PTM datasets such as acetylation, ubiquitination and methylation. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1820-x) contains supplementary material, which is available to authorized users.

While adipogenesis is controlled by a cascade of transcription factors, the global gene expression profiles in the early phase of adipogenesis are not well defined. Using microarray analysis of gene expression in 3T3-L1 cells, we have identified evidence for the activity of 2,568 genes during the early phase of adipocyte differentiation. One of these, the ISL1 gene, was of interest since its expression was markedly upregulated 1 h after initiation of differentiation, with a subsequent rapid decline. Overexpression of ISL1 at early times during adipocyte differentiation but not at later times was found to profoundly inhibit differentiation. This was accompanied by moderate downregulation of peroxisome proliferator-activated receptor γ (PPARγ) levels, substantial downregulation of PPARγ downstream genes, and downregulation of bone morphogenetic protein 4 (BMP4) levels in preadipocytes. Readdition of BMP4 overcame the inhibitory effect of ISL1 on the expression of PPARγ but not aP2, a gene downstream of PPARγ, and BMP4 also partially rescued ISL1 inhibition of adipogenesis, an effect which is additive with rosiglitazone. These results suggest that ISL1 is intimately involved in early regulation of adipogenesis, modulating PPARγ expression and activity via BMP4-dependent and -independent mechanisms. Our time course gene expression survey sets the stage for further studies to explore other early and immediate regulators.

In patients with metastatic melanoma, the identification and validation of accurate prognostic biomarkers will assist rational treatment planning. Studies based on "-omics" technologies have focussed on a single high-throughput data type such as gene or microRNA transcripts. Occasionally, these features have been evaluated in conjunction with limited clinico-pathologic data. With the increased availability of multiple data types, there is a pressing need to tease apart which of these sources contain the most valuable prognostic information. We evaluated and integrated several data types derived from the same tumor specimens in AJCC stage III melanoma patients-gene, protein, and microRNA expression as well as clinical, pathologic and mutation information-to determine their relative impact on prognosis. We used classification frameworks based on pre-validation and bootstrap multiple imputation to compare the prognostic power of each data source, both individually as well as integratively. We found that the prognostic utility of clinico-pathologic information was not out-performed by any of the various "-omics" platforms. Rather, a combination of clinico-pathologic variables and mRNA expression data performed best. Furthermore, a patient-based classification analysis revealed that the prognostic accuracy of various data types was not the same for different patients. This indicates that ongoing development in the individualized evaluation of melanoma patients must take account of the value of both traditional and novel "-omics" measurements.

Background The large-scale sequencing of 5' cap enriched cDNA promises to reveal the diversity of transcription initiation across entire genomes. The process of transcription is noisy, and there is often no single, exact start site. This creates the need for a fast and simple method of identifying transcription start peaks based on this type of data. Due to both biological and technical noise, many of the peaks seen are not real transcription initiation events. Classification of the observed peaks is an essential filtering step in the discovery of genuine initiation locations. Results We develop a two-stage approach consisting of a fast and simple algorithm based on a sliding window with Poisson null distribution for detecting the genomic locations of peaks, followed by a linear support vector machine classifier to distinguish between peaks which represent the initiation of transcription and peaks that do not. Comparison of classification performance to the best existing method based on whole genome segmentation showed comparable precision and improved recall. Internal features, which are intrinsic to the data and require no further experiments, had high precision and recall rates. Addition of pooled external data or matched RNA sequencing data resulted in gains of recall with equivalent precision. Conclusions The Poisson sliding window model is an effective and fast way of taking the peak neighbourhood into account, and finding statistically significant peaks over a range of transcript expression values. It is orders of magnitude faster than doing whole genome segmentation. The support vector classification scheme has better precision and recall than existing methods. Integrating additional datasets is shown to provide minor gains in recall, in comparison to using only the cap-sequencing data.

Motivation: Mass spectrometry-based iTRAQ protein quantification is a high-throughput assay for determining relative protein expressions and identifying disease biomarkers. Processing and analysis of these large and complex data involves a number of distinct components and it is desirable to have a pipeline to efficiently integrate these together. To date, there are limited public available comprehensive analysis pipelines for iTRAQ data and many of these existing pipelines have limited visualization tools and no convenient interfaces with downstream analyses. We have developed a new open source comprehensive iTRAQ analysis pipeline, OCAP, integrating a wavelet-based preprocessing algorithm which provides better peak picking, a new quantification algorithm and a suite of visualizsation tools. OCAP is mainly developed in C++ and is provided as a standalone version (OCAP_standalone) as well as an R package. The R package (OCAP) provides the necessary interfaces with downstream statistical analysis. Availability: OCAP is freely available and can be downloaded at http://www.maths.usyd.edu.au/u/penghao Contact: penghao.wang@sydney.edu.au

The rapid increase in the use of microarray studies has generated many questions on how to plan and design experiments that will effectively utilize this technology. Investigators often require answers to questions relating to microarray platforms, RNA samples, options for replication, allocation of samples to arrays, sample sizes, appropriate downstream analysis, and many others. Careful consideration of these issues is critical to ensure the efficiency and reliability of the actual microarray experiments, and will assist in enhancing interpretability of the experimental results.

The rapid increase in the use of microarray studies has generated many questions on how to plan and design experiments that will effectively utilize this technology. Investigators often require answers to questions relating to microarray platforms, RNA samples, options for replication, allocation of samples to arrays, sample sizes, appropriate downstream analysis, and many others. Careful consideration of these issues is critical to ensure the efficiency and reliability of the actual microarray experiments, and will assist in enhancing interpretability of the experimental results.

Microarray technology is widely employed for studying the molecular mechanisms underlying complex diseases. However, analyses of individual diseases or models of diseases frequently yield extensive lists of differentially expressed genes with uncertain relationships to disease pathogenesis.To compare gene expression changes in a heterogeneous set of lung disease models in order to identify common gene expression changes seen in diverse forms of lung pathology, as well as relatively small subsets of genes likely to be involved in specific pathophysiological processes.We profiled lung gene expression in 12 mouse models of infection, allergy, and lung injury. A linear model was used to estimate transcript expression changes for each model, and hierarchical clustering was used to compare expression patterns between models. Selected expression changes were verified by quantitative polymerase chain reaction.A total of 24 transcripts, including many involved in inflammation and immune activation, were differentially expressed in a substantial majority (9 or more) of the models. Expression patterns distinguished three groups of models: (1) bacterial infection (n = 5), with changes in 89 transcripts, including many related to nuclear factor-kappaB signaling, cytokines, chemokines, and their receptors; (2) bleomycin-induced diseases (n = 2), with changes in 53 transcripts, including many related to matrix remodeling and Wnt signaling; and (3) T helper cell type 2 (allergic) inflammation (n = 5), with changes in 26 transcripts, including many encoding epithelial secreted molecules, ion channels, and transporters.This multimodel dataset highlights novel genes likely involved in various pathophysiological processes and will be a valuable resource for the investigation of molecular mechanisms underlying lung disease pathogenesis.

Insulin resistance (IR) is one of the earliest predictors of type 2 diabetes. However, diagnosis of IR is limited. High fat fed mouse models provide key insights into IR. We hypothesized that early features of IR are associated with persistent changes in gene expression (GE) and endeavored to (a) develop novel methods for improving signal:noise in analysis of human GE using mouse models; (b) identify a GE motif that accurately diagnoses IR in humans; and (c) identify novel biology associated with IR in humans. We integrated human muscle GE data with longitudinal mouse GE data and developed an unbiased three-level cross-species analysis platform (single gene, gene set, and networks) to generate a gene expression motif (GEM) indicative of IR. A logistic regression classification model validated GEM in three independent human data sets (n=115). This GEM of 93 genes substantially improved diagnosis of IR compared with routine clinical measures across multiple independent data sets. Individuals misclassified by GEM possessed other metabolic features raising the possibility that they represent a separate metabolic subclass. The GEM was enriched in pathways previously implicated in insulin action and revealed novel associations between β-catenin and Jak1 and IR. Functional analyses using small molecule inhibitors showed an important role for these proteins in insulin action. This study shows that systems approaches for identifying molecular signatures provides a powerful way to stratify individuals into discrete metabolic groups. Moreover, we speculate that the β-catenin pathway may represent a novel biomarker for IR in humans that warrant future investigation. A new molecular profile for insulin resistance (IR) has been identified by integrating gene expression (GE) profiles from two different species. Prof James from the University of Sydney and colleagues from Australia and the US sought an easier way to diagnose IR–one of the earliest predictors of type 2 diabetes. Generating a generic GE-based signature of the condition is difficult because, though gene activity changes, they vary among individuals. To solve this problem, James’s team combined GE patterns from mice fed high-fat diets and humans to generate a signature of IR including 93 different genes. When examined in 115 individuals from three different studies, this gene collection was able to better discriminate between insulin resistant and insulin sensitive patients compared to standard measures. Further investigation of this geneset could prove useful in the clinic.

Background The most widely used amplification method for microarray analysis of gene expression uses T7 RNA polymerase-driven in vitro transcription (IVT) to produce complementary RNA (cRNA) that can be hybridized to arrays. However, multiple rounds of amplification are required when assaying very small amounts of starting RNA. Moreover, certain cRNA-DNA mismatches are more stable than the analogous cDNA-DNA mismatches and this might increase non-specific hybridization. We sought to determine whether a recently developed linear isothermal amplification method (ribo-SPIA) that produces single stranded cDNA would offer advantages over traditional IVT-based methods for microarray-based analyses of transcript expression. Results A single round of ribo-SPIA amplification produced sufficient sscDNA for hybridizations when as little as 5 ng of starting total RNA was used. Comparisons of probe set signal intensities obtained from replicate amplifications showed consistently high correlations (r = 0.99). We compared gene expression in two different human RNA samples using ribo-SPIA. Compared with one round IVT, ribo-SPIA had a larger dynamic range and correlated better with quantitative PCR results even though we used 1000-fold less starting RNA. The improved dynamic range was associated with decreases in hybridization to mismatch control probes. Conclusion The use of amplified sscDNA may offer substantial advantages over IVT-based amplification methods, especially when very limited amounts of starting RNA are available. The use of sscDNA targets instead of cRNA targets appears to improve hybridization specificity.

Background Large scale single cell transcriptome profiling has exploded in recent years and has enabled unprecedented insight into the behavior of individual cells. Identifying genes with high levels of expression using data from single cell RNA sequencing can be useful to characterize very active genes and cells in which this occurs. In particular single cell RNA-Seq allows for cell-specific characterization of high gene expression, as well as gene coexpression. Results We offer a versatile modeling framework to identify transcriptional states as well as structures of coactivation for different neuronal cell types across multiple datasets. We employed a gamma-normal mixture model to identify active gene expression across cells, and used these to characterize markers for olfactory sensory neuron cell maturity, and to build cell-specific coactivation networks. We found that combined analysis of multiple datasets results in more known maturity markers being identified, as well as pointing towards some novel genes that may be involved in neuronal maturation. We also observed that the cell-specific coactivation networks of mature neurons tended to have a higher centralization network measure than immature neurons. Conclusion Integration of multiple datasets promises to bring about more statistical power to identify genes and patterns of interest. We found that transforming the data into active and inactive gene states allowed for more direct comparison of datasets, leading to identification of maturity marker genes and cell-specific network observations, taking into account the unique characteristics of single cell transcriptomics data. Electronic supplementary material The online version of this article (doi:10.1186/s12918-016-0370-4) contains supplementary material, which is available to authorized users.