Next Article in Journal
Neuronal Nicotinic Acetylcholine Receptor Modulators from Cone Snails
Next Article in Special Issue
Cnidarian Interaction with Microbial Communities: From Aid to Animal’s Health to Rejection Responses
Previous Article in Journal
Smenamide A Analogues. Synthesis and Biological Activity on Multiple Myeloma Cells
Previous Article in Special Issue
Proteomic Analyses of the Unexplored Sea Anemone Bunodactis verrucosa
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Marine Drugs
Volume 16
Issue 6
10.3390/md16060207
marinedrugs-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Holo-Transcriptome of the Zoantharian Protopalythoa variabilis (Cnidaria: Anthozoa): A Plentiful Source of Enzymes for Potential Application in Green Chemistry, Industrial and Pharmaceutical Biotechnology

by
Jean-Étienne R. L. Morlighem
1,2,
Chen Huang
3,
Qiwen Liao
3,
Paula Braga Gomes
4,
Carlos Daniel Pérez
5,*,
Álvaro Rossan De Brandão Prieto-da-Silva
6,
Simon Ming-Yuen Lee
3 and
Gandhi Rádis-Baptista
1,2,*
1
Northeast Biotechnology Network (RENORBIO), Post-Graduation Program in Biotechnology, Federal University of Ceará, Fortaleza 60440-900, Brazil
2
Laboratory of Biochemistry and Biotechnology, Institute for Marine Sciences, Federal University of Ceará, Fortaleza 60165-081, Brazil
3
State Key Laboratory of Quality Research in Chinese Medicine and Institute of Chinese Medical Sciences, University of Macau, Macau 519020, China
4
Department of Biology, Federal Rural University of Pernambuco, Recife 52171-900, Brazil
5
Academic Center in Vitória, Federal University of Pernambuco, Vitória de Santo Antão 50670-901, Pernambuco, Brazil
6
Laboratory of Genetics, Butantan Institute, São Paulo 05503-900, Brazil
*
Authors to whom correspondence should be addressed.
Mar. Drugs 2018, 16(6), 207; https://doi.org/10.3390/md16060207
Submission received: 1 May 2018 / Revised: 5 June 2018 / Accepted: 8 June 2018 / Published: 13 June 2018
(This article belongs to the Special Issue Jellyfish and Polyps: Cnidarians as Sustainable Resources for Biotechnological Applications and Bioprospecting)
Browse Figures
Versions Notes

Abstract

:
Marine invertebrates, such as sponges, tunicates and cnidarians (zoantharians and scleractinian corals), form functional assemblages, known as holobionts, with numerous microbes. This type of species-specific symbiotic association can be a repository of myriad valuable low molecular weight organic compounds, bioactive peptides and enzymes. The zoantharian Protopalythoa variabilis (Cnidaria: Anthozoa) is one such example of a marine holobiont that inhabits the coastal reefs of the tropical Atlantic coast and is an interesting source of secondary metabolites and biologically active polypeptides. In the present study, we analyzed the entire holo-transcriptome of P. variabilis, looking for enzyme precursors expressed in the zoantharian-microbiota assemblage that are potentially useful as industrial biocatalysts and biopharmaceuticals. In addition to hundreds of predicted enzymes that fit into the classes of hydrolases, oxidoreductases and transferases that were found, novel enzyme precursors with multiple activities in single structures and enzymes with incomplete Enzyme Commission numbers were revealed. Our results indicated the predictive expression of thirteen multifunctional enzymes and 694 enzyme sequences with partially characterized activities, distributed in 23 sub-subclasses. These predicted enzyme structures and activities can prospectively be harnessed for applications in diverse areas of industrial and pharmaceutical biotechnology.

1. Introduction

The worldwide demand for useful enzymes is continuously increasing in various industries, including the pharmaceutical sector, green chemistry, fine chemicals and basic and applied biomedical research. The market for enzymes that are more catalytically efficient than currently used enzymes, more environmentally friendly, and have potential use as drugs in the pharmaceutical industry, is foreseen to reach over USD 17 billion, by 2024 [1]. Additionally, the global market for enzymes that have applications in molecular biology and analytical kits is estimated to reach approximately USD 12.5 billion by 2021 [2]. In green chemistry, biocatalysts are used in environmentally benign chemical synthesis, in particular, halogenation and transamination, and the future of enzyme-catalyzed reactions is foreseen to rely on the enantio (selective) synthesis and kinetic resolution of chemicals, preferentially from renewable sources [3]. The development of eco-friendly and greener biocatalyzed processes aims to achieve better environmental factor values, i.e., mass of waste per unit of product [4]. Therefore, enzymes with improved properties, including physical and chemical stability, high conversion rates, change of substrate specificity, stereoselectivity, independence of cofactors for activity, ability to produce new chemicals, possibility to assemble multi-enzyme complex, among others, are desirable. Enzymes with these attributes and properties can be prospected by screening microbiomes [5,6] and/or by in vitro mutagenesis followed by direct evolution [7,8].
With the advent of “omics” sciences, the search for novel enzymes has progressed primarily by means of high throughput metagenomics assays of uncultured environmental samples [9]. In addition, the untapped biodiversity of distinct marine environments has become a hot target source for a wide variety of bioactive compounds, particularly, enzymes, that are halotolerant, thermostable and adapted to a range of pressures and substrate selectivity [10]. Some recent examples of marine-derived enzymes include the new flavin-dependent halogenase from a marine sponge metagenome [11] and several new α-amylases isolated from a sea anemone microbial community [12]. In two recent articles, the potential of marine biomes, particularly microbiomes, was emphasized, especially with regard to marine enzymes [13,14]. The successful bioprospection of natural compounds from marine species has been reported primarily from invertebrates of the phyla Porifera and Cnidaria [15]. Most of these compounds, which were secondary metabolites that were made into chemotherapeutics or drug-leads, were originally isolated from symbiotic microorganisms, rather than their host, in the holobiont assemblage [16]. Sponges, cnidarians, tunicates and other marine invertebrates can harbor a great diversity of microbial symbionts [17,18,19]. Despite the microbial origin of many compounds derived from the holobiont assemblage, the own coral tissues comprise unique resources of diverse chemicals with distinct pharmacological activities, such as anti-inflammatory and anti-proliferative activities [20,21,22]. Therefore, prospecting biopharmaceuticals from unusual marine species, as sources of unique enzymes, focusing particularly on improved and novel biocatalysts, is also warranted.
Our local, potential of discovery, biodiverse sites, i.e., the Brazilian large marine ecosystems—namely, the Brazilian Shelves, which are rich in coral reefs and are marine biodiversity hotspots, have been largely under-studied, particularly concerning to the discovery of novel enzymes and even natural products [15,23]. Recently, we conducted whole RNA sequencing of the anthozoan Protopalythoa variabilis (Cnidaria: Anthozoa) and described a repertoire of bioactive peptides with cnidarian toxin features in its transcriptome [24]. P. variabilis is a zoantharian species that preferentially grows in the shallow, warm water of the Atlantic coast and is present in abundance along the coastal reefs of the Brazilian Northeast [25]. Colonies of P. variabilis, like other anthozoans, harbor symbiotic zooxanthellae and a consortium of other microorganisms [26,27,28]. In the present work, we investigated the expressed enzymatic content of the P. variabilis holo-transcriptome. This species of zoantharian belongs to the phylum Cnidaria, one of the first groups that diverged from the Bilateria and is positioned at the base of Metazoa. The holo-transcriptomic analysis of P. variabilis extends the prospection of marine organisms for biotechnological studies and biopharmaceutical applications, which are less numerous than those of terrestrial and microbial origin [29]. In contrast to other analytical strategies reported in most recent articles, we focused our investigation on the entire transcriptome of P. variabilis, from a holobiont perspective, searching for enzyme precursors expressed in the zoantharian-microbiota assemblage that resulted in the prediction of numerous enzyme sequences relevant to biotechnology and green chemistry. These putative enzymes include oxidoreductases, transferases, hydrolases, lyases, isomerases and ligases, which have potential applications in several industrial fields, such as the production of pharmaceutics and fine chemicals, bioconversion and biopolymers, and green chemistry, to mention a few. Moreover, as exemplified in the present work, a single zoantharian species can be viewed as a species-specific repository of a unique collection of marine enzymes.

2. Results and Discussion

2.1. Biodiversity in the P. variabilis Holobiont

Coral reefs are niches for different life forms, ranging from small fishes and crustaceans to associated and endosymbiotic microbial communities. Unraveling the biodiversity of a zoantharian holobiont is interesting not only from the ecological point of view, but also essential for the comprehension of the interconnected metabolic pathways, which ultimately depend on the symbiotic interactions and their enzymatic activities. To verify the overall diversity of life forms in the P. variabilis assemblage, we used three gene sequences commonly accepted for barcoding in species identification: the mitochondrial 16S rRNA gene for the taxonomic identification of bacteria and archaea, the mitochondrial cytochrome C oxidase subunit I (COI) gene for animals, algae, and dinoflagellates, and RuBisCo (rbcL) gene for plant and microalgae. The results of the species identification in the holobiont are shown in Figure 1A and the species are listed in Supplementary Table S1. In addition to the presence of the most common dinoflagellate algae of the genus Symbiodinium, involved in mutualistic symbiosis with cnidarians, the largest majority of the holobiont community identified is composed of uncultured species of cyanobateria and proteobacteria. Interestingly, four COI sequences found in the P. variabilis holobiont transcriptome had their best hits against terrestrial flying insects (Endopterygota), seemingly suggesting the interaction of P. variabilis with a group of organisms lacking known genetic information that belong to the phylum Arthropoda.

2.2. Annotation and Functional Prediction

