Extended local similarity analysis (eLSA) of microbial community and other time series data with replicates
- Li C Xia^{1},
- Joshua A Steele^{2, 3},
- Jacob A Cram^{3},
- Zoe G Cardon^{4},
- Sheri L Simmons^{5},
- Joseph J Vallino^{4},
- Jed A Fuhrman^{3} and
- Fengzhu Sun^{1}Email author
https://doi.org/10.1186/1752-0509-5-S2-S15
© Xia et al; licensee BioMed Central Ltd. 2011
Published: 14 December 2011
Abstract
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.
Keywords
Background
In recent years, advances in microbial molecular technologies, such as next generation sequencing and molecular profiling, have enabled researchers to spatially and temporally characterize natural microbial communities without laboratory cultivation [1]. However, to reveal existing symbiotic relationships and microbe-environment interactions, it is necessary to mine and analyze temporal and spatial co-occurrence association patterns of organisms within these new datasets [2, 3]. Time series data, in particular, are receiving increased attention, since not only ordinary associations, but also other local and potentially time-delayed associations can be inferred from these datasets. Here local association indicates that the association only occurs in a subinterval of the time of interest, and time-delayed association indicates that there is a time lag for the response of one organism to the change in another organism. The rapid accrual of time series data is not limited to the microbial ecology field. Progress in high-throughput low-cost experimental technologies has also brought such changes to gene transcription and translation studies. Thus, while the subjects may vary, the association network we build from local and potentially time-delayed association patterns will likely pave the way to a better understanding of these systems.
To analyze microbial community and other data under various conditions, researchers typically use techniques such as Pearson’s Correlation Coefficient (PCC), principal component analysis (PCA), multi-dimensional scaling (MDS), discriminant function analysis (DFA) and canonical correlation analysis (CCA) [4–8]. Although these analytic methods yield interesting patterns, they generally analyze the data throughout the whole time interval of interest without considering potential local and time-delayed associations. We are specifically interested in discovering local and potentially time-delayed associations. Such associations have been shown to play important roles in understanding gene expression dynamics and the association of organisms in microbial communities [9–12].
To understand local and time-delayed associations, we originally designed a Local Similarity Analysis (LSA) for time series data measured typically at successive and equal time intervals without replicates [11]. Studies adopting the original LSA technique have shown interesting and novel discoveries for microbial community datasets. To name a few, Paver et al. [10] successfully applied LSA to study glycolate-utilizing bacterial and phytoplankton associations, while Shade et al. [13] used LSA to discover bacterial association dynamics during lake mixing.
Since biological experiments are often associated with many potential sources of noise, repeated measurements (replicates) are usually carried out in order to better assess inherent uncertainties of the quantities of interest [14]. Furthermore, data emerging from such experiments are typically analyzed by mean effect or by the development of profiles where variability is not properly accounted for [15]. Temporal and spatial data with replicates are being generated in Dr. Cardon’s laboratory and others. The lack of support for replicated data in the original LSA program has prevented its application to these new datasets. With replicates, it is possible to evaluate the variation of and to give a bootstrap confidence interval for the local similarity (LS) score as defined in Ruan et al.[11]. Furthermore, the original LSA is restricted by the low computing efficiency of the R language, as well as poor handling of missing values. In order to improve upon these issues and make the technique more accessible to the scientific community, we developed an extended LSA technique, named eLSA, and implemented it as a C++ extension to Python.
Briefly, given time series data of two factors and a user-constrained delay limit, eLSA finds the configuration of the data that yields the highest local similarity (LS) score, which is a type of similarity metric. For example, within a delay limit of two units, the first time spot of one series might be aligned to the third time spot of the other series, thus maximizing their LS. For a dataset of many factors, eLSA is applied to each pairwise combination of factors in the dataset. Candidate associations are then evaluated statistically by a permutation test, which calculates the p-value which is the proportion of scores exceeding the original LS score after shuffling the first series and re-evaluating the LS score many times, and by the false discovery rate (FDR q-value), which is used to correct multiple comparisons. Researchers can use eLSA to detect undirected associations, i.e., association patterns without time delays, and directed associations, where the change of one factor may temporally lead or follow another factor.
The organization of the paper is as follows. In the “Methods” section, we describe the LSA algorithm for calculating LS score with replicates, data normalization, estimation of confidence interval for the LS score, and testing the statistical significance of a LS score. In the “Results” section, we first show the efficacy of eLSA by simulations, then describe briefly the pipeline of eLSA, and finally apply the pipeline to analyze a microbiological dataset and a gene expression dataset. The paper concludes with some discussion and conclusions.
Methods
Pearson’s correlation coefficient-based analysis
where , , and are the means of X and Y, respectively. The statistical significance of r is tested by the fact that follows a t-distribution (degree of freedom: v = n – 2, mean: 0 and variance v / (v – 2)) when m = 1. For a pair of non-replicated series where m = 1, PCC is a straightforward and powerful method to test and identify linear relationship between two bivariate normally distributed random variables. It is widely adopted in the literature but with limitations. Specifically, when the real relationships are more complex, for example, the association between the two factors only occurs in a subinterval of the region of interest or the change of one factor has a time-delay in response to the change of another factor. Several methods, including the original LSA method, have been proposed to overcome such difficulties [11, 16].
Local similarity analysis with replicates
- (1)
For i, j in {1,2,…,n}^{2}:
- (2)
For i, j in {1,2,…,n}^{2} with |i – j|≤ D:
P_{i+1,j+1} = max{0,P_{i,j} + S_{ XY }[F(X_{ i }),F(Y_{ j })]} and
- (3)
P_{ max }(X,Y) = max_{1≤i,j≤n}P_{i,j} and
S_{ sgn }(X,Y) = sgn[P_{ max }(X,Y) – N_{ max }(X,Y)].
The S_{ max }(X,Y) obtained is the maximum local similarity score possible for all configurations of m-replicated time series X and Y within time-delay D. In this extended algorithm, the scalars x_{ i } ’s and y_{ i } ’ s from the non-replicated series in Ruan et al.[11] are replaced by vector functions F(X_{ i })’s and F(Y_{ j })’s to handle data with replicates. Alternatively, we can also consider F(X_{ i })’s and F(Y_{ j })’s as the same input data for the original algorithm in Ruan et al.[11], except that they are F-transformed data. In addition, this extended LSA framework easily accommodates the original version of LSA without replicates using m = 1 as a special case.
Different ways of summarizing the replicate data
Notice that the only additional component we introduced in the eLSA algorithm is the function F. Many reports have suggested different possible forms for F, and several computational methods have been proposed for summarizing the additional information available from replicates, including the simple average method (abbreviated as ‘simple’) and the Standard Deviation (SD)-weighted average method (abbreviated as ‘SD’), and the multivariate correlation coefficient method [17–19]. However, the result of the multivariate correlation coefficient method from Zhu et al.[17] can be shown to be the same as the ‘simple’ method. Therefore, in eLSA, we used the first two methods. We also propose the use of median in place of average and Median Absolute Deviation (MAD) in place of SD when robust statistics are needed to handle outliers [20]. The corresponding methods are named simple median method (abbreviated as ‘Med’) and MAD-weighted median method (abbreviated as ‘MAD’), respectively.
The ‘simple’ method is, in spirit, to take the mean profiles to represent the replicated series. In practice, we take F to be the simple average of repeated measurements: . The ‘SD’ method, on the other hand, takes the standard deviation of the replicates into account. Here we take F to be the replicate average divided by its standard deviation (SD): . Importantly, this method utilizes the variability information available, and, as such, it is claimed to be better than the ‘simple’ method in estimating the true correlation [18]. However, in order for the ‘SD’ method to be effective, a relatively large number of replicates, m, are needed, e.g., m ≥ 5. For a small number of replicates, the ‘SD’ method may not work well since the standard deviation may not be reliably estimated. Further, if we replace average with median and SD with MAD, we obtain the ‘Med’ method: F(X_{ i }) = Median(X_{ i }) and the ‘MAD’ method: , where MAD(X_{ i }) = Median(|X_{ i } – Median(X_{ i })|). The two transformations have similar properties as their corresponding average and SD versions, but they are more robust.
Bootstrap confidence interval for the LS score
With replicate data, researchers can study the variation of quantities of interest and to give their confidence intervals. Due to the complexity of calculating the LS score, the probability distribution of the LS score is hard to study theoretically. Thus, we resort to bootstrap to give a bootstrap confidence interval (CI) for the LS score. Bootstrap is a re-sampling method for studying the variation of an estimated quantity based on available sample data [21]. In this study, we use bootstrap to estimate a confidence interval for the LS score. For a given type I error α, the 1 – α confidence interval is the estimated range that covers the true value with probability 1 – α. Thus, for a given number, B, of bootstraps, we construct the bootstrap sample set , where each and are samples with replacement from X_{ i } and Y_{ j }, respectively. The rest of the calculation is the same as that used for the original data, and we obtain . Without the loss of generality, we suppose that these values are sorted in ascending order: . Then, a 1 – α bootstrap CI of S_{ max } can be estimated by , as suggested by Efron et al.[21].
Data normalization
where Φ is the cumulative distribution function of the standard normal distribution. We will take Z = Z_{[1:n]} obtained through the above procedure as the normalization of X. Therefore, the normalization steps are taken after the F-transformation.
Permutation test to evaluate the statistical significance of LSA association
where I(·) is the indicator function. With large enough number of permutations, we can evaluate the p-value to any desired accuracy.
False discovery rate (FDR) estimation
In most biological studies, a large number of factors need to be considered. If there are T factors, there will be eLSA pairwise calculations, representing its quadratic growth in T. In order to avoid many falsely declared associated pairs of factors, we need to correct for multiple testing. Many methods have been developed to correct for multiple testing and here we use the method by Storey et al.[23] to address this issue. In particular, we report the q-value, Q, for each pair of factors. The q-value for a pair of factors is the proportion of false positives incurred when that particular pair of factors is declared significant.
Computation complexity and implementation
For a single pair of time series, the time complexity for calculating the LS score using the dynamic programming algorithm is O(n), where n is the number of time points. The estimation of the bootstrap confidence interval for the LS score using B bootstraps will need O(Bn) calculations. The estimation of statistical significance for each pair of factors using L permutations will need O(Ln) calculations. Thus, the number of calculations for a full analysis of each pair of factors will be O(BLn). If there are a total of T factors, there are a total of pairs of factors that need to be compared. Thus, the number of calculations for a full analysis of T factors will be in the order of O(T^{2}BLn), which can be computationally intensive.
In summary, the internal support for replicates and the use of CI estimates are the two major methodological enhancements to LSA. The eLSA software, however, also incorporates other new features, such as faster permutation and false discovery rate evaluations and more options to handle missing values. Other implementation details are available from the software documentation.
Results
Simulations and benchmarks
We generated simulated data to show the efficacy of eLSA in capturing time-dependent association patterns, such as time-delayed associations and associations within a subinterval. We also studied the difference between the eLSA inference using the simple average (referred to as ‘simple’) method, the SD-weighted average method (referred to as ‘SD’), the median (referred to as “Med”) method, and the MAD (referred to as ‘MAD’) method.
Time-delayed association
We see that the two series closely follows each other if we shift the Y series three units toward right. In this particular example, the PCC is -0.258 (P=0.272) while the LS score using ‘simple’ averaging method is 0.507 with a p-value of 0.006. We did 1000 bootstraps and the 95% bootstrap confidence interval for this particular example is (0.448, 0.549). Therefore, this time-delayed association is only found significant by the eLSA analysis.
Association within a subinterval
In this case, we assume X and Y are positively associated within a subinterval and not associated in other regions. In our simulation, we generate 20 time spots of the two series by sampling (X_{ j },Y_{ j }) from a bivariate-normal distribution with mean μ = 0 and covariance matrix where ρ = 0.8 for 6 ≤ j ≤ 15, and ρ = 0 for j ≤ 5 or 16 ≤ j ≤ 20. The generated (X_{ j },Y_{ j })’s are further perturbed m times by a measurement disturbance ε_{ ij } : N(0, 0.01) to obtain the m-replicated series. One generated series are shown in Figure 1b for a typical simulation with m = 5.
We can see the two series mostly closely follow each other within the intended subinterval 6 ≤ j ≤ 15. In this particular example, the PCC is 0.258 (P=0.272) while the LS score using ‘simple’ averaging method is 0.428 with a p-value of 0.028. We did 1000 bootstraps and the 95% bootstrap confidence interval is (0.404, 0.446). This pattern is again uniquely captured by the eLSA analysis. In real applications, there are many other possibilities that two factors are associated without a significant Pearson or Spearman’s correlation coefficient. The eLSA can capture these associations as long as their LS score can be maximized through dynamically enumerating their configurations.
Different summarizing function
Mean and standard error of the estimated LS score
m=1 | m=5 | m=10 | m=15 | m=20 | ||||||
---|---|---|---|---|---|---|---|---|---|---|
F-function | mean | se. | mean | se. | mean | se. | mean | se. | mean | se. |
‘simple’ | .495 | .078 | .495 | .085 | .491 | .088 | .493 | .076 | .496 | .091 |
‘SD’ | na. | na. | .332 | .127 | .391 | .124 | .412 | .119 | .435 | .109 |
‘Med’ | .495 | .078 | .490 | .090 | .490 | .090 | .490 | .083 | .498 | .083 |
‘MAD’ | na. | na. | .494 | .115 | .302 | .128 | .325 | .129 | .371 | .119 |
Running time comparison
We benchmarked the running time performance of the new eLSA implementation and the old R script. For a dataset of 72 time series each with 35 time points, we tried eLSA analysis with 100 bootstraps, 1000 permutations and a delay limit of 3. It took the old script 20462 seconds to finish the computation while the new C++ program used 2054 seconds, which is about 9 times faster. Meanwhile, the new implementation also reduces the memory consumption and increases input/output efficiency. The benchmark is carried out on a “Dell, PE1950, Xeon E5420, 2.5GHz, 12010MB RAM” computing node.
The eLSA analysis pipeline
F-transformation and data normalization
The eLSA tool accepts a matrix file where each row is a time series for one factor. It fills up missing data by a user-specified method. Zero to third order spline-based methods and the nearest neighbour method as implemented in the Scipy (http://www.scipy.org) interpolation module are available. It then transforms the data by the user-specified F function and normalizes the F-transformed data by the normal score transformation following Li et al.[22] (see Methods).
Local similarity scoring
Local similarity analysis calculates the highest similarity score between any pair of factors. Users can specify parameters, including, for example, the maximum shifts allowed. Local Similarity score is calculated using the eLSA dynamic programming algorithm (see Methods).
Permutation test
The statistical significance, the p-value, of LS score is evaluated using a permutation test. Briefly, eLSA randomly shuffles the components of the original time series and recalculates the LS score for the pairs. The p-value is approximated by the fraction of permutation scores that are larger (in absolute value) than the original score. Confidence interval for a given LS score is also found by bootstrapping from the replicated data. Finally, users can obtain significant eLSA association results by the combined use of p-value and FDR q-value thresholds as their filtering criteria.
Association network construction
Using only the significant associations, users can construct a partially directed association network. Generally, for two factors X and Y, if the time interval [s_{1},t_{1}] in X and [s_{2},t_{2}] in Y have the highest LS and s_{1} <s_{2}, we can infer that X leads Y; in other words, X possibly activates Y. In network visualization software (e.g., Cytoscape [24]), one can use arrows to directionally indicate these lead patterns (i.e., X to Y, if X leads Y; otherwise undirected, if no direction is inferred). One can also use lines to indicate association types (solid, if X is positively associated with Y; otherwise dashed). Following these rules, one can build a partially directed association network based on eLSA results.
Microbial community data analysis
As an immediate application, we applied the eLSA pipeline to a set of real microbial community time series data. This San Pedro Ocean Time Series (SPOTs) dataset, originally reported in Steele et al.[2] and Countway et al.[25], was collected following a biological feature (i.e. the chlorophyll maximum depth) off the coast of Southern California. The bacterial community was analyzed using the ARISA [4] technique and the protistan community was analyzed using the T-RFLP [26] technique. The dataset is composed of monthly sampled data from September 2000 to March 2004, including 40 time points without replicates. We analyzed the dataset with a delay limit of 3 months and 1000 permutations to evaluate the statistical significance of the LSA score. In this dataset, the factor names, including the operational taxonomic units and environmental factors, are previously defined by Steele et al. [2].
Significant associations found in real datasets
Found by eLSA | Found by PCC | Found by both | Found by eLSA | Found by PCC | Found by both | ||
---|---|---|---|---|---|---|---|
Dataset | # of factors | P ≤ 0.01 | P ≤ 0.01 | P ≤ 0.01 | P ≤ 0.05 | P ≤ 0.05 | P ≤ 0.05 |
Q ≤ 0.01 | Q ≤ 0.01 | Q ≤ 0.01 | Q ≤ 0.05 | Q ≤ 0.05 | Q ≤ 0.05 | ||
Microbial | 515 | 1643 | 3237 | 293 | 2804 | 4242 | 658 |
C. elegans | 446 | 42532 | 56605 | 39114 | 57991 | 71799 | 54201 |
Top LS scores from the microbial community data
X | Y | LS | Xs | Ys | Len | D | P | PCC | Ppcc | Q | Qpcc |
---|---|---|---|---|---|---|---|---|---|---|---|
Euk239 | Euk269 | 0.82 | 1 | 1 | 40 | 0 | 0 | 0.09 | 0.59 | 0.02 | 1.00 |
Bac609 | Bac675 | 0.77 | 1 | 2 | 39 | -1 | 0 | 0.14 | 0.41 | 0.00 | 1.00 |
Euk381 | Euk462 | 0.77 | 1 | 1 | 40 | 0 | 0 | 0.44 | 0.00 | 0.02 | 0.11 |
Euk583 | Bac989 | 0.68 | 2 | 1 | 39 | 1 | 0 | 0.30 | 0.06 | 0.02 | 0.73 |
Euk229 | Euk339 | 0.57 | 1 | 2 | 39 | -1 | 0 | 0.05 | 0.77 | 0.02 | 1.00 |
Euk97 | boxy | -0.62 | 15 | 15 | 21 | 0 | 0 | -0.42 | 0.01 | 0.00 | 0.17 |
Euk98 | boxy | -0.62 | 15 | 15 | 21 | 0 | 0 | -0.42 | 0.01 | 0.00 | 0.17 |
Euk109 | boxy | -0.62 | 15 | 15 | 21 | 0 | 0 | -0.42 | 0.01 | 0.00 | 0.17 |
Euk112 | boxy | -0.62 | 15 | 15 | 21 | 0 | 0 | -0.42 | 0.01 | 0.00 | 0.17 |
Euk116 | boxy | -0.62 | 15 | 15 | 21 | 0 | 0 | -0.42 | 0.01 | 0.00 | 0.17 |
Gene expression data analysis
Although LSA had its roots grounded in microbial community analysis, the technique can be readily applied to other biological time series data, such as replicated gene expression time series data from microarray and RNA-Seq experiments [27–29]. Here we show an example of applying eLSA to the dauer exit gene expression profile time series data of 446 genes from a C elegans study. The result of the original study suggests that the 446 genes under investigation have similar kinetics in both the dauer exit and the L1 starvation time course [30]. Here we use the dauer exit time series data consisting of 12 hourly time spots, each with four replicates. We analyzed the dataset with a delay limit of 3 hours and with 1000 permutations and 100 bootstraps.
The results are summarized in Table 2. Comparing the C. elegans results to those of the microbial community, we see that gene-gene associations in this network are much denser, since a smaller number of genes end up with a much larger, rather than smaller, number of eLSA significant associations (e.g. 2804 versus 57991 for Q ≤ 0.05 and P ≤ 0.05, see Table 2). Also different is that about 93% of these associations are found by PCC analysis as well. The high congruence between PCC and eLSA analysis may be due to the fact that about 90% of the eLSA findings are without delays, which thus are also amenable to PCC analysis.
Because these genes do not change expression level in both dauer exit and L1 starvation conditions, they are considered as common feeding response genes [30]. However, it is not clear whether they are correlated with each other in expression profiles under the dauer exit condition. To study this, we combined all eLSA and PCC significant associations with Q ≤ 0.05 and P ≤ 0.05, and found the average degree of the resulting association network is around 169, while that of previous microbial community data is around 12. Such high average degree for C. elegans genes shows the high similarity of their expression profiles, which also reflects their intimate functional coordination along the process. Therefore, our result suggests those feeding response genes are likely to be co-expressed under the dauer exit condition.
Top LS scores from the C. elegans gene-expression data
X | Y | LS | lowCI | upCI | Xs | Ys | Len | D | P | PCC | Ppcc | Q | Qpcc |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
48087 | 27993 | 0.53 | 0.41 | 0.61 | 1 | 2 | 11 | -1 | 0.00 | 0.56 | 0.06 | 0.00 | 0.01 |
32607 | 51986 | 0.52 | 0.41 | 0.61 | 2 | 1 | 10 | 1 | 0.01 | 0.51 | 0.09 | 0.00 | 0.01 |
29504 | 48087 | 0.52 | 0.40 | 0.61 | 2 | 1 | 11 | 1 | 0.00 | 0.41 | 0.18 | 0.00 | 0.03 |
23193 | 27993 | 0.51 | 0.41 | 0.59 | 1 | 2 | 11 | -1 | 0.00 | 0.48 | 0.11 | 0.00 | 0.02 |
29494 | 30208 | 0.51 | 0.39 | 0.61 | 2 | 1 | 11 | 1 | 0.00 | 0.58 | 0.05 | 0.00 | 0.01 |
27993 | 53694 | -0.55 | -0.62 | -0.44 | 2 | 1 | 11 | 1 | 0.00 | -0.53 | 0.08 | 0.00 | 0.01 |
436287 | 53694 | -0.54 | -0.62 | -0.44 | 2 | 1 | 11 | 1 | 0.01 | -0.55 | 0.06 | 0.00 | 0.01 |
48941 | 53694 | -0.52 | -0.61 | -0.42 | 2 | 1 | 11 | 1 | 0.00 | -0.38 | 0.22 | 0.00 | 0.03 |
29494 | 22857 | -0.52 | -0.61 | -0.41 | 2 | 1 | 11 | 1 | 0.00 | -0.49 | 0.10 | 0.00 | 0.02 |
29494 | 436727 | -0.52 | -0.61 | -0.40 | 2 | 1 | 11 | 1 | 0.01 | -0.55 | 0.06 | 0.00 | 0.01 |
Discussion and conclusions
The eLSA technique extends LSA to time series data with replicates. This will help investigators better utilize the available information from their sample replicates and assist them in more effective and reliable hypothesis generation of time-dependent associations. In addition, a bootstrap framework is developed to estimate the confidence interval for the LS score. We also provided flexible missing value options and integrated efficient multiple testing control methods for the new eLSA technique. Using the microbial community and gene expression datasets, we demonstrated that eLSA uniquely captures additional time-dependent associations, including local and time-delayed association patterns, when compared to ordinary correlation methods, such as PCC. In this paper, we described the applications of our method with the time series data. Actually, the eLSA can be applied to any type of data with some gradients, including the response to different levels of treatments, temperature, humidity, or spatial distributions.
Currently, we use permutation test to assess the statistical significance of LS scores and bootstrap re-sampling to estimate the confidence interval of LS score. Both the permutation test and bootstrap methods are time consuming if high precise determination of statistical significance or confidence interval is desired. Theoretical developments on the distribution of the LS score are needed to eliminate or mitigate the computational burden required for these processes, and would be interesting topics for future studies. There is also a minimum sample number requirement for eLSA analysis. We suggest the sample number to be greater than 5+D, where D is the desired delay limit, since shifting and trimming by eLSA will further reduce the effective sample number and result in lower statistical power.
Declarations
Acknowledgements
The authors would like to thank Cheryl Chow, Rohan Sachdeva, Barbara Campbell, Anders Andersson and Stefan Bertilsson for testing the eLSA software packages and web services and providing valuable suggestions. We thank Jun Zhao of PIBBS at University of Southern California for helpful discussion of C. elegans dataset analysis. We thank an anonymous reviewer for suggesting the “Med” and the “MAD” approaches. We also thank the Molecular Computational Biology Program at University of Southern California for providing computing resources. This research is partially supported by the National Science Foundation (NSF) DMS-1043075 and OCE 1136818.
This article has been published as part of BMC Systems Biology Volume 5 Supplement 2, 2011: 22nd International Conference on Genome Informatics: Systems Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/1752-0509/5?issue=S2.
Authors’ Affiliations
References
- Fuhrman JA: Microbial community structure and its functional implications. Nature. 2009, 459: 193-199. 10.1038/nature08058.View ArticlePubMedGoogle Scholar
- Steele JA, Countway PD, Xia L, Vigil PD, Beman JM, Kim DY, Chow CE, Sachdeva R, Jones AC, Schwalbach MS: Marine bacterial, archaeal and protistan association networks reveal ecological linkages. ISME J. 2011, 5: 1414-1425. 10.1038/ismej.2011.24.PubMed CentralView ArticlePubMedGoogle Scholar
- Chaffron S, Rehrauer H, Pernthaler J, von Mering C: A global network of coexisting microbes from environmental and whole-genome sequence data. Genome Res. 2010, 20: 947-959. 10.1101/gr.104521.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Fisher MM, Triplett EW: Automated approach for ribosomal intergenic spacer analysis of microbial diversity and its application to freshwater bacterial communities. Appl Environ Microbiol. 1999, 65: 4630-4636.PubMed CentralPubMedGoogle Scholar
- Stepanauskas R, Moran MA, Bergamaschi BA, Hollibaugh JT: Covariance of bacterioplankton composition and environmental variables in a temperate delta system. Aquat Microb Ecol. 2003, 31: 85-98.View ArticleGoogle Scholar
- Van Mooy BAS, Devol AH, Keil RG: Relationship between bacterial community structure, light, and carbon cycling in the eastern subarctic North Pacific. Limnology and Oceanography. 2004, 49: 1056-1062. 10.4319/lo.2004.49.4.1056.View ArticleGoogle Scholar
- Yannarell AC, Triplett EW: Geographic and environmental sources of variation in lake bacterial community composition. Appl Environ Microbiol. 2005, 71: 227-239. 10.1128/AEM.71.1.227-239.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Yannarell AC, Triplett EW: Within- and between-lake variability in the composition of bacterioplankton communities: investigations using multiple spatial scales. Appl Environ Microbiol. 2004, 70: 214-223. 10.1128/AEM.70.1.214-223.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Li X, Rao S, Jiang W, Li C, Xiao Y, Guo Z, Zhang Q, Wang L, Du L, Li J: Discovery of Time-Delayed Gene Regulatory Networks based on temporal gene expression profiling. BMC Bioinformatics. 2006, 7: 26-10.1186/1471-2105-7-26.PubMed CentralView ArticlePubMedGoogle Scholar
- Paver SF, Kent AD: Temporal patterns in glycolate-utilizing bacterial community composition correlate with phytoplankton population dynamics in humic lakes. Microb Ecol. 2010, 60: 406-418. 10.1007/s00248-010-9722-6.View ArticlePubMedGoogle Scholar
- Ruan Q, Dutta D, Schwalbach MS, Steele JA, Fuhrman JA, Sun F: Local similarity analysis reveals unique associations among marine bacterioplankton species and environmental factors. Bioinformatics. 2006, 22: 2532-2538. 10.1093/bioinformatics/btl417.View ArticlePubMedGoogle Scholar
- Wang G, Yin L, Zhao Y, Mao K: Efficiently mining time-delayed gene expression patterns. IEEE Trans Syst Man Cybern B Cybern. 2010, 40: 400-411.View ArticlePubMedGoogle Scholar
- Shade A, Chiu CY, McMahon KD: Differential bacterial dynamics promote emergent community robustness to lake mixing: an epilimnion to hypolimnion transplant experiment. Environ Microbiol. 2010, 12: 455-466. 10.1111/j.1462-2920.2009.02087.x.View ArticlePubMedGoogle Scholar
- Lee ML, Kuo FC, Whitmore GA, Sklar J: Importance of replication in microarray gene expression studies: statistical methods and evidence from repetitive cDNA hybridizations. Proc Natl Acad Sci U S A. 2000, 97: 9834-9839. 10.1073/pnas.97.18.9834.PubMed CentralView ArticlePubMedGoogle Scholar
- Nguyen TT, Almon RR, DuBois DC, Jusko WJ, Androulakis IP: Importance of replication in analyzing time-series gene expression data: corticosteroid dynamics and circadian patterns in rat liver. BMC Bioinformatics. 2010, 11: 279-10.1186/1471-2105-11-279.PubMed CentralView ArticlePubMedGoogle Scholar
- Balasubramaniyan R, Hullermeier E, Weskamp N, Kamper J: Clustering of gene expression data using a local shape-based similarity measure. Bioinformatics. 2005, 21: 1069-1077. 10.1093/bioinformatics/bti095.View ArticlePubMedGoogle Scholar
- Zhu D, Li Y, Li H: Multivariate correlation estimator for inferring functional relationships from replicated genome-wide data. Bioinformatics. 2007, 23: 2298-2305. 10.1093/bioinformatics/btm328.View ArticlePubMedGoogle Scholar
- Yao J, Chang C, Salmi ML, Hung YS, Loraine A, Roux SJ: Genome-scale cluster analysis of replicated microarrays using shrinkage correlation coefficient. BMC Bioinformatics. 2008, 9: 288-10.1186/1471-2105-9-288.PubMed CentralView ArticlePubMedGoogle Scholar
- Littell RC, Pendergast J, Natarajan R: Modelling covariance structure in the analysis of repeated measures data. Stat Med. 2000, 19: 1793-1819. 10.1002/1097-0258(20000715)19:13<1793::AID-SIM482>3.0.CO;2-Q.View ArticlePubMedGoogle Scholar
- Venables WN, Ripley BD: Modern Applied Statistics with S. 1997, Springer, 4View ArticleGoogle Scholar
- Efron B, Tibshirani R: An Introduction to the Bootstrap. 1998, Boca Raton; London: Chapman & Hall/CRCGoogle Scholar
- Li KC: Genome-wide coexpression dynamics: theory and application. Proc Natl Acad Sci U S A. 2002, 99: 16875-16880. 10.1073/pnas.252466999.PubMed CentralView ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.PubMed CentralView ArticlePubMedGoogle Scholar
- Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, Christmas R, Avila-Campilo I, Creech M, Gross B: Integration of biological networks and gene expression data using Cytoscape. Nat Protoc. 2007, 2: 2366-2382. 10.1038/nprot.2007.324.PubMed CentralView ArticlePubMedGoogle Scholar
- Countway PD, Vigil PD, Schnetzer A, Moorthi SD, Caron DA: Seasonal analysis of protistan community structure and diversity at the USC Microbial Observatory (San Pedro Channel, North Pacific Ocean). Limnology and Oceanography. 2010, 55: 2381-2396. 10.4319/lo.2010.55.6.2381.View ArticleGoogle Scholar
- Vigil P, Countway PD, Rose J, Lonsdale DJ, Gobler CJ, Caron DA: Rapid shifts in dominant taxa among microbial eukaryotes in estuarine ecosystems. Aquat Microb Ecol. 2008, 54: 83-100.View ArticleGoogle Scholar
- Bar-Joseph Z: Analyzing time series gene expression data. Bioinformatics. 2004, 20: 2493-2503. 10.1093/bioinformatics/bth283.View ArticlePubMedGoogle Scholar
- Tai YC, Speed TP: On gene ranking using replicated microarray time course data. Biometrics. 2009, 65: 40-51. 10.1111/j.1541-0420.2008.01057.x.View ArticlePubMedGoogle Scholar
- Tai YC, Speed TP: A multivariate empirical Bayes statistic for replicated microarray time course data. Ann Stat. 2006, 34: 2387-2412. 10.1214/009053606000000759.View ArticleGoogle Scholar
- Wang J, Kim SK: Global analysis of dauer gene expression in Caenorhabditis elegans. Development. 2003, 130: 1621-1634. 10.1242/dev.00363.View ArticlePubMedGoogle Scholar
- Kapp LD, Lorsch JR: The molecular mechanics of eukaryotic translation. Annu Rev Biochem. 2004, 73: 657-704. 10.1146/annurev.biochem.73.030403.080419.View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.