Ranko Gacesa

KEGG orthology-based annotation of the predicted proteome of Acropora digitifera: ZoophyteBase - an open access and searchable database of a coral genome

BackgroundContemporary coral reef research has firmly established that a genomic approach is urgently needed to better understand the effects of anthropogenic environmental stress and global climate change on coral holobiont interactions. Here we present KEGG orthology-based annotation of the complete genome sequence of the scleractinian coral Acropora digitifera and provide the first comprehensive view of the genome of a reef-building coral by applying advanced bioinformatics.DescriptionSequences from the KEGG database of protein function were used to construct hidden Markov models. These models were used to search the predicted proteome of A. digitifera to establish complete genomic annotation. The annotated dataset is published in ZoophyteBase, an open access format with different options for searching the data. A particularly useful feature is the ability to use a Google-like search engine that links query words to protein attributes. We present features of the annotation that underpin the molecular structure of key processes of coral physiology that include (1) regulatory proteins of symbiosis, (2) planula and early developmental proteins, (3) neural messengers, receptors and sensory proteins, (4) calcification and Ca2+-signalling proteins, (5) plant-derived proteins, (6) proteins of nitrogen metabolism, (7) DNA repair proteins, (8) stress response proteins, (9) antioxidant and redox-protective proteins, (10) proteins of cellular apoptosis, (11) microbial symbioses and pathogenicity proteins, (12) proteins of viral pathogenicity, (13) toxins and venom, (14) proteins of the chemical defensome and (15) coral epigenetics.ConclusionsWe advocate that providing annotation in an open-access searchable database available to the public domain will give an unprecedented foundation to interrogate the fundamental molecular structure and interactions of coral symbiosis and allow critical questions to be addressed at the genomic level based on combined aspects of evolutionary, developmental, metabolic, and environmental perspectives.

Draft Genome Sequence of Streptomyces rapamycinicus Strain NRRL 5491, the Producer of the Immunosuppressant Rapamycin

ABSTRACT Streptomyces rapamycinicus strain NRRL 5491 produces the important drug rapamycin. It has a large genome of 12.7 Mb, of which over 3 Mb consists of 48 secondary metabolite biosynthesis clusters.

Predicting substrate specificity of adenylation domains of nonribosomal peptide synthetases and other protein properties by latent semantic indexing

Successful genome mining is dependent on accurate prediction of protein function from sequence. This often involves dividing protein families into functional subtypes (e.g., with different substrates). In many cases, there are only a small number of known functional subtypes, but in the case of the adenylation domains of nonribosomal peptide synthetases (NRPS), there are >500 known substrates. Latent semantic indexing (LSI) was originally developed for text processing but has also been used to assign proteins to families. Proteins are treated as ‘‘documents’’ and it is necessary to encode properties of the amino acid sequence as ‘‘terms’’ in order to construct a term-document matrix, which counts the terms in each document. This matrix is then processed to produce a document-concept matrix, where each protein is represented as a row vector. A standard measure of the closeness of vectors to each other (cosines of the angle between them) provides a measure of protein similarity. Previous work encoded proteins as oligopeptide terms, i.e. counted oligopeptides, but used no information regarding location of oligopeptides in the proteins. A novel tokenization method was developed to analyze information from multiple alignments. LSI successfully distinguished between two functional subtypes in five well-characterized families. Visualization of different ‘‘concept’’ dimensions allows exploration of the structure of protein families. LSI was also used to predict the amino acid substrate of adenylation domains of NRPS. Better results were obtained when selected residues from multiple alignments were used rather than the total sequence of the adenylation domains. Using ten residues from the substrate binding pocket performed better than using 34 residues within 8xa0Å of the active site. Prediction efficiency was somewhat better than that of the best published method using a support vector machine.

Rising levels of atmospheric oxygen and evolution of Nrf2

In mammals, the master transcription regulator of antioxidant defences is provided by the Nrf2 protein. Phylogenetic analyses of Nrf2 sequences are used here to derive a molecular clock that manifests persuasive evidence that Nrf2 orthologues emerged, and then diverged, at two time points that correlate with well-established geochemical and palaeobiological chronologies during progression of the ‘Great Oxygenation Event’. We demonstrate that orthologues of Nrf2 first appeared in fungi around 1.5u2009Ga during the Paleoproterozoic when photosynthetic oxygen was being absorbed into the oceans. A subsequent significant divergence in Nrf2 is seen during the split between fungi and the Metazoa approximately 1.0–1.2u2009Ga, at a time when oceanic ventilation released free oxygen to the atmosphere, but with most being absorbed by methane oxidation and oxidative weathering of land surfaces until approximately 800u2009Ma. Atmospheric oxygen levels thereafter accumulated giving rise to metazoan success known as the Cambrian explosion commencing at ~541u2009Ma. Atmospheric O2 levels then rose in the mid Paleozoic (359–252u2009Ma), and Nrf2 diverged once again at the division between mammals and non-mammalian vertebrates during the Permian-Triassic boundary (~252u2009Ma). Understanding Nrf2 evolution as an effective antioxidant response may have repercussions for improved human health.

Bacterial diversity of polluted surface sediments in the northern Adriatic Sea

Samples were collected from sea sediments at seven sites in the northern Adriatic Sea that included six sites next to industrial complexes and one from a tourist site (recreational beach). The samples were assayed for alkanes and polycyclic aromatic hydrocarbons. The composition of the hydrocarbon samples suggested that industrial pollution was present in most cases. A sample from one site was also grown aerobically under crude oil enrichment in order to evaluate the response of indigenous bacterial populations to crude oil exposure. Analysis of 16S rRNA gene sequences showed varying microbial biodiversity depending on the level of pollution--ranging from low (200 detected genera) to high (1000+ genera) biodiversity, with lowest biodiversity observed in polluted samples. This indicated that there was considerable biodiversity in all sediment samples but it was severely restricted after exposure to crude oil selection pressure. Phylogenetic analysis of putative alkB genes showed high evolutionary diversity of the enzymes in the samples and suggested great potential for bioremediation and bioprospecting. The first systematic analysis of bacterial communities from sediments of the northern Adriatic Sea is presented, and it will provide a baseline assessment that may serve as a reference point for ecosystem changes and hydrocarbon degrading potential--a potential that could soon gain importance due to plans for oil exploitation in the area.

Gene duplications are extensive and contribute significantly to the toxic proteome of nematocysts isolated from Acropora digitifera (Cnidaria: Anthozoa: Scleractinia)

BackgroundGene duplication followed by adaptive selection is a well-accepted process leading to toxin diversification in venoms. However, emergent genomic, transcriptomic and proteomic evidence now challenges this role to be at best equivocal to other processess . Cnidaria are arguably the most ancient phylum of the extant metazoa that are venomous and such provide a definitive ancestral anchor to examine the evolution of this trait.MethodsHere we compare predicted toxins from the translated genome of the coral Acropora digitifera to putative toxins revealed by proteomic analysis of soluble proteins discharged from nematocysts, to determine the extent to which gene duplications contribute to venom innovation in this reef-building coral species. A new bioinformatics tool called HHCompare was developed to detect potential gene duplications in the genomic data, which is made freely available (https://github.com/rgacesa/HHCompare).ResultsA total of 55 potential toxin encoding genes could be predicted from the A. digitifera genome, of which 36 (65xa0%) had likely arisen by gene duplication as evinced using the HHCompare tool and verified using two standard phylogeny methods. Surprisingly, only 22xa0% (12/55) of the potential toxin repertoire could be detected following rigorous proteomic analysis, for which only half (6/12) of the toxin proteome could be accounted for as peptides encoded by the gene duplicates. Biological activities of these toxins are dominatedby putative phospholipases and toxic peptidases.ConclusionsGene expansions in A. digitifera venom are the most extensive yet described in any venomous animal, and gene duplication plays a significant role leading to toxin diversification in this coral species. Since such low numbers of toxins were detected in the proteome, it is unlikely that the venom is evolving rapidly by prey-driven positive natural selection. Rather we contend that the venom has a defensive role deterring predation or harm from interspecific competition and overgrowth by fouling organisms. Factors influencing translation of toxin encoding genes perhaps warrants more profound experimental consideration.

Evolutionary concepts in natural products discovery: what actinomycetes have taught us

Actinomycetes are a very important source of natural products for the pharmaceutical industry and other applications. Most of the strains belong to Streptomyces or related genera, partly because they are particularly amenable to growth in the laboratory and industrial fermenters. It is unlikely that chemical synthesis can fulfil the needs of the pharmaceutical industry for novel compounds so there is a continuing need to find novel natural products. An evolutionary perspective can help this process in several ways. Genome mining attempts to identify secondary metabolite biosynthetic clusters in DNA sequences, which are likely to produce interesting chemical entities. There are often technical problems in assembling the DNA sequences of large modular clusters in genome and metagenome projects, which can be overcome partially using information about the evolution of the domain sequences. Understanding the evolutionary mechanisms of modular clusters should allow simulation of evolutionary pathways in the laboratory to generate novel compounds.

Databases of the thiotemplate modular systems (CSDB) and their in silico recombinants (r-CSDB)

Modular biosynthetic clusters are responsible for the synthesis of many important pharmaceutical products. They include polyketide synthases (PKS clusters), non-ribosomal synthetases (NRPS clusters), and mixed clusters (containing both PKS and NRPS modules). The ClustScan database (CSDB) contains highly annotated descriptions of 170 clusters. The database has a hierarchical organization, which allows easy extraction of DNA and protein sequences of polypeptides, modules, and domains as well as an organization of the annotation so as to be able to predict the product chemistry to view it or export it in a standard SMILES format. The recombinant ClustScan database contains information about predicted recombinants between PKS clusters. The recombinants are generated by modeling homologous recombination and are associated with annotation and prediction of product chemistry automatically generated by the model. The database contains over 20,000 recombinants and is a resource for in silico approaches to detecting promising new compounds. Methods are available to construct the corresponding recombinants in the laboratory.

Bioinformatics analyses provide insight into distant homology of the Keap1-Nrf2 pathway.

An essential requirement for the evolution of early eukaryotic life was the development of effective means to protect against metabolic oxidative stress and exposure to environmental toxicants. In present-day mammals, the master transcription factor Nrf2 regulates basal level homeostasis and inducible expression of numerous detoxifying and antioxidant genes. To examine early evolution of the Keap1-Nrf2 pathway, we present bioinformatics analyses of distant homology of mammalian Keap1 and Nrf2 proteins across the Kingdoms of Life. Software written for this analysis is made freely available on-line. Furthermore, utilizing protein modeling and virtual screening methods, we demonstrate potential for Nrf2 activation by competitive inhibition of its binding to Keap1, specifically by UV-protective fungal mycosporines and marine mycosporine-like amino acids (MAAs). We contend that coevolution of Nrf2-activating secondary metabolites by fungi and other extant microbiota may provide prospective compound leads for the design of new therapeutics to target activation of the human Keap1-Nrf2 pathway for treating degenerative diseases of ageing.

O-Methyltransferase Is Shared between the Pentose Phosphate and Shikimate Pathways and Is Essential for Mycosporine-Like Amino Acid Biosynthesis in Anabaena variabilis ATCC 29413

The parent core structure of mycosporine‐like amino acids (MAAs) is 4‐deoxygadusol, which, in cyanobacteria, is derived from conversion of the pentose phosphate pathway intermediate sedoheptulose 7‐phosphate by the enzymes 2‐epi‐5‐epivaliolone synthase (EVS) and O‐methyltransferase (OMT). Yet, deletion of the EVS gene from Anabaena variabilis ATCC 29413 was shown to have little effect on MAA production, thus suggesting that its biosynthesis is not exclusive to the pentose phosphate pathway. Herein, we report how, using pathway‐specific inhibitors, we demonstrated unequivocally that MAA biosynthesis occurs also via the shikimate pathway. In addition, complete in‐frame gene deletion of the OMT gene from A. variabilis ATCC 29413 reveals that, although biochemically distinct, the pentose phosphate and shikimate pathways are inextricably linked to MAA biosynthesis in this cyanobacterium. Furthermore, proteomic data reveal that the shikimate pathway is the predominate route for UV‐induced MAA biosynthesis.