Skip to main content

Phylogenetic reconstruction using secondary structures of Internal Transcribed Spacer 2 (ITS2, rDNA): finding the molecular and morphological gap in Caribbean gorgonian corals



Most phylogenetic studies using current methods have focused on primary DNA sequence information. However, RNA secondary structures are particularly useful in systematics because they include characteristics, not found in the primary sequence, that give "morphological" information. Despite the number of recent molecular studies on octocorals, there is no consensus opinion about a region that carries enough phylogenetic resolution to solve intrageneric or close species relationships. Moreover, intrageneric morphological information by itself does not always produce accurate phylogenies; intra-species comparisons can reveal greater differences than intra-generic ones. The search for new phylogenetic approaches, such as by RNA secondary structure analysis, is therefore a priority in octocoral research.


Initially, twelve predicted RNA secondary structures were reconstructed to provide the basic information for phylogenetic analyses; they accorded with the 6 helicoidal ring model, also present in other groups of corals and eukaryotes. We obtained three similar topologies for nine species of the Caribbean gorgonian genus Eunicea (candelabrum corals) with two sister taxa as outgroups (genera Plexaura and Pseudoplexaura) on the basis of molecular morphometrics of ITS2 RNA secondary structures only, traditional primary sequence analyses and maximum likelihood, and a Bayesian analysis of the combined data. The latter approach allowed us to include both primary sequence and RNA molecular morphometrics; each data partition was allowed to have a different evolution rate. In addition, each helix was partitioned as if it had evolved at a distinct rate. Plexaura flexuosa was found to group within Eunicea; this was best supported by both the molecular morphometrics and combined analyses. We suggest Eunicea flexuosa (Lamouroux, 1821) comb. nov., and we present a new species description including Scanning Electron Microscopy (SEM) images of morphological characteristics (sclerites). Eunicea flexuosa, E. pallida, E. laxispica and E. mammosa formed a separate clade in the molecular phylogenies, and were reciprocally monophyletic with respect to other Eunicea (Euniceopsis subgenus, e.g. E. tourneforti and E. laciniata) in the molecular morphometrics tree, with the exception of E. fusca. Moreover, we suggest a new diagnostic character for Eunicea, also present in E. flexuosa: middle layer sclerites > 1 mm in length.


ITS2 was a reliable sequence for intrageneric studies in gorgonian octocorals because of the amount of phylogenetic signal, and was corroborated against morphological characters separating Eunicea from Plexaura. The ITS2 RNA secondary structure approach to phylogeny presented here did not rely on alignment methods such as INDELS, but provided clearly homologous characters for partition analysis and RNA molecular morphometrics. These approaches support the divergence of Eunicea flexuosa comb. nov. from the outgroup Plexaura, although it has been considered part of this outgroup for nearly two centuries because of morphological resemblance.


Millions of RNA sequences (as DNA) have been deposited in public databases such as Genbank (NCBI), and over one hundred thousand of these correspond to genes such as 18S that are widely used in phylogenetic studies. Although functional RNA sequences fold themselves to form secondary structures that depend sensitively on the primary sequence, only a handful of phylogenetic studies have taken advantage of this secondary information. Most phylogenetic studies using current methods have focused on primary sequence information, but RNA secondary structures are particularly useful in systematics because they include characteristics, not found in the primary sequence, that give "morphological" information for reconstructing the complete tree of life [1]. Ribosomal gene secondary structures have conserved zones that make it easier to reconstruct the structures of unknown homologous RNA sequences by direct comparison. Although molecular phylogenetic inference relies heavily on single-copy genes, ribosomal genes that show concerted evolution still offer one of the best alternatives in the genome. The aim of this paper is to demonstrate the value of RNA secondary structures in improving and enhancing phylogenetic inference for a group of marine invertebrates.

