MBGD update 2018: microbial genome database based on hierarchical orthology relations covering closely related and distantly related comparisons

The complete genome sequence of unculturable Mycoplasma faucium obtained through clinical metagenomic next-generation sequencing

Introduction

Diagnosing Mycoplasma faucium poses challenges, and it's unclear if its rare isolation is due to infrequent occurrence or its fastidious nutritional requirements.

Methods

This study analyzes the complete genome sequence of M. faucium, obtained directly from the pus of a sternum infection in a lung transplant patient using metagenomic sequencing.

Results

Genome analysis revealed limited therapeutic options for the M. faucium infection, primarily susceptibility to tetracyclines. Three classes of mobile genetic elements were identified: two new insertion sequences, a new prophage (phiUMCG-1), and a species-specific variant of a mycoplasma integrative and conjugative element (MICE). Additionally, a Type I Restriction-Modification system was identified, featuring 5'-terminally truncated hsdS pseudogenes with overlapping repeats, indicating the potential for forming alternative hsdS variants through recombination.

Conclusion

This study represents the first-ever acquisition of a complete circularized bacterial genome directly from a patient sample obtained from invasive infection of a primary sterile site using culture-independent, PCR-free clinical metagenomics.

InParanoiDB 9: Ortholog Groups for Protein Domains and Full-Length Proteins

Prediction of orthologs is an important bioinformatics pursuit that is frequently used for inferring protein function and evolutionary analyses. The InParanoid database is a well known resource of ortholog predictions between a wide variety of organisms. Although orthologs have historically been inferred at the level of full-length protein sequences, many proteins consist of several independent protein domains that may be orthologous to domains in other proteins in a way that differs from the full-length protein case. To be able to capture all types of orthologous relations, conventional full-length protein orthologs can be complemented with orthologs inferred at the domain level. We here present InParanoiDB 9, covering 640 species and providing orthologs for both protein domains and full-length proteins. InParanoiDB 9 was built using the faster InParanoid-DIAMOND algorithm for orthology analysis, as well as Domainoid and Pfam to infer orthologous domains. InParanoiDB 9 is based on proteomes from 447 eukaryotes, 158 bacteria and 35 archaea, and includes over one billion predicted ortholog groups. A new website has been built for the database, providing multiple search options as well as visualization of groups of orthologs and orthologous domains. This release constitutes a major upgrade of the InParanoid database in terms of the number of species as well as the new capability to operate on the domain level. InParanoiDB 9 is available at https://inparanoidb.sbc.su.se/.

MADGAN:A microbe-disease association prediction model based on generative adversarial networks

Researches have demonstrated that microorganisms are indispensable for the nutrition transportation, growth and development of human bodies, and disorder and imbalance of microbiota may lead to the occurrence of diseases. Therefore, it is crucial to study relationships between microbes and diseases. In this manuscript, we proposed a novel prediction model named MADGAN to infer potential microbe-disease associations by combining biological information of microbes and diseases with the generative adversarial networks. To our knowledge, it is the first attempt to use the generative adversarial network to complete this important task. In MADGAN, we firstly constructed different features for microbes and diseases based on multiple similarity metrics. And then, we further adopted graph convolution neural network (GCN) to derive different features for microbes and diseases automatically. Finally, we trained MADGAN to identify latent microbe-disease associations by games between the generation network and the decision network. Especially, in order to prevent over-smoothing during the model training process, we introduced the cross-level weight distribution structure to enhance the depth of the network based on the idea of residual network. Moreover, in order to validate the performance of MADGAN, we conducted comprehensive experiments and case studies based on databases of HMDAD and Disbiome respectively, and experimental results demonstrated that MADGAN not only achieved satisfactory prediction performances, but also outperformed existing state-of-the-art prediction models.

ProPan: a comprehensive database for profiling prokaryotic pan-genome dynamics

