Translational Innovation in Medicine and Complexity

Summary Based on population genomic and environmental data, genomewide ecological association studies aim at detecting allele frequencies that exhibit significant statistical association with ecological gradients. Ecological association studies can provide lists of genetic polymorphisms that are potentially involved in local adaptation to environmental conditions through natural selection. Here, we present the R package LEA that enables users to run ecological association studies from the R command line. The package can perform analyses of population structure and genome scans for adaptive alleles from large genomic data sets. It derives advantages from R programming functionalities to adjust significance values for multiple testing issues and to visualize results. This note also illustrates the main steps of ecological association studies and the typical use of LEA for analysing data sets based on R commands.

The Escherichia coli species represents one of the best-studied model organisms, but also encompasses a variety of commensal and pathogenic strains that diversify by high rates of genetic change. We uniformly (re-) annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Within the approximately 18,000 families of orthologous genes, we found approximately 2,000 common to all strains. Although recombination rates are much higher than mutation rates, we show, both theoretically and using phylogenetic inference, that this does not obscure the phylogenetic signal, which places the B2 phylogenetic group and one group D strain at the basal position. Based on this phylogeny, we inferred past evolutionary events of gain and loss of genes, identifying functional classes under opposite selection pressures. We found an important adaptive role for metabolism diversification within group B2 and Shigella strains, but identified few or no extraintestinal virulence-specific genes, which could render difficult the development of a vaccine against extraintestinal infections. Genome flux in E. coli is confined to a small number of conserved positions in the chromosome, which most often are not associated with integrases or tRNA genes. Core genes flanking some of these regions show higher rates of recombination, suggesting that a gene, once acquired by a strain, spreads within the species by homologous recombination at the flanking genes. Finally, the genome's long-scale structure of recombination indicates lower recombination rates, but not higher mutation rates, at the terminus of replication. The ensuing effect of background selection and biased gene conversion may thus explain why this region is A+T-rich and shows high sequence divergence but low sequence polymorphism. Overall, despite a very high gene flow, genes co-exist in an organised genome.

DNA barcoding should provide rapid, accurate and automatable species identifications by using a standardized DNA region as a tag. Based on sequences available in GenBank and sequences produced for this study, we evaluated the resolution power of the whole chloroplast trnL (UAA) intron (254-767 bp) and of a shorter fragment of this intron (the P6 loop, 10-143 bp) amplified with highly conserved primers. The main limitation of the whole trnL intron for DNA barcoding remains its relatively low resolution (67.3% of the species from GenBank unambiguously identified). The resolution of the P6 loop is lower (19.5% identified) but remains higher than those of existing alternative systems. The resolution is much higher in specific contexts such as species originating from a single ecosystem, or commonly eaten plants. Despite the relatively low resolution, the whole trnL intron and its P6 loop have many advantages: the primers are highly conserved, and the amplification system is very robust. The P6 loop can even be amplified when using highly degraded DNA from processed food or from permafrost samples, and has the potential to be extensively used in food industry, in forensic science, in diet analyses based on feces and in ancient DNA studies.

The R package pcadapt performs genome scans to detect genes under selection based on population genomic data. It assumes that candidate markers are outliers with respect to how they are related to population structure. Because population structure is ascertained with principal component analysis, the package is fast and works with large-scale data. It can handle missing data and pooled sequencing data. By contrast to population-based approaches, the package handle admixed individuals and does not require grouping individuals into populations. Since its first release, pcadapt has evolved in terms of both statistical approach and software implementation. We present results obtained with robust Mahalanobis distance, which is a new statistic for genome scans available in the 2.0 and later versions of the package. When hierarchical population structure occurs, Mahalanobis distance is more powerful than the communality statistic that was implemented in the first version of the package. Using simulated data, we compare pcadapt to other computer programs for genome scans (BayeScan, hapflk, OutFLANK, sNMF). We find that the proportion of false discoveries is around a nominal false discovery rate set at 10% with the exception of BayeScan that generates 40% of false discoveries. We also find that the power of BayeScan is severely impacted by the presence of admixed individuals whereas pcadapt is not impacted. Last, we find that pcadapt and hapflk are the most powerful in scenarios of population divergence and range expansion. Because pcadapt handles next-generation sequencing data, it is a valuable tool for data analysis in molecular ecology.

