ProDeGe: a computational protocol for fully automated decontamination of genomes

Single amplified genomes and genomes assembled from metagenomes have enabled the exploration of uncultured microorganisms at an unprecedented scale. However, both these types of products are plagued by contamination. Since these genomes are now being generated in a high-throughput manner and sequences from them are propagating into public databases to drive novel scientific discoveries, rigorous quality controls and decontamination protocols are urgently needed. Here, we present ProDeGe (Protocol for fully automated Decontamination of Genomes), the first computational protocol for fully automated decontamination of draft genomes. ProDeGe classifies sequences into two classes—clean and contaminant—using a combination of homology and feature-based methodologies. On average, 84% of sequence from the non-target organism is removed from the data set (specificity) and 84% of the sequence from the target organism is retained (sensitivity). The procedure operates successfully at a rate of ~0.30 CPU core hours per megabase of sequence and can be applied to any type of genome sequence.Recent technological advancements have enabled the large-scale sampling of genomes from uncultured microbial taxa, through the high-throughput sequencing of single amplified genomes (SAGs; Rinke et al., 2013; Swan et al., 2013) and assembly and binning of genomes from metagenomes (GMGs; Cuvelier et al., 2010; Sharon and Banfield, 2013). The importance of these products in assessing community structure and function has been established beyond doubt (Kalisky and Quake, 2011). Multiple Displacement Amplification (MDA) and sequencing of single cells has been immensely successful in capturing rare and novel phyla, generating valuable references for phylogenetic anchoring. However, efforts to conduct MDA and sequencing in a high-throughput manner have been heavily impaired by contamination from DNA introduced by the environmental sample, as well as introduced during the MDA or sequencing process (Woyke et al., 2011; Engel et al., 2014; Field et al., 2014). Similarly, metagenome binning and assembly often carries various errors and artifacts depending on the methods used (Nielsen et al., 2014). Even cultured isolate genomes have been shown to lack immunity to contamination with other species (Parks et al., 2014; Mukherjee et al., 2015). As sequencing of these genome product types rapidly increases, contaminant sequences are finding their way into public databases as reference sequences. It is therefore extremely important to define standardized and automated protocols for quality control and decontamination, which would go a long way towards establishing quality standards for all microbial genome product types.Current procedures for decontamination and quality control of genome sequences in single cells and metagenome bins are heavily manual and can consume hours/megabase when performed by expert biologists. Supervised decontamination typically involves homology-based inspection of ribosomal RNA sequences and protein coding genes, as well as visual analysis of k-mer frequency plots and guanine–cytosine content (Clingenpeel, 2015). Manual decontamination is also possible through the software SmashCell (Harrington et al., 2010), which contains a tool for visual identification of contaminants from a self-organizing map and corresponding U-matrix. Another existing software tool, DeconSeq (Schmieder and Edwards, 2011), automatically removes contaminant sequences, however, the contaminant databases are required input. The former lacks automation, whereas the latter requires prior knowledge of contaminants, rendering both applications impractical for high-throughput decontamination.Here, we introduce ProDeGe, the first fully automated computational protocol for decontamination of genomes. ProDeGe uses a combination of homology-based and sequence composition-based approaches to separate contaminant sequences from the target genome draft. It has been pre-calibrated to discard at least 84% of the contaminant sequence, which results in retention of a median 84% of the target sequence. The standalone software is freely available at http://prodege.jgi-psf.org//downloads/src and can be run on any system that has Perl, R (R Core Team, 2014), Prodigal (Hyatt et al., 2010) and NCBI Blast (Camacho et al., 2009) installed. A graphical viewer allowing further exploration of data sets and exporting of contigs accompanies the web application for ProDeGe at http://prodege.jgi-psf.org, which is open to the wider scientific community as a decontamination service (Supplementary Figure S1).The assembly and corresponding NCBI taxonomy of the data set to be decontaminated are required inputs to ProDeGe (Figure 1a). Contigs are annotated with genes following which, eukaryotic contamination is removed based on homology of genes at the nucleotide level using the eukaryotic subset of NCBI''s Nucleotide database as the reference. For detecting prokaryotic contamination, a curated database of reference contigs from the set of high-quality genomes within the Integrated Microbial Genomes (IMG; Markowitz et al., 2014) system is used as the reference. This ensures that errors in public reference databases due to poor quality of sequencing, assembly and annotation do not negatively impact the decontamination process. Contigs determined as belonging to the target organism based on nucleotide level homology to sequences in the above database are defined as ‘Clean'', whereas those aligned to other organisms are defined as ‘Contaminant''. Contigs whose origin cannot be determined based on alignment are classified as ‘Undecided''. Classified clean and contaminated contigs are used to calibrate the separation in the subsequent 5-mer based binning module, which classifies undecided contigs as ‘Clean'' or ‘Contaminant'' using principal components analysis (PCA) of 5-mer frequencies. This parameter can also be specified by the user. When data sets do not have taxonomy deeper than phylum level, or a single confident taxonomic bin cannot be detected using sequence alignment, solely 9-mer based binning is used due to more accurate overall classification. In the absence of a user-defined cutoff, a pre-calibrated cutoff for 80% or more specificity separates the clean contigs from contaminated sequences in the resulting PCA of the 9-mer frequency matrix. Details on ProDeGe''s custom database, evaluation of the performance of the system and exploration of the parameter space to calibrate ProDeGe for a high accurate classification rate are provided in the Supplementary Material.Open in a separate windowFigure 1(a) Schematic overview of the ProDeGe engine. (b) Features of data sets used to validate ProDeGe: SAGs from the Arabidopsis endophyte sequencing project, MDM project, public data sets found in IMG but not sequenced at the JGI, as well as genomes from metagenomes. All the data and results can be found in Supplementary Table S3.The performance of ProDeGe was evaluated using 182 manually screened SAGs (Figure 1b,Supplementary Table S1) from two studies whose data sets are publicly available within the IMG system: genomes of 107 SAGs from an Arabidopsis endophyte sequencing project and 75 SAGs from the Microbial Dark Matter (MDM) project* (only 75/201 SAGs from the MDM project had 1:1 mapping between contigs in the unscreened and the manually screened versions, hence these were used; Rinke et al., 2013). Manual curation of these SAGs demonstrated that the use of ProDeGe prevented 5311 potentially contaminated contigs in these data sets from entering public databases. Figure 2a demonstrates the sensitivity vs specificity plot of ProDeGe results for the above data sets. Most of the data points in Figure 2a cluster in the top right of the box reflecting a median retention of 89% of the clean sequence (sensitivity) and a median rejection of 100% of the sequence of contaminant origin (specificity). In addition, on average, 84% of the bases of a data set are accurately classified. ProDeGe performs best when the target organism has sequenced homologs at the class level or deeper in its high-quality prokaryotic nucleotide reference database. If the target organism''s taxonomy is unknown or not deeper than domain level, or there are few contigs with taxonomic assignments, a target bin cannot be assessed and thus ProDeGe removes contaminant contigs using sequence composition only. The few samples in Figure 2a that demonstrate a higher rate of false positives (lower specificity) and/or reduced sensitivity typically occur when the data set contains few contaminant contigs or ProDeGe incorrectly assumes that the largest bin is the target bin. Some data sets contain a higher proportion of contamination than target sequence and ProDeGe''s performance can suffer under this condition. However, under all other conditions, ProDeGe demonstrates high speed, specificity and sensitivity (Figure 2). In addition, ProDeGe demonstrates better performance in overall classification when nucleotides are considered than when contigs are considered, illustrating that longer contigs are more accurately classified (Supplementary Table S1).Open in a separate windowFigure 2ProDeGe accuracy and performance scatterplots of 182 manually curated single amplified genomes (SAGs), where each symbol represents one SAG data set. (a) Accuracy shown by sensitivity (proportion of bases confirmed ‘Clean'') vs specificity (proportion of bases confirmed ‘Contaminant'') from the Endophyte and Microbial Dark Matter (MDM) data sets. Symbol size reflects input data set size in megabases. Most points cluster in the top right of the plot, showing ProDeGe''s high accuracy. Median and average overall results are shown in Supplementary Table S1. (b) ProDeGe completion time in central processing unit (CPU) core hours for the 182 SAGs. ProDeGe operates successfully at an average rate of 0.30 CPU core hours per megabase of sequence. Principal components analysis (PCA) of a 9-mer frequency matrix costs more computationally than PCA of a 5-mer frequency matrix used with blast-binning. The lack of known taxonomy for the MDM data sets prevents blast-binning, thus showing longer finishing times than the endophyte data sets, which have known taxonomy for use in blast-binning.All SAGs used in the evaluation of ProDeGe were assembled using SPAdes (Bankevich et al., 2012). In-house testing has shown that reads assembled with SPAdes from different strains or even slightly divergent species of the same genera may be combined into the same contig (Personal communications, KT and Robert Bowers). Ideally, the DNA in a well that gets sequenced belongs to a single cell. In the best case, contaminant sequences need to be at least from a different species to be recognized as such by the homology-based screening stage. In the absence of closely related sequenced organisms, contaminant sequences need to be at least from a different genus to be recognized as such by the composition-based screening stage (Supplementary Material). Thus, there is little risk of ProDeGe separating sequences from clonal populations or strains. We have found species- and genus-level contamination in MDA samples to be rare.To evaluate the quality of publicly available uncultured genomes, ProDeGe was used to screen 185 SAGs and 14 GMGs (Figure 1b). Compared with CheckM (Parks et al., 2014), a tool which calculates an estimate of genome sequence contamination using marker genes, ProDeGe generally marks a higher proportion of sequence as ‘Contaminant'' (Supplementary Table S2). This is because ProDeGe has been calibrated to perform at high specificity levels. The command line version of ProDeGe allows users to conduct their own calibration and specify a user-defined distance cutoff. Further, CheckM only outputs the proportion of contamination, but ProDeGe actually labels each contig as ‘Clean'' or ‘Contaminant'' during the process of automated removal.The web application for ProDeGe allows users to export clean and contaminant contigs, examine contig gene calls with their corresponding taxonomies, and discover contig clusters in the first three components of their k-dimensional space. Non-linear approaches for dimensionality reduction of k-mer vectors are gaining popularity (van der Maaten and Hinton, 2008), but we observed no systematic advantage of using t-Distributed Stochastic Neighbor Embedding over PCA (Supplementary Figure S2).ProDeGe is the first step towards establishing a standard for quality control of genomes from both cultured and uncultured microorganisms. It is valuable for preventing the dissemination of contaminated sequence data into public databases, avoiding resulting misleading analyses. The fully automated nature of the pipeline relieves scientists of hours of manual screening, producing reliably clean data sets and enabling the high-throughput screening of data sets for the first time. ProDeGe, therefore, represents a critical component in our toolkit during an era of next-generation DNA sequencing and cultivation-independent microbial genomics.     

