BAGET 2.0: an updated web tool for the effortless retrieval of prokaryotic gene context and sequence

Motivation

The retrieval of a single gene sequence and context from completely sequenced bacterial and archaeal genomes constitutes an intimidating task for the wet bench biologist. Existing web-based genome browsers are either too complex for routine use or only provide a subset of the available prokaryotic genomes.

Results

We have developed BAGET 2.0 (Bacterial and Archaeal Gene Exploration Tool), an updated web service granting access in just three mouse clicks to the sequence and synteny of any gene from completely sequenced bacteria and archaea. User-provided annotated genomes can be processed as well. BAGET 2.0 relies on a local database updated on a daily basis.

Availability and implementation

BAGET 2.0 befits all current browsers such as Chrome, Firefox, Edge, Opera and Safari. Internet Explorer 11 is supported. BAGET 2.0 is freely accessible at https://archaea.i2bc.paris-saclay.fr/baget/

Comparative genomics reveals conserved positioning of essential genomic clusters in highly rearranged Thermococcales chromosomes

The genomes of the 21 completely sequenced Thermococcales display a characteristic high level of rearrangements. As a result, the prediction of their origin and termination of replication on the sole basis of chromosomal DNA composition or skew is inoperative. Using a different approach based on biologically relevant sequences, we were able to determine oriC position in all 21 genomes. The position of dif, the site where chromosome dimers are resolved before DNA segregation could be predicted in 19 genomes. Computation of the core genome uncovered a number of essential gene clusters with a remarkably stable chromosomal position across species, in sharp contrast with the scrambled nature of their genomes. The active chromosomal reorganization of numerous genes acquired by horizontal transfer, mainly from mobile elements, could explain this phenomenon.

•

Thermococcales chromosomal landmarks were uncovered using biologically relevant sequences.

•

Core genomes procedures predict integration of mobile elements on Thermococcales chromosomes.

•

Thermococcales genomes are highly rearranged but core clusters positions remain invariable.

•

Thermococcales core genes are more expressed and predominantly encoded on the leading strand.

Operator recognition by the ROK transcription factor family members, NagC and Mlc

NagC and Mlc, paralogous members of the ROK family of proteins with almost identical helix-turn-helix DNA binding motifs, specifically regulate genes for transport and utilization of N-acetylglucosamine and glucose. We previously showed that two amino acids in a linker region outside the canonical helix-turn-helix motif are responsible for Mlc site specificity. In this work we identify four amino acids in the linker, which are required for recognition of NagC targets. These amino acids allow Mlc and NagC to distinguish between a C/G and an A/T bp at positions ±11 of the operators. One linker position, glycine in NagC and arginine in Mlc, corresponds to the major specificity determinant for the two proteins. In certain contexts it is possible to switch repression from Mlc-style to NagC-style, by interchanging this glycine and arginine. Secondary determinants are supplied by other linker positions or the helix-turn-helix motif. A wide genomic survey of unique ROK proteins shows that glycine- and arginine-rich sequences are present in the linkers of nearly all ROK family repressors. Conserved short sequence motifs, within the branches of the ROK evolutionary tree, suggest that these sequences could also be involved in operator recognition in other ROK family members.

Thermococcus nautili sp. nov., a hyperthermophilic archaeon isolated from a hydrothermal deep-sea vent

Thermococcus nautili, strain 30-1T (formerly reported as Thermococcus nautilus), was isolated from a hydrothermal chimney sample collected from the East Pacific Rise at a depth of 2633 m on the ‘La chainette PP57’ area. Cells were motile, irregular cocci with a polar tuft of flagella (0.8–1.5 µm) and divided by constriction. The micro-organism grew optimally at 87.5 °C (range 55–95 °C), at pH 7 (range pH 4–9) and with 2 % NaCl (range 1–4 %). Doubling time was 64 min in Zillig’s broth medium under optimal conditions. Growth was strictly anaerobic. It grew preferentially in the presence of elemental sulfur or cystine, which are reduced to H2S, on complex organic substrates such as yeast extract, tryptone, peptone, Casamino acids and casein. Slow growth was observed on starch and pyruvate. Strain 30-1T was resistant to chloramphenicol and tetracyclin (at 100 µg ml−1) but sensitive to kanamycin and rifampicin. The G+C content of the genomic DNA was 54 mol%. Strain 30-1T harboured three plasmids named pTN1, pTN2 and pTN3 and produced membrane vesicles that incorporate pTN1 and pTN3. As determined by 16S rRNA gene sequence analysis, strain 30-1T is related most closely to Thermococcus sp. AM4 (99.3 % similarity) and Thermococcus gammatolerans DSM 15229T (99.2 %). DNA–DNA hybridization values (in silico) with these two closest relatives were below the threshold value of 70 % (33 % with Thermococcus sp. AM4 and 32 % with T. gammatolerans DSM 15229T) and confirmed that strain 30-1 represents a novel species. On the basis of the data presented, strain 30-1T is considered to represent a novel species of the genus Thermococcus , for which the name Thermococcus nautili sp. nov. is proposed. The type strain is 30-1T ( = CNCM 4275 = JCM 19601).

