The InterPro protein families database: the classification resource after 15 years

Morphologically Different Pectobacterium brasiliense Bacteriophages PP99 and PP101: Deacetylation of O-Polysaccharide by the Tail Spike Protein of Phage PP99 Accompanies the Infection

Soft rot caused by numerous species of Pectobacterium and Dickeya is a serious threat to the world production of potatoes. The application of bacteriophages to combat bacterial infections in medicine, agriculture, and the food industry requires the selection of comprehensively studied lytic phages and the knowledge of their infection mechanism for more rational composition of therapeutic cocktails. We present the study of two bacteriophages, infective for the Pectobacterium brasiliense strain F152. Podoviridae PP99 is a representative of the genus Zindervirus, and Myoviridae PP101 belongs to the still unclassified genomic group. The structure of O-polysaccharide of F152 was established by sugar analysis and 1D and 2D NMR spectroscopy:

→ 4)-α-D-Manp6Ac-(1→ 2)-α-D-Manp-(1→ 3)-β-D-Galp-(1→ 3↑1α-l-6dTalpAc0−2

The recombinant tail spike protein of phage PP99, gp55, was shown to deacetylate the side chain talose residue of bacterial O-polysaccharide, thus providing the selective attachment of the phage to the cell surface. Both phages demonstrate lytic behavior, thus being prospective for therapeutic purposes.

The EngCP endo α-N-acetylgalactosaminidase is a virulence factor involved in Clostridium perfringens gas gangrene infections

Clostridium perfringens is the causative agent of human clostridial myonecrosis; the major toxins involved in this disease are α-toxin and perfringolysin O. The RevSR two-component regulatory system has been shown to be involved in regulating virulence in a mouse myonecrosis model. Previous microarray and RNAseq analysis of a revR mutant implied that factors other than the major toxins may play a role in virulence. The RNAseq data showed that the expression of the gene encoding the EngCP endo α-N-acetylgalactosaminidase (CPE0693) was significantly down-regulated in a revR mutant. Enzymes from this family have been identified in several Gram-positive pathogens and have been postulated to contribute to their virulence. In this study, we constructed an engCP mutant of C. perfringens and showed that it was significantly less virulent than its wild-type parent strain. Virulence was restored by complementation in trans with the wild-type engCP gene. We also demonstrated that purified EngCP was able to hydrolyse α-dystroglycan derived from C2C12 mouse myotubes. However, EngCP had little effect on membrane permeability in mice, suggesting that EngCP may play a role other than the disruption of the structural integrity of myofibres. Glycan array analysis indicated that EngCP could recognise structures containing the monosaccharide N-acetlygalactosamine at 4C, but could recognise structures terminating in galactose, glucose and N-acetylglucosamine under conditions where EngCP was enzymatically active. In conclusion, we have obtained evidence that EngCP is required for virulence in C. perfringens and, although classical exotoxins are important for disease, we have now shown that an O-glycosidase also plays an important role in the disease process.

SpinachBase: a central portal for spinach genomics

Spinach (Spinacia oleracea L.) is a nutritious vegetable enriched with many essential minerals and vitamins. A reference spinach genome has been recently released, and additional spinach genomic resources are being rapidly developed. Therefore, there is an urgent need of a central database to store, query, analyze and integrate various resources of spinach genomic data. To this end, we developed SpinachBase (http://spinachbase.org), which provides centralized public accesses to genomic data as well as analytical tools to assist research and breeding in spinach. The database currently stores the spinach reference genome sequence, and sequences and comprehensive functional annotations of protein-coding genes predicted from the genome. The database also contains gene expression profiles derived from RNA-Seq experiments as well as highly co-expressed genes and genetic variants called from transcriptome sequences of 120 cultivated and wild Spinacia accessions. Biochemical pathways have been predicted from spinach protein-coding genes and are available through a pathway database (SpinachCyc) within SpinachBase. SpinachBase provides a suite of analysis and visualization tools including a genome browser, sequence similarity searches with BLAST, functional enrichment and functional classification analyses and functions to query and retrieve gene sequences and annotations.

High-resolution structure of intramolecularly proteolyzed human mucin-1 SEA domain

SEA domains are ubiquitous in large proteins associated with highly glycosylated environments. Certain SEA domains undergo intramolecular proteolysis involving a nucleophilic attack of a serine hydroxyl group on the preceding glycine carbonyl. The mucin-1 (MUC1) SEA domain has been extensively investigated as a model of intramolecular proteolysis. Since neither a general base, a general acid, nor an oxyanion hole could be identified in MUC1 SEA, it has been suggested that proteolysis is accelerated by a non-planarity of the scissile peptide bond imposed by protein folding. A reactant distorted peptide bond has been also invoked to explain the autoproteolysis of several unrelated proteins. However, the only evidence of peptide distortion in MUC1 SEA stems from molecular dynamic simulations of the reactant modeled upon a single NMR structure of the cleaved product. We report the first high-resolution X-ray structure of cleaved MUC1 SEA. Structural comparison with uncleaved SEA domains suggests that the number of residues evolutionarily inserted in the cleaved loop of MUC1 SEA precludes the formation of a properly hydrogen-bonded beta turn. By sequence analysis, we show that this conformational frustration is shared by all known cleaved SEA domains. In addition, alternative conformations of the uncleaved precursor could be modeled in which the scissile peptide bond is planar. The implications of these structures for autoproteolysis are discussed in the light of the previous research on autoproteolysis.

Whole Genome Sequencing and Assembly of the Asian Honey Bee Apis dorsata

The Asian honey bee (Apis dorsata) is distinct from its more widely distributed cousin Apis mellifera by a few key characteristics. Most prominently, A. dorsata, nest in the open by forming a colony clustered around the honeycomb, whereas A. mellifera nest in concealed cavities. Additionally, the worker and reproductive castes are all of the same size in A. dorsata. In order to investigate these differences, we performed whole genome sequencing of A. dorsata using a hybrid Oxford Nanopore and Illumina approach. The 223 Mb genome has an N50 of 35 kb with the largest scaffold of 302 kb. We have found that there are many genes in the dorsata genome that are distinct from other hymenoptera and also large amounts of transposable elements, and we suggest some candidate genes for A. dorsata's exceptional level of defensive aggression.

Learning important features from multi-view data to predict drug side effects

The problem of drug side effects is one of the most crucial issues in pharmacological development. As there are many limitations in current experimental and clinical methods for detecting side effects, a lot of computational algorithms have been developed to predict side effects with different types of drug information. However, there is still a lack of methods which could integrate heterogeneous data to predict side effects and select important features at the same time. Here, we propose a novel computational framework based on multi-view and multi-label learning for side effect prediction. Four different types of drug features are collected and graph model is constructed from each feature profile. After that, all the single view graphs are combined to regularize the linear regression functions which describe the relationships between drug features and side effect labels. L1 penalties are imposed on the regression coefficient matrices in order to select features relevant to side effects. Additionally, the correlations between side effect labels are also incorporated into the model by graph Laplacian regularization. The experimental results show that the proposed method could not only provide more accurate prediction for side effects but also select drug features related to side effects from heterogeneous data. Some case studies are also supplied to illustrate the utility of our method for prediction of drug side effects.

The developmental phosphoproteome of Haemonchus contortus

Protein phosphorylation plays essential roles in many cellular processes. Despite recent progress in the genomics, transcriptomics and proteomics of socioeconomically important parasitic nematodes, there is scant phosphoproteomic data to underpin molecular biological discovery. Here, using the phosphopeptide enrichment-based LC-MS/MS and data-independent acquisition (DIA) quantitation, we characterised the first developmental phosphoproteome of the parasitic nematode Haemonchus contortus - one of the most pathogenic parasites of ruminant livestock. Totally, 1804 phosphorylated proteins with 4406 phosphorylation sites ('phosphosites') from different developmental stages/sexes were identified. Bioinformatic analyses of quantified 'phosphosites' exhibited distinctive stage- and sex-specific patterns during development, and identified a subset of phosphoproteins proposed to play crucial roles in processes such as spindle positioning, signal transduction and kinase activity. A sequence-based comparison of the phosphoproteome of H. contortus with those of two free-living nematode species (Caenorhabditis elegans and Pristionchus pacificus) suggested a limited number of common protein phosphorylation events among these species. Our findings infer active roles for protein phosphorylation in the adaptation of a parasitic nematode to a constantly changing external environment. The phosphoproteomic data set for H. contortus provides a basis to better understand phosphorylation and associated biological processes (e.g., regulation of signal transduction), and might enable the discovery of novel anthelmintic targets. SIGNIFICANCE: Here, we report the first phosphoproteome for a socioeconomically parasitic nematode (Haemonchus contortus). This phosphoproteome exhibits distinctive patterns during development, suggesting active roles of post-translational modification in the parasite's adaptation to changing environments within and outside of the host animal. This work sheds a light on the developmental phosphorylation in a parasitic nematode, and could enable the discovery of novel interventions against major pathogens.

Systematic analysis of linker histone PTM hotspots reveals phosphorylation sites that modulate homologous recombination and DSB repair

Double strand-breaks (DSBs) of genomic DNA caused by ionizing radiation or mutagenic chemicals are a common source of mutation, recombination, chromosomal aberration, and cell death. Linker histones are DNA packaging proteins with established roles in chromatin compaction, gene transcription, and in homologous recombination (HR)-mediated DNA repair. Using a machine-learning model for functional prioritization of eukaryotic post-translational modifications (PTMs) in combination with genetic and biochemical experiments with the yeast linker histone, Hho1, we discovered that site-specific phosphorylation sites regulate HR and HR-mediated DSB repair. Five total sites were investigated (T10, S65, S141, S173, and S174), ranging from high to low function potential as determined by the model. Of these, we confirmed S173/174 are phosphorylated in yeast by mass spectrometry and found no evidence of phosphorylation at the other sites. Phospho-nullifying mutations at these two sites results in a significant decrease in HR-mediated DSB repair templated either with oligonucleotides or a homologous chromosome, while phospho-mimicing mutations have no effect. S65, corresponding to a mammalian phosphosite that is conserved in yeast, exhibited similar effects. None of the mutations affected base- or nucleotide-excision repair, nor did they disrupt non-homologous end joining or RNA-mediated repair of DSBs when sequence heterology between the break and repair template strands was low. More extensive analysis of the S174 phospho-null mutant revealed that its repression of HR and DSB repair is proportional to the degree of sequence heterology between DSB ends and the HR repair template. Taken together, these data demonstrate the utility of machine learning for the discovery of functional PTM hotspots, reveal linker histone phosphorylation sites necessary for HR and HR-mediated DSB repair, and provide insight into the context-dependent control of DNA integrity by the yeast linker histone Hho1.

TGFB1-Mediated Gliosis in Multiple Sclerosis Spinal Cords Is Favored by the Regionalized Expression of HOXA5 and the Age-Dependent Decline in Androgen Receptor Ligands

In multiple sclerosis (MS) patients with a progressive form of the disease, spinal cord (SC) functions slowly deteriorate beyond age 40. We previously showed that in the SC of these patients, large areas of incomplete demyelination extend distance away from plaque borders and are characterized by a unique progliotic TGFB1 (Transforming Growth Factor Beta 1) genomic signature. Here, we attempted to determine whether region- and age-specific physiological parameters could promote the progression of SC periplaques in MS patients beyond age 40. An analysis of transcriptomics databases showed that, under physiological conditions, a set of 10 homeobox (HOX) genes are highly significantly overexpressed in the human SC as compared to distinct brain regions. Among these HOX genes, a survey of the human proteome showed that only HOXA5 encodes a protein which interacts with a member of the TGF-beta signaling pathway, namely SMAD1 (SMAD family member 1). Moreover, HOXA5 was previously found to promote the TGF-beta pathway. Interestingly, SMAD1 is also a protein partner of the androgen receptor (AR) and an unsupervised analysis of gene ontology terms indicates that the AR pathway antagonizes the TGF-beta/SMAD pathway. Retrieval of promoter analysis data further confirmed that AR negatively regulates the transcription of several members of the TGF-beta/SMAD pathway. On this basis, we propose that in progressive MS patients, the physiological SC overexpression of HOXA5 combined with the age-dependent decline in AR ligands may favor the slow progression of TGFB1-mediated gliosis. Potential therapeutic implications are discussed.

Identification of a metallothionein gene in honey bee Apis mellifera and its expression profile in response to Cd, Cu and Pb exposure

Metallothioneins are ubiquitous proteins important in metal homeostasis and detoxification. However, they have not previously been identified in honey bees or other Hymenoptera, where metallothioneins could be of ecophysiological and ecotoxicological significance. Better understanding of the molecular responses to stress induced by toxic metals could contribute to honey bee conservation. In addition, honey bee metallothionein could represent a biomarker for monitoring environmental quality. Here we identify and characterize a metallothionein gene in Apis mellifera (AmMT). AmMT is 1,680 bp long and encodes a 48 amino acids protein with 15 cysteines and no aromatic residues. A metal response element upstream of the start codon, coupled with numerous cis-regulatory elements indicate the functional context of AmMT. Molecular modelling predicts several transition metal binding sites, and comparative phylogenetic analysis revealed five putative metallothionein proteins in three other hymenoptera species. AmMT was characterized by cloning the full-length coding sequence of the putative metallothionein. Recombinant AmMT was found to increase metal tolerance upon overexpression in Escherichia coli supplemented with Cd, Cu or Pb. Finally, in laboratory tests on honey bees, gene expression profiles showed a dose-dependant relationship between Cd, Cu and Pb concentrations present in food and AmMT expression, while field experiments showed induction of AmMT in bees from an industrial site compared to those from an urban area. These studies suggest that AmMT has metal binding properties in agreement with a possible role in metal homeostasis. Further functional and structural characterization of metallothionein in honey bees and other Hymenoptera are necessary.

Time-coursed transcriptome analysis identifies key expressional regulation in growth cessation and dormancy induced by short days in Paulownia

Maintaining the viability of the apical shoot is critical for continued vertical growth in plants. Terminal shoot of tree species Paulownia cannot regrow in subsequent years. The short day (SD) treatment leads to apical growth cessation and dormancy. To understand the molecular basis of this, we further conducted global RNA-Seq based transcriptomic analysis in apical shoots to check regulation of gene expression. We obtained ~219 million paired-end 125-bp Illumina reads from five time-courses and de novo assembled them to yield 49,054 unigenes. Compared with the untreated control, we identified 1540 differentially expressed genes (DEGs) which were found to involve in 116 metabolic pathways. Expression of 87% of DEGs exhibited switch-on or switch-off pattern, indicating key roles in growth cessation. Most DEGs were enriched in the biological process of gene ontology categories and at later treatment stages. The pathways of auxin and circadian network were most affected and the expression of associated DEGs was characterised. During SD induction, auxin genes IAA, ARF and SAURs were down-regulated and circadian genes including PIF3 and PRR5 were up-regulated. PEPC in photosynthesis was constitutively upregulated, suggesting a still high CO₂ concentrating activity; however, the converting CO₂ to G3P in the Calvin cycle is low, supported by reduced expression of GAPDH encoding the catalysing enzyme for this step. This indicates a de-coupling point in the carbon fixation. The results help elucidate the molecular mechanisms for SD inducing dormancy and cessation in apical shoots.

HeteroDualNet: A Dual Convolutional Neural Network With Heterogeneous Layers for Drug-Disease Association Prediction via Chou's Five-Step Rule

Identifying new treatments for existing drugs can help reduce drug development costs and explore novel indications of drugs. The prediction of associations between drugs and diseases is challenging because their similarities and relations are complicated and non-linear. We propose a HeteroDualNet model to address this issue. Firstly, three types of matrices are extracted to represent intra-drug similarities, intra-disease similarity and drug-disease associations. The intra-drug similarities consider three drug features and a newly introduced drug-related disease correlation. Secondly, an embedding mechanism is proposed to integrate these matrices in a heterogenous drug-disease association layer (hetero-layer). Further, a neighbouring heterogeneous layer (hetero-layer-N) is constructed to incorporate the biological premise that similar drugs can often treat related diseases. Finally, a dual convolutional neural network is built with hetero-layer and hetero-layer-N as two branches to learn from characteristics of drug-disease and the relations of their neighbours simultaneously. HeteroDualNet outperformed the other four methods in comparison over a public dataset of 763 drugs and 681 diseases in terms of Areas Under the Curves of Receiver Operating Characteristics and Precision-Recall, and recall rate at top k. Case study of five drugs further proved the capacity of HeteroDualNet in finding reliable disease candidates of drugs as validated by database records or literature. Our findings show that the embedded heterogenous layers of original and neighbouring drug-disease representations in a dual neural network improved the association prediction performance.

TreeGrafter: phylogenetic tree-based annotation of proteins with Gene Ontology terms and other annotations

Summary

TreeGrafter is a new software tool for annotating protein sequences using pre-annotated phylogenetic trees. Currently, the tool provides annotations to Gene Ontology (GO) terms, and PANTHER family and subfamily. The approach is generalizable to any annotations that have been made to internal nodes of a reference phylogenetic tree. TreeGrafter takes each input query protein sequence, finds the best matching homologous family in a library of pre-calculated, pre-annotated gene trees, and then grafts it to the best location in the tree. It then annotates the sequence by propagating annotations from ancestral nodes in the reference tree. We show that TreeGrafter outperforms subfamily HMM scoring for correctly assigning subfamily membership, and that it produces highly specific annotations of GO terms based on annotated reference phylogenetic trees. This method will be further integrated into InterProScan, enabling an even broader user community.

Availability and implementation

TreeGrafter is freely available on the web at https://github.com/pantherdb/TreeGrafter, including as a Docker image.

Supplementary information

Supplementary data are available at Bioinformatics online.

Modeling post-translational modifications and cancer-associated mutations that impact the heterochromatin protein 1α-importin α heterodimers

Heterochromatin protein 1α (HP1α) is a protein that mediates cancer-associated processes in the cell nucleus. Proteomic experiments, reported here, demonstrate that HP1α complexes with importin α (IMPα), a protein necessary for its nuclear transport. This data is congruent with Simple Linear Motif (SLiM) analyses that identify an IMPα-binding motif within the linker that joins the two globular domains of this protein. Using molecular modeling and dynamics simulations, we develop a model of the IMPα-HP1α complex and investigate the impact of phosphorylation and genomic variants on their interaction. We demonstrate that phosphorylation of the HP1α linker likely regulates its association with IMPα, which has implications for HP1α access to the nucleus, where it functions. Cancer-associated genomic variants do not abolish the interaction of HP1α but instead lead to rearrangements where the variant proteins maintain interaction with IMPα, but with less specificity. Combined, this new mechanistic insight bears biochemical, cell biological, and biomedical relevance.

Armadillo-like helical domain containing-4 is dynamically expressed in both the first and second heart fields

Armadillo repeat and Armadillo-like helical domain containing proteins form a large family with diverse and fundamental functions in many eukaryotes. Herein we investigated the spatiotemporal expression pattern of Armadillo-like helical domain containing 4 (or Armh4) as an uncharacterized protein coding mouse gene, within the mouse embryo during the initial stages of heart morphogenesis. We found Armh4 is initially expressed in both first heart field as well as the second heart field progenitors and subsequently within predominantly their cardiomyocyte derivatives. Armh4 expression is initially cardiac-restricted in the developing embryo and is expressed in second heart field subpharyngeal mesoderm prior to cardiomyocyte differentiation, but Armh4 diminishes as the embryonic heart matures into the fetal heart. Armh4 is subsequently expressed in craniofacial structures and neural crest-derived dorsal root and trigeminal ganglia. Whereas lithium chloride-induced stimulation of Wnt/β-catenin signaling elevated Armh4 expression in both second heart field subpharyngeal mesodermal progenitors and outflow tract, right ventricle and atrial cardiomyocytes, neither a systemic loss of Islet-1 nor an absence of cardiac neural crest cells had any effect upon Armh4 expression. These results confirm that Wnt/β-catenin-responsive Armh4 is a useful specific biomarker of the FHF and SHF cardiomyocyte derivatives only.

Light-Regulation of Tryptophan Synthase by Combining Protein Design and Enzymology

The spatiotemporal control of enzymes by light is of growing importance for industrial biocatalysis. Within this context, the photo-control of allosteric interactions in enzyme complexes, common to practically all metabolic pathways, is particularly relevant. A prominent example of a metabolic complex with a high application potential is tryptophan synthase from Salmonella typhimurium (TS), in which the constituting TrpA and TrpB subunits mutually stimulate each other via a sophisticated allosteric network. To control TS allostery with light, we incorporated the unnatural amino acid o-nitrobenzyl-O-tyrosine (ONBY) at seven strategic positions of TrpA and TrpB. Initial screening experiments showed that ONBY in position 58 of TrpA (aL58ONBY) inhibits TS activity most effectively. Upon UV irradiation, ONBY decages to tyrosine, largely restoring the capacity of TS. Biochemical characterization, extensive steady-state enzyme kinetics, and titration studies uncovered the impact of aL58ONBY on the activities of TrpA and TrpB and identified reaction conditions under which the influence of ONBY decaging on allostery reaches its full potential. By applying those optimal conditions, we succeeded to directly light-activate TS(aL58ONBY) by a factor of ~100. Our findings show that rational protein design with a photo-sensitive unnatural amino acid combined with extensive enzymology is a powerful tool to fine-tune allosteric light-activation of a central metabolic enzyme complex.

Nature's recyclers: anaerobic microbial communities drive crude biomass deconstruction

Microbial communities within anaerobic ecosystems have evolved to degrade and recycle carbon throughout the earth. A number of strains have been isolated from anaerobic microbial communities, which are rich in carbohydrate active enzymes (CAZymes) to liberate fermentable sugars from crude plant biomass (lignocellulose). However, natural anaerobic communities host a wealth of microbial diversity that has yet to be harnessed for biotechnological applications to hydrolyze crude biomass into sugars and value-added products. This review highlights recent advances in 'omics' techniques to sequence anaerobic microbial genomes, decipher microbial membership, and characterize CAZyme diversity in anaerobic microbiomes. With a focus on the herbivore rumen, we further discuss methods to discover new CAZymes, including those found within multi-enzyme fungal cellulosomes. Emerging techniques to characterize the interwoven metabolism and spatial interactions between anaerobes are also reviewed, which will prove critical to developing a predictive understanding of anaerobic communities to guide in microbiome engineering.

Genomic Analyses Provide Insights Into the Evolutionary History and Genetic Diversity of Auricularia Species

Species in the genus Auricularia play important roles for people’s food and nutrition especially Auricularia cornea and A. heimuer. To understand their evolutionary history, genome structure, and population-level genetic variation, we performed a high-quality genome sequencing of Auricularia cornea and the corresponding comparative genomic analysis. The genome size of A. cornea was similar to Auricularia subglabra, but 1.5 times larger than that of A. heimuer. Several factors were responsible for genome size variation including gene numbers, repetitive elements, and gene lengths. Phylogenomic analysis revealed that the estimated divergence time between A. heimuer and other Auricularia is ∼79.1 million years ago (Mya), while the divergence between A. cornea and A. subglabra occurred in ∼54.8 Mya. Population genomic analysis also provided insight into the demographic history of A. cornea and A. heimuer, indicating that their populations fluctuated over time with global climate change during Marine Isotope Stage 5-2. Moreover, despite the highly similar external morphologies of A. cornea and A. heimuer, their genomic properties were remarkably different. The A. cornea genome only shared 14% homologous syntenic blocks with A. heimuer and possessed more genes encoding carbohydrate-active enzymes and secondary metabolite biosynthesis proteins. The cross-taxa transferability rates of simple sequence repeat (SSR) and insertion or deletion (InDel) markers within the genus Auricularia were also lower than that previously observed for species within the same genus. Taken together, these results indicate a high level of genetic differentiation between these two Auricularia species. Consequently, our study provides new insights into the genomic evolution and genetic differentiation of Auricularia species that will facilitate future genetic breeding.

Collection of antimicrobial peptides database and its derivatives: Applications and beyond

Collection of antimicrobial peptides (CAMP), CAMPSign, and ClassAMP are open-access resources that have been developed to enhance research on antimicrobial peptides (AMPs). Comprehensive information on AMPs and machine learning-based predictive models are made available for users through these resources. As of date, CAMP_R3 has 10,247 sequences, 757 structures, and 114 family-specific signatures of AMPs along with associated tools for AMP sequence and structure analysis. CAMPSign uses family-specific sequence conservation, in the form of patterns and hidden Markov models for identification of AMPs. ClassAMP can be used to classify AMPs as antibacterial, antifungal, or antiviral based on sequence information. Here we describe CAMP and its derivatives and illustrate, with a few examples, the contribution of these online resources to the advancement of our current understanding of AMPs.

Reference gene and small RNA data from multiple tissues of Davidia involucrata Baill

Davidia involucrata Baill. is a rare plant endemic to China. Its exclusive evolutionary position and specific floral organs endow it with a high research value. However, a lack of genomic resources has constrained the study of D. involucrata functional genomics. Here, we report D. involucrata transcriptome reads from different floral tissues pooled from six individuals at two developmental stages using Illumina HiSeq technology and the construction of a high-quality reference gene set containing a total of 104,463 unigenes with an N50 of 1,693 bp and 48,529 high-quality coding sequences. The transcriptome data exhibited 89.24% full-length completeness with respect to the benchmarking universal single-copy (BUSCO) dataset and a PLAZA CoreGF weighted score of 98.85%. In total, 65,534 (62.73%) unigenes were functionally annotated, including 58 transcription factor families and 44,327 simple sequence repeats (SSRs). In addition, 96 known and 112 novel miRNAs were identified in the parallel small RNA sequencing of each sample. All these high-quality data could provide a valuable annotated gene set for subsequent studies of D. involucrata.

A draft genome and transcriptome of common milkweed (Asclepias syriaca) as resources for evolutionary, ecological, and molecular studies in milkweeds and Apocynaceae

Milkweeds (Asclepias) are used in wide-ranging studies including floral development, pollination biology, plant-insect interactions and co-evolution, secondary metabolite chemistry, and rapid diversification. We present a transcriptome and draft nuclear genome assembly of the common milkweed, Asclepias syriaca. This reconstruction of the nuclear genome is augmented by linkage group information, adding to existing chloroplast and mitochondrial genomic resources for this member of the Apocynaceae subfamily Asclepiadoideae. The genome was sequenced to 80.4× depth and the draft assembly contains 54,266 scaffolds ≥1 kbp, with N50 = 3,415 bp, representing 37% (156.6 Mbp) of the estimated 420 Mbp genome. A total of 14,474 protein-coding genes were identified based on transcript evidence, closely related proteins, and ab initio models, and 95% of genes were annotated. A large proportion of gene space is represented in the assembly, with 96.7% of Asclepias transcripts, 88.4% of transcripts from the related genus Calotropis, and 90.6% of proteins from Coffea mapping to the assembly. Scaffolds covering 75 Mbp of the Asclepias assembly formed 11 linkage groups. Comparisons of these groups with pseudochromosomes in Coffea found that six chromosomes show consistent stability in gene content, while one may have a long history of fragmentation and rearrangement. The progesterone 5β-reductase gene family, a key component of cardenolide production, is likely reduced in Asclepias relative to other Apocynaceae. The genome and transcriptome of common milkweed provide a rich resource for future studies of the ecology and evolution of a charismatic plant family.

Characterization of Vibrio cholerae's Extracellular Nuclease Xds

The Gram-negative bacterium Vibrio cholerae encodes two nucleases, Dns and Xds, which play a major role during the human pathogen's lifecycle. Dns and Xds control three-dimensional biofilm formation and bacterial detachment from biofilms via degradation of extracellular DNA and thus contribute to the environmental, inter-epidemic persistence of the pathogen. During intestinal colonization the enzymes help evade the innate immune response, and therefore promote survival by mediating escape from neutrophil extracellular traps. Xds has the additional function of degrading extracellular DNA down to nucleotides, which are an important nutrient source for V. cholerae. Thus, Xds is a key enzyme for survival fitness during distinct stages of the V. cholerae lifecycle and could be a potential therapeutic target. This study provides detailed information about the enzymatic properties of Xds using purified protein in combination with a real time nuclease activity assay. The data define an optimal buffer composition for Xds activity as 50 mM Tris/HCl pH 7, 100 mM NaCl, 10 mM MgCl2, and 20 mM CaCl2. Moreover, maximal activity was observed using substrate DNA with low GC content and ambient temperatures of 20-25°C. In silico analysis and homology modeling predicted an exonuclease domain in the C-terminal part of the protein. Biochemical analyses with truncated variants and point mutants of Xds confirm that the C-terminal region is sufficient for nuclease activity. We also find that residues D787 and H837 within the predicted exonuclease domain are key to formation of the catalytic center.

Genome-wide identification and characterization of the Hsp70 gene family in allopolyploid rapeseed (Brassica napus L.) compared with its diploid progenitors

Heat shock protein 70 (Hsp70) plays an essential role in plant growth and development, as well as stress response. Rapeseed (Brassica napus L.) originated from recently interspecific hybridization between Brassica rapa and Brassica oleracea. In this study, a total of 47 Hsp70 genes were identified in B. napus (A_nA_nC_nC_n genome), including 22 genes from A_n subgenome and 25 genes from C_n subgenome. Meanwhile, 29 and 20 Hsp70 genes were explored in B. rapa (A_rA_r genome) and B. oleracea (C_oC_o genome), respectively. Based on phylogenetic analysis, 114 Hsp70 proteins derived from B. napus, B. rapa, B. oleracea and Arabidopsis thaliana, were divided into 6 subfamilies containing 16 A_r-A_n and 11 C_o-C_n reliable orthologous pairs. The homology and synteny analysis indicated whole genome triplication and segmental duplication may be the major contributor for the expansion of Hsp70 gene family. Intron gain of BnHsp70 genes and domain loss of BnHsp70 proteins also were found in B. napus, associating with intron evolution and module evolution of proteins after allopolyploidization. In addition, transcriptional profiles analyses indicated that expression patterns of most BnHsp70 genes were tissue-specific. Moreover, Hsp70 orthologs exhibited different expression patterns in the same tissue and C_n subgenome biased expression was observed in leaf. These findings contribute to exploration of the evolutionary adaptation of polyploidy and will facilitate further application of BnHsp70 gene functions.

An immunoproteomic approach to identify antigenic proteins in Nocardia farcinica IFM 10152

Nocardia farcinica is the etiological agent of nocardiosis, leading to serious pulmonary or systemic infections. To uncover virulence factors and early diagnostic markers, secreted proteins of N. farcinica IFM 10152 were analyzed using an immunoproteome-based approach. A total of 5 proteins were identified by matrix-assisted laser desorption (MALDI-TOF-MS). Bioinformatic analyses showed that the identified proteins were involved in defense against the host innate immune system and required for pathogenesis. All proteins were expressed in E. coli and antigenicity was analyzed with Western blot. To our knowledge, these proteins with antigenicity were identified for the first time in N. farcinica and they may help elucidate the pathogenesis underlying Nocardia and provide potential future diagnostic markers.

Plant evolution and environmental adaptation unveiled by long-read whole-genome sequencing of Spirodela

Aquatic plants have to adapt to the environments distinct from where land plants grow. A critical aspect of adaptation is the dynamics of sequence repeats, not resolved in older sequencing platforms due to incomplete and fragmented genome assemblies from short reads. Therefore, we used PacBio long-read sequencing of the Spirodela polyrhiza genome, reaching a 44-fold increase of contiguity with an N50 (a median of contig lengths) of 831 kb and filling 95.4% of gaps left from the previous version. Reconstruction of repeat regions indicates that sequentially nested long terminal repeat (LTR) retrotranspositions occur early in monocot evolution, featured with both prokaryote-like gene-rich regions and eukaryotic repeat islands. Protein-coding genes are reduced to 18,708 gene models supported by 492,435 high-quality full-length PacBio complementary DNA (cDNA) sequences. Different from land plants, the primitive architecture of Spirodela's adventitious roots and lack of lateral roots and root hairs are consistent with dispensable functions of nutrient absorption. Disease-resistant genes encoding antimicrobial peptides and dirigent proteins are expanded by tandem duplications. Remarkably, disease-resistant genes are not only amplified, but also highly expressed, consistent with low levels of 24-nucleotide (nt) small interfering RNA (siRNA) that silence the immune system of land plants, thereby protecting Spirodela against a wide spectrum of pathogens and pests. The long-read sequence information not only sheds light on plant evolution and adaptation to the environment, but also facilitates applications in bioenergy and phytoremediation.

Bioinformatic and Molecular Analysis of Inverse Autotransporters from Escherichia coli

Escherichia coli is one of the most prevalent facultative anaerobes of the human gut. E. coli normally exists as a harmless commensal but can also cause disease following the acquisition of genes that enhance its pathogenicity. Adhesion is an important first step in colonization of the host and is mediated by an array of cell surface components. In E. coli, these include a family of adhesins secreted by the type V secretion system. Here, we identified and characterized new proteins from an emerging subclass of the type V secretion system known as the inverse autotransporters (IATs). We found that IAT-encoding genes are present in a wide range of strains and showed that three novel IATs were localized on the E. coli cell surface and mediated biofilm formation. Overall, this study provides new insight into the prevalence, function, and regulation of IATs in E. coli.

Proteins secreted by the type V secretion system possess multiple functions, including the capacity to mediate adhesion, aggregation, and biolfilm formation. The type V secretion system can be divided into five subclasses, one of which is the type Ve system. Proteins of the type Ve secretion system are also referred to as inverse autotransporters (IATs). In this study, we performed an in silico analysis of 126 completely sequenced Escherichia coli genomes available in the NCBI database and identified several distinct IAT-encoding gene families whose distribution varied throughout the E. coli phylogeny. The genes included three characterized IATs (intimin, fdeC, and yeeJ) and four uncharacterized IATs (here named iatA, iatB, iatC, and iatD). The four iat genes were cloned from the completely sequenced environmental E. coli strain SMS-3-5 and characterized. Three of these IAT proteins (IatB, IatC, and IatD) were expressed at the cell surface and possessed the capacity to mediate biofilm formation in a recombinant E. coli K-12 strain. Further analysis of the iatB gene, which showed a unique association with extraintestinal E. coli strains, suggested that its regulation is controlled by the LeuO global regulator. Overall, this study provides new data describing the prevalence, sequence variation, domain structure, function, and regulation of IATs found in E. coli.

IMPORTANCEEscherichia coli is one of the most prevalent facultative anaerobes of the human gut. E. coli normally exists as a harmless commensal but can also cause disease following the acquisition of genes that enhance its pathogenicity. Adhesion is an important first step in colonization of the host and is mediated by an array of cell surface components. In E. coli, these include a family of adhesins secreted by the type V secretion system. Here, we identified and characterized new proteins from an emerging subclass of the type V secretion system known as the inverse autotransporters (IATs). We found that IAT-encoding genes are present in a wide range of strains and showed that three novel IATs were localized on the E. coli cell surface and mediated biofilm formation. Overall, this study provides new insight into the prevalence, function, and regulation of IATs in E. coli.

Gene prioritization using Bayesian matrix factorization with genomic and phenotypic side information

Motivation

Most gene prioritization methods model each disease or phenotype individually, but this fails to capture patterns common to several diseases or phenotypes. To overcome this limitation, we formulate the gene prioritization task as the factorization of a sparsely filled gene-phenotype matrix, where the objective is to predict the unknown matrix entries. To deliver more accurate gene-phenotype matrix completion, we extend classical Bayesian matrix factorization to work with multiple side information sources. The availability of side information allows us to make non-trivial predictions for genes for which no previous disease association is known.

Results

Our gene prioritization method can innovatively not only integrate data sources describing genes, but also data sources describing Human Phenotype Ontology terms. Experimental results on our benchmarks show that our proposed model can effectively improve accuracy over the well-established gene prioritization method, Endeavour. In particular, our proposed method offers promising results on diseases of the nervous system; diseases of the eye and adnexa; endocrine, nutritional and metabolic diseases; and congenital malformations, deformations and chromosomal abnormalities, when compared to Endeavour.

Availability and implementation

The Bayesian data fusion method is implemented as a Python/C++ package: https://github.com/jaak-s/macau. It is also available as a Julia package: https://github.com/jaak-s/BayesianDataFusion.jl. All data and benchmarks generated or analyzed during this study can be downloaded at https://owncloud.esat.kuleuven.be/index.php/s/UGb89WfkZwMYoTn.

Supplementary information

Supplementary data are available at Bioinformatics online.

Identification of histone methylation modifiers and their expression patterns during somatic embryogenesis in Hevea brasiliensis

Histone methylation plays a crucial role in various biological processes, from heterochromatin formation to transcriptional regulation. Currently, no information is available regarding histone methylation modifiers in the important rubber-producing plant Hevea brasiliensis. Here, we identified 47 histone methyltransferase (HMT) genes and 25 histone demethylase (HDM) genes as possible members of the histone methylation modifiers in the rubber tree genome. According to the structural features of HMT and HDM, the HbHMTs were classified into two groups (HbPRMs and HbSDGs), the HbHDMs have two groups (HbLSDs and HbJMJs). Expression patterns were analyzed in five different tissues and at different phases of somatic embryogenesis. HbSDG10, 21, 25, 33, HbJMJ2, 18, 20 were with high expression at different phases of somatic embryogenesis. HbSDG10,14, 20, 21, 33 and HbPRMT4 were expressed highly in anther, HbSDG14, 20, 21, 22, 23, 33, 35 and HbPRMT1 HbJMJ7 and HbLSD1, 2, 3, 4 showed high expression levels in callus. HbSDG1, 7, 10, 13, 14, 18, 19, 21, 22, 23, 35, HbPRMT1, 8, HbJMJ5, 7, 11, 16, 20 and HbLSD2, 3, 4 were expressed highly in somatic embryo. HbSDG10, 21, 25, 33, HbLSD2, 3 were expressed highly in bud of regenerated plant. The analyses reveal that HbHMTs and HbHDMs exhibit different expression patterns at different phases during somatic embryogenesis, implying that some HbHMTs and HbHDMs play important roles during somatic embryogenesis. This study provide fundamental information for further studies on histone methylation in Hevea brasiliensis.

Computational and experimental elucidation of Plasmodium falciparum phosphoethanolamine methyltransferase inhibitors: Pivotal drug target

INTRODUCTION

Plasmodium falciparum synthesizes phosphatidylcholine for the membrane development through serine decarboxylase-phosphoethanolamine methyltransferase pathway for growth in human host. Phosphoethanolamine-methyltransferase (PfPMT) is a crucial enzyme for the synthesis of phosphocholine which is a precursor for phosphatidylcholine synthesis and is considered as a pivotal drug target as it is absent in the host. The inhibition of PfPMT may kill malaria parasite and hence is being considered as potential target for rational antimalarial drug designing.

METHODS

In this study, we have used computer aided drug designing (CADD) approaches to establish potential PfPMT inhibitors from Asinex compound library virtually screened for ADMET and the docking affinity. The selected compounds were tested for in-vitro schizonticidal, gametocidal and cytotoxicity activity. Nontoxic compounds were further studied for PfPMT enzyme specificity and antimalarial efficacy for P. berghei in albino mice model.

RESULTS

Our results have identified two nontoxic PfPMT competitive inhibitors ASN.1 and ASN.3 with better schizonticidal and gametocidal activity which were found to inhibit PfPMT at IC50 1.49μM and 2.31μM respectively. The promising reduction in parasitaemia was found both in orally (50 & 10 mg/kg) and intravenous (IV) (5& 1 mg/kg) however, the better growth inhibition was found in intravenous groups.

CONCLUSION

We report that the compounds containing Pyridinyl-Pyrimidine and Phenyl-Furan scaffolds as the potential inhibitors of PfPMT and thus may act as promising antimalarial inhibitor candidates which can be further optimized and used as leads for template based antimalarial drug development.

Transcriptome profiling of Gossypium arboreum during fiber initiation and the genome-wide identification of trihelix transcription factors

Cotton fiber initiation is the first step in fiber development, and it determines the yield. Here, genome-wide transcriptome profiling of Gossypium arboreum was performed to determine the molecular basis of cotton fiber initiation. A comparison of the transcriptomes of fiber-bearing ovules at -0.5, 0, 0.5, 1, 1.5, 2, 2.5 and 3 d post-anthesis detected 12,049 differentially expressed genes that mainly participated in ribosome, carbon metabolism and amino acid biosynthesis pathways. Genes encoding alcohol dehydrogenase 1 and hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase, involving in fatty acid degradation and flavonoid biosynthesis, were enriched. Furthermore, 1049 differentially expressed transcription factors were identified. Among these, 17 were trihelix family transcription factors, which play important roles in plant development and responses to biotic and abiotic stresses. In total, 52 full-length trihelix genes, named as GaGTs, were identified in G. arboreum and located in 12 of the 13 cotton chromosomes. Transcriptomic data and a quantitative real-time PCR analysis indicated that several GaGTs were significantly induced during fiber initiation in G. arboreum. Thus, the genome-wide comprehensive analysis of gene expression in G. arboreum fiber initiation will serve as a useful resource for unraveling the functions of specific genes. The phylogenetic relationships and expression analyses of the G. arboreum trihelix genes established a solid foundation for future comprehensive functional analyses of the GaGTs.

Structural and Functional Characterization of YdjI, an Aldolase of Unknown Specificity in Escherichia coli K12

The ydj gene cluster is found in 80% of sequenced Escherichia coli genomes and other closely related species in the human microbiome. On the basis of the annotations of the enzymes located in this cluster, it is expected that together they catalyze the catabolism of an unknown carbohydrate. The focus of this investigation is on YdjI, which is in the ydj gene cluster of E. coli K-12. It is predicted to be a class II aldolase of unknown function. Here we describe a structural and functional characterization of this enzyme. YdjI catalyzes the hydrogen/deuterium exchange of the pro-S hydrogen at C3 of dihydroxyacetone phosphate (DHAP). In the presence of DHAP, YdjI catalyzes an aldol condensation with a variety of aldo sugars. YdjI shows a strong preference for higher-order (seven-, eight-, and nine-carbon) monosaccharides with specific hydroxyl stereochemistries and a negatively charged terminus (carboxylate or phosphate). The best substrate is l-arabinuronic acid with an apparent k_cat of 3.0 s^-1. The product, l-glycero-l-galacto-octuluronate-1-phosphate, has a k_cat/K_m value of 2.1 × 10³ M^-1 s^-1 in the retro-aldol reaction with YdjI. This is the first recorded synthesis of l-glycero-l-galacto-octuluronate-1-phosphate and six similar carbohydrates. The crystal structure of YdjI, determined to a nominal resolution of 1.75 Å (Protein Data Bank entry 6OFU ), reveals unusual positions for two arginine residues located near the active site. Computational docking was utilized to distinguish preferable binding orientations for l-glycero-l-galacto-octuluronate-1-phosphate. These results indicate a possible alternative binding orientation for l-glycero-l-galacto-octuluronate-1-phosphate compared to that observed in other class II aldolases, which utilize shorter carbohydrate molecules.

Tetrahymena Comparative Genomics Database (TCGD): a community resource for Tetrahymena

Ciliates are a large and diverse group of unicellular organisms characterized by having the following two distinct type of nuclei within a single cell: micronucleus (MIC) and macronucleus (MAC). Although the genomes of several ciliates in different groups have been sequenced, comparative genomics data for multiple species within a ciliate genus are not yet available. Here we collected the genome information and comparative genomics analysis results for 10 species in the Tetrahymena genus, including the previously sequenced model organism Tetrahymena thermophila and 9 newly sequenced species, and constructed a genus-level comparative analysis platform, the Tetrahymena Comparative Genomics Database (TCGD). Genome sequences, transcriptomic data, gene models, functional annotation, ortholog groups and synteny maps were built into this database and a user-friendly interface was developed for searching, visualizing and analyzing these data. In summary, the TCGD (http://ciliate.ihb.ac.cn) will be an important and useful resource for the ciliate research community.

PamulDB: a comprehensive genomic resource for the study of human- and animal-pathogenic Pasteurella multocida

Pasteurella multocida can infect a wide range of host, including humans and animals of economic importance. Genomics studies on the pathogen have produced a large amount of omics data, which are deposited in GenBank but lacks a dedicated and comprehensive resource for further analysis and integration so that need to be brought together centrally in a coherent and systematic manner. Here we have collected the genomic data for 176 P. multocida strains that are categorized into 11 host groups and 9 serotype groups, and developed the open-access P. multocida Database (PamulDB) to make this resource readily available. The PamulDB implements and integrates Chado for genome data management, Drupal for web content management, and bioinformatics tools like NCBI BLAST, HMMER, PSORTb and OrthoMCL for data analysis. All the P. multocida genomes have been further annotated for search and analysis of homologous sequence, phylogeny, gene ontology, transposon, protein subcellular localization and secreted protein. Transcriptomic data of P. multocida are also selectively adopted for gene expression analysis. The PamulDB has been developing and improving to better aid researchers with identifying and classifying of pathogens, dissecting mechanisms of the pathogen infection and host response.

Functional Characterization of YdjH, a Sugar Kinase of Unknown Specificity in Escherichia coli K12

The ydj gene cluster is annotated to catalyze the catabolism of an unknown carbohydrate. Previously, YdjI, a class II aldolase, was shown to catalyze the retro-aldol cleavage of l-glycero-l-galacto-octuluronate-1-phosphate into DHAP and l-arabinuronate. In this report, the functional characterization of YdjH is presented. YdjH catalyzes the phosphorylation of 2-keto-monosaccharides at the C1 hydroxyl group with a substrate profile significantly more stringent than that of YdjI. Similar to YdjI, YdjH shows a strong preference for higher-order monosaccharides (seven to nine carbons) with a carboxylate terminus. The best substrate was determined to be l-glycero-l-galacto-octuluronate, yielding l-glycero-l-galacto-octuluronate-1-phosphate with a k_cat of 16 s^-1 and a k_cat/K_m of 2.1 × 10⁴ M^-1 s^-1. This is apparently the first reported example of kinase activity with eight-carbon monosaccharides. Two crystal structures of YdjH were previously determined to 2.15 and 1.8 Å resolution (Protein Data Bank entries 3H49 and 3IN1 ). We present an analysis of the active site layout and use computational docking to identify potential key residues in the binding of l-glycero-l-galacto-octuluronate.

Why do eukaryotic proteins contain more intrinsically disordered regions?

Intrinsic disorder is more abundant in eukaryotic than prokaryotic proteins. Methods predicting intrinsic disorder are based on the amino acid sequence of a protein. Therefore, there must exist an underlying difference in the sequences between eukaryotic and prokaryotic proteins causing the (predicted) difference in intrinsic disorder. By comparing proteins, from complete eukaryotic and prokaryotic proteomes, we show that the difference in intrinsic disorder emerges from the linker regions connecting Pfam domains. Eukaryotic proteins have more extended linker regions, and in addition, the eukaryotic linkers are significantly more disordered, 38% vs. 12-16% disordered residues. Next, we examined the underlying reason for the increase in disorder in eukaryotic linkers, and we found that the changes in abundance of only three amino acids cause the increase. Eukaryotic proteins contain 8.6% serine; while prokaryotic proteins have 6.5%, eukaryotic proteins also contain 5.4% proline and 5.3% isoleucine compared with 4.0% proline and ≈ 7.5% isoleucine in the prokaryotes. All these three differences contribute to the increased disorder in eukaryotic proteins. It is tempting to speculate that the increase in serine frequencies in eukaryotes is related to regulation by kinases, but direct evidence for this is lacking. The differences are observed in all phyla, protein families, structural regions and type of protein but are most pronounced in disordered and linker regions. The observation that differences in the abundance of three amino acids cause the difference in disorder between eukaryotic and prokaryotic proteins raises the question: Are amino acid frequencies different in eukaryotic linkers because the linkers are more disordered or do the differences cause the increased disorder?

Intrinsic disorder is essential for various functions in eukaryotic cells and is a signature of eukaryotic proteins. Here, we try to understand the origin of the difference in disorder between eukaryotic and prokaryotic proteins. We show that eukaryotic proteins contain more extended linker regions and that these linker regions are significantly more disordered. Further, we show, for the first time, that the difference in disorder originates from a systematic difference in amino acid frequencies between eukaryotic and prokaryotic proteins. Three amino acids contribute to the difference in disorder; serine and proline are more abundant in eukaryotic linkers, while isoleucine is less frequent. These shifts in frequencies are observed in all phyla, protein families, structural regions and type of protein but are most pronounced in disordered and linker regions. It is tempting to speculate that the increase in serine frequencies in eukaryotes is related to regulation by kinases, but direct evidence for this is lacking. Anyhow the widespread of the shifts in abundance indicates that the differences are ancient and caused be some yet not fully understood selective difference acting on eukaryotic and prokaryotic proteins.

PMut: a web-based tool for the annotation of pathological variants on proteins, 2017 update

We present here a full update of the PMut predictor, active since 2005 and with a large acceptance in the field of predicting Mendelian pathological mutations. PMut internal engine has been renewed, and converted into a fully featured standalone training and prediction engine that not only powers PMut web portal, but that can generate custom predictors with alternative training sets or validation schemas. PMut Web portal allows the user to perform pathology predictions, to access a complete repository of pre-calculated predictions, and to generate and validate new predictors. The default predictor performs with good quality scores (MCC values of 0.61 on 10-fold cross validation, and 0.42 on a blind test with SwissVar 2016 mutations). The PMut portal is freely accessible at http://mmb.irbbarcelona.org/PMut. A complete help and tutorial is available at http://mmb.irbbarcelona.org/PMut/help.

The Bologna Annotation Resource (BAR 3.0): improving protein functional annotation

BAR 3.0 updates our server BAR (Bologna Annotation Resource) for predicting protein structural and functional features from sequence. We increase data volume, query capabilities and information conveyed to the user. The core of BAR 3.0 is a graph-based clustering procedure of UniProtKB sequences, following strict pairwise similarity criteria (sequence identity ≥40% with alignment coverage ≥90%). Each cluster contains the available annotation downloaded from UniProtKB, GO, PFAM and PDB. After statistical validation, GO terms and PFAM domains are cluster-specific and annotate new sequences entering the cluster after satisfying similarity constraints. BAR 3.0 includes 28 869 663 sequences in 1 361 773 clusters, of which 22.2% (22 241 661 sequences) and 47.4% (24 555 055 sequences) have at least one validated GO term and one PFAM domain, respectively. 1.4% of the clusters (36% of all sequences) include PDB structures and the cluster is associated to a hidden Markov model that allows building template-target alignment suitable for structural modeling. Some other 3 399 026 sequences are singletons. BAR 3.0 offers an improved search interface, allowing queries by UniProtKB-accession, Fasta sequence, GO-term, PFAM-domain, organism, PDB and ligand/s. When evaluated on the CAFA2 targets, BAR 3.0 largely outperforms our previous version and scores among state-of-the-art methods. BAR 3.0 is publicly available and accessible at http://bar.biocomp.unibo.it/bar3.

Genome of the African cassava whitefly Bemisia tabaci and distribution and genetic diversity of cassava-colonizing whiteflies in Africa

The whitefy Bemisia tabaci, a species complex consisting of many morphologically indistinguishable species divided into distinct clades, is one of the most globally important agricultural pests and plant virus vectors. Cassava-colonizing B. tabaci transmits viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). Half of all cassava plants in Africa are affected by these viral diseases, resulting in annual production losses of more than US$ 1 billion. Here we report the draft genome of the cassava whitefly B. tabaci Sub-Saharan Africa - East and Central Africa (SSA-ECA), the super-abundant population that has been associated with the rapid spread of viruses causing the pandemics of CMD and CBSD. The SSA-ECA genome assembled from Illumina short reads has a total size of 513.7 Mb and a scaffold N50 length of 497 kb, and contains 15,084 predicted protein-coding genes. Phylogenetic analysis suggests that SSA-ECA diverged from MEAM1 around 5.26 million years ago. A comprehensive genetic analysis of cassava-colonizing B. tabaci in Africa was also conducted, in which a total of 243 whitefly specimens were collected from 18 countries representing all major cassava-growing regions in the continent and genotyped using NextRAD sequencing. Population genomic analyses confirmed the existence of six major populations linked by gene flow and inferred the distribution patterns of these populations across the African continent. The genome of SSA-ECA and the genetic findings provide valuable resources and guidance to facilitate whitefly research and the development of strategies to control cassava viral diseases spread by whiteflies.

The OST-complex as target for RNAi-based pest control in Nilaparvata lugens

RNAi-based pest control strategies are emerging as environment friendly and species-specific alternatives for the use of conventional pesticides. Because N-glycosylation is important for many biological processes, such as growth and development, the early steps of protein N-glycosylation are promising targets for an RNAi-based pest control strategy. Through injection of dsRNAs, the expression of the catalytic subunits of the oligosaccharyl transferase complex was efficiently silenced in nymphs of the notorious rice pest insect Nilaparvata lugens. Silencing of both STT3 isoforms resulted in a high mortality of the N. lugens nymphs. However, our data reveals the occurrence of a functional redundancy between the two isoforms when silencing only one of the isoforms. These observations confirm the potential to use the early genes in the N-glycosylation pathway as targets for an RNAi-based pest control strategy. In addition, the existence of a functional redundancy between the two STT3 isoforms presents a factor which one must take into account when designing RNAi-based approaches.

Time-resolved dissection of the molecular crosstalk driving Fusarium head blight in wheat provides new insights into host susceptibility determinism

Fungal plant diseases are controlled by a complex molecular dialogue that involves pathogen effectors able to manipulate plant susceptibility factors at the earliest stages of the interaction. By probing the wheat-Fusarium graminearum pathosystem, we profiled the coregulations of the fungal and plant proteins shaping the molecular responses of a 96-hr-long infection's dynamics. Although no symptoms were yet detectable, fungal biomass swiftly increased along with an extremely diverse set of secreted proteins and candidate effectors supposed to target key plant organelles. Some showed to be early accumulated during the interaction or already present in spores, otherwise stored in germinating spores and detectable in an in vitro F. graminearum exudate. Wheat responses were swiftly set up and were evidenced before any visible symptom. Significant wheat protein abundance changes co-occurred along with the accumulation of putative secreted fungal proteins and predicted effectors. Regulated wheat proteins were closely connected to basal cellular processes occurring during spikelet ontogeny, and particular coregulation patterns were evidenced between chloroplast proteins and fungal proteins harbouring a predicted chloroplast transit peptide. The described plant and fungal coordinated responses provide a resourceful set of data and expand our understanding of the wheat-F. graminearum interaction.

Co-complex protein membership evaluation using Maximum Entropy on GO ontology and InterPro annotation

Motivation

Protein–protein interactions (PPI) play a crucial role in our understanding of protein function and biological processes. The standardization and recording of experimental findings is increasingly stored in ontologies, with the Gene Ontology (GO) being one of the most successful projects. Several PPI evaluation algorithms have been based on the application of probabilistic frameworks or machine learning algorithms to GO properties. Here, we introduce a new training set design and machine learning based approach that combines dependent heterogeneous protein annotations from the entire ontology to evaluate putative co-complex protein interactions determined by empirical studies.

Results

PPI annotations are built combinatorically using corresponding GO terms and InterPro annotation. We use a S.cerevisiae high-confidence complex dataset as a positive training set. A series of classifiers based on Maximum Entropy and support vector machines (SVMs), each with a composite counterpart algorithm, are trained on a series of training sets. These achieve a high performance area under the ROC curve of ≤0.97, outperforming go2ppi—a previously established prediction tool for protein-protein interactions (PPI) based on Gene Ontology (GO) annotations.

Availability and implementation

https://github.com/ima23/maxent-ppi

Supplementary information

Supplementary data are available at Bioinformatics online.

Metagenomic next-generation sequencing of samples from pediatric febrile illness in Tororo, Uganda

Febrile illness is a major burden in African children, and non-malarial causes of fever are uncertain. In this retrospective exploratory study, we used metagenomic next-generation sequencing (mNGS) to evaluate serum, nasopharyngeal, and stool specimens from 94 children (aged 2–54 months) with febrile illness admitted to Tororo District Hospital, Uganda. The most common microbes identified were Plasmodium falciparum (51.1% of samples) and parvovirus B19 (4.4%) from serum; human rhinoviruses A and C (40%), respiratory syncytial virus (10%), and human herpesvirus 5 (10%) from nasopharyngeal swabs; and rotavirus A (50% of those with diarrhea) from stool. We also report the near complete genome of a highly divergent orthobunyavirus, tentatively named Nyangole virus, identified from the serum of a child diagnosed with malaria and pneumonia, a Bwamba orthobunyavirus in the nasopharynx of a child with rash and sepsis, and the genomes of two novel human rhinovirus C species. In this retrospective exploratory study, mNGS identified multiple potential pathogens, including 3 new viral species, associated with fever in Ugandan children.

CAPRI enables comparison of evolutionarily conserved RNA interacting regions

RNA-protein complexes play essential regulatory roles at nearly all levels of gene expression. Using in vivo crosslinking and RNA capture, we report a comprehensive RNA-protein interactome in a metazoan at four levels of resolution: single amino acids, domains, proteins and multisubunit complexes. We devise CAPRI, a method to map RNA-binding domains (RBDs) by simultaneous identification of RNA interacting crosslinked peptides and peptides adjacent to such crosslinked sites. CAPRI identifies more than 3000 RNA proximal peptides in Drosophila and human proteins with more than 45% of them forming new interaction interfaces. The comparison of orthologous proteins enables the identification of evolutionary conserved RBDs in globular domains and intrinsically disordered regions (IDRs). By comparing the sequences of IDRs through evolution, we classify them based on the type of motif, accumulation of tandem repeats, conservation of amino acid composition and high sequence divergence.

Comprehensive characterisation of RNA-protein interactions requires different levels of resolution. Here, the authors present an integrated mass spectrometry-based approach that allows them to define the Drosophila RNA-protein interactome from the level of multisubunit complexes down to the RNA-binding amino acid.

Laying the foundation for genomically-based risk assessment in chronic myeloid leukemia

Outcomes for patients with chronic myeloid leukemia (CML) have substantially improved due to advances in drug development and rational treatment intervention strategies. Despite these significant advances there are still unanswered questions on patient management regarding how to more reliably predict treatment failure at the time of diagnosis and how to select frontline tyrosine kinase inhibitor (TKI) therapy for optimal outcome. The BCR-ABL1 transcript level at diagnosis has no established prognostic impact and cannot guide frontline TKI selection. BCR-ABL1 mutations are detected in ~50% of TKI resistant patients but are rarely responsible for primary resistance. Other resistance mechanisms are largely uncharacterized and there are no other routine molecular testing strategies to facilitate the evaluation and further stratification of TKI resistance. Advances in next-generation sequencing technology has aided the management of a growing number of other malignancies, enabling the incorporation of somatic mutation profiles in diagnosis, classification, and prognostication. A largely unexplored area in CML research is whether expanded genomic analysis at diagnosis, resistance, and disease transformation can enhance patient management decisions, as has occurred for other cancers. The aim of this article is to review publications that reported mutated cancer-associated genes in CML patients at various disease phases. We discuss the frequency and type of such variants at initial diagnosis and at the time of treatment failure and transformation. Current limitations in the evaluation of mutants and recommendations for future reporting are outlined. The collective evaluation of mutational studies over more than a decade suggests a limited set of cancer-associated genes are indeed recurrently mutated in CML and some at a relatively high frequency. Genomic studies have the potential to lay the foundation for improved diagnostic risk classification according to clinical and genomic risk, and to enable more precise early identification of TKI resistance.

Combining learning and constraints for genome-wide protein annotation

BACKGROUND

The advent of high-throughput experimental techniques paved the way to genome-wide computational analysis and predictive annotation studies. When considering the joint annotation of a large set of related entities, like all proteins of a certain genome, many candidate annotations could be inconsistent, or very unlikely, given the existing knowledge. A sound predictive framework capable of accounting for this type of constraints in making predictions could substantially contribute to the quality of machine-generated annotations at a genomic scale.

RESULTS

We present OCELOT, a predictive pipeline which simultaneously addresses functional and interaction annotation of all proteins of a given genome. The system combines sequence-based predictors for functional and protein-protein interaction (PPI) prediction with a consistency layer enforcing (soft) constraints as fuzzy logic rules. The enforced rules represent the available prior knowledge about the classification task, including taxonomic constraints over each GO hierarchy (e.g. a protein labeled with a GO term should also be labeled with all ancestor terms) as well as rules combining interaction and function prediction. An extensive experimental evaluation on the Yeast genome shows that the integration of prior knowledge via rules substantially improves the quality of the predictions. The system largely outperforms GoFDR, the only high-ranking system at the last CAFA challenge with a readily available implementation, when GoFDR is given access to intra-genome information only (as OCELOT), and has comparable or better results (depending on the hierarchy and performance measure) when GoFDR is allowed to use information from other genomes. Our system also compares favorably to recent methods based on deep learning.

Unprecedented reorganization of holocentric chromosomes provides insights into the enigma of lepidopteran chromosome evolution

Chromosome evolution presents an enigma in the mega-diverse Lepidoptera. Most species exhibit constrained chromosome evolution with nearly identical haploid chromosome counts and chromosome-level gene collinearity among species more than 140 million years divergent. However, a few species possess radically inflated chromosomal counts due to extensive fission and fusion events. To address this enigma of constraint in the face of an exceptional ability to change, we investigated an unprecedented reorganization of the standard lepidopteran chromosome structure in the green-veined white butterfly (Pieris napi). We find that gene content in P. napi has been extensively rearranged in large collinear blocks, which until now have been masked by a haploid chromosome number close to the lepidopteran average. We observe that ancient chromosome ends have been maintained and collinear blocks are enriched for functionally related genes suggesting both a mechanism and a possible role for selection in determining the boundaries of these genome-wide rearrangements.

Phytoplasma-induced Changes in the Acetylome and Succinylome of Paulownia tomentosa Provide Evidence for Involvement of Acetylated Proteins in Witches' Broom Disease

Lysine acetylation and succinylation are post-translational modifications of proteins that have been shown to play roles in plants response to pathogen infection. Phytoplasma infection can directly alter multiple metabolic processes in the deciduous plant Paulownia and lead to Paulownia witches' broom (PaWB) disease, the major cause of Paulownia mortality worldwide. However, the extent and function of lysine aceylation and succinylation during phytoplasma infection have yet to be explored. Here, we investigated the changes in the proteome, acetylome, and succinylome of phytoplasma-infected Paulownia tomentosa seedlings using quantitative mass spectrometry. In total, we identified 8963 proteins, 2893 acetylated proteins (5558 acetylation sites), and 1271 succinylated proteins (1970 succinylation sites), with 425 (533 sites) simultaneously acetylated and succinylated. Comparative analysis revealed that 276 proteins, 546 acetylated proteins (741 acetylation sites) and 5 succinylated proteins (5 succinylation sites) were regulated in response to phytoplasma infection, suggesting that acetylation may be more important than succinylation in PaWB. Enzymatic assays showed that acetylation of specific sites in protochlorophyllide reductase and RuBisCO, key enzymes in chlorophyll and starch biosynthesis, respectively, modifies their activity in phytoplasma-infected seedlings. On the basis of these results, we propose a model to elucidate the molecular mechanism of responses to PaWB and offer a resource for functional studies on the effects of acetylation on protein function.

A novel thermostable and halophilic thioredoxin reductase from the Red Sea Atlantis II hot brine pool

The highly extreme conditions of the lower convective layer in the Atlantis II (ATII) Deep brine pool of the Red Sea make it an ideal environment for the search for novel enzymes that can function under extreme conditions. In the current study, we isolated a novel sequence of a thioredoxin reductase (TrxR) enzyme from the metagenomic dataset established from the microbial community that resides in the lower convective layer of Atlantis II. The gene was cloned, expressed and characterized for redox activity, halophilicity, and thermal stability. The isolated thioredoxin reductase (ATII-TrxR) was found to belong to the high-molecular-weight class of thioredoxin reductases. A search for conserved domains revealed the presence of an extra domain (Crp) in the enzyme sequence. Characterization studies of ATII-TrxR revealed that the enzyme was halophilic (maintained activity at 4 M NaCl), thermophilic (optimum temperature was 65°C) and thermostable (60% of its activity was retained at 70°C). Additionally, the enzyme utilized NADH in addition to NADPH as an electron donor. In conclusion, a novel thermostable and halophilic thioredoxin reductase has been isolated with a unique sequence that adapts to the harsh conditions of the brine pools making this protein a good candidate for biological research and industrial applications.