Approximately 40% of the unigenes (51,792 sequences) identified in this study matched to entries in at least one of the four public protein databases, corresponding to “predicted or “annotated” hits (89% of all contigs) and “hypothetical” or “putative hits” (11%).
Supplementary Figures S1 and S2 present the complete statistical analysis of the sequence annotations. Compared to other de novo transcriptome assembly studies of cnidarian species, the initial level of BLAST annotations observed here was in the expected range [30,31,32,33]. Based on the selected BLASTx matches for the annotations, the taxonomic distribution of species from which the predicted protein sequences originated is depicted in Figure 1B. As can be observed, the three most representative species belong to the phylum Cnidaria, class Anthozoa, subclass Hexacorallia (namely, Exaiptasia pallida, Nematostella vectensis, and Acropora digitifera), which together accounted for 21,852 matched contigs in the annotation process for the P. variabilis holobiont transcripts (~42% of all annotated transcripts).

2.2.1. Functional Classification of P. variabilis Predicted Enzymes

After a GO-slim reduction, a total of 29,866 unigenes (~58% of all annotated transcripts) were classified within 35 different ontological categories (Supplementary Figure S3A). The functional annotation of the P. variabilis holo-transcriptome, which returned GO terms for almost 30,000 unigenes, was substantially higher than that observed in the majority of previous studies but collectively consistent regarding the gene distribution within the categories of cellular or metabolic processes and functions.

2.2.2. Assignment of Metabolic Pathways of the P. variabilis Predicted Enzymes

The KEGG pathway mapping analysis of the annotated unigenes retrieved 135 pathway maps (Supplementary Figure S3B), which were similar to the GO annotations in the category of “Biological Process”. However, the coverage differed greatly (Supplementary Table S2). These maps are available in Supplementary File S1.
From these analyses, a relatively high number of contigs were determined to encode enzymes involved in the metabolism of terpenoids and polyketides (Supplementary Table S2). This finding is interesting from the viewpoint of drug discovery, since numerous important bioactive secondary metabolites used in therapy, such as taxol, steroids and macrolide antibiotics, prostaglandins and tetracyclines, are synthesized in these biosynthetic routes. Indeed, the C-15 epimer of prostaglandin A2 (PGA2) and related structures have been identified and isolated from the Caribbean coral Plexaura homomalla, which produces and accumulates such compound [34].

2.3. Biotechnologically Relevant Enzymes from the P. variabilis Holo-Transcriptome

The identification of industrially relevant enzyme biocatalysts from environmental samples has been concretized by means of metagenomic analysis [9,35,36], a combination of metagenomics and metaproteomics technologies [37], and the combined application of transcriptomics and proteomics [38]. We based our search on analyzing the P. variabilis holo-transcriptome, providing an insight into the unexplored enzymes and related biosynthesis pathways expressed in this zoantharian-microbiota assemblage. Among all predicted proteins, 771 predicted enzymes classes with complete Enzyme Commission (EC) numbers were identified from the P. variabilis holo-transcriptome, representing a cumulative number of 6571 unigenes, with 606 of the predicted enzymes classes (5385 unigenes) mapped to KEGG pathways (Supplementary Table S3). Breaking down the enzymes by classes, 22% of them were oxidoreductases—EC 1 (170 enzymes for 806 unigenes), 35% were transferases—EC 2 (266 enzymes for 1208 unigenes), 21% were hydrolases—EC 3 (162 enzymes for 3792 unigenes), 10% were lyases—EC 4 (78 enzymes for 280 unigenes), 5% were isomerases—EC 5 (38 enzymes for 247 unigenes)), and 7% were ligases—EC 6 (57 enzymes for 238 unigenes).
Based on these results, we focused our study on enzymes with recognized relevance in chemical and pharmaceutical industries, as well as for green chemistry. For instance, transaminases (or aminotransferases) are stereoselective, cofactor-recycling enzymes that catalyze the transfer or exchange of an amino group from an amino-containing substrate to an acceptor molecule, resulting in the synthesis of chiral amino acids and amines. According to Supplementary Table S3, seven contigs encoding transaminases retrieved from the P. variabilis holo-transcriptome mapped mainly to antibiotic and amino acid biosynthesis pathways. The search for and development of optimized and high-performance transaminases that exhibit catalytic polyvalence (activity with a wide range of substrates), reactional stability and the possibility to form multi-enzyme complexes, are in demand for applications in green organic chemistry and the production of fine chemicals, food additives and pharmaceuticals [39]. In the following sections, additional selected examples of predicted enzymatic activities from the P. variabilis holobiont are presented, classified according to the major commercial and industrial application and discussed.

2.3.1. Relevance in the Treatment of Rare Diseases and Other Biomedical and Clinical Applications

Several classes of hydrolases, such as galactosidases and proteases, are used in therapy and the preparation of biopharmaceuticals. For instance, fibrinolytic enzymes are used clinically as thrombolytic agents to treat myocardial infarction, asparaginase and arginine deaminase are used to treat leukemia and solid tumors, and α- and β-galactosidases are used to treat lysosomal storage disease [40]. Moreover, naïve or structure-guided directly evolved glycohydrolases have been tested in trials to remove sugar residues from the surface of erythrocytes to obtain universal blood [41]. Supplementary Table S4 presents, a list of predicted P. variabilis enzymes for which counterparts are used in the treatment of rare metabolic diseases, cancer, and for experimental blood production and organ transplantation. Several putative enzymes in the P. variabilis holo-transcriptome are homologous to therapeutic enzymes currently used to treat lysosomal storage disorders (LSDs). LSDs encompass a large number of metabolic diseases, primarily characterized by a lack of hydrolases and defects in the degradation of glycoproteins, glycolipid, glycosaminoglycan and glycogen [42]. Importantly, replacement therapy with human recombinant enzymes has been proven to be effective treatment in clinical and pre-clinical settings [40,42,43].
In the fields of green chemistry, industrial and pharmaceutical biotechnology, glucosidases have been used to prepare glucopolymers of polyvinyl saccharides, such as poly(glucosyl)-acrylates, which function as drug carrier systems and non-ionic polymeric surfactants. For example, Kloosterman and colleagues (2014) [44] utilized β-glucosidase to synthesize the monomers 4-(β-glucosyloxy)-butyl acrylate, 2-(β-glucosyloxy)-ethyl acrylate and methacrylate, as a means to avoid multiple reaction steps, thereby precluding the formation of undesirable isomers. β-glucosidases were also predicted from the P. variabilis holo-transcriptome (Supplementary Table S4).
Two important therapeutic enzymes used in cancer therapy are l-asparaginase and arginine deaminase, for which putative homologues were also identified in this study (Supplementary Table S4). l-asparaginase catalyzes the conversion of l-asparagine into l-aspartate, preventing cancer cell survival in patients with lymphoblastic leukemia, while arginine deaminase catalytically removes, by deamination, excess adenosine molecules in the circulation of patients with cancer, thereby reducing the toxicity to the immune system caused by elevated adenosine levels [45]. In the preparation of universal O-type blood, the selective removal of α-GalNAc and α-Gal residues of the A and B oligosaccharide antigens on the surface of red blood cells can be achieved by α-galactosidases and α-N-acetylgalactosaminidases, respectively [46]. Eleven sequences from the P. variabilis holo-transcriptome mapped to a putative β-N-acetylhexosaminidase (EC: 3.2.1.52). One of the eleven sequences is closely related to the clade Cnidaria, whereas another belongs to the protist-algae group, and three others are in the archaea-bacteria clade. The last one, Unigene28224, fits in between the protist-algae and archaea-bacteria groups (Supplementary Figure S4A). Despite the divergence of the sequences, the proton donor glutamic acid residues in the catalytic site were found to be conserved across all P. variabilis sequences (Supplementary Figure S4B).
Many enzymes and therapeutic proteins used in clinical and experimental clinical trials exist in PEGylated forms, in which serum stability and half-life are usually increased, while the immunological response is decreased. The PEGylation of proteins can potentially be achieved by biocatalysis using transglutaminases, which carry out an acyl transfer reaction [47]. As noted on Supplementary Table S4, a transglutaminase (EC 2.3.2.13) was found that may be used to catalyze the covalent transfer of the PEG moiety to therapeutic enzymes and proteins. Transglutaminases are also useful in other biotech industries, such as food processing, biopolymer production and leather and wool treatment [48].

2.3.2. Relevance in Colorant, Aromas, Flavor, Fragrance, Cosmetic, and Hygienic Industries

Another group of polymer-degrading hydrolases, for which homologues are used in the fine chemicals industry, comprises glycosidases, alpha-amylase, polygalacturonase, beta-glucosidase, 1,4-alpha-glucosidase and cellulase (Supplementary Table S5). Moreover, in the industry of fine chemicals, oxidoreductases are widely used biocatalysts, and oxygenases (mono- and dioxygenases and peroxidases) are important enzymes for the modification of terpenoids. In corals, a diverse array of diterpenoids has been reported [49], with activities including anti-inflammatory, antifouling and antifeedant, anti-infective (antimicrobial, antiviral, anti-parasite), anticancer and cytotoxic effects. Thus, different species of marine cnidarian holobionts appear to be interesting sources not only for terpenoids themselves but also for enzymes involved in their biosynthesis. The holo-transcriptome of P. variabilis has revealed some of these predicted enzyme sequences.

2.3.3. Relevance in Agrochemical, Food and Feed Industries

In the P. variabilis holo-transcriptome, many predicted enzymatic activities were also identified that may be relevant to the food industry (Supplementary Table S6). Hydrolases and oxidoreductases have emerged as preferred biocatalysts on an industrial scale for the generation of chirality and enantio (selective) kinetic resolutions of chemicals, especially from renewable sources [50].
From an ecological viewpoint, the expression of chitinases and enzymes related to chitin degradation, including chitodextrinase, in the P. variabilis holobiont is suggestive of a complimentary mechanism of self-defense in zoantharians, equipping the organisms to cope with disease-causing agents. Chitinases have known anti-fungal [51] and algicidal properties, contributing to nutrient cycling in marine environments [49]. In pharmaceutical biotechnology, chitinases are useful for preparing chitosan composites for numerous applications, ranging from biosensors, tissue engineering and drug delivery systems to nanoarrays and lab-on-chip devices [52]. Moreover, the coupled reactions of N-acetylhexosaminidase (a chitinase) and deacetylases can produce, a valuable nutraceutical supplement, d-Glucosamine, with high yields from polymeric chitin in a proof-of-concept environmentally friendly dual-catalysis process [53]. Interestingly, a predicted deacetylase (N-acetyl-d-glucosamine-6-phosphate deacetylase, EC 3.5.1.25) was found in the holo-transcriptome of P. variabilis (Supplementary Table S3).
Again, from an ecological viewpoint, the presence of transcriptional precursors encoding enzymes for the biosynthesis of herbicide intermediates, in the P. variabilis holo-transcriptome, can be seen as a possible way for this anthozoan species to restrain macroalgae overgrowth, as seen in corals [54]. Hence, considering the marine environment in which these enzymes presumably work, the possibility to obtain naturally evolved salt-tolerant biocatalysts for industrial bioprocesses is high.

2.3.4. Relevance in Bioconversion and Biopolymer Synthesis

In the P. variabilis holo-transcriptome, some, but not all, industrially useful putative glycosidases that could potentially be applied to the processing of lignocellulose were found (Supplementary Table S7). Moreover, enzymatic activities involved in the production of precursors of biopolymers, such as recyclable bioplastic, were also observed. Examples of such enzymes include, 3-oxoacyl-ACP reductase (E.C. 1.1.1.100), acetyl-CoA C-acyltransferase (2.3.1.16) and enoyl-CoA hydratase (4.2.1.17), which catalyzes the biosynthesis of polyhydroxyalkanate (PHA) precursors. The predicted enzymes epoxide hydrolase (3.3.2.9), nitrile hydratase (4.2.1.84) and γ-glutamyltransferase (2.3.2.13) are important for bioconversion in the renewable energy industry and for the synthesis of fine chemicals, in green chemistry and bioremediation, as will be discussed later. Finally, the peroxidases identified in this zoantharian holo-transcriptome could be further evaluated for applications in the conversion and biosynthesis of phenol- and vinyl-containing polymers, accordingly to the current use of enzymes of these catalytic classes.

2.3.5. Relevance in the Cleaning and Recovery Industries

The identification of enzymatic activities in a given microbiome is useful for estimating the capacity of the microbes to recover an environmental site. The disclosed enzymatic activities, serve as indicators of the bioremediation potential for a given set of contaminants, and they also indicate potential enzymatic catalysts that may be isolated for downstream processing. In this regard, among the enzymes listed in Supplementary Table S8, a very interesting enzyme is glutathione γ-glutamylcysteinyltransferase (E.C. 2.3.15.2), also known as phytochelatin synthase, which catalyzes the synthesis of phytochelatins. Phytochelatins are cysteine-rich peptides responsible for the chelation and sequestration of essential microelements (e.g., copper and zinc) and toxic heavy metals (e.g., cadmium, lead and mercury). Members of the phytochelatin family are biosynthesized from glutathione and are structurally characterized by n repetitions of the γ-GluCys dipeptide followed by a terminal Gly, i.e., (γ-GluCys)n-Gly [55]. The genes encoding enzymes for phytochelatin synthesis have a recognized wide phylogenetic distribution, indicating the importance of coping with metal and metalloid (arsenic, selenium and silver) detoxification across species [55,56]. Thus, the expression of phytochelatin synthase in the P. variabilis holo-transcriptome could be related to detoxification of heavy metals and metal homeostasis in zoantharians. From the biotechnological point of view, immobilized phytochelatin synthase has been utilized to prepare highly stable cadmium-sulfur (CdS) nanocrystals of tunable sizes with optical and electrical properties [57]. Prepared in this way, nanocrystals are useful in the fabrication of semiconductor quantum dots (QDs) and for application as fluorescent tags in biological systems for molecular imaging. Another potential technological use for phytochelatin synthases is related to bioremediation; engineered bacteria have been designed that overexpress a mutated enzyme in combination with a cadmium protein transporter, resulting in Cd accumulation [58].
Last but not least, in this group of enzymes relevant to the recovery and processing industry, epoxide hydrolase and nitrile hydratase homologous were found. Epoxide hydrolases (EC 3.3.2.3) are cofactor-independent biocatalysts that convert epoxides to the corresponding diols, and epoxide hydrolases of microbial origin are useful for producing enantioselective intermediates with vicinal diols in the synthesis of fine chemicals and pharmaceuticals, such as (S)-Ibuprofen, (R)-Mevanolactone and (R)-Eliprodil [59,60]. Nitrile hydratases are Fe- or Co-type metalloenzymes that convert nitriles (organo-cyanides) into industrially valuable amides, such as acrylamides, from acrylonitriles. Tons of acrylamide are used as coagulators, soil conditioners and additives in the paper industry, as well as adhesives, paint components and agents for petroleum recovery [61,62]. Additionally, wild-type and engineered nitrile hydratases that act on myriad cyanide-containing substrates can be applied in bioremediation, wastewater treatment and even for the development of cyanide biosensors [63]. These examples provide a glimpse of the biotechnologically significant enzymes found in the holo-transcriptome of P. variabilis that may be used as biocatalysts with improved activities and selectivity over those currently in use. Moreover, one can speculate on the eco-physiological role of these enzymes in the cnidarian-holobiont assemblage.

2.3.6. Relevance in Molecular Biology and Analytical Applications

In Supplementary Table S9, homologs of enzymes useful for manipulating nucleic acids found in the P. variabilis holo-transcriptome are listed. Among these, two enzymes—deoxyribonuclease I and N-acetylglucosaminidase—have also been used to treat cystic fibrosis and lysosomal storage disease, respectively.
In recombinant DNA technology, deoxyribonuclease (DNase) I catalyzes the hydrolysis of phosphodiester bonds and cleaves double-stranded (ds) and single-stranded (ss) DNA in a sequence-nonspecific manner. DNase I is used, for instance, to make “nick translations” and in DNase I footprinting—a technique used to study the interaction of ligands (proteins and drugs) with DNA [64]. Other predicted endonucleases found in the P. variabilis holo-transcriptome include DNase II and IV; DNase II is a lysosomal “acidic” DNase that preferentially cleaves ssDNA, whereas DNase IV (flap endonuclease-1) is a multifunctional enzyme that cleaves 5′-ssDNA flaps of DNA or RNA. DNase IV has been investigated as a marker of disease risk, since it is involved in DNA metabolism, genomic stability and apoptosis [65].
Several types of predicted ribonucleases were also found, such as ribonucleases (RNases) H and T (Supplementary Table S9). RNase H is a monomeric enzyme that degrades RNA in RNA:DNA heteroduplexes and is useful for the synthesis of complementary DNA (cDNA); RNase III works in multicomponent assemblies to bind and cleave dsRNAs, especially in the processing of dsRNA maturation and the decay of coding and noncoding RNAs, such as miRNAs and siRNAs [66].
DNA-dependent DNA polymerases catalyze the biosynthesis of polydeoxyribonucleotides; diverse applications for DNA polymerases exist, including nucleotide sequencing, in vitro synthesis of the second cDNA strand, DNA amplification and the preparation of DNA hybridization probes [67]. RNA polymerases (RNA-dependent DNA polymerases) are enzymes that transcribe gene sequences into the corresponding RNAs. RNA polymerases are used to prepare hybridization probes and in vitro-transcribed RNA to direct a high-level of expression of cloned genes, as well as to produce capture reagents for RNA-binding proteins and as antisense probes [68]. In the P. variabilis zoantharian holo-transcriptome, predicted and structurally conserved RNA-dependent RNA polymerase (reverse transcriptase) sequences were retrieved. Because reverse transcriptase catalyzes the synthesis of DNA from an RNA template, the uses for this type of enzyme include the preparation of cDNA and the molecular design of inhibitors of retrovirus replication [69]. DNA and RNA ligases catalyze the formation of 3′ → 5′ phosphodiester bonds in nucleic acid (DNA and RNA) molecules. These enzymes have a range of uses in recombinant DNA technology, from elongation or circularization of dsDNA, in the case of DNA ligases, to 3′-labeling of RNA, DNA 5′-tailing of DNA and production of elongated molecules, e.g., in cDNA cloning, in the case of RNA ligases [70]. Transcripts encoding both types of ligases were found to be expressed in the zoantharian holo-transcriptome in this study.

2.4. Prediction of Enzymes with Two or More Activities

In searching for different types of enhanced and unique marine biocatalysts, we wondered whether the P. variabilis holobiont transcriptome contains enzymes with multiple activities. A widely known example of a dual catalytic enzyme is RuBisCO (ribulose-1,5-bisphosphate carboxylase/oxygenase, EC 4.1.1.39), which catalyzes both CO2 fixation in the dark phase of photosynthesis and carbon oxidation in the process of photorespiration [71]. Another interestingly example of a dual catalytic enzyme from plants is the enzyme hydroxycinnamoyl-Coenzyme A:quinate hydroxycinnamoyl transferase (HQT, EC 2.3.1.99), which catalyzes the formation of the strong phenolic antioxidants, chlorogenic and dicaffeoylquinic acids, which are useful as phytonutrients in foods and as pharmaceuticals. It has been demonstrated that in addition to the transesterification of caffeoyl-CoA with quinic acid to produce chlorogenic acid, HQT can form dicaffeoylquinic acid via its chlorogenate:chlorogenate transferase activity [72]. In biotechnology, multi-enzyme systems have been designed with different assembly strategies to mimic natural enzyme complexes and pathways, with the aim of improving catalytic efficiency [73]. We devised an annotation iterative process that resulted in the prediction of 13 putative enzymes with two to three enzymatic activities (Table 1). Six of these enzymes have activities that are partially identified and are related to more than one metabolic pathway (CL12403.contig2, CL2444.contig1, unigene12818, unigene14615, unigene32504, as well as unigene33780), and another four enzymes with dual activities were positioned in a single metabolic pathway. Enzymes that were predicted to possess dual catalytic activity are also shown in Supplementary Figure S5; these were found as a result of our stringent analysis. Some putative dual-activity enzymes might have been missed; however, this strategy proved to be effective, as it included an initial convenient and rapid data mining and screening approach of bi-functional biocatalysts from this zoantharian holo-transcriptome.

