Modeling the temporal dynamics of the gut microbial community in adults and infants
There is increasing recognition that the human gut microbiome is a contributor to many aspects of human physiology and health including obesity , non-alcoholic fatty liver disease , inflammatory diseases , cancer , metabolic diseases , aging , and neurodegenerative disorders 1–14 ., This suggests that the human gut microbiome may play important roles in the diagnosis , treatment , and ultimately prevention of human disease ., These applications require an understanding of the temporal variability of the microbiota over the lifespan of an individual particularly since we now recognize that our microbiota is highly dynamic , and that the mechanisms underlying these changes are linked to ecological resilience and host health 15–17 ., Due to the lack of data and insufficient methodology , we currently have major gaps in our understanding of fundamental mechanisms related to the temporal behavior of the microbiome ., Critically , we currently do not have a clear characterization of how and why our gut microbiome varies in time , and whether these dynamics are consistent across humans ., It is also unclear whether we can define ‘stable’ or ‘healthy’ dynamics as opposed to ‘abnormal’ or ‘unhealthy’ dynamics , which could potentially reflect an underlying health condition or an environmental factor affecting the individual , such as antibiotics exposure or diet ., Moreover , there is no consensus as to whether the gut microbial community structure varies continuously or jumps between discrete community states , and whether or not these states are shared across individuals 18 , 19 ., Notably , recent work 20 suggests that the human gut microbiome composition is dominated by environmental factors rather than by host genetics , emphasizing the dynamic nature of this ecosystem ., The need for understanding the temporal dynamics of the microbiome and its interaction with host attributes have led to a rise in longitudinal studies that record the temporal variation of microbial communities in a wide range of environments , including the human gut microbiome ., These time series studies are enabling increasingly comprehensive analyses of how the microbiome changes over time , which are in turn beginning to provide insights into fundamental questions about microbiome dynamics 16 , 17 , 21 ., One of the most fundamental questions that still remains unanswered is to what degree the microbial community in the gut is deterministically dependent on its initial composition ( e . g . , microbial composition at birth ) ., More generally , it is unknown to what degree the microbial composition of the gut at a given time determines the microbial composition at a later time ., Additionally , there is only preliminary evidence of the long-term effects of early life events on the gut microbial community composition , and it is currently unclear whether these long-term effects traverse through a predefined set of potential trajectories 21 , 22 ., To address these questions , it is important to quantify the dependency of the microbial community at a given time on past community composition 23 , 24 ., This task has been previously studied in theoretical settings ., Specifically , the generalized Lotka-Volterra family of models infer changes in community composition through defined species-species or species-resource interaction terms , and are popular for describing internal ecological dynamics ., Recently , a few methods that rely on deterministic regularized model fitting using generalized Lotka-Volterra equations have been proposed ( e . g . , 25–27 ) ., Nonetheless , the importance of pure autoregressive factors ( a stochastic process in which future values are a function of the weighted sum of past values ) in driving gut microbial dynamics is , as yet , unclear ., Other approaches that utilize the full potential of longitudinal data , can often reveal insights about the autoregressive nature of the microbiome ., These include , for example , the sparse vector autoregression ( sVAR ) model , ( Gibbons et al . 24 ) , which assumes linear dynamics and is built around an autoregressive type of model , ARIMA Poisson ( Ridenhour et al . 28 ) , which assumes log-linear dynamics and suggests modeling the read counts along time using Poisson regression , and TGP-CODA ( Aijo et al . 2018 29 ) , which uses a Bayesian probabilistic model that combines a multinomial distribution with Gaussian processes ., Particularly , Gibbons et al . 24 , uses the sparse vector autoregression ( sVAR ) model to show evidence that the human gut microbial community has two dynamic regimes: autoregressive and non-autoregressive ., The autoregressive regime includes taxa that are affected by the community composition at previous time points , while the non-autoregressive regime includes taxa that their appearance in a specific time is random and or does not depend on the previous time points ., In this paper , we show that previous studies substantially underestimate the autoregressive component of the gut microbiome ., In order to quantify the dependency of taxa on past composition of the microbial community , we introduce Microbial community Temporal Variability Linear Mixed Model ( MTV-LMM ) , a ready-to-use scalable framework that can simultaneously identify and predict the dynamics of hundreds of time-dependent taxa across multiple hosts ., MTV-LMM is based on a linear mixed model , a heavily used tool in statistical genetics and other areas of genomics 30 , 31 ., Using MTV-LMM we introduce a novel concept we term ‘time-explainability’ , which corresponds to the fraction of temporal variance explained by the microbiome composition at previous time points ., Using time-explainability researchers can select the microorganisms whose abundance can be explained by the community composition at previous time points in a rigorous manner ., MTV-LMM has a few notable advantages ., First , unlike the sVAR model and the Bayesian approach proposed by Aijo et al . 29 , MTV-LMM models all the individual hosts simultaneously , thus leveraging the information across an entire population while adjusting for the host’s effect ( e . g , . host’s genetics or environment ) ., This provides MTV-LMM an increased power to detect temporal dependencies , as well as the ability to quantify the consistency of dynamics across individuals ., The Poisson regression method suggested by Ridenhour et al . 28 also utilizes the information from all individuals , but does not account for the individual effects , which may result in an inflated autoregressive component ., Second , MTV-LMM is computationally efficient , allowing it to model the dynamics of a complex ecosystem like the human gut microbiome by simultaneously evaluating the time-series of hundreds of taxa , across multiple hosts , in a timely manner ., Other methods , ( e . g . , TGP-CODA 29 , MDSINE 26 etc . ) can model only a small number of taxa ., Third , MTV-LMM can serve as a feature selection method , selecting only the taxa affected by the past composition of the microbiome ., The ability to identify these time-dependent taxa is crucial when fitting a time series model to study the microbial community temporal dynamics ., Finally , we demonstrate that MTV-LMM can serve as a standalone prediction model that outperforms commonly used models by an order of magnitude in predicting the taxa abundance ., We applied MTV-LMM to synthetic data , as suggested by Ajio et al . 2018 29 as well as to three real longitudinal studies of the gut microbiome ( David et al . 17 , Caporaso et al . 16 , and DIABIMMUNE 21 ) ., These datasets contain longitudinal abundance data using 16S rRNA gene sequencing ., Nonetheless , MTV-LMM is agnostic to the sequencing data type ( i . e . , 16s rRNA or shotgun sequencing ) ., Using MTV-LMM we find that in contrast to previous reports , a considerable portion of microbial taxa , in both infants and adults , display temporal structure that is predictable using the previous composition of the microbial community ., Moreover , we show that , on average , the time-explainability is an order of magnitude larger than previously estimated for these datasets ., We begin with an informal description of the main idea and utility of MTV-LMM ., A more comprehensive description can be found in the Methods ., MTV-LMM is motivated by our assumption that the temporal changes in the abundance of taxa are a time-homogeneous high-order Markov process ., MTV-LMM models the transitions of this Markov process by fitting a sequential linear mixed model ( LMM ) to predict the relative abundance of taxa at a given time point , given the microbial community composition at previous time points ., Intuitively , the linear mixed model correlates the similarity between the microbial community composition across different time points with the similarity of the taxa abundance at the next time points ., MTV-LMM is making use of two types of input data: ( 1 ) continuous relative abundance of focal taxa j at previous time points and ( 2 ) quantile-binned relative abundance of the rest of the microbial community at previous time points ., The output of MTV-LMM is prediction of continuous relative abundance , for each taxon , at future time points ., In order to apply linear mixed models , MTV-LMM generates a temporal kinship matrix , which represents the similarity between every pair of samples across time , where a sample is a normalization of taxa abundances at a given time point for a given individual ( see Methods ) ., When predicting the abundance of taxa j at time t , the model uses both the global state of the entire microbial community in the last q time points , as well as the abundance of taxa j in the previous p time points ., The parameters p and q are determined by the user , or can be determined using a cross-validation approach; a more formal description of their role is provided in the Methods ., MTV-LMM has the advantage of increased power due to a low number of parameters coupled with an inherent regularization mechanism , similar in essence to the widely used ridge regularization , which provides a natural interpretation of the model ., We evaluated MTV-LMM by testing its accuracy in predicting the abundance of taxa at a future time point using real time series data ., Such evaluation will mitigate overfitting , since the future data points are held out from the algorithm ., To measure accuracy on real data , we used the squared Pearson correlation coefficient between estimated and observed relative abundance along time , per taxon ., In addition we validated MTV-LMM using synthetic data , illustrating realistic dynamics and abundance distribution , as suggested by Aijo et al . 2018 29 ., Following 29 , we evaluate the performance of the model using the ‘estimation-error’ , defined to be the Euclidean distance between estimated and observed relative abundance , per time point ( see Supplementary Information S1 Note ) ., We used real time series data from three different datasets , each composed of longitudinal abundance data ., These three datasets are David et al . 17 ( 2 adult donors—DA , DB—average 250 time points per individual ) , Caporaso et al . 16 ( 2 adult donors—M3 , F4—average 231 time points per individual ) , and the DIABIMMUNE dataset 21 ( 39 infant donors—average 28 time points per individual ) ., In these datasets , the temporal parameters p and q were estimated using a validation set , and ranged from 0 to 3 ., See Methods for further details ., We compared the results of MTV-LMM to common approaches that are widely used for temporal microbiome modeling , namely the AR ( 1 ) model ( see Methods ) , the sparse vector autoregression model sVAR 24 , the ARIMA Poisson regression 28 and TGP-CODA 29 ., Overall , MTV-LMM’s prediction accuracy is higher than AR’s ( Supplementary Information S1 Table ) and significantly outperforms both the sVAR method and the Poisson regression across all datasets , using real time-series data ( Fig 1 ) ., In addition , since TGP-CODA can not be fully applied to these real datasets ( due to scalability limitations ) , we used synthetic data , considering a scenario of 200 taxa and 70 time points with realistic dynamics and abundance distribution , as suggested by the authors of this method ., Similarly to the real data , MTV-LMM significantly outperforms all the compared methods ( Supplementary Information S1 Fig ) ., We applied MTV-LMM to the DIABIMMUNE infant dataset and estimated the species-species association matrix across all individuals , using 1440 taxa that passed a preliminary screening according to temporal presence-absence patterns ( see Methods ) ., We found that most of these effects are close to zero , implying a sparse association pattern ., Next , we applied a principal component analysis ( PCA ) to the estimated species-species associations and found a strong phylogenetic structure ( PerMANOVA P-value = 0 . 001 ) suggesting that closely related species have similar association patterns within the microbial community ( Fig 2 ) ., These findings are supported by Thompson et al . 32 , who suggested that ecological interactions are phylogenetically conserved , where closely related species interact with similar partners ., Gomez et al . 33 tested these assumptions on a wide variety of hosts and found that generalized interactions can be evolutionary conserved ., We note that the association matrix estimated by MTV-LMM should be interpreted with caution since the number of possible associations is quadratic in the number of species , and it is , therefore , unfeasible to infer with high accuracy all the associations ., However , we can still aggregate information across species or higher taxonomic levels to uncover global patterns of the microbial composition dynamics ( e . g . , principal component analysis ) ., In order to address the fundamental question regarding the gut microbiota temporal variation , we quantify its autoregressive component ., Namely , we quantify to what degree the abundance of different taxa can be inferred based on the microbial community composition at previous time points ., In statistical genetics , the fraction of phenotypic variance explained by genetic factors is called heritability and is typically evaluated under an LMM framework 30 ., Intuitively , linear mixed models estimate heritability by measuring the correlation between the genetic similarity and the phenotypic similarity of pairs of individuals ., We used MTV-LMM to define an analogous concept that we term time-explainability , which corresponds to the fraction of temporal variance explained by the microbiome composition at previous time points ., In order to highlight the effect of the microbial community , we next estimated the time-explainability of taxa in each dataset , using the parameters q = 1 , p = 0 ., The resulting model corresponds to the formula: taxat = microbiome community ( t−1 ) + individual effect ( t−1 ) + unknown effects ., Of the taxa we examined , we identified a large portion of them to have a statistically significant time-explainability component across datasets ., Specifically , we found that over 85% of the taxa included in the temporal kinship matrix are significantly explained by the time-explainability component , with estimated time-explainability average levels of 23% in the DIABIMMUNE infant dataset ( sd = 15% ) , 21% in the Caporaso et al . ( 2011 ) dataset ( sd = 15% ) and 14% in the David el al . dataset ( sd = 10% ) ( Fig 3 , Supplementary Information S2 Fig ) ., Notably , we found that higher time explanability is associated with higher prediction accuracy ( Supplementary Information S3 Fig ) ., As a secondary analysis , we aggregated the time-explainability by taxonomic order , and found that in some orders ( non-autoregressive orders ) all taxa are non-autoregressive , while in others ( mixed orders ) we observed the presence of both autoregressive and non-autoregressive taxa ( Fig 4 , Supplementary Information S4 Fig ) , where an autoregressive taxa have a statistically significant time-explainability component ., Particularly , in the DIABIMMUNE infant data set , there are 7244 taxa , divided into 55 different orders ., However , the taxa recognized by MTV-LMM as autoregressive ( 1387 out of 7244 ) are represented in only 19 orders out of the 55 ., The remaining 36 orders do not include any autoregressive taxa ., Unlike the autoregressive organisms , these non-autoregressive organisms carry a strong phylogenetic structure ( t-test p-value < 10−16 ) , that may indicate a niche/habitat filtering ., This observation is consistent with the findings of Gibbons et al . 24 , who found a strong phylogenetic structure in the non-autoregressive organisms in the adult microbiome ., Notably , across all datasets , there is no significant correlation between the order dominance ( number of taxa in the order ) and the magnitude of its time-explainability component ( median Pearson r = 0 . 12 ) ., For example , in the DIABIMMUNE data set , the proportion of autoregressive taxa within the 19 mixed orders varies between 2% and 75% , where the average is approximately 20% ., In the most dominant order , Clostridiales ( representing 68% of the taxa ) , approximately 20% of the taxa are autoregressive and the average time-explainability is 23% ., In the second most dominant order , Bacteroidales , approximately 35% of the taxa are autoregressive and the average time-explainability is 31% ., In the Bifidobacteriales order , approximately 75% of the taxa are autoregressive , and the average time-explainability is 19% ( Fig 4 ) ., We hypothesize that the large fraction of autoregressive taxa in the Bifidobacteriales order , specifically in the infants dataset , can be partially attributed to the finding made by 34 , according to which some sub-species in this order appear to be specialized in the fermentation of human milk oligosaccharides and thus can be detected in infants but not in adults ., This emphasizes the ability of MTV-LMM to identify taxa that have prominent temporal dynamics that are both habitat and host-specific ., As an example of MTV-LMM’s ability to differentiate autoregressive from non-autoregressive taxa within the same order , we examined Burkholderiales , a relatively rare order ( less than 2% of the taxa in the data ) with 76 taxa overall , where only 19 of which were recognized as autoregressive by MTV-LMM ., Indeed , by examining the temporal behavior of each non-autoregressive taxa in this order , we witnessed abrupt changes in abundance over time , where the maximal number of consecutive time points with abundance greater than 0 is very small ., On the other hand , in the autoregressive taxa , we witnessed a consistent temporal behavior , where the maximal number of consecutive time points with abundance greater than 0 is well over 10 ( Supplementary Information S5 Fig ) ., The colonization of the human gut begins at birth and is characterized by a succession of microbial consortia 35–38 , where the diversity and richness of the microbiome reach adult levels in early childhood ., A longitudinal study has recently been used to show that infant gut microbiome begins transitioning towards an adult-like community after weaning 39 ., This observation is validated using our infant longitudinal data set ( DIABIMMUNE ) by applying PCA to the temporal kinship matrix ( Fig 5 ) ., Our analysis reveals that the first principal component ( accounting for 26% of the overall variability ) is associated with time ., Specifically , there is a clear clustering of the time samples from the first nine months of an infant’s life and the rest of the time samples ( months 10 − 36 ) which may be correlated to weaning ., As expected , we find a strong autoregressive component in an infant microbiome , which is highly associated with temporal variation across individuals ., By applying PCA to the temporal kinship matrix , we demonstrate that there is high similarity in the microbial community composition of infants at least in the first 9 months ., This similarity increases the power of our algorithm and thus helps MTV-LMM to detect autoregressive taxa ., In contrast to the infant microbiome , the adult microbiome is considered relatively stable 16 , 40 , but with considerable variation in the constituents of the microbial community between individuals ., Specifically , it was previously suggested that each individual adult has a unique gut microbial signature 41–43 , which is affected , among others factors , by environmental factors 20 and host lifestyle ( i . e . , antibiotics consumption , high-fat diets 17 etc . ) ., In addition , 17 showed that over the course of one year , differences between individuals were much larger than variation within individuals ., This observation was validated in our adult datasets ( David et al . and Caporaso et al . ) by applying PCA to the temporal kinship matrices ., In both David et al . and Caporaso et al . , the first principal component , which accounts for 61% and 43% of the overall variation respectively , is associated with the individual’s identity ( Fig 6 ) ., Using MTV-LMM we observed that despite the large similarity along time within adult individuals , there is also a non-negligible autoregressive component in the adult microbiome ., The fraction of variance explained by time across individuals can range from 6% up to 79% for different taxa ., These results shed more light on the temporal behavior of taxa in the adult microbiome , as opposed to that of infants , which are known to be highly affected by time 39 ., MTV-LMM uses a linear mixed model ( see 44 for a detailed review ) , a natural extension of standard linear regression , for the prediction of time series data ., We describe the technical details of the linear mixed model below ., We assume that the relative abundance levels of focal taxa j at time point t depend on a linear combination of the relative abundance levels of the microbial community at previous time points ., We further assume that temporal changes in relative abundance levels , in taxa j , are a time-homogeneous high-order Markov process ., We model the transitions of this Markov process using a linear mixed model , where we fit the p previous time points of taxa j as fixed effects and the q previous time points of the rest of the microbial community as random effects ., p and q are the temporal parameters of the model ., For simplicity of exposition , we present the generative linear mixed model that motivates the approach taken in MTV-LMM in two steps ., In the first step we model the microbial dynamics in one individual host ., In the second step we extend our model to N individuals , while accounting for the hosts’ effect ., We first describe the model assuming there is only one individual ., Consider a microbial community of m taxa measured at T equally spaced time points ., We get as input an m × T matrix M , where Mjt represents the relative-abundance levels of taxa j at time point t ., Let yj = ( Mj , p+1 , … , MjT ) t be a ( T −, p ) × 1 vector of taxa j relative abundance , across T − p time points starting at time point p + 1 and ending at time point T . Let Xj be a ( T −, p ) × ( p + 1 ) matrix of p + 1 covariates , comprised of an intercept vector as well as the first p time lags of taxa j ( i . e . , the relative abundance of taxa j in the p time points prior to the one predicted ) ., Formally , for k = 1 we have X t k j = 1 , and for 1 < k ≤ p + 1 we have X t k j = M j , t - k + 1 for t ≥ k ., For simplicity of exposition and to minimize the notation complexity , we assume for now that p = 1 . Let W be an ( T −, q ) × q ⋅ m normalized relative abundance matrix , representing the first q time lags of the microbial community ., For simplicity of exposition we describe the model in the case q = 1 , and then Wtj = Mjt ( in the more general case , we have Wtj = M⌈j/q⌉ , t− ( j mod, q ) , where p , q ≤ T − 1 ) ., With these notations , we assume the following linear model:, y j = X j β j + W u j + ϵ j , ( 1 ), where uj and ϵj are independent random variables distributed as uj∼ N ( 0 m , σ u j 2 I m ) and ϵ j ∼ N ( 0 T - 1 , σ ϵ j 2 I T - 1 ) ., The parameters of the model are βj ( fixed effects ) , σ u j 2 , and σ ϵ j 2 . We note that environmental factors known to be correlated with taxa abundance levels ( e . g . , diet , antibiotic usage 17 , 20 ) can be added to the model as fixed linear effects ( i . e . , added to the matrix Xj ) ., Given the high variability in the relative abundance levels , along with our desire to efficiently capture the effects of multiple taxa in the microbial community on each focal taxa j , we represent the microbial community input data ( matrix M ) using its quantiles ., Intuitively , we would like to capture the information as to whether a taxa is present or absent , or potentially introduce a few levels ( i . e . , high , medium , and low abundance ) ., To this end , we use the quantiles of each taxa to transform the matrix M into a matrix M ˜ , where M ˜ j t ∈ { 0 , 1 , 2 } depending on whether the abundance level is low ( below 25% quantile ) , medium , or high ( above 75% quantile ) ., We also tried other normalization strategies , including quantile normalization , which is typically used in gene expression eQTL analysis 45 , 46 , and the results were qualitatively similar ( see Supplementary Information S6 Fig ) ., We subsequently replace the matrix W by a matrix W ˜ , which is constructed analogously to W , but using M ˜ instead of M . Notably , both the fixed effect ( the relative abundance of yj at previous time points ) and the output of MTV-LMM are the continuous relative abundance ., The random effects are quantile-binned relative abundance of the rest of the microbial community at previous time points ( matrix W ˜ ) ., Thus , our model can now be described as, y j = X j β j + W ˜ u j + ϵ j ( 2 ) So far , we described the model assuming we have time series data from one individual ., We next extend the model to the case where time series data is available from multiple individuals ., In this case , we assume that the relative abundance levels of m taxa , denoted as the microbial community , have been measured at T time points across N individuals ., We assume the input consists of N matrices , M1 , … , MN , where matrix Mi corresponds to individual i , and it is of size m × T . Therefore , the outcome vector yj is now an n × 1 vector , composed of N blocks , where n = ( T − 1 ) N , and block i corresponds to the time points of individual, i . Formally , y k j = M j , ( k m o d ( T - 1 ) ) ⌈ k / ( T - 1 ) ⌉ ., Similarly , we define Xj and W ˜ as block matrices , with N different blocks , where corresponds to individual, i . When applied to multiple individuals , Model ( 2 ) may overfit to the individual effects ( e . g . , due to the host genetics and or environment ) ., In other words , since our goal is to model the changes in time , we need to condition these changes in time on the individual effects , that are unwanted confounders for our purposes ., We therefore construct a matrix H by randomly permuting the rows of each block matrix i in W ˜ , where the permutation is conducted only within the same individual ., Formally , we apply permutation πi ∈ ST−1 on the rows of each block matrix i , Mi , corresponding to individual i , where ST−1 is the set of all permutations of ( T − 1 ) elements ., In each πi , we are simultaneously permuting the entire microbial community ., Hence , matrix H corresponds to the data of each one of the individuals , but with no information about the time ( since the data was shuffled across the different time points ) ., With this addition , our final model is given by, y j = X j β j + W ˜ u j + H r + ϵ j , ( 3 ), where u j ∼ N ( 0 m , σ u j 2 I m ) and ϵ j ∼ N ( 0 n , σ ϵ j 2 I n ) , and r ∼ N ( 0 m , σ r 2 I m ) ., It is easy to verify that an equivalent mathematical representation of model 3 can be given by, y j ∼ N ( X j β j , σ A R j 2 K 1 + σ i n d 2 K 2 + σ ϵ j 2 I ) , ( 4 ), where σ A R j 2 = m σ u j 2 , K 1 = 1 m W ˜ W ˜ T , σ i n d 2 = m σ r 2 , K 2 = 1 m H H T . We will refer to K1 as the temporal kinship matrix , which represents the similarity between every pair of samples across time ( i . e . , represents the cross-correlation structure of the data ) ., We note that for the simplicity of exposition , we assumed so far that each sample has the same number of time points T , however in practice the number of samples may vary between the different individuals ., It is easy to extend the above model to the case where individual i has Ti time points , however the notations become cumbersome; the implementation of MTV-LMM , however takes into account a variable number of time points across the different individuals ., Once the distribution of yj is specified , one can proceed to estimate the fixed effects βj and the variance of the random effects using maximum likelihood approaches ., One common approach for estimating variance components is known as restricted maximum likelihood ( REML ) ., We followed the procedure described in the GCTA software package 47 , under ‘GREML analysis’ , originally developed for genotype data , and re-purposed it for longitudinal microbiome data ., GCTA implements the restricted maximum likelihood method via the average information ( AI ) algorithm ., Specifically , we performed a restricted maximum likelihood analysis using the function “–reml” followed by the option “–mgrm” ( reflects multiple variance components ) to estimate the variance explained by the microbial community at previous time points ., To predict the random effects by the BLUP ( best linear unbiased prediction ) method we use “–reml-pred-rand” ., This option is actually to predict the total temporal effect ( called “breeding value” in animal genetics ) of each time point attributed by the aggregated effect of the taxa used to estimate the temporal kinship matrix ., In both functions , to represent yj ( the abundance of taxa j at the next time point ) , we use the option “–pheno” ., For a detailed description see Supplementary Information S3 Note ., We define the term time-explainability , denoted as χ , to be the temporal variance explained by the microbial community in the previous time points ., Formally , for taxa j we define, χ j = σ A R j 2 σ A R j 2 + σ i n d 2 + σ ϵ j 2, The time-explainability was estimated with GCTA , using the temporal kinship matrix ., In order to measure the accuracy of time-explainability estimation , the average confidence interval width was estimated by computing the confidence interval widths for all autoregressive taxa and averaging the results ., Additionally , we adjust the time-explainability P-values for multiple comparisons using the Benjamini-Hochberg method 48 ., We now turn to the task of predicting y t j using the taxa abundance in time t − 1 ( or more generally in the last few time points ) ., Using our model notation , we are given xj and w ˜ , the covariates associated with a newly observed time point t in taxa j , and we would like to predict y t j with the greatest possible accuracy ., For a simple linear regression model , the answer is simply taking the covariate vector x and multiplying it by the estimated coefficients β ^ : y ^ t j = x T β ^ ., This practice yields unbiased estimates ., However , when attempting prediction in the linear mixed model case , things are not so simple ., One could adopt the same approach , but since the effects of the random components are not directly estimated , the vector of covariates w ˜ will not contribute directly to the predicted value of y t j , and will only affect the variance of the prediction , resulting in an unbiased but inefficient estimate ., Instead , one can use the correlation between the realized values of W ˜u , to attempt a better guess at the realization of w ˜ u for the new sample ., This is achieved by computing the distribution of the outcome of the new sample conditional on the full dataset , by using the following property of the multivariate 
Given the highly dynamic and complex nature of the human gut microbial community , the ability to identify and predict time-dependent compositional patterns of microbes is crucial to our understanding of the structure and functions of this ecosystem ., One factor that could affect such time-dependent patterns is microbial interactions , wherein community composition at a given time point affects the microbial composition at a later time point ., However , the field has not yet settled on the degree of this effect ., Specifically , it has been recently suggested that only a minority of taxa depend on the microbial composition in earlier times ., To address the issue of identifying and predicting temporal microbial patterns we developed a new model , MTV-LMM ( Microbial Temporal Variability Linear Mixed Model ) , a linear mixed model for the prediction of microbial community temporal dynamics ., MTV-LMM can identify time-dependent microbes ( i . e . , microbes whose abundance can be predicted based on the previous microbial composition ) in longitudinal studies , which can then be used to analyze the trajectory of the microbiome over time ., We evaluated the performance of MTV-LMM on real and synthetic time series datasets , and found that MTV-LMM outperforms commonly used methods for microbiome time series modeling ., Particularly , we demonstrate that the effect of the microbial composition in previous time points on the abundance of taxa at later time points is underestimated by a factor of at least 10 when applying previous approaches ., Using MTV-LMM , we demonstrate that a considerable portion of the human gut microbiome , both in infants and adults , has a significant time-dependent component that can be predicted based on microbiome composition in earlier time points ., This suggests that microbiome composition at a given time point is a major factor in defining future microbiome composition and that this phenomenon is considerably more common than previously reported for the human gut microbiome .
The ability to characterize and predict temporal trajectories of the microbial community in the human gut is crucial to our understanding of the structure and functions of this ecosystem ., In this study we develop MTV-LMM , a method for modeling time-series microbial community data ., Using MTV-LMM we find that in contrast to previous reports , a considerable portion of microbial taxa in both infants and adults display temporal structure that is predictable using the previous composition of the microbial community ., In reaching this conclusion we have adopted a number of concepts common in statistical genetics for use with longitudinal microbiome studies ., We introduce concepts such as time-explainability and the temporal kinship matrix , which we believe will be of use to other researchers studying microbial dynamics , through the framework of linear mixed models ., In particular we find that the association matrix estimated by MTV-LMM reveals known phylogenetic relationships and that the temporal kinship matrix uncovers known temporal structure in infant microbiome and inter-individual differences in adult microbiome ., Finally , we demonstrate that MTV-LMM significantly outperforms commonly used methods for temporal modeling of the microbiome , both in terms of its prediction accuracy as well as in its ability to identify time-dependent taxa .
taxonomy, children, ecology and environmental sciences, microbiome, community structure, statistics, microbiology, multivariate analysis, age groups, phylogenetics, data management, non-coding rna, mathematics, infants, cellular structures and organelles, microbial genomics, families, research and analysis methods, computer and information sciences, medical microbiology, mathematical and statistical techniques, principal component analysis, evolutionary systematics, ribosomes, people and places, community ecology, biochemistry, rna, ribosomal rna, cell biology, ecology, nucleic acids, genetics, biology and life sciences, population groupings, physical sciences, genomics, evolutionary biology, statistical methods