The IMG/M data management and analysis system v.6.0: new tools and advanced capabilities

Phylogeny and Evolutionary History of Respiratory Complex I Proteins in Melainabacteria

The evolution of oxygenic photosynthesis was one of the most transformative evolutionary events in Earth's history, leading eventually to the oxygenation of Earth's atmosphere and, consequently, the evolution of aerobic respiration. Previous work has shown that the terminal electron acceptors (complex IV) of aerobic respiration likely evolved after the evolution of oxygenic photosynthesis. However, complex I of the respiratory complex chain can be involved in anaerobic processes and, therefore, may have pre-dated the evolution of oxygenic photosynthesis. If so, aerobic respiration may have built upon respiratory chains that pre-date the rise of oxygen in Earth's atmosphere. The Melainabacteria provide a unique opportunity to examine this hypothesis because they contain genes for aerobic respiration but likely diverged from the Cyanobacteria before the evolution of oxygenic photosynthesis. Here, we examine the phylogenies of translated complex I sequences from 44 recently published Melainabacteria metagenome assembled genomes and genomes from other Melainabacteria, Cyanobacteria, and other bacterial groups to examine the evolutionary history of complex I. We find that complex I appears to have been present in the common ancestor of Melainabacteria and Cyanobacteria, supporting the idea that aerobic respiration built upon respiratory chains that pre-date the evolution of oxygenic photosynthesis and the rise of oxygen.

Development of a time-series shotgun metagenomics database for monitoring microbial communities at the Pacific coast of Japan

Although numerous metagenome, amplicon sequencing-based studies have been conducted to date to characterize marine microbial communities, relatively few have employed full metagenome shotgun sequencing to obtain a broader picture of the functional features of these marine microbial communities. Moreover, most of these studies only performed sporadic sampling, which is insufficient to understand an ecosystem comprehensively. In this study, we regularly conducted seawater sampling along the northeastern Pacific coast of Japan between March 2012 and May 2016. We collected 213 seawater samples and prepared size-based fractions to generate 454 subsets of samples for shotgun metagenome sequencing and analysis. We also determined the sequences of 16S rRNA (n = 111) and 18S rRNA (n = 47) gene amplicons from smaller sample subsets. We thereafter developed the Ocean Monitoring Database for time-series metagenomic data ( http://marine-meta.healthscience.sci.waseda.ac.jp/omd/ ), which provides a three-dimensional bird's-eye view of the data. This database includes results of digital DNA chip analysis, a novel method for estimating ocean characteristics such as water temperature from metagenomic data. Furthermore, we developed a novel classification method that includes more information about viruses than that acquired using BLAST. We further report the discovery of a large number of previously overlooked (TAG)n repeat sequences in the genomes of marine microbes. We predict that the availability of this time-series database will lead to major discoveries in marine microbiome research.

Fecal Methylmercury Correlates With Gut Microbiota Taxa in Pacific Walruses (Odobenus rosmarus divergens)

OBJECTIVES

Methylmercury metabolism was investigated in Pacific walruses (Odobenus rosmarus divergens) from St. Lawrence Island, Alaska, United States.

METHODS

Total mercury and methylmercury concentrations were measured in fecal samples and paired colon samples (n = 16 walruses). Gut microbiota composition and diversity were determined using 16S rRNA gene sequencing. Associations between fecal and colon mercury and the 24 most prevalent gut microbiota taxa were investigated using linear models.

RESULTS

In fecal samples, the median values for total mercury, methylmercury, and %methylmercury (of total mercury) were 200 ng/g, 4.7 ng/g, and 2.5%, respectively, while in colon samples, the median values for the same parameters were 28 ng/g, 7.8 ng/g, and 26%, respectively. In fecal samples, methylmercury was negatively correlated with one Bacteroides genus, while members of the Oscillospirales order were positively correlated with both methylmercury and %methylmercury (of total mercury). In colon samples, %methylmercury (of total mercury) was negatively correlated with members of two genera, Romboutsia and Paeniclostridium.

CONCLUSIONS

Median %methylmercury (of total mercury) was 10 times higher in the colon compared to the fecal samples, suggesting that methylmercury was able to pass through the colon into systemic circulation. Fecal total mercury and/or methylmercury concentrations in walruses were comparable to some human studies despite differences in seafood consumption rates, suggesting that walruses excreted less mercury. There are no members (at this time) of the Oscillospirales order which are known to contain the genes to methylate mercury, suggesting the source of methylmercury in the gut was from diet and not in vivo methylation.

Metagenome Sequencing to Explore Phylogenomics of Terrestrial Cyanobacteria

Cyanobacteria are ubiquitous microorganisms with crucial ecosystem functions, yet most knowledge of their biology relates to aquatic taxa. We have constructed metagenomes for 50 taxonomically well-characterized terrestrial cyanobacterial cultures. These data will support phylogenomic studies of evolutionary relationships and gene content among these unique algae and their aquatic relatives.

Description of Nocardioides piscis sp. nov., Sphingomonas piscis sp. nov. and Sphingomonas sinipercae sp. nov., isolated from the intestine of fish species Odontobutis interrupta (Korean spotted sleeper) and Siniperca scherzeri (leopard mandarin fish)

A polyphasic taxonomic approach was used to characterize three novel bacterial strains, designated as HDW12AT, HDW-15BT, and HDW15CT, isolated from the intestine of fish species Odontobutis interrupta or Siniperca scherzeri. All isolates were obligate aerobic, non-motile bacteria, and grew optimally at 30°C. Phylogenetic analysis based on 16S rRNA sequences revealed that strain HDW12AT was a member of the genus Nocardioides, and closely related to Nocardioides allogilvus CFH 30205T (98.9% sequence identities). Furthermore, strains HDW15BT and HDW15CT were members of the genus Sphingomonas, and closely related to Sphingomonas lutea JS5T and Sphingomonas sediminicola Dae 20T (97.1% and 97.9% sequence identities), respectively. Strain HDW12AT contained MK-8 (H4), and strains HDW15BT and HDW15CT contained Q-10 as the respiratory quinone. Major polar lipid components of strain HDW12AT were diphosphatidylglycerol, phosphatidylglycerol, and phosphatidylinositol, and those of strains HDW15BT and HDW15CT were sphingoglycolipid, diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, and phosphatidylcholine. The G + C content of strains HDW12AT, HDW15BT, and HDW15CT were 69.7, 63.3, and 65.5%, respectively. The results of phylogenetic, phenotypic, chemotaxonomic, and genotypic analyses suggest that strain HDW12AT represents a novel species within the genus Nocardioides, and strains HDW15BT and HDW15CT represent two novel species within the genus Sphingomonas. We propose the names Nocardioides piscis for strain HDW12AT (= KACC 21336T = KCTC 49321T = JCM 33670T), Sphingomonas piscis for strain HDW15BT (= KACC 21341T = KCTC 72588T = JCM 33738T), and Sphingomonas sinipercae for strain HDW15CT (= KACC 21342T = KCTC 72589T = JCM 33739T).

Improved protein structure prediction by deep learning irrespective of co-evolution information

Predicting the tertiary structure of a protein from its primary sequence has been greatly improved by integrating deep learning and co-evolutionary analysis, as shown in CASP13 and CASP14. We describe our latest study of this idea, analyzing the efficacy of network size and co-evolution data and its performance on both natural and designed proteins. We show that a large ResNet (convolutional residual neural networks) can predict structures of correct folds for 26 out of 32 CASP13 free-modeling (FM) targets and L/5 long-range contacts with precision over 80%. When co-evolution is not used ResNet still can predict structures of correct folds for 18 CASP13 FM targets, greatly exceeding previous methods that do not use co-evolution either. Even with only primary sequence ResNet can predict structures of correct folds for all tested human-designed proteins. In addition, ResNet may fare better for the designed proteins when trained without co-evolution than with co-evolution. These results suggest that ResNet does not simply denoise co-evolution signals, but instead may learn important protein sequence-structure relationship. This has important implications on protein design and engineering especially when co-evolutionary data is unavailable.

DOE JGI Metagenome Workflow

The DOE Joint Genome Institute (JGI) Metagenome Workflow performs metagenome data processing, including assembly; structural, functional, and taxonomic annotation; and binning of metagenomic data sets that are subsequently included into the Integrated Microbial Genomes and Microbiomes (IMG/M) (I.-M. A. Chen, K. Chu, K. Palaniappan, A. Ratner, et al., Nucleic Acids Res, 49:D751–D763, 2021, https://doi.org/10.1093/nar/gkaa939) comparative analysis system and provided for download via the JGI data portal (https://genome.jgi.doe.gov/portal/). This workflow scales to run on thousands of metagenome samples per year, which can vary by the complexity of microbial communities and sequencing depth. Here, we describe the different tools, databases, and parameters used at different steps of the workflow to help with the interpretation of metagenome data available in IMG and to enable researchers to apply this workflow to their own data. We use 20 publicly available sediment metagenomes to illustrate the computing requirements for the different steps and highlight the typical results of data processing. The workflow modules for read filtering and metagenome assembly are available as a workflow description language (WDL) file (https://code.jgi.doe.gov/BFoster/jgi_meta_wdl). The workflow modules for annotation and binning are provided as a service to the user community at https://img.jgi.doe.gov/submit and require filling out the project and associated metadata descriptions in the Genomes OnLine Database (GOLD) (S. Mukherjee, D. Stamatis, J. Bertsch, G. Ovchinnikova, et al., Nucleic Acids Res, 49:D723–D733, 2021, https://doi.org/10.1093/nar/gkaa983).

IMPORTANCE The DOE JGI Metagenome Workflow is designed for processing metagenomic data sets starting from Illumina fastq files. It performs data preprocessing, error correction, assembly, structural and functional annotation, and binning. The results of processing are provided in several standard formats, such as fasta and gff, and can be used for subsequent integration into the Integrated Microbial Genomes and Microbiomes (IMG/M) system where they can be compared to a comprehensive set of publicly available metagenomes. As of 30 July 2020, 7,155 JGI metagenomes have been processed by the DOE JGI Metagenome Workflow. Here, we present a metagenome workflow developed at the JGI that generates rich data in standard formats and has been optimized for downstream analyses ranging from assessment of the functional and taxonomic composition of microbial communities to genome-resolved metagenomics and the identification and characterization of novel taxa. This workflow is currently being used to analyze thousands of metagenomic data sets in a consistent and standardized manner.

Genomics, Exometabolomics, and Metabolic Probing Reveal Conserved Proteolytic Metabolism of Thermoflexus hugenholtzii and Three Candidate Species From China and Japan

Thermoflexus hugenholtzii JAD2^T, the only cultured representative of the Chloroflexota order Thermoflexales, is abundant in Great Boiling Spring (GBS), NV, United States, and close relatives inhabit geothermal systems globally. However, no defined medium exists for T. hugenholtzii JAD2^T and no single carbon source is known to support its growth, leaving key knowledge gaps in its metabolism and nutritional needs. Here, we report comparative genomic analysis of the draft genome of T. hugenholtzii JAD2^T and eight closely related metagenome-assembled genomes (MAGs) from geothermal sites in China, Japan, and the United States, representing “Candidatus Thermoflexus japonica,” “Candidatus Thermoflexus tengchongensis,” and “Candidatus Thermoflexus sinensis.” Genomics was integrated with targeted exometabolomics and ¹³C metabolic probing of T. hugenholtzii. The Thermoflexus genomes each code for complete central carbon metabolic pathways and an unusually high abundance and diversity of peptidases, particularly Metallo- and Serine peptidase families, along with ABC transporters for peptides and some amino acids. The T. hugenholtzii JAD2^T exometabolome provided evidence of extracellular proteolytic activity based on the accumulation of free amino acids. However, several neutral and polar amino acids appear not to be utilized, based on their accumulation in the medium and the lack of annotated transporters. Adenine and adenosine were scavenged, and thymine and nicotinic acid were released, suggesting interdependency with other organisms in situ. Metabolic probing of T. hugenholtzii JAD2^T using ¹³C-labeled compounds provided evidence of oxidation of glucose, pyruvate, cysteine, and citrate, and functioning glycolytic, tricarboxylic acid (TCA), and oxidative pentose-phosphate pathways (PPPs). However, differential use of position-specific ¹³C-labeled compounds showed that glycolysis and the TCA cycle were uncoupled. Thus, despite the high abundance of Thermoflexus in sediments of some geothermal systems, they appear to be highly focused on chemoorganotrophy, particularly protein degradation, and may interact extensively with other microorganisms in situ.

Bioinformatic Prediction of an tRNASec Gene Nested inside an Elongation Factor SelB Gene in Alphaproteobacteria

In bacteria, selenocysteine (Sec) is incorporated into proteins via the recoding of a particular codon, the UGA stop codon in most cases. Sec-tRNASec is delivered to the ribosome by the Sec-dedicated elongation factor SelB that also recognizes a Sec-insertion sequence element following the codon on the mRNA. Since the excess of SelB may lead to sequestration of Sec-tRNASec under selenium deficiency or oxidative stress, the expression levels of SelB and tRNASec should be regulated. In this bioinformatic study, I analyzed the Rhizobiales SelB species because they were annotated to have a non-canonical C-terminal extension. I found that the open reading frame (ORF) of diverse Alphaproteobacteria selB genes includes an entire tRNASec sequence (selC) and overlaps with the start codon of the downstream ORF. A remnant tRNASec sequence was found in the Sinorhizobium melilotiselB genes whose products have a shorter C-terminal extension. Similar overlapping traits were found in Gammaproteobacteria and Nitrospirae. I hypothesized that once the tRNASec moiety is folded and processed, the expression of the full-length SelB may be repressed. This is the first report on a nested tRNA gene inside a protein ORF in bacteria.

Unexpected distribution of the 4-formylaminooxyvinylglycine (FVG) biosynthetic pathway in Pseudomonas and beyond

The biological herbicide and antibiotic 4-formylaminooxyvinylglycine (FVG) was originally isolated from several rhizosphere-associated strains of Pseudomonas fluorescens. Biosynthesis of FVG is dependent on the gvg biosynthetic gene cluster in P. fluorescens. In this investigation, we used comparative genomics to identify strains with the genetic potential to produce FVG due to presence of a gvg gene cluster. These strains primarily belong to two groups of Pseudomonas, P. fluorescens and P. syringae, however, a few strains with the gvg cluster were found outside of Pseudomonas. Mass spectrometry confirmed that all tested strains of the P. fluorescens species group produced FVG. However, P. syringae strains did not produce FVG under standard conditions. Several lines of evidence regarding the transmission of the gvg cluster including a robust phylogenetic analysis suggest that it was introduced multiple times through horizontal gene transfer within the Pseudomonas lineage as well as in select lineages of Thiomonas, Burkholderia and Pantoea. Together, these data broaden our understanding of the evolution and diversity of FVG biosynthesis. In the course of this investigation, additional gene clusters containing only a subset of the genes required to produce FVG were identified in a broad range of bacteria, including many non-pseudomonads.

High Potential for Biomass-Degrading Enzymes Revealed by Hot Spring Metagenomics

Enzyme stability and activity at elevated temperatures are important aspects in biotechnological industries, such as the conversion of plant biomass into biofuels. In order to reduce the costs and increase the efficiency of biomass conversion, better enzymatic processing must be developed. Hot springs represent a treasure trove of underexplored microbiological and protein chemistry diversity. Herein, we conduct an exploratory study into the diversity of hot spring biomass-degrading potential. We describe the taxonomic diversity and carbohydrate active enzyme (CAZyme) coding potential in 71 publicly available metagenomic datasets from 58 globally distributed terrestrial geothermal features. Through taxonomic profiling, we detected a wide diversity of microbes unique to varying temperature and pH ranges. Biomass-degrading enzyme potential included all five classes of CAZymes and we described the presence or absence of genes encoding 19 glycosyl hydrolases hypothesized to be involved with cellulose, hemicellulose, and oligosaccharide degradation. Our results highlight hot springs as a promising system for the further discovery and development of thermo-stable biomass-degrading enzymes that can be applied toward generation of renewable biofuels. This study lays a foundation for future research to further investigate the functional diversity of hot spring biomass-degrading enzymes and their potential utility in biotechnological processing.

The 2021 Nucleic Acids Research database issue and the online molecular biology database collection

The 2021 Nucleic Acids Research database Issue contains 189 papers spanning a wide range of biological fields and investigation. It includes 89 papers reporting on new databases and 90 covering recent changes to resources previously published in the Issue. A further ten are updates on databases most recently published elsewhere. Seven new databases focus on COVID-19 and SARS-CoV-2 and many others offer resources for studying the virus. Major returning nucleic acid databases include NONCODE, Rfam and RNAcentral. Protein family and domain databases include COG, Pfam, SMART and Panther. Protein structures are covered by RCSB PDB and dispersed proteins by PED and MobiDB. In metabolism and signalling, STRING, KEGG and WikiPathways are featured, along with returning KLIFS and new DKK and KinaseMD, all focused on kinases. IMG/M and IMG/VR update in the microbial and viral genome resources section, while human and model organism genomics resources include Flybase, Ensembl and UCSC Genome Browser. Cancer studies are covered by updates from canSAR and PINA, as well as newcomers CNCdatabase and Oncovar for cancer drivers. Plant comparative genomics is catered for by updates from Gramene and GreenPhylDB. The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). The NAR online Molecular Biology Database Collection has been substantially updated, revisiting nearly 1000 entries, adding 90 new resources and eliminating 86 obsolete databases, bringing the current total to 1641 databases. It is available at https://www.oxfordjournals.org/nar/database/c/.