Molecular footprint of drug-selective pressure in a human immunodeficiency virus transmission chain

Known human immunodeficiency virus (HIV) transmission histories are invaluable models for investigating the evolutionary and transmission dynamics of the virus and to assess the accuracy of phylogenetic reconstructions. Here we have characterized an HIV-1 transmission chain consisting of nine infected patients, almost all of whom were treated with antiviral drugs at later stages of infection. Partial pol and env gp41 regions of the HIV genome were directly sequenced from plasma viral RNA for at least one sample from each patient. Phylogenetic analyses in pol using likelihood methods inferred an evolutionary history not fully compatible with the known transmission history. This could be attributed to parallel evolution of drug resistance mutations resulting in the incorrect clustering of multidrug-resistant virus. On the other hand, a fully compatible phylogenetic tree was reconstructed from the env sequences. We were able to identify and quantify the molecular footprint of drug-selective pressure in pol using maximum likelihood inference under different codon substitution models. An increased fixation rate of mutations in the HIV population of the multidrug-resistant patient was demonstrated using molecular clock modeling. We show that molecular evolutionary analyses, guided by a known transmission history, can reveal the presence of confounding factors like natural selection and caution should be taken when accurate descriptions of HIV evolution are required.status: publishe

Randomness Increases Order in Biological Evolution

n this text, we revisit part of the analysis of anti-entropy in Bailly and Longo (2009} and develop further theoretical reflections. In particular, we analyze how randomness, an essential component of biological variability, is associated to the growth of biological organization, both in ontogenesis and in evolution. This approach, in particular, focuses on the role of global entropy production and provides a tool for a mathematical understanding of some fundamental observations by Gould on the increasing phenotypic complexity along evolution. Lastly, we analyze the situation in terms of theoretical symmetries, in order to further specify the biological meaning of anti-entropy as well as its strong link with randomness

Gene expression profiles among murine strains segregate with distinct differences in the progression of radiation-induced lung disease.

Molecular mechanisms underlying development of acute pneumonitis and/or late fibrosis following thoracic irradiation remain poorly understood. Here, we hypothesize that heterogeneity in disease progression and phenotypic expression of radiation-induced lung disease (RILD) across murine strains presents an opportunity to better elucidate mechanisms driving tissue response toward pneumonitis and/or fibrosis. Distinct differences in disease progression were observed in age- and sex-matched CBA/J, C57L/J and C57BL/6J mice over 1 year after graded doses of whole-thorax lung irradiation (WTLI). Separately, comparison of gene expression profiles in lung tissue 24 h post-exposure demonstrated \u3e5000 genes to be differentially expressed (P\u3c0.01; \u3etwofold change) between strains with early versus late onset of disease. An immediate divergence in early tissue response between radiation-sensitive and -resistant strains was observed. In pneumonitis-prone C57L/J mice, differentially expressed genes were enriched in proinflammatory pathways, whereas in fibrosis-prone C57BL/6J mice, genes were enriched in pathways involved in purine and pyrimidine synthesis, DNA replication and cell division. At 24 h post-WTLI, different patterns of cellular damage were observed at the ultrastructural level among strains but microscopic damage was not yet evident under light microscopy. These data point toward a fundamental difference in patterns of early pulmonary tissue response to WTLI, consistent with the macroscopic expression of injury manifesting weeks to months after exposure. Understanding the mechanisms underlying development of RILD might lead to more rational selection of therapeutic interventions to mitigate healthy tissue damage

Predator and detritivore niche width helps to explain biocomplexity of experimental detritus-based food webs in four aquatic and terrestrial ecosystems

In the study of food webs, the existence and explanation of recurring patterns, such as the scale invariance of linkage density, predator–prey ratios and mean chain length, constitute long-standing issues. Our study focused on litter-associated food webs and explored the influence of detritivore and predator niche width (as d13C range) on web topological structure. To compare patterns within and between aquatic and terrestrial ecosystems and take account of intra-habitat variability, we constructed 42 macroinvertebrate patch-scale webs in four different habitats (lake, lagoon, beech forest and cornfield), using an experimental approach with litterbags. The results suggest that although web differences exist between ecosystems, patterns are more similar within than between aquatic and terrestrial web types. In accordance with optimal foraging theory, we found that the niche width of predators and prey increased with the number of predators and prey taxa as a proportion of total taxa in the community. The tendency was more marked in terrestrial ecosystems and can be explained by a lower per capita food level than in aquatic ecosystems, particularly evident for predators. In accordance with these results, the number of links increased with the number of species but with a significantly sharper regression slope for terrestrial ecosystems. As a consequence, linkage density, which was found to be directly correlated to niche width, increased with the total number of species in terrestrial webs, whereas it did not change significantly in aquatic ones, where connectance scaled negatively with the total number of species. In both types of ecosystem, web robustness to rare species removal increased with connectance and the niche width of predators. In conclusion, although limited to litter-associated macroinvertebrate assemblages, this study highlights structural differences and similarities between aquatic and terrestrial detrital webs, providing field evidence of the central role of niche width in determining the structure of detritus-based food webs and posing foraging optimisation constraints as a general mechanistic explanation of food web complexity differences within and between ecosystem types

