Functional and evolutionary significance of unknown genes from uncultivated taxa.

Many of the Earth's microbes remain uncultured and understudied, limiting our understanding of the functional and evolutionary aspects of their genetic material, which remain largely overlooked in most metagenomic studies1. Here we analysed 149,842 environmental genomes from multiple habitats2-6 and compiled a curated catalogue of 404,085 functionally and evolutionarily significant novel (FESNov) gene families exclusive to uncultivated prokaryotic taxa. All FESNov families span multiple species, exhibit strong signals of purifying selection and qualify as new orthologous groups, thus nearly tripling the number of bacterial and archaeal gene families described to date. The FESNov catalogue is enriched in clade-specific traits, including 1,034 novel families that can distinguish entire uncultivated phyla, classes and orders, probably representing synapomorphies that facilitated their evolutionary divergence. Using genomic context analysis and structural alignments we predicted functional associations for 32.4% of FESNov families, including 4,349 high-confidence associations with important biological processes. These predictions provide a valuable hypothesis-driven framework that we used for experimental validatation of a new gene family involved in cell motility and a novel set of antimicrobial peptides. We also demonstrate that the relative abundance profiles of novel families can discriminate between environments and clinical conditions, leading to the discovery of potentially new biomarkers associated with colorectal cancer. We expect this work to enhance future metagenomics studies and expand our knowledge of the genetic repertory of uncultivated organisms.

Towards the biogeography of prokaryotic genes.

Microbial genes encode the majority of the functional repertoire of life on earth. However, despite increasing efforts in metagenomic sequencing of various habitats1-3, little is known about the distribution of genes across the global biosphere, with implications for human and planetary health. Here we constructed a non-redundant gene catalogue of 303 million species-level genes (clustered at 95% nucleotide identity) from 13,174 publicly available metagenomes across 14 major habitats and use it to show that most genes are specific to a single habitat. The small fraction of genes found in multiple habitats is enriched in antibiotic-resistance genes and markers for mobile genetic elements. By further clustering these species-level genes into 32 million protein families, we observed that a small fraction of these families contain the majority of the genes (0.6% of families account for 50% of the genes). The majority of species-level genes and protein families are rare. Furthermore, species-level genes, and in particular the rare ones, show low rates of positive (adaptive) selection, supporting a model in which most genetic variability observed within each protein family is neutral or nearly neutral.

eggNOG 6.0: enabling comparative genomics across 12 535 organisms.

The eggNOG (evolutionary gene genealogy Non-supervised Orthologous Groups) database is a bioinformatics resource providing orthology data and comprehensive functional information for organisms from all domains of life. Here, we present a major update of the database and website (version 6.0), which increases the number of covered organisms to 12 535 reference species, expands functional annotations, and implements new functionality. In total, eggNOG 6.0 provides a hierarchy of over 17M orthologous groups (OGs) computed at 1601 taxonomic levels, spanning 10 756 bacterial, 457 archaeal and 1322 eukaryotic organisms. OGs have been thoroughly annotated using recent knowledge from functional databases, including KEGG, Gene Ontology, UniProtKB, BiGG, CAZy, CARD, PFAM and SMART. eggNOG also offers phylogenetic trees for all OGs, maximising utility and versatility for end users while allowing researchers to investigate the evolutionary history of speciation and duplication events as well as the phylogenetic distribution of functional terms within each OG. Furthermore, the eggNOG 6.0 website contains new functionality to mine orthology and functional data with ease, including the possibility of generating phylogenetic profiles for multiple OGs across species or identifying single-copy OGs at custom taxonomic levels. eggNOG 6.0 is available at http://eggnog6.embl.de.

GeCoViz: genomic context visualisation of prokaryotic genes from a functional and evolutionary perspective.