2.5. Annotation of Novel Predicted Enzyme Sequences with Partial EC Number

In addition to the prediction of enzymes with activities that are fully characterized, groups of enzymes that are of general interest comprise expressed sequences that are completely new. These, as analyzed herein, could not be mapped to a specific and detailed catalytic reaction, i.e., they comprise enzymatic precursors with an EC number lacking the fourth categorization numbers. Indeed, several predicted enzymes were found with a sequence similarity that were reasonably close to be classified within a known sub-subclass (EC with at least three numbers), but, still, distinct enough to be completely mapped in given reactional group of characterized enzymes. Predicted precursors that fit in this category could hypothetically point to isozymes with already described reactions but working on different substrates, with distinct kinetic parameters, in distinct catalytic conditions or, even, comprising a totally novel catalytic reaction. Based on these facts, from the P. variabilis holo-transcriptome, additional 694 predicted enzyme sequences were found with incomplete EC numbers (unknown fourth serial digit), distributed into 23 sub-subclasses (Table 2), with a large majority related to the class of hydrolases. In this study, based on these preceding findings, we focused our further analysis on two sub-subclasses that P. variabilis sequences are grouped in distinct new clades, representing new structures, namely, cysteine dioxygenases (EC:1.13.11.20) and carboxypeptidases A, B, A2, and U (ECs:3.4.17.1, 3.4.17.2, 3.4.17.15, 3.4.17.20). I Initially, all predicted P. variabilis sequences mapped to a sub-subclass were evaluated by phylogenetic inference with corresponding counterparts representing each species for each enzyme in this sub-subclass. It was found that eleven P. variabilis sequences were related to the cysteine dioxygenases while forming an out-clade, and five sequences with the carboxypeptidases A, B, A2, and U forming a distinct clade all together. Afterwards, a second tree was inferred to confirm that, even if they are related, the P. variabilis sequences form a distinct clade, showing their uniqueness (Figure 2 and Figure 3). Cysteine dioxygenase is a key enzyme in the synthesis of taurine, an important compound, product of cysteine metabolism, that is used in functional foods, as well as in pharmaceutical and cosmetic industries. Several patents are granted to the production of taurine by fermentation methods with transgenic microorganisms. Carboxypeptidases are generally applied either in research or in the pharmaceutical industries, but also found some application in food industries, as exemplified by the use of the carboxylase A in baking industry [74].

3. Materials and Methods

3.1. Biological Sample

Specimens of Protopalythoa variabilis (Duerden, 1898) were collected in coastal reefs of Porto de Galinhas, Pernambuco, Brazil (8°30′20″ S, 35°00′34″ W) during low tide. A voucher specimen (GPA 181) was identified by us and was kept at the cnidarian collection of the Anthozoan Research Group (GPA) at the Academic Center of Vitória, Federal University of Pernambuco (Brazil). This species has been mentioned by some authors as Palythoa variabilis, in reason of a proposed synonymy for the genera Palythoa and Protopalythoa [75]. However, the issue of distinctive genera is not completely solved yet (see, for example, [76]), despite new molecular phylogenetic approach, based on the universal target-enrichment baits, has been recently developed to help resolve long-standing controversial relationships in the class Anthozoa [77]. Thus, until an extensive revision of the group is not definitively resolved with morphological and molecular data precisely combined, with inclusion of species of both genera, the binomial nomenclature Protopalythoa variabilis is used herein, as for decades.

3.2. RNA Library Construction and Origin of Zoantharian RNA Sequences

The RNA isolation, library preparation and the transcriptome assembly of the P. variabilis holobiont were performed as described in one of our previous articles [24]. RNA sequencing was performed using the Illumina HiSeq 2500 platform. Reads were cleaned up before the de novo transcriptome assembly using the Trinity method for transcriptome reconstruction [78]. This Transcriptome Shotgun Assembly (TSA) project was deposited in DDBJ/EMBL/GenBank under the accession GCVI00000000, associated with the BioProject PRJNA279783 and biosample SAMN03450566. The statistics of the RNA sequencing and contig assembling are summarized in the Supplementary Table S1.

3.3. Assessment of the Biodiversity Composition of the P. variabilis Holobiont

P. variabilis transcripts with a high similarity to 16S rRNA, COI, and rbcL sequences were identified using BLASTn with an E-value limit of 1E-40. The closest related species were characterized by homology search of these transcripts against the NCBI nr database. Only the first hits in concordance with the selected gene markers aforementioned were retained for species composition identification in the holobiont assemblage.

3.4. Sequence Annotations for Enzyme Precursors in the P. variabilis Holobiont

The unigenes from the P. variabilis transcriptome were investigated for structural enzyme homology using BLASTx (BLAST+ suite, version 2.5.0) [79] with a fixed E-value of 1E-5 against four public protein databases: the NCBI non-redundant (nr) protein database (accessed from October to November 2016), the Clusters of Orthologous Groups (COG) database, version 2003–2014, the UniProtKB/Swiss-Prot database, downloaded in October 2016), and the EuKaryotic Orthologous Groups (KOG) database, version 2003. Basic statistics of the sequence annotations are described in Supplementary Figures S1 and S2. Species information of the selected annotations was extracted from the BLAST output files to discern the taxonomic distribution.

3.5. Gene Ontology, Enzyme Codes and KEGG Pathway Assignments

The Blast2GO software, version 4.0.2 [80], was used for the subsequent steps under default parameters to carry out the InterProScan protein domain analysis, followed by the Gene Ontology (GO) annotation. The annotations were then subjected to a generic GO-slim reduction, prior to the mapping of the Enzyme Commission codes (EC) and KEGG pathways. The GO annotations chart was plotted using WEGO [81], while further KEGG pathways information was retrieved through the KEGG BRITE hierarchies site [82].

3.6. Sequence Alignment and Phylogenetic Inference

Multiple sequence alignments of predicted enzyme sequences were performed using MUSCLE v3.8.31 [83] or Kalign v2.9.0b2 [84], depending on whether the enzyme domain/motif was a single structural unit or multiple-domain repeats, respectively. The phylogenetic tree was inferred using MEGA 7.0.26 [85] based on the LG+G+I model, with a bootstrap test of 500 replicates and edited using FigTree v 1.4.3. The amino acid sequence identities and similarities were determined using Jalview 2.9.0b2 [86].

3.7. Prediction of Enzymes with Two or More Activities

An alternative annotation strategy was used to identify potential enzymes with two or more enzymatic activities. To favor the multiple domains (hits) discovery, we first performed a BLASTx search against the UniProtKB/Swiss-Prot database, keeping an E-value of 1E-5 with the additional “culling_limit” option set to “1”. The selected set was reduced to only sequences with at least two non-overlapping hits on the same frame. During a second round of selection, only sequences with predicted protein products containing regions found by BLAST were kept. Finally, an InterProScan protein domain analysis [87] was performed to map each sequence to its GO annotation and its associated EC codes and KEGG pathways using the external EC2GO [88] and KEGG2GO mapping databases [89,90]. Only sequences with two or more ECs were kept, and data on the enzyme’s substrate(s) and product(s) were retrieved from KEGG Enzyme [82].

3.8. Analysis of Predicted Enzyme with Partial EC Number

Results of the annotation generated during the Blast2GO analysis were used to extract a list of the transcripts with a partial EC number. For each of these sequences, the corresponding protein sequence is predicted and subsequently submitted to InterProScan 5 via RESTful service 80 to validate the presence of a catalytic domain corresponding to the partial EC. These steps were automatized using a python script, depending on the external EC2GO mapping database [88]. To perform the observed distance analysis between the predicted P. variabilis and the known enzyme sequences, the available collections of sequences corresponding to each EC of interest were retrieved from the BRENDA website (www.brenda-enzymes.org, release 2018.1) [91] and sequences with a length shorter or higher to two population standard deviation were removed. The multiple sequence alignment was done using Kalign v2.04 77, and the analysis was computed with MEGA 7.0.26 78 using the Neighbor-Joining method and the p-distance method, removing missing ambiguous positions for each sequence pair, and with a bootstrap test of 500 replicates. The cladograms were edited in the Interactive Tree Of Life website [92].

4. Conclusions