Inference of individual ancestry coefficients, which is important for population genetic and association studies, is commonly performed using computer-intensive likelihood algorithms. With the availability of large population genomic data sets, fast versions of likelihood algorithms have attracted considerable attention. Reducing the computational burden of estimation algorithms remains, however, a major challenge. Here, we present a fast and efficient method for estimating individual ancestry coefficients based on sparse nonnegative matrix factorization algorithms. We implemented our method in the computer program sNMF and applied it to human and plant data sets. The performances of sNMF were then compared to the likelihood algorithm implemented in the computer program ADMIXTURE. Without loss of accuracy, sNMF computed estimates of ancestry coefficients with runtimes ∼10-30 times shorter than those of ADMIXTURE.

Adaptation to local environments often occurs through natural selection acting on a large number of loci, each having a weak phenotypic effect. One way to detect these loci is to identify genetic polymorphisms that exhibit high correlation with environmental variables used as proxies for ecological pressures. Here, we propose new algorithms based on population genetics, ecological modeling, and statistical learning techniques to screen genomes for signatures of local adaptation. Implemented in the computer program "latent factor mixed model" (LFMM), these algorithms employ an approach in which population structure is introduced using unobserved variables. These fast and computationally efficient algorithms detect correlations between environmental and genetic variation while simultaneously inferring background levels of population structure. Comparing these new algorithms with related methods provides evidence that LFMM can efficiently estimate random effects due to population history and isolation-by-distance patterns when computing gene-environment correlations, and decrease the number of false-positive associations in genome scans. We then apply these models to plant and human genetic data, identifying several genes with functions related to development that exhibit strong correlations with climatic gradients.

Plants generate effective responses to infection by recognizing both conserved and variable pathogen-encoded molecules. Pathogens deploy virulence effector proteins into host cells, where they interact physically with host proteins to modulate defense. We generated an interaction network of plant-pathogen effectors from two pathogens spanning the eukaryote-eubacteria divergence, three classes of Arabidopsis immune system proteins, and ~8000 other Arabidopsis proteins. We noted convergence of effectors onto highly interconnected host proteins and indirect, rather than direct, connections between effectors and plant immune receptors. We demonstrated plant immune system functions for 15 of 17 tested host proteins that interact with effectors from both pathogens. Thus, pathogens from different kingdoms deploy independently evolved virulence proteins that interact with a limited set of highly connected cellular hubs to facilitate their diverse life-cycle strategies.

Abstract On the basis of simulated data, this study compares the relative performances of the Bayesian clustering computer programsstructure,geneland,geneclustand a new program namedtess. While these four programs can detect population genetic structure from multilocus genotypes, only the last three ones include simultaneous analysis from geographical data. The programs are compared with respect to their abilities to infer the number of populations, to estimate membership probabilities, and to detect genetic discontinuities and clinal variation. The results suggest that combining analyses usingtessandstructureoffers a convenient way to address inference of spatial population structure.

Summary 1. Many recent statistical applications involve inference under complex models, where it is computationally prohibitive to calculate likelihoods but possible to simulate data. Approximate Bayesian computation (ABC) is devoted to these complex models because it bypasses the evaluation of the likelihood function by comparing observed and simulated data. 2. We introduce the R package ‘abc’ that implements several ABC algorithms for performing parameter estimation and model selection. In particular, the recently developed nonlinear heteroscedastic regression methods for ABC are implemented. The ‘abc’ package also includes a cross‐validation tool for measuring the accuracy of ABC estimates and to calculate the misclassification probabilities when performing model selection. The main functions are accompanied by appropriate summary and plotting tools. 3. R is already widely used in bioinformatics and several fields of biology. The R package ‘abc’ will make the ABC algorithms available to a large number of R users. ‘abc’ is a freely available R package under the GPL license, and it can be downloaded at http://cran.r‐project.org/web/packages/abc/index.html .

CGD is an immunodeficiency caused by deletions or mutations in genes that encode subunits of the leukocyte NADPH oxidase complex. Normally, assembly of the NADPH oxidase complex in phagosomes of certain phagocytic cells leads to a "respiratory burst", essential for the clearance of phagocytosed micro-organisms. CGD patients lack this mechanism, which leads to life-threatening infections and granuloma formation. However, a clear picture of the clinical course of CGD is hampered by its low prevalence (approximately 1:250,000). Therefore, extensive clinical data from 429 European patients were collected and analyzed. Of these patients 351 were males and 78 were females. X-linked (XL) CGD (gp91(phox) deficient) accounted for 67% of the cases, autosomal recessive (AR) inheritance for 33%. AR-CGD was diagnosed later in life, and the mean survival time was significantly better in AR patients (49.6 years) than in XL CGD (37.8 years), suggesting a milder disease course in AR patients. The disease manifested itself most frequently in the lungs (66% of patients), skin (53%), lymph nodes (50%), gastrointestinal tract (48%) and liver (32%). The most frequently cultured micro-organisms per episode were Staphylococcus aureus (30%), Aspergillus spp. (26%), and Salmonella spp. (16%). Surprisingly, Pseudomonas spp. (2%) and Burkholderia cepacia (<1%) were found only sporadically. Lesions induced by inoculation with BCG occurred in 8% of the patients. Only 71% of the patients received antibiotic maintenance therapy, and 53% antifungal prophylaxis. 33% were treated with gamma-interferon. 24 patients (6%) had received a stem cell transplantation. The most prominent reason of death was pneumonia and pulmonary abscess (18/84 cases), septicemia (16/84) and brain abscess (4/84). These data provide further insight in the clinical course of CGD in Europe and hopefully can help to increase awareness and optimize the treatment of these patients.

Epistatic interactions between mutations play a prominent role in evolutionary theories. Many studies have found that epistasis is widespread, but they have rarely considered beneficial mutations. We analyzed the effects of epistasis on fitness for the first five mutations to fix in an experimental population of Escherichia coli. Epistasis depended on the effects of the combined mutations--the larger the expected benefit, the more negative the epistatic effect. Epistasis thus tended to produce diminishing returns with genotype fitness, although interactions involving one particular mutation had the opposite effect. These data support models in which negative epistasis contributes to declining rates of adaptation over time. Sign epistasis was rare in this genome-wide study, in contrast to its prevalence in an earlier study of mutations in a single gene.

The Mediterranean diet has been linked to a number of health benefits, including reduced mortality risk and lower incidence of cardiovascular disease. Definitions of the Mediterranean diet vary across some settings, and scores are increasingly being employed to define Mediterranean diet adherence in epidemiological studies. Some components of the Mediterranean diet overlap with other healthy dietary patterns, whereas other aspects are unique to the Mediterranean diet. In this forum article, we asked clinicians and researchers with an interest in the effect of diet on health to describe what constitutes a Mediterranean diet in different geographical settings, and how we can study the health benefits of this dietary pattern.

Fall detection of the elderly is a major public health problem. Thus it has generated a wide range of applied research and prompted the development of telemonitoring systems to enable the early diagnosis of fall conditions. This article is a survey of systems, algorithms and sensors, for the automatic early detection of the fall of elderly persons. It points out the difficulty to compare the performances of the different systems due to the lack of a common framework. It then proposes a procedure for this evaluation.

Enzymatic fuel cells use enzymes to produce energy from bioavailable substrates. However, such biofuel cells are limited by the difficult electrical wiring of enzymes to the electrode. Here we show the efficient wiring of enzymes in a conductive pure carbon nanotube matrix for the fabrication of a glucose biofuel cell (GBFC). Glucose oxidase and laccase were respectively incorporated in carbon nanotube disks by mechanical compression. The characterization of each bioelectrode shows an open circuit potential corresponding to the redox potential of the respective enzymes, and high current densities for glucose oxidation and oxygen reduction. The mediatorless GBFC delivers a high power density up to 1.3 mW cm−2 and an open circuit voltage of 0.95 V. Moreover, the GBFC remains stable for 1 month and delivers 1 mW cm−2 power density under physiological conditions (5×10−3 mol l−1 glucose, pH 7). To date, these values are the best performances obtained for a GBFC. Glucose biofuel cells can be used to produce clean energy from renewable sources, but their use is limited by poor stability and low power output. In this study, bioelectrodes are fabricated using carbon nanotubes and the resulting biofuel cells have improved stability and power.

Thierry Wirth, Philip Supply, Stefan Niemann and colleagues analyze 4,987 Mycobacterium tuberculosis strains of the Beijing lineage isolated from 99 countries. They report whole-genome sequencing of 110 representative strains, characterize global population structure and reconstruct the evolutionary history of this lineage. Mycobacterium tuberculosis strains of the Beijing lineage are globally distributed and are associated with the massive spread of multidrug-resistant (MDR) tuberculosis in Eurasia. Here we reconstructed the biogeographical structure and evolutionary history of this lineage by genetic analysis of 4,987 isolates from 99 countries and whole-genome sequencing of 110 representative isolates. We show that this lineage initially originated in the Far East, from where it radiated worldwide in several waves. We detected successive increases in population size for this pathogen over the last 200 years, practically coinciding with the Industrial Revolution, the First World War and HIV epidemics. Two MDR clones of this lineage started to spread throughout central Asia and Russia concomitantly with the collapse of the public health system in the former Soviet Union. Mutations identified in genes putatively under positive selection and associated with virulence might have favored the expansion of the most successful branches of the lineage.

BACKGROUND: The authors conducted a patient-based survey of practices to fully describe the assessment and the management of pain and sedation of a large cohort of mechanically ventilated patients during their first week of intensive care unit (ICU) stay. METHODS: A total of 1,381 adult patients were included in a prospective, observational study in 44 ICUs in France. Pain and sedation assessment, analgesic and sedative use, and analgesic management during procedural pain were collected on days 2, 4, and 6 of the ICU stay. RESULTS: The observed rates of assessment on day 2 for sedation (43%) and analgesia (42%) were significantly smaller than that of use of sedatives (72%) and opioids (90%), also noted on days 4 and 6. The use of protocols/guidelines for sedation/analgesia in the ICU reduced the proportion of patients who were treated, although not evaluated. A large proportion of assessed patients were in a deep state of sedation (40-50%). Minor changes in the dosages of the main prescribed agents for sedation (midazolam, propofol) and analgesia (sufentanil, fentanyl, morphine, remifentanil) were found across 6 days of the patient's ICU stay. Procedural pain was specifically managed for less than 25% of patients; during those procedures, the proportion of patients with pain significantly increased from the baseline pain evaluation. CONCLUSIONS: Excessively deep states of sedation and a lack of analgesia during painful procedures must be prevented. To facilitate systematic pain and sedation assessment and to adjust daily drug dosages accordingly, it seems crucial to promote educational programs and elaboration of protocols/guidelines in the ICU.