Evaluation of genetic association of the INK4 locus with primary open angle glaucoma in East Indian population

INK4 locus at chromosome 9p21 has been reported to be associated with primary open angle glaucoma (POAG) and its subtypes along with the associated optic disc parameters across the populations of European, Japanese and African ancestries. The locus encodes three tumor suppressor genes namely CDKN2A, ARF, CDKN2B and a long non-coding RNA CDKN2B-AS1 (also known as ANRIL). Here, we report association study of 34 SNPs from INK4 locus with POAG in a population of Indo-European ancestry from the eastern part of India (350 patients and 354 controls). With 81% power to detect genetic association we observed only nominal association of rs1011970 (uncorrected p = 0.048) with POAG and rs10120688 (uncorrected p = 0.048) in patients without a high intra-ocular pressure (IOP<21 mm of Hg) compared to controls. This study, in contrast to the previous reports, suggests lack of significant genetic association of INK4 locus with POAG in East Indian population which needs to be replicated in larger studies in diverse world populations

OpenAlea: Scientific Workflows Combining Data Analysis and Simulation

International audienceAnalyzing biological data (e.g., annotating genomes, assembling NGS data...) may involve very complex and inter-linked steps where several tools are combined together. Scientific workflow systems have reached a level of maturity that makes them able to support the design and execution of such in-silico experiments, and thus making them increasingly popular in the bioinformatics community. However, in some emerging application domains such as system biology, developmental biology or ecology, the need for data analysis is combined with the need to model complex multi-scale biological systems, possibly involving multiple simulation steps. This requires the scientific work-flow to deal with retro-action to understand and predict the relationships between structure and function of these complex systems. OpenAlea (openalea.gforge.inria.fr) is the only scientific workflow system able to uniformly address the problem, which made it successful in the scientific community. One of its main originality is to introduce higher-order dataflows as a means to uniformly combine classical data analysis with modeling and simulation. In this demonstration paper, we provide for the first time the description of the OpenAlea system involving an original combination of features. We illustrate the demonstration on a high-throughput workflow in phenotyping, phenomics, and environmental control designed to study the interplay between plant architecture and climatic change

Specialized dynamical properties of promiscuous residues revealed by simulated conformational ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein-protein interaction prediction and design methods. © 2013 American Chemical Society

Species-level classification of the vaginal microbiome

Background The application of next-generation sequencing to the study of the vaginal microbiome is revealing the spectrum of microbial communities that inhabit the human vagina. High-resolution identification of bacterial taxa, minimally to the species level, is necessary to fully understand the association of the vaginal microbiome with bacterial vaginosis, sexually transmitted infections, pregnancy complications, menopause, and other physiological and infectious conditions. However, most current taxonomic assignment strategies based on metagenomic 16S rDNA sequence analysis provide at best a genus-level resolution. While surveys of 16S rRNA gene sequences are common in microbiome studies, few well-curated, body-site-specific reference databases of 16S rRNA gene sequences are available, and no such resource is available for vaginal microbiome studies. Results We constructed the Vaginal 16S rDNA Reference Database, a comprehensive and non-redundant database of 16S rDNA reference sequences for bacterial taxa likely to be associated with vaginal health, and we developed STIRRUPS, a new method that employs the USEARCH algorithm with a curated reference database for rapid species-level classification of 16S rDNA partial sequences. The method was applied to two datasets of V1-V3 16S rDNA reads: one generated from a mock community containing DNA from six bacterial strains associated with vaginal health, and a second generated from over 1,000 mid-vaginal samples collected as part of the Vaginal Human Microbiome Project at Virginia Commonwealth University. In both datasets, STIRRUPS, used in conjunction with the Vaginal 16S rDNA Reference Database, classified more than 95% of processed reads to a species-level taxon using a 97% global identity threshold for assignment. Conclusions This database and method provide accurate species-level classifications of metagenomic 16S rDNA sequence reads that will be useful for analysis and comparison of microbiome profiles from vaginal samples. STIRRUPS can be used to classify 16S rDNA sequence reads from other ecological niches if an appropriate reference database of 16S rDNA sequences is available