Variability and bias in microbiome metagenomic sequencing: an interlaboratory study comparing experimental protocols

Several studies have documented the significant impact of methodological choices in microbiome analyses. The myriad of methodological options available complicate the replication of results and generally limit the comparability of findings between independent studies that use differing techniques and measurement pipelines. Here we describe the Mosaic Standards Challenge (MSC), an international interlaboratory study designed to assess the impact of methodological variables on the results. The MSC did not prescribe methods but rather asked participating labs to analyze 7 shared reference samples (5 × human stool samples and 2 × mock communities) using their standard laboratory methods. To capture the array of methodological variables, each participating lab completed a metadata reporting sheet that included 100 different questions regarding the details of their protocol. The goal of this study was to survey the methodological landscape for microbiome metagenomic sequencing (MGS) analyses and the impact of methodological decisions on metagenomic sequencing results. A total of 44 labs participated in the MSC by submitting results (16S or WGS) along with accompanying metadata; thirty 16S rRNA gene amplicon datasets and 14 WGS datasets were collected. The inclusion of two types of reference materials (human stool and mock communities) enabled analysis of both MGS measurement variability between different protocols using the biologically-relevant stool samples, and MGS bias with respect to ground truth values using the DNA mixtures. Owing to the compositional nature of MGS measurements, analyses were conducted on the ratio of Firmicutes: Bacteroidetes allowing us to directly apply common statistical methods. The resulting analysis demonstrated that protocol choices have significant effects, including both bias of the MGS measurement associated with a particular methodological choices, as well as effects on measurement robustness as observed through the spread of results between labs making similar methodological choices. In the analysis of the DNA mock communities, MGS measurement bias was observed even when there was general consensus among the participating laboratories. This study was the result of a collaborative effort that included academic, commercial, and government labs. In addition to highlighting the impact of different methodological decisions on MGS result comparability, this work also provides insights for consideration in future microbiome measurement study design.

A novel sequencing-based vaginal health assay combining self-sampling, HPV detection and genotyping, STI detection, and vaginal microbiome analysis

The composition of the vaginal microbiome, including both the presence of pathogens involved in sexually transmitted infections (STI) as well as commensal microbiota, has been shown to have important associations for a woman's reproductive and general health. Currently, healthcare providers cannot offer comprehensive vaginal microbiome screening, but are limited to the detection of individual pathogens, such as high-risk human papillomavirus (hrHPV), the predominant cause of cervical cancer. There is no single test on the market that combines HPV, STI, and microbiome screening. Here, we describe a novel inclusive vaginal health assay that combines self-sampling with sequencing-based HPV detection and genotyping, vaginal microbiome analysis, and STI-associated pathogen detection. The assay includes genotyping and detection of 14 hrHPV types, 5 low-risk HPV types (lrHPV), as well as the relative abundance of 31 bacterial taxa of clinical importance, including Lactobacillus, Sneathia, Gardnerella, and 3 pathogens involved in STI, with high sensitivity, specificity, and reproducibility. For each of these taxa, reference ranges were determined in a group of 50 self-reported healthy women. The HPV sequencing portion of the test was evaluated against the digene High-Risk HPV HC2 DNA test. For hrHPV genotyping, agreement was 95.3% with a kappa of 0.804 (601 samples); after removal of samples in which the digene hrHPV probe showed cross-reactivity with lrHPV types, the sensitivity and specificity of the hrHPV genotyping assay were 94.5% and 96.6%, respectively, with a kappa of 0.841. For lrHPV genotyping, agreement was 93.9% with a kappa of 0.788 (148 samples), while sensitivity and specificity were 100% and 92.9%, respectively. This novel assay could be used to complement conventional cervical cancer screening, because its self-sampling format can expand access among women who would otherwise not participate, and because of its additional information about the composition of the vaginal microbiome and the presence of pathogens.

VarQ: A Tool for the Structural and Functional Analysis of Human Protein Variants

Understanding the functional effect of Single Amino acid Substitutions (SAS), derived from the occurrence of single nucleotide variants (SNVs), and their relation to disease development is a major issue in clinical genomics. Despite the existence of several bioinformatic algorithms and servers that predict if a SAS is pathogenic or not, they give little or no information at all on the reasons for pathogenicity prediction and on the actual predicted effect of the SAS on the protein function. Moreover, few actual methods take into account structural information when available for automated analysis. Moreover, many of these algorithms are able to predict an effect that no necessarily translates directly into pathogenicity. VarQ is a bioinformatic pipeline that incorporates structural information for the detailed analysis and prediction of SAS effect on protein function. It is an online tool which uses UniProt id and automatically analyzes known and user provided SAS for their effect on protein activity, folding, aggregation and protein interactions, among others. We show that structural information, when available, can improve the SAS pathogenicity diagnosis and more important explain its causes. We show that VarQ is able to correctly reproduce previous analysis of RASopathies related mutations, saving extensive and time consuming manual curation. VarQ assessment was performed over a set of previously manually curated RASopathies (diseases that affects the RAS/MAPK signaling pathway) related variants, showing its ability to correctly predict the phenotypic outcome and its underlying cause. This resource is available online at http://varq.qb.fcen.uba.ar/. Supporting Information & Tutorials may be found in the webpage of the tool.

Next generation sequencing panel based on single molecule molecular inversion probes for detecting genetic variants in children with hypopituitarism

BACKGROUND

Congenital Hypopituitarism is caused by genetic and environmental factors. Over 30 genes have been implicated in isolated and/or combined pituitary hormone deficiency. The etiology remains unknown for up to 80% of the patients, but most cases have been analyzed by limited candidate gene screening. Mutations in the PROP1 gene are the most common known cause, and the frequency of mutations in this gene varies greatly by ethnicity. We designed a custom array to assess the frequency of mutations in known hypopituitarism genes and new candidates, using single molecule molecular inversion probes sequencing (smMIPS).

METHODS

We used this panel for the first systematic screening for causes of hypopituitarism in children. Molecular inversion probes were designed to capture 693 coding exons of 30 known genes and 37 candidate genes. We captured genomic DNA from 51 pediatric patients with CPHD (n = 43) or isolated GH deficiency (IGHD) (n = 8) and their parents and conducted next generation sequencing.

RESULTS

We obtained deep coverage over targeted regions and demonstrated accurate variant detection by comparison to whole-genome sequencing in a control individual. We found a dominant mutation GH1, p.R209H, in a three-generation pedigree with IGHD.

CONCLUSIONS

smMIPS is an efficient and inexpensive method to detect mutations in patients with hypopituitarism, drastically limiting the need for screening individual genes by Sanger sequencing.

A quantitative model for oxygen uptake and release in a family of hemeproteins

MOTIVATION

Hemeproteins have many diverse functions that largely depend on the rate at which they uptake or release small ligands, like oxygen. These proteins have been extensively studied using either simulations or experiments, albeit only qualitatively and one or two proteins at a time.

RESULTS

We present a physical-chemical model, which uses data obtained exclusively from computer simulations, to describe the uptake and release of oxygen in a family of hemeproteins, called truncated hemoglobins (trHbs). Through a rigorous statistical analysis we demonstrate that our model successfully recaptures all the reported experimental oxygen association and dissociation kinetic rate constants, thus allowing us to establish the key factors that determine the rates at which these hemeproteins uptake and release oxygen. We found that internal tunnels as well as the distal site water molecules control ligand uptake, whereas oxygen stabilization by distal site residues controls ligand release. Because these rates largely determine the functions of these hemeproteins, these approaches will also be important tools in characterizing the trHbs members with unknown functions.

CONTACT

lboechi@ic.fcen.uba.ar

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evolutionary and Functional Relationships in the Truncated Hemoglobin Family

Predicting function from sequence is an important goal in current biological research, and although, broad functional assignment is possible when a protein is assigned to a family, predicting functional specificity with accuracy is not straightforward. If function is provided by key structural properties and the relevant properties can be computed using the sequence as the starting point, it should in principle be possible to predict function in detail. The truncated hemoglobin family presents an interesting benchmark study due to their ubiquity, sequence diversity in the context of a conserved fold and the number of characterized members. Their functions are tightly related to O2 affinity and reactivity, as determined by the association and dissociation rate constants, both of which can be predicted and analyzed using in-silico based tools. In the present work we have applied a strategy, which combines homology modeling with molecular based energy calculations, to predict and analyze function of all known truncated hemoglobins in an evolutionary context. Our results show that truncated hemoglobins present conserved family features, but that its structure is flexible enough to allow the switch from high to low affinity in a few evolutionary steps. Most proteins display moderate to high oxygen affinities and multiple ligand migration paths, which, besides some minor trends, show heterogeneous distributions throughout the phylogenetic tree, again suggesting fast functional adaptation. Our data not only deepens our comprehension of the structural basis governing ligand affinity, but they also highlight some interesting functional evolutionary trends.

The N-terminal pre-A region of Mycobacterium tuberculosis 2/2HbN promotes NO-dioxygenase activity

