InterPro in 2019: improving coverage, classification and access to protein sequence annotations

Biofilm and Pathogenesis-Related Proteins in the Foodborne P. fluorescens ITEM 17298 With Distinctive Phenotypes During Cold Storage

In food chain, Pseudomonas spp. cause spoilage by reducing shelf life of fresh products, especially during cold storage, with a high economic burden for industries. However, recent studies have shed new light on health risks occurring when they colonize immunocompromised patient tissues. Likewise to P. aeruginosa, they exhibit antibiotic resistance and biofilm formation, responsible for their spread and persistence in the environment. Biofilm formation might be induced by environmental stresses, such as temperature fluctuations causing physiological and metabolic changes exacerbating food spoilage (by protease and pigment synthesis), and the production of adhesion molecules, chemotactic or underestimated virulence factors. In order to provide a new insight into phenotypic biodiversity of Pseudomonas spoilers isolated from cold stored cheese, in this work 19 Pseudomonas spp. were investigated for biofilm, pigments, exopolysaccharide production and motility at low temperature. Only nine strains showed these phenotypic traits and the blue pigmenting cheese strain P. fluorescens ITEM 17298 was the most distinctive. In addition, this strain decreased the survival probability of infected Galleria mellonella larvae, showing, for the first time, a pathogenic potential. Genomic and proteomic analyses performed on the ITEM 17298 planktonic cells treated or not with lactoferrin derived antibiofilm peptides allowed to reveal specific biofilm related-pathways as well as proteins involved in pathogenesis. Indeed, several genes were found related to signaling system by cGMP-dependent protein kinases, cellulose, rhamnolipid and alginate synthesis, antibiotic resistance, adhesion and virulence factors. The proteome of the untreated ITEM 17298, growing at low temperature, showed that most of the proteins associated with biofilm regulation, pigmentation motility, antibiotic resistance and pathogenecity were repressed, or decreased their levels in comparison to that of the untreated cultures. Thus, the results of this work shed light on the complex pathways network allowing psychrotrophic pseudomonads to adapt themselves to food-refrigerated conditions and enhance their spoilage. In addition, the discovery of virulence factors and antibiotic resistance determinants raises some questions about the need to deeper investigate these underestimated bacteria in order to increase awareness and provide input to update legislation on their detection limits in foods.

Bovine Genome Database: new annotation tools for a new reference genome