Despite progress towards the resolution of octocoral phylogeny using mitochondrial gene sequences [2], there is no consensus about a sequence region that provides enough resolution for closely-related species relationships. The most variable regions known so far in octocorals are the Internal Transcribed Spacers (ITSs), which have provided good resolution among a few soft coral species that would be very difficult to distinguish morphologically because their characters show excessive homoplasy [3, 4]. However, the ITS region shows saturation and in some cases many INDELS (insertions-deletions) in an alignment, which can produce misleading results for high taxonomic ranks [5]. Opportunely, analyses of the predicted RNA secondary structures of some genes such as 16S [4] and ITS2 [6] have provided an alternative way of increasing the phylogenetic signal of RNA sequences for octocorals. Here, we studied the predicted RNA secondary structures of ITS2 to resolve close-species relationships within the gorgonian genus Eunicea (candelabrum corals).

ITS2 (Internal Transcribed Spacer 2) sequences, found in the tandem arrays of the nuclear ribosomal RNA between the 5.8S and 28S genes, have not been considered useful for molecular systematics in certain invertebrate groups mainly because of their excessive INDELS and/or intragenomic variation [7, 8]. Nonetheless, they seem to afford a promising nuclear region for lower eukaryotes, particularly at the intra-familial level and down to closely-related species, where other known regions are nearly invariant [912]. In scleractinian corals, except for the genus Acropora, ITS sequences give reliable phylogenetic information especially if the predicted RNA secondary structures are compared [6]. These sequences vary considerably within octocoral genera, providing information for classifying a recently-described species, Alaskagorgia aleutiana Sánchez and Cairns, within the Plexauridae; this was uncertain on the basis of morphological characters, but the results were consistent with mitochondrial DNA analysis [13, 5]. Variation was also found within the Gorgoniidae family; ITS2 secondary structures retained and classified different gorgonian genera [5]. Perhaps they are among the most variable genomic regions in these organisms, while mitochondrial DNA and other nuclear sequences are highly conserved [10, 11, 14]. Unfortunately, ITS alignment always produces multiple and variable INDELS depending on the different gap opening and/or extension penalties (e.g. ClustalW), which make phylogenetic reconstruction somewhat unreliable.

The known predicted ITS2 RNA secondary structure is a model with a common core in most eukaryotes [15]. It is a large-scale marker that is not limited to a specific taxonomic level and has even been used in mega-systematics analysis [15]. Among other advantages, molecular morphometrics makes it possible to identify characters informative about parsimony that are not found in the primary sequence [16, 17], and to correct primary alignments on the basis of secondary structure information [3]. Could the predicted ITS2 RNA secondary structures help to overcome the great difficulties in the taxonomy of lower invertebrates such as octocorals?

Caribbean octocorals are very abundant in coral reefs with up to 60 species in a single location, and many species and genera show few or no morphological differences [18]. Gorgonian octocorals from the genus Eunicea Lamouroux, 1816 (Cnidaria; Octocorallia; Alcyonacea [=Gorgonacea]; Plexauridae), known as "sea candelabrum", are endemic to the Caribbean region (Tropical Western Atlantic Ocean excluding the Brazilian region). There are 13 valid species, although more than 30 species have been described in, or assigned to, Eunicea [18, 19]. No phylogeny for the Eunicea genus has been constructed on the basis of morphological or molecular characters, and this could be important for testing the monophyly of the internal subgenera [19] and other morphological characters that incur the risk of confusion with Plexaura octocorals [18]. In this study, in order to differentiate between Eunicea and Plexaura gorgonian corals using molecular and morphological data, we present a phylogeny based on the predicted secondary structures of the ITS2 sequences (molecular morphometrics) and contrast it with traditional molecular systematics based on primary DNA alignment and Bayesian analysis.


The ITS2 nuclear region ranged from 252 bp in Eunicea laciniata to a minimum length of 185 bp in Plexaura kuna. Twelve new predicted RNA secondary structures were constructed for this region, all following the 6 helicoidal ring model except for E. pallida, which lacked the sixth helix (Fig. 1I). Two major common features were present in these structures: (1) 5' G-C, R-R, G-C, G-C, bulge on helix IV and (2) 5' GUGC, bulge, CAAGG with its complementary pair base on helix V (Fig. 1).

Figure 1
figure 1