Synteny conservation analysis is a well-established methodology to investigate the potential functional role of unknown prokaryotic genes. However, bioinformatic tools to reconstruct and visualise genomic contexts usually depend on slow computations, are restricted to narrow taxonomic ranges, and/or do not allow for the functional and interactive exploration of neighbouring genes across different species. Here, we present GeCoViz, an online resource built upon 12 221 reference prokaryotic genomes that provides fast and interactive visualisation of custom genomic regions anchored by any target gene, which can be sought by either name, orthologous group (KEGGs, eggNOGs), protein domain (PFAM) or sequence. To facilitate functional and evolutionary interpretation, GeCoViz allows to customise the taxonomic scope of each analysis and provides comprehensive annotations of the neighbouring genes. Interactive visualisation options include, among others, the scaled representations of gene lengths and genomic distances, and on the fly calculation of synteny conservation of neighbouring genes, which can be highlighted based on custom thresholds. The resulting plots can be downloaded as high-quality images for publishing purposes. Overall, GeCoViz offers an easy-to-use, comprehensive, fast and interactive web-based tool for investigating the genomic context of prokaryotic genes, and is freely available at https://gecoviz.cgmlab.org.

A protease-mediated mechanism regulates the cytochrome c6/plastocyanin switch in Synechocystis sp. PCC 6803.

After the Great Oxidation Event (GOE), iron availability was greatly decreased, and photosynthetic organisms evolved several alternative proteins and mechanisms. One of these proteins, plastocyanin, is a type I blue-copper protein that can replace cytochrome c6 as a soluble electron carrier between cytochrome b6f and photosystem I. In most cyanobacteria, expression of these two alternative proteins is regulated by copper availability, but the regulatory system remains unknown. Herein, we provide evidence that the regulatory system is composed of a BlaI/CopY-family transcription factor (PetR) and a BlaR-membrane protease (PetP). PetR represses petE (plastocyanin) expression and activates petJ (cytochrome c6), while PetP controls PetR levels in vivo. Using whole-cell extracts, we demonstrated that PetR degradation requires both PetP and copper. Transcriptomic analysis revealed that the PetRP system regulates only four genes (petE, petJ, slr0601, and slr0602), highlighting its specificity. Furthermore, the presence of petE and petRP in early branching cyanobacteria indicates that acquisition of these genes could represent an early adaptation to decreased iron bioavailability following the GOE.

Structure-based analyses of Salmonella RcsB variants unravel new features of the Rcs regulon.

RcsB is a transcriptional regulator that controls expression of numerous genes in enteric bacteria. RcsB accomplishes this role alone or in combination with auxiliary transcriptional factors independently or dependently of phosphorylation. To understand the mechanisms by which RcsB regulates such large number of genes, we performed structural studies as well as in vitro and in vivo functional studies with different RcsB variants. Our structural data reveal that RcsB binds promoters of target genes such as rprA and flhDC in a dimeric active conformation. In this state, the RcsB homodimer docks the DNA-binding domains into the major groove of the DNA, facilitating an initial weak read-out of the target sequence. Interestingly, comparative structural analyses also show that DNA binding may stabilize an active conformation in unphosphorylated RcsB. Furthermore, RNAseq performed in strains expressing wild-type or several RcsB variants provided new insights into the contribution of phosphorylation to gene regulation and assign a potential role of RcsB in controlling iron metabolism. Finally, we delimited the RcsB box for homodimeric active binding to DNA as the sequence TN(G/A)GAN4TC(T/C)NA. This RcsB box was found in promoter, intergenic and intragenic regions, facilitating both increased or decreased gene transcription.

Correction to: Genome analysis of Salmonella enterica subsp. diarizonae isolates from invasive human infections reveals enrichment of virulence-related functions in lineage ST1256.

An amendment to this paper has been published and can be accessed via the original article.

Genome analysis of Salmonella enterica subsp. diarizonae isolates from invasive human infections reveals enrichment of virulence-related functions in lineage ST1256.

BACKGROUND: Salmonella enterica subsp. diarizonae (IIIb) is frequently isolated from the environment, cold-blooded reptiles, sheep and humans; however only a few studies describe the isolation of this subspecies from invasive human infections. The factors contributing to this unusual behavior are currently unknown. RESULTS: We report here the genome features of two diarizonae strains, SBO13 and SBO27, isolated from endocervical tissue collected post-abortion and from cerebrospinal fluid of a newborn child, respectively, in the city of Santa Cruz, Bolivia. Although isolated six years apart, SBO27 in 2008 and SBO13 in 2014, both strains belong to the same sequence type 1256 (ST1256) and show a high degree of genome conservation sharing more than 99% of their genes, including the conservation of a ~ 10 kb plasmid. A prominent feature of the two genomes is the presence of 24 genomic islands (GIs), in addition to 10 complete Salmonella pathogenicity islands (SPI) and fragments of SPI-7, a pathogenicity island first reported in the human-adapted serovar Typhi. Some of the GIs identified in SBO13 and SBO27 harbor genes putatively encoding auto-transporters involved in adhesion, lipopolysaccharide modifying enzymes, putative toxins, pili-related proteins, efflux pumps, and several putative membrane cation transport related-genes, among others. These two Bolivian isolates also share genes encoding the type-III secretion system effector proteins SseK2, SseK3 and SlrP with other diarizonae sequence types (ST) mainly-associated with infections in humans. The sseK2, sseK3 and slrP genes were either absent or showing frameshift mutations in a significant proportion of genomes from environmental diarizonae isolates. CONCLUSIONS: The comparative genomic study of two diarizonae strains isolated in Bolivia from human patients uncovered the presence of many genes putatively related to virulence. The statistically-significant acquisition of a unique combination of these functions by diarizonae strains isolated from humans may have impacted the ability of these isolates to successfully infect the human host.

Identification of the direct regulon of NtcA during early acclimation to nitrogen starvation in the cyanobacterium Synechocystis sp. PCC 6803.

In cyanobacteria, nitrogen homeostasis is maintained by an intricate regulatory network around transcription factor NtcA. Although mechanisms controlling NtcA activity appear to be well understood, its regulon remains poorly defined. To determine the NtcA regulon during the early stages of nitrogen starvation for the model cyanobacterium Synechocystis sp. PCC 6803, we performed chromatin immunoprecipitation, followed by sequencing (ChIP-seq), in parallel with transcriptome analysis (RNA-seq). Through combining these methods, we determined 51 genes activated and 28 repressed directly by NtcA. In addition to genes associated with nitrogen and carbon metabolism, a considerable number of genes without current functional annotation were among direct targets providing a rich reservoir for further studies. The NtcA regulon also included eight non-coding RNAs, of which Ncr1071, Syr6 and NsiR7 were experimentally validated, and their putative targets were computationally predicted. Surprisingly, we found substantial NtcA binding associated with delayed expression changes indicating that NtcA can reside in a poised state controlled by other factors. Indeed, a role of PipX as modulating factor in nitrogen regulation was confirmed for selected NtcA-targets. We suggest that the indicated poised state of NtcA enables a more differentiated response to nitrogen limitation and can be advantageous in native habitats of Synechocystis.

The CopRS two-component system is responsible for resistance to copper in the cyanobacterium Synechocystis sp. PCC 6803.

Photosynthetic organisms need copper for cytochrome oxidase and for plastocyanin in the fundamental processes of respiration and photosynthesis. However, excess of free copper is detrimental inside the cells and therefore organisms have developed homeostatic mechanisms to tightly regulate its acquisition, sequestration, and efflux. Herein we show that the CopRS two-component system (also known as Hik31-Rre34) is essential for copper resistance in Synechocystis sp. PCC 6803. It regulates expression of a putative heavy-metal efflux-resistance nodulation and division type copper efflux system (encoded by copBAC) as well as its own expression (in the copMRS operon) in response to the presence of copper in the media. Mutants in this two-component system or the efflux system render cells more sensitive to the presence of copper in the media and accumulate more intracellular copper than the wild type. Furthermore, CopS periplasmic domain is able to bind copper, suggesting that CopS could be able to detect copper directly. Both operons (copMRS and copBAC) are also induced by the photosynthetic inhibitor 2,5-dibromo-3-methyl-6-isopropyl-p-benzoquinone but this induction requires the presence of copper in the media. The reduced response of two mutant strains to copper, one lacking plastocyanin and a second one impaired in copper transport to the thylakoid, due to the absence of the P(I)-type ATPases PacS and CtaA, suggests that CopS can detect intracellular copper. In addition, a tagged version of CopS with a triple HA epitope localizes to both the plasma and the thylakoid membranes, suggesting that CopS could be involved in copper detection in both the periplasm and the thylakoid lumen.

Extracellular vesicles as an alternative copper-secretion mechanism in bacteria.

Metal homeostasis is fundamental for optimal performance of cell metabolic pathways. Over the course of evolution, several systems emerged to warrant an intracellular metal equilibrium. When exposed to growth-challenging copper concentrations, Gram-negative bacteria quickly activate copper-detoxification mechanisms, dependent on transmembrane-protein complexes and metallochaperones that mediate metal efflux. Here, we show that vesiculation is also a common bacterial response mechanism to high copper concentrations, and that extracellular vesicles (EVs) play a role in transporting copper. We present evidence that bacteria from different ecological niches release copious amounts of EVs when exposed to copper. Along with the activation of the classical detoxification systems, we demonstrate that copper-stressed cells of the cyanobacterium Synechocystis sp. PCC6803 release EVs loaded with the copper-binding metallochaperone CopM. Under standard growth conditions, CopM-loaded EVs could also be isolated from a Synechocystis strain lacking a functional TolC-protein, which we characterize here as exhibiting a copper-sensitive phenotype. Analyses of Synechocystis tolC-mutant's EVs isolated from cells cultivated under standard conditions indicated the presence of copper therein, in significantly higher levels as compared to those from the wild-type. Altogether, these results suggest that release of EVs in bacteria represent a novel copper-secretion mechanism, shedding light into alternative mechanisms of bacterial metal resistance.

ChIP-seq Experiment and Data Analysis in the Cyanobacterium Synechocystis sp. PCC 6803.

Nitrogen is an essential nutrient for all living organisms. In cyanobacteria, a group of oxygenic photosynthetic bacteria, nitrogen homeostasis is maintained by an intricate regulatory network around the transcription factor NtcA. Although mechanisms controlling NtcA activity appear to be well understood, the sets of genes under its control (i.e., its regulon) remain poorly defined. In this protocol, we describe the procedure for chromatin immunoprecipitation using NtcA antibodies, followed by DNA sequencing analysis (ChIP-seq) during early acclimation to nitrogen starvation in the cyanobacterium Synechocystis sp. PCC 6803 (hereafter Synechocystis). This protocol can be extended to analyze any DNA-binding protein in cyanobacteria for which suitable antibodies exist.

Draft Genome Sequence of Salmonella enterica subsp. enterica Serovar Infantis Strain SPE101, Isolated from a Chronic Human Infection.

We report a 4.99-Mb draft genome sequence of Salmonella enterica subsp. enterica serovar Infantis strain SPE101, isolated from feces of a 5-month-old breast-fed female showing diarrhea associated with severe dehydration and malnutrition. The infection prolonged for 6 months despite antibiotic treatment.

Ni interferes in the Cu-regulated transcriptional switch petJ/petE in Synechocystis sp. PCC 6803.

Plastocyanin (petE) plays an essential role in photosynthesis as an electron carrier between cytochrome b6 f and photosystem I, and in some cyanobacteria it can be replaced by the haem-containing protein, cytochrome c6 (petJ). In Synechocystis sp. PCC 6803, transcription of petE and petJ is activated and repressed, respectively, by Cu. Here, we show that Ni can act similarly to Cu in inducing petE and repressing petJ, thus leading to a partial switch between cytochrome c6 and plastocyanin. Transcription of these genes is only altered by Ni in Cu-depleted medium, and none of the Ni-dependent transcription factors described in Synechocystis, NrsR and InrS seem to be involved in this regulation. Finally, we show that plastocyanin is essential for growth under conditions of excess Ni.

Extracellular Proteins: Novel Key Components of Metal Resistance in Cyanobacteria?

Metals are essential for all living organisms and required for fundamental biochemical processes. However, when in excess, metals can turn into highly-toxic agents able to disrupt cell membranes, alter enzymatic activities, and damage DNA. Metal concentrations are therefore tightly controlled inside cells, particularly in cyanobacteria. Cyanobacteria are ecologically relevant prokaryotes that perform oxygenic photosynthesis and can be found in many different marine and freshwater ecosystems, including environments contaminated with heavy metals. As their photosynthetic machinery imposes high demands for metals, homeostasis of these micronutrients has been widely studied in cyanobacteria. So far, most studies have focused on how cells are capable of controlling their internal metal pools, with a strong bias toward the analysis of intracellular processes. Ultrastructure, modulation of physiology, dynamic changes in transcription and protein levels have been studied, but what takes place in the extracellular environment when cells are exposed to an unbalanced metal availability remains largely unknown. The interest in studying the subset of proteins present in the extracellular space has only recently begun and the identification and functional analysis of the cyanobacterial exoproteomes are just emerging. Remarkably, metal-related proteins such as the copper-chaperone CopM or the iron-binding protein FutA2 have already been identified outside the cell. With this perspective, we aim to raise the awareness that metal-resistance mechanisms are not yet fully known and hope to motivate future studies assessing the role of extracellular proteins on bacterial metal homeostasis, with a special focus on cyanobacteria.

The Transcriptional Landscape of the Photosynthetic Model Cyanobacterium Synechocystis sp. PCC6803.

Cyanobacteria exhibit a great capacity to adapt to different environmental conditions through changes in gene expression. Although this plasticity has been extensively studied in the model cyanobacterium Synechocystis sp. PCC 6803, a detailed analysis of the coordinated transcriptional adaption across varying conditions is lacking. Here, we report a meta-analysis of 756 individual microarray measurements conducted in 37 independent studies-the most comprehensive study of the Synechocystis transcriptome to date. Using stringent statistical evaluation, we characterized the coordinated adaptation of Synechocystis' gene expression on systems level. Evaluation of the data revealed that the photosynthetic apparatus is subjected to greater changes in expression than other cellular components. Nevertheless, network analyses indicated a significant degree of transcriptional coordination of photosynthesis and various metabolic processes, and revealed the tight co-regulation of components of photosystems I, II and phycobilisomes. Detailed inspection of the integrated data led to the discovery a variety of regulatory patterns and novel putative photosynthetic genes. Intriguingly, global clustering analyses suggested contrasting transcriptional response of metabolic and regulatory genes stress to conditions. The integrated Synechocystis transcriptome can be accessed and interactively analyzed via the CyanoEXpress website (http://cyanoexpress.sysbiolab.eu).

CopM is a novel copper-binding protein involved in copper resistance in Synechocystis sp. PCC 6803.

Copper resistance system in the cyanobacterium Synechocystis sp. PCC 6803 comprises two operons, copMRS and copBAC, which are expressed in response to copper in the media. copBAC codes for a heavy-metal efflux-resistance nodulation and division (HME-RND) system, while copMRS codes for a protein of unknown function, CopM, and a two-component system CopRS, which controls the expression of these two operons. Here, we report that CopM is a periplasmic protein able to bind Cu(I) with high affinity (KD ~3 × 10(-16) ). Mutants lacking copM showed a sensitive copper phenotype similar to mutants affected in copB, but lower than mutants of the two-component system CopRS, suggesting that CopBAC and CopM constitute two independent resistance mechanisms. Moreover, constitutive expression of copM is able to partially suppress the copper sensitivity of the copR mutant strain, pointing out that CopM per se is able to confer copper resistance. Furthermore, constitutive expression of copM was able to reduce total cellular copper content of the copR mutant to the levels determined in the wild-type (WT) strain. Finally, CopM was localized not only in the periplasm but also in the extracellular space, suggesting that CopM can also prevent copper accumulation probably by direct copper binding outside the cell.

The unfolded protein response and its potential role in Huntington's disease elucidated by a systems biology approach.

