Robust estimation of microbial diversity in theory and in practice

Quantifying diversity is of central importance for the study of structure, function and evolution of microbial communities. The estimation of microbial diversity has received renewed attention with the advent of large-scale metagenomic studies. Here, we consider what the diversity observed in a sample tells us about the diversity of the community being sampled. First, we argue that one cannot reliably estimate the absolute and relative number of microbial species present in a community without making unsupported assumptions about species abundance distributions. The reason for this is that sample data do not contain information about the number of rare species in the tail of species abundance distributions. We illustrate the difficulty in comparing species richness estimates by applying Chao's estimator of species richness to a set of in silico communities: they are ranked incorrectly in the presence of large numbers of rare species. Next, we extend our analysis to a general family of diversity metrics ("Hill diversities"), and construct lower and upper estimates of diversity values consistent with the sample data. The theory generalizes Chao's estimator, which we retrieve as the lower estimate of species richness. We show that Shannon and Simpson diversity can be robustly estimated for the in silico communities. We analyze nine metagenomic data sets from a wide range of environments, and show that our findings are relevant for empirically-sampled communities. Hence, we recommend the use of Shannon and Simpson diversity rather than species richness in efforts to quantify and compare microbial diversity.Comment: To be published in The ISME Journal. Main text: 16 pages, 5 figures. Supplement: 16 pages, 4 figure

A model for transition of 5 '-nuclease domain of DNA polymerase I from inert to active modes

Bacteria contain DNA polymerase I (PolI), a single polypeptide chain consisting of similar to 930 residues, possessing DNA-dependent DNA polymerase, 3'-5' proofreading and 5'-3' exonuclease (also known as flap endonuclease) activities. PolI is particularly important in the processing of Okazaki fragments generated during lagging strand replication and must ultimately produce a double-stranded substrate with a nick suitable for DNA ligase to seal. PolI's activities must be highly coordinated both temporally and spatially otherwise uncontrolled 5'-nuclease activity could attack a nick and produce extended gaps leading to potentially lethal double-strand breaks. To investigate the mechanism of how PolI efficiently produces these nicks, we present theoretical studies on the dynamics of two possible scenarios or models. In one the flap DNA substrate can transit from the polymerase active site to the 5'-nuclease active site, with the relative position of the two active sites being kept fixed; while the other is that the 5'-nuclease domain can transit from the inactive mode, with the 5'-nuclease active site distant from the cleavage site on the DNA substrate, to the active mode, where the active site and substrate cleavage site are juxtaposed. The theoretical results based on the former scenario are inconsistent with the available experimental data that indicated that the majority of 5'-nucleolytic processing events are carried out by the same PolI molecule that has just extended the upstream primer terminus. By contrast, the theoretical results on the latter model, which is constructed based on available structural studies, are consistent with the experimental data. We thus conclude that the latter model rather than the former one is reasonable to describe the cooperation of the PolI's polymerase and 5'-3' exonuclease activities. Moreover, predicted results for the latter model are presented

Stromal Down-Regulation of Macrophage CD4/CCR5 Expression and NF-κB Activation Mediates HIV-1 Non-Permissiveness in Intestinal Macrophages

Tissue macrophages are derived exclusively from blood monocytes, which as monocyte-derived macrophages support HIV-1 replication. However, among human tissue macrophages only intestinal macrophages are non-permissive to HIV-1, suggesting that the unique microenvironment in human intestinal mucosa renders lamina propria macrophages non-permissive to HIV-1. We investigated this hypothesis using blood monocytes and intestinal extracellular matrix (stroma)-conditioned media (S-CM) to model the exposure of newly recruited monocytes and resident macrophages to lamina propria stroma, where the cells take up residence in the intestinal mucosa. Exposure of monocytes to S-CM blocked up-regulation of CD4 and CCR5 expression during monocyte differentiation into macrophages and inhibited productive HIV-1 infection in differentiated macrophages. Importantly, exposure of monocyte-derived macrophages simultaneously to S-CM and HIV-1 also inhibited viral replication, and sorted CD4+ intestinal macrophages, a proportion of which expressed CCR5+, did not support HIV-1 replication, indicating that the non-permissiveness to HIV-1 was not due to reduced receptor expression alone. Consistent with this conclusion, S-CM also potently inhibited replication of HIV-1 pseudotyped with vesicular stomatitis virus glycoprotein, which provides CD4/CCR5-independent entry. Neutralization of TGF-β in S-CM and recombinant TGF-β studies showed that stromal TGF-β inhibited macrophage nuclear translocation of NF-κB and HIV-1 replication. Thus, the profound inability of intestinal macrophages to support productive HIV-1 infection is likely the consequence of microenvironmental down-regulation of macrophage HIV-1 receptor/coreceptor expression and NF-κB activation

Placental syncytiotrophoblast constitutes a major barrier to vertical transmission of Listeria monocytogenes.

Listeria monocytogenes is an important cause of maternal-fetal infections and serves as a model organism to study these important but poorly understood events. L. monocytogenes can infect non-phagocytic cells by two means: direct invasion and cell-to-cell spread. The relative contribution of each method to placental infection is controversial, as is the anatomical site of invasion. Here, we report for the first time the use of first trimester placental organ cultures to quantitatively analyze L. monocytogenes infection of the human placenta. Contrary to previous reports, we found that the syncytiotrophoblast, which constitutes most of the placental surface and is bathed in maternal blood, was highly resistant to L. monocytogenes infection by either internalin-mediated invasion or cell-to-cell spread. Instead, extravillous cytotrophoblasts-which anchor the placenta in the decidua (uterine lining) and abundantly express E-cadherin-served as the primary portal of entry for L. monocytogenes from both extracellular and intracellular compartments. Subsequent bacterial dissemination to the villous stroma, where fetal capillaries are found, was hampered by further cellular and histological barriers. Our study suggests the placenta has evolved multiple mechanisms to resist pathogen infection, especially from maternal blood. These findings provide a novel explanation why almost all placental pathogens have intracellular life cycles: they may need maternal cells to reach the decidua and infect the placenta

PCB pollution continues to impact populations of orcas and other dolphins in European waters

Organochlorine (OC) pesticides and the more persistent polychlorinated biphenyls (PCBs) have well-established dose-dependent toxicities to birds, fish and mammals in experimental studies, but the actual impact of OC pollutants on European marine top predators remains unknown. Here we show that several cetacean species have very high mean blubber PCB concentrations likely to cause population declines and suppress population recovery. In a large pan-European meta-analysis of stranded (n = 929) or biopsied (n = 152) cetaceans, three out of four species:- striped dolphins (SDs), bottlenose dolphins (BNDs) and killer whales (KWs) had mean PCB levels that markedly exceeded all known marine mammal PCB toxicity thresholds. Some locations (e.g. western Mediterranean Sea, south-west Iberian Peninsula) are global PCB "hotspots" for marine mammals. Blubber PCB concentrations initially declined following a mid-1980s EU ban, but have since stabilised in UK harbour porpoises and SDs in the western Mediterranean Sea. Some small or declining populations of BNDs and KWs in the NE Atlantic were associated with low recruitment, consistent with PCB-induced reproductive toxicity. Despite regulations and mitigation measures to reduce PCB pollution, their biomagnification in marine food webs continues to cause severe impacts among cetacean top predators in European seas

Potentials of Mean Force for Protein Structure Prediction Vindicated, Formalized and Generalized

Understanding protein structure is of crucial importance in science, medicine and biotechnology. For about two decades, knowledge based potentials based on pairwise distances -- so-called "potentials of mean force" (PMFs) -- have been center stage in the prediction and design of protein structure and the simulation of protein folding. However, the validity, scope and limitations of these potentials are still vigorously debated and disputed, and the optimal choice of the reference state -- a necessary component of these potentials -- is an unsolved problem. PMFs are loosely justified by analogy to the reversible work theorem in statistical physics, or by a statistical argument based on a likelihood function. Both justifications are insightful but leave many questions unanswered. Here, we show for the first time that PMFs can be seen as approximations to quantities that do have a rigorous probabilistic justification: they naturally arise when probability distributions over different features of proteins need to be combined. We call these quantities reference ratio distributions deriving from the application of the reference ratio method. This new view is not only of theoretical relevance, but leads to many insights that are of direct practical use: the reference state is uniquely defined and does not require external physical insights; the approach can be generalized beyond pairwise distances to arbitrary features of protein structure; and it becomes clear for which purposes the use of these quantities is justified. We illustrate these insights with two applications, involving the radius of gyration and hydrogen bonding. In the latter case, we also show how the reference ratio method can be iteratively applied to sculpt an energy funnel. Our results considerably increase the understanding and scope of energy functions derived from known biomolecular structures

Metabolomics demonstrates divergent responses of two Eucalyptus species to water stress

Past studies of water stress in Eucalyptus spp. generally highlighted the role of fewer than five “important” metabolites, whereas recent metabolomic studies on other genera have shown tens of compounds are affected. There are currently no metabolite profiling data for responses of stress-tolerant species to water stress. We used GC–MS metabolite profiling to examine the response of leaf metabolites to a long (2 month) and severe (Ψpredawn < −2 MPa) water stress in two species of the perennial tree genus Eucalyptus (the mesic Eucalyptus pauciflora and the semi-arid Eucalyptus dumosa). Polar metabolites in leaves were analysed by GC–MS and inorganic ions by capillary electrophoresis. Pressure–volume curves and metabolite measurements showed that water stress led to more negative osmotic potential and increased total osmotically active solutes in leaves of both species. Water stress affected around 30–40% of measured metabolites in E. dumosa and 10–15% in E. pauciflora. There were many metabolites that were affected in E. dumosa but not E. pauciflora, and some that had opposite responses in the two species. For example, in E. dumosa there were increases in five acyclic sugar alcohols and four low-abundance carbohydrates that were unaffected by water stress in E. pauciflora. Re-watering increased osmotic potential and decreased total osmotically active solutes in E. pauciflora, whereas in E. dumosa re-watering led to further decreases in osmotic potential and increases in total osmotically active solutes. This experiment has added several extra dimensions to previous targeted analyses of water stress responses in Eucalyptus, and highlights that even species that are closely related (e.g. congeners) may respond differently to water stress and re-waterin