Increased Frequency of T Follicular Helper Cells and Elevated Interleukin-27 Plasma Levels in Patients with Pemphigus

Pemphigus is an autoimmune disease in which IgG auto-antibodies (auto-ab) against the desmosomal cadherins desmoglein (Dsg) 3 and Dsg1 cause loss of epidermal keratinocyte adhesion. Aim of this study was to investigate cytokines derived from antigen-presenting cells (APC) and their relation to CD4⁺ T cell subpopulations and to the auto-ab response in pemphigus. In this regard, patients with pemphigus were compared to patients with myasthenia gravis (MG), an unrelated auto-ab–mediated autoimmune disease, and healthy controls. In pemphigus and MG, the plasma concentrations of the APC-derived immunomodulatory cytokine IL-27 were highly increased. Strikingly, IL-27 strongly correlated with Dsg-specific IgG auto-ab titers. T helper (Th) 17 cells were augmented in both pemphigus and MG patients while T follicular helper (Tfh) cells, which are essential in providing B cell help, were increased only in pemphigus along with increasing plasma concentrations of IL-21, a cytokine produced by Th17 and Tfh cells. Moreover, we could detect Dsg3-specific autoreactive T cells producing IL-21 upon ex vivo stimulation with Dsg3. These findings suggest that IL-27 and IL-21-producing T cells, are involved in the pathogenesis of pemphigus. The further characterization of IL-21-producing T cells and of the role of IL-27 will lead to a more defined understanding of the auto-ab response in pemphigus.     