Marine invertebrates with associated microbiota form complex holobiont assemblages, which are attractive sources of biologically active organic compounds and (poly-)peptides, including enzymes. Marine enzymes have a high potential to be applied in green organic synthesis and in pharmaceutical and industrial biotechnology. The search for improved biocatalysts can be carried out using different strategies, such as screening a huge number of environmental samples, pursuing enzyme engineering, mining genomic and proteomic data, or a combination of more than one approach. The process of data mining transcriptomes has some advantages over genomic analysis; the most obvious advantage is that only enzymes that are expressed in a given environmental context are retrieved, including enzymes not completely characterized and with unknow enzyme-substrate specificity. This is particularly advantageous in the case of the marine assemblages of microbionts that form species-specific holobionts, from which the purification of enzymes with a high yield may be a concern. Thus, once identified, the cloning and the recombinant production of desirable marine biocatalysts can be structure-guided and based on the nature of expressed transcripts.
According to data reported in this work, the zoantharian P. variabilis expresses a variety of putative enzymes that could potentially be converted into biotechnologically useful biocatalysts and biopharmaceuticals. This holo-transcriptomic data demonstrates that a single holobiont assemblage comprises a unique repository of relevant biotechnological enzymes. Finally, the integrative analyses of this holo-transcriptome point to a valuable marine resource for the discovery of improved enzymes with applications in green chemistry, industrial and pharmaceutical biotechnology.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-3397/16/6/207/s1. Figure S1: Summary of the unigenes mapped to the public databases. Figure S2: Characteristics of the sequence homology search results. Figure S3: GO and KEGG pathways assignments. Figure S4: Maximum Likelihood (ML) phylogenetic tree of the predicted P. variabilis beta-N-acetylhexosaminidase (3.2.1.52) and their closest homologous sequences. Figure S5: Prediction of enzymes with two activities closely positioned in given metabolic pathway. Table S1: Species from the BLAST search of the barcode sequences from the P. variabilis holobiont transcriptome against the NCBI nr database. Table S2: KEGG pathways mapping summary. Table S3: List of predicted enzymes in P. variabilis. Table S4: List of enzymatic activities with relevance in treatment of rare diseases and other pharmaceutical fine chemicals predicted in Protopalythoa variabilis holo-transcriptome. Table S5: List of enzymatic activities with relevance in colorant, aromas, flavor, fragrance, cosmetic and hygienic industries predicted in Protopalythoa variabilis holo-transcriptome. Table S6: List of enzymatic activities with relevance in agrochemical, and food and feed industries predicted in Protopalythoa variabilis holo-transcriptome. Table S7: List of enzymatic activities with relevance in bioconversion and biopolymer synthesis predicted in Protopalythoa variabilis holo-transcriptome. Table S8: List of enzymatic activities with relevance in other industries predicted in Protopalythoa variabilis holo-transcriptome. Table S9: List of enzymatic activities with relevance in molecular biology and analytical applications predicted in Protopalythoa variabilis holo-transcriptome. Table S10: RNA-sequencing and assembling statistics. File S1: KEGG pathways.

Author Contributions

Conceived and designed the experiments: J.-É.R.L.M. and G.R.-B. Performed the experiments: J.-É.R.L.M., C.H. and Q.L. Analyzed the data: J.-É.R.L.M., A.R.B.P.-S. and G.R.-B. Contributed reagents, materials and/or analysis tools: S.M.-Y.L. and G.R.-B. Field trip missions, species identification and deposit collection: P.B.G. and C.D.P. wrote the manuscript with the contribution from all other authors: J.-É.R.L.M. and G.R.-B. All authors have read and approved the final article.

Funding

Financial support for G.R.-B. and research at the Institute for Marine Sciences, Federal University of Ceará, was provided by the Brazilian National Council for Scientific and Technological Development, CNPq (Marine Biotechnology Network—Proc. 408835/2013-3), the Ministry of Science, Technology and Innovation (MCTI), Brasília, DF, Brazil. At the University Federal of Pernambuco, financial support was CNPq (Proc. 408934/2013-1) granted to C.D.P. Research at the University of Macau was supported by grants from the Science and Technology Development Fund (FDCT) of Macao SAR and Research Committee, University of Macau, Macau, China.

Acknowledgments