In contrast to the NADPH oxidases Nox1 and Nox2, which generate superoxide (O(2)(·-)), Nox4 produces hydrogen peroxide (H(2)O(2)). We constructed chimeric proteins and mutants to address the protein region that specifies which reactive oxygen species is produced. Reactive oxygen species were measured with luminol/horseradish peroxidase and Amplex Red for H(2)O(2) versus L-012 and cytochrome c for O(2)(·-). The third extracytosolic loop (E-loop) of Nox4 is 28 amino acids longer than that of Nox1 or Nox2. Deletion of E-loop amino acids only present in Nox4 or exchange of the two cysteines in these stretches switched Nox4 from H(2)O(2) to O(2)(·-) generation while preserving expression and intracellular localization. In the presence of an NO donor, the O(2)()-producing Nox4 mutants, but not wild-type Nox4, generated peroxynitrite, excluding artifacts of the detection system as the apparent origin of O(2)(·-). In Cos7 cells, in which Nox4 partially localizes to the plasma membrane, an antibody directed against the E-loop decreased H(2)O(2) but increased O(2)(·-) formation by Nox4 without affecting Nox1-dependent O(2)(·-) formation. The E-loop of Nox4 but not Nox1 and Nox2 contains a highly conserved histidine that could serve as a source for protons to accelerate spontaneous dismutation of superoxide to form H(2)O(2). Mutation of this but not of four other conserved histidines also switched Nox4 from H(2)O(2) to O(2)(·-) formation. Thus, H(2)O(2) formation is an intrinsic property of Nox4 that involves its E-loop. The structure of the E-loop may hinder O(2)(·-) egress and/or provide a source for protons, allowing dismutation to form H(2)O(2).

Most source separation algorithms are based on a model of stationary sources. However, it is a simple matter to take advantage of possible nonstationarities of the sources to achieve separation. This paper develops novel approaches in this direction based on the principles of maximum likelihood and minimum mutual information. These principles are exploited by efficient algorithms in both the off-line case (via a new joint diagonalization procedure) and in the on-line case (via a Newton-like procedure). Some experiments showing the good performance of our algorithms and evidencing an interesting feature of our methods are presented: their ability to achieve a kind of super-efficiency. The paper concludes with a discussion contrasting separating methods for non-Gaussian and nonstationary models and emphasizing that, as a matter of fact, "what makes the algorithms work" is-strictly speaking-not the nonstationarity itself but rather the property that each realization of the source signals has a time-varying envelope.

The history of click-speaking Khoe-San, and African populations in general, remains poorly understood. We genotyped ~2.3 million single-nucleotide polymorphisms in 220 southern Africans and found that the Khoe-San diverged from other populations ≥100,000 years ago, but population structure within the Khoe-San dated back to about 35,000 years ago. Genetic variation in various sub-Saharan populations did not localize the origin of modern humans to a single geographic region within Africa; instead, it indicated a history of admixture and stratification. We found evidence of adaptation targeting muscle function and immune response; potential adaptive introgression of protection from ultraviolet light; and selection predating modern human diversification, involving skeletal and neurological development. These new findings illustrate the importance of African genomic diversity in understanding human evolutionary history.

By 2050, about one third of the French population will be over 65. Our laboratory's current research focuses on the monitoring of elderly people at home, to detect a loss of autonomy as early as possible. Our aim is to quantify criteria such as the international activities of daily living (ADL) or the French Autonomie Gerontologie Groupes Iso-Ressources (AGGIR) scales, by automatically classifying the different ADL performed by the subject during the day. A Health Smart Home is used for this. Our Health Smart Home includes, in a real flat, infrared presence sensors (location), door contacts (to control the use of some facilities), temperature and hygrometry sensor in the bathroom, and microphones (sound classification and speech recognition). A wearable kinematic sensor also informs postural transitions (using pattern recognition) and walk periods (frequency analysis). This data collected from the various sensors are then used to classify each temporal frame into one of the ADL that was previously acquired (seven activities: hygiene, toilet use, eating, resting, sleeping, communication, and dressing/undressing). This is done using support vector machines. We performed a 1-h experimentation with 13 young and healthy subjects to determine the models of the different activities, and then we tested the classification algorithm (cross validation) with real data.