The Role of Nature and Nurture for Individual Differences in Primary Emotional Systems: Evidence from a Twin Study

The present study investigated for the first time the relative importance of genetics and environment on individual differences in primary emotionality as measured with the Affective Neuroscience Personality Scales (ANPS) by means of a twin-sibling study design. In N = 795 participants (n = 303 monozygotic twins, n = 172 dizygotic twins and n = 267 non-twin full siblings), moderate to strong influences of genetics on individual differences in these emotional systems are observed. Lowest heritability estimates are presented for the SEEKING system (33%) and highest for the PLAY system (69%). Further, multivariate genetic modeling was applied to the data showing that associations among the six ANPS scales were influences by both, a genetic as well as an environmental overlap between them. In sum, the study underlines the usefulness of the ANPS for biologically oriented personality psychology research.     

Attentional Bias towards Positive Emotion Predicts Stress Resilience

There is extensive evidence for an association between an attentional bias towards emotionally negative stimuli and vulnerability to stress-related psychopathology. Less is known about whether selective attention towards emotionally positive stimuli relates to mental health and stress resilience. The current study used a modified Dot Probe task to investigate if individual differences in attentional biases towards either happy or angry emotional stimuli, or an interaction between these biases, are related to self-reported trait stress resilience. In a nonclinical sample (N = 43), we indexed attentional biases as individual differences in reaction time for stimuli preceded by either happy or angry (compared to neutral) face stimuli. Participants with greater attentional bias towards happy faces (but not angry faces) reported higher trait resilience. However, an attentional bias towards angry stimuli moderated this effect: The attentional bias towards happy faces was only predictive for resilience in those individuals who also endorsed an attentional bias towards angry stimuli. An attentional bias towards positive emotional stimuli may thus be a protective factor contributing to stress resilience, specifically in those individuals who also endorse an attentional bias towards negative emotional stimuli. Our findings therefore suggest a novel target for prevention and treatment interventions addressing stress-related psychopathology.     

