On the possible use of hydraulic force to assist with building the step pyramid of saqqara.

The Step Pyramid of Djoser in Saqqara, Egypt, is considered the oldest of the seven monumental pyramids built about 4,500 years ago. From transdisciplinary analysis, it was discovered that a hydraulic lift may have been used to build the pyramid. Based on our mapping of the nearby watersheds, we show that one of the unexplained massive Saqqara structures, the Gisr el-Mudir enclosure, has the features of a check dam with the intent to trap sediment and water. The topography beyond the dam suggests a possible ephemeral lake west of the Djoser complex and water flow inside the 'Dry Moat' surrounding it. In the southern section of the moat, we show that the monumental linear rock-cut structure consisting of successive, deep compartments combines the technical requirements of a water treatment facility: a settling basin, a retention basin, and a purification system. Together, the Gisr el-Mudir and the Dry Moat's inner south section work as a unified hydraulic system that improves water quality and regulates flow for practical purposes and human needs. Finally, we identified that the Step Pyramid's internal architecture is consistent with a hydraulic elevation mechanism never reported before. The ancient architects may have raised the stones from the pyramid centre in a volcano fashion using the sediment-free water from the Dry Moat's south section. Ancient Egyptians are famous for their pioneering and mastery of hydraulics through canals for irrigation purposes and barges to transport huge stones. This work opens a new line of research: the use of hydraulic force to erect the massive structures built by Pharaohs.

MODEL ASSISTED VARIABLE CLUSTERING: MINIMAX-OPTIMAL RECOVERY AND ALGORITHMS.

The problem of variable clustering is that of estimating groups of similar components of a p-dimensional vector X = (X 1, … , X p ) from n independent copies of X. There exists a large number of algorithms that return data-dependent groups of variables, but their interpretation is limited to the algorithm that produced them. An alternative is model-based clustering, in which one begins by defining population level clusters relative to a model that embeds notions of similarity. Algorithms tailored to such models yield estimated clusters with a clear statistical interpretation. We take this view here and introduce the class of G-block covariance models as a background model for variable clustering. In such models, two variables in a cluster are deemed similar if they have similar associations will all other variables. This can arise, for instance, when groups of variables are noise corrupted versions of the same latent factor. We quantify the difficulty of clustering data generated from a G-block covariance model in terms of cluster proximity, measured with respect to two related, but different, cluster separation metrics. We derive minimax cluster separation thresholds, which are the metric values below which no algorithm can recover the model-defined clusters exactly, and show that they are different for the two metrics. We therefore develop two algorithms, COD and PECOK, tailored to G-block covariance models, and study their minimax-optimality with respect to each metric. Of independent interest is the fact that the analysis of the PECOK algorithm, which is based on a corrected convex relaxation of the popular K-means algorithm, provides the first statistical analysis of such algorithms for variable clustering. Additionally, we compare our methods with another popular clustering method, spectral clustering. Extensive simulation studies, as well as our data analyses, confirm the applicability of our approach.

Capitalizing on opportunistic data for monitoring relative abundances of species.

With the internet, a massive amount of information on species abundance can be collected by citizen science programs. However, these data are often difficult to use directly in statistical inference, as their collection is generally opportunistic, and the distribution of the sampling effort is often not known. In this article, we develop a general statistical framework to combine such "opportunistic data" with data collected using schemes characterized by a known sampling effort. Under some structural assumptions regarding the sampling effort and detectability, our approach makes it possible to estimate the relative abundance of several species in different sites. It can be implemented through a simple generalized linear model. We illustrate the framework with typical bird datasets from the Aquitaine region in south-western France. We show that, under some assumptions, our approach provides estimates that are more precise than the ones obtained from the dataset with a known sampling effort alone. When the opportunistic data are abundant, the gain in precision may be considerable, especially for rare species. We also show that estimates can be obtained even for species recorded only in the opportunistic scheme. Opportunistic data combined with a relatively small amount of data collected with a known effort may thus provide access to accurate and precise estimates of quantitative changes in relative abundance over space and/or time.

Hybridization within Saccharomyces Genus Results in Homoeostasis and Phenotypic Novelty in Winemaking Conditions.