Compared with conventional comparative genomics, the recent studies in pan-genomics have provided further insights into species genomic dynamics, taxonomy and identification, pathogenicity and environmental adaptation. To better understand genome characteristics of species of interest and to fully excavate key metabolic and resistant genes and their conservations and variations, here we present ProPan (https://ngdc.cncb.ac.cn/propan), a public database covering 23 archaeal species and 1,481 bacterial species (in a total of 51,882 strains) for comprehensively profiling prokaryotic pan-genome dynamics. By analyzing and integrating these massive datasets, ProPan offers three major aspects for the pan-genome dynamics of the species of interest: 1) the evaluations of various species' characteristics and composition in pan-genome dynamics; 2) the visualization of map association, the functional annotation and presence/absence variation for all contained species' gene clusters; 3) the typical characteristics of the environmental adaptation, including resistance genes prediction of 126 substances (biocide, antimicrobial drug and metal) and evaluation of 31 metabolic cycle processes. Besides, ProPan develops a very user-friendly interface, flexible retrieval and multi-level real-time statistical visualization. Taken together, ProPan will serve as a weighty resource for the studies of prokaryotic pan-genome dynamics, taxonomy and identification as well as environmental adaptation.

Nitric Oxide Production and Effects in Group B Streptococcus Chorioamnionitis

Intrauterine infection, or chorioamnionitis, due to group B Streptococcus (GBS) is a common cause of miscarriage and preterm birth. To cause chorioamnionitis, GBS must bypass maternal-fetal innate immune defenses including nitric oxide (NO), a microbicidal gas produced by nitric oxide synthases (NOS). This study examined placental NO production and its role in host-pathogen interactions in GBS chorioamnionitis. In a murine model of ascending GBS chorioamnionitis, placental NOS isoform expression quantified by RT-qPCR revealed a four-fold expression increase in inducible NOS, no significant change in expression of endothelial NOS, and decreased expression of neuronal NOS. These NOS expression results were recapitulated ex vivo in freshly collected human placental samples that were co-incubated with GBS. Immunohistochemistry of wild type C57BL/6 murine placentas with GBS chorioamnionitis demonstrated diffuse inducible NOS expression with high-expression foci in the junctional zone and areas of abscess. Pregnancy outcomes between wild type and inducible NOS-deficient mice did not differ significantly although wild type dams had a trend toward more frequent preterm delivery. We also identified possible molecular mechanisms that GBS uses to survive in a NO-rich environment. In vitro exposure of GBS to NO resulted in dose-dependent growth inhibition that varied by serovar. RNA-seq on two GBS strains with distinct NO resistance phenotypes revealed that both GBS strains shared several detoxification pathways that were differentially expressed during NO exposure. These results demonstrate that the placental immune response to GBS chorioamnionitis includes induced NO production and indicate that GBS activates conserved stress pathways in response to NO exposure.

Predicting Microbe-Disease Association Based on Multiple Similarities and LINE Algorithm

Numerous microbes have been found to have vital impacts on human health through affecting biological processes. Therefore, exploring potential associations between microbes and diseases will promote the understanding and diagnosis of diseases. In this study, we present a novel computational model, named MSLINE, to infer potential microbe-disease associations by integrating Multiple Similarities and Large-scale Information Network Embedding (LINE) based on known associations. Specifically, on the basis of known microbe-disease associations from the Human Microbe-Disease Association Database, we first increase the known associations by collecting proven associations from existing literatures. We then construct a microbe-disease heterogeneous network (MDHN) by integrating known associations and multiple similarities (including Gaussian interaction profile kernel similarity, microbe function similarity, disease semantic similarity and disease-symptom similarity). After that, we implement random walk and LINE algorithm on MDHN to learn its structure information. Finally, we score the microbe-disease associations according to the structure information for every nodes. In the Leave-one-out cross validation and 5-fold cross validation, MSLINE performs better compared to other existing methods. Moreover, case studies of different diseases proved that MSLINE could predict the potential microbe-disease associations efficiently.

Emerging frontiers in microbe-mediated pesticide remediation: Unveiling role of omics and In silico approaches in engineered environment

The overuse of pesticides for augmenting agriculture productivity always comes at the cost of environment, biodiversity, and human health and has put the land, water, and environmental footprints under severe threat throughout the globe. Underpinning and maximizing the microbiome functions in pesticide-contaminated environments has become a prerequisite for a sustainable environment and resilient agriculture. It is imperative to elucidate the metabolic network of the microbial communities and environmental variables at the contaminated site to predict the best strategy for remediation and soil microbe-pesticide interactions. High throughput next-generation sequencing and in silico analysis allow us to identify and discern the members and characteristics of core microbiomes at the contaminated site. Integration of modern high throughput multi-omics investigations and informatics pipelines provide novel approaches and pathways to capitalize on the core microbiomes for enhancing environmental functioning and mitigation. The role of eco-genomics tools in visualising the microbial network, taxonomy, functional potential, and environmental variables in contaminated habitats is discussed in this review. The integrated role of the potential microbe identification as individual or consortia, mechanistic approach for pesticide degradation, identification of responsible enzymes/genes, and in silico approach is emphasized for the prospects of the area.

Bioinformatic tools support decision-making in plant disease management

Food loss due to pathogens is a major concern in agriculture, requiring the need for advanced disease detection and prevention measures to minimize pathogen damage to plants. Novel bioinformatic tools have opened doors for the low-cost rapid identification of pathogens and prevention of disease. The number of these tools is growing fast and a comprehensive and comparative summary of these resources is currently lacking. Here, we review all current bioinformatic tools used to identify the mechanisms of pathogen pathogenicity, plant resistance protein identification, and the detection and treatment of plant disease. We compare functionality, data volume, data sources, performance, and applicability of all tools to provide a comprehensive toolbox for researchers in plant disease management.

Complete Genome Sequence of Ferrigenium kumadai An22, a Microaerophilic Iron-Oxidizing Bacterium Isolated from a Paddy Field Soil

Ferrigenium kumadai An22^T (= JCM 30584^T = NBRC 112974^T = ATCC TSD-51^T) is a microaerophilic iron oxidizer isolated from paddy field soil and belongs to the family Gallionellaceae. Here, we report the complete genome sequence of F. kumadai An22^T, which was obtained from the hybrid data of Oxford Nanopore long-read and Illumina short-read sequencing.

Analysis of a photosynthetic cyanobacterium rich in internal membrane systems via gradient profiling by sequencing (Grad-seq)

Although regulatory small RNAs have been reported in photosynthetic cyanobacteria, the lack of clear RNA chaperones involved in their regulation poses a conundrum. Here, we analyzed the full complement of cellular RNAs and proteins using gradient profiling by sequencing (Grad-seq) in Synechocystis 6803. Complexes with overlapping subunits such as the CpcG1-type versus the CpcL-type phycobilisomes or the PsaK1 versus PsaK2 photosystem I pre(complexes) could be distinguished, supporting the high quality of this approach. Clustering of the in-gradient distribution profiles followed by several additional criteria yielded a short list of potential RNA chaperones that include an YlxR homolog and a cyanobacterial homolog of the KhpA/B complex. The data suggest previously undetected complexes between accessory proteins and CRISPR-Cas systems, such as a Csx1-Csm6 ribonucleolytic defense complex. Moreover, the exclusive association of either RpoZ or 6S RNA with the core RNA polymerase complex and the existence of a reservoir of inactive sigma-antisigma complexes is suggested. The Synechocystis Grad-seq resource is available online at https://sunshine.biologie.uni-freiburg.de/GradSeqExplorer/ providing a comprehensive resource for the functional assignment of RNA-protein complexes and multisubunit protein complexes in a photosynthetic organism.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

GIMICA: host genetic and immune factors shaping human microbiota

Besides the environmental factors having tremendous impacts on the composition of microbial community, the host factors have recently gained extensive attentions on their roles in shaping human microbiota. There are two major types of host factors: host genetic factors (HGFs) and host immune factors (HIFs). These factors of each type are essential for defining the chemical and physical landscapes inhabited by microbiota, and the collective consideration of both types have great implication to serve comprehensive health management. However, no database was available to provide the comprehensive factors of both types. Herein, a database entitled 'Host Genetic and Immune Factors Shaping Human Microbiota (GIMICA)' was constructed. Based on the 4257 microbes confirmed to inhabit nine sites of human body, 2851 HGFs (1368 single nucleotide polymorphisms (SNPs), 186 copy number variations (CNVs), and 1297 non-coding ribonucleic acids (RNAs)) modulating the expression of 370 microbes were collected, and 549 HIFs (126 lymphocytes and phagocytes, 387 immune proteins, and 36 immune pathways) regulating the abundance of 455 microbes were also provided. All in all, GIMICA enabled the collective consideration not only between different types of host factor but also between the host and environmental ones, which is freely accessible without login requirement at: https://idrblab.org/gimica/.

COG database update: focus on microbial diversity, model organisms, and widespread pathogens

The Clusters of Orthologous Genes (COG) database, also referred to as the Clusters of Orthologous Groups of proteins, was created in 1997 and went through several rounds of updates, most recently, in 2014. The current update, available at https://www.ncbi.nlm.nih.gov/research/COG, substantially expands the scope of the database to include complete genomes of 1187 bacteria and 122 archaea, typically, with a single genome per genus. In addition, the current version of the COGs includes the following new features: (i) the recently deprecated NCBI's gene index (gi) numbers for the encoded proteins are replaced with stable RefSeq or GenBank\ENA\DDBJ coding sequence (CDS) accession numbers; (ii) COG annotations are updated for >200 newly characterized protein families with corresponding references and PDB links, where available; (iii) lists of COGs grouped by pathways and functional systems are added; (iv) 266 new COGs for proteins involved in CRISPR-Cas immunity, sporulation in Firmicutes and photosynthesis in cyanobacteria are included; and (v) the database is made available as a web page, in addition to FTP. The current release includes 4877 COGs. Future plans include further expansion of the COG collection by adding archaeal COGs (arCOGs), splitting the COGs containing multiple paralogs, and continued refinement of COG annotations.

Lipid Pathway Databases with a Focus on Algae

Pathways of lipid biosynthesis are highly complex and have been established in model organisms such as Arabidopsis thaliana and Chlamydomonas reinhardtii, whereas in other organisms, we need bioinformatic tools to map individual enzymes onto reference pathways. In this chapter, we explain representative tools that are useful in identifying algal orthologs of lipid biosynthetic enzymes and finding new enzymes that are possibly involved in the pathway of interest. All descriptions in this chapter refer to in silico (i.e., computer-based) methods rather than laboratory experiments.

Large-scale production of enzymes for biotechnology uses

Enzymes are biocatalysts that speed up the chemical reaction to obtain the final valuable product/s. Biotechnology has revolutionized the use of traditional enzymes to be applicable in industries such as food, beverage, personal and household care, agriculture, bioenergy, pharmaceutical, and various other segments. With respect to the exponential growth of enzymes in biotech industries, it becomes important to highlight the advancements and impact of enzyme technology over recent years. In this review article, we discuss the existing and emerging production approaches, applications, developments, and global need for enzymes. Special emphasis is given to the predominantly utilized hydrolytic microbial enzymes in industrial bioprocesses.

Helicobacter pylori diversification during chronic infection within a single host generates sub-populations with distinct phenotypes

Helicobacter pylori chronically infects the stomach of approximately half of the world's population. Manifestation of clinical diseases associated with H. pylori infection, including cancer, is driven by strain properties and host responses; and as chronic infection persists, both are subject to change. Previous studies have documented frequent and extensive within-host bacterial genetic variation. To define how within-host diversity contributes to phenotypes related to H. pylori pathogenesis, this project leverages a collection of 39 clinical isolates acquired prospectively from a single subject at two time points and from multiple gastric sites. During the six years separating collection of these isolates, this individual, initially harboring a duodenal ulcer, progressed to gastric atrophy and concomitant loss of acid secretion. Whole genome sequence analysis identified 1,767 unique single nucleotide polymorphisms (SNPs) across isolates and a nucleotide substitution rate of 1.3x10-4 substitutions/site/year. Gene ontology analysis identified cell envelope genes among the genes with excess accumulation of nonsynonymous SNPs (nSNPs). A maximum likelihood tree based on genetic similarity clusters isolates from each time point separately. Within time points, there is segregation of subgroups with phenotypic differences in bacterial morphology, ability to induce inflammatory cytokines, and mouse colonization. Higher inflammatory cytokine induction in recent isolates maps to shared polymorphisms in the Cag PAI protein, CagY, while rod morphology in a subgroup of recent isolates mapped to eight mutations in three distinct helical cell shape determining (csd) genes. The presence of subgroups with unique genetic and phenotypic properties suggest complex selective forces and multiple niches within the stomach during chronic infection.

The Quest for Orthologs benchmark service and consensus calls in 2020

The identification of orthologs—genes in different species which descended from the same gene in their last common ancestor—is a prerequisite for many analyses in comparative genomics and molecular evolution. Numerous algorithms and resources have been conceived to address this problem, but benchmarking and interpreting them is fraught with difficulties (need to compare them on a common input dataset, absence of ground truth, computational cost of calling orthologs). To address this, the Quest for Orthologs consortium maintains a reference set of proteomes and provides a web server for continuous orthology benchmarking (http://orthology.benchmarkservice.org). Furthermore, consensus ortholog calls derived from public benchmark submissions are provided on the Alliance of Genome Resources website, the joint portal of NIH-funded model organism databases.

Predicting human microbe-disease associations via graph attention networks with inductive matrix completion

MOTIVATION

human microbes play a critical role in an extensive range of complex human diseases and become a new target in precision medicine. In silico methods of identifying microbe-disease associations not only can provide a deep insight into understanding the pathogenic mechanism of complex human diseases but also assist pharmacologists to screen candidate targets for drug development. However, the majority of existing approaches are based on linear models or label propagation, which suffers from limitations in capturing nonlinear associations between microbes and diseases. Besides, it is still a great challenge for most previous methods to make predictions for new diseases (or new microbes) with few or without any observed associations.

RESULTS

in this work, we construct features for microbes and diseases by fully exploiting multiply sources of biomedical data, and then propose a novel deep learning framework of graph attention networks with inductive matrix completion for human microbe-disease association prediction, named GATMDA. To our knowledge, this is the first attempt to leverage graph attention networks for this important task. In particular, we develop an optimized graph attention network with talking-heads to learn representations for nodes (i.e. microbes and diseases). To focus on more important neighbours and filter out noises, we further design a bi-interaction aggregator to enforce representation aggregation of similar neighbours. In addition, we combine inductive matrix completion to reconstruct microbe-disease associations to capture the complicated associations between diseases and microbes. Comprehensive experiments on two data sets (i.e. HMDAD and Disbiome) demonstrated that our proposed model consistently outperformed baseline methods. Case studies on two diseases, i.e. asthma and inflammatory bowel disease, further confirmed the effectiveness of our proposed model of GATMDA.

AVAILABILITY

python codes and data set are available at: https://github.com/yahuilong/GATMDA.

CONTACT

luojiawei@hnu.edu.cn.

Contents, Construction Methods, Data Resources, and Functions Comparative Analysis of Bacteria Databases

Many bacterial-related databases are developed to meet the researchers' needs of analysis and search for a number of bacterial information. However, these databases have different data resources, construction methods, data formats, and analysis tools. It's difficult for researchers to select appropriate databases and analysis tools to promote their researches. In the paper, we compared the contents, construction methods, data sources, update frequency, scope and scale of data, analysis tools, and features of nine famous bacterial databases: CARD, EffectiveDB, MBGD, MPD, PATRCI, PHI-base, VFDB, gcMeta and SILVA, and help researchers to better make better use of these databases. In addition, we also hope this review can help researchers develop a more comprehensive database and better tools to meet the needs of researchers.

Advances and Applications in the Quest for Orthologs

Gene families evolve by the processes of speciation (creating orthologs), gene duplication (paralogs), and horizontal gene transfer (xenologs), in addition to sequence divergence and gene loss. Orthologs in particular play an essential role in comparative genomics and phylogenomic analyses. With the continued sequencing of organisms across the tree of life, the data are available to reconstruct the unique evolutionary histories of tens of thousands of gene families. Accurate reconstruction of these histories, however, is a challenging computational problem, and the focus of the Quest for Orthologs Consortium. We review the recent advances and outstanding challenges in this field, as revealed at a symposium and meeting held at the University of Southern California in 2017. Key advances have been made both at the level of orthology algorithm development and with respect to coordination across the community of algorithm developers and orthology end-users. Applications spanned a broad range, including gene function prediction, phylostratigraphy, genome evolution, and phylogenomics. The meetings highlighted the increasing use of meta-analyses integrating results from multiple different algorithms, and discussed ongoing challenges in orthology inference as well as the next steps toward improvement and integration of orthology resources.

A hands-on introduction to querying evolutionary relationships across multiple data sources using SPARQL

The increasing use of Semantic Web technologies in the life sciences, in particular the use of the Resource Description Framework (RDF) and the RDF query language SPARQL, opens the path for novel integrative analyses, combining information from multiple data sources. However, analyzing evolutionary data in RDF is not trivial, due to the steep learning curve required to understand both the data models adopted by different RDF data sources, as well as the equivalent SPARQL constructs required to benefit from this data - in particular, recursive property paths. In this article, we provide a hands-on introduction to querying evolutionary data across several data sources that publish orthology information in RDF, namely: The Orthologous MAtrix (OMA), the European Bioinformatics Institute (EBI) RDF platform, the Database of Orthologous Groups (OrthoDB) and the Microbial Genome Database (MBGD). We present four protocols in increasing order of complexity. In these protocols, we demonstrate through SPARQL queries how to retrieve pairwise orthologs, homologous groups, and hierarchical orthologous groups. Finally, we show how orthology information in different data sources can be compared, through the use of federated SPARQL queries.

A hands-on introduction to querying evolutionary relationships across multiple data sources using SPARQL

The increasing use of Semantic Web technologies in the life sciences, in particular the use of the Resource Description Framework (RDF) and the RDF query language SPARQL, opens the path for novel integrative analyses, combining information from multiple sources. However, analyzing evolutionary data in RDF is not trivial, due to the steep learning curve required to understand both the data models adopted by different RDF data sources, as well as the SPARQL query language. In this article, we provide a hands-on introduction to querying evolutionary data across multiple sources that publish orthology information in RDF, namely: The Orthologous MAtrix (OMA), the European Bioinformatics Institute (EBI) RDF platform, the Database of Orthologous Groups (OrthoDB) and the Microbial Genome Database (MBGD). We present four protocols in increasing order of complexity. In these protocols, we demonstrate through SPARQL queries how to retrieve pairwise orthologs, homologous groups, and hierarchical orthologous groups. Finally, we show how orthology information in different sources can be compared, through the use of federated SPARQL queries.