Predicted ITS2 RNA secondary structures for 9 species of the genus Eunicea and their structure formation enthalpies according to MFOLD: (A) Eunicea flexuosa, numbering represents characters used for molecular morphometry [See additional file 2], ΔG = -58.96 kcal/mole; (B) Eunicea laciniata, ΔG = -104.17 kcal/mole; (C) Eunicea sp. 1, ΔG = -104.17 kcal/mole; (D) Eunicea sp. 2, ΔG = -71.45 kcal/mole; (E) E. laxispica, ΔG = -62.67 kcal/mole; (F) E. fusca, ΔG = -65.96 kcal/mole; (G) E. tourneforti, ΔG = -69.67 kcal/mole; (H) E. mammosa, ΔG = -67.66 kcal/mole; (I) E. pallida, ΔG = -69.53 kcal/mole.

Three different phylogenetic analyses were obtained from the ITS2 sequence alignment. The molecular morphometrics analysis relied on 42 parsimony-informative characters and resulted in one most parsimonious tree using ordered Wagner parsimony (Length = 270; Retention index = 0.57). Both bootstrapping and decay indices showed strong support for Plexaura and Eunicea as reciprocally monophyletic groups (95% and 9, respectively) with the exception of E. flexuosa (new combination from Plexaura flexuosa, see details below), which was clearly a derived species of Eunicea (Fig. 2A). The topology based only on the predicted RNA secondary structure of the ITS2 region resolved most relationships among the species studied (Fig. 2A). Maximum likelihood analysis using traditional primary sequence alignment retained one tree with the assumed model (TVMef+G) selected by AIC. The tree searches had six substitution rates (A-C, 1.2431; A-G, 3.0064; A-T, 1.0809; C-G, 0.3376; C-T, 3.0064; and G-T, 1.0000), not assuming the proportion of invariant sites, and a gamma-shape parameter (1.413); node support was examined by 1000 bootstrap replicates (Fig. 2B). Bayesian analysis of the same alignment retained the same topology and supported the same branches (data not shown). This last topology differed in two main respects from the molecular morphometrics tree: the placement of E. mammosa within the outgroup P. kuna Lasker, Cofforth and Kim, with low bootstrap and Bayesian clade credibility values (< 60%), and of E. fusca Duchassaing and Michelotti as a derived instead of a basal species (Fig. 2A–B). The Bayesian inference in which helices were to evolve independently at different rates in a mixed matrix secondary structure was similar to the molecular morphometrics tree though less resolved within Eunicea (Fig. 2C). Eunicea flexuosa (Lamouroux), E. pallida Garcia-Parrado and Alcolado, E. laxispica (Lamarck, 1815) and E. mammosa Lamouroux formed a separate clade in the molecular phylogenies (except for E. mammosa in the maximum likelihood tree), most probably the Eunicea subgenus [19], and were reciprocally monophyletic with respect to other Eunicea (Euniceopsis subgenus : E. laciniata Duchassaing & Michelotti, E. tourneforti Milne Edwards and Haime, and Eunicea spp. 1–2), except for E. fusca in the molecular morphometrics tree (Fig. 2). The predicted ITS2 RNA secondary structure models showed that E. mammosa, E. pallida, E. laxispica and E. flexuosa exhibited a long helix IV (Fig. 1A, E, H, I) that could be a sinapomorphic feature of the Eunicea subgenus. Plexaura spp., in contrast, lacked that helix [see Additional file 1].

Figure 2
figure 2

Phylogenetic analyses of Eunicea spp. A. Most parsimonious phylogram using molecular morphometrics data [See additional file 2]. Numbers represent bootstrap values of 1,000 replicates/decay index (Bremer support). B. Maximum Likelihood phylogram of the primary ITS2 alignment with the assumed model (TVMef+G) selected by AIC; tree searches had six substitution rates (A-C, 1.2431; A-G, 3.0064; A-T, 1.0809; C-G, 0.3376; C-T, 3.0064; and G-T, 1.0000), an assumed proportion of invariant sites (= 0), and a gamma-shape parameter (1.4130). Bayesian probabilities and bootstrap values of 1,000 replicates are separated by a slash. C. Hypothetical Bayesian analysis phylogeny based on the combined predicted secondary structure and secondary sequence alignment data (e.g. including secondary structure partitions) of the ITS2 region. Inset photo: Eunicea flexuosa (Lamouroux, 1821) comb. nov., San Salvador, Bahamas, 12 m (1999).

Systematic considerations

Genus Eunicea Lamouroux,1816

Diagnosis: Plexaurids with spindles of the middle ring layer averaging >> 1 mm in length and greater than 5 mm in some species (e.g. Fig. 3A–B). The lengths of the middle layer spindles in Plexaura are noticeably greater than the surface layer club sclerites in Eunicea (e.g. Fig. 3A–D).

Figure 3
figure 3

SEM images of candelabrum coral sclerites. A-D. Sclerites (ornamented spindles, scale 0.1 mm) from the middle layer: A. Eunicea flexuosa (Lamouroux,1821) comb. nov. (scale 0.1 mm). B. E. laciniata Duchassaing & Michelotti C. Plexaura kuna Lasker et al. (1996). D. Eunicea sp. 1. E-F. Club sclerites from the coenenchyme surface. E. P. kuna (scale 0.1 mm). F. E. flexuosa (scale 0.04 mm).

Eunicea flexuosa(Lamouroux,1821) comb. nov.

(Fig. 3)

Plexaura flexuosa Lamouroux 1821: 135 [20] (Cuba)-Holotype: Mus. Hist. Nat. Paris, Paris.

Plexaura flexuosa: Bayer 1961:104 [19].


Eunicea with colonies exhibiting many forms ranging from bushy, candelabrum, and dense stands of short branches. Dichotomous pattern of branching (Fig. 2 inset). Polyp aperture as a small prominence of the lower lip or a well-developed calyx with prominent lower lip. Polyp armature, without collar, consisting of ornate rods between 0.05 and 0.3 mm long. Axial layer with capstans and spindle sclerites with diverse arrangements and ornamentations; purple coloration, 0.01–0.26 mm long. Middle layer sclerites robust spindles up to 4.5 mm (Fig. 3A). Surface layer sclerites with foliated club sclerites exhibiting fused lobules in the head, occasionally serrated; 0.22 mm long (Fig. 3F).

Distribution and habitat

Widespread in the Tropical Western Atlantic. All kinds of reef or rocky environments with hard ground and some water movement (moderate to rough) between 0.5 and 30 m depth.

Material examined

ICN-MHN-CO-076 (Instituto de Ciencias Naturales, U. Nacional, Colombia) (U-58), Isla Narsa, Capurganá, Gulf of Urabá, Caribbean coast of Colombia, col. J.A. Sánchez, October 2, 1995, slope edge 12 m. ICN-MHN-CO-087 (F-19) Isla Fuerte, bajo "el bobito", Caribbean coast of Colombia, col. J.A. Sánchez, October 3, 1995, marginal terrace 10 m. ICN-MHN-CO-089 (Pfl-R Tayrona National Park, Arrecifes, "La Piscinita", Santa Marta, Colombian Caribbean, col. J.A. Sánchez, rock barrier 1 m. ICN-MHN-CO-092 (AR-1), Tayrona National Park, Arrecifes, "La Piscinita", Santa Marta, Colombian Caribbean, col. J.A. Sánchez, rock barrier 1 m. USNM (Smithsonian Institution, National Museum of Natural History, USA) 97680 (12), Tesoro Island, Rosario Islands, Caribbean coast of Colombia, col. J.A. Sánchez & A. Ramirez, September 16, 1992; fore-reef terrace, dead stands of Acropora palmata 7 m; INV-CNID-382, (AR-3), Tayrona National Park, Arrecifes, "La Piscinita", Santa Marta, Colombian Caribbean, col. J.A. Sánchez, rock barrier 1 m. USNM 51811, North Atlantic Ocean; Bermuda, Somerset Island, 5–6 m prof. col. Tucker T., 1 September 1960. USNM 14380, North Atlantic Ocean, Bahamas, New Providence, Island, Nassau, East end of Island, col. Nye W. Jr., 1886, Albatross R/V.


Traditionally, morphological identification of Eunicea and Plexaura species has been based mostly on microscopic characteristics [19]. Morphological variation within species can overlap with variation among species, so molecular identification has become an important alternative for octocorals in general. ITS2 has proven to be a good molecular marker, resolving octocoral species at different taxonomic levels [5]. In particular, the ITS2 sequences and predicted RNA secondary structures have afforded a new view of the phylogeny of the Eunicea genus. This approach supports the divergence of Eunicea flexuosa comb. nov. from Plexaura, although it has been considered part of this outgroup for nearly two centuries because of morphological resemblance. Nonetheless, the morphological observation of two subgenera in Eunicea [19] seemed congruent with the molecular results.

The ITS2 RNA secondary structures yielded homologous models that grouped conserved species features in the Eunicea subgenus (e.g. long helix IV), though it is important to note that these molecular structures can also produce homoplasies because of selective pressures acting on particular helix features such as the 6-helicoidal ring core, conserved among most eukaryotes. ITS2 structures are important for proper pre-rRNA maturation; any change in their helix conformation can affect the process of 25S production [21]. Each helix was partitioned as if it had evolved at a different rate and substitution bias, because they do not all have the same role in this dynamic conformation [22]. This approach resulted in a supported Bayesian inference phylogram very similar to the molecular morphometrics approach. The major difference between the topologies obtained by the Bayesian (e.g. molecular morphometrics using partitions) and primary sequence (e.g maximum likelihood) methodologies, was the placement of E. mammosa in and out of Eunicea, which could be due to the lack of phylogenetic signal and the likely saturation of ITS2 sequences after multiple alignment (e.g. [5]).

The combined molecular-morphological Bayesian approach data matrix has not been widely used, mainly because until recently the methods for combining datasets were only applied to parsimony-based analysis, which has been shown to be insufficient when the combined data imply different evolutionary scenarios [23]. This type of analysis, combining molecular data with "typical" morphological characteristics, has sometimes failed to yield a satisfactory model of character evolution [24]. Here, we combined the corrected primary alignment (e.g. aligning sequences from each helix only) with the molecular morphometrics data. Traditionally, some problems in explaining changes in character state arise from assuming that each character has exactly the same state at a particular time, in opposition to punctuated equilibria. In molecular morphometrics, the possibility of change is based in the number of nucleotides present in a particular helix, which can be readily identified. The other issue, corresponding to the comparison between character states, was resolved using only the correction from the parsimonious-informative sites [24]. Nevertheless, because there were inconsistencies in the placement of a few species (E. mammosa and E. fusca), this study needs to be extended, hopefully including all Eunicea representatives, in order to gain a better understanding of the systematics of this group as well as the evolution of their predicted ITS2 RNA secondary structures.


The molecular study of this genus, which is the most species-rich among shallow-water Caribbean gorgonian corals, is just beginning. Consequently, ITS2 RNA secondary structure analysis could be a valuable tool for distinguishing new species and completing Eunicea systematics: morphological phylogenetic reconstruction of octocorals is very difficult (e.g. [13]), and ITS2 secondary structure contains more information than the usual primary sequence alignment [5, 26].


ITS2 sequences

Sequences were obtained using primers designed by [26], which target the region between the 5' end of the 5.8S and the 3' end of the 28S ribosomal genes, containing the complete ITS2 (5.8S-436: 5'-AGCATGTCTGTCTGAGTGTTGG-3' and 28S-663: 5'-GGGTAATCTTGCCTGATCTGAG-3': numbers relate to the Alcyonium digitatum sequence, Genbank Acc. No. AF262347: [2]). DNA was extracted from material preserved in 95% ethanol using the DNeasy kit (Qiagen). Template DNA for sequencing was obtained from a combination of two PCR reactions containing 1 μl DNA template (1:50 dilution of genomic DNA extract), 2 units Taq polymerase (Promega), 3 μl 10X Buffer (Promega), MgCl2, 0.15 μM dNTP mixture and 0.16 μM of each primer in 56 μl total volume (completed with doubly-distilled water). PCR conditions were: one initial period of 2.0 min at 94°C followed by 30 cycles of 30 s at 94°C, 45 s at 56°C, 1.0 45 s at 72°C, and a final extension step for 5.0 min at 72°C. Genes were purified using the Edge Biosystems kit and sequenced using BigDye 3.1 (AB 3100, capillary electrophoresis automated sequencer). Consensus sequences were obtained automatically by assembling the two complementary DNA chromatograms (Sequencher software). Sequence information is available in Genbank (Accession numbers EF 490973-84).

Predicted ITS2 RNA secondary structures

Secondary structures for nine octocoral species were reconstructed by aligning their sequences (using Bioedit) [25] with homologous structures already published [26]. In addition, manual alignment was performed by visual homology to construct a Dedicated Comparative Sequence Editor (DCSE) format in order to perform different phylogenetic analyses. The DCSE format uses square brackets ([, ]) to delimit each helix, braces ({, }) for bulges and loops, and hyphens for gaps in the alignment; terminal loops appear separated by inverse brackets (] and [) (see Additional file 1). The acquired structures with restrictions and constraints were submitted in MFOLD [27]. RNA was folded at a fixed temperature of 37°C, and the structure chosen from different output files was the desired 6-helicoidal ring or the one with the highest negative free energy if various similar structures were obtained.

Phylogenetic analyses

Sixteen sequences of the ITS2 region of nine Eunicea species (Eunicea sp. 1, E. fusca, E. mammosa, E. tourneforti, E. laxispica, E. laciniata, E.pallida, Eunicea sp. 2 and E.flexuosa) and three corresponding to Pseudoplexaura crucis Bayer, Plexaura homomalla (Esper) and Plexaura kuna were analyzed. The manually-obtained alignment with secondary structure information was used to construct a matrix for cladistic analysis based on molecular morphometrics using geometrical features and base numbering [16]. Thus, secondary structure helices were numbered and treated like characters that vary depending on their base number. Bulges and internal loops were assigned the number of the helix followed by a letter; ending helix loops and separation segments were designated "i". Loop complements, helix complements and separation segments were indicated with apostrophes (e.g. loop 4a with its complement 4a'). The nucleotides were counted for each character and the number scaled as discrete character states ranging from 0 to 8 or by subtracting the lesser count from the rest (each nucleotide addition counted as a new character, 0 being the character with the lowest nucleotide number) [See additional file 2]. A semi-exhaustive search using ordered (Wagner) maximum parsimony and the branch and bound algorithm was executed in PAUP* [29]. Branch support was given by estimating the decay index (Bremer support) and using 1000 bootstrap replicates in PAUP* [30].

Initially, the sequences were aligned in BioEdit [25] using ClustalW multiple alignment [28] with the default gap and extension penalties used by this program. The primary alignment was submitted to the best-fit model from ModelTest [31] obtained by the Akaike Information Criterion (AIC) in order to perform a search with the branch and bound algorithm in PAUP* [29]. INDELS were treated as missing data and not used in any of the analyses. Bayesian inference of phylogeny was done using MrBayes [32], Bayesian-estimated likelihood (settings according to MrModeltest), 10 million Monte Carlo Markov chain generations (Bayesian-Monte Carlo simulation by MrBayes sampling every 100 simulation, burn-in 10000).

A third phylogenetic analysis was carried using the Bayesian approach with combined datasets [33] using Mrbayes 3.1 [32]. In this approach, each data partition is allowed to have a different evolution rate (using the option prset ratepr = variable). The molecular morphometrics data were treated as DATATYPE = standard mode and no evolutionary model was selected, but a different rate of change for each helix was set (γ). In addition, each helix was partitioned as if it had evolved at a different rate. The model for the sequence dataset was searched using the software Mrmodeltest. Two separate chains were run with 10 million generations each to test for the convergence of the parameter estimates from two different departure points of the search.

Morphological characters

Scanning Electron Microscope (SEM) characters were established following the methods from [13] for the number of specimens given in the "systematics considerations" section.


  1. Caetano-Anollés G: Tracing the evolution of RNA structure in ribosomes. Nucleic Acids Res. 2002, 30: 2575-2587. 10.1093/nar/30.11.2575.

    Article  PubMed Central  PubMed  Google Scholar 

  2. McFadden CS, Donahue R, Hadland BK, Weston R: A molecular phylogenetic analysis of reproductive trait evolution in the soft coral genus Alcyonium. Evolution. 2001, 55: 54-67.

    Article  CAS  PubMed  Google Scholar 

  3. Sánchez JA, Lasker HR, Taylor DJ: Phylogenetic analyses among octocorals (Cnidaria) according to mitochondrial and nuclear DNA sequences (lsu-rRNA 16S, and ssu-rRNA 18S) support two convergent clades of branching gorgonians. Mol Phylogenet Evol. 2003, 29: 31-42. 10.1016/S1055-7903(03)00090-3. (DOI:10.1016/S1055-7903(03)00090-3)

    Article  Google Scholar 

  4. Sánchez JA, Mcfadden CS, France SC, Lasker HR: Molecular Phylogenetic analyses of shallow-water Caribbean octocorals. Marine Biol. 2003, 142: 975-987. (DOI10.1007/s00227-003-1018-7)

    Google Scholar 

  5. Aguilar C, Sánchez JA: Contribution of ITS2 sequences and predicted RNA secondary structures to octocoral systematics. Bull Marine Sci.

  6. Chen CA, Chang CC, Wei NV, Chen CH, Lein YT, Lin HE, Dai CF, Wallace CC: Secondary structure and phylogenetic utility of the ribosomal internal transcribed spacer 2 (ITS2) in Scleractinian Corals. Zool Stud. 2004, 43: 759-771.

    CAS  Google Scholar 

  7. Fabry S, Kohler A, Coleman AW: Intraspecies analysis: Comparison of ITS sequence data and gene intron sequence data with breeding data for a worldwide collection of Gonium pectorale. J Mol Evol. 1999, 48: 94-101. 10.1007/PL00006449.

    Article  CAS  PubMed  Google Scholar 

  8. Harris DJ, Crandall KA: Intragenomic variation within ITS1 and ITS2 of freshwater crayfishes (Decapoda: Cambaridae): Implications for phylogenetic and microsatellite studies. Mol Biol Evol. 2000, 17: 284-291.

    Article  CAS  PubMed  Google Scholar 

  9. Knowlton N: Molecular genetic analyses of species boundaries in the sea. Hydrobiología. 2000, 420: 73-90. 10.1023/A:1003933603879.

    Article  CAS  Google Scholar 

  10. Berntson EA, Bayer FM, McArthur AG, France SG: Phylogenetic relationships within the Octocorallia (Cnidaria : Anthozoa) based on nuclear 18S rRNA sequences. Marine Biol. 2001, 138: 235-246. 10.1007/s002270000457.

    Article  CAS  Google Scholar 

  11. LaJeunesse TC: Investigating the biodiversity, ecology, and phylogeny of endosymbiotic dinoflagellates in the genus Symbiodinium using the ITS region: in search of a "species" level marker. J Phycol. 2001, 37: 866-880. 10.1046/j.1529-8817.2001.01031.x.

    Article  CAS  Google Scholar 

  12. Shearer TL, van Oppen MJH, Romano SL, Wörheide G: Slow mitochondrial DNA sequence evolution in the Anthozoa (Cnidaria). Mol Ecol. 2002, 11: 2475-2487. 10.1046/j.1365-294X.2002.01652.x.

    Article  CAS  PubMed  Google Scholar 

  13. Sánchez JA, Cairns SD: An unusual new gorgonian coral from the Aleutian Islands, North Pacific, Alaska. Zool Med. 2004, 78: 265-274.

    Google Scholar 

  14. France SC, Hoover LL: DNA sequences of the mitochondrial COI gene have low levels of divergence among deep-sea octocorals (Cnidaria: Anthozoa). Hydrobiologia. 2002, 471: 149-155. 10.1023/A:1016517724749.

    Article  Google Scholar 

  15. Schultz J, Maisel S, Gerlach D, Müller T, Wolf M: A common core of secondary structure of the internal transcribed spacer 2 (ITS2) throughout the Eukaryota. RNA. 2005, 11: 361-364. 10.1261/rna.7204505.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Billoud B, Guerrucci MA, Masselot M, Deutsch JS: Cirripede phylogeny using a novel approach: molecular morphometrics. Mol Biol Evol. 2000, 17: 1435-1445.

    Article  CAS  PubMed  Google Scholar 

  17. Swain TD, Taylor DJ: Structural rRNA characters support monophyly of raptorial limbs and paraphyly of limb specialization in water fleas. Proc Roy Soc London B. 2003, 270: 887-896. 10.1098/rspb.2002.2297.

    Article  CAS  Google Scholar 

  18. Sánchez JA, Wirshing H: A field key to the identification of zooxanthellate octocorals from the Caribbean and Western Atlantic. Caribbean J Sci. 2005, 41: 508-522.

    Google Scholar 

  19. Bayer FM: The shallow water Octocorallia of the West Indian region. Stud Fauna Curaçao. 1961, 12: 1-373.

    Google Scholar 

  20. Lamouroux JVF: Exposition méthodique des genres de l'ordre des Polypiers, avec description et celles des principales espèces, figurées dans 84 planches; les 63 premières appartenant a l'Histoire des Zoophytes d'Ellis et Solander Paris. 1821-

  21. Cote CA, Peculis BA: Role of the ITS2-proximal item and evidence for indirect recognition of processing sites in pre-rRNA processing in yeast. Nuc Acids Res. 2001, 29: 2106-2116. 10.1093/nar/29.10.2106.

    Article  CAS  Google Scholar 

  22. Cote CA, Greer CL, Peculis BA: Dynamic conformational model for the role of ITS2 in pre-rRNA processing in yeast. RNA. 2002, 8: 786-797. 10.1017/S1355838202023063.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Pollock DD, Bruno WJ: Assessing an unknown evolutionary process: Effect of increasing site-specific knowledge through taxon addition. Mol Biol Evol. 2000, 17: 1854-1858.

    Article  CAS  PubMed  Google Scholar 

  24. Lewis PO: A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001, 50: 913-925. 10.1080/106351501753462876.

    Article  CAS  PubMed  Google Scholar 

  25. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp. 1999, 41: 95-98.

    CAS  Google Scholar 

  26. Aguilar C, Sánchez JA: Phylogenetic hypotheses of gorgoniid octocorals according to ITS2 and their predicted RNA secondary structures. Mol Phylogenet Evol. 2007, 43: 774-786. 10.1016/j.ympev.2006.11.005. doi:10.1016/j.ympev.2006.11.005.

    Article  CAS  PubMed  Google Scholar 

  27. Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucl Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (and other methods). Version 4.0b10. 2002, Sunderland, MA: Sinauer Associates

    Google Scholar 

  30. Hillis DM, Bull JJ: An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst Biol. 1993, 42: 182-192. 10.2307/2992540.

    Article  Google Scholar 

  31. Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.

    Article  CAS  PubMed  Google Scholar 

  32. Huelsenbeck JP, Ronquist F: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.

    Article  CAS  PubMed  Google Scholar 

  33. Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL: Bayesian phylogenetic analysis of combined data. Syst Biol. 2004, 53: 47-67. 10.1080/10635150490264699.

    Article  PubMed  Google Scholar 

Download references


Special thanks to the Universidad de los Andes (Faculty of Sciences and Department of Biological Sciences) for resources, Howard Lasker (SUNY, University at Buffalo) for field collections in the Bahamas, Klaus Ruetzler (CCRE, Smithsonian) for field collections in Carrie Bow Cay, Belize, and a postdoctoral fellowship (J.A. Sánchez, 2002–2003) from the Smithsonian Institution-NMNH (Liz Zimmer, Ken Schallop, Lee Weight and Stephen D. Cairns) for initial molecular analyses. Colleagues and students from BIOMMAR provided helpful discussions of the paper.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Juan A Sánchez.

Additional information

Authors' contributions

AG performed most of the phylogenetic analyses and had a major role in writing the manuscript. CA reconstructed the ITS2 RNA secondary structure, collaborated in the phylogenetic analyses and had a major role in writing the manuscript. JAS collected all the specimens studied, obtained the DNA sequences in the laboratory and performed the SEM analyses of sclerites; he also contributed to write the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: DCSE Alignment. Sequences alignment in DCSE format showing secondary structures of fourteen species of octocorals. (DOC 60 KB)

Additional file 2: Molecular morphometrics. Molecular morphometrics species/characters matrix (DOC 144 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Grajales, A., Aguilar, C. & Sánchez, J.A. Phylogenetic reconstruction using secondary structures of Internal Transcribed Spacer 2 (ITS2, rDNA): finding the molecular and morphological gap in Caribbean gorgonian corals. BMC Evol Biol 7, 90 (2007).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: