Multivariate correlation estimator for inferring functional relationships from replicated genome-wide data

Identifying local associations in biological time series: algorithms, statistical significance, and applications

Local associations refer to spatial-temporal correlations that emerge from the biological realm, such as time-dependent gene co-expression or seasonal interactions between microbes. One can reveal the intricate dynamics and inherent interactions of biological systems by examining the biological time series data for these associations. To accomplish this goal, local similarity analysis algorithms and statistical methods that facilitate the local alignment of time series and assess the significance of the resulting alignments have been developed. Although these algorithms were initially devised for gene expression analysis from microarrays, they have been adapted and accelerated for multi-omics next generation sequencing datasets, achieving high scientific impact. In this review, we present an overview of the historical developments and recent advances for local similarity analysis algorithms, their statistical properties, and real applications in analyzing biological time series data. The benchmark data and analysis scripts used in this review are freely available at http://github.com/labxscut/lsareview.

Cluster analysis of replicated alternative polyadenylation data using canonical correlation analysis

BACKGROUND

Alternative polyadenylation (APA) has emerged as a pervasive mechanism that contributes to the transcriptome complexity and dynamics of gene regulation. The current tsunami of whole genome poly(A) site data from various conditions generated by 3' end sequencing provides a valuable data source for the study of APA-related gene expression. Cluster analysis is a powerful technique for investigating the association structure among genes, however, conventional gene clustering methods are not suitable for APA-related data as they fail to consider the information of poly(A) sites (e.g., location, abundance, number, etc.) within each gene or measure the association among poly(A) sites between two genes.

RESULTS

Here we proposed a computational framework, named PASCCA, for clustering genes from replicated or unreplicated poly(A) site data using canonical correlation analysis (CCA). PASCCA incorporates multiple layers of gene expression data from both the poly(A) site level and gene level and takes into account the number of replicates and the variability within each experimental group. Moreover, PASCCA characterizes poly(A) sites in various ways including the abundance and relative usage, which can exploit the advantages of 3' end deep sequencing in quantifying APA sites. Using both real and synthetic poly(A) site data sets, the cluster analysis demonstrates that PASCCA outperforms other widely-used distance measures under five performance metrics including connectivity, the Dunn index, average distance, average distance between means, and the biological homogeneity index. We also used PASCCA to infer APA-specific gene modules from recently published poly(A) site data of rice and discovered some distinct functional gene modules. We have made PASCCA an easy-to-use R package for APA-related gene expression analyses, including the characterization of poly(A) sites, quantification of association between genes, and clustering of genes.

CONCLUSIONS

By providing a better treatment of the noise inherent in repeated measurements and taking into account multiple layers of poly(A) site data, PASCCA could be a general tool for clustering and analyzing APA-specific gene expression data. PASCCA could be used to elucidate the dynamic interplay of genes and their APA sites among various biological conditions from emerging 3' end sequencing data to address the complex biological phenomenon.

Time-resolved proteome profiling of normal lung development

Biochemical networks mediating normal lung morphogenesis and function have important implications for ameliorating morbidity and mortality in premature infants. Although several transcript-level studies have examined normal lung development, corresponding protein-level analyses are lacking. Here we performed proteomics analysis of murine lungs from embryonic to early adult ages to identify the molecular networks mediating normal lung development. We identified 8,932 proteins, providing a deep and comprehensive view of the lung proteome. Analysis of the proteomics data revealed discrete modules and the underlying regulatory and signaling network modulating their expression during development. Our data support the cell proliferation that characterizes early lung development and highlight responses of the lung to exposure to a nonsterile oxygen-rich ambient environment and the important role of lipid (surfactant) metabolism in lung development. Comparison of dynamic regulation of proteomic and recent transcriptomic analyses identified biological processes under posttranscriptional control. Our study provides a unique proteomic resource for understanding normal lung formation and function and can be freely accessed at Lungmap.net.

Uncovering robust patterns of microRNA co-expression across cancers using Bayesian Relevance Networks

Co-expression networks have long been used as a tool for investigating the molecular circuitry governing biological systems. However, most algorithms for constructing co-expression networks were developed in the microarray era, before high-throughput sequencing-with its unique statistical properties-became the norm for expression measurement. Here we develop Bayesian Relevance Networks, an algorithm that uses Bayesian reasoning about expression levels to account for the differing levels of uncertainty in expression measurements between highly- and lowly-expressed entities, and between samples with different sequencing depths. It combines data from groups of samples (e.g., replicates) to estimate group expression levels and confidence ranges. It then computes uncertainty-moderated estimates of cross-group correlations between entities, and uses permutation testing to assess their statistical significance. Using large scale miRNA data from The Cancer Genome Atlas, we show that our Bayesian update of the classical Relevance Networks algorithm provides improved reproducibility in co-expression estimates and lower false discovery rates in the resulting co-expression networks. Software is available at www.perkinslab.ca.

A generalized dSpliceType framework to detect differential splicing and differential expression events using RNA-Seq

Transcriptomes are routinely compared in term of a list of differentially expressed genes followed by functional enrichment analysis. Due to the technology limitations of microarray, the molecular mechanisms of differential expression is poorly understood. Using RNA-seq data, we propose a generalized dSpliceType framework to systematically investigate the synergistic and antagonistic effects of differential splicing and differential expression. We applied the method to two public RNA-seq data sets and compared the transcriptomes between treatment and control conditions. The generalized dSpliceType detects and prioritizes a list of genes that are differentially expressed and/or spliced. In particular, the multivariate dSpliceType is among the fist to utilize sequential dependency of normalized base-wise read coverage signals and capture biological variability among replicates using a multivariate statistical model. We compared dSpliceType with two other methods in terms of five most common types of differential splicing events between two conditions using RNA-Seq. dSpliceType is free, available from http://dsplicetype.sourceforge.net/.

CorSig: a general framework for estimating statistical significance of correlation and its application to gene co-expression analysis

With the rapid increase of omics data, correlation analysis has become an indispensable tool for inferring meaningful associations from a large number of observations. Pearson correlation coefficient (PCC) and its variants are widely used for such purposes. However, it remains challenging to test whether an observed association is reliable both statistically and biologically. We present here a new method, CorSig, for statistical inference of correlation significance. CorSig is based on a biology-informed null hypothesis, i.e., testing whether the true PCC (ρ) between two variables is statistically larger than a user-specified PCC cutoff (τ), as opposed to the simple null hypothesis of ρ = 0 in existing methods, i.e., testing whether an association can be declared without a threshold. CorSig incorporates Fisher's Z transformation of the observed PCC (r), which facilitates use of standard techniques for p-value computation and multiple testing corrections. We compared CorSig against two methods: one uses a minimum PCC cutoff while the other (Zhu's procedure) controls correlation strength and statistical significance in two discrete steps. CorSig consistently outperformed these methods in various simulation data scenarios by balancing between false positives and false negatives. When tested on real-world Populus microarray data, CorSig effectively identified co-expressed genes in the flavonoid pathway, and discriminated between closely related gene family members for their differential association with flavonoid and lignin pathways. The p-values obtained by CorSig can be used as a stand-alone parameter for stratification of co-expressed genes according to their correlation strength in lieu of an arbitrary cutoff. CorSig requires one single tunable parameter, and can be readily extended to other correlation measures. Thus, CorSig should be useful for a wide range of applications, particularly for network analysis of high-dimensional genomic data.

SOFTWARE AVAILABILITY

A web server for CorSig is provided at http://202.127.200.1:8080/probeWeb. R code for CorSig is freely available for non-commercial use at http://aspendb.uga.edu/downloads.

Graphical identification of cancer-associated gene subnetworks based on small proteomics data sets

Proteomics is a rapidly emerging frontier in post-genomics medicine and biology, but the quantitative analysis and validation of proteomic data are in need of further improvements. Before selecting potential candidate proteomic biomarkers, it is important to understand the broader context of how biological processes are regulated under different conditions or in different phenotypes. The enrichment of proteomic data consists of extracting as much biological meaning as possible from curated, pathway-based, functional protein interaction networks. Currently, most of the enrichment tools are intended for microarray data and require parametric data, whereas proteomic data are often nonparametric. In this study, we aimed to select a suite of interactive tools that can enrich proteomic results with a graphical overview. This facilitated diagnosis and interpretation prior to further analysis. From a list of proteins, a network was constructed using a map of the most severely disrupted biological process, and the disease entity was then identified on the basis of clinical data. Taken together, this graphical and interactive method ranks potential proteins via functional analysis in order to improve the choice of biomarkers for validation with the following advantages: 1) It adds neighbor proteins that are not selected by mass spectrometry analysis, but could in fact be key proteins; 2) pinpoints the biological process most often involved; and 3) predicts the most likely disease on the basis of clinical data.

Assessing numerical dependence in gene expression summaries with the jackknife expression difference

Statistical methods to test for differential expression traditionally assume that each gene's expression summaries are independent across arrays. When certain preprocessing methods are used to obtain those summaries, this assumption is not necessarily true. In general, the erroneous assumption of dependence results in a loss of statistical power. We introduce a diagnostic measure of numerical dependence for gene expression summaries from any preprocessing method and discuss the relative performance of several common preprocessing methods with respect to this measure. Some common preprocessing methods introduce non-trivial levels of numerical dependence. The issue of (between-array) dependence has received little if any attention in the literature, and researchers working with gene expression data should not take such properties for granted, or they risk unnecessarily losing statistical power.

Extended local similarity analysis (eLSA) of microbial community and other time series data with replicates

Background

The increasing availability of time series microbial community data from metagenomics and other molecular biological studies has enabled the analysis of large-scale microbial co-occurrence and association networks. Among the many analytical techniques available, the Local Similarity Analysis (LSA) method is unique in that it captures local and potentially time-delayed co-occurrence and association patterns in time series data that cannot otherwise be identified by ordinary correlation analysis. However LSA, as originally developed, does not consider time series data with replicates, which hinders the full exploitation of available information. With replicates, it is possible to understand the variability of local similarity (LS) score and to obtain its confidence interval.

Results

We extended our LSA technique to time series data with replicates and termed it extended LSA, or eLSA. Simulations showed the capability of eLSA to capture subinterval and time-delayed associations. We implemented the eLSA technique into an easy-to-use analytic software package. The software pipeline integrates data normalization, statistical correlation calculation, statistical significance evaluation, and association network construction steps. We applied the eLSA technique to microbial community and gene expression datasets, where unique time-dependent associations were identified.

Conclusions

The extended LSA analysis technique was demonstrated to reveal statistically significant local and potentially time-delayed association patterns in replicated time series data beyond that of ordinary correlation analysis. These statistically significant associations can provide insights to the real dynamics of biological systems. The newly designed eLSA software efficiently streamlines the analysis and is freely available from the eLSA homepage, which can be accessed at http://meta.usc.edu/softs/lsa.

Effects of a 2-y dietary weight-loss intervention on cholesterol metabolism in moderately obese men

BACKGROUND

Long-term dietary weight loss results in complex metabolic changes. However, its effect on cholesterol metabolism in obese subjects is still unclear.

OBJECTIVE

We assessed the effects of 2 y of weight loss achieved with various diet regimens on phytosterols (markers of intestinal cholesterol absorption), lanosterol (marker of de novo cholesterol synthesis), and changes in apolipoprotein concentrations.

DESIGN

We conducted the 2-y Dietary Intervention Randomized Controlled Trial (DIRECT-a study of low-fat, Mediterranean, and low-carbohydrate diets). We assessed circulating phytosterol and lanosterol concentrations and their ratios to cholesterol and apolipoproteins A-I and B-100 in 90 DIRECT participants at 0, 6, and 24 mo.

RESULTS

We observed a significant upregulation of the markers of cholesterol absorption (campesterol: +16.8%, P < 0.001) and a downregulation of the markers of cholesterol synthesis (lanosterol: -16.5%, P = 0.008) during the active weight-loss phase (first 6 mo, weight loss of 5%, 6%, and 10% in the 3 diet groups, respectively), followed by a rebound (campesterol: -6.2%, P = 0.045; lanosterol: +43.7%, P < 0.001) during the next 18 mo (weight gain of 1%, 1%, and 2% in the 3 diet groups, respectively). HDL cholesterol continuously increased during the study (17.0%, P < 0.001), whereas LDL cholesterol remained constant. At the end of the 24-mo follow-up period, campesterol (P < 0.001) and lanosterol (P = 0.016) amounts were significantly higher than baseline values. The concentration of apolipoprotein B-100 correlated with cholesterol metabolism (ρ = 0.299 and P = 0.020 for lanosterol; ρ = -0.105 and NS for campesterol), and the homeostasis model assessment of insulin resistance correlated with lanosterol (ρ = 0.09, P = 0.001).

CONCLUSIONS

Long-term weight loss is related to a characteristic response suggestive of altered cholesterol and apolipoprotein metabolism. Various diets have a similar effect on these effects. DIRECT is registered at clinicaltrials.gov as NCT00160108.

A generalized multivariate approach to pattern discovery from replicated and incomplete genome-wide measurements