A unique defense mechanisms by which Mycobacterium tuberculosis protects itself from nitrosative stress is based on the O2 -dependent NO-dioxygenase (NOD) activity of truncated hemoglobin 2/2HbN (Mt2/2HbN). The NOD activity largely depends on the efficiency of ligand migration to the heme cavity through a two-tunnel (long and short) system; recently, it was also correlated with the presence at the Mt2/2HbN N-terminus of a short pre-A region, not conserved in most 2/2HbNs, whose deletion results in a drastic reduction of NO scavenging. In the present study, we report the crystal structure of Mt2/2HbN-ΔpreA, lacking the pre-A region, at a resolution of 1.53 Å. We show that removal of the pre-A region results in long range effects on the protein C-terminus, promoting the assembly of a stable dimer, both in the crystals and in solution. In the Mt2/2HbN-ΔpreA dimer, access of heme ligands to the short tunnel is hindered. Molecular dynamics simulations show that the long tunnel branch is the only accessible pathway for O2 -ligand migration to/from the heme, and that the gating residue Phe(62)E15 partly restricts the diameter of the tunnel. Accordingly, kinetic measurements indicate that the kon value for peroxynitrite isomerization by Mt2/2HbN-ΔpreA-Fe(III) is four-fold lower relative to the full-length protein, and that NO scavenging by Mt2/2HbN-ΔpreA-Fe(II)-O2 is reduced by 35-fold. Therefore, we speculate that Mt2/2HbN evolved to host the pre-A region as a mechanism for preventing dimerization, thus reinforcing the survival of the microorganism against the reactive nitrosative stress in macrophages.

DATABASE

Coordinates and structure factors have been deposited in the Protein Data Bank under accession number 5AB8.

Ligand uptake modulation by internal water molecules and hydrophobic cavities in hemoglobins

Internal water molecules play an active role in ligand uptake regulation, since displacement of retained water molecules from protein surfaces or cavities by incoming ligands can promote favorable or disfavorable effects over the global binding process. Detection of these water molecules by X-ray crystallography is difficult given their positional disorder and low occupancy. In this work, we employ a combination of molecular dynamics simulations and ligand rebinding over a broad time range to shed light into the role of water molecules in ligand migration and binding. Computational studies on the unliganded structure of the thermostable truncated hemoglobin from Thermobifida fusca (Tf-trHbO) show that a water molecule is in the vicinity of the iron heme, stabilized by WG8 with the assistance of YCD1, exerting a steric hindrance for binding of an exogenous ligand. Mutation of WG8 to F results in a significantly lower stabilization of this water molecule and in subtle dynamical structural changes that favor ligand binding, as observed experimentally. Water is absent from the fully hydrophobic distal cavity of the triple mutant YB10F-YCD1F-WG8F (3F), due to the lack of residues capable of stabilizing it nearby the heme. In agreement with these effects on the barriers for ligand rebinding, over 97% of the photodissociated ligands are rebound within a few nanoseconds in the 3F mutant case. Our results demonstrate the specific involvement of water molecules in shaping the energetic barriers for ligand migration and binding.

H-bonding networks of the distal residues and water molecules in the active site of Thermobifida fusca hemoglobin

The ferric form of truncated hemoglobin II from Thermobifida fusca (Tf-trHb) and its triple mutant WG8F-YB10F-YCD1F at neutral and alkaline pH, and in the presence of CN(-) have been characterized by resonance Raman spectroscopy, electron paramagnetic resonance spectroscopy, and molecular dynamics simulations. Tf-trHb contains three polar residues in the distal site, namely TrpG8, TyrCD1 and TyrB10. Whereas TrpG8 can act as a potential hydrogen-bond donor, the tyrosines can act as donors or acceptors. Ligand binding in heme-containing proteins is determined by a number of factors, including the nature and conformation of the distal residues and their capability to stabilize the heme-bound ligand via hydrogen-bonding and electrostatic interactions. Since both the RR Fe-OH(-) and Fe-CN(-) frequencies are very sensitive to the distal environment, detailed information on structural variations has been obtained. The hydroxyl ligand binds only the WT protein giving rise to two different conformers. In form 1 the anion is stabilized by H-bonds with TrpG8, TyrCD1 and a water molecule, in turn H-bonded to TyrB10. In form 2, H-bonding with TyrCD1 is mediated by a water molecule. Unlike the OH(-) ligand, CN(-) binds both WT and the triple mutant giving rise to two forms with similar spectroscopic characteristics. The overall results clearly indicate that H-bonding interactions both with distal residues and water molecules are important structural determinants in the active site of Tf-trHb. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Chronic pain in a victim of attempted homicide

This Medical Psychology Forum focuses on the problem of chronic pain in a patient with a neurologic deficit incurred when he became the victim of an attempted homicide. Murderous anger, wishes for revenge and intolerance of helplessness and passivity proved closely intertwined with the patient's awareness of pain and his capacity to bear it. The use of a dorsal column stimulator to “mask” the pain and the role of personality factors in his adjustment to disability are discussed.