New Acylated Flavonol Glycosides and a Phenolic Profile of Pritzelago alpina,a Forgotten Edible Alpine Plant

Thirteen acylated flavonoid glycosides, 1 – 13 , including eleven new congeners, 3 – 13 , were isolated from the aerial parts of Pritzelago alpina (Brassicaceae) by a combination of column chromatography on Sephadex LH‐20, and preparative and semi‐preparative HPLC. The structures were established by extensive NMR and MS experiments in combination with acid hydrolysis and sugar analysis by GC/MS. The new compounds were shown to be kaempferol and quercetin glycosides acylated for most of them by a branched short chain fatty acid or a hydroxycinnamic acid residue on the sugar portion. As shown by a HPLC‐DAD analysis of a MeOH extract, these compounds are the main phenolic constituents in the aerial parts of the plant.     

Xylene‐Bridged Phosphaviologen Oligomers and Polymers as High‐Performance Electrode‐Modifiers for Li‐Ion Batteries

Lithium‐ion batteries are one of the most common forms of energy storage devices used in society today. Due to the inherent limitations of conventional Li‐ion batteries, organic materials have surfaced as potentially suitable electrode alternatives with improved performance and sustainability. Viologens and phosphaviologens in particular, are strong electron‐accepting materials with excellent kinetic properties, making them suitable candidates for battery applications. In this paper, new polymeric species of the latter moieties are reported that lead to improved electrode stability and device performance. The performance of the phosphaviologen is further enhanced through the utilization of both redox steps, allowing for good performance proof‐of‐concept hybrid organic/Li‐ion batteries. This opens the potential for more sustainable and improved battery performance for use in current energy applications.     

Low‐Overpotential Electrocatalytic Water Splitting with Noble‐Metal‐Free Nanoparticles Supported in a sp3 N‐Rich Flexible COF

Covalent organic frameworks (COFs) are crystalline organic polymers with tunable structures. Here, a COF is prepared using building units with highly flexible tetrahedral sp³ nitrogens. This flexibility gives rise to structural changes which generate mesopores capable of confining very small (<2 nm sized) non‐noble‐metal‐based nanoparticles (NPs). This nanocomposite shows exceptional activity toward the oxygen‐evolution reaction from alkaline water with an overpotential of 258 mV at a current density of 10 mA cm^?2. The overpotential observed in the COF‐nanoparticle system is the best in class, and is close to the current record of ≈200 mV for any noble‐metal‐free electrocatalytic water splitting system—the Fe–Co–Ni metal‐oxide‐film system. Also, it possesses outstanding kinetics (Tafel slope of 38.9 mV dec^?1) for the reaction. The COF is able to stabilize such small‐sized NP in the absence of any capping agent because of the COF–Ni(OH)₂ interactions arising from the N‐rich backbone of the COF. Density‐functional‐theory modeling of the interaction between the hexagonal Ni(OH)₂ nanosheets and the COF shows that in the most favorable configuration the Ni(OH)₂ nanosheets are sandwiched between the sp³ nitrogens of the adjacent COF layers and this can be crucial to maximizing their synergistic interactions.     

