Next-Generation High-Throughput Functional Annotation of Microbial Genomes

High-Content Screening and Computational Prediction Reveal Viral Genes That Suppress the Innate Immune Response

Suppression of the host innate immune response is a critical aspect of viral replication. Upon infection, viruses may introduce one or more proteins that inhibit key immune pathways, such as the type I interferon pathway. However, the ability to predict and evaluate viral protein bioactivity on targeted pathways remains challenging and is typically done on a single-virus or -gene basis. Here, we present a medium-throughput high-content cell-based assay to reveal the immunosuppressive effects of viral proteins. To test the predictive power of our approach, we developed a library of 800 genes encoding known, predicted, and uncharacterized human virus genes. We found that previously known immune suppressors from numerous viral families such as Picornaviridae and Flaviviridae recorded positive responses. These include a number of viral proteases for which we further confirmed that innate immune suppression depends on protease activity. A class of predicted inhibitors encoded by Rhabdoviridae viruses was demonstrated to block nuclear transport, and several previously uncharacterized proteins from uncultivated viruses were shown to inhibit nuclear transport of the transcription factors NF-κB and interferon regulatory factor 3 (IRF3). We propose that this pathway-based assay, together with early sequencing, gene synthesis, and viral infection studies, could partly serve as the basis for rapid in vitro characterization of novel viral proteins. IMPORTANCE Infectious diseases caused by viral pathogens exacerbate health care and economic burdens. Numerous viral biomolecules suppress the human innate immune system, enabling viruses to evade an immune response from the host. Despite our current understanding of viral replications and immune evasion, new viral proteins, including those encoded by uncultivated viruses or emerging viruses, are being unearthed at a rapid pace from large-scale sequencing and surveillance projects. The use of medium- and high-throughput functional assays to characterize immunosuppressive functions of viral proteins can advance our understanding of viral replication and possibly treatment of infections. In this study, we assembled a large viral-gene library from diverse viral families and developed a high-content assay to test for inhibition of innate immunity pathways. Our work expands the tools that can rapidly link sequence and protein function, representing a practical step toward early-stage evaluation of emerging and understudied viruses.

Approaches to Unmask Functioning of the Uncultured Microbial Majority From Extreme Habitats on the Seafloor

Researchers have recognized the potential of enzymes and metabolic pathways hidden among the unseen majority of Earth’s microorganisms for decades now. Most of the microbes expected to colonize the seafloor and its subsurface are currently uncultured. Thus, their ability and contribution to element cycling remain enigmatic. Given that the seafloor covers ∼70% of our planet, this amounts to an uncalled potential of unrecognized metabolic properties and interconnections catalyzed by this microbial dark matter. Consequently, a tremendous black box awaits discovery of novel enzymes, catalytic abilities, and metabolic properties in one of the largest habitats on Earth. This mini review summarizes the current knowledge of cultivation-dependent and -independent techniques applied to seafloor habitats to unravel the role of the microbial dark matter. It highlights the great potential that combining microbiological and biogeochemical data from in situ experiments with molecular tools has for providing a holistic understanding of bio-geo-coupling in seafloor habitats and uses hydrothermal vent systems as a case example.

Construction of Aerobic/Anaerobic-Substrate-Induced Gene Expression Procedure for Exploration of Metagenomes From Subseafloor Sediments

Substrate-induced gene expression (SIGEX) is a high-throughput promoter-trap method. It is a function-based metagenomic screening tool that relies on transcriptional activation of a reporter gene green fluorescence protein (gfp) by a metagenomic DNA library upon induction with a substrate. However, its use is limited because of the relatively small size of metagenomic DNA libraries and incompatibility with screening metagenomes from anaerobic environments. In this study, these limitations of SIGEX were addressed by fine-tuning metagenome DNA library construction protocol and by using Evoglow, a green fluorescent protein that forms a chromophore even under anaerobic conditions. Two metagenomic libraries were constructed for subseafloor sediments offshore Shimokita Peninsula (Pacific Ocean) and offshore Joetsu (Japan Sea). The library construction protocol was improved by (a) eliminating short DNA fragments, (b) applying topoisomerase-based high-efficiency ligation, (c) optimizing insert DNA concentration, and (d) column-based DNA enrichment. This led to a successful construction of metagenome DNA libraries of approximately 6 Gbp for both samples. SIGEX screening using five aromatic compounds (benzoate, 3-chlorobenzoate, 3-hydroxybenzoate, phenol, and 2,4-dichlorophenol) under aerobic and anaerobic conditions revealed significant differences in the inducible clone ratios under these conditions. 3-Chlorobenzoate and 2,4-dichlorophenol led to a higher induction ratio than that for the other non-chlorinated aromatic compounds under both aerobic and anaerobic conditions. After the further screening of induced clones, a clone induced by 3-chlorobenzoate only under anaerobic conditions was isolated and characterized. The clone harbors a DNA insert that encodes putative open reading frames of unknown function. Previous aerobic SIGEX attempts succeeded in the isolation of gene fragments from anaerobes. This study demonstrated that some gene fragments require a strict in vivo reducing environment to function and may be potentially missed when screened by aerobic induction. The newly developed anaerobic SIGEX scheme will facilitate functional exploration of metagenomes from the anaerobic biosphere.

Prokaryotic Genome Annotation

In the last decade, the high-throughput and relatively low cost of short-read sequencing technologies have revolutionized prokaryotic genomics. This has led to an exponential increase in the number of bacterial and archaeal genome sequences available, as well as corresponding increase of genome assembly and annotation tools developed. Together, these hardware and software technologies have given scientists unprecedented options to study their chosen microbial systems without the need for large teams of bioinformaticists or supercomputing facilities. While these analysis tools largely fall into only a few categories, each may have different requirements, caveats and file formats, and some may be rarely updated or even abandoned. And so, despite the apparent ease in sequencing and analyzing a prokaryotic genome, it is no wonder that the budding genomicist may quickly find oneself overwhelmed. Here, we aim to provide the reader with an overview of genome annotation and its most important considerations, as well as an easy-to-follow protocol to get started with annotating a prokaryotic genome.

Manual Annotation Studio (MAS): a collaborative platform for manual functional annotation of viral and microbial genomes

Background

Functional genome annotation is the process of labelling functional genomic regions with descriptive information. Manual curation can produce higher quality genome annotations than fully automated methods. Manual annotation efforts are time-consuming and complex; however, software can help reduce these drawbacks.

Results

We created Manual Annotation Studio (MAS) to improve the efficiency of the process of manual functional annotation prokaryotic and viral genomes. MAS allows users to upload unannotated genomes, provides an interface to edit and upload annotations, tracks annotation history and progress, and saves data to a relational database. MAS provides users with pertinent information through a simple point and click interface to execute and visualize results for multiple homology search tools (blastp, rpsblast, and HHsearch) against multiple databases (Swiss-Prot, nr, CDD, PDB, and an internally generated database). MAS was designed to accept connections over the local area network (LAN) of a lab or organization so multiple users can access it simultaneously. MAS can take advantage of high-performance computing (HPC) clusters by interfacing with SGE or SLURM and data can be exported from MAS in a variety of formats (FASTA, GenBank, GFF, and excel).

Conclusions

MAS streamlines and provides structure to manual functional annotation projects. MAS enhances the ability of users to generate, interpret, and compare results from multiple tools. The structure that MAS provides can improve project organization and reduce annotation errors. MAS is ideal for team-based annotation projects because it facilitates collaboration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08029-8.

RNA atlas of human bacterial pathogens uncovers stress dynamics linked to infection

Bacterial processes necessary for adaption to stressful host environments are potential targets for new antimicrobials. Here, we report large-scale transcriptomic analyses of 32 human bacterial pathogens grown under 11 stress conditions mimicking human host environments. The potential relevance of the in vitro stress conditions and responses is supported by comparisons with available in vivo transcriptomes of clinically important pathogens. Calculation of a probability score enables comparative cross-microbial analyses of the stress responses, revealing common and unique regulatory responses to different stresses, as well as overlapping processes participating in different stress responses. We identify conserved and species-specific 'universal stress responders', that is, genes showing altered expression in multiple stress conditions. Non-coding RNAs are involved in a substantial proportion of the responses. The data are collected in a freely available, interactive online resource (PATHOgenex).