The Bovine Genome Database (BGD) (http://bovinegenome.org) has been the key community bovine genomics database for more than a decade. To accommodate the increasing amount and complexity of bovine genomics data, BGD continues to advance its practices in data acquisition, curation, integration and efficient data retrieval. BGD provides tools for genome browsing (JBrowse), genome annotation (Apollo), data mining (BovineMine) and sequence database searching (BLAST). To augment the BGD genome annotation capabilities, we have developed a new Apollo plug-in, called the Locus-Specific Alternate Assembly (LSAA) tool, which enables users to identify and report potential genome assembly errors and structural variants. BGD now hosts both the newest bovine reference genome assembly, ARS-UCD1.2, as well as the previous reference genome, UMD3.1.1, with cross-genome navigation and queries supported in JBrowse and BovineMine, respectively. Other notable enhancements to BovineMine include the incorporation of genomes and gene annotation datasets for non-bovine ruminant species (goat and sheep), support for multiple assemblies per organism in the Regions Search tool, integration of additional ontologies and development of many new template queries. To better serve the research community, we continue to focus on improving existing tools, developing new tools, adding new datasets and encouraging researchers to use these resources.

The DNA methylation landscape of giant viruses

DNA methylation is an important epigenetic mark that contributes to various regulations in all domains of life. Giant viruses are widespread dsDNA viruses with gene contents overlapping the cellular world that also encode DNA methyltransferases. Yet, virtually nothing is known about the methylation of their DNA. Here, we use single-molecule real-time sequencing to study the complete methylome of a large spectrum of giant viruses. We show that DNA methylation is widespread, affecting 2/3 of the tested families, although unevenly distributed. We also identify the corresponding viral methyltransferases and show that they are subject to intricate gene transfers between bacteria, viruses and their eukaryotic host. Most methyltransferases are conserved, functional and under purifying selection, suggesting that they increase the viruses' fitness. Some virally encoded methyltransferases are also paired with restriction endonucleases forming Restriction-Modification systems. Our data suggest that giant viruses' methyltransferases are involved in diverse forms of virus-pathogens interactions during coinfections.

modEnrichr: a suite of gene set enrichment analysis tools for model organisms

High-throughput experiments produce increasingly large datasets that are difficult to analyze and integrate. While most data integration approaches focus on aligning metadata, data integration can be achieved by abstracting experimental results into gene sets. Such gene sets can be made available for reuse through gene set enrichment analysis tools such as Enrichr. Enrichr currently only supports gene sets compiled from human and mouse, limiting accessibility for investigators that study other model organisms. modEnrichr is an expansion of Enrichr for four model organisms: fish, fly, worm and yeast. The gene set libraries within FishEnrichr, FlyEnrichr, WormEnrichr and YeastEnrichr are created from the Gene Ontology, mRNA expression profiles, GeneRIF, pathway databases, protein domain databases and other organism-specific resources. Additionally, libraries were created by predicting gene function from RNA-seq co-expression data processed uniformly from the gene expression omnibus for each organism. The modEnrichr suite of tools provides the ability to convert gene lists across species using an ortholog conversion tool that automatically detects the species. For complex analyses, modEnrichr provides API access that enables submitting batch queries. In summary, modEnrichr leverages existing model organism databases and other resources to facilitate comprehensive hypothesis generation through data integration.

CDD/SPARCLE: the conserved domain database in 2020

As NLM's Conserved Domain Database (CDD) enters its 20th year of operations as a publicly available resource, CDD curation staff continues to develop hierarchical classifications of widely distributed protein domain families, and to record conserved sites associated with molecular function, so that they can be mapped onto user queries in support of hypothesis-driven biomolecular research. CDD offers both an archive of pre-computed domain annotations as well as live search services for both single protein or nucleotide queries and larger sets of protein query sequences. CDD staff has continued to characterize protein families via conserved domain architectures and has built up a significant corpus of curated domain architectures in support of naming bacterial proteins in RefSeq. These architecture definitions are available via SPARCLE, the Subfamily Protein Architecture Labeling Engine. CDD can be accessed at https://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml.

The ABCD database: a repository for chemically defined antibodies

The ABCD (for AntiBodies Chemically Defined) database is a repository of sequenced antibodies, integrating curated information about the antibody and its antigen with cross-links to standardized databases of chemical and protein entities. It is freely available to the academic community, accessible through the ExPASy server (https://web.expasy.org/abcd/). The ABCD database aims at helping to improve reproducibility in academic research by providing a unique, unambiguous identifier associated to each antibody sequence. It also allows to determine rapidly if a sequenced antibody is available for a given antigen.

Crystal structure of the Enterococcus faecalis α-N-acetylgalactosaminidase, a member of the glycoside hydrolase family 31

Glycoside hydrolases catalyze the hydrolysis of glycosidic linkages in carbohydrates. The glycoside hydrolase family 31 (GH31) contains α-glucosidase, α-xylosidase, α-galactosidase, and α-transglycosylase. Recent work has expanded the diversity of substrate specificity of GH31 enzymes, and α-N-acetylgalactosaminidases (αGalNAcases) belonging to GH31 have been identified in human gut bacteria. Here, we determined the first crystal structure of a truncated form of GH31 αGalNAcase from the human gut bacterium Enterococcus faecalis. The enzyme has a similar fold to other reported GH31 enzymes and an additional fibronectin type 3-like domain. Additionally, the structure in complex with N-acetylgalactosamine reveals that conformations of the active site residues, including its catalytic nucleophile, change to recognize the ligand. Our structural analysis provides insight into the substrate recognition and catalytic mechanism of GH31 αGalNAcases.

MicrobioLink: An Integrated Computational Pipeline to Infer Functional Effects of Microbiome-Host Interactions

Microbiome–host interactions play significant roles in health and in various diseases including autoimmune disorders. Uncovering these inter-kingdom cross-talks propels our understanding of disease pathogenesis and provides useful leads on potential therapeutic targets. Despite the biological significance of microbe–host interactions, there is a big gap in understanding the downstream effects of these interactions on host processes. Computational methods are expected to fill this gap by generating, integrating, and prioritizing predictions—as experimental detection remains challenging due to feasibility issues. Here, we present MicrobioLink, a computational pipeline to integrate predicted interactions between microbial and host proteins together with host molecular networks. Using the concept of network diffusion, MicrobioLink can analyse how microbial proteins in a certain context are influencing cellular processes by modulating gene or protein expression. We demonstrated the applicability of the pipeline using a case study. We used gut metaproteomic data from Crohn’s disease patients and healthy controls to uncover the mechanisms by which the microbial proteins can modulate host genes which belong to biological processes implicated in disease pathogenesis. MicrobioLink, which is agnostic of the microbial protein sources (bacterial, viral, etc.), is freely available on GitHub.

The human DEPhOsphorylation Database DEPOD: 2019 update

The human Dephosphorylation Database (DEPOD) is a manually curated resource that harbors human phosphatases, their protein and non-protein substrates, dephosphorylation sites and the associated signaling pathways. We report here an update to DEPOD by integrating and/or linking to annotations from 69 other open access databases including disease associations, phosphorylating kinases, protein interactions, and also genome browsers. We also provide tools to visualize protein interactions, protein structures, phosphorylation networks, evolutionary conservation of proteins, dephosphorylation sites, and short linear motifs within various proteins. The updated version of DEPOD contains 254 human phosphatases, 336 protein and 83 non-protein substrates, and 1215 manually curated phosphatase-substrate relationships. In addition, we have improved the data access as all the data in DEPOD can now be easily downloaded in a user-friendly format. With multiple significant improvements, DEPOD continues serving as a key resource for research on phosphatase-kinase networks.

Database URL: www.depod.org

MGnify: the microbiome analysis resource in 2020

MGnify (http://www.ebi.ac.uk/metagenomics) provides a free to use platform for the assembly, analysis and archiving of microbiome data derived from sequencing microbial populations that are present in particular environments. Over the past 2 years, MGnify (formerly EBI Metagenomics) has more than doubled the number of publicly available analysed datasets held within the resource. Recently, an updated approach to data analysis has been unveiled (version 5.0), replacing the previous single pipeline with multiple analysis pipelines that are tailored according to the input data, and that are formally described using the Common Workflow Language, enabling greater provenance, reusability, and reproducibility. MGnify's new analysis pipelines offer additional approaches for taxonomic assertions based on ribosomal internal transcribed spacer regions (ITS1/2) and expanded protein functional annotations. Biochemical pathways and systems predictions have also been added for assembled contigs. MGnify's growing focus on the assembly of metagenomic data has also seen the number of datasets it has assembled and analysed increase six-fold. The non-redundant protein database constructed from the proteins encoded by these assemblies now exceeds 1 billion sequences. Meanwhile, a newly developed contig viewer provides fine-grained visualisation of the assembled contigs and their enriched annotations.

Structural analysis of the meiosis-related protein MS5 reveals non-canonical papain enhancement by cystatin-like folds

MS5 is a meiosis-related protein belonging to the Brassicaceae-specific domain of unknown function family and characterized by the MS5 superfamily domain (MSD). In this study, we elucidated the three-dimensional crystal structure and potential biochemical function of the MSD. It was observed that the MSD adopts a cystatin-like fold, mainly consisting of a central α-helix and four- or five-stranded antiparallel β-sheets that wrap around it. However, unlike cystatins, which inhibit cysteine proteases, the MSD displayed allosteric activation of papain. We believe that our study provides insight into novel mechanisms of proteolytic enzyme regulation and may serve as a basis for functional studies of the MS5 family proteins in plants.

MaGenDB: a functional genomics hub for Malvaceae plants

Malvaceae is a family of flowering plants containing many economically important plant species including cotton, cacao and durian. Recently, the genomes of several Malvaceae species have been decoded, and many omics data were generated for individual species. However, no integrative database of multiple species, enabling users to jointly compare and analyse relevant data, is available for Malvaceae. Thus, we developed a user-friendly database named MaGenDB (http://magen.whu.edu.cn) as a functional genomics hub for the plant community. We collected the genomes of 13 Malvaceae species, and comprehensively annotated genes from different perspectives including functional RNA/protein element, gene ontology, KEGG orthology, and gene family. We processed 374 sets of diverse omics data with the ENCODE pipelines and integrated them into a customised genome browser, and designed multiple dynamic charts to present gene/RNA/protein-level knowledge such as dynamic expression profiles and functional elements. We also implemented a smart search system for efficiently mining genes. In addition, we constructed a functional comparison system to help comparative analysis between genes on multiple features in one species or across closely related species. This database and associated tools will allow users to quickly retrieve large-scale functional information for biological discovery.

Defining Critical Genes During Spherule Remodeling and Endospore Development in the Fungal Pathogen, Coccidioides posadasii

Coccidioides immitis and C. posadasii are soil dwelling dimorphic fungi found in North and South America. Inhalation of aerosolized asexual conidia can result in asymptomatic, acute, or chronic respiratory infection. In the United States there are approximately 350,000 new infections per year. The Coccidioides genus is the only known fungal pathogen to make specialized parasitic spherules, which contain endospores that are released into the host upon spherule rupture. The molecular determinants involved in this key step of infection remain largely elusive as 49% of genes are hypothetical with unknown function. An attenuated mutant strain C. posadasii Δcts2/Δard1/Δcts3 in which chitinase genes 2 and 3 were deleted was previously created for vaccine development. This strain does not complete endospore development, which prevents completion of the parasitic lifecycle. We sought to identify pathways active in the wild-type strain during spherule remodeling and endospore formation that have been affected by gene deletion in the mutant. We compared the transcriptome and volatile metabolome of the mutant Δcts2/Δard1/Δcts3 to the wild-type C735. First, the global transcriptome was compared for both isolates using RNA sequencing. The raw reads were aligned to the reference genome using TOPHAT2 and analyzed using the Cufflinks package. Genes of interest were screened in an in vivo model using NanoString technology. Using solid phase microextraction (SPME) and comprehensive two-dimensional gas chromatography - time-of-flight mass spectrometry (GC × GC-TOFMS) volatile organic compounds (VOCs) were collected and analyzed. Our RNA-Seq analyses reveal approximately 280 significantly differentially regulated transcripts that are either absent or show opposite expression patterns in the mutant compared to the parent strain. This suggests that these genes are tied to networks impacted by deletion and may be critical for endospore development and/or spherule rupture in the wild-type strain. Of these genes, 14 were specific to the Coccidioides genus. We also found that the wild-type and mutant strains differed significantly in their production versus consumption of metabolites, with the mutant displaying increased nutrient scavenging. Overall, our results provide the first targeted list of key genes that are active during endospore formation and demonstrate that this approach can define targets for functional assays in future studies.

The Biotechnological Potential of the Marine Diatom Skeletonema dohrnii to the Elevated Temperature and pCO2 Concentration

Marine diatoms are promising candidates for biotechnological applications, since they contain high-value compounds, naturally. To facilitate the production of these compounds, stress conditions are often preferable; however, challenges remain with respect to maximizing a metabolic potential for the large-scale cultivation. Here, we sequenced the transcriptome of diatom Skeletonema dohrnii under the actual (21 °C, 400 ppm) and elevated (25 °C, 1000 ppm) temperature and pCO₂ condition. Results indicated that cells grown at higher temperature and pCO₂ showed increasing growth rate, pigment composition, and biochemical productivity as did the expression of chlorophyll, carotenoid and bioactive compound related genes or transcripts. Furthermore, performing de novo transcriptome, we identified 32,884 transcript clusters and found 10,974 of them were differentially expressed between these two conditions. Analyzing the functions of differentially expressed transcripts, we found many of them involved in core metabolic and biosynthesis pathways, including chlorophyll metabolism, carotenoid, phenylpropanoid, phenylalanine and tyrosine, and flavonoid biosynthesis was upregulated. Moreover, we here demonstrated that utilizing a unique bio-fixation ability, S. dohrnii is capable of suppressing central carbon metabolism to promote lipid productivity, fatty acid contents and other bioactive compounds under high temperature and pCO₂ treatment. Our study suggests that this S. dohrnii species could be a potential candidate for wide-scale biotechnological applications under elevated temperature and CO₂ conditions.

Characterization of the Porcine CLEC12A and Analysis of Its Expression on Blood Dendritic Cell Subsets

CLEC12A has been proposed as a suitable target for delivering antigen to dendritic cells (DCs) to enhance vaccine efficacy both in human and mouse. In this study, we have characterized the porcine homolog of CLEC12A (poCLEC12A). Using new monoclonal antibodies (mAb), raised against its ectodomain, poCLEC12A was found to be expressed on alveolar macrophages, blood conventional type 1 and type 2 DCs and plasmacytoid DCs, but not on monocytes, T cells, B cells or NK cells, in contrast to its human and murine homologs. Western blot analysis showed that in alveolar macrophages this receptor is expressed both as a monomer and a dimer. After binding to DCs, anti- poCLEC12A mAb was efficiently internalized. No significant changes were observed in TNFα or IFNα secretion by plasmacytoid DCs stimulated with either CpGs (ODN2216) or polyinosinic-polycytidylic acid (poly I:C), upon incubation with mAb. These results provide the basis for future investigations aimed to assess the ability of anti-poCLEC12A mAbs to improve vaccine efficacy by targeting antigen to DCs.

NetGO: improving large-scale protein function prediction with massive network information

Automated function prediction (AFP) of proteins is of great significance in biology. AFP can be regarded as a problem of the large-scale multi-label classification where a protein can be associated with multiple gene ontology terms as its labels. Based on our GOLabeler—a state-of-the-art method for the third critical assessment of functional annotation (CAFA3), in this paper we propose NetGO, a web server that is able to further improve the performance of the large-scale AFP by incorporating massive protein-protein network information. Specifically, the advantages of NetGO are threefold in using network information: (i) NetGO relies on a powerful learning to rank framework from machine learning to effectively integrate both sequence and network information of proteins; (ii) NetGO uses the massive network information of all species (>2000) in STRING (other than only some specific species) and (iii) NetGO still can use network information to annotate a protein by homology transfer, even if it is not contained in STRING. Separating training and testing data with the same time-delayed settings of CAFA, we comprehensively examined the performance of NetGO. Experimental results have clearly demonstrated that NetGO significantly outperforms GOLabeler and other competing methods. The NetGO web server is freely available at http://issubmission.sjtu.edu.cn/netgo/.

INGA 2.0: improving protein function prediction for the dark proteome

Our current knowledge of complex biological systems is stored in a computable form through the Gene Ontology (GO) which provides a comprehensive description of genes function. Prediction of GO terms from the sequence remains, however, a challenging task, which is particularly critical for novel genomes. Here we present INGA 2.0, a new version of the INGA software for protein function prediction. INGA exploits homology, domain architecture, interaction networks and information from the ‘dark proteome’, like transmembrane and intrinsically disordered regions, to generate a consensus prediction. INGA was ranked in the top ten methods on both CAFA2 and CAFA3 blind tests. The new algorithm can process entire genomes in a few hours or even less when additional input files are provided. The new interface provides a better user experience by integrating filters and widgets to explore the graph structure of the predicted terms. The INGA web server, databases and benchmarking are available from URL: https://inga.bio.unipd.it/.

Genome-wide association study leads to novel genetic insights into resistance to Aspergillus flavus in maize kernels

BACKGROUND

Fungus infection in staple grains affects the food storage and threatens food security. The Aspergillus flavus is known to infect multiple grains and produce mycotoxin Aflatoxin B1, which is mutagenic, teratogenic and causes immunosuppression in animals. However, the molecular mechanism of maize resistance to A. flavus is largely unknown.

RESULTS

Here we used corn kernels to investigate resistance genes to A. flavus using genome-wide association study (GWAS) of 313 inbred lines. We characterized the resistance levels of kernels after inoculating with A. flavus. The GWAS with 558,529 SNPs identified four associated loci involving 29 candidate genes that were linked to seed development, resistance or infection, and involved in signal pathways, seed development, germination, dormancy, epigenetic modification, and antimicrobial activity. In addition, a few candidate genes were also associated with several G-protein signaling and phytohormones that might involve in synergistic work conferring different resistance during seed development. Expression of 16 genes out of 29 during kernel development was also associated with resistance levels.

CONCLUSIONS

We characterized the resistance levels of 313 maize kernels after inoculating with A. flavus, and found four associated loci and 16 candidate maize genes. The expressed 16 genes involved in kernel structure and kernel composition most likely contribute to mature maize kernels' resistance to A. flavus, and in particular, in the development of pericarp. The linked candidate genes could be experimentally transformed to validate and manipulate fungal resistance. Thus this result adds value to maize kernels in breeding programs.

Clemen C, Eichinger L. Functional Characterisation of the Autophagy ATG12~5/16 Complex in Dictyostelium discoideum

Macroautophagy, a highly conserved and complex intracellular degradative pathway, involves more than 20 core autophagy (ATG) proteins, among them the hexameric ATG12~5/16 complex, which is part of the essential ubiquitin-like conjugation systems in autophagy. Dictyostelium discoideumatg5 single, atg5/12 double, and atg5/12/16 triple gene knock-out mutant strains displayed similar defects in the conjugation of ATG8 to phosphatidylethanolamine, development, and cell viability upon nitrogen starvation. This implies that ATG5, 12 and 16 act as a functional unit in canonical autophagy. Macropinocytosis of TRITC dextran and phagocytosis of yeast were significantly decreased in ATG5¯ and ATG5¯/12¯ and even further in ATG5¯/12¯/16¯ cells. In contrast, plaque growth on Klebsiella aerogenes was about twice as fast for ATG5¯ and ATG5¯/12¯/16¯ cells in comparison to AX2, but strongly decreased for ATG5¯/12¯ cells. Along this line, phagocytic uptake of Escherichia coli was significantly reduced in ATG5¯/12¯ cells, while no difference in uptake, but a strong increase in membrane association of E. coli, was seen for ATG5¯ and ATG5¯/12¯/16¯ cells. Proteasomal activity was also disturbed in a complex fashion, consistent with an inhibitory activity of ATG16 in the absence of ATG5 and/or ATG12. Our results confirm the essential function of the ATG12~5/16 complex in canonical autophagy, and furthermore are consistent with autophagy-independent functions of the complex and its individual components. They also strongly support the placement of autophagy upstream of the ubiquitin-proteasome system (UPS), as a fully functional UPS depends on autophagy.

Heavy Metal Stress-Associated Proteins in Rice and Arabidopsis: Genome-Wide Identification, Phylogenetics, Duplication, and Expression Profiles Analysis

Heavy metal exposure is a serious environmental stress in plants. However, plants have evolved several strategies to improve their heavy metal tolerance. Heavy metal-associated proteins (HMPs) participate in heavy metal detoxification. Here, we identified 46 and 55 HMPs in rice and Arabidopsis, respectively, and named them OsHMP 1–46 and AtHMP 1–55 according to their chromosomal locations. The HMPs from both plants were divided into six clades based on the characteristics of their heavy metal-associated domains (HMA). The HMP gene structures and motifs varied greatly among the different classifications. The HMPs had high collinearity and were segmentally duplicated. A cis-element analysis revealed that the HMPs may be regulated by different transcription factors. An expression profile analysis disclosed that only eight OsHMPs were constitutive in rice tissues. Of these, the expression of OsHMP37 was far higher than that of the other seven genes while OsHMP28 was expressed exclusively in the roots. For Arabidopsis, nine AtHMPs presented with very high transcript levels in all organs. Most of the selected OsHMPs were differentially expressed in various tissues under different heavy metal stresses. Only OsHMP09, OsHMP18, and OsHMP22 showed higher expression levels in all tissues under different heavy metal stresses. In contrast, most of the selected AtHMPs had nearly constant expression levels in different tissues under various heavy metal stresses. The AtHMP20, AtHMP23, AtHMP25, AtHMP31, AtHMP35, AtHMP46 expression levels under different heavy metal stresses were higher in the leaves and roots. The foregoing discoveries elucidated HMP evolution in monocotyledonous and dicotyledonous plants and may helpful functionally characterize HMPs in the future.

Structural and Proteomic Characterization of the Initiation of Giant Virus Infection

Since their discovery, giant viruses have expanded our understanding of the principles of virology. Due to their gargantuan size and complexity, little is known about the life cycles of these viruses. To answer outstanding questions regarding giant virus infection mechanisms, we set out to determine biomolecular conditions that promote giant virus genome release. We generated four infection intermediates in Samba virus (Mimivirus genus, lineage A) as visualized by cryoelectron microscopy (cryo-EM), cryoelectron tomography (cryo-ET), and scanning electron microscopy (SEM). Each of these four intermediates reflects similar morphology to a stage that occurs in vivo. We show that these genome release stages are conserved in other mimiviruses. Finally, we identified proteins that are released from Samba and newly discovered Tupanvirus through differential mass spectrometry. Our work revealed the molecular forces that trigger infection are conserved among disparate giant viruses. This study is also the first to identify specific proteins released during the initial stages of giant virus infection.

The Whole-Genome Sequencing and Hybrid Assembly of Mytilus coruscus

The hard-shelled mussel (Mytilus coruscus) is an economically important shellfish that has been cultivated for the last decade. Due to over-exploitation, most mussel stocks have dramatically declined. Efforts to study this species' natural distribution, genetics, breeding, and cultivation have been hindered by the lack of a high-quality reference genome. To address this, we produced a hybrid high-quality reference genome of M. coruscus using a long-read platform to assemble the genome and short-read, high-quality technology to accurately correct for sequence errors. The genome was assembled into 10,484 scaffolds, a total length of 1.90 Gb, and a scaffold N50 of 898 kb. Ab initio annotation of the M. coruscus genome assembly identified a total of 42,684 genes. This accurate reference genome of M. coruscus provides an essential resource with the advantage of enabling the genome-scale selective breeding of M. coruscus. More importantly, it will also help in deciphering the speciation and local adaptation of the Mytilus species.

The Pheromone Module SteC-MkkB-MpkB-SteD-HamE Regulates Development, Stress Responses and Secondary Metabolism in Aspergillus fumigatus

In order for eukaryotes to efficiently detect and respond to environmental stimuli, a myriad of protein signaling pathways are utilized. An example of highly conserved signaling pathways in eukaryotes are the mitogen-activated protein kinase (MAPK) pathways. In fungi, MAPK pathways have been shown to regulate a diverse array of biological processes, such as asexual and sexual development, stress responses and the production of secondary metabolites (SMs). In the model fungus Aspergillus nidulans, a MAPK pathway known as the pheromone module is utilized to regulate both development and SM production. This signaling cascade consists of the three kinases SteC, MkkB, and MpkB, as well as the SteD adaptor protein and the HamE scaffold. In this study, homologs of each of these proteins have been identified in the opportunistic human pathogen A. fumigatus. By performing epitope tagging and mass spectrometry experiments, we have shown that these proteins form a pentameric complex, similar to what is observed in A. nidulans. This complex has been shown to assemble in the cytoplasm and MpkB enters the nucleus, where it would presumably interact with various transcription factors. Pheromone module mutant strains exhibit drastic reductions in asexual sporulation, vegetative growth rate and production of SMs, such as gliotoxin. Mutants also display increased sensitivity to cell wall and oxidative stress agents. Overall, these data provide evidence of the existence of a conserved MAP kinase signaling pathway in Aspergillus species and suggest that this pathway is critical for the regulation of fungal development and secondary metabolism.

Type II and type IV toxin-antitoxin systems show different evolutionary patterns in the global Klebsiella pneumoniae population

The Klebsiella pneumoniae species complex includes important opportunistic pathogens which have become public health priorities linked to major hospital outbreaks and the recent emergence of multidrug-resistant hypervirulent strains. Bacterial virulence and the spread of multidrug resistance have previously been linked to toxin-antitoxin (TA) systems. TA systems encode a toxin that disrupts essential cellular processes, and a cognate antitoxin which counteracts this activity. Whilst associated with the maintenance of plasmids, they also act in bacterial immunity and antibiotic tolerance. However, the evolutionary dynamics and distribution of TA systems in clinical pathogens are not well understood. Here, we present a comprehensive survey and description of the diversity of TA systems in 259 clinically relevant genomes of K. pneumoniae. We show that TA systems are highly prevalent with a median of 20 loci per strain. Importantly, these toxins differ substantially in their distribution patterns and in their range of cognate antitoxins. Classification along these properties suggests different roles of TA systems and highlights the association and co-evolution of toxins and antitoxins.

Alarming India-wide phenomenon of antifungal resistance in dermatophytes: A multicentre study

BACKGROUND

An alarming increase in recalcitrant dermatophytosis has been witnessed in India over the past decade. Drug resistance may play a major role in this scenario.

OBJECTIVES

The aim of the present study was to determine the prevalence of in vitro resistance to terbinafine, itraconazole and voriconazole in dermatophytes, and to identify underlying mutations in the fungal squalene epoxidase (SQLE) gene.

PATIENTS/METHODS

We analysed skin samples from 402 patients originating from eight locations in India. Fungi were identified by microbiological and molecular methods, tested for antifungal susceptibility (terbinafine, itraconazole, voriconazole), and investigated for missense mutations in SQLE.

RESULTS

Trichophyton (T.) mentagrophytes internal transcribed spacer (ITS) Type VIII was found in 314 (78%) samples. Eighteen (5%) samples harboured species identified up to the T interdigitale/mentagrophytes complex, and T rubrum was detected in 19 (5%) samples. 71% of isolates were resistant to terbinafine. The amino acid substitution Phe397Leu in the squalene epoxidase of resistant T mentagrophytes was highly prevalent (91%). Two novel substitutions in resistant Trichophyton strains, Ser395Pro and Ser443Pro, were discovered. The substitution Ala448Thr was found in terbinafine-sensitive and terbinafine-resistant isolates but was associated with increased MICs of itraconazole and voriconazole.

CONCLUSIONS

The high frequencies of terbinafine resistance in dermatophytes are worrisome and demand monitoring and further research. Squalene epoxidase substitutions between Leu393 and Ser443 could serve as markers of resistance in the future.

Multi-faceted nature of the tRNA isopentenyltransferase

Transfer RNA isopentenylation an adenine 37 position (A37) is a universal modification known in prokaryotes and eukaryotes. A set of highly homologous enzymes catalyse a series of reactions, leading to tRNA modifications, aimed to increase adaptation to environmental condition through the control of translation efficiency and reading frame maintenance. Transfer RNA-isopentenylation-related (TI-related) functions are well studied in bacteria, mitochondria of yeast and human, but completely unexplored in plants. Transfer RNA-isopentenylation-unrelated (TI-unrelated) functions participate in adaptation to environmental stresses via the regulation of sterol metabolism, gene silencing/suppression and amyloid fibrils formation. TI-unrelated functions are mostly studied in yeast. Finally, the degradation of A37-modified tRNA releases a set of bioactive compounds known as cis-cytokinins. Although all organisms are able to produce cis-cytokinins, its physiological role is still a matter of debates. For several species of bacteria and fungi, cis-cytokinins are known to play a crucial role in pathogenesis. In mammalian and human models cis-cytokinins have tumour-suppressing and anti-inflammation effects. This review aims to summarise current knowledge of the TI-related and TI-unrelated functions and main bioactive by-products of isopentenylated tRNA degradation.

Wheat PP2C-a10 regulates seed germination and drought tolerance in transgenic Arabidopsis

A wheat protein phosphatase PP2C-a10, which interacted with TaDOG1L1 and TaDOG1L4, promoted seed germination and decreased drought tolerance of transgenic Arabidopsis. Seed dormancy and germination are critical to plant fitness. DELAY OF GERMINATION 1 (DOG1) is a quantitative trait locus for dormancy in Arabidopsis thaliana. Some interactions between DOG1 and the type 2C protein phosphatases (PP2Cs) have been reported in Arabidopsis. However, the research on molecular functions and regulations of DOG1Ls and group A PP2Cs in wheat (Triticum aestivum. L), an important crop plant, is rare. In this study, the whole TaDOG1L family was identified. Expression analysis revealed that TaDOG1L2, TaDOG1L4 and TaDOG1L-N2 specially expressed in wheat grains, while others displayed distinct expression patterns. Yeast two-hybrid analysis of TaDOG1Ls and group A TaPP2Cs revealed interaction patterns differed from those in Arabidopsis, and TaDOG1L1 and TaDOG1L4 interacted with TaPP2C-a10. The qRT-PCR analysis showed that TaPP2C-a10 exhibited the highest transcript level in wheat grains. Further investigation showed that ectopic expression of TaPP2C-a10 in Arabidopsis promoted seed germination and decreased sensitivity to ABA during germination stage. Additionally, TaPP2C-a10 transgenic Arabidopsis exhibited decreased tolerance to drought stress. Finally, the phylogenetic analysis indicated that TaPP2C-a10 gene was conserved in angiosperm during evolutionary process. Overall, our results reveal the role of TaPP2C-a10 in seed germination and abiotic stress response, as well as the functional diversity of TaDOG1L family.

Developmental Mechanisms of Fleshy Fruit Diversity in Rosaceae

Rosaceae (the rose family) is an economically important family that includes species prized for high-value fruits and ornamentals. The family also exhibits diverse fruit types, including drupe (peach), pome (apple), drupetum (raspberry), and achenetum (strawberry). Phylogenetic analysis and ancestral fruit-type reconstruction suggest independent evolutionary paths of multiple fleshy fruit types from dry fruits. A recent whole genome duplication in the Maleae/Pyreae tribe (with apple, pear, hawthorn, and close relatives; referred to as Maleae here) may have contributed to the evolution of pome fruit. MADS-box genes, known to regulate floral organ identity, are emerging as important regulators of fruit development. The differential competence of floral organs to respond to fertilization signals may explain the different abilities of floral organs to form fleshy fruit. Future comparative genomics and functional studies in closely related Rosaceae species with distinct fruit types will test hypotheses and provide insights into mechanisms of fleshy fruit diversity. These efforts will be facilitated by the wealth of genome data and resources in Rosaceae.

ECMPride: prediction of human extracellular matrix proteins based on the ideal dataset using hybrid features with domain evidence

Extracellular matrix (ECM) proteins play an essential role in various biological processes in multicellular organisms, and their abnormal regulation can lead to many diseases. For large-scale ECM protein identification, especially through proteomic-based techniques, a theoretical reference database of ECM proteins is required. In this study, based on the experimentally verified ECM datasets and by the integration of protein domain features and a machine learning model, we developed ECMPride, a flexible and scalable tool for predicting ECM proteins. ECMPride achieved excellent performance in predicting ECM proteins, with appropriate balanced accuracy and sensitivity, and the performance of ECMPride was shown to be superior to the previously developed tool. A new theoretical dataset of human ECM components was also established by applying ECMPride to all human entries in the SwissProt database, containing a significant number of putative ECM proteins as well as the abundant biological annotations. This dataset might serve as a valuable reference resource for ECM protein identification.

Hyd ubiquitinates the NF-κB co-factor Akirin to operate an effective immune response in Drosophila

The Immune Deficiency (IMD) pathway in Drosophila melanogaster is activated upon microbial challenge with Gram-negative bacteria to trigger the innate immune response. In order to decipher this nuclear factor κB (NF-κB) signaling pathway, we undertook an in vitro RNAi screen targeting E3 ubiquitin ligases specifically and identified the HECT-type E3 ubiquitin ligase Hyperplastic discs (Hyd) as a new actor in the IMD pathway. Hyd mediated Lys⁶³ (K63)-linked polyubiquitination of the NF-κB cofactor Akirin was required for efficient binding of Akirin to the NF-κB transcription factor Relish. We showed that this Hyd-dependent interaction was required for the transcription of immunity-related genes that are activated by both Relish and Akirin but was dispensable for the transcription of genes that depend solely on Relish. Therefore Hyd is key in NF-κB transcriptional selectivity downstream of the IMD pathway. Drosophila depleted of Akirin or Hyd failed to express the full set of genes encoding immune-induced anti-microbial peptides and succumbed to immune challenges. We showed further that UBR5, the mammalian homolog of Hyd, was also required downstream of the NF-κB pathway for the activation of Interleukin 6 (IL6) transcription by LPS or IL-1β in cultured human cells. Our findings link the action of an E3 ubiquitin ligase to the activation of immune effector genes, deepening our understanding of the involvement of ubiquitination in inflammation and identifying a potential target for the control of inflammatory diseases.

Ubiquitination has been recently identified in pathogenesis and progression of various diseases where inflammation is critical. NF-κB transcription factors are key actors in the transcriptional cascade leading to inflammation as they activate genes with pro- or anti-inflammatory activities. The similarity between the immune pathways in flies and mammals makes Drosophila melanogaster an excellent model to study the innate response. Accordingly, we decided to identify E3 ubiquitin-ligases involved in the regulation of NF-κB pathway, using Drosophila as a model system. A RNAi based screen in immortalized embryonic macrophage-like Drosophila cells points to the HECT-E3 ubiquitin ligase Hyd as a new regulator of the Immune-deficiency (IMD) NF-κB pathway, activated after Gram-negative immune challenge.

More precisely, we showed that Hyd acts at the level of Akirin, an evolutionarily conserved player in the NF-κB pathway, required for the transcription of pro-inflammatory genes, but not for the NF-κB-dependent genes contributing to the down-regulation of inflammation.

In addition, we could show that the human homologue of Hyd (UBR5) acts genetically at the level of human AKIRIN2, pointing to a unique dichotomy between Hyd/Akirin-dependent and -independent gene activation, allowing for the decoupling activation and resolution of inflammation. These results identified UBR5 as a putative target for anti-inflammatory compounds.

Honey-bee-associated prokaryotic viral communities reveal wide viral diversity and a profound metabolic coding potential

This study uses viral-like particle purification and subsequent unbiased genome sequencing to identify prokaryotic viruses associated with Apis mellifera. Interestingly, bacteriophages found in honey bees show a high diversity and span different viral taxa. This diversity sharply contrasts with the state-of-the-art knowledge on the relatively simple bee bacterial microbiome. The identification of multiple auxiliary metabolic genes suggests that these bacteriophages possess the coding potential to intervene in essential microbial pathways related to health and possibly also to disease. This study sheds light on a neglected part of the bee microbiota and opens avenues of in vivo research on the interaction of bacteriophages with their bacterial host, which likely has strongly underappreciated consequences on bee health.

Honey bees (Apis mellifera) produce an enormous economic value through their pollination activities and play a central role in the biodiversity of entire ecosystems. Recent efforts have revealed the substantial influence that the gut microbiota exert on bee development, food digestion, and homeostasis in general. In this study, deep sequencing was used to characterize prokaryotic viral communities associated with honey bees, which was a blind spot in research up until now. The vast majority of the prokaryotic viral populations are novel at the genus level, and most of the encoded proteins comprise unknown functions. Nevertheless, genomes of bacteriophages were predicted to infect nearly every major bee-gut bacterium, and functional annotation and auxiliary metabolic gene discovery imply the potential to influence microbial metabolism. Furthermore, undiscovered genes involved in the synthesis of secondary metabolic biosynthetic gene clusters reflect a wealth of previously untapped enzymatic resources hidden in the bee bacteriophage community.

Novel Target Exploration from Hypothetical Proteins of Klebsiella pneumoniae MGH 78578 Reveals a Protein Involved in Host-Pathogen Interaction

The opportunistic pathogen Klebsiella pneumoniae is a causative agent of several hospital-acquired infections. It has become resistant to a wide range of currently available antibiotics, leading to high mortality rates among patients; this has further led to a demand for novel therapeutic intervention to treat such infections. Using a series of in silico analyses, the present study aims to explore novel drug/vaccine candidates from the hypothetical proteins of K. pneumoniae. A total of 540 proteins were found to be hypothetical in this organism. Analysis of these 540 hypothetical proteins revealed 30 pathogen-specific proteins essential for pathogen survival. A motifs/domain family analysis, similarity search against known proteins, gene ontology, and protein–protein interaction analysis of the shortlisted 30 proteins led to functional assignment for 17 proteins. They were mainly cataloged as enzymes, lipoproteins, stress-induced proteins, transporters, and other proteins (viz., two-component proteins, skeletal proteins and toxins). Among the annotated proteins, 16 proteins, located in the cytoplasm, periplasm, and inner membrane, were considered as potential drug targets, and one extracellular protein was considered as a vaccine candidate. A druggability analysis indicated that the identified 17 drug/vaccine candidates were “novel”. Furthermore, a host–pathogen interaction analysis of these identified target candidates revealed a betaine/carnitine/choline transporters (BCCT) family protein showing interactions with five host proteins. Structure prediction and validation were carried out for this protein, which could aid in structure-based inhibitor design.

Aspartate aminotransferase Rv3722c governs aspartate-dependent nitrogen metabolism in Mycobacterium tuberculosis

Gene rv3722c of Mycobacterium tuberculosis is essential for in vitro growth, and encodes a putative pyridoxal phosphate-binding protein of unknown function. Here we use metabolomic, genetic and structural approaches to show that Rv3722c is the primary aspartate aminotransferase of M. tuberculosis, and mediates an essential but underrecognized role in metabolism: nitrogen distribution. Rv3722c deficiency leads to virulence attenuation in macrophages and mice. Our results identify aspartate biosynthesis and nitrogen distribution as potential species-selective drug targets in M. tuberculosis.

A community proposal to integrate structural bioinformatics activities in ELIXIR (3D-Bioinfo Community)

Structural bioinformatics provides the scientific methods and tools to analyse, archive, validate, and present the biomolecular structure data generated by the structural biology community. It also provides an important link with the genomics community, as structural bioinformaticians also use the extensive sequence data to predict protein structures and their functional sites. A very broad and active community of structural bioinformaticians exists across Europe, and 3D-Bioinfo will establish formal platforms to address their needs and better integrate their activities and initiatives. Our mission will be to strengthen the ties with the structural biology research communities in Europe covering life sciences, as well as chemistry and physics and to bridge the gap between these researchers in order to fully realize the potential of structural bioinformatics. Our Community will also undertake dedicated educational, training and outreach efforts to facilitate this, bringing new insights and thus facilitating the development of much needed innovative applications e.g. for human health, drug and protein design. Our combined efforts will be of critical importance to keep the European research efforts competitive in this respect.

Here we highlight the major European contributions to the field of structural bioinformatics, the most pressing challenges remaining and how Europe-wide interactions, enabled by ELIXIR and its platforms, will help in addressing these challenges and in coordinating structural bioinformatics resources across Europe. In particular, we present recent activities and future plans to consolidate an ELIXIR 3D-Bioinfo Community in structural bioinformatics and propose means to develop better links across the community. These include building new consortia, organising workshops to establish data standards and seeking community agreement on benchmark data sets and strategies. We also highlight existing and planned collaborations with other ELIXIR Communities and other European infrastructures, such as the structural biology community supported by Instruct-ERIC, with whom we have synergies and overlapping common interests.

Fold-Change-Specific Enrichment Analysis (FSEA): Quantification of Transcriptional Response Magnitude for Functional Gene Groups

Gene expression profiling data contains more information than is routinely extracted with standard approaches. Here we present Fold-Change-Specific Enrichment Analysis (FSEA), a new method for functional annotation of differentially expressed genes from transcriptome data with respect to their fold changes. FSEA identifies Gene Ontology (GO) terms, which are shared by the group of genes with a similar magnitude of response, and assesses these changes. GO terms found by FSEA are fold-change-specifically (e.g., weakly, moderately, or strongly) affected by a stimulus under investigation. We demonstrate that many responses to abiotic factors, mutations, treatments, and diseases occur in a fold-change-specific manner. FSEA analyses suggest that there are two prevailing responses of functionally-related gene groups, either weak or strong. Notably, some of the fold-change-specific GO terms are invisible by classical algorithms for functional gene enrichment, Singular Enrichment Analysis (SEA), and Gene Set Enrichment Analysis (GSEA). These are GO terms not enriched compared to the genome background but strictly regulated by a factor within specific fold-change intervals. FSEA analysis of a cancer-related transcriptome suggested that the gene groups with a tightly coordinated response can be the valuable source to search for possible regulators, markers, and therapeutic targets in oncogenic processes. Availability and Implementation: FSEA is implemented as the FoldGO Bioconductor R package and a web-server.

Deep conservation of prion-like composition in the eukaryotic prion-former Pub1/Tia1 family and its relatives

Pub1 protein is an important RNA-binding protein functional in stress granule assembly in budding yeast Saccharomyces cerevisiae and, as its co-ortholog Tia1, in humans. It is unique among proteins in evidencing prion-like aggregation in both its yeast and human forms. Previously, we noted that Pub1/Tia1 was the only protein linked to human disease that has prion-like character and and has demonstrated such aggregation in both species. Thus, we were motivated to probe further into the evolution of the Pub1/Tia1 family (and its close relative Nam8 and its orthologs) to gain a picture of how such a protein has evolved over deep evolutionary time since the last common ancestor of eukaryotes. Here, we discover that the prion-like composition of this protein family is deeply conserved across eukaryotes, as is the prion-like composition of its close relative Nam8/Ngr1. A sizeable minority of protein orthologs have multiple prion-like domains within their sequences (6-20% depending on criteria). The number of RNA-binding RRM domains is conserved at three copies over >86% of the Pub1 family (>71% of the Nam8 family), but proteins with just one or two RRM domains occur frequently in some clades, indicating that these are not due to annotation errors. Overall, our results indicate that a basic scaffold comprising three RNA-binding domains and at least one prion-like region has been largely conserved since the last common ancestor of eukaryotes, providing further evidence that prion-like aggregation may be a very ancient and conserved phenomenon for certain specific proteins.

Pathological, Morphological, Cytogenomic, Biochemical and Molecular Data Support the Distinction between Colletotrichum cigarro comb. et stat. nov. and Colletotrichum kahawae

The genus Colletotrichum has witnessed tremendous variations over the years in the number of species recognized, ranging from 11 to several hundreds. Host-specific fungal species, once the rule, are now the exception, with polyphagous behavior regarded as normal in this genus. The species Colletotrichum kahawae was created to accommodate the pathogens that have the unique ability to infect green developing coffee berries causing the devastating Coffee Berry Disease in Africa, but its close phylogenetic relationship to a polyphagous group of fungi in the C. gloeosporioides species complex led some researchers to regard these pathogens as members of a wider species. In this work we combine pathological, morphological, cytogenomic, biochemical, and molecular data of a comprehensive set of phylogenetically-related isolates to show that the Coffee Berry Disease pathogen forms a separate species, C. kahawae, and also to assign the closely related fungi, previously in C. kahawae subsp. cigarro, to a new species, C. cigarro comb. et stat. nov. This taxonomic clarification provides an opportunity to link phylogeny and functional biology, and additionally enables a much-needed tool for plant pathology and agronomy, associating exclusively C. kahawae to the Coffee Berry Disease pathogen.

Dynamic Surface Proteomes of Allergenic Fungal Conidia

Fungal spores and hyphal fragments play an important role as allergens in respiratory diseases. In this study, we performed trypsin shaving and secretome analyses to identify the surface-exposed proteins and secreted/shed proteins of Aspergillus fumigatus conidia, respectively. We investigated the surface proteome under different conditions, including temperature variation and germination. We found that the surface proteome of resting A. fumigatus conidia is not static but instead unexpectedly dynamic, as evidenced by drastically different surface proteomes under different growth conditions. Knockouts of two abundant A. fumigatus surface proteins, ScwA and CweA, were found to function only in fine-tuning the cell wall stress response, implying that the conidial surface is very robust against perturbations. We then compared the surface proteome of A. fumigatus to other allergy-inducing molds, including Alternaria alternata, Penicillium rubens, and Cladosporium herbarum, and performed comparative proteomics on resting and swollen conidia, as well as secreted proteins from germinating conidia. We detected 125 protein ortholog groups, including 80 with putative catalytic activity, in the extracellular region of all four molds, and 42 nonorthologous proteins produced solely by A. fumigatus. Ultimately, this study highlights the dynamic nature of the A. fumigatus conidial surface and provides targets for future diagnostics and immunotherapy.

Alienimonas chondri sp. nov., a novel planctomycete isolated from the biofilm of the red alga Chondrus crispus

The phylum Planctomycetes comprises bacteria with peculiar and very unique characteristics among prokaryotes. In marine environments, macroalgae biofilms are well known for harboring planctomycetal diversity. Here, we describe a novel isolate obtained from the biofilm of the red alga Chondrus crispus collected at a rocky beach in Porto, Portugal. The novel strain LzC2^T is motile, rosette-forming with spherical- to ovoid-shaped cells. LzC2^T forms magenta- to pinkish-colored colonies in M13 and M14 media. Transmission and scanning electron microscopy observations showed a division by polar and lateral budding. Mother cells are connected to the daughter cells by a tubular neck-like structure. The strain requires salt for growth. Vitamins are not required for growth. Optimal growth occurs from 15 to 30°C and within a pH range from 5.5 to 10.0. Major fatty acids are anteiso-C15:0 (54.2%) and iso-C15:0 (19.5%). Phosphatidylglycerol, diphosphatidylglycerol and an unidentified glycolipid represent the main lipids and menaquinone 6 (MK-6) is the only quinone present. 16S rRNA gene-based phylogenetic analysis supports the affiliation to the phylum Planctomycetes and family Planctomycetaceae, with Alienimonas as the closest relative. Strain LzC2^T shares 97% 16S rRNA gene sequence similarity with Alienimonas californiensis. LzC2^T has a genome size of 5.3 Mb and a G+C content of 68.3%. Genotypic and phenotypic comparison with the closest relatives strongly suggest that LzC2^T (=CECT 30038^T=LMG XXXT) is a new species of the genus Alienimonas, for which we propose the name Alienimonas chondri sp. nov., represented by LzC2^T as type strain. 16S rRNA gene accession number: GenBank=MN757873.1. Genome accession number: GenBank=WTPX00000000.

Bioinformatics analysis of diversity in bacterial glycan chain-termination chemistry and organization of carbohydrate-binding modules linked to ABC transporters

The structures of bacterial cell surface glycans are remarkably diverse. In spite of this diversity, the general strategies used for their assembly are limited. In one of the major processes, found in both Gram-positive and Gram-negative bacteria, the glycan is polymerized in the cytoplasm on a polyprenol lipid carrier and exported from the cytoplasm by an ATP-binding cassette (ABC) transporter. The ABC transporter actively participates in determining the chain length of the glycan substrate, which impacts functional properties of the glycoconjugate products. A subset of these systems employs an additional elaborate glycan capping strategy that dictates the size distribution of the products. The hallmarks of prototypical capped glycan systems are a chain-terminating enzyme possessing a coiled-coil molecular ruler and an ABC transporter possessing a carbohydrate-binding module, which recognizes the glycan cap. To date, detailed investigations are limited to a small number of prototypes, and here, we used our current understanding of these processes for a bioinformatics census of other examples in available genome sequences. This study not only revealed additional instances of existing terminators but also predicted new chemistries as well as systems that diverge from the established prototypes. These analyses enable some new functional hypotheses and offer a roadmap for future research.

Genome analysis of cellulose and hemicellulose degrading Micromonospora sp. CP22

In this study, a bacterial strain CP22 with ability to produce cellulase, xylanase and mannanase was isolated from the oil palm compost. Based on the 16S rRNA gene analysis, the strain was affiliated to genus Micromonospora. To further investigate genes that are related to cellulose and hemicellulose degradation, the genome of strain CP22 was sequenced, annotated and analyzed. The de novo assembled genome of strain CP22 featured a size of 5,856,203 bp with G + C content of 70.84%. Detailed genome analysis on lignocellulose degradation revealed a total of 60 genes consisting of 47 glycoside hydrolase domains and 16 carbohydrate esterase domains predicted to be involved in cellulolytic and hemicellulolytic deconstruction. Particularly, 20 genes encode for cellulases (8 endoglucanases, 3 exoglucanases and 9 β-glucosidases) and 40 genes encode for hemicellulases (15 endo-1,4-β-xylanase, 3 β-xylosidase, 3 α-arabinofuranosidase, 10 acetyl xylan esterase, 6 polysaccharide deacetylase, 1 β-mannanase, 1 β-mannosidase and 1 α-galactosidase). Thirty-two genes encoding carbohydrate-binding modules (CBM) from six different families (CBM2, CBM4, CBM6, CBM9, CBM13 and CBM22) were present in the genome of strain CP22. These CBMs were found in 27 cellulolytic and hemicellulolytic genes, indicating their potential role in enhancing the substrate-binding capability of the enzymes. CBM2 and CBM13 are the major CBMs present in cellulases and hemicellulases (xylanases and mannanases), respectively. Moreover, a GH10 xylanase was found to contain 3 CBMs (1 CBM9 and 2 CBM22) and these CBMs were reported to bind specifically to xylan. This genome-based analysis could facilitate the exploration of this strain for lignocellulosic biomass degradation.

Antibiotics from Haloarchaea: What Can We Learn from Comparative Genomics?

The knowledge of antibiotics produced by Archaea (archaeocins) is still limited. So far, only two types of archaeocins are known: (i) sulfolobicins, produced by the extremely thermophilic Sulfolobus spp. and (ii) haloarcheocins, produced by halophilic archaea. Haloarcheocins were first discovered in the 1980s, but most of their characterisation was solely based on supernatant-based assays. Only a few were successfully purified and sequenced, and even fewer have a proposed biosynthetic model. Furthermore, their mode of action, ecological role and biotechnological potential are still to be explored. Haloarcheocin C8 (HalC8) is the best well-characterised haloarcheocin. We applied an approach of comparative genomics in order to go a step further in the knowledge of their biosynthetic clusters as well as the clusters encoding HalC8-like peptides. These peptides can be classified, at least, into 4 different clades, and there is low gene conservation between them. However, the putative function of some proteins is conserved. These include uncharacterized major facilitator superfamily proteins, transmembrane peptides, DNA-binding transcriptional regulators and proteins with extracellular domains. Our analysis reinforces the association of these proteins with HalC8/HalC8-like biosynthesis. Their functionality is unknown, and, in an era where it is known that haloarchaea are not confined to high salt habitats, the advance in the knowledge of their specialised metabolites will be imperative.

High expression in maize pollen correlates with genetic contributions to pollen fitness as well as with coordinated transcription from neighboring transposable elements

In flowering plants, gene expression in the haploid male gametophyte (pollen) is essential for sperm delivery and double fertilization. Pollen also undergoes dynamic epigenetic regulation of expression from transposable elements (TEs), but how this process interacts with gene expression is not clearly understood. To explore relationships among these processes, we quantified transcript levels in four male reproductive stages of maize (tassel primordia, microspores, mature pollen, and sperm cells) via RNA-seq. We found that, in contrast with vegetative cell-limited TE expression in Arabidopsis pollen, TE transcripts in maize accumulate as early as the microspore stage and are also present in sperm cells. Intriguingly, coordinate expression was observed between highly expressed protein-coding genes and their neighboring TEs, specifically in mature pollen and sperm cells. To investigate a potential relationship between elevated gene transcript level and pollen function, we measured the fitness cost (male-specific transmission defect) of GFP-tagged coding sequence insertion mutations in over 50 genes identified as highly expressed in the pollen vegetative cell, sperm cell, or seedling (as a sporophytic control). Insertions in seedling genes or sperm cell genes (with one exception) exhibited no difference from the expected 1:1 transmission ratio. In contrast, insertions in over 20% of vegetative cell genes were associated with significant reductions in fitness, showing a positive correlation of transcript level with non-Mendelian segregation when mutant. Insertions in maize gamete expressed2 (Zm gex2), the sole sperm cell gene with measured contributions to fitness, also triggered seed defects when crossed as a male, indicating a conserved role in double fertilization, given the similar phenotype previously demonstrated for the Arabidopsis ortholog GEX2. Overall, our study demonstrates a developmentally programmed and coordinated transcriptional activation of TEs and genes in pollen, and further identifies maize pollen as a model in which transcriptomic data have predictive value for quantitative phenotypes.

Non-intrinsic ATP-binding cassette proteins ABCI19, ABCI20 and ABCI21 modulate cytokinin response at the endoplasmic reticulum in Arabidopsis thaliana

The non-intrinsic ABC proteins ABCI20 and ABCI21 are induced by light under HY5 regulation, localize to the ER, and ameliorate cytokinin-driven growth inhibition in young Arabidopsis thaliana seedlings. The plant ATP-binding cassette (ABC) I subfamily (ABCIs) comprises heterogeneous proteins containing any of the domains found in other ABC proteins. Some ABCIs are known to function in basic metabolism and stress responses, but many remain functionally uncharacterized. ABCI19, ABCI20, and ABCI21 of Arabidopsis thaliana cluster together in a phylogenetic tree, and are suggested to be targets of the transcription factor ELONGATED HYPOCOTYL 5 (HY5). Here, we reveal that these three ABCIs are involved in modulating cytokinin responses during early seedling development. The ABCI19, ABCI20 and ABCI21 promoters harbor HY5-binding motifs, and ABCI20 and ABCI21 expression was induced by light in a HY5-dependent manner. abci19 abci20 abci21 triple and abci20 abci21 double knockout mutants were hypersensitive to cytokinin in seedling growth retardation assays, but did not show phenotypic differences from the wild type in either control medium or auxin-, ABA-, GA-, ACC- or BR-containing media. ABCI19, ABCI20, and ABCI21 were expressed in young seedlings and the three proteins interacted with each other, forming a large protein complex at the endoplasmic reticulum (ER) membrane. These results suggest that ABCI19, ABCI20, and ABCI21 fine-tune the cytokinin response at the ER under the control of HY5 at the young seedling stage.

Grad-seq in a Gram-positive bacterium reveals exonucleolytic sRNA activation in competence control

RNA–protein interactions are the crucial basis for many steps of bacterial gene expression, including post‐transcriptional control by small regulatory RNAs (sRNAs). In stark contrast to recent progress in the analysis of Gram‐negative bacteria, knowledge about RNA–protein complexes in Gram‐positive species remains scarce. Here, we used the Grad‐seq approach to draft a comprehensive landscape of such complexes in Streptococcus pneumoniae, in total determining the sedimentation profiles of ~ 88% of the transcripts and ~ 62% of the proteins of this important human pathogen. Analysis of in‐gradient distributions and subsequent tag‐based protein capture identified interactions of the exoribonuclease Cbf1/YhaM with sRNAs that control bacterial competence for DNA uptake. Unexpectedly, the nucleolytic activity of Cbf1 stabilizes these sRNAs, thereby promoting their function as repressors of competence. Overall, these results provide the first RNA/protein complexome resource of a Gram‐positive species and illustrate how this can be utilized to identify new molecular factors with functions in RNA‐based regulation of virulence‐relevant pathways.

The shell matrix of the european thorny oyster, Spondylus gaederopus: microstructural and molecular characterization

Molluscs, the largest marine phylum, display extraordinary shell diversity and sophisticated biomineral architectures. However, mineral-associated biomolecules involved in biomineralization are still poorly characterised. We report the first comprehensive structural and biomolecular study of Spondylus gaederopus, a pectinoid bivalve with a peculiar shell texture. Used since prehistoric times, this is the best-known shell of Europe's cultural heritage. We find that Spondylus microstructure is very poor in mineral-bound organics, which are mostly intercrystalline and concentrated at the interface between structural layers. Using high-resolution liquid chromatography tandem mass spectrometry (LC-MS/MS) we characterized several shell protein fractions, isolated following different bleaching treatments. Several peptides were identified as well as six shell proteins, which display features and domains typically found in biomineralized tissues, including the prevalence of intrinsically disordered regions. It is very likely that these sequences only partially represent the full proteome of Spondylus, considering the lack of genomics data for this genus and the fact that most of the reconstructed peptides do not match with any known shell proteins, representing consequently lineage-specific sequences. This work sheds light onto the shell matrix involved in the biomineralization in spondylids. Our proteomics data suggest that Spondylus has evolved a shell-forming toolkit, distinct from that of other better studied pectinoids - fine-tuned to produce shell structures with high mechanical properties, while limited in organic content. This study therefore represents an important milestone for future studies on biomineralized skeletons and provides the first reference dataset for forthcoming molecular studies of Spondylus archaeological artifacts.

Klebsiella pneumoniae type VI secretion system-mediated microbial competition is PhoPQ controlled and reactive oxygen species dependent

Klebsiella pneumoniae is recognized as an urgent threat to human health due to the increasing isolation of multidrug resistant strains. Hypervirulent strains are a major concern due to their ability to cause life-threating infections in healthy hosts. The type VI secretion system (T6SS) is widely implicated in microbial antagonism, and it mediates interactions with host eukaryotic cells in some cases. In silico search for genes orthologous to T6SS component genes and T6SS effector genes across 700 K. pneumoniae genomes shows extensive diversity in T6SS genes across the K. pneumoniae species. Temperature, oxygen tension, pH, osmolarity, iron levels, and NaCl regulate the expression of the T6SS encoded by a hypervirulent K. pneumoniae strain. Polymyxins and human defensin 3 also increase the activity of the T6SS. A screen for regulators governing T6SS uncover the correlation between the transcription of the T6SS and the ability to kill E. coli prey. Whereas H-NS represses the T6SS, PhoPQ, PmrAB, Hfq, Fur, RpoS and RpoN positively regulate the T6SS. K. pneumoniae T6SS mediates intra and inter species bacterial competition. This antagonism is only evident when the prey possesses an active T6SS. The PhoPQ two component system governs the activation of K. pneumoniae T6SS in bacterial competitions. Mechanistically, PhoQ periplasmic domain, and the acid patch within, is essential to activate K. pneumoniae T6SS. Klebsiella T6SS also mediates anti-fungal competition. We have delineated the contribution of each of the individual VgrGs in microbial competition and identified VgrG4 as a T6SS effector. The DUF2345 domain of VgrG4 is sufficient to intoxicate bacteria and yeast. ROS generation mediates the antibacterial effects of VgrG4, and the antitoxin Sel1E protects against the toxic activity of VgrG4. Our findings provide a better understanding of the regulation of the T6SS in bacterial competitions, and place ROS as an early event in microbial competition.

Klebsiella pneumoniae has been singled out as an “urgent threat to human health” due to extremely drug resistant strains. Numerous studies investigate the molecular mechanisms underlying antibiotic resistance in K. pneumoniae, while others dissect the virulence strategies of this pathogen. However, there is still limited knowledge on the fitness of Klebsiella in the environment, and, particularly, the competition of Klebsiella with other species. Here, we demonstrate that Klebsiella exploits the type VI secretion system (T6SS) nanoweapon to kill bacterial competitors and fungi. K. pneumoniae perceives T6SS attacks from bacterial competitors, resulting in retaliation against the aggressive cell. The perception of the attack involved the sensor PhoPQ and led to the up-regulation of the T6SS. We identified one of the toxins deployed by the T6SS to antagonize other microbes and revealed how Klebsiella protects itself from this toxin. Our findings provide a better understanding of the T6SS role in microbial competition and uncover new aspects on how bacteria regulate T6SS-mediated microbial antagonism.

Unraveling 1,4-Butanediol Metabolism in Pseudomonas putida KT2440

Plastics, in all forms, are a ubiquitous cornerstone of modern civilization. Although humanity undoubtedly benefits from the versatility and durability of plastics, they also cause a tremendous burden for the environment. Bio-upcycling is a promising approach to reduce this burden, especially for polymers that are currently not amenable to mechanical recycling. Wildtype P. putida KT2440 is able to grow on 1,4-butanediol as sole carbon source, but only very slowly. Adaptive laboratory evolution (ALE) led to the isolation of several strains with significantly enhanced growth rate and yield. Genome re-sequencing and proteomic analysis were applied to characterize the genomic and metabolic basis of efficient 1,4-butanediol metabolism. Initially, 1,4-butanediol is oxidized to 4-hydroxybutyrate, in which the highly expressed dehydrogenase enzymes encoded within the PP_2674-2680 ped gene cluster play an essential role. The resulting 4-hydroxybutyrate can be metabolized through three possible pathways: (i) oxidation to succinate, (ii) CoA activation and subsequent oxidation to succinyl-CoA, and (iii) beta oxidation to glycolyl-CoA and acetyl-CoA. The evolved strains were both mutated in a transcriptional regulator (PP_2046) of an operon encoding both beta-oxidation related genes and an alcohol dehydrogenase. When either the regulator or the alcohol dehydrogenase is deleted, no 1,4-butanediol uptake or growth could be detected. Using a reverse engineering approach, PP_2046 was replaced by a synthetic promotor (14g) to overexpress the downstream operon (PP_2047-2051), thereby enhancing growth on 1,4-butanediol. This work provides a deeper understanding of microbial 1,4-butanediol metabolism in P. putida, which is also expandable to other aliphatic alpha-omega diols. It enables the more efficient metabolism of these diols, thereby enabling biotechnological valorization of plastic monomers in a bio-upcycling approach.

Wikidata as a knowledge graph for the life sciences

Wikidata is a community-maintained knowledge base that has been assembled from repositories in the fields of genomics, proteomics, genetic variants, pathways, chemical compounds, and diseases, and that adheres to the FAIR principles of findability, accessibility, interoperability and reusability. Here we describe the breadth and depth of the biomedical knowledge contained within Wikidata, and discuss the open-source tools we have built to add information to Wikidata and to synchronize it with source databases. We also demonstrate several use cases for Wikidata, including the crowdsourced curation of biomedical ontologies, phenotype-based diagnosis of disease, and drug repurposing.

Draft genomes of female and male turbot Scophthalmus maximus

Turbot (Scophthalmus maximus) is a commercially important flatfish species in aquaculture. It has a drastic sexual dimorphism, with females growing faster than males. In the present study, we sequenced and de novo assembled female and male turbot genomes. The assembled female genome was 568 Mb (scaffold N50, 6.2 Mb, BUSCO 97.4%), and the male genome was 584 Mb (scaffold N50, 5.9 Mb, BUSCO 96.6%). Using two genetic maps, we anchored female scaffolds representing 535 Mb onto 22 chromosomes. Annotation of the female anchored genome identified 87.8 Mb transposon elements and 20,134 genes. We identified 17,936 gene families, of which 369 gene families were flatfish specific. Phylogenetic analysis showed that the turbot, Japanese flounder and Chinese tongue sole form a clade that diverged from other teleosts approximately 78 Mya. This report of female and male turbot draft genomes and annotated genes provides a new resource for identifying sex determination genes, elucidating the evolution of adaptive traits in flatfish and developing genetic techniques to increase the sustainability of turbot aquaculture.