Arsenic as an Endocrine Disruptor: Arsenic Disrupts Retinoic Acid Receptor–and Thyroid Hormone Receptor–Mediated Gene Regulation and Thyroid Hormone–Mediated Amphibian Tail Metamorphosis

Background: Chronic exposure to excess arsenic in drinking water has been strongly associated with increased risks of multiple cancers, diabetes, heart disease, and reproductive and developmental problems in humans. We previously demonstrated that As, a potent endocrine disruptor at low, environmentally relevant levels, alters steroid signaling at the level of receptor-mediated gene regulation for all five steroid receptors. Objectives: The goal of this study was to determine whether As can also disrupt gene regulation via the retinoic acid (RA) receptor (RAR) and/or the thyroid hormone (TH) receptor (TR) and whether these effects are similar to previously observed effects on steroid regulation. Methods and results: Human embryonic NT2 or rat pituitary GH3 cells were treated with 0.01–5 μM sodium arsenite for 24 hr, with or without RA or TH, respectively, to examine effects of As on receptor-mediated gene transcription. At low, noncytotoxic doses, As significantly altered RAR-dependent gene transcription of a transfected RAR response element–luciferase construct and the native RA-inducible cytochrome P450 CYP26A gene in NT2 cells. Likewise, low-dose As significantly altered expression of a transfected TR response element–luciferase construct and the endogenous TR-regulated type I deiodinase (DIO1) gene in a similar manner in GH3 cells. An amphibian ex vivo tail metamorphosis assay was used to examine whether endocrine disruption by low-dose As could have specific pathophysiologic consequences, because tail metamorphosis is tightly controlled by TH through TR. TH-dependent tail shrinkage was inhibited in a dose-dependent manner by 0.1– 4.0 μM As. Conclusions: As had similar effects on RAR- and TR-mediated gene regulation as those previously observed for the steroid receptors, suggesting a common mechanism or action. Arsenic also profoundly affected a TR-dependent developmental process in a model animal system at very low concentrations. Because RAR and TH are critical for both normal human development and adult function and their dysregulation is associated with many disease processes, disruption of these hormone receptor–dependent processes by As is also potentially relevant to human developmental problems and disease risk

Does the lateral intercondylar ridge disappear in ACL deficient patients?

The aim of this study was to determine whether there is a difference in the presence of the lateral intercondylar ridge and the lateral bifurcate ridge between patients with sub-acute and chronic ACL injuries. We hypothesized that the ridges would be present less often with chronic ACL deficiency. Twenty-five patients with a chronic ACL injury were matched for age and gender to 25 patients with a sub-acute ACL injury. The lateral intercondylar ridge and lateral bifurcate ridge were scored as either present, absent, or indeterminate due to insufficient visualization by three blinded observers. The kappa for the three observers was .61 for the lateral intercondylar ridge and .58 for the lateral bifurcate ridge. The lateral intercondylar ridge was present in 88% of the sub-acute patients and 88% of the chronic patients. The lateral bifurcate ridge was present in 48% of the sub-acute and 48% of the chronic patients. This matched-pairs case–control study was unable to show a difference in the presence of the femoral bony ridges between patients with acute and chronic ACL injuries. The authors would suggest looking for the ridges as a landmark of the native ACL insertion site during ACL reconstruction in both acute and chronic ACL injuries

A coarse-to-fine approach to prostate boundary segmentation in ultrasound images

BACKGROUND: In this paper a novel method for prostate segmentation in transrectal ultrasound images is presented. METHODS: A segmentation procedure consisting of four main stages is proposed. In the first stage, a locally adaptive contrast enhancement method is used to generate a well-contrasted image. In the second stage, this enhanced image is thresholded to extract an area containing the prostate (or large portions of it). Morphological operators are then applied to obtain a point inside of this area. Afterwards, a Kalman estimator is employed to distinguish the boundary from irrelevant parts (usually caused by shadow) and generate a coarsely segmented version of the prostate. In the third stage, dilation and erosion operators are applied to extract outer and inner boundaries from the coarsely estimated version. Consequently, fuzzy membership functions describing regional and gray-level information are employed to selectively enhance the contrast within the prostate region. In the last stage, the prostate boundary is extracted using strong edges obtained from selectively enhanced image and information from the vicinity of the coarse estimation. RESULTS: A total average similarity of 98.76%(± 0.68) with gold standards was achieved. CONCLUSION: The proposed approach represents a robust and accurate approach to prostate segmentation