A Mycobacterial Systems Resource for the Research Community

Functional characterization of bacterial proteins lags far behind the identification of new protein families. This is especially true for bacterial species that are more difficult to grow and genetically manipulate than model systems such as Escherichia coli and Bacillus subtilis To facilitate functional characterization of mycobacterial proteins, we have established a Mycobacterial Systems Resource (MSR) using the model organism Mycobacterium smegmatis This resource focuses specifically on 1,153 highly conserved core genes that are common to many mycobacterial species, including Mycobacterium tuberculosis, in order to provide the most relevant information and resources for the mycobacterial research community. The MSR includes both biological and bioinformatic resources. The biological resource includes (i) an expression plasmid library of 1,116 genes fused to a fluorescent protein for determining protein localization; (ii) a library of 569 precise deletions of nonessential genes; and (iii) a set of 843 CRISPR-interference (CRISPRi) plasmids specifically targeted to silence expression of essential core genes and genes for which a precise deletion was not obtained. The bioinformatic resource includes information about individual genes and a detailed assessment of protein localization. We anticipate that integration of these initial functional analyses and the availability of the biological resource will facilitate studies of these core proteins in many Mycobacterium species, including the less experimentally tractable pathogens M. abscessus, M. avium, M. kansasii, M. leprae, M. marinum, M. tuberculosis, and M. ulceransIMPORTANCE Diseases caused by mycobacterial species result in millions of deaths per year globally, and present a substantial health and economic burden, especially in immunocompromised patients. Difficulties inherent in working with mycobacterial pathogens have hampered the development and application of high-throughput genetics that can inform genome annotations and subsequent functional assays. To facilitate mycobacterial research, we have created a biological and bioinformatic resource (https://msrdb.org/) using Mycobacterium smegmatis as a model organism. The resource focuses specifically on 1,153 proteins that are highly conserved across the mycobacterial genus and, therefore, likely perform conserved mycobacterial core functions. Thus, functional insights from the MSR will apply to all mycobacterial species. We believe that the availability of this mycobacterial systems resource will accelerate research throughout the mycobacterial research community.

Novel modulators of p53-signaling encoded by unknown genes of emerging viruses

The p53 transcription factor plays a key role both in cancer and in the cell-intrinsic response to infections. The ORFEOME project hypothesized that novel p53-virus interactions reside in hitherto uncharacterized, unknown, or hypothetical open reading frames (orfs) of human viruses. Hence, 172 orfs of unknown function from the emerging viruses SARS-Coronavirus, MERS-Coronavirus, influenza, Ebola, Zika (ZIKV), Chikungunya and Kaposi Sarcoma-associated herpesvirus (KSHV) were de novo synthesized, validated and tested in a functional screen of p53 signaling. This screen revealed novel mechanisms of p53 virus interactions and two viral proteins KSHV orf10 and ZIKV NS2A binding to p53. Originally identified as the target of small DNA tumor viruses, these experiments reinforce the notion that all viruses, including RNA viruses, interfere with p53 functions. These results validate this resource for analogous systems biology approaches to identify functional properties of uncharacterized viral proteins, long non-coding RNAs and micro RNAs.

New viruses are constantly emerging. The ORFEOME project was based on the hypothesis that every virus, regardless of its molecular makeup and biology should encode functions that intersect the p53 signaling network, since p53 guards the cell from genomic insults, of which depositing a foreign, viral nucleic acid is one. The result of the ORFEOME screen of proteins without any known function, of predicted open reading frames and of suspected non-coding RNAs is the identification of two viral proteins that interact with p53. The first one, orf10, is encoded by Kaposi Sarcoma-associated herpesvirus and the second one, NS2A, is encoded by the Zika virus.

Microbiome definition re-visited: old concepts and new challenges

The field of microbiome research has evolved rapidly over the past few decades and has become a topic of great scientific and public interest. As a result of this rapid growth in interest covering different fields, we are lacking a clear commonly agreed definition of the term "microbiome." Moreover, a consensus on best practices in microbiome research is missing. Recently, a panel of international experts discussed the current gaps in the frame of the European-funded MicrobiomeSupport project. The meeting brought together about 40 leaders from diverse microbiome areas, while more than a hundred experts from all over the world took part in an online survey accompanying the workshop. This article excerpts the outcomes of the workshop and the corresponding online survey embedded in a short historical introduction and future outlook. We propose a definition of microbiome based on the compact, clear, and comprehensive description of the term provided by Whipps et al. in 1988, amended with a set of novel recommendations considering the latest technological developments and research findings. We clearly separate the terms microbiome and microbiota and provide a comprehensive discussion considering the composition of microbiota, the heterogeneity and dynamics of microbiomes in time and space, the stability and resilience of microbial networks, the definition of core microbiomes, and functionally relevant keystone species as well as co-evolutionary principles of microbe-host and inter-species interactions within the microbiome. These broad definitions together with the suggested unifying concepts will help to improve standardization of microbiome studies in the future, and could be the starting point for an integrated assessment of data resulting in a more rapid transfer of knowledge from basic science into practice. Furthermore, microbiome standards are important for solving new challenges associated with anthropogenic-driven changes in the field of planetary health, for which the understanding of microbiomes might play a key role. Video Abstract.

Metal-ion-induced expression of gene fragments from subseafloor micro-organisms in the Kumano forearc basin, Nankai Trough

AIMS

Using substrate-induced gene-expression (SIGEX) screening on subseafloor sediment samples from the Nankai Trough, Japan, we identified gene fragments showing an induction response to metal ions.

METHODS AND RESULTS

Environmental DNA libraries in Escherichia coli host cells were tested by the addition of metal ions (Ni²⁺ , Co²⁺ , Ga³⁺ or Mo⁶⁺ ), followed by cell sorting of clones exhibiting green fluorescence upon co-expression of green fluorescence protein downstream of the inserted gene fragments. One clone displayed Ni²⁺ -specific induction, three clones displayed Ga³⁺ -specific induction and three clones displayed an induction response to multiple metal ions. DNA sequence analysis showed that a variety of genes showed induction responses in the screened clones.

CONCLUSIONS

Using the SIGEX approach, we retrieved gene fragments with no previously identified response to metal ions that exhibited metal-ion-induced expression. This method has the potential to promote exploration of gene function through gene-induction response.

SIGNIFICANCE AND IMPACT OF THE STUDY

We successfully linked gene-induction response with sequence information for gene fragments of previously unknown function. The SIGEX-based approach exhibited the potential to identify genetic function in unknown gene pools from the deep subseafloor biosphere, as well as novel genetic components for future biotechnological applications.

Molecular control of gene expression by Brucella BaaR, an IclR-type transcriptional repressor

The general stress response sigma factor σE1 directly and indirectly regulates the transcription of dozens of genes that influence stress survival and host infection in the zoonotic pathogen Brucella abortus Characterizing the functions of σE1-regulated genes therefore would contribute to our understanding of B. abortus physiology and infection biology. σE1 indirectly activates transcription of the IclR family regulator Bab2_0215, but the function of this regulator remains undefined. Here, we present a structural and functional characterization of Bab2_0215, which we have named B rucella adipic acid-activated regulator (BaaR). We found that BaaR adopts a classic IclR-family fold and directly represses the transcription of two operons with predicted roles in carboxylic acid oxidation. BaaR binds two sites on chromosome II between baaR and a divergently transcribed hydratase/dehydrogenase (acaD2), and it represses transcription of both genes. We identified three carboxylic acids (adipic acid, tetradecanedioic acid, and ϵ-aminocaproic acid) and a lactone (ϵ-caprolactone) that enhance transcription from the baaR and acaD2 promoters. However, neither the activating acids nor caprolactone enhanced transcription by binding directly to BaaR. Induction of baaR transcription by adipic acid required the gene bab2_0213, which encodes a major facilitator superfamily transporter, suggesting that Bab2_0213 transports adipic acid across the inner membrane. We conclude that a suite of structurally related organic molecules activate transcription of genes repressed by BaaR. Our study provides molecular-level understanding of a gene expression program in B. abortus that is downstream of σE1.

Genome-wide discovery of novel M1T1 group A streptococcal determinants important for fitness and virulence during soft-tissue infection

The Group A Streptococcus remains a significant human pathogen causing a wide array of disease ranging from self-limiting to life-threatening invasive infections. Epithelium (skin or throat) colonization with progression to the subepithelial tissues is the common step in all GAS infections. Here, we used transposon-sequencing (Tn-seq) to define the GAS 5448 genetic requirements for in vivo fitness in subepithelial tissue. A near-saturation transposon library of the M1T1 GAS 5448 strain was injected subcutaneously into mice, producing suppurative inflammation at 24 h that progressed to prominent abscesses with tissue necrosis at 48 h. The library composition was monitored en masse by Tn-seq and ratios of mutant abundance comparing the output (12, 24 and 48 h) versus input (T₀) mutant pools were calculated for each gene. We identified a total of 273 subcutaneous fitness (scf) genes with 147 genes (55 of unknown function) critical for the M1T1 GAS 5448 fitness in vivo; and 126 genes (53 of unknown function) potentially linked to in vivo fitness advantage. Selected scf genes were validated in competitive subcutaneous infection with parental 5448. Two uncharacterized genes, scfA and scfB, encoding putative membrane-associated proteins and conserved among Gram-positive pathogens, were further characterized. Defined scfAB mutants in GAS were outcompeted by wild type 5448 in vivo, attenuated for lesion formation in the soft tissue infection model and dissemination to the bloodstream. We hypothesize that scfAB play an integral role in enhancing adaptation and fitness of GAS during localized skin infection, and potentially in propagation to other deeper host environments.

The WHO ranks the Group A Streptococcus (GAS) in the top 10 leading causes of morbidity and mortality from infectious diseases worldwide. GAS is a strict human pathogen causing both benign superficial infections as well as life-threatening invasive diseases. All GAS infections begin by colonization of an epithelium (throat or skin) followed by propagation into subepithelial tissues. The genetic requirements for M1T1 GAS 5448 within this niche were interrogated by in vivo transposon sequencing (Tn-seq), identifying 273 subcutaneous fitness (scf) genes with 108 of those previously of “unknown function”. Two yet uncharacterized genes, scfA and scfB, were shown to be critical during GAS 5448 soft tissue infection and dissemination into the bloodstream. Thus, this study improves the functional annotation of the GAS genome, providing new insights into GAS pathophysiology and enhancing the development of novel GAS therapeutics.

MERS-CoV Accessory ORFs Play Key Role for Infection and Pathogenesis

While dispensable for viral replication, coronavirus (CoV) accessory open reading frame (ORF) proteins often play critical roles during infection and pathogenesis. Utilizing a previously generated mutant, we demonstrate that the absence of all four Middle East respiratory syndrome CoV (MERS-CoV) accessory ORFs (deletion of ORF3, -4a, -4b, and -5 [dORF3-5]) has major implications for viral replication and pathogenesis. Importantly, attenuation of the dORF3-5 mutant is primarily driven by dysregulated host responses, including disrupted cell processes, augmented interferon (IFN) pathway activation, and robust inflammation. In vitro replication attenuation also extends to in vivo models, allowing use of dORF3-5 as a live attenuated vaccine platform. Finally, examination of ORF5 implicates a partial role in modulation of NF-κB-mediated inflammation. Together, the results demonstrate the importance of MERS-CoV accessory ORFs for pathogenesis and highlight them as potential targets for surveillance and therapeutic treatments moving forward.

The initial emergence and periodic outbreaks of MERS-CoV highlight a continuing threat posed by zoonotic pathogens to global public health. In these studies, mutant virus generation demonstrates the necessity of accessory ORFs in regard to MERS-CoV infection and pathogenesis. With this in mind, accessory ORF functions can be targeted for both therapeutic and vaccine treatments in response to MERS-CoV and related group 2C coronaviruses. In addition, disruption of accessory ORFs in parallel may offer a rapid response platform to attenuation of future emergent strains based on both SARS- and MERS-CoV accessory ORF mutants.