Development of an embryonic skeletogenic mesenchyme lineage in a sea cucumber reveals the trajectory of change for the evolution of novel structures in echinoderms
© McCauley et al.; licensee BioMed Central Ltd. 2012
Received: 18 June 2012
Accepted: 22 June 2012
Published: 9 August 2012
The mechanisms by which the conserved genetic “toolkit” for development generates phenotypic disparity across metazoans is poorly understood. Echinoderm larvae provide a great resource for understanding how developmental novelty arises. The sea urchin pluteus larva is dramatically different from basal echinoderm larval types, which include the auricularia-type larva of its sister taxon, the sea cucumbers, and the sea star bipinnaria larva. In particular, the pluteus has a mesodermally-derived larval skeleton that is not present in sea star larvae or any outgroup taxa. To understand the evolutionary origin of this structure, we examined the molecular development of mesoderm in the sea cucumber, Parastichopus parvimensis.
By comparing gene expression in sea urchins, sea cucumbers and sea stars, we partially reconstructed the mesodermal regulatory state of the echinoderm ancestor. Surprisingly, we also identified expression of the transcription factor alx1 in a cryptic skeletogenic mesenchyme lineage in P. parvimensis. Orthologs of alx1 are expressed exclusively within the sea urchin skeletogenic mesenchyme, but are not expressed in the mesenchyme of the sea star, which suggests that alx1 + mesenchyme is a synapomorphy of at least sea urchins and sea cucumbers. Perturbation of Alx1 demonstrates that this protein is necessary for the formation of the sea cucumber spicule. Overexpression of the sea star alx1 ortholog in sea urchins is sufficient to induce additional skeleton, indicating that the Alx1 protein has not evolved a new function during the evolution of the larval skeleton.
The proposed echinoderm ancestral mesoderm state is highly conserved between the morphologically similar, but evolutionarily distant, auricularia and bipinnaria larvae. However, the auricularia, but not bipinnaria, also develops a simple skelotogenic cell lineage. Our data indicate that the first step in acquiring these novel cell fates was to re-specify the ancestral mesoderm into molecularly distinct territories. These new territories likely consisted of only a few cells with few regulatory differences from the ancestral state, thereby leaving the remaining mesoderm to retain its original function. The new territories were then free to take on a new fate. Partitioning of existing gene networks was a necessary pre-requisite to establish novelty in this system.
KeywordsSea cucumber Echinoderm Evolution of novelty Co-option Skeletogenesis Alx1
Gene regulatory network
Whole mount in situ hybridization.
It has been known for many years that most animals rely on the use of the same basic set of regulatory genes during development . The great challenge is to determine how novel structures are specified by orthologous genes. Gene regulatory networks (GRNs) for novel structures can presumably be built “from scratch” (also termed de novo construction), or be re-wired from existing GRNs, which can include the co-option and fusion of GRN subcircuits. In particular, co-option, that is, the reuse of ancestral regulatory programs in novel contexts, might be a major mechanism through which novelties arise during evolution (reviewed in ). For example, this process is thought to have occurred in the formation of butterfly eyespots, patterning appendage outgrowth, the formation of novel insect appendages, and in the evolution of gnathostome oral teeth, among many others [3–8]. Relevant to this study, it has also been suggested that the sea urchin larval skeleton might have arisen via co-option of an adult echinoderm skeletogenic program . However, it is essentially unknown how an ancestral GRN can be rewired to accommodate the acquisition of a co-opted regulatory program in order to produce a viable new state.
The larvae of echinoderms provide a rich source of morphological variation for understanding the evolution of novelty. Broadly, echinoderms have two types of feeding larvae: the pluteus-like larvae of sea urchins and brittle stars, and the auricularia-like larvae of sea cucumbers and sea stars . Interestingly, most sea lilies, which are the earliest-branching echinoderms , have a secondarily-derived non-feeding larva, but at least one species is known to go through an auricularia-like stage . Due to the similarity between the hemichordate tornaria larva and the echinoderm auricularia-type larva, this is considered to be the basal form , possibly to the whole deuterostome clade. Intriguingly, larval morphology does not track with phylogenetic relationships among the echinoderms. Sea urchins and sea cucumbers are sister taxa, while brittle stars either form the outgroup to this clade or are the sister group to sea stars, though the exact phylogeny remains unclear due to the rapid divergence of these classes [11, 14, 15]. Sea lilies are thought to be the most basal echinoderms. Thus, the derived plutei of sea urchins and brittle stars do not group together, but instead are scattered among the ancestral auricularia-like larvae in the echinoderm phylogenetic tree.
There are many differences between echinoderm larval forms, but perhaps the most dramatic and obvious is the larval skeleton that is found in the sea urchin plutei that is not at all present in sea star larvae. The larval skeleton provides the structure that gives the sea urchin the dramatic pluteus morphology. The larval skeleton in modern sea urchins (non-irregular euechinoids, hereafter referred to as simply “sea urchins”) forms from the skeletogenic mesenchyme (SM) that arises from micromeres at the vegetal-most pole (reviewed in ). The remaining mesoderm goes on to produce mesenchymal blastocoelar cells, as well as pigment cells, epithelial coeloms and muscle [17, 18]. The specification of the mesodermal territories is extraordinarily well known in sea urchins, and the SM in particular has, perhaps, the best characterized GRN for any developmental system . Each territory expresses a unique complement of transcription factors during blastula stages, though there is considerable overlap in the expression of individual factors among these lineages. alx1, erg, ets1, foxn2/3, hex, tbr and tgif are expressed in the SM lineage [20–27]. The dorsally-restricted pigment cell precursors express gcm, foxn2/3 and gata4/5/6[26, 28, 29], while the ventrally-localized presumptive blastocoelar cells express erg, gata1/2/3, foxn2/3, ets1 and gata4/5/6[25, 26, 29, 30]. Each of these mesodermal cell types goes on to exhibit distinct behavior during gastrulation. The SM ingresses into the blastocoel at the onset of gastrulation and migrates to form a characteristic ring at the dorsal side of the embryo, which terminates in two ventro-lateral clusters. During late gastrulation, cells of the SM fuse and begin secreting the larval skeleton. The pigment cells are the next population of mesenchyme to ingress from the tip of the archenteron during early gastrulation; they migrate to the dorsal ectoderm, into which they intercalate and begin secreting pigment. Blastocoelar cells do not ingress until the late gastrula stage .
Previous studies have revealed that orthologs of many of the regulatory genes expressed in the sea urchin mesodermal territories are expressed in the presumptive mesoderm of sea stars [32–34], which will give rise to blastocoelar cells and epithelial mesoderm, but not an SM lineage. We, therefore, undertook an analysis of the expression of orthologous transcription factors in the embryos of a phylogenetic intermediate, the sea cucumber, Parastichopus parvimensis, to better correlate changes in regulatory gene expression with morphological novelty. This study is the first comprehensive gene expression analysis performed in any species of sea cucumber and identifies a highly conserved pleisiomorphic mesodermal regulatory state. We show for the first time that sea cucumbers develop a mesodermally-derived skeletogenic mesenchyme from the vegetal pole of blastulae, although the resulting larval skeleton is much less morphologically complex than that of sea urchins. Our results suggest that the evolution of larval skeletogenic cells in echinoderms occurred in a step-wise manner, with just minor changes in the regulatory state of the mesodermal precursors occurring prior to subdivisions of the ancestral gene regulatory network subcircuits.
Results and discussion
The development of the sea cucumber, Parastichopus parvimensis
Although no sea cucumber has been developed as a model organism, several have been the subject of embryological studies (reviewed in ). Several papers reported the development of Stichopus (now referred to as Apostichopus, Parastichopus and Stichopus) species [35–37]. In particular, the early development of S. tremulus was previously described . Early development in S. tremulus is very similar to what we observe in P. parvimensis (shown in Additional file 1). Early cleavage is equal and little cell-cell adhesion is seen between blastomeres. Blastulae are formed by 16 hours after fertilization and hatch from the fertilization envelope around 26 hours. Prior to gastrulation, the embryos elongate along the animal-vegetal axis, and a thickening is observed at the vegetal pole, termed the vegetal plate. Mesenchyme cells ingress from the vegetal plate before invagination of the archenteron occurs. At the mid-gastrula stage around 48 hours of development, the mesenchyme has migrated. By 72 hours of development, the mouth has formed, and the embryo is now an early auricularia larva. Presumptive muscle cells are seen associated with the foregut, and a thickened ciliary band is evident in the oral hood and looping above the anus. A small skeletal element is situated at the posterior end of the larva, under the posterior coleom.
The ancestral echinoderm mesodermal regulatory state
Therefore, the sea cucumber, like the sea star, forms a central domain of mesodermally fated cells in which orthologs of erg, ets1, gata4/5/6, foxn2/3, tbr, tgif and gata1/2/3 are co-expressed. As such, these transcription factors, at least, represent the core set of regulatory genes that has been maintained in a single mesodermal domain of blastulae since sea cucumbers and sea stars last shared a common ancestor some 450 million years ago . This suggests that there has been strong selection for this mesodermal regulatory state in the auricularia-like larvae of these two classes of echinoderms. This pleisiomorphic mesoderm is a mixture of the regulatory states that are now segregated into separate mesodermal lineages in modern sea urchins. Of course, transcription factors uniquely expressed in the mesoderm of one organism or another have been identified and more are likely to be found.
The sea cucumber, like the sea urchin, has an SM lineage that arises from the mesoderm during early development
Exclusion of TF expression from the sea cucumber SM during gastrulation reveals additional similarities with the sea urchin skeletogenic lineage
The evolution of novelty
A comparison of the regulatory state of embryonic mesenchyme in three classes of echinoderm embryos, those of sea urchins, sea cucumbers and sea stars, provides an extraordinary opportunity to understand the trajectory of change associated with a novel morphology. The sea urchin ingresses skeletogenic, pigment and blastoceolar cells as mesenchyme. Like the sea urchin, the phylogenetically intermediate sea cucumber also ingresses a skeletogenic lineage that is distinct from blastocoelar cells. Thus, the skeletogenic mesenchyme, defined as a population of mesenchyme expressing alx1 that secretes the larval skeleton, is a synapomorphy of at least sea cucumbers and sea urchins. In the sea cucumber SM, however, this lineage secretes only a very small skeletal element in contrast to the dramatic skeleton that defines the morphology of the sea urchin pluteus larva.
An understanding of the evolutionary trajectory of echinoderm larval skeleton formation could be strengthened by a molecular study of the embryology of this process in brittle stars. If the third scenario is upheld, parsimony also supports a single origin of the echinoderm larval skeleton and, therefore, a monophyletic grouping of sea urchins, sea cucumbers and brittle stars, with the sea star as the outgroup that forms no skeleton. Molecular phylogenetic analysis finds equal support for this grouping as for there being a separate brittle star + sea star clade [11, 14, 15].
The other important, and unexpected, observation to arise from this comparison is that there is a loss of gene expression in each of the sea urchin mesodermal lineages relative to the ancestral mesoderm (Figure 5). The sea urchin uses distinct GRNs, linked in part by intercellular signaling, to drive the development of its skeletogenic, pigment and blastocoelar cells. Each of these sub-populations retains some, but not all, features of the ancestral network. Thus, partitioning of an ancestral GRN may have been an important mechanism for allowing evolutionary change.
Culturing and microinjection of embryos
Adult animals were collected off the southern California coast, USA. Spawning was induced in adult sea cucumbers by intra-coelomic injection of 100nM NGLWY-amide followed by heat shock in room temperature sea water after the method of Kato et al.. Freshly shed eggs were mixed with dilute sperm to fertilize, and embryos were cultured in artificial seawater at 15°C until the desired stage. Sea star and sea urchin embryos were prepared according to standard protocols [31, 43] and cultured at 15°C until the desired stage. Zygotes were injected as previously described for sea urchin and sea star [43, 44]. Sea cucumber microinjection was performed following protocols for sea stars.
Cloning sea cucumber and sea star transcription factor orthologs
RNA was isolated from embryos at blastula, gastrula and larva stages using the Total Mammalian RNA Miniprep kit (Sigma, St. Louis, MO, USA). First strand cDNA synthesis was performed with the iScript Select cDNA Synthesis kit (BioRad, Hercules, CA, USA). Degenerate primers were designed against conserved domains and PCR was performed under non-stringent conditions. Once partial sequences were obtained, 5′ and 3′ RACE was performed according to directions (GeneRacer; Invitrogen, Carlsbad, CA, USA). PCR primer sequences are available on request. Sequences were deposited in GenBank under the following accession numbers: PmAlx1, GenBank:JQ740823; PpAlx1, GenBank:JQ740824; PpErg, GenBank:JQ740825; PpFoxa, GenBank:JQ740826; PpFoxn2/3, GenBank:JQ740827; PpGata1/2/3, GenBank:JQ740828; PpGata4/5/6, GenBank:JQ740829; PpGcm, GenBank:JQ740830; PpTbr, GenBank:JQ740831; PpTgif, GenBank:JQ740832.
Whole mount in situ hybridization (WMISH)
Antisense digoxigenin-labeled RNA probes were generated against genes of interest. WMISH was performed as previously described for sea stars . Embryos were photographed using DIC optics on a Leica DMI4000B at 200× magnification using the Leica Application Suite software (Leica, Wetzlar, Germany).
Two-color in situ hybridization
Antisense dinitrophenol-labeled RNA probes were generated against PpAlx1, PpFoxa, PpGcm, and PmEts1. Embryos were fixed as for WMISH, and two-color WMISH was carried out as described by Yankura et al.. In some cases, embryos were counterstained with 1 μM DAPI. Embryos were photographed on a LSM 510 Meta/UV DuoScan Inverted Spectral Confocal Microscope with the ZEN 2009 Imaging suite (Carl Zeiss, Thornwood, NY, USA).
RNA overexpression of PmAlx1 and RFP mRNA in sea urchins and knockdown of PpAlx1
Full-length PmAlx1 and mCherry, an RFP derivative, were cloned into the pCS2+ vector. The resulting plasmids were linearized with HpaI and NotI, respectively. mRNA was synthesized using the mMessage mMachine kit (Ambion; Grand Island, NY, USA). The PpAlx1 morpholino antisense oligonucleotide (MASO) 5′ GGCTCACAAAGTCTGAAATAATCAT 3′ was designed by GeneTool, LLC (Philomath, O, USA). The mRNAs and MASO were injected into embryos in 200 mM KCl containing 0.5 mg/mL fluorescent rhodamine tracer.
The authors would like to thank Kristen Yankura and Alys Cheatle for helpful commentary and discussion. We also thank Hiroshi Wada for discussion of Alx1/4 expression among echinoderms and the helpful comments of three anonymous reviewers. Animals were collected by Patrick Leahy and Peter Halmay or Marinus Inc. under permit number SC-11478 from California Department of Fish and Game to VFH. This research was funded by NSF grant IOS-0844948 and the Winters Foundation.
- Carroll SB, Grenier JK, Weatherbee SD: From DNA to Diversity: Molecular Genetics and the Evolution of Animal Design. 2nd edition. Malden, Massachusetts: Blackwell Publishing Group; 2005.
- True JR, Carroll SB: Gene co-option in physiological and morphological evolution. Annu Rev Cell Dev Biol 2002, 18:53–80.PubMedView Article
- Simonnet F, Moczek AP: Conservation and diversification of gene function during mouthpart development in Onthophagus beetles. Evol Dev 2011, 13:280–289.PubMedView Article
- Fraser GJ, Hulsey CD, Bloomquist RF, Uyesugi K, Manley NR, Streelman JT: An ancient gene network is co-opted for teeth on old and new jaws. PLoS Biol 2009, 7:e1000031.View Article
- Bowsher J, Nijhout H: Partial co-option of the appendage patterning pathway in the development of abdominal appendages in the sepsid fly Themira biloba . Dev Genes Evol 2009, 219:577–587.PubMedView Article
- Janssen R, Eriksson BJ, Budd GE, Akam M, Prpic N-M: Gene expression patterns in an onychophoran reveal that regionalization predates limb segmentation in pan-arthropods. Evol Dev 2010, 12:363–372.PubMedView Article
- Saenko S, Brakefield P, Beldade P: Single locus affects embryonic segment polarity and multiple aspects of an adult evolutionary novelty. BMC Biol 2011, 8:111.View Article
- Lemons D, Fritzenwanker JH, Gerhart J, Lowe CJ, McGinnis W: Co-option of an anteroposterior head axis patterning system for proximodistal patterning of appendages in early bilaterian evolution. Dev Biol 2010, 344:358–362.PubMedView Article
- Gao F, Davidson EH: Transfer of a large gene regulatory apparatus to a new developmental address in echinoid evolution. Proc Natl Acad Sci U S A 2008, 105:6091–6096.PubMedView Article
- Hyman LH: The Invertebrates: Echinodermata. Volume IV. New York: McGraw-Hill; 1955.
- Wada H, Satoh N: Phylogenetic relationships among extant classes of echinoderms, as inferred from sequences of 18 S rDNA, coincide with relationships deduced from the fossil record. J Mol Evol 1994, 38:41–49.PubMedView Article
- Nakano H, Hibino T, Oji T, Hara Y, Amemiya S: Larval stages of a living sea lily (stalked crinoid echinoderm). Nature 2003, 421:158–160.PubMedView Article
- Nielsen MG, Popodi E, Minsuk S, Raff RA: Evolutionary convergence in Otx expression in the pentameral adult rudiment in direct-developing sea urchins. Dev Genes Evol 2003, 213:73–82.PubMed
- Pisani D, Feuda R, Peterson KJ, Smith AB: Resolving phylogenetic signal from noise when divergence is rapid: a new look at the old problem of echinoderm class relationships. Mol Phylogenet Evol 2012, 62:27–34.PubMedView Article
- Janies DA, Voight JR, Daly M: Echinoderm phylogeny including xyloplax, a progenetic asteroid. Syst Biol 2011, 60:420–438.PubMedView Article
- Wray GA: The evolution of cell lineage in echinoderms. Am Zool 1994, 34:353–363.
- Tanaka S, Dan K: Study of the lineage and cell cycle of small micromeres in embryos of the sea urchin, Hemicentrotus pulcherrimus . Dev Growth Differ 1990, 32:145–156.View Article
- Ruffins SW, Ettensohn CA: A fate map of the vegetal plate of the sea urchin ( Lytechinus variegtus ) mesenchyme blastula. Development 1996, 122:253–263.PubMed
- Oliveri P, Tu Q, Davidson EH: Global regulatory logic for specification of an embryonic cell lineage. Proc Natl Acad Sci U S A 2008, 105:5955–5962.PubMedView Article
- Croce J, Lhomond G, Lozano J-C, Gache C: ske-T, a T-box gene expressed in the skeletogenic mesenchyme lineage of the sea urchin embryo. MechDev 2001, 107:159–162.
- Ettensohn CA, Illies MR, Oliveri P, De Jong DL: Alx1, a member of the Cart1/Alx3/Alx4 subfamily of Paired-class homeodomain proteins, is an essential component of the gene network controlling skeletogenic fate specification in the sea urchin embryo. Development 2003, 130:2917–2928.PubMedView Article
- Howard-Ashby M, Materna SC, Brown CT, Chen L, Cameron RA, Davidson EH: Identification and characterization of homeobox transcription factor genes in Strongylocentrotus purpuratus, and their expression in embryonic development. Dev Biol 2006, 300:74–89.PubMedView Article
- Kurokawa D, Kitajima T, Mitsunaga-Nakatsubo K, Amemiya S, Shimada H, Akasaka K: HpEts, an ets-related transcription factor implicated in primary mesenchyme cell differentiation in the sea urchin embryo. Mech Dev 1999, 80:41–52.PubMedView Article
- Rho HK, McClay DR: The control of foxN2/3 expression in sea urchin embryos and its function in the skeletogenic gene regulatory network. Development 2011, 138:937–945.PubMedView Article
- Rizzo F, Fernandez-Serra M, Squarzoni P, Archimandritis A, Arnone MI: Identification and developmental expression of the ets gene family in the sea urchin ( Strongylocentrotus purpuratus ). Dev Biol 2006, 300:35–48.PubMedView Article
- Tu Q, Brown T, Davidson EH, Oliveri P: Sea urchin Forkhead gene family: Phylogeny and embryonic expression. Dev Biol 2006, 300:49–62.PubMedView Article
- Zhu X, Mahairas G, Illies M, Cameron RA, Davidson EH, Ettensohn CA: A large-scale analysis of mRNAs expressed by primary mesenchyme cells of the sea urchin embryo. Development 2001, 128:2615–2627.PubMed
- Ransick A, Davidson EH: cis-regulatory processing of Notch signaling input to the sea urchin glial cells missing gene during mesoderm specification. Dev Biol 2006, 297:587–602.PubMedView Article
- Lee PY, Davidson EH: Expression of Spgatae, the Strongylocentrotus purpuratus ortholog of vertebrate GATA4/5/6 factors. Gene Expr Patterns 2004, 5:161–165.PubMedView Article
- Duboc V, Lapraz F, Saudemont A, Bessodes N, Mekpoh F, Haillot E, Quirin M, Lepage T: Nodal and BMP2/4 pattern the mesoderm and endoderm during development of the sea urchin embryo. Development 2010, 137:223–235.PubMedView Article
- Cameron RA, Fraser SE, Britten RJ, Davidson EH: Macromere cell fates during sea urchin development. Development 1991, 113:1085–1091.PubMed
- Hinman VF, Davidson EH: Evolutionary plasticity of developmental gene regulatory network architecture. Proc Natl Acad Sci U S A 2007, 104:19404–19409.PubMedView Article
- Hinman VF, Nguyen AT, Cameron RA, Davidson EH: Developmental gene regulatory network architecture across 500 million years of echinoderm evolution. Proc Natl Acad Sci U S A 2003, 100:13356–13361.PubMedView Article
- McCauley BS, Weideman EP, Hinman VF: A conserved gene regulatory network subcircuit drives different developmental fates in the vegetal pole of highly divergent echinoderm embryos. Dev Biol 2010, 340:200–208.PubMedView Article
- Holland ND: Electron microscopy of the sea cucumber, Stichopus tremulous, from unfertilized egg through hatched blastula. Acta Zool 1981, 62:89–112.View Article
- Smiley S: Metamorphosis of Stichopus califronicus (Echinodermata: Holothuroidea) and its phylogenetic implications. Biol Bull 1986, 171:611–631.View Article
- Shoguchi E, Harada Y, Numakunai T, Satoh N: Expression of the Otx gene in the ciliary bands during sea cucumber embryogenesis. Genesis 2000, 27:58.PubMedView Article
- Oliveri P, Walton KD, Davidson EH, McClay DR: Repression of mesodermal fate by foxa, a key endoderm regulator of the sea urchin embryo. Development 2006, 133:4173–4181.PubMedView Article
- Koga H, Matsubara M, Fujitani H, Miyamoto N, Komatsu M, Kiyomoto M, Akasaka K, Wada H: Functional evolution of Ets in echinoderms with focus on the evolution of echinoderm larval skeletons. Dev Genes Evol 2010, 220:107–115.PubMedView Article
- Hinman VF, Davidson EH: Expression of a gene encoding a Gata transcription factor during embryogenesis of the starfish Asterina miniata . Gene Expr Patterns 2003, 3:419–422.PubMedView Article
- Ettensohn CA, Kitazawa C, Cheers MS, Leonard JD, Sharma T: Gene regulatory networks and developmental plasticity in the early sea urchin embryo: alternative deployment of the skeletogenic gene regulatory network. Development 2007, 134:3077–3087.PubMedView Article
- Kato S, Tsurumaru S, Taga M, Yamane T, Shibata Y, Kaoru O, Fujiwara A, Yamano K, Yoshikuni M: Neuronal peptides induce oocyte maturation and spawning of a sea cucumber, Apostichopus japonicus . Dev Biol 2009, 326:169–176.PubMedView Article
- Hinman VF, Nguyen AT, Davidson EH: Expression and function of a starfish Otx ortholog, Am Otx: a conserved role for Otx proteins in endoderm development that predates divergence of the eleutherozoa. Mech Dev 2003, 120:1165–1176.PubMedView Article
- Cheers MS, Ettensohn CA: Rapid microinjection of fertilized eggs. Methods Cell Biol 2004, 74:287–310.PubMedView Article
- Yankura KA, Martik ML, Jennings CK, Hinman VF: Uncoupling of complex regulatory patterning during evolution of larval development in echinoderms. BMC Biol 2010, 8:143.PubMedView Article
- Altschul SF, Wootton JC, Gertz EM, Agarwala R, Morgulis A, Schaffer AA, Yu Y-K: Protein database searches using compositionally-adjusted substitution matrices. FEBS Journal 2005,272(20):5101–5109.PubMedView Article
- Anisimova M, Gascuel O: Approximate likelihood ratio test for branches: a fast, accurate and powerful alternative. Systematic Biology 2006,55(4):539–552.PubMedView Article
- Dereeper A, Guignon V, Blanc G, Audic S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M, et al.: Phylogeny.fr: Robust phylogenetic analysis for the non-specialist. Nucleic Acids Research 2008,36(suppl 2):W465-W469.PubMedView Article
- Guindon SP, Gascuel O: A simple, fast, and accurate algorithm to estimate large Phylogenies by maximum likelihood. Systematic Biology 2003,52(5):696–704.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.