Identification of candidate genes related to calanolide biosynthesis by transcriptome sequencing of Calophyllum brasiliense (Calophyllaceae)
© The Author(s). 2016
Received: 10 February 2016
Accepted: 28 July 2016
Published: 15 August 2016
Calophyllum brasiliense is highlighted as an important resource of calanolides, which are dipyranocoumarins that inhibit the reverse transcriptase of human immunodeficiency virus type 1 (HIV-1 RT). Despite having great medicinal importance, enzymes involved in calanolide, biosynthesis and the pathway itself, are still largely unknown. Additionally, no genomic resources exist for this plant species.
In this work, we first analyzed the transcriptome of C. brasiliense leaves, stem, and roots using a RNA-seq strategy, which provided a dataset for functional gene mining. According to the structures of the calanolides, putative biosynthetic pathways were proposed. Finally, candidate unigenes in the transcriptome dataset, potentially involved in umbelliferone and calanolide (angular pyranocoumarin) biosynthetic pathways, were screened using mainly homology-based BLAST and phylogenetic analyses.
The unigene dataset that was generated in this study provides an important resource for further molecular studies of C. brasiliense, especially for functional analysis of candidate genes involved in the biosynthetic pathways of linear and angular pyranocoumarins.
Calophyllum sp. (Calophyllaceae) is a large group of tropical trees with more than 180–200 species . Currently, some species of this genus have aroused great interest in the scientific community due to their promising phytochemical aspects. In Mexico, the most widely distributed species among the eight found in America is Calophyllum brasiliense Cambes , which grows in tropical rain forests from Brazil to northwest of Mexico . This species contains a large number and variety of secondary metabolites including flavonoids, triterpenes, coumarins, chromones, and xanthones , some of which exhibit interesting anti-leishmanial, anti-bacterial, anti-cancer, anti-parasitic, and anti-viral properties [4, 5]. Two chemotypes have been classified according to their geographical origin. Chemotype 1 (CTP 1), which grows in Sierra de Santa Marta, State of Veracruz, Mexico, produces mammea type coumarins with high in vitro cytotoxic activity against human tumor cells and antibacterial properties against Staphyloccoccus aureus, S. epidermidis and Bacillus subtilis . Meanwhile, chemotype 2 (CTP 2) grows in San Andres Tuxtla, State of Veracruz, Mexico, and produces calanolides, a series of tetracyclic dipyranocoumarins that exhibit an inhibitory effect against the reverse transcriptase of the human immunodeficiency virus type 1 (HIV-1 RT) [2, 7]. There are three different calanolides (A, B and C) that have been found in C. brasiliense and exhibit a significant inhibition on replication of the HIV-1 virus. Interestingly, these bioactive compounds show no toxicity to MT2 human lymphocytes . Additional studies have shown that a high dose of B and C calanolides causes an increased number of spleen megakaryocytes and no alteration of hepatocytes . Calanolide A, which possesses the highest inhibition of viral replication, has been synthesized and has been reported to have similar actions to the natural product [9, 10]. This molecule is in fact in clinical development as a novel therapeutic agent against HIV-1 infection [11, 12]. In plants, calanolides can be detected mainly in leaves (from CTP2), even if they come from seedlings of C. brasiliense that were germinated from seed and grown in a greenhouse . Calanolides can be also detected in plant callus , cell suspension cultures, and leaves from 12-month-old plants that were in vitro regenerated from the young, nodal-stems of C. brasiliense plants .
The metabolic pathways in the biosynthesis of calanolides involve multiple and complex series of enzymatic reactions in which L-phenylalanine and trans-cinnamic acids can be considered as primary precursors. It is important to emphasize that some intermediates such as 7-hydroxycoumarin (umbelliferone) have interesting properties such as anti-fungal [15, 16] and insecticidal  activities, and they have also been shown to have strong inhibitory activity on proliferation of human bladder carcinoma E-J cell lines ( cited in ).
Despite the great pharmacological importance of C. brasiliense, the genomic basis of the synthesis and function of metabolic compounds such as calanolides remains poorly understood. Here, we present the first report of a complete transcriptome analysis of C.brasiliense. The main goal of this study was to characterize the transcriptome of C.brasiliense (CT2) for future gene identification and functional genomics studies of this species. We carried out de novo transcriptome sequencing and assembly of RNA libraries derived from terminal leaves, stems, and roots that come from in vitro regenerated C. brasiliense seedlings. We provide annotation to public databases and categorize the transcripts into biological functions and pathways. In addition, calanolide biosynthetic pathways are suggested, and based on the homologies of some genes, we propose some of them to be promising candidates for future analyses of the calanolide biosynthetic pathway.
Results and discussions
Sequencing and assembly
A total of 16,842,368 paired-end reads (2x150) were generated (5,276,841 for leaves, 5,000,558 for stem and 6,240,602 for roots). Prior to the assembly process, the paired reads were trimmed and/or merged together using the SeqPrep pipeline (see methods for more details). A de novo assembly was generated using Oases , a Bruijn graph-based assembler designed as an extension of Velvet  mainly used to assemble short-read sequences derived from transcriptomics data. Velvet/Oases produced a total of 61,620 contigs ranged from 0.1 to 10 kb, with an average length of 547.28 bp (Additional file 1). The GC contents of the contig set was approximately 44.7 %, which is similar to the GC content of the coding regions from other species within the Malpighiales order (reviewed in ). The N50 of these contigs was also estimated and resulted in a moderately high value of 867 bp. A fairly large number (40,727) of assembled contigs (40,727, which represents a 66.01 % of the total), were between 200 bp and 500 bp in length, indicating the presence of assembled fragments. For practical purposes, in the presented work, all contigs from the dataset are referred to as unigenes. The BLASTx algorithm  was used to annotate the unigenes based on the traditional top-BLAST-hit annotation method. As references, a collection of protein databases including the Arabidopsis thaliana (Arabidopsis) and plant RefSeq databases were used for this purpose. A significant value (e-value) of 10−5 was applied as threshold in the BLASTx similarity searches (Additional file 2: Table S1).
Based on the Arabidopsis top hits, Gene Ontology (GO) annotations for the C. brasiliense unigenes were obtained. WEGO software  was used to perform GO functional classification into the three major categories (biological process, molecular function, and cellular components). Among the unigenes with Arabidopsis hits, 42,090 (68.30 %) were assigned to different gene ontology categories with a total of 367,994 functional terms, of which 103,865 are unique. Biological processes comprised the majority of the functional terms (178,629; 48.54 %), followed by cellular components (95,428; 25.93 %) and molecular functions (93,937; 25.52 %) (Additional file 3: Figure S1; see also Additional file 4: Table S2). Top-ranked categories of GO biological processes were the sub-categories corresponding to cellular (27,090 unigenes) and metabolic (24,653 unigenes) processes. Interestingly, response to stimulus (14,101 unigenes) and biological regulation (12,646 unigenes) were also prominently represented among GO biological processes categories. In addition to functional annotation based on GOs, C. brasiliense unigenes were classified based on metabolic pathways available and described in Kyoto Encyclopedia of Genes and Genomes (KEGG). KEGG Automatic Annotation Server (KAAS; ) was used to assign to C.brasiliense unigenes the KEGG Orthology (KO) codes and enzyme commission (EC) numbers. KO codes were assigned to 4,881 unigenes, of which 1,733 could be associated to specific EC numbers related to 226 different metabolic pathways (Additional file 2: Table S1).
Gene expression profiles of C. brasiliense organs
In order to gain insight into the organ-function connection, the top 10 organ-specific unigenes were surveyed. In leaves, three C. brasiliense unigenes (UN36044, UN28345 and UN34582), which are homologous to Arabidopsis members of the glycine-aspartic acid-serine-leucine motif lipase/ hydrolase (GDSL lipase) family, were highly expressed as well as some homologs (UN09544 and UN13106) of 3-ketoacyl-CoA synthase proteins (KSC). Consistently, members of the GDSL lipase gene family, such as AT5G33370 and AT3G04290, are co-regulated with genes involved in cutin biosynthesis as CER6/KCS6/CUT1 (AT1G68530), which has a dominant role in the elongation of very-long-chain fatty acids for cuticular wax synthesis [27–29]. The lipoxygenase/peroxygenase pathway is also involved in biosynthesis of cutin monomers ; this could explain the expression profile of the UN00603, which was identified as a leaf-specific unigene and annotated as homologous to chloroplast lipoxygenase LOX2 (AT3G45140), an enzyme required for wound-induced jasmonic acid accumulation in Arabidopsis . The presence of all of these genes in the leaf transcriptome was expected considering that epicuticular waxes are produced either exclusively during leaf development and expansion, or during the entire lifetime of the leaf .
Regarding stems, a homolog of fasciclin-like arabinogalactan protein FLA12 (AT5G60490; unigene UN35075) exhibited one of the highest expression levels as a stem-specific gene, which is consistent with previous reports showing that the expression of some members of the FLA gene family are correlated with the onset of secondary-wall cellulose synthesis in Arabidopsis stems, and with wood formation in the stems and branches of trees. This data suggests that unigene UN35075 may play a biological role in C. brasiliense stem development . Additionally, genes encoding enzymes related to monolignol biosynthesis, such as phenylalanine ammonia-lyase (PAL1; UN01310), caffeoyl-CoA 3-O-methyltransferase (CCoAOMT; UN21637 and UN12250), and 4-coumarate: CoA ligase (4CL; UN01988), were identified as preferentially expressed in stems, although transcripts could be detected in all three organs sampled. Lignin, which plays a crucial role in conducting tissue in plant stems, is synthesized from the oxidative coupling of monolignols . In addition, 16 unigenes homologous to Arabidopsis IRX proteins were also classified as specifically or preferentially expressed in stems; this was expected considering that the irregular xylem (irx) mutant is characterized by a reduction in cellulose in stem tissue . Finally, two C. barsiliense unigenes (UN02173 and UN03226), which were homologous to members of the family of high affinity phosphorous transporters (PHT1), were identified as highly expressed only in roots, as well as the UN03833 unigene, a homolog of ARSK (AT2G26290), a poorly characterized gene encoding a root-specific kinase .
In order to validate the expression profiles obtained by normalized read counts, RT-qPCR was performed using nine chosen genes. All genes evaluated showed RT-qPCR expression profiles in complete agreement with the profiles derived from read counts analyses (Fig. 1).
Functional annotation of preferentially expressed genes
The biosynthesis of umbelliferone, a key precursor in the calanolides formation
Coumarins are synthesized in plants via the shikimate pathway, in which phenylalanine is an end product that also gives rise to the aromatic amino acids tyrosine and tryptophan and other small molecules such as flavonoids and hydroxycinnamic acid conjugates . Successive para- and ortho- hydroxylation of trans-cinnamate (conjugate base of trans-cinnamic acid) leads to the formation of coumarin via 2-coumarate, or via 4-coumarate, to the formation of hydroxycoumarins such as umbelliferone (7-hydroxycoumarin). Other hydroxycoumarins lacking oxygenation at C-7 also share the trans-cinnamic acid as its precursor. According to EC numbers assigned to C. brasiliense unigenes, with only one exception (glutamate-prephenate aminotransferase; EC: 184.108.40.206), homologs from all enzymes required for the formation of L-phenylalanine via the shikimate pathway were identified (Additional file 3: Figure S3).
Similar approaches to those described above were used to identify the remaining C. brasiliense genes involved in umbelliferone biosynthesis. A total of three unigenes were identified as homologous to the only trans-cinnamate 4-hydroxylase (AT2G30490; C4H) in the Arabidopsis genome (Additional file 3: Figure S4 and Additional file 7). C4H is a plant-specific cytochrome P450 (PF00067) that catalyzes the second step of the multibranched phenylpropanoid pathway . Regarding 4-coumarate CoA ligase [4CL; EC:220.127.116.11], a total of six unigenes were detected as homologous to these proteins; however, only in three of them (UN01603, UN01988 and UN01725), were complete open reading frames identified (Additional file 8). The motif/domain searches revealed that both the AMP-binding C-terminal domain (PF13193) and common AMP-binding central domain (PF00501) are present in the translated sequences corresponding to these C. brasiliense unigenes (Additional file 3: Figure S5). Commonly, most angiosperms encode to a small family of 4CL (e.g., seven members in case of Arabidopsis). These enzymes are involved in the last step of the general phenylpropanoid pathway, and in addition to using 4-coumarate as substrate, they also convert p-coumaric acid, ferulic acid, caffeic acid and 5-OH-ferulic acid with different catalytic efficiency [43, 44].
Finally, the bi-directional best hit (BBH) method was used to identify a homolog of 4-coumaroyl 2′-hydroxylase [EC:1.14.11.-] from Ruta graveolens. The 4-coumaroyl 2′-hydroxylase isolated/characterized from R. graveolens (Accession JF799117.1) is the only enzyme that has been specifically assigned to coumarin synthesis, and to a lesser extent this enzyme also accepts 4-coumaroyl-CoA to produce umbelliferone . The C. brasiliense unigene UN02124, in which the cytochrome P450 conserved domain PF00067 is present, showed 87.2 % identity with the published protein from R. graveolens (Additional file 3: Figure S6 and Additional file 9). The final step of the in vivo pathway for the synthesis of umbelliferone involves a trans-cis isomerization followed by a subsequent lactonization of the 2, 4-dihydroxy-cinnamoyl-CoA that closes the side chain, and this reaction occurs spontaneously (Fig. 3).
Considering that we were able to identify homologs to all genes involved in the umbelliferone biosynthetic pathway via 4-coummarate, and due the absence of several transcripts that potentially encode for enzymes such as cinnamate 2-hydroxylase (EC: 18.104.22.168), 2-coumarate O-β-glucosyltransferase (EC: 22.214.171.124), 2-coumarate β-D-glucoside isomerase (EC: 5.2.1.-) and coumarinic acid glucoside β-glucosidase (EC:126.96.36.199), which were first characterized in Melilotus alba [46, 47] were they are all involved in coumarin biosynthetic pathway via 2-coumarate, we suggest that in C. brasiliense, umbelliferone is synthesized via the 4-coumarate pathway. This was expected, because in contrast with mammals, only in a few plant species (e.g. Catharanthus roseus and Conium maculatum) has it been suggested that enzymes capable of carrying out the hydroxylation of coumarin in C-7 to produce umbelliferone exist (reviewed in ).
Analysis of putative candidate genes involved in the calanolide (angular pyranocoumarins) biosynthetic pathway
Previous studies conducted in Pimpinella magna and Pastinaca sativa plants could be considered the first experimental evidence to prove that linear and angular furanocuoumarins are derived from umbelliferone, prenylated at either the C-6 (leading to the formation of demethylsuberonosin; DMS) or C-8 position (osthenol), respectively . In later years, additional investigations revealed that the cyclization of demethylsuberosin leads to (+)-marmesin formation, taking place through an enzymatic reaction that occurs in the presence of NADPH and molecular oxygen . These ‘mamersin synthases’ have been identified as cytochrome P450 monooxygenases in Petroselinum crispum and Ammi maju plant species (reviewed in ). The range of reactions catalyzed by P450s include the epoxidation of olefins by insertion of an ‘oxen’ , and the reactive product of this reaction often inactivates the enzyme by alkylation of the prosthetic heme group . However, no intermediate was released from mamersin synthase reaction, and it is likely that the 7-hydroxyl group of demethylsuberosin delocalizes the double bond electrons and favors the instantaneous cyclization to the dihydrofuranocoumarin (Fig. 5a). Model mechanisms have been proposed for the reactions mediated by catalytic action of the P450 enzymes, and one of these mechanisms consists of primary interaction of the catalytic P450 oxo-derivative, formed by heterolytic cleavage of the oxygen-oxygen bond in the ferric-hydroperoxy species, with aliphatic double bonds [52, 54]. This mechanism of reaction is compatible with such a cyclization of demethylsuberosin to (+)-marmesin, avoiding the formation of an intermediate epoxide (Fig. 5a). The formation of (−)-columbianetin from osthenol is catalyzed in an analogous fashion and the subsequent activities of both enzymes, angelicin and psoralen synthases, have been supported experimentally [55, 56]. It appears feasible that the synthesis of linear and angular pyranocoumarins such xanthyletin , graveolone , or seselin  may be produced concomitantly with furanocoumarins in a very similar way as from demethylsuberosin or osthenol, respectively (Fig. 5b).
Considering the remarkable structural similarities between furanocoumarins and pyranocoumarins (linear and angular), and the fact that both classes of compounds are derived from the same precursors (DMS and osthenol), it is possible to hypothesize that the enzymes involved in these biosynthetic pathways might share common ancestry. This hypothesis is consistent with the observation that furanocoumarins and pyranocuoumarins do not usually coexist in the same species. Even when many unknowns about biosynthesis and the function of furanocoumarins remain unresolved, it is clear that our knowledge about closely related compounds, such as pyranocoumarins, appears to be even less.
Concerning furanocoumarins, it has been reported that they are produced by a wide variety of plants in response to pathogen or herbivore attack. They are activated by ultraviolet light and can be highly toxic to certain vertebrate and invertebrate herbivores due to their integration into DNA, which contributes to rapid cell death [49, 62]. It has been suggested that linear and angular furanocoumarins are the results of a co-evolutionary process between plant and insects. Plant-insect interaction studies reveal that linear furanocoumarins are more toxic than angular ones; however, angular structures apparently produce a synergistic effect when they are combined with linear ones . When mixed, linear and angular furanocoumarins result in a combination that is more difficult for insects to detoxify . Apparently, during evolution, angular furanocoumarins appeared later than linear ones, a hypothesis that finds support based on the observation that angular furanocoumarins are always found concomitantly with linear structures, while linear types can be found alone (reviewed in [49, 63]).
During recent years, many enzymes involved in furanocoumarin biosynthesis have been described at the molecular level, including three P450 monooxygenases (psoralen-, angelicin-, and (+)-marmesin- synthases [55, 56]), bergaptol O-methyltransferase , as well as the key enzyme, a prenyltransferase ) involved in the critical step leading to precursors synthesis of linear and angular furanocoumarins (DMS and osthenol respectively; see Figs. 5 and 6). Additional P450-dependent enzymatic steps have remained unresolved, but the participation of this class of enzymes has been suggested in many steps of furanocoumarin biosynthesis. (Figs. 5 and 6, see  for an extended revision). Furanocoumarins are produced by a wide variety of plants in response to pathogen or herbivore attack.
The biosynthesis of pyranocoumarins has not yet been investigated. However, considering that psoralen synthase shows 70 % identity with angelisin synthase, and that its participation in linear and angular furanocoumarins biosynthesis has been previously demonstrated [55, 56], and based on the assumption that enzymes involved in pyranocoumarin biosynthesis might share a common ancestor of unknown functionality (perhaps as a result of gene duplications and subsequent molecular evolution), we used these and other enzymes involved in furanocoumarin biosynthesis as sequence references to identify homologs in the C. brasiliense unigenes set.
First, the prenyltransferase (PT) identified in parsley (Petroselinum crispum) , the angelicin synthase (CYP71AJ4) from Pastinaca sativa, the psoralen synthase (CYP71AJ1) identified from Ammi majus and their orthologs (CYP71AJ2 and CYP71AJ3 from Apium graveolens and Pastinaca sativa, respectively) [55, 56], were used as references in tBLASTn similarity searches (e-value 10−6) against the C. brasiliense unigene dataset. Additionally, the poorly characterized CYP82H1 isolated from Ammi majus, which increases its expression levels in plant-fungi interactions and accompanying furanocoumarin biosynthesis, was also included . Sequences of unigenes that upon translation show at least 20 % identity over a resides window that represents a complete coding sequence (CDS), or at least 70 % of the homologous protein, were retained for future analysis.
A total of six PT-like sequences were identified. The subsequent comparison against PFAM databases confirmed the presence of the UbiA prenyltransferase domain (PF01040). The unigene UN01964 (with 35 % shared identity) was resolved in the same clade as the parsley prenyltransferase (Additional file 3: Figure S8 and Additional file 10). Considering these results, we propose that UN01964 should be considered a leading candidate to encode a putative prenyl transferase involved in biosynthesis of the precursors that lead to synthesis of linear and angular pyranocoumarins.
Regarding the subsequent steps in pyranocoumarin biosynthesis, a total of 34 unigenes were identified as homologous to the previously characterized P450 monooxygeneases involved in the furanocoumarin biosynthetic pathway (stringency levels: coverage ≥ 70 %, identity ≥ 20 %). Three major clades were recognized on the phylogenetic tree, and one of these included only three C. brasiliense unigenes. The second clade brought together all CYP71AJ proteins that were included in phylogenetic analyses and a total of fourteen C. brasiliense unigenes. Finally, the third clade included the CYP82H1 P450 monooxygenase and the remaining C. brasiliense CYP-like sequences that were identified (Additional file 3: Figure S9, Additional file 11). The CYP71-related clade comprises two sister clades named as classes I and II, respectively. The translated Class I unigenes showed a percent identity that ranged from 31 to 46 % with respect to CYP71AJ proteins, while class II members were on average ~10 % less similar (Additional file 4: Table S8). Low percent similarities at the protein level were expected considering that the presence of furanocoumarins has not been reported in C. brasiliense, which is instead capable of producing pyranocoumarins (linear and angular, including the calonolide compounds) in young seedlings (mainly at leaves) and in callus cultures . According to the expression profiles of these genes, 44 % were identified as preferentially expressed in some of the organs sampled. UN02363, UN03063, UN04124, and UN02841 were selected as preferentially expressed in leaves. With only one exception (UN02841), these genes were classified as CYP71-related; two of them (UN02363 and UN03063) grouped in class II and the other one into class I (UN04124). In addition, UN04124 possesses 35.6 and 35.8 % identity with the psoralen and angelicin synthases of Pastinaca sativa. Altogether, this data suggests that UN04124 can be considered a prime candidate for involvement in angular and/or linear pyranocoumarin biosynthesis.
The unigene dataset generated in this study provides an important resource for further molecular studies of C. brasiliense, especially for characterizing candidate genes in the biosynthetic pathways of linear and angular pyranocoumarins. Using appropriate approaches, a series of candidate genes were identified. Consecutive analyses were conducted to determine their corresponding expression patterns and phylogenetic relationships. The candidate genes identified in C. brasiliense transcriptome that were suggested to possibly be involved in the biosynthesis of calanolides (angular coumarins), could be cloned and characterized in further studies. Additional bioinformatic analyses could be conducted in order to reduce the number of candidate proteins that may catalyze specific reactions at particular steps. We suggest that at least for P450 monooxygenase enzymes, docking and molecular dynamics analyses could be performed in order to reduce the number of candidate genes involved in the calanolide biosynthetic pathway.
Previously, Bernabé et al., , have shown that callus cultures or young plants of C. brasiliense grown in a greenhouse are capable of producing calanolides. In order to guarantee the production of calanolides, in vitro seedlings were regenerated from the young nodal-stem of the same plants (chemotype 2) used previously by Bernabé et al. These plants were germinated from seeds collected in San Andrés Tuxtla, State of Veracruz. The plants were grown in a green house during 5–6 months approximately. Young nodal-stem were collected from the plants and used as the source of explants. Leaves, stems, and roots were collected from in vitro regenerated 12-month-old plants. According to the manufacturer’s instructions, total RNAs were isolated with TRIzol® Reagent (Life technologies). RNAs isolated from the three organs sampled, were re-purified with the RNeasy kit (Qiagen) and treated with RNase-free DNase I (Invitrogen) in order to remove the DNA residues. The quality and purity of RNAs were assessed with OD260/230 ratio by using the NanoDrop 2000 (Thermo Fisher). RNA integrity was evaluated by RNA integrity number (RIN) using an Agilent 2100 Bioanalyzer (Agilent Technologies). Only RNAs with RIN values ≥ 8.5 were used from library generation.
cDNA library preparation and sequencing (RNA-seq)
cDNA preparation, library construction, and sequencing were performed at the Genomic Services laboratory, LANGEBIO-CINVESTAV, Mexico by using the MiSeq™ platform according to the manufacturer’s instructions (Illumina, San Diego, CA). Briefly, poly (A) RNA was isolated from 20 μg of total RNA using Sera-mag Magnetic Oligo (dT) Beads (Illumina). To avoid priming bias when synthesizing cDNA, the purified mRNA was first fragmented into small pieces. Then the double-stranded cDNA was synthesized using the SuperScript Double-Stranded cDNA Synthesis kit (Invitrogen, Camarillo, CA) with random hexamer (N6) primers (Illumina). The synthesized cDNA was subjected to end-repair and phosphorylation using T4 DNA polymerase, Klenow DNA polymerase and T4 PNK. These repaired cDNA fragments were 3′ adenylated using Klenow Exo- (3′ to 5′ exo minus, Illumina). Illumina paired-end adapters were ligated to the ends of these 3′-adenylated cDNA fragments, using specific barcodes for each sample (roots, stem, and leaves). Fifteen rounds of PCR amplification were performed to enrich the purified cDNA template using PCR Primer PE 1.0 and PE 2.0 (Illumina) with Phusion DNA Polymerase. The cDNA library was constructed with a fragment length range of 278 bp (±0.5 SD). Finally, after validating on an Agilent Technologies 2100 Bioanalyzer using the Agilent DNA 1000 chip kit, cDNA libraries were sequenced on a paired-end (2x150) flow cell using Illumina MiSeq sequencer. Files containing sequence reads and quality scores were deposited in the Short Read Archive of the National Center for Biotechnology Information (NCBI) [Accession number SRP079249].
Data filtering and de novo assembly
Forward and reverse read pairs (generated by Illumina-MiSeq) were merged to form a single “longer-reads” using the SeqPrep pipeline (https://github.com/jstjohn/SeqPrep), with default parameters (a quality score cutoff of phred 33, a minimum merged read length of 15 bp and no mismatches in the overlapping region). Paired-end reads that did not overlap were trimmed with a sliding window approach (window size 10 bases, shift 1 base). The FASTX-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/index.html) was used to this purpose. Reads were discarded if they were smaller than 30 bases after trimming, orphan reads were also removed in order to keep pairs only. Velvet assembler using the Oases module  was used for sequence assembly. A unigenes set from C. brasiliense was generated considering only resulting contigs with a minimum size of 100 bp.
Annotation of C. brasiliense unigenes
To annotate unigenes obtained by de novo assembly, we performed sequence similarity searches using the BLASTx algorithm (e-value 10−5) on Arabidopsis thaliana (TAIR v11), and other plant proteins (NCBI; ftp://ftp.ncbi.nlm.nih.gov/refseq/release/plant/) databases. Top protein matches from Arabidopsis or additional plant proteins were assigned to each of the C. brasiliense unigenes. The gene ontology (GO) functional classes and pathways for each C. brasiliense unigene were assigned based on Arabidopsis GO SLIM and pathway annotation (ftp://ftp.arabidopsis.org/home/tair/Ontologies/). The data were statistically analyzed using the WEGO software  which is a useful tool for plotting GO annotation results. Additionally, the unigenes were also analyzed using the KEGG Automatic Annotation Server (KAAS; http://www.genome.jp/tools/kaas/) to provide annotations of KEGG Orthology (KO) codes. The bi-directional best hit (BBH) method was used. Enzyme Commission (EC) numbers were also assigned based on the annotations extracted from Kyoto Encyclopedia of Genes and Genomes (KEGG) and were cross-checked with orthologous gene annotation projections from Plant Metabolic Network (PMN; http://www.plantcyc.org/), and if available, corresponding EC-formatted MetaCyc  cross-references were added.
Expression profile analysis of C. brasiliense transcriptome
After assembling of the C. brasiliense transcriptome, every RNA-seq library was separately aligned to the generated transcriptome assembly using Bowtie . Counting of alignments was done using RSEM . Differential expression statistical analysis was done using the FPKM (fragments per Kilo bases of contigs for per million mapped reads) values and statistical method described by Sketel . Briefly, all clusters were submitted to a log likelihood ratio statistics which tends asymptotically to a χ2 distribution as described by Stekel. It is based on a single statistical test to describe the extent to which a gene is differentially expressed between libraries. This method permits in any number of libraries to identify differential expressed genes.
In order to identify the coding sequence in their correct open reading frame, eight of the unigenes of C. brasiliense identified as differentially expressed genes were aligned versus their Arabidopsis homologues. Using the SeaView program, protein-coding nucleotide sequences were aligned based on their corresponding amino acid translations (Additional file 12). Gene-specific primer pairs (Additional file 4: Table S9), were designed using the Primer3 v.0.4.0 web tool (http://bioinfo.ut.ee/primer3-0.4.0/primer3/) and then used in RT-qPCR assays.
A total of 10 ug of total RNA was reverse transcribed using SuperScript® III Reverse Transcriptase (Life Technologies) according to the manufacturer’s instructions. qPCR of selected genes was carried out through SYBR green chemistry (Applied Biosystem) on a real time thermal cycler (AB7500, Applied Biosystem). UBQ11 (UN18770) was used as an internal control. The thermal cycling program was set to 95 °C for 5 min, 40 cycles of 95 °C for 30 s, 60 °C for 30 s, and 72 °C for 1 min. All reactions were run in duplicates of three biological replicates.
A maximum likelihood framework was used in phylogenetic analyses, which were performed with SeaView v2.4 software . The alignments and phylogenetic analysis were drived by SeaView using Muscle  and PhyML  programs, respectively. Topology, branch lengths and, equilibrium frequencies, were optimized and the PhyML option was used under LG (Le and Gascuel) model . The starting tree was determined using BioNJ, and both nearest-neighbor interchange (NNI) and subtree pruning and regrafting (SPR) algorithms for tree searching were used. Branch robustness was analyzed by approximate likelihood-ratio test (aLRT) .
The C. brasiliense’s 3D protein structures were modeled by the rigid body grouping method, using the Swiss-Model workspace (http://swissmodel.expasy.org/) [73, 74]. This server is used to align the target sequences and template structure available in the Protein Data Bank (PDB). Once the template has been selected the 3D structure of the target sequences can be modeled. Templates used for modeling were 1w27.1 for phenylalanine ammonia lyases (PAL) and 3a9u.1 for 4-coumarate: CoA ligases (4CL). Each model generated was checked for various parameters that include Z, GMQE (Global Model Quality Estimation) and QMEAN (Qualitative Model Energy ANalysis) scores to assess the accuracy of the model. The modeled C. brasiliense proteins were superimposed onto their corresponding homologues structures using the SWISS-PDB viewer v4.1.0 program .
We would like to thank Victor A. Albert for his positive and relevant comments which improved the quality of this manuscript. We also thank Victor’s support for language editing. Finally, we wish to thank Aldo Segura-Cabrera, Ana-Luisa Kiel-Martínez and Ofelia Ferrera-Rodríguez for fruitful discussions and valuable suggestions.
This work was supported by Consejo Nacional de Ciencia y Tecnología (CONACyT), grant 223323 (EIL).
Availability of data and materials
All supporting data are included within the article or in the additional files.
HBG-R wrote the paper with significant contributions by EI-L, HBG-R, FC-S, AB-A, AG-A, JLO-R, AA-S, EV and EI-L collected data. HBG-R, AG-A, JLO-R, AA-S, EV and EI-L analyzed data. EI-L conceived of and led the study. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Stevens PF. A revision of the Old World species of Calophyllum (Guttiferae). Journal of the Arnold Arboretum. 1980;61:117–699.View ArticleGoogle Scholar
- Huerta-Reyes M, Basualdo Mdel C, Lozada L, Jimenez-Estrada M, Soler C, Reyes-Chilpa R. HIV-1 inhibition by extracts of Clusiaceae species from Mexico. Biol Pharm Bull. 2004;27(6):916–20.View ArticlePubMedGoogle Scholar
- Pennington TD, Sarukhán J. Arboles tropicales de México: manual para la identificación de capo de los principales. Mexico: Instituto Nacional de Investigaciones Forestales; 1968.Google Scholar
- Bernabé-Antonio A, Álvarez-Berber LP, Cruz-Sosa F. Biological Importance of Phytochemicals from Calophyllum brasiliense Cambess. Annual Research & Review in Biology. 2014;4(10):1502–17.View ArticleGoogle Scholar
- Brenzan MA, Santos AO, Nakamura CV, Filho BP, Ueda-Nakamura T, Young MC, Correa AG, Junior JA, Morgado-Diaz JA, Cortez DA. Effects of (−) mammea A/BB isolated from Calophyllum brasiliense leaves and derivatives on mitochondrial membrane of Leishmania amazonensis. Phytomedicine : international journal of phytotherapy and phytopharmacology. 2012;19(3–4):223–30.View ArticleGoogle Scholar
- Reyes-Chilpa R, Estrada-Muniz E, Apan TR, Amekraz B, Aumelas A, Jankowski CK, Vazquez-Torres M. Cytotoxic effects of mammea type coumarins from Calophyllum brasiliense. Life Sci. 2004;75(13):1635–47.View ArticlePubMedGoogle Scholar
- Huerta-Reyes M, Basualdo Mdel C, Abe F, Jimenez-Estrada M, Soler C, Reyes-Chilpa R. HIV-1 inhibitory compounds from Calophyllum brasiliense leaves. Biol Pharm Bull. 2004;27(9):1471–5.View ArticlePubMedGoogle Scholar
- García-Zebadúa JC, Magos-Guerrero GA, Mumbrú-Massip M, Estrada-Muñoz E, Contreras-Barrios MA, Huerta-Reyes M, Campos-Lara MG, Jiménez-Estrada M, Reyes-Chilpa R. Inhibition of HIV-1 reverse transcriptase, toxicological and chemical profile of Calophyllum brasiliense extracts from Chiapas, Mexico. Fitoterapia. 2011;82(7):1027–34.View ArticleGoogle Scholar
- Flavin MT, Rizzo JD, Khilevich A, Kucherenko A, Sheinkman AK, Vilaychack V, Lin L, Chen W, Greenwood EM, Pengsuparp T, et al. Synthesis, chromatographic resolution, and anti-human immunodeficiency virus activity of (+/−)-calanolide A and its enantiomers. J Med Chem. 1996;39(6):1303–13.View ArticlePubMedGoogle Scholar
- Flavin MT, Xu ZQ, Rizzo JD, Kucherenko A, Khilevich A, Sheinkman AK, Wilaychack V, Lin L, Chen W, Boulanger WA. Method for the preparation of (±)-calanolide a and intermediates thereof. In. Google Patents; 1996.Google Scholar
- Cragg GM, Newman DJ. Plants as a source of anti-cancer and anti-HIV agents. Ann Appl Biol. 2003;143(2):127–33.View ArticleGoogle Scholar
- Kashman Y, Gustafson KR, Fuller RW, Cardellina 2nd JH, McMahon JB, Currens MJ, Buckheit Jr RW, Hughes SH, Cragg GM, Boyd MR. The calanolides, a novel HIV-inhibitory class of coumarin derivatives from the tropical rainforest tree, Calophyllum lanigerum. J Med Chem. 1992;35(15):2735–43.View ArticlePubMedGoogle Scholar
- Bernabé-Antonio A, Estrada-Zúñiga ME, Buendía-González L, Reyes-Chilpa R, Chávez-Ávila VM, Cruz-Sosa F. Production of anti-HIV-1 calanolides in a callus culture of Calophyllum brasiliense (Cambes). Plant Cell Tiss Organ Cult. 2010;103(1):33–40.View ArticleGoogle Scholar
- Cisneros-Torres D. Análisis de compuestos activos en extractos de hoja y de cultivo de células en suspensión de Calophyllum brasiliense Cambess y evaluación de la actividad anti-inflamatoria. México: Universidad Autónoma Metropolitana - Iztapalapa; 2015.Google Scholar
- Bai X-n, Liang W, Cheng J, Ma L-Q, Liu Y-B, Shi G-l, Wang Y-N, Gu J-C. Inhibitory effect and antifunal mechanism of umbelliferone on plant pathogenic fungi. In: Zhu E, Sambath S, editors. Information Technology and Agricultural Engineering. Berlin Heidelberg: Springer; 2012. p. 693–702.View ArticleGoogle Scholar
- Montagner C, De Souza SM, Groposoa C, Delle Monache F, Smania EF, Smania Jr A. Antifungal activity of coumarins. Zeitschrift fur Naturforschung C, Journal of biosciences. 2008;63(1–2):21–8.PubMedGoogle Scholar
- Zhang G, Wang Y, Xu H, Wu G, Zhao S. Isolation and identification of extraction of Stellera chamaejasme. Anhui Nongxueyuan Xuebao. 2000;27(4):345–7.Google Scholar
- Yang XW, Xu B, Ran FX, Wang RQ, Wu J, Cui JR. Inhibitory effects of 11 coumarin compounds against growth of human bladder carcinoma cell line E-J in vitro. Zhong xi yi jie he xue bao = Journal of Chinese integrative medicine. 2007;5(1):56–60.View ArticlePubMedGoogle Scholar
- Schulz MH, Zerbino DR, Vingron M, Birney E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012;28(8):1086–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18(5):821–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Ibarra-Laclette E, Albert VA, Herrera-Estrella A, Herrera-Estrella L. Is GC bias in the nuclear genome of the carnivorous plant Utricularia driven by ROS-based mutation and biased gene conversion? Plant Signal Behav. 2011;6(11):1631–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.View ArticlePubMedPubMed CentralGoogle Scholar
- 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–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.View ArticlePubMedGoogle Scholar
- Stekel DJ, Git Y, Falciani F. The comparison of gene expression from multiple cDNA libraries. Genome Res. 2000;10(12):2055–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Fiebig A, Mayfield JA, Miley NL, Chau S, Fischer RL, Preuss D. Alterations in CER6, a gene identical to CUT1, differentially affect long-chain lipid content on the surface of pollen and stems. Plant Cell. 2000;12(10):2001–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Joubes J, Raffaele S, Bourdenx B, Garcia C, Laroche-Traineau J, Moreau P, Domergue F, Lessire R. The VLCFA elongase gene family in Arabidopsis thaliana: phylogenetic analysis, 3D modelling and expression profiling. Plant Mol Biol. 2008;67(5):547–66.View ArticlePubMedGoogle Scholar
- Kannangara R, Branigan C, Liu Y, Penfield T, Rao V, Mouille G, Hofte H, Pauly M, Riechmann JL, Broun P. The transcription factor WIN1/SHN1 regulates Cutin biosynthesis in Arabidopsis thaliana. Plant Cell. 2007;19(4):1278–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Blée E, Schuber F. Biosynthesis of cutin monomers: involvement of a lipoxygenase/peroxygenase pathway. Plant J. 1993;4(1):113–23.View ArticleGoogle Scholar
- Bell E, Creelman RA, Mullet JE. A chloroplast lipoxygenase is required for wound-induced jasmonic acid accumulation in Arabidopsis. Proc Natl Acad Sci U S A. 1995;92(19):8675–9.View ArticlePubMedPubMed CentralGoogle Scholar
- MacMillan CP, Mansfield SD, Stachurski ZH, Evans R, Southerton SG. Fasciclin-like arabinogalactan proteins: specialization for stem biomechanics and cell wall architecture in Arabidopsis and Eucalyptus. Plant J. 2010;62(4):689–703.View ArticlePubMedGoogle Scholar
- Vanholme R, Demedts B, Morreel K, Ralph J, Boerjan W. Lignin biosynthesis and structure. Plant Physiol. 2010;153(3):895–905.View ArticlePubMedPubMed CentralGoogle Scholar
- Turner SR, Somerville CR. Collapsed xylem phenotype of Arabidopsis identifies mutants deficient in cellulose deposition in the secondary cell wall. Plant Cell. 1997;9(5):689–701.View ArticlePubMed CentralGoogle Scholar
- Hwang I, Goodman HM. An Arabidopsis thaliana root-specific kinase homolog is induced by dehydration, ABA, and NaCl. Plant J. 1995;8(1):37–43.View ArticlePubMedGoogle Scholar
- Dewick PM. The shikimate pathway: aromatic amino acids and phenylpropanoids. In: Medicinal Natural Products. Chichester: John Wiley & Sons, Ltd; 2009. p. 137–86.View ArticleGoogle Scholar
- Hyun MW, Yun YH, Kim JY, Kim SH. Fungal and plant phenylalanine ammonia-lyase. Mycobiology. 2011;39(4):257–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Birney E, Clamp M, Durbin R. GeneWise and Genomewise. Genome Res. 2004;14(5):988–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Gouy M, Guindon S, Gascuel O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27(2):221–4.View ArticlePubMedGoogle Scholar
- Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, et al. The Pfam protein families database. Nucleic Acids Res. 2004;32(Database issue):D138–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Wanner LA, Li G, Ware D, Somssich IE, Davis KR. The phenylalanine ammonia-lyase gene family in Arabidopsis thaliana. Plant Mol Biol. 1995;27(2):327–38.View ArticlePubMedGoogle Scholar
- Mizutani M, Ohta D, Sato R. Isolation of a cDNA and a genomic clone encoding cinnamate 4-hydroxylase from Arabidopsis and its expression manner in planta. Plant Physiol. 1997;113(3):755–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Ehlting J, Buttner D, Wang Q, Douglas CJ, Somssich IE, Kombrink E. Three 4-coumarate:coenzyme A ligases in Arabidopsis thaliana represent two evolutionarily divergent classes in angiosperms. Plant J. 1999;19(1):9–20.View ArticlePubMedGoogle Scholar
- Hamberger B, Hahlbrock K. The 4-coumarate: CoA ligase gene family in Arabidopsis thaliana comprises one rare, sinapate-activating and three commonly occurring isoenzymes. Proc Natl Acad Sci U S A. 2004;101(7):2209–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Vialart G, Hehn A, Olry A, Ito K, Krieger C, Larbat R, Paris C, Shimizu B, Sugimoto Y, Mizutani M, et al. A 2-oxoglutarate-dependent dioxygenase from Ruta graveolens L. exhibits p-coumaroyl CoA 2′-hydroxylase activity (C2′H): a missing step in the synthesis of umbelliferone in plants. Plant J. 2012;70(3):460–70.View ArticlePubMedGoogle Scholar
- Gestetner B, Conn EE. The 2-hydroxylation of trans-cinnamic acid by chloroplasts from Melilotus alba Desr. Arch Biochem Biophys. 1974;163(2):617–24.View ArticlePubMedGoogle Scholar
- Poulton JE, McRee DE, Conn EE. Intracellular localization of two enzymes involved in coumarin biosynthesis in Melilotus alba. Plant Physiol. 1980;65(2):171–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Robbins MP. Biochemistry of plant secondary metabolism. Annual Plant Reviews, Volume 2. Edited by Michael Wink. Eur J Plant Pathol. 2000;106(5):487.View ArticleGoogle Scholar
- Bourgaud F, Hehn A, Larbat R, Doerper S, Gontier E, Kellner S, Matern U. Biosynthesis of coumarins in plants: a major pathway still to be unravelled for cytochrome P450 enzymes. Phytochem Rev. 2006;5(2):293–308.View ArticleGoogle Scholar
- Harborne JB. The natural coumarins: occurrence, chemistry and biochemistry (Book). Plant Cell Environ. 1982;5(6):435–6.View ArticleGoogle Scholar
- Matern U, Strasser H, Wendorff H, Hamerski D. CHAPTER 1 - Coumarins and Furanocoumarins. In: Vasil FCK, editor. Phytochemicals in Plant Cell Cultures. Oxford: Academic; 1988. p. 3–21.View ArticleGoogle Scholar
- Wink M. Introduction: Biochemistry, physiology and ecological functions of secondary metabolites. In: Annual Plant Reviews Volume 40: Biochemistry of Plant Secondary Metabolism. Chichester: Wiley-Blackwell; 2010. p. 1–19.View ArticleGoogle Scholar
- Bolwell GP, Bozak K, Zimmerlin A. Plant cytochrome P450. Phytochemistry. 1994;37(6):1491–506.View ArticlePubMedGoogle Scholar
- Halkier BA. Catalytic reactivities and structure/function relationships of cytochrome P450 enzymes. Phytochemistry. 1996;43(1):1–21.View ArticleGoogle Scholar
- Larbat R, Kellner S, Specker S, Hehn A, Gontier E, Hans J, Bourgaud F, Matern U. Molecular cloning and functional characterization of psoralen synthase, the first committed monooxygenase of furanocoumarin biosynthesis. J Biol Chem. 2007;282(1):542–54.View ArticlePubMedGoogle Scholar
- Larbat R, Hehn A, Hans J, Schneider S, Jugdé H, Schneider B, Matern U, Bourgaud F. Isolation and functional characterization of CYP71AJ4 encoding for the first P450 monooxygenase of angular furanocoumarin biosynthesis. J Biol Chem. 2009;284(8):4776–85.View ArticleGoogle Scholar
- Khan A, Kunesch G, Chuilon S, Ravisé A. Structure and biological activity of xanthyletin : a new phytoalexin of CITRUS. Fruits. 1985;40(12):807–11.Google Scholar
- Beier RC, Ivie GW, Oertli EH. Linear furanocoumarins and graveolone from the common herb parsley. Phytochemistry. 1994;36(4):869–72.View ArticleGoogle Scholar
- Tomer E, Goren R, Monselise SP. Isolation and identification of seselin in Citrus roots. Phytochemistry. 1969;8(7):1315–6.View ArticleGoogle Scholar
- Borges F, Roleira F, Milhazes N, Santana L, Uriarte E. Simple coumarins and analogues in medicinal chemistry: occurrence, synthesis and biological activity. Curr Med Chem. 2005;12(8):887–916.View ArticlePubMedGoogle Scholar
- Poulsen S-A, Davis R. Natural products that inhibit carbonic anhydrase. In: Frost SC, McKenna R, editors. Carbonic Anhydrase: Mechanism, Regulation, Links to Disease, and Industrial Applications. Netherlands: Springer; 2014. p. 325–47.View ArticleGoogle Scholar
- Berenbaum MR, Zangerl AR. Chemical phenotype matching between a plant and its insect herbivore. Proc Natl Acad Sci U S A. 1998;95(23):13743–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Karamat F, Olry A, Munakata R, Koeduka T, Sugiyama A, Paris C, Hehn A, Bourgaud F, Yazaki K. A coumarin-specific prenyltransferase catalyzes the crucial biosynthetic reaction for furanocoumarin formation in parsley. Plant J. 2014;77(4):627–38.View ArticlePubMedGoogle Scholar
- Hehmann M, Lukacin R, Ekiert H, Matern U. Furanocoumarin biosynthesis in Ammi majus L. Cloning of bergaptol O-methyltransferase. European journal of biochemistry / FEBS. 2004;271(5):932–40.View ArticlePubMedGoogle Scholar
- Larbat R. Contribution à l’étude des P450 impliqués dans la biosynthèse des furocoumarines. France: Unité Mixte de Recherche Agronomie et Environnement (UMR INPL-INRA); 2006.Google Scholar
- Caspi R, Foerster H, Fulcher CA, Kaipa P, Krummenacker M, Latendresse M, Paley S, Rhee SY, Shearer AG, Tissier C, et al. The MetaCyc Database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nucleic Acids Res. 2008;36 suppl 1:D623–31.PubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.View ArticlePubMedPubMed CentralGoogle Scholar
- Li B, Dewey C. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC bioinformatics. 2011;12(1):323.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52(5):696–704.View ArticlePubMedGoogle Scholar
- Le SQ, Gascuel O. An improved general amino acid replacement matrix. Mol Biol Evol. 2008;25(7):1307–20.View ArticlePubMedGoogle Scholar
- Anisimova M, Gascuel O. Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol. 2006;55(4):539–52.View ArticlePubMedGoogle Scholar
- Arnold K, Bordoli L, Kopp J, Schwede T. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006;22(2):195–201.View ArticlePubMedGoogle Scholar
- Kiefer F, Arnold K, Kunzli M, Bordoli L, Schwede T. The SWISS-MODEL Repository and associated resources. Nucleic Acids Res. 2009;37(Database issue):D387–92.View ArticlePubMedGoogle Scholar
- Guex N, Peitsch MC. SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis. 1997;18(15):2714–23.View ArticlePubMedGoogle Scholar
- Nischal L, Mohsin M, Khan I, Kardam H, Wadhwa A, Abrol YP, Iqbal M, Ahmad A. Identification and comparative analysis of microRNAs associated with low-N tolerance in rice genotypes. PLoS ONE. 2012;7(12), e50261.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramawat KG, Dass S, Mathur M. The chemical diversity of bioactive molecules and therapeutic potential of medicinal plants. In: Ramawat KG, editor. Herbal Drugs: Ethnomedicine to Modern Medicine. Berlin Heidelberg: Springer; 2009. p. 7–32.View ArticleGoogle Scholar