MGcV: the microbial genomic context viewer for comparative genome analysis

Background

Conserved gene context is used in many types of comparative genome analyses. It is used to provide leads on gene function, to guide the discovery of regulatory sequences, but also to aid in the reconstruction of metabolic networks. We present the Microbial Genomic context Viewer (MGcV), an interactive, web-based application tailored to strengthen the practice of manual comparative genome context analysis for bacteria.

Results

MGcV is a versatile, easy-to-use tool that renders a visualization of the genomic context of any set of selected genes, genes within a phylogenetic tree, genomic segments, or regulatory elements. It is tailored to facilitate laborious tasks such as the interactive annotation of gene function, the discovery of regulatory elements, or the sequence-based reconstruction of gene regulatory networks. We illustrate that MGcV can be used in gene function annotation by visually integrating information on prokaryotic genes, like their annotation as available from NCBI with other annotation data such as Pfam domains, sub-cellular location predictions and gene-sequence characteristics such as GC content. We also illustrate the usefulness of the interactive features that allow the graphical selection of genes to facilitate data gathering (e.g. upstream regions, ID’s or annotation), in the analysis and reconstruction of transcription regulation. Moreover, putative regulatory elements and their corresponding scores or data from RNA-seq and microarray experiments can be uploaded, visualized and interpreted in (ranked-) comparative context maps. The ranked maps allow the interpretation of predicted regulatory elements and experimental data in light of each other.

Conclusion

MGcV advances the manual comparative analysis of genes and regulatory elements by providing fast and flexible integration of gene related data combined with straightforward data retrieval. MGcV is available at http://mgcv.cmbi.ru.nl.

SyntTax: a web server linking synteny to prokaryotic taxonomy

Background

The study of the conservation of gene order or synteny constitutes a powerful methodology to assess the orthology of genomic regions and to predict functional relationships between genes. The exponential growth of microbial genomic databases is expected to improve synteny predictions significantly. Paradoxically, this genomic data plethora, without information on organisms relatedness, could impair the performance of synteny analysis programs.

Results

In this work, I present SyntTax, a synteny web service designed to take full advantage of the large amount or archaeal and bacterial genomes by linking them through taxonomic relationships. SyntTax incorporates a full hierarchical taxonomic tree allowing intuitive access to all completely sequenced prokaryotes. Single or multiple organisms can be chosen on the basis of their lineage by selecting the corresponding rank nodes in the tree. The synteny methodology is built upon our previously described Absynte algorithm with several additional improvements.

Conclusions

SyntTax aims to produce robust syntenies by providing prompt access to the taxonomic relationships connecting all completely sequenced microbial genomes. The reduction in redundancy offered by lineage selection presents the benefit of increasing accuracy while reducing computation time. This web tool was used to resolve successfully several conserved complex gene clusters described in the literature. In addition, particular features of SyntTax permit the confirmation of the involvement of the four components constituting the E. coli YgjD multiprotein complex responsible for tRNA modification. By analyzing the clustering evolution of alternative gene fusions, new proteins potentially interacting with this complex could be proposed. The web service is available at http://archaea.u-psud.fr/SyntTax.

Junker: an intergenic explorer for bacterial genomes

In the past few decades, scientists from all over the world have taken a keen interest in novel functional units such as small regulatory RNAs, small open reading frames, pseudogenes, transposons, integrase binding attB/attP sites, repeat elements within the bacterial intergenic regions (IGRs) and in the analysis of those “junk” regions for genomic complexity. Here we have developed a web server, named Junker, to facilitate the in-depth analysis of IGRs for examining their length distribution, four-quadrant plots, GC percentage and repeat details. Upon selection of a particular bacterial genome, the physical genome map is displayed as a multiple loci with options to view any loci of interest in detail. In addition, an IGR statistics module has been created and implemented in the web server to analyze the length distribution of the IGRs and to understand the disordered grouping of IGRs across the genome by generating the four-quadrant plots. The proposed web server is freely available at the URL http://pranag.physics.iisc.ernet.in/junker/.

Absynte: a web tool to analyze the evolution of orthologous archaeal and bacterial gene clusters

SUMMARY

Absynte (Archaeal and Bacterial Synteny Explorer) is a web-based service designed to display local syntenies in completely sequenced prokaryotic chromosomes. The genomic contexts are determined with a multiple center star clustering topology on the basis of a user-provided protein sequence and all (or a set of) chromosomes from the publicly available archaeal and bacterial genomes. The results consist in a dynamic web page where a consistent color-coding permits a rapid visual evaluation of the relative positioning of genes with similar sequences within the synteny. Each gene composing the synteny can be further queried interactively using either local or remote databases. Absynte results can be exported in .CSV or high-resolution, .PDF formats for printing, archival, further editing or publication purposes. Performance, real-time computation, user-friendliness and daily database updates constitute the principal advantages of Absynte over similar web services.

AVAILABILITY

http://archaea.u-psud.fr/absynte

CONTACT

jacques.oberto@igmors.u-psud.fr.

FITBAR: a web tool for the robust prediction of prokaryotic regulons

Background

The binding of regulatory proteins to their specific DNA targets determines the accurate expression of the neighboring genes. The in silico prediction of new binding sites in completely sequenced genomes is a key aspect in the deeper understanding of gene regulatory networks. Several algorithms have been described to discriminate against false-positives in the prediction of new binding targets; however none of them has been implemented so far to assist the detection of binding sites at the genomic scale.

Results

FITBAR (Fast Investigation Tool for Bacterial and Archaeal Regulons) is a web service designed to identify new protein binding sites on fully sequenced prokaryotic genomes. This tool consists in a workbench where the significance of the predictions can be compared using different statistical methods, a feature not found in existing resources. The Local Markov Model and the Compound Importance Sampling algorithms have been implemented to compute the P-value of newly discovered binding sites. In addition, FITBAR provides two optimized genomic scanning algorithms using either log-odds or entropy-weighted position-specific scoring matrices. Other significant features include the production of a detailed genomic context map for each detected binding site and the export of the search results in spreadsheet and portable document formats. FITBAR discovery of a high affinity Escherichia coli NagC binding site was validated experimentally in vitro as well as in vivo and published.

Conclusions

FITBAR was developed in order to allow fast, accurate and statistically robust predictions of prokaryotic regulons. This feature constitutes the main advantage of this web tool over other matrix search programs and does not impair its performance. The web service is available at http://archaea.u-psud.fr/fitbar.

Qri7/OSGEPL, the mitochondrial version of the universal Kae1/YgjD protein, is essential for mitochondrial genome maintenance

Yeast Qri7 and human OSGEPL are members of the orthologous Kae1(OSGEP)/YgjD protein family, the last class of universally conserved proteins without assigned function. Phylogenetic analyses indicate that the eukaryotic Qri7(OSGEPL) proteins originated from bacterial YgjD proteins. We have recently shown that the archaeal Kae1 protein is a DNA-binding protein that exhibits apurinic endonuclease activity in vitro. We show here that the Qri7/OSGEPL proteins localize in mitochondria and are involved in mitochondrial genome maintenance in two model eukaryotic organisms, Saccharomyces cerevisiae and Caenorhabditis elegans. Furthermore, S. cerevisiae Qri7 complements the loss of the bacterial YgjD protein in Escherichia coli, suggesting that Qri7/OSGEPL and YgjD proteins have retained similar functions in modern organisms. We suggest to name members of the Kae1(OSGEP)/YgjD family UGMP, for Universal Genome Maintenance Proteins.

The HU regulon is composed of genes responding to anaerobiosis, acid stress, high osmolarity and SOS induction

BACKGROUND

The Escherichia coli heterodimeric HU protein is a small DNA-bending protein associated with the bacterial nucleoid. It can introduce negative supercoils into closed circular DNA in the presence of topoisomerase I. Cells lacking HU grow very poorly and display many phenotypes.

METHODOLOGY/PRINCIPAL FINDINGS

We analyzed the transcription profile of every Escherichia coli gene in the absence of one or both HU subunits. This genome-wide in silico transcriptomic approach, performed in parallel with in vivo genetic experimentation, defined the HU regulon. This large regulon, which comprises 8% of the genome, is composed of four biologically relevant gene classes whose regulation responds to anaerobiosis, acid stress, high osmolarity, and SOS induction.

CONCLUSIONS/SIGNIFICANCE

The regulation a large number of genes encoding enzymes involved in energy metabolism and catabolism pathways by HU explains the highly pleiotropic phenotype of HU-deficient cells. The uniform chromosomal distribution of the many operons regulated by HU strongly suggests that the transcriptional and nucleoid architectural functions of HU constitute two aspects of a unique protein-DNA interaction mechanism.