Toward a Low‐Cost Artificial Leaf: Driving Carbon‐Based and Bifunctional Catalyst Electrodes with Solution‐Processed Perovskite Photovoltaics

Molecular hydrogen can be generated renewably by water splitting with an “artificial‐leaf device”, which essentially comprises two electrocatalyst electrodes immersed in water and powered by photovoltaics. Ideally, this device should operate efficiently and be fabricated with cost‐efficient means using earth‐abundant materials. Here, a lightweight electrocatalyst electrode, comprising large surface‐area NiCo₂O₄ nanorods that are firmly anchored onto a carbon–paper current collector via a dense network of nitrogen‐doped carbon nanotubes is presented. This electrocatalyst electrode is bifunctional in that it can efficiently operate as both anode and cathode in the same alkaline solution, as quantified by a delivered current density of 10 mA cm^?2 at an overpotential of 400 mV for each of the oxygen and hydrogen evolution reactions. By driving two such identical electrodes with a solution‐processed thin‐film perovskite photovoltaic assembly, a wired artificial‐leaf device is obtained that features a Faradaic H₂ evolution efficiency of 100%, and a solar‐to‐hydrogen conversion efficiency of 6.2%. A detailed cost analysis is presented, which implies that the material‐payback time of this device is of the order of 100 days.     

Quantitative assessment of methyl‐esterification and other side reactions in a standard propionylation protocol for detection of histone modifications

Histone modifications play an important role in regulating chromatin stability and gene expression, but to date, investigating them remains challenging. In order to obtain peptides suitable for MS‐based analysis, chemical derivatization of N‐terminus and lysine residues by propionic anhydride is commonly performed. Several side reactions (methyl‐esterification, amidation, solvolysis, overpropionylation, and missed propionylation) during propionylation protocols have been described, yet their relative abundances remain vague. Because methyl‐esterification could interfere with correct interpretation of the modification pattern, it is essential to take measures to avoid it. Here we present in‐depth quantitative analyses of methyl‐esterification and the other side reactions in a standard propionylation protocol containing methanol, and when replacing methanol with isopropanol or acetonitrile. We show that the use of alternative solvents can eliminate methyl‐esterification and that even though other side reactions are not prevented, their contribution can be kept relatively small. We also show that replacing methanol can be of importance also in other proteomics methods, such as mixed cation exchange, using methanol under acidic conditions.     

Poromicromechanics reveals that physiological bone strains induce osteocyte-stimulating lacunar pressure

Mechanical loads which are macroscopically acting onto bony organs, are known to influence the activities of biological cells located in the pore spaces of bone, in particular so the signaling and production processes mediated by osteocytes. The exact mechanisms by which osteocytes are actually able to “feel” the mechanical loading and changes thereof, has been the subject of numerous studies, and, while several hypotheses have been brought forth over time, this topic has remained a matter of debate. Relaxation times reported in a recent experimental study of Gardinier et al. (Bone 46(4):1075–1081, 2010) strongly suggest that the lacunar pores are likely to experience, during typical physiological load cycles, not only fluid transport, but also undrained conditions. The latter entail the buildup of lacunar pore pressures, which we here quantify by means of a thorough multiscale modeling approach. In particular, the proposed model is based on classical poroelasticity theory, and able to account for multiple pore spaces. First, the model reveals distinct nonlinear dependencies of the resulting lacunar (and vascular) pore pressures on the underlying bone composition, highlighting the importance of a rigorous multiscale approach for appropriate computation of the aforementioned pore pressures. Then, the derived equations are evaluated for macroscopic (uniaxial as well as hydrostatic) mechanical loading of physiological magnitude. The resulting model-predicted pore pressures agree very well with the pressures that have been revealed, by means of in vitro studies, to be of adequate magnitude for modulating the responses of biological cells, including osteocytes. This underlines that osteocytes may respond to many types of loading stimuli at the same time, in particular so to fluid flow and hydrostatic pressure.