Despite its biotechnological interest, hybridization, which can result in hybrid vigor, has not commonly been studied or exploited in the yeast genus. From a diallel design including 55 intra- and interspecific hybrids between Saccharomyces cerevisiae and S. uvarum grown at two temperatures in enological conditions, we analyzed as many as 35 fermentation traits with original statistical and modeling tools. We first showed that, depending on the types of trait--kinetics parameters, life-history traits, enological parameters and aromas -, the sources of variation (strain, temperature and strain * temperature effects) differed in a large extent. Then we compared globally three groups of hybrids and their parents at two growth temperatures: intraspecific hybrids S. cerevisiae * S. cerevisiae, intraspecific hybrids S. uvarum * S. uvarum and interspecific hybrids S. cerevisiae * S. uvarum. We found that hybridization could generate multi-trait phenotypes with improved oenological performances and better homeostasis with respect to temperature. These results could explain why interspecific hybridization is so common in natural and domesticated yeast, and open the way to applications for wine-making.

The spatial distribution of Mustelidae in France.

We estimated the spatial distribution of 6 Mustelidae species in France using the data collected by the French national hunting and wildlife agency under the "small carnivorous species logbooks" program. The 1500 national wildlife protection officers working for this agency spend 80% of their working time traveling in the spatial area in which they have authority. During their travels, they occasionally detect dead or living small and medium size carnivorous animals. Between 2002 and 2005, each car operated by this agency was equipped with a logbook in which officers recorded information about the detected animals (species, location, dead or alive, date). Thus, more than 30000 dead or living animals were detected during the study period. Because a large number of detected animals in a region could have been the result of a high sampling pressure there, we modeled the number of detected animals as a function of the sampling effort to allow for unbiased estimation of the species density. For dead animals -- mostly roadkill -- we supposed that the effort in a given region was proportional to the distance traveled by the officers. For living animals, we had no way to measure the sampling effort. We demonstrated that it was possible to use the whole dataset (dead and living animals) to estimate the following: (i) the relative density -- i.e., the density multiplied by an unknown constant -- of each species of interest across the different French agricultural regions, (ii) the sampling effort for living animals for each region, and (iii) the relative detection probability for various species of interest.

Yeast proteome variations reveal different adaptive responses to grape must fermentation.

Saccharomyces cerevisiae and S. uvarum are two domesticated species of the Saccharomyces sensu stricto clade that diverged around 100 Ma after whole-genome duplication. Both have retained many duplicated genes associated with glucose fermentation and are characterized by the ability to achieve grape must fermentation. Nevertheless, these two species differ for many other traits, indicating that they underwent different evolutionary histories. To determine how the evolutionary histories of S. cerevisiae and S. uvarum are mirrored on the proteome, we analyzed the genetic variability of the proteomes of domesticated strains of these two species by quantitative mass spectrometry. Overall, 445 proteins were quantified. Massive variations of protein abundances were found, that clearly differentiated the two species. Abundance variations in specific metabolic pathways could be related to phenotypic traits known to discriminate the two species. In addition, proteins encoded by duplicated genes were shown to be differently recruited in each species. Comparing the strain differentiation based on the proteome variability to those based on the phenotypic and genetic variations further revealed that the strains of S. uvarum and some strains of S. cerevisiae displayed similar fermentative performances despite strong proteomic and genomic differences. Altogether, these results indicate that the ability of S. cerevisae and S. uvarum to complete grape must fermentation arose through different evolutionary roads, involving different metabolic pathways and duplicated genes.

Including shared peptides for estimating protein abundances: a significant improvement for quantitative proteomics.

Inferring protein abundances from peptide intensities is the key step in quantitative proteomics. The inference is necessarily more accurate when many peptides are taken into account for a given protein. Yet, the information brought by the peptides shared by different proteins is commonly discarded. We propose a statistical framework based on a hierarchical modeling to include that information. Our methodology, based on a simultaneous analysis of all the quantified peptides, handles the biological and technical errors as well as the peptide effect. In addition, we propose a practical implementation suitable for analyzing large data sets. Compared to a method based on the analysis of one protein at a time (that does not include shared peptides), our methodology proved to be far more reliable for estimating protein abundances and testing abundance changes. The source codes are available at http://pappso.inra.fr/bioinfo/all_p/.

Graph selection with GGMselect.

Applications on inference of biological networks have raised a strong interest in the problem of graph estimation in high-dimensional Gaussian graphical models. To handle this problem, we propose a two-stage procedure which first builds a family of candidate graphs from the data, and then selects one graph among this family according to a dedicated criterion. This estimation procedure is shown to be consistent in a high-dimensional setting, and its risk is controlled by a non-asymptotic oracle-like inequality. The procedure is tested on a real data set concerning gene expression data, and its performances are assessed on the basis of a large numerical study. The procedure is implemented in the R-package GGMselect available on the CRAN.