Estimation of pairwise correlation from incomplete and replicated molecular profiling data is an ubiquitous problem in pattern discovery analysis, such as clustering and networking. However, existing methods solve this problem by ad hoc data imputation, followed by aveGation coefficient type approaches, which might annihilate important patterns present in the molecular profiling data. Moreover, these approaches do not consider and exploit the underlying experimental design information that specifies the replication mechanisms. We develop an Expectation-Maximization (EM) type algorithm to estimate the correlation structure using incomplete and replicated molecular profiling data with a priori known replication mechanism. The approach is sufficiently generalized to be applicable to any known replication mechanism. In case of unknown replication mechanism, it is reduced to the parsimonious model introduced previously. The efficacy of our approach was first evaluated by comprehensively comparing various bivariate and multivariate imputation approaches using simulation studies. Results from real-world data analysis further confirmed the superior performance of the proposed approach to the commonly used approaches, where we assessed the robustness of the method using data sets with up to 30 percent missing values.

Cross-platform analysis of global microRNA expression technologies

Background

Although analysis of microRNAs (miRNAs) by DNA microarrays is gaining in popularity, these new technologies have not been adequately validated. We examined within and between platform reproducibility of four miRNA array technologies alongside TaqMan PCR arrays.

Results

Two distinct pools of reference materials were selected in order to maximize differences in miRNA content. Filtering for miRNA that yielded signal above background revealed 54 miRNA probes (matched by sequence) across all platforms. Using this probeset as well as all probes that were present on an individual platform, within-platform analyses revealed Spearman correlations of >0.9 for most platforms. Comparing between platforms, rank analysis of the log ratios of the two reference pools also revealed high correlation (range 0.663-0.949). Spearman rank correlation and concordance correlation coefficients for miRNA arrays against TaqMan qRT-PCR arrays were similar for all of the technologies. Platform performances were similar to those of previous cross-platform exercises on mRNA and miRNA microarray technologies.

Conclusions

These data indicate that miRNA microarray platforms generated highly reproducible data and can be recommended for the study of changes in miRNA expression.

Effects of scanning sensitivity and multiple scan algorithms on microarray data quality

Background

Maximizing the utility of DNA microarray data requires optimization of data acquisition through selection of an appropriate scanner setting. To increase the amount of useable data, several approaches have been proposed that incorporate multiple scans at different sensitivities to reduce the quantification error and to minimize effects of saturation. However, no direct comparison of their efficacy has been made. In the present study we compared individual scans at low, medium and high sensitivity with three methods for combining data from multiple scans (either 2-scan or 3-scan cases) using an actual dataset comprising 40 technical replicates of a reference RNA standard.

Results

Of the individual scans, the low scan exhibited the lowest background signal, the highest signal-to-noise ratio, and equivalent reproducibility to the medium and high scans. Most multiple scan approaches increased the range of probe intensities compared to the individual scans, but did not increase the dynamic range (the proportion of useable data). Approaches displayed striking differences in the background signal and signal-to-noise ratio. However, increased probe intensity range and improved signal-to-noise ratios did not necessarily correlate with improved reproducibility. Importantly, for one multiple scan method that combined 3 scans, reproducibility was significantly improved relative to individual scans and all other multiple scan approaches. The same method using 2 scans yielded significantly lower reproducibility, attributable to a lack-of-fit of the statistical model.

Conclusions

Our data indicate that implementation of a suitable multiple scan approach can improve reproducibility, but that model validation is critical to ensure accurate estimates of probe intensity.

A statistical framework for consolidating "sibling" probe sets for Affymetrix GeneChip data

Background

Affymetrix GeneChip typically contains multiple probe sets per gene, defined as sibling probe sets in this study. These probe sets may or may not behave similar across treatments. The most appropriate way of consolidating sibling probe sets suitable for analysis is an open problem. We propose the Analysis of Variance (ANOVA) framework to decide which sibling probe sets can be consolidated.

Results

The ANOVA model allows us to separate the sibling probe sets into two types: those behave similarly across treatments and those behave differently across treatments. We found that consolidation of sibling probe sets of the former type results in large increase in the number of differentially expressed genes under various statistical criteria. The approach to selecting sibling probe sets suitable for consolidating is implemented in R language and freely available from .

Conclusion

Our ANOVA analysis of sibling probe sets provides a statistical framework for selecting sibling probe sets for consolidation. Consolidating sibling probe sets by pooling data from each greatly improves the estimates of a gene expression level and results in identification of more biologically relevant genes. Sibling probe sets that do not qualify for consolidation may represent annotation errors or other artifacts, or may correspond to differentially processed transcripts of the same gene that require further analysis.