Huntington ´s disease (HD) is a progressive, neurodegenerative disease with a fatal outcome. Although the disease-causing gene (huntingtin) has been known for over 20 years, the exact mechanisms leading to neuronal cell death are still controversial. One potential mechanism contributing to the massive loss of neurons observed in the brain of HD patients could be the unfolded protein response (UPR) activated by accumulation of misfolded proteins in the endoplasmic reticulum (ER). As an adaptive response to counter-balance accumulation of un- or misfolded proteins, the UPR upregulates transcription of chaperones, temporarily attenuates new translation, and activates protein degradation via the proteasome. However, persistent ER stress and an activated UPR can also cause apoptotic cell death. Although different studies have indicated a role for the UPR in HD, the evidence remains inconclusive. Here, we present extensive bioinformatic analyses that revealed UPR activation in different experimental HD models based on transcriptomic data. Accordingly, we have identified 53 genes, including RAB5A, HMGB1, CTNNB1, DNM1, TUBB, TSG101, EEF2, DYNC1H1, SLC12A5, ATG5, AKT1, CASP7 and SYVN1 that provide a potential link between UPR and HD. To further elucidate the potential role of UPR as a disease-relevant process, we examined its connection to apoptosis based on molecular interaction data, and identified a set of 40 genes including ADD1, HSP90B1, IKBKB, IKBKG, RPS3A and LMNB1, which seem to be at the crossroads between these two important cellular processes. Remarkably, we also found strong correlation of UPR gene expression with the length of the polyglutamine tract of Huntingtin, which is a critical determinant of age of disease onset in human HD patients pointing to the UPR as a promising target for therapeutic intervention. The study is complemented by a newly developed web-portal called UPR-HD (http://uprhd.sysbiolab.eu) that enables visualization and interactive analysis of UPR-associated gene expression across various HD models.

Global transcriptional profiles of the copper responses in the cyanobacterium Synechocystis sp. PCC 6803.

Copper is an essential element involved in fundamental processes like respiration and photosynthesis. However, it becomes toxic at high concentration, which has forced organisms to control its cellular concentration. We have recently described a copper resistance system in the cyanobacterium Synechocystis sp. PCC 6803, which is mediated by the two-component system, CopRS, a RND metal transport system, CopBAC and a protein of unknown function, CopM. Here, we report the transcriptional responses to copper additions at non-toxic (0.3 µM) and toxic concentrations (3 µM) in the wild type and in the copper sensitive copR mutant strain. While 0.3 µM copper slightly stimulated metabolism and promoted the exchange between cytochrome c6 and plastocyanin as soluble electron carriers, the addition of 3 µM copper catalyzed the formation of ROS, led to a general stress response and induced expression of Fe-S cluster biogenesis genes. According to this, a double mutant strain copRsufR, which expresses constitutively the sufBCDS operon, tolerated higher copper concentration than the copR mutant strain, suggesting that Fe-S clusters are direct targets of copper toxicity in Synechocystis. In addition we have also demonstrated that InrS, a nickel binding transcriptional repressor that belong to the CsoR family of transcriptional factor, was involved in heavy metal homeostasis, including copper, in Synechocystis. Finally, global gene expression analysis of the copR mutant strain suggested that CopRS only controls the expression of copMRS and copBAC operons in response to copper.

Metals in cyanobacteria: analysis of the copper, nickel, cobalt and arsenic homeostasis mechanisms.

Traces of metal are required for fundamental biochemical processes, such as photosynthesis and respiration. Cyanobacteria metal homeostasis acquires an important role because the photosynthetic machinery imposes a high demand for metals, making them a limiting factor for cyanobacteria, especially in the open oceans. On the other hand, in the last two centuries, the metal concentrations in marine environments and lake sediments have increased as a result of several industrial activities. In all cases, cells have to tightly regulate uptake to maintain their intracellular concentrations below toxic levels. Mechanisms to obtain metal under limiting conditions and to protect cells from an excess of metals are present in cyanobacteria. Understanding metal homeostasis in cyanobacteria and the proteins involved will help to evaluate the use of these microorganisms in metal bioremediation. Furthermore, it will also help to understand how metal availability impacts primary production in the oceans. In this review, we will focus on copper, nickel, cobalt and arsenic (a toxic metalloid) metabolism, which has been mainly analyzed in model cyanobacterium Synechocystis sp. PCC 6803.