Taxonomic classification method for metagenomics based on core protein families with Core-Kaiju.

Characterizing species diversity and composition of bacteria hosted by biota is revolutionizing our understanding of the role of symbiotic interactions in ecosystems. Determining microbiomes diversity implies the assignment of individual reads to taxa by comparison to reference databases. Although computational methods aimed at identifying the microbe(s) taxa are available, it is well known that inferences using different methods can vary widely depending on various biases. In this study, we first apply and compare different bioinformatics methods based on 16S ribosomal RNA gene and shotgun sequencing to three mock communities of bacteria, of which the compositions are known. We show that none of these methods can infer both the true number of taxa and their abundances. We thus propose a novel approach, named Core-Kaiju, which combines the power of shotgun metagenomics data with a more focused marker gene classification method similar to 16S, but based on emergent statistics of core protein domain families. We thus test the proposed method on various mock communities and we show that Core-Kaiju reliably predicts both number of taxa and abundances. Finally, we apply our method on human gut samples, showing how Core-Kaiju may give more accurate ecological characterization and a fresh view on real microbiomes.

Upscaling Statistical Patterns from Reduced Storage in Social and Life Science Big Datasets.

Recent technological and computational advances have enabled the collection of data at an unprecedented rate. On the one hand, the large amount of data suddenly available has opened up new opportunities for new data-driven research but, on the other hand, it has brought into light new obstacles and challenges related to storage and analysis limits. Here, we strengthen an upscaling approach borrowed from theoretical ecology that allows us to infer with small errors relevant patterns of a dataset in its entirety, although only a limited fraction of it has been analysed. In particular we show that, after reducing the input amount of information on the system under study, by applying our framework it is still possible to recover two statistical patterns of interest of the entire dataset. Tested against big ecological, human activity and genomics data, our framework was successful in the reconstruction of global statistics related to both the number of types and their abundances while starting from limited presence/absence information on small random samples of the datasets. These results pave the way for future applications of our procedure in different life science contexts, from social activities to natural ecosystems.

The distance decay of similarity in tropical rainforests. A spatial point processes analytical formulation.

In this paper we are concerned with the analytical description of the change in floristic composition (species turnover) with the distance between two plots of a tropical rainforest due to the clustering of the individuals of the different species. We describe the plant arrangement by a superposition of spatial point processes and in this framework we introduce an analytical function which represents the average spatial density of the Sørensen similarity between two infinitesimal plots at distance r. We see that the decay in similarity with the distance is essentially described by the pair correlation function of the superposed process and that it is governed by the most abundant species. We test our analytical model with empirical data obtained for the Barro Colorado Island and Pasoh rainforests. To this end we adopt the statistical estimator for the pair correlation function in Shimatani (2001) and we design a novel one for the Sørensen similarity. Furthermore, we test our analytical formula by modeling the forest study area with Neyman-Scott point processes. We conclude comparing the advantages of our approach with other ones existing in literature.

Synthetic eco-evolutionary dynamics in simple molecular environment.

The understanding of eco-evolutionary dynamics, and in particular the mechanism of coexistence of species, is still fragmentary and in need of test bench model systems. To this aim we developed a variant of SELEX in vitro selection to study the evolution of a population of ∼1015 single-strand DNA oligonucleotide 'individuals'. We begin with a seed of random sequences which we select via affinity capture from ∼1012 DNA oligomers of fixed sequence ('resources') over which they compete. At each cycle ('generation'), the ecosystem is replenished via PCR amplification of survivors. Massive parallel sequencing indicates that across generations the variety of sequences ('species') drastically decreases, while some of them become populous and dominate the ecosystem. The simplicity of our approach, in which survival is granted by hybridization, enables a quantitative investigation of fitness through a statistical analysis of binding energies. We find that the strength of individual resource binding dominates the selection in the first generations, while inter- and intra-individual interactions become important in later stages, in parallel with the emergence of prototypical forms of mutualism and parasitism.

Seeing the forest for the trees through metabolic scaling.

We demonstrate that when power scaling occurs for an individual tree and in a forest, there is great resulting simplicity notwithstanding the underlying complexity characterizing the system over many size scales. Our scaling framework unifies seemingly distinct trends in a forest and provides a simple yet promising approach to quantitatively understand a bewilderingly complex many-body system with imperfectly known interactions. We show that the effective dimension, D tree , of a tree is close to 3, whereas a mature forest has D forest approaching 1. We discuss the energy equivalence rule and show that the metabolic rate-mass relationship is a power law with an exponent D/(D + 1) in both cases leading to a Kleiber's exponent of 3/4 for a tree and 1/2 for a forest. Our work has implications for understanding carbon sequestration and for climate science.

Upscaling human activity data: A statistical ecology approach.

Big data require new techniques to handle the information they come with. Here we consider four datasets (email communication, Twitter posts, Wikipedia articles and Gutenberg books) and propose a novel statistical framework to predict global statistics from random samples. More precisely, we infer the number of senders, hashtags and words of the whole dataset and how their abundances (i.e. the popularity of a hashtag) change through scales from a small sample of sent emails per sender, posts per hashtag and word occurrences. Our approach is grounded on statistical ecology as we map inference of human activities into the unseen species problem in biodiversity. Our findings may have applications to resource management in emails, collective attention monitoring in Twitter and language learning process in word databases.

Upscaling species richness and abundances in tropical forests.

The quantification of tropical tree biodiversity worldwide remains an open and challenging problem. More than two-fifths of the number of worldwide trees can be found either in tropical or in subtropical forests, but only ≈0.000067% of species identities are known. We introduce an analytical framework that provides robust and accurate estimates of species richness and abundances in biodiversity-rich ecosystems, as confirmed by tests performed on both in silico-generated and real forests. Our analysis shows that the approach outperforms other methods. In particular, we find that upscaling methods based on the log-series species distribution systematically overestimate the number of species and abundances of the rare species. We finally apply our new framework on 15 empirical tropical forest plots and quantify the minimum percentage cover that should be sampled to achieve a given average confidence interval in the upscaled estimate of biodiversity. Our theoretical framework confirms that the forests studied are comprised of a large number of rare or hyper-rare species. This is a signature of critical-like behavior of species-rich ecosystems and can provide a buffer against extinction.

Update on psychological functioning in adults with congenital heart disease: a systematic review.

The population of adults with congenital heart disease is increasing due to advancements in cardiology and cardiac surgery. Many patients face medical complications and psychosocial difficulties; however, it is not yet clear whether there is a direct relationship between medical status and the psychological functioning of these patients. This systematic review of the relevant literature is an attempt to: provide a comparison between the population of adults with congenital heart disease, the healthy reference population and similar cardiac populations when it comes to psychological functioning; explore the relationship between medical status/cardiac condition and psychological functioning; and identify the predictors of psychological distress in this population.