Our gratitude to the Program on Toxinology (Issue 2010), the Coordination for the Improvement of Higher Education Personnel (CAPES), the Ministry of Education of the Federal Government of Brazil, Brasília, DF, Brazil. J-ÉRM was a doctoral fellowship recipient from CAPES.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Enzymes Market by Type (Industrial, Specialty), by Product (Carbohydrases, Proteases, Lipases, Polymerases & Nucleases), by Application (Food & Beverages, Detergents, Animal Feed, Textile, Paper & Pulp, Nutraceutical, Personal Care & Cosmetics, Wastewater, Research & Biotechnology, Diagnostics, Biocatalyst) and Segment Forecasts to 2024. Available online: https://www.grandviewresearch.com/industry-analysis/enzymes-industry (accessed on 13 June 2018).
  2. Molecular Biology Enzymes and Kits & Reagents Market Analysis by Product (Kits and Reagents, Enzymes), by Application (PCR, Sequencing, Cloning, Epigenetics), by End User, and Segment Forecasts, 2018–2025. Available online: https://www.grandviewresearch.com/industry-analysis/molecular-biology-enzymes-kits-reagents-market (accessed on 13 June 2018).
  3. Itoh, T.; Hanefeld, U. Enzyme catalysis in organic synthesis. Green Chem. 2017, 19, 331–332. [Google Scholar] [CrossRef]
  4. Wenda, S.; Illner, S.; Mell, A.; Kragl, U. Industrial biotechnology-the future of green chemistry? Green Chem. 2011, 13, 3007–3047. [Google Scholar] [CrossRef]
  5. Yamaguchi, S. The quest for industrial enzymes from microorganisms. Biosci. Biotechnol. Biochem. 2017, 81, 54–58. [Google Scholar] [CrossRef] [PubMed]
  6. Shimizu, S.; Ogawa, J.; Kataoka, M.; Kobayashi, M. Screening of novel microbial enzymes for the production of biologically and chemically useful compounds. Adv. Biochem. Eng. Biotechnol. 1997, 58, 45–87. [Google Scholar] [PubMed]
  7. Wang, J.B.; Li, G.; Reetz, M.T. Enzymatic site-selectivity enabled by structure-guided directed evolution. Chem. Commun. 2017, 53, 3916–3928. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Hibbert, E.G.; Dalby, P.A. Directed evolution strategies for improved enzymatic performance. Microb. Cell Fact. 2005, 4, 29. [Google Scholar] [CrossRef] [PubMed]
  9. Madhavan, A.; Sindhu, R.; Parameswaran, B.; Sukumaran, R.K.; Pandey, A. Metagenome Analysis: A Powerful Tool for Enzyme Bioprospecting. Appl. Biochem. Biotechnol. 2017, 183, 636–651. [Google Scholar] [CrossRef] [PubMed]
  10. Trincone, A. Marine biocatalysts: Enzymatic features and applications. Mar. Drugs 2011, 9, 478–499. [Google Scholar] [CrossRef] [PubMed]
  11. Smith, D.R.M.; Uria, A.R.; Helfrich, E.J.N.; Milbredt, D.; van Pee, K.H.; Piel, J.; Goss, R.J.M. An Unusual Flavin-Dependent Halogenase from the Metagenome of the Marine Sponge Theonella swinhoei WA. ACS Chem. Biol. 2017, 12, 1281–1287. [Google Scholar] [CrossRef] [PubMed]
  12. Sarian, F.D.; Janecek, S.; Pijning, T.; Ihsanawati; Nurachman, Z.; Radjasa, O.K.; Dijkhuizen, L.; Natalia, D.; van der Maarel, M.J. A new group of glycoside hydrolase family 13 alpha-amylases with an aberrant catalytic triad. Sci. Rep. 2017, 7, 44230. [Google Scholar] [CrossRef] [PubMed]
  13. Lima, R.N.; Porto, A.L. Recent Advances in Marine Enzymes for Biotechnological Processes. Adv. Food Nutr. Res. 2016, 78, 153–192. [Google Scholar] [PubMed]
  14. Trincone, A. Enzymatic Processes in Marine Biotechnology. Mar. Drugs 2017, 15, 93. [Google Scholar] [CrossRef] [PubMed]
  15. Leal, M.C.; Puga, J.; Serodio, J.; Gomes, N.C.; Calado, R. Trends in the discovery of new marine natural products from invertebrates over the last two decades—Where and what are we bioprospecting? PLoS ONE 2012, 7, e30580. [Google Scholar] [CrossRef] [PubMed]
  16. Lindequist, U. Marine-Derived Pharmaceuticals—Challenges and Opportunities. Biomol. Ther. 2016, 24, 561–571. [Google Scholar] [CrossRef] [PubMed]
  17. Ryu, T.; Seridi, L.; Moitinho-Silva, L.; Oates, M.; Liew, Y.J.; Mavromatis, C.; Wang, X.; Haywood, A.; Lafi, F.F.; Kupresanin, M.; et al. Hologenome analysis of two marine sponges with different microbiomes. BMC Genomics 2016, 17, 158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Macintyre, L.; Zhang, T.; Viegelmann, C.; Martinez, I.J.; Cheng, C.; Dowdells, C.; Abdelmohsen, U.R.; Gernert, C.; Hentschel, U.; Edrada-Ebel, R. Metabolomic tools for secondary metabolite discovery from marine microbial symbionts. Mar. Drugs 2014, 12, 3416–3448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Dunlap, W.C.; Battershill, C.N.; Liptrot, C.H.; Cobb, R.E.; Bourne, D.G.; Jaspars, M.; Long, P.F.; Newman, D.J. Biomedicinals from the phytosymbionts of marine invertebrates: A molecular approach. Methods 2007, 42, 358–376. [Google Scholar] [CrossRef] [PubMed]
  20. Huang, C.Y.; Sung, P.J.; Uvarani, C.; Su, J.H.; Lu, M.C.; Hwang, T.L.; Dai, C.F.; Wu, S.L.; Sheu, J.H. Glaucumolides A and B, Biscembranoids with New Structural Type from a Cultured Soft Coral Sarcophyton glaucum. Sci. Rep. 2015, 5, 15624. [Google Scholar] [CrossRef] [PubMed]
  21. Huang, C.Y.; Su, J.H.; Liaw, C.C.; Sung, P.J.; Chiang, P.L.; Hwang, T.L.; Dai, C.F.; Sheu, J.H. Bioactive Steroids with Methyl Ester Group in the Side Chain from a Reef Soft Coral Sinularia brassica Cultured in a Tank. Mar. Drugs 2017, 15, 280. [Google Scholar] [CrossRef] [PubMed]
  22. Hsiao, T.H.; Sung, C.S.; Lan, Y.H.; Wang, Y.C.; Lu, M.C.; Wen, Z.H.; Wu, Y.C.; Sung, P.J. New Anti-Inflammatory Cembranes from the Cultured Soft Coral Nephthea columnaris. Mar. Drugs 2015, 13, 3443–3453. [Google Scholar] [CrossRef] [PubMed]
  23. Ferrer, M.; Martinez-Martinez, M.; Bargiela, R.; Streit, W.R.; Golyshina, O.V.; Golyshin, P.N. Estimating the success of enzyme bioprospecting through metagenomics: Current status and future trends. Microb. Biotechnol. 2016, 9, 22–34. [Google Scholar] [CrossRef] [PubMed]
  24. Huang, C.; Morlighem, J.R.; Zhou, H.; Lima, E.P.; Gomes, P.B.; Cai, J.; Lou, I.; Perez, C.D.; Lee, S.M.; Radis-Baptista, G. The Transcriptome of the Zoanthid Protopalythoa variabilis (Cnidaria, Anthozoa) Predicts a Basal Repertoire of Toxin-like and Venom-Auxiliary Polypeptides. Genome Biol. Evol. 2016, 8, 3045–3064. [Google Scholar] [CrossRef] [PubMed]
  25. Rabelo, E.F.; Soares, M.D.O.; Bezerra, L.E.A.; Matthews-Cascon, H. Distribution pattern of zoanthids (Cnidaria: Zoantharia) on a tropical reef. Mar. Biol. Res. 2015, 11, 584–592. [Google Scholar] [CrossRef]
  26. Rabelo, E.; Rocha, L.; Colares, G.; Araújo Bomfim, T.; Lúcia Rodrigues Nogueira, V.; Katzenberger, M.; Matthews-Cascon, H.; Melo, V. Symbiodinium Diversity Associated with Zoanthids (Cnidaria: Hexacorallia) in Northeastern Brazil; Springer: Dordrecht, The Netherlands, 2015. [Google Scholar]
  27. Sun, W.; Zhang, F.; He, L.; Li, Z. Pyrosequencing reveals diverse microbial community associated with the zoanthid Palythoa australiae from the South China Sea. Microb. Ecol. 2014, 67, 942–950. [Google Scholar] [CrossRef] [PubMed]
  28. Chu, N.D.; Vollmer, S.V. Caribbean corals house shared and host-specific microbial symbionts over time and space. Environ. Microbiol. Rep. 2016, 8, 493–500. [Google Scholar] [CrossRef] [PubMed]
  29. Newman, D.J.; Cragg, G.M.; Snader, K.M. The influence of natural products upon drug discovery. Nat. Prod. Rep. 2000, 17, 215–234. [Google Scholar] [CrossRef] [PubMed]
  30. Urbarova, I.; Karlsen, B.O.; Okkenhaug, S.; Seternes, O.M.; Johansen, S.D.; Emblem, A. Digital marine bioprospecting: Mining new neurotoxin drug candidates from the transcriptomes of cold-water sea anemones. Mar. Drugs 2012, 10, 2265–2279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Anderson, D.A.; Walz, M.E.; Weil, E.; Tonellato, P.; Smith, M.C. RNA-Seq of the Caribbean reef-building coral Orbicella faveolata (Scleractinia-Merulinidae) under bleaching and disease stress expands models of coral innate immunity. PeerJ 2016, 4, e1616. [Google Scholar] [CrossRef] [PubMed]
  32. Ponce, D.; Brinkman, D.L.; Potriquet, J.; Mulvenna, J. Tentacle Transcriptome and Venom Proteome of the Pacific Sea Nettle, Chrysaora fuscescens (Cnidaria: Scyphozoa). Toxins 2016, 8, 102. [Google Scholar] [CrossRef] [PubMed]
  33. Hasegawa, Y.; Watanabe, T.; Takazawa, M.; Ohara, O.; Kubota, S. De Novo Assembly of the Transcriptome of Turritopsis, a Jellyfish that Repeatedly Rejuvenates. Zool. Sci. 2016, 33, 366–371. [Google Scholar] [CrossRef] [PubMed]
  34. Valmsen, K.; Jarving, I.; Boeglin, W.E.; Varvas, K.; Koljak, R.; Pehk, T.; Brash, A.R.; Samel, N. The origin of 15R-prostaglandins in the Caribbean coral Plexaura homomalla: Molecular cloning and expression of a novel cyclooxygenase. Proc. Natl. Acad. Sci. USA 2001, 98, 7700–7705. [Google Scholar] [CrossRef] [PubMed]
  35. Uria, A.R.; Zilda, D.S. Metagenomics-Guided Mining of Commercially Useful Biocatalysts from Marine Microorganisms. Adv. Food Nutr. Res. 2016, 78, 1–26. [Google Scholar] [PubMed]
  36. Popovic, A.; Tchigvintsev, A.; Tran, H.; Chernikova, T.N.; Golyshina, O.V.; Yakimov, M.M.; Golyshin, P.N.; Yakunin, A.F. Metagenomics as a Tool for Enzyme Discovery: Hydrolytic Enzymes from Marine-Related Metagenomes. Adv. Exp. Med. Biol. 2015, 883, 1–20. [Google Scholar] [PubMed]
  37. Sukul, P.; Schakermann, S.; Bandow, J.E.; Kusnezowa, A.; Nowrousian, M.; Leichert, L.I. Simple discovery of bacterial biocatalysts from environmental samples through functional metaproteomics. Microbiome 2017, 5, 28. [Google Scholar] [CrossRef] [PubMed]
  38. Sturmberger, L.; Wallace, P.W.; Glieder, A.; Birner-Gruenberger, R. Synergism of proteomics and mRNA sequencing for enzyme discovery. J. Biotechnol. 2016, 235, 132–138. [Google Scholar] [CrossRef] [PubMed]
  39. Guo, F.; Berglund, P. Transaminase biocatalysis: Optimization and application. Green Chem. 2017, 19, 333–360. [Google Scholar] [CrossRef]
  40. McGrath, B.M.; Walsh, G. Directory of Therapeutic Enzymes; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2005. [Google Scholar]
  41. Sun, Z.; Ilie, A.; Reetz, M.T. Towards the Production of Universal Blood by Structure-Guided Directed Evolution of Glycoside Hydrolases. Angew. Chem. 2015, 54, 9158–9160. [Google Scholar] [CrossRef] [PubMed]
  42. Greiner-Tollersrud, O.K.; Berg, T. Lysosomal Storage Disorders. In Madame Curie Bioscience Database [Internet]; Landes Bioscience: Austin, TX, USA, 2013. [Google Scholar]
  43. Desnick, R.J.; Schuchman, E.H. Enzyme replacement therapy for lysosomal diseases: Lessons from 20 years of experience and remaining challenges. Annu. Rev. Genomics Hum. Genet. 2012, 13, 307–335. [Google Scholar] [CrossRef] [PubMed]
  44. Kloosterman, W.M.J.; Roest, S.; Priatna, S.R.; Stavila, E.; Loos, K. Chemo-enzymatic synthesis route to poly(glucosyl-acrylates) using glucosidase from almonds. Green Chem. 2014, 16, 1837–1846. [Google Scholar] [CrossRef]
  45. Vellard, M. The enzyme as drug: Application of enzymes as pharmaceuticals. Curr. Opin. Biotechnol. 2003, 14, 444–450. [Google Scholar] [CrossRef]
  46. Olsson, M.L.; Clausen, H. Modifying the red cell surface: Towards an ABO-universal blood supply. Br. J. Haematol. 2008, 140, 3–12. [Google Scholar] [CrossRef] [PubMed]
  47. Dozier, J.K.; Distefano, M.D. Site-Specific PEGylation of Therapeutic Proteins. Int. J. Mol. Sci. 2015, 16, 25831–25864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Zhang, D.; Zhu, Y.; Chen, J. Microbial transglutaminase production: Understanding the mechanism. Biotechnol. Genet. Eng. Rev. 2010, 26, 205–222. [Google Scholar] [CrossRef] [PubMed]
  49. Li, Y.; Lei, X.; Zhu, H.; Zhang, H.; Guan, C.; Chen, Z.; Zheng, W.; Fu, L.; Zheng, T. Chitinase producing bacteria with direct algicidal activity on marine diatoms. Sci. Rep. 2016, 6, 21984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Hollmann, F.; Arends, I.W.C.E.; Holtmann, D. Enzymatic reductions for the chemist. Green Chem. 2011, 13, 2285–2314. [Google Scholar] [CrossRef]
  51. Han, Y.; Yang, B.; Zhang, F.; Miao, X.; Li, Z. Characterization of antifungal chitinase from marine Streptomyces sp. DA11 associated with South China Sea sponge Craniella australiensis. Mar. Biotechnol. 2009, 11, 132–140. [Google Scholar] [CrossRef] [PubMed]
  52. Hein, S.; Wang, K.; Stevens, W.F.; Kjems, J. Chitosan composites for biomedical applications: Status, challenges and perspectives. Mater. Sci. Technol. 2008, 24, 1053–1061. [Google Scholar] [CrossRef]
  53. Lyu, Y.; Laborda, P.; Huang, K.; Cai, Z.P.; Wang, M.; Aimin, L.; Doherty, C.; Liu, L.; Flitsch, S.; Voglmeir, J. Highly Efficient and Selective Biocatalytic Production of Glucosamine from Chitin. Green Chem. 2016, 19, 527–535. [Google Scholar]
  54. Shearer, T.L.; Snell, T.W.; Hay, M.E. Gene expression of corals in response to macroalgal competitors. PLoS ONE 2014, 9, e114525. [Google Scholar] [CrossRef] [PubMed]
  55. Cobbett, C.; Goldsbrough, P. Phytochelatins and metallothioneins: Roles in heavy metal detoxification and homeostasis. Annu. Rev. Plant Biol. 2002, 53, 159–182. [Google Scholar] [CrossRef] [PubMed]
  56. Bundy, J.G.; Kille, P. Metabolites and metals in Metazoa—What role do phytochelatins play in animals? Metallomics 2014, 6, 1576–1582. [Google Scholar] [CrossRef] [PubMed]
  57. Liu, F.; Kang, S.H.; Lee, Y.-I.; Choa, Y.-H.; Mulchandani, A.; Myung, N.V.; Chen, W. Enzyme mediated synthesis of phytochelatin-capped CdS nanocrystals. Appl. Phys. Lett. 2010, 97, 123703. [Google Scholar] [CrossRef]
  58. Kang, S.H.; Singh, S.; Kim, J.Y.; Lee, W.; Mulchandani, A.; Chen, W. Bacteria metabolically engineered for enhanced phytochelatin production and cadmium accumulation. Appl. Environ. Microbiol. 2007, 73, 6317–6320. [Google Scholar] [CrossRef] [PubMed]
  59. Archelas, A.; Furstoss, R. Synthetic applications of epoxide hydrolases. Curr. Opin. Chem. Biol. 2001, 5, 112–119. [Google Scholar] [CrossRef]
  60. Van Loo, B.; Kingma, J.; Arand, M.; Wubbolts, M.G.; Janssen, D.B. Diversity and biocatalytic potential of epoxide hydrolases identified by genome analysis. Appl. Environ. Microbiol. 2006, 72, 2905–2917. [Google Scholar] [CrossRef] [PubMed]
  61. Yamada, H.; Shimizu, S.; Kobayashi, M. Hydratases involved in nitrile conversion: Screening, characterization and application. Chem. Rec. 2001, 1, 152–161. [Google Scholar] [CrossRef] [PubMed]
  62. Yamada, H.; Kobayashi, M. Nitrile hydratase and its application to industrial production of acrylamide. Biosci. Biotechnol. Biochem. 1996, 60, 1391–1400. [Google Scholar] [CrossRef] [PubMed]
  63. Park, J.M.; Trevor Sewell, B.; Benedik, M.J. Cyanide bioremediation: The potential of engineered nitrilases. Appl. Microbiol. Biotechnol. 2017, 101, 3029–3042. [Google Scholar] [CrossRef] [PubMed]
  64. Eun, H.-M. 3-Nucleases. In Enzymology Primer for Recombinant DNA Technology; Academic Press: San Diego, CA, USA, 1996; pp. 145–232. [Google Scholar]
  65. Kim, I.; Hur, N.W.; Shin, H.D.; Park, B.L.; Cheong, H.S.; Bae, S.C. Associations of DNase IV polymorphisms with autoantibodies in patients with systemic lupus erythematosus. Rheumatology 2008, 47, 996–999. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Nicholson, A.W. Ribonuclease III mechanisms of double-stranded RNA cleavage. Wiley Interdiscip. Rev. RNA 2014, 5, 31–48. [Google Scholar] [CrossRef] [PubMed]
  67. Eun, H.-M. 6-DNA Polymerases. In Enzymology Primer for Recombinant DNA Technology; Academic Press: San Diego, CA, USA, 1996; pp. 345–489. [Google Scholar]
  68. Eun, H.-M. 7-RNA Polymerases. In Enzymology Primer for Recombinant DNA Technology; Academic Press: San Diego, CA, USA, 1996; pp. 491–565. [Google Scholar]
  69. Bachman, J. Reverse-transcription PCR (RT-PCR). Methods Enzymol. 2013, 530, 67–74. [Google Scholar] [PubMed]
  70. Eun, H.-M. 2-Ligases. In Enzymology Primer for Recombinant DNA Technology; Academic Press: San Diego, CA, USA, 1996; pp. 109–144. [Google Scholar]
  71. Spreitzer, R.J.; Salvucci, M.E. Rubisco: Structure, regulatory interactions, and possibilities for a better enzyme. Annu. Rev. Plant Biol. 2002, 53, 449–475. [Google Scholar] [CrossRef] [PubMed]
  72. Moglia, A.; Lanteri, S.; Comino, C.; Hill, L.; Knevitt, D.; Cagliero, C.; Rubiolo, P.; Bornemann, S.; Martin, C. Dual catalytic activity of hydroxycinnamoyl-coenzyme A quinate transferase from tomato allows it to moonlight in the synthesis of both mono- and dicaffeoylquinic acids. Plant Physiol. 2014, 166, 1777–1787. [Google Scholar] [CrossRef] [PubMed]
  73. Wang, S.Z.; Zhang, Y.H.; Ren, H.; Wang, Y.L.; Jiang, W.; Fang, B.S. Strategies and perspectives of assembling multi-enzyme systems. Crit. Rev. Biotechnol. 2017, 37, 1024–1037. [Google Scholar] [CrossRef] [PubMed]
  74. Kupski, L.; Queiroz, M.I.; Badiale-Furlong, E. Application of carboxypeptidase A to a baking process to mitigate contamination of wheat flour by ochratoxin A. Process Biochem. 2018, 64, 248–254. [Google Scholar] [CrossRef]
  75. Reimer, J.D.; Ono, S.; Takishita, K.; Tsukahara, J.; Maruyama, T. Molecular evidence suggesting species in the zoanthid genera Palythoa and Protopalythoa (Anthozoa: Hexacorallia) are congeneric. Zool. Sci. 2006, 23, 87–94. [Google Scholar] [CrossRef] [PubMed]
  76. Swain, T.D. Revisiting the phylogeny of Zoanthidea (Cnidaria: Anthozoa): Staggered alignment of hypervariable sequences improves species tree inference. Mol. Phyl. Evol. 2018, 118, 1–12. [Google Scholar] [CrossRef] [PubMed]
  77. Quattrini, A.M.; Faircloth, B.C.; Duenas, L.F.; Bridge, T.C.L.; Brugler, M.R.; Calixto-Botia, I.F.; DeLeo, D.M.; Foret, S.; Herrera, S.; Lee, S.M.Y.; et al. Universal target-enrichment baits for anthozoan (Cnidaria) phylogenomics: New approaches to long-standing problems. Mol. Ecol. Resour. 2018, 18, 281–295. [Google Scholar] [CrossRef] [PubMed]
  78. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protocols 2013, 8, 1494–1512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [PubMed]
  80. Conesa, A.; Gotz, S.; Garcia-Gomez, J.M.; Terol, J.; Talon, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [PubMed]
  81. Ye, J.; Fang, L.; Zheng, H.; Zhang, Y.; Chen, J.; Zhang, Z.; Wang, J.; Li, S.; Li, R.; Bolund, L.; et al. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34, W293–W297. [Google Scholar] [CrossRef] [PubMed]
  82. Kanehisa, M.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016, 44, D457–D462. [Google Scholar] [CrossRef] [PubMed]
  83. Edgar, R.C. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004, 5, 113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Lassmann, T.; Frings, O.; Sonnhammer, E.L. Kalign2: High-performance multiple alignment of protein and nucleotide sequences allowing external features. Nucleic Acids Res. 2009, 37, 858–865. [Google Scholar] [CrossRef] [PubMed]
  85. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  86. Waterhouse, A.M.; Procter, J.B.; Martin, D.M.; Clamp, M.; Barton, G.J. Jalview Version 2—A multiple sequence alignment editor and analysis workbench. Bioinformatics 2009, 25, 1189–1191. [Google Scholar] [CrossRef] [PubMed]
  87. Jones, P.; Binns, D.; Chang, H.Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [PubMed]
  88. Hill, D.P.; Davis, A.P.; Richardson, J.E.; Corradi, J.P.; Ringwald, M.; Eppig, J.T.; Blake, J.A. PROGRAM DESCRIPTION: Strategies for Biological Annotation of Mammalian Systems: Implementing Gene Ontologies in Mouse Genome Informatics. Genomics 2001, 74, 121–128. [Google Scholar] [CrossRef] [PubMed]
  89. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene Ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25. [Google Scholar] [CrossRef] [PubMed]
  90. The Gene Ontology Consortium. Expansion of the Gene Ontology knowledgebase and resources. Nucleic Acids Res. 2017, 45, D331–D338. [Google Scholar]
  91. Placzek, S.; Schomburg, I.; Chang, A.; Jeske, L.; Ulbrich, M.; Tillack, J.; Schomburg, D. BRENDA in 2017: New perspectives and new tools in BRENDA. Nucleic Acids Res. 2017, 45, D380–D388. [Google Scholar] [CrossRef] [PubMed]
  92. Letunic, I.; Bork, P. Interactive tree of life (iTOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016, 44, W242–W245. [Google Scholar] [CrossRef] [PubMed]
Marinedrugs 16 00207 g001a 550Marinedrugs 16 00207 g001b 550
Figure 1. Biodiversity composition and taxonomic classification of unigenes from the P. variabilis holo-transcriptome. (A) Number of different species identified using the 16S rRNA, COI, and rbcL barcode sequences. (B) Taxonomic classification of unigenes from the P. variabilis holo-transcriptome after BLASTx analysis. Inserted Box: the distribution of top-hit species in the P. variabilis holobiont with the cnidarian species in dark red, the haptophytes and Stramenopiles in purple, the alveolates in green, and the green algae in gray.
Figure 1. Biodiversity composition and taxonomic classification of unigenes from the P. variabilis holo-transcriptome. (A) Number of different species identified using the 16S rRNA, COI, and rbcL barcode sequences. (B) Taxonomic classification of unigenes from the P. variabilis holo-transcriptome after BLASTx analysis. Inserted Box: the distribution of top-hit species in the P. variabilis holobiont with the cnidarian species in dark red, the haptophytes and Stramenopiles in purple, the alveolates in green, and the green algae in gray.
Marinedrugs 16 00207 g001aMarinedrugs 16 00207 g001b
Marinedrugs 16 00207 g002 550
Figure 2. Cladogram depiction of the P. variabilis predicted sequences related to cysteine dioxygenase (EC:1.13.11.20) members. Tree based on the distance between the protein sequences of cysteine dioxygenase (EC:1.13.11.20) from 106 species and eleven P. variabilis predicted enzymes with an incompletely annotated EC:1.13.11 activity. Twenty catechol 2,3-dioxygenase (EC:1.13.11.2) sequences were used as outgroup. Only bootstrap values greater than 50% are shown at the branch points, in blue, purple, or red color for values comprised between 50–69%, 70–89%, and 90–100% respectively. Enzyme activities are indicated at the name base by circles colored as indicated in the legend.table.
Figure 2. Cladogram depiction of the P. variabilis predicted sequences related to cysteine dioxygenase (EC:1.13.11.20) members. Tree based on the distance between the protein sequences of cysteine dioxygenase (EC:1.13.11.20) from 106 species and eleven P. variabilis predicted enzymes with an incompletely annotated EC:1.13.11 activity. Twenty catechol 2,3-dioxygenase (EC:1.13.11.2) sequences were used as outgroup. Only bootstrap values greater than 50% are shown at the branch points, in blue, purple, or red color for values comprised between 50–69%, 70–89%, and 90–100% respectively. Enzyme activities are indicated at the name base by circles colored as indicated in the legend.table.
Marinedrugs 16 00207 g002
Marinedrugs 16 00207 g003 550
Figure 3. Cladogram depiction of the P. variabilis predicted sequences related to carboxypeptidase A, B, A2, and U (EC:3.4.17.1, 2, 15, and 20) members. Tree based on the distance between the protein sequences of carboxypeptidase A, B, A2, and U (EC:3.4.17.1, 2, 15, and 20) from 40, 19, 4, and 5 species respectively and five P. variabilis predicted enzymes with an incompletely annotated EC:3.4.17 activity. Sequences from 24 carboxypeptidase Taq (EC:3.4.17.19) sequences were used as outgroup. Tree legend same as in Figure 2.
Figure 3. Cladogram depiction of the P. variabilis predicted sequences related to carboxypeptidase A, B, A2, and U (EC:3.4.17.1, 2, 15, and 20) members. Tree based on the distance between the protein sequences of carboxypeptidase A, B, A2, and U (EC:3.4.17.1, 2, 15, and 20) from 40, 19, 4, and 5 species respectively and five P. variabilis predicted enzymes with an incompletely annotated EC:3.4.17 activity. Sequences from 24 carboxypeptidase Taq (EC:3.4.17.19) sequences were used as outgroup. Tree legend same as in Figure 2.
Marinedrugs 16 00207 g003
Table
Table 1. List of predicted multi-functional enzymes from the P. variabilis holo-transcriptome with dual catalytic activities.
Table 1. List of predicted multi-functional enzymes from the P. variabilis holo-transcriptome with dual catalytic activities.
UnigeneECsActivitiesSubstrateProduct
CL12403.Contig11.5.1.20methylenetetrahydrofolate reductase (NAD(P)H)5-methyltetrahydrofolate5,10-methylenetetrahydrofolate
2.1.2.1glycine hydroxymethyltransferase5,10-methylenetetrahydrofolatetetrahydrofolate
CL12403.Contig21.5.1.20methylenetetrahydrofolate reductase (NAD(P)H)5-methyltetrahydrofolate5,10-methylenetetrahydrofolate
2transferase??
CL2444.Contig14.2.1.11phosphopyruvate hydratase2-phospho-d-glyceratephosphoenolpyruvate
5.3.1.1triose-phosphate isomerased-glyceraldehyde 3-phosphateglycerone phosphate
Unigene128181.3.99.1succinate dehydrogenasesuccinatefumarate
1.8.4oxidoreductase, acting on a sulfur group of donors??
5.3.4.1protein disulfide isomerase-S-S- bonds rearrangementnone
Unigene146152.4.2.30NAD ADP-ribosyltransferaseNAD+nicotinamide
2.4.2.31NAD(P)+protein-arginine ADP-ribosyltransferaseNAD+ & protein-l-argininenicotinamide & protein-omega-N-(ADP-d-ribosyl)-l-arginine
Unigene280092.5.1.21farnesyl-diphosphate farnesyltransferase(2E,6E)-farnesyl diphosphatesqualene
1.14oxidoreductase??
Unigene325042.1.1O-methyltransferase??
2.6.1.42branched-chain-amino-acid transaminaseL-leucine4-methyl-2-oxopentanoate
Unigene337802.7.7.9UTP:glucose-1-phosphate uridylyltransferasealpha-d-glucose 1-phosphateUDP-glucose
5.3.1.1triose-phosphate isomerased-glyceraldehyde 3-phosphateglycerone phosphate
Unigene348073.4.21serine-type endopeptidase??
3.4.24metalloendopeptidase??
Unigene389182.7.2.3phosphoglycerate kinase3-phospho-d-glycerate3-phospho-d-glyceroyl phosphate
4.2.1.11phosphopyruvate hydratase2-phospho-d-glyceratephosphoenolpyruvate
Unigene524682.1.1methyltransferase??
2.5.1transferase, transferring alkyl or aryl groups??
Unigene95622.1.2.113-methyl-2-oxobutanoate hydroxymethyltransferase5,10-methylenetetrahydrofolate & 3-methyl-2-oxobutanoatetetrahydrofolate & 2-dehydropantoate
6.3.2.1pantoate-beta-alanine ligase(R)-pantoate (R)-pantothenate
Unigene98042.5.1.21farnesyl-diphosphate farnesyltransferase(2E,6E)-farnesyl diphosphatesqualene
5.3.3.2isopentenyl-diphosphate delta-isomeraseisopentenyl diphosphatedimethylallyl diphosphate
Table
Table 2. List of predicted enzymes with incomplete Enzyme Commission numbers.
Table 2. List of predicted enzymes with incomplete Enzyme Commission numbers.
ClassSub-Subclass EC NumberSub-Subclass Principal Enzyme Type(s)Sub-Subclass Known Enzyme EntriesP. variabilis Number Sequences
Oxidoreductases1.13.11dioxygenase8026
1.14.11dioxygenase, hydroxylase, demethylase569
1.14.12dioxygenase231
1.14.13monooxygenase, hydroxylase2351
1.14.16monooxygenase73
1.14.17monooxygenase313
1.14.19desaturase515
Transferases2.7.10protein-tyrosine kinase272
2.7.11protein-serine/threonine kinase33100
Hydrolases3.1.13exoribonuclease53
3.1.21endodeoxyribonuclease91
3.4.11aminopeptidase2325
3.4.13dipeptidase127
3.4.15peptidyl-dipeptidase41
3.4.16serine-type carboxypeptidase410
3.4.17metallocarboxypeptidase2025
3.4.19omega peptidase129
3.4.21serine endopeptidase100136
3.4.22cysteine endopeptidase5820
3.4.23aspartic endopeptidase4170
3.4.24metalloendopeptidase83120
3.4.25threonine endopeptidase234
Lyases4.1.99carbon-carbon lyases163

Share and Cite

MDPI and ACS Style

R. L. Morlighem, J.-É.; Huang, C.; Liao, Q.; Braga Gomes, P.; Daniel Pérez, C.; De Brandão Prieto-da-Silva, Á.R.; Ming-Yuen Lee, S.; Rádis-Baptista, G. The Holo-Transcriptome of the Zoantharian Protopalythoa variabilis (Cnidaria: Anthozoa): A Plentiful Source of Enzymes for Potential Application in Green Chemistry, Industrial and Pharmaceutical Biotechnology. Mar. Drugs 2018, 16, 207. https://doi.org/10.3390/md16060207

AMA Style

R. L. Morlighem J-É, Huang C, Liao Q, Braga Gomes P, Daniel Pérez C, De Brandão Prieto-da-Silva ÁR, Ming-Yuen Lee S, Rádis-Baptista G. The Holo-Transcriptome of the Zoantharian Protopalythoa variabilis (Cnidaria: Anthozoa): A Plentiful Source of Enzymes for Potential Application in Green Chemistry, Industrial and Pharmaceutical Biotechnology. Marine Drugs. 2018; 16(6):207. https://doi.org/10.3390/md16060207

Chicago/Turabian Style

R. L. Morlighem, Jean-Étienne, Chen Huang, Qiwen Liao, Paula Braga Gomes, Carlos Daniel Pérez, Álvaro Rossan De Brandão Prieto-da-Silva, Simon Ming-Yuen Lee, and Gandhi Rádis-Baptista. 2018. "The Holo-Transcriptome of the Zoantharian Protopalythoa variabilis (Cnidaria: Anthozoa): A Plentiful Source of Enzymes for Potential Application in Green Chemistry, Industrial and Pharmaceutical Biotechnology" Marine Drugs 16, no. 6: 207. https://doi.org/10.3390/md16060207

APA Style

R. L. Morlighem, J. -É., Huang, C., Liao, Q., Braga Gomes, P., Daniel Pérez, C., De Brandão Prieto-da-Silva, Á. R., Ming-Yuen Lee, S., & Rádis-Baptista, G. (2018). The Holo-Transcriptome of the Zoantharian Protopalythoa variabilis (Cnidaria: Anthozoa): A Plentiful Source of Enzymes for Potential Application in Green Chemistry, Industrial and Pharmaceutical Biotechnology. Marine Drugs, 16(6), 207. https://doi.org/10.3390/md16060207

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop