The Hawaiian Rhodophyta Biodiversity Survey (2006-2010): a summary of principal findings
© Sherwood et al; licensee BioMed Central Ltd. 2010
Received: 27 July 2010
Accepted: 22 November 2010
Published: 22 November 2010
The Hawaiian red algal flora is diverse, isolated, and well studied from a morphological and anatomical perspective, making it an excellent candidate for assessment using a combination of traditional taxonomic and molecular approaches. Acquiring and making these biodiversity data freely available in a timely manner ensures that other researchers can incorporate these baseline findings into phylogeographic studies of Hawaiian red algae or red algae found in other locations.
A total of 1,946 accessions are represented in the collections from 305 different geographical locations in the Hawaiian archipelago. These accessions represent 24 orders, 49 families, 152 genera and 252 species/subspecific taxa of red algae. One order of red algae (the Rhodachlyales) was recognized in Hawaii for the first time and 196 new island distributional records were determined from the survey collections. One family and four genera are reported for the first time from Hawaii, and multiple species descriptions are in progress for newly discovered taxa. A total of 2,418 sequences were generated for Hawaiian red algae in the course of this study - 915 for the nuclear LSU marker, 864 for the plastidial UPA marker, and 639 for the mitochondrial COI marker. These baseline molecular data are presented as neighbor-joining trees to illustrate degrees of divergence within and among taxa. The LSU marker was typically most conserved, followed by UPA and COI. Phylogenetic analysis of a set of concatenated LSU, UPA and COI sequences recovered a tree that broadly resembled the current understanding of florideophyte red algal relationships, but bootstrap support was largely absent above the ordinal level. Phylogeographic trends are reported here for some common taxa within the Hawaiian Islands and include examples of those with, as well as without, intraspecific variation.
The UPA and COI markers were determined to be the most useful of the three and are recommended for inclusion in future algal biodiversity surveys. Molecular data for the survey provide the most extensive assessment of Hawaiian red algal diversity and, in combination with the morphological/anatomical and distributional data collected as part of the project, provide a solid baseline data set for future studies of the flora. The data are freely available via the Hawaiian Algal Database (HADB), which was designed and constructed to accommodate the results of the project. We present the first DNA sequence reference collection for a tropical Pacific seaweed flora, whose value extends beyond Hawaii since many Hawaiian taxa are shared with other tropical areas.
The Rhodophyta, or red algae of Hawaii represent one of the best studied red algal floras worldwide from a morphological perspective , making it an ideal flora for molecular characterization. Endemism estimates for the marine (19.5%) and freshwater (20.0%) algal floras of the islands are based on species lists from compiled traditional taxonomic surveys, and the potential for unrecognized genetic diversity is high [1, 2]. It is almost certain that these estimates will change as molecular techniques are applied on a more routine basis to the Hawaiian flora and in other areas of the tropical Pacific. Until this survey, only a small proportion of the published studies of Hawaiian Rhodophyta employed molecular approaches, although those few resulted in substantial taxonomic revision to the Hawaiian representatives of the genera Akalaphycus, Galaxaura, Ganonema and Liagora [3–5]. Additionally, a new red algal order (the Pihiellales) was described from Hawaiian material as a result of molecular phylogenetic analyses .
There has also been a bias toward the study of large, conspicuous members of the Hawaiian seaweed flora . The red algae, however, display a vast array of morphologies and reproductive strategies , and the breadth of this diversity is often overlooked in general survey efforts. It is certain that more concentrated study of the less conspicuous turf and crustose morphologies, as well as the epiphytic and parasitic representatives, will yield unrecognized diversity [1, 8]. Indeed, examination of a number of red algal taxa using molecular techniques has revealed cryptic diversity that was previously unrecognized based on morphology and anatomy, e.g. [9, 10].
Red algal systematics and taxonomy are fraught with a number of issues, often making identification difficult. Since the end of the 19th century, the classification of red algae has been largely based on female reproductive anatomy and post-fertilization events [7, 11, 12]. Later additions to the suite of taxonomic characters included the ultrastructure of pit connections between cells  and the use of DNA sequence data [14, 15]. Current red algal systematics usually employs a blend of biological, morphological and evolutionary species concepts, depending on the subgroup under study. For example, although the ability to interbreed and produce fertile offspring is not usually explicitly tested in phycological studies, some researchers have shown the biological species concept to be applicable [16, 17]. Morphology continues to play a key role in red algal systematics, although these data are now typically viewed in an evolutionary context, and the inclusion of molecular data has become commonplace in the last two decades.
The traditional system of relying exclusively on reproductive characters for taxonomic purposes has several disadvantages. These characters can be evaluated only for those taxa that undergo sexual reproduction, which may be encountered rarely in some taxa. Furthermore, a number of red algae from Hawaii are believed to be strictly asexual, including the genera Stylonema, Chroodactylon, Erythrocladia, Erythrotrichia, Compsopogon and Hildenbrandia [1, 18, 19]. For others, only asexual stages have been found in Hawaii, even though they are known to be sexually reproductive in other geographical regions (e.g. Peyssonnelia ). Some taxa are also highly seasonal in their timing of reproduction, meaning that collections must be made in a narrow window of the year to observe the necessary characters (e.g. Acrosymphyton, Dudresnaya and Naccaria ). Still other taxa rarely undergo sexual reproduction, with no apparent seasonal pattern (e.g. Acrochaetium ).
Heteromorphic taxa (those with morphologically distinct gametophyte and sporophyte life history phases) present another challenge for red algal systematics and, in Hawaii, include members of the Bangiales, Batrachospermales, Bonnemaisoniales, Gigartinales and some members of the Nemaliales [1, 2, 20]. Until recently, researchers have relied on culturing methods to link heteromorphic counterparts of the same taxon, and this approach has had profound impacts on the taxonomy of red algae [21, 22]. A number of taxa in the Hawaiian Islands are known for which these two phases are spatially separated and morphologically distinguishable, and collecting often results in the recovery of only one of these phases. Often the gametophyte is the more conspicuous of the two generations, as an upright phase, while the tetrasporophyte is an overlooked filamentous or crustose stage [1, 23]. In cases like these, DNA sequences associated with each morphological specimen can reveal molecular identity of morphologically distinct life history stages.
Biodiversity research is based on the ability to apply taxonomic names to collections. This is often not simple for the red algae, for reasons described above. DNA barcoding is a more recent development in biodiversity studies that has the potential to allow greater discrimination among cryptic species than traditional taxonomic methods. DNA barcodes are short, species-specific DNA sequences from one or more genetic loci, allowing for differentiation of individuals at the species level . The purpose of a DNA barcode is to assign unknown individuals to species and to enhance the discovery of new species [25, 26], both of which are fundamental aims of this project. One of the major components of the present study is the acquisition and analysis of short DNA marker sequences from the nuclear (partial 28 S rRNA gene, or LSU ), plastidial (partial 23 S rRNA gene, Universal Plastid Amplicon, or UPA [28, 29]) and mitochondrial (partial cytochrome c oxidase I gene, or COI ) genomes of Hawaiian red algae to assess intraspecific variation and to provide a DNA sequence framework for the Hawaiian Rhodophyta.
The Hawaiian Rhodophyta Biodiversity Survey (2006-2010) was undertaken to document and archive the marine, freshwater and terrestrial red algae of the Hawaiian Islands with an unprecedented depth of analysis and coverage by vouchering morphological samples and total genomic DNA extracts and undertaking an assessment of DNA sequence diversity for as many samples as possible. Compilations of freely available taxonomic data sets such as this fulfil the mandate of the "Biodiversity Commons" . Resulting archived collections are invaluable for studies into the phylogenetic relationships among taxa, species-level descriptions of the flora, preliminary population genetic analyses, and biogeographic studies within the Hawaiian Islands and among other tropical Pacific island systems.
The Hawaiian archipelago extends in a northwesterly direction for almost 2,400 km, includes a series of small islands, atolls and seamounts , and is well known as a center of endemism and biological uniqueness . In 2006 the Northwestern Hawaiian Islands were declared the Papahanaumokuakea National Marine Monument, and this large portion of the archipelago is understudied for most algae. An imminent need exists to document the biodiversity of this unique island chain in the face of threats from ongoing alien and/or invasive species, algal blooms, climate change and corresponding sea-level rise [33, 34]. Characterizing native biodiversity provides essential baseline data for future detection of algal introductions, and for monitoring changes in assemblages at particular sites. To this end, we undertook a series of collecting expeditions and analyses to document the Hawaiian red algal diversity, and here present a summary of our findings.
Diversity of sampling locations
Summary of Rhodophyta collections
The vast majority of accessions analyzed were from marine and brackish locations (1,928), while 16 were from freshwater habitats and two were from terrestrial locations. A total of 24 orders, 49 families, 152 genera and 252 species/subspecific taxa of red algae are represented in the collections (species counts do not include those identified only to genus as "sp."). Additional file 2 presents a checklist of taxa collected, identified and analyzed as part of the biodiversity survey, including new island and state records.
New records reported here do not include most of the coralline red algae (orders Corallinales and Sporolithales), which are being studied in detail by other researchers. In our new records, we included specimens from the Bernice P. Bishop Museum (BISH) that were not reported in Abbott's 1999  assessment of distributional records, and also compared our distributional data to publications since this treatise and to those not covered therein. Included are those for freshwater red algae [35–39], descriptions of more recently recognized species of marine algae (e.g. Dasya atropurpurea ; the order Pihiellales ) and reports of previously unrecognized taxa (e.g. marine species of the genus Hildenbrandia ).
One new order of red algae (the Rhodachlyales) was recorded from our survey collections, although a number of additional new ordinal names are applied compared to those listed in Abbott  due to nomenclatural revision and inclusion of previously unreported taxa (i.e. the orders Acrosymphytales, Batrachospermales, Colaconematales, Compsopogonales, Halymeniales, Hildenbrandiales, Nemastomatales, Peyssonneliales, Pihiellales, Plocamiales, Rhodachlyales, Sporolithales, Stylonematales and Thoreales). The family Rhodachlyaceae is reported for the first time, as are the genera Chroothece, Madagascaria, Pseudobangia and Rhodachlya. Species designations for most of these newly reported genera are still being assessed, but it is likely that many of them represent undescribed species. Many new island distributional records (196 in total) were determined from the survey collections (Additional file 2).
Abbott  listed 343 species of marine red algae in the Hawaiian Islands, of which 218 were also identified from our biodiversity survey collections. The remaining 125 taxa, which included large numbers of species in the genera Colaconema, Ceramium/Gayliella, Herposiphonia and Polysiphonia/Neosiphonia, were not identified during our survey; however, these are primarily diminutive taxa that require a great deal of specialized morphological expertise for identification, and were not surveyed as readily as the more conspicuous macrophyte taxa. It is also possible that some of these taxa are not permanent members of the Hawaiian red algal flora since attempts to collect from reported localities have not yielded examples of these taxa.
Summary of DNA sequence data
DNA sequence diversity of Hawaiian Rhodophyta
Some trends of marker utility are evident. Overall, a pattern of #LSU sequences > #UPA sequences > #COI sequences was observed in a number of instances, including for the Acrosymphytales, Gigartinales, Nemastomatales, Peyssonneliales, Plocamiales and Pihiellales (Figures 16, 17 and 18, Additional file 3), the Gelidiales, Gracilariales and Halymeniales (Figures 19, 20 and 21, Additional file 3), and many of the Ceramiales (Figures 28, 29, 30, 31, 32 and 33, Additional file 3). In contrast to this, the large number of COI sequences for the Bonnemaisoniales (Figure 12) relative to LSU (Figure 10) and UPA (Figure 11) is due to targeting of the COI region for a mtDNA-based phylogeographic study . Overall, amplification and sequencing success was much higher for the LSU and UPA markers than for the COI marker, based on the laboratory protocols employed in this study.
A number of taxa were represented by multiple clusters of sequences according to one or more of the markers, including Hydrolithon reinboldii (Corallinales) (Figures 13, 14 and 15), Grateloupia filicina (Halymeniales) (Figures 19, 20 and 21), Dichotomaria marginata (Nemaliales) (Figures 22 and 24), Champia parvula (Rhodymeniales) (Figure 27), Spyridia filamentosa (Ceramiales) (Figures 28, 29 and 30), Martensia flabelliformis and M. fragilis (Ceramiales) (Figures 31, 32 and 33), Laurencia crustiformans and L. mcdermidiae (Ceramiales) (Figures 35 and 36), and Amansia glomerata (Ceramiales) (Figures 37, 38 and 39). These clusters may be due to cryptic or incipient speciation, or phylogeographic structure, and these taxa should be investigated through detailed anatomical and morphological comparison of specimens (and comparison to the type specimens) as well as increased sampling for molecular analyses, in order to clarify which of these processes accounts for the diversity observed in the NJ trees.
Other taxa were identified that will require taxonomic attention to reconcile patterns of molecular diversity with currently applied taxonomic names. For example, the genera Ahnfeltiopsis (Gigartinales) (Figures 16 and 18), Peyssonnelia (Peyssonneliales) (Figures 16 and 17) and Gelidium (Gelidiales) (Figures 19 and 20) all have clusters of sequences for at least one marker with multiple species names applied. Along these same lines, some pairs of taxa were shown to possibly be conspecific according to all three molecular markers employed here (e.g. Gracilaria eppihippisora /G. salicornia, and G. coronopifolia /G. dotyi [Figures 19 and 20]).
Phylogenetic analysis of concatenated UPA, LSU and COI sequences
Phylogeographic analyses of widespread taxa
A. glomerata is widely distributed throughout the Hawaiian archipelago . Statistical parsimony analysis yielded three distinct networks; one containing NWHI samples (French Frigate Shoals), and two containing MHI samples (Figure 41a; ).
This species is sought after as a food source in the Hawaiian culture . Statistical parsimony analysis of COI sequences for this taxon yielded two networks formed at the 95% connection limit; one that contains three samples from the south shore of the island of Oahu (which possibly represent a recent introduction), belonging to Lineage 4 (sensu [43, 45]), and the remainder forming a large network that is broadly divided into Lineage 1 (widespread throughout the MHI and Nihoa) and Lineage 2 (NWHI and MHI in distribution) (Figure 41b; ).
A single network was formed by statistical parsimony analysis containing both MHI and NWHI samples (Figure 41c). The NWHI samples were closely related to one another and separated from MHI samples. Additionally, greater nucleotide diversity was present for MHI samples, especially those from the island of Oahu.
Two networks were formed by statistical parsimony analysis, both of which contained only MHI samples (Figure 41d). Samples from the island of Molokai formed two distinct haplotypes in the most populated network, while samples from the other islands (Kauai, Oahu, Maui and Hawaii) were represented in the most geographically diverse haplotypes (Figure 41d).
Hawaiian Rhodophyta with low intraspecific divergence
Coelothrix irregularis (Rhodymeniales)
Oahu, Molokai, Maui, Hawaii
Lithophyllum insipidum (Corallinales)
Oahu, Lanai, Hawaii
Haliptilon subulatum (Corallinales)
Kauai, Oahu, Molokai, Hawaii
Tolypiocladia glomerulata (Ceramiales)
Oahu, Molokai, Maui, Hawaii
Polysiphonia howei (Ceramiales)
Oahu, Maui, Hawaii
Acanthophora pacifica (Ceramiales)
Oahu, Molokai, Maui, Kahoolawe, Hawaii
Laurencia majuscula (Ceramiales)
Laysan, French Frigate Shoals, Kauai, Oahu, Molokai, Lanai
The Hawaiian Rhodophyta Biodiversity Survey yielded only a few new taxonomic records for the Hawaiian Islands, but many new distributional records. The flora has been remarkably well-studied from a morphological and anatomical perspective, most notably by Abbott , and many of her designations are supported by our molecular investigations. However, many taxa have also been revealed to be in need of in-depth taxonomic study based on our molecular survey data due to instances of suspected cryptic or incipient speciation (e.g. Amansia glomerata and Spyridia filamentosa; see results for full list), or conspecificity (e.g. Gracilaria eppihippisora and G. salicornia). As such, our survey data provide a unique opportunity to flag Hawaiian red algal taxa for further study. Although identification of these taxonomic issues, rather than full resolution, was the goal of the Hawaiian Rhodophyta Biodiversity Survey, several of these flagged taxa have been investigated in detail under the auspices of the project (e.g. [8, 10, 43]), and we hope that other researchers will be able to make use of the data for similar purposes, for it is detailed taxonomic studies that will ultimately lead to a floristic revision of the Hawaiian red algae. Also, additional taxonomic expertise for some of the lineages of red algae investigated during the project (in particular, for many of the Ceramiales), will be needed for taxonomic resolution. However, despite adding relatively few new taxa to the species list, we have substantially updated the reported range of many taxa, with 196 new island-level distributional records included from our sampling records.
The vast majority of accessions analyzed were from marine and brackish locations (1,928, or 99%), while 16 were from freshwater habitats and two were from terrestrial locations. This breakdown is close to the accepted consensus that, worldwide, approximately 97% of red algal species are marine .
The three molecular markers chosen for diversity assessment of our red algal collections (nuclear LSU, mitochondrial COI and plastid UPA) differ markedly in both their degree of conservation and ease of amplification/sequencing, based on the protocols that were applied to our survey samples. Since the Hawaiian Rhodophyta Biodiversity Survey encompassed an evolutionarily diverse set of samples, and we aimed to generate a large number of sequences (i.e. >2,000), we were constrained by the need for minimal variations of protocol for each marker (including primer design, PCR cocktail components/concentrations and amplification conditions). The LSU and UPA markers were sequenced without modification (i.e. according to the primers and conditions described in  and  respectively), and these first two markers yielded the highest and second highest overall number of sequences, respectively. In contrast, the COI marker required the design of an additional reverse primer, and yielded the lowest number of sequences. Ease of data acquisition and performance characteristics of each of the three markers are further examined and discussed in , as well as below.
All clean sequences (those with clearly readable chromatograms that could be assembled with confidence) were retained in HADB, regardless of whether or not they corresponded to the taxon under study; contaminant sequences (e.g. from epiphytes) were re-labelled as such and given new accession numbers. All three markers yielded some contaminant sequences, but the UPA marker yielded the highest number (30), followed by LSU (15) and COI (7).
The LSU marker was the simplest for which to reliably obtain clean sequence data; only a single pair of primers was needed, and most samples were successfully sequenced on the first attempt. The disadvantage of the LSU marker, however, was a lack of reliable species-level resolution (e.g. Figure 21, Nemaliales LSU sequence data show no divergence between some members of the genera Akalaphycus, Ganonema, Izziella, Liagora and Trichogloea). A recent investigation of the characteristics and phylogenetic "performance" of the LSU, UPA and COI markers based on this same data set confirmed that the LSU marker is the most conserved of the three, that it was the easiest for which to obtain sequence data, and illustrated that florideophyte red algal sequence data for the marker became saturated at F84 distances of approximately 0.30 (although there was variation when individual orders were examined) . The LSU marker was included in the survey to provide representation from the nuclear genome (both the UPA and COI are organellar markers). However, given that the level of resolution for the LSU marker differs so strongly from the other two markers (Figure 4), it would be difficult to recommend its inclusion in future biodiversity surveys.
The UPA marker, which was first described in 2006  and tested for a selection of algal lineages in 2007 , can be amplified and sequenced for most cyanobacterial and plastid-containing algal lineages, with only a few genus-level exceptions noted to date (e.g. Cladophora, Symbiodinium). This universality is exceptional among algal primers, and lends great potential for the use of this marker in assessment of environmental samples that contain multiple algal lineages . For this reason, we believe it is important to continue to generate framework data for this marker for as many algal representatives as possible. Recent work has shown that the UPA marker falls between in the COI and LSU markers in terms of ease of data acquisition, that it is quite conserved (although slightly less so than LSU), and that saturation is not an issue for any group of the florideophyte red algae based on UPA sequences .
The COI marker was somewhat more difficult to amplify and sequence across red algal lineages, and two combinations of primers were used to obtain the 639 sequences reported here (although new primers have been designed recently that are generally more successful; G.W. Saunders, pers. comm.). Others have reported requiring up to seven primer combinations to obtain COI sequences for samples within a single family of red algae . However, the COI marker is widely recognized as the DNA barcode for red algae [9, 41, 49], and as such, provides a useful sequence data set due to the immediate comparative power of COI sequences against the 11,925 red algal DNA barcode sequences in the Barcode of Life Data Systems (BOLD) v.2.5 database (accessed 05 November 2010 ). This trade-off between universal amplification and resolution has been noted in the past [42, 49]. We recently reported that these COI florideophyte red algal sequence data were the most difficult of the three markers to obtain (i.e. low ease of acquisition), but that the COI marker is the least conserved (and is especially variable at the third codon position), and that saturation of the entire data set occurs at F84 distances of approximately 0.11 (although this number varies substantially depending on which order is being examined) .
Maximum likelihood analysis of the concatenated LSU+UPA+COI sequences yielded a phylogenetic tree that included support for most of the florideophyte red algal orders and families included in the data set (i.e. 3/17 orders and 4/37 families were not recovered as monophyletic; Figure 40). A number of explanations may be responsible for the lack of full resolution in the phylogeny, including insufficient sequence data [47, 51], insufficient taxonomic representation, or unresolved identifications of some specimens. Verbruggen et al. , in their recent research on the data requirements for resolving the red algal tree of life, indicated that up to 2.8 × 105 nucleotides of sequence data may be required to resolve some regions of the phylogeny, which is orders of magnitude more data than that provided by our concatenated three-marker data set. We  preliminarily evaluated the phylogenetic performance of the three markers with the same data set by assessing changes in neighbour-joining bootstrap support with the addition of marker data (progressing from UPA, to UPA+LSU, to UPA+LSU+COI), and also concluded that the fully concatenated phylogeny was supported largely at the family and ordinal level, with higher level relationships not resolved. Nonetheless, the sequences of the three markers will make a valuable addition to the red algal species assignment framework of data, and in concatenated form are useful for examining phylogenetic relationships at lower taxonomic levels.
Considering the taxa studied using statistical parsimony analysis (Figure 41), the debate arises as to whether observed breaks in networks represent strong phylogeographic signal or species boundaries. It has been suggested that the 95% connection limit under statistical parsimony, which is commonly used in phylogeographic studies, is a useful criterion for indicating biological species boundaries . Applying this criterion to our analyses suggests that cryptic species may be common in the Hawaiian red algal flora. Only one of the four taxa investigated formed a single haplotype network at the 95% connection limit (Dasya iridescens), while the remaining taxa formed two (Asparagopsis taxiformis and Dichotomaria marginata) or three networks (Amansia glomerata). If this pattern holds true for a number of other Hawaiian red algal taxa (which will only be determined as more taxa are studied in depth), then there may be many more species than currently recognized in the flora. Sorting out phylogeographic signal from species boundaries requires an understanding of the limits of phenotypic plasticity for each case. Most of our identifications were made, at least preliminarily, based on morphology and anatomy. We largely relied on the detailed descriptions of Abbott  but supplemented these with more recent or taxon-specific studies that were relevant to each taxon (a list of this additional literature can be found in ). Nonetheless, lately there have been a number of cases in the phycological literature where cryptic species complexes have been revealed through molecular data, and in a number of these instances other characters (morphological, anatomical, life history, physiological, etc.) have later been identified that allowed recognition of the species on grounds other than molecular data, e.g. [49, 53]. More detailed phylogeographic, morphological/anatomical and life history analyses of suspected species complexes will be necessary for confirmation, but the DNA sequences from the present survey provide a valuable starting point for identification of these phenomena, and pave the way for taxonomic revision of the floristic checklist.
Interestingly, not all common Hawaiian Rhodophyta species have notable intraspecific sequence variation. Eight taxa representing three different orders (Table 1) were reported to have little to no sequence divergence for the three markers sequenced. Although these taxa require stronger sampling for confirmation and resolution of these trends, one important finding from the survey is that phylogeographic patterns are likely to vary, even for taxa considered to be native to the Hawaiian Islands, and that histories of individual taxa need to be studied (e.g. time since colonization or speciation, colonization order of the islands, interisland dispersal vectors, introduction events from other regions) to obtain a fuller understanding of intraspecific molecular patterns. It does not seem to be the case that increased within-taxon sampling will always reveal obvious phylogeographic patterns. This difference could be due to many factors, including the length of time since colonization of the species (or speciation, in the case of endemic taxa), ease of dispersal, dispersal vectors and type of life history (i.e. whether there is an easily spread life history stage, such as a floating, filamentous form), habitat and physicochemical parameter tolerances.
Ekrem et al.  emphasized that DNA barcode data will only ever be as useful as the reference library of sequences against which comparisons are made. This is a compelling argument for support of biodiversity surveys of entire floras, such as this one, that attempt to account for the diversity of a region with both morphological and molecular sequence characterizations, and that make these data available to others in easily-accessible formats for comparisons. Zhang et al. , however, rightly point out that sample sizes for most DNA barcode studies are insufficient to assess genetic diversity of a species, and the level of intraspecific sampling for our project is no exception. Given the limited time and funds for the survey, we aimed to strike a balance between representation of as many taxa as possible and sufficient coverage of common taxa to assess intraspecific diversity. We emphasize that the patterns discussed here are largely preliminary, and that a more thorough examination is warranted to reveal the full story.
The value of the Hawaiian Rhodophyta Biodiversity Survey is multifaceted: through this project we have generated baseline molecular diversity data for much of the flora, revealed taxa requiring further systematic investigation (both in terms of possible cryptic species that have not been previously recognized, and morphologically/molecularly similar taxa that may be synonymous), begun to elucidate phylogeographic patterns for Hawaiian red algae, revised the species checklist for Hawaiian Rhodophyta, and reported a number of new records (both taxonomic and distributional). All of these data (our own collections plus those from BISH), as well as photographs and micrographs, specimen details and collecting site information, are available via the publicly accessible Hawaiian Algal Database , and the data have also been shared with the Global Biodiversity Information Facility  and GenBank . Voucher specimens and the DNA extract library are available to other researchers, enhancing the long-term value of the survey collections. The value of the survey will also increase as other studies make use of the data and our initial observations are developed into additional projects.
Here we describe the principal findings of the Hawaiian Rhodophyta Biodiversity Survey, a four-year effort to inventory and assess the species of red algae from the Hawaiian Islands. Numerous new taxonomic and distributional records resulted from the collections, which allowed for an updated taxonomic checklist of Hawaiian red algae. Three molecular markers (one each from the nucleus, plastid and mitochondrion) were sequenced for numerous collections, and the data from each marker compared to establish baseline molecular diversity profiles for the flora. A number of cases were highlighted where the molecular data helped to identify taxa or groups of taxa in need of systematic revision, illustrating the value of constructing this baseline molecular diversity data set. In the case of widespread and common species, the sequence data were also helpful in beginning to elucidate phylogeographic patterns, which differed substantially depending on the taxon. A maximum likelihood analysis based on concatenated sequence data for 17 red algal orders yielded a tree topology that is largely consistent with current understanding of red algal phylogeny, but with lower support at deeper nodes, indicating that the DNA barcode-like markers have phylogenetic utility at and below the ordinal level. The value of each of the markers for biodiversity surveys was also assessed, and recommendations made for inclusion or exclusion in future surveys of algal floras.
Samples were field-collected or subsampled from accessioned material at the Bernice P. Bishop Museum's Herbarium Pacificum (BISH). Annual collecting expeditions to the islands of Kauai, Oahu, Molokai, Maui, Lanai and Hawaii were undertaken to survey nearshore coastal areas and streams. Samples from the Main Hawaiian Islands of Niihau and Kahoolawe were processed as they were opportunistically obtained, as were those from the Northwestern Hawaiian Islands. We aimed to collect three specimens per species, although there was variation around this number due to rare or, conversely, extremely widespread taxa.
Sample processing and analysis
Field-collected samples were processed as described in Sherwood . Samples were examined using light microscopy (Olympus BX-41 compound microscope [Olympus, Center Valley, PA, USA] and a Zeiss SteREO Discovery.v12 stereomicroscope, Carl Zeiss Microimaging Inc., Göttingen, Germany). Identifications were made using the most current sources, with  forming the foundation for identifications and more recent primary literature being added where appropriate [3–6, 40]. Material analyzed from BISH was labelled as per the assigned identification on the herbarium sheet, except in cases of obvious mis-identification or where microscopic examination revealed the sample to belong to a different taxon (herbarium sheets were subsequently annotated to reflect these changes).
DNA was extracted from field-collected and herbarium material following the methods described in . All DNA extracts are housed in the Hawaiian Algal DNA Library at the University of Hawaii at Manoa, which is maintained at -80°C in the Sherwood Laboratory. Duplicate DNA extracts of material from herbarium sheets housed at BISH are maintained at Bishop Museum's Pacific Center for Molecular Biodiversity. Three short DNA markers were amplified and sequenced for as many of the collections as possible; these markers were selected to represent different genomes within the cell (i.e. nucleus - part of the 28 S rRNA gene (LSU), plastid - the Universal Plastid Amplicon (UPA) [28, 29], and mitochondrion - the cytochrome oxidase I barcode region (COI) ). PCR primers and cycling conditions used for each of the markers are referenced in , with information on additional LSU and COI primer derivation as follows.
We aligned the entire GenBank content for Rhodophyta LSU (>650 sequences ranging from ~800-2700 bp) and designed primers to amplify a central portion where most of the data overlapped. These primers, nu28SF (5'-GGAATCCGCYAAGGAGTGTG-3'; Tm = 56.9°C) and nu28SR (5'-TGCCGACTTCCCTTACCTGC-3'; Tm = 59.7°C) produced an amplicon of variable length; 580-650 bp depending on the taxon. Cycling conditions for this amplicon were as follows: 94°C for 2 min, 40 cycles of 94°C for 20 s, 55°C for 30 s, 72°C for 50 s, and a 5 min final extension at 72°C.
Data storage and dissemination
All relevant project data (photographs and micrographs, locality information, DNA sequences, identifications, voucher type and archival information) were entered into the Hawaiian Algal Database (HADB); a database that was specially designed and constructed for the project to allow organization and accessibility of the data. HADB is internet accessible at http://algae.manoa.hawaii.edu and is fully described in .
Mapping of GPS data
Collection locality data (GPS points) were downloaded from HADB and plotted on a map of the Hawaiian Islands using GPS Visualizer . All coastal, freshwater and terrestrial MHI points were included, while those from the NWHI were excluded since most were only at island-scale resolution.
Construction of DNA sequence diversity trees
Sequences for each marker were downloaded from HADB and aligned using Clustal X . Neighbor-joining trees of sequences for the UPA, COI and LSU markers were constructed based on the maximum-likelihood composite model in MEGA v.4 . Separate trees were constructed for each marker, and for each red algal order (or groups of related orders when sample numbers were low enough to still visualize taxon labels on trees), or by family (in the case of the order Ceramiales, which accounts for approximately half of the Hawaiian red algal diversity). Summary labels were added to the trees to indicate the following: a) previously identified lineages of species from the literature (in these cases the citations are given in the figure legends), b) groups of sequences that, to the best of our knowledge, belong to the same species although some sequence divergence is present (e.g. Chroodactylon "ornatum"; Figure 5), c) groups of sequences that do not have a consistent taxonomy but nonetheless likely represent a single taxon (e.g. Platoma sp.; Figure 16), and d) groups of sequences that do or do not have a consistent taxonomy and likely represent multiple species (e.g. Peyssonnelia spp.; Figures 16 and 17).
Phylogenetic analysis of concatenated sequence data
All accessions for which sequences of all three markers were available (n = 284) were concatenated using the software package APE  with R . Fast evolving loops of LSU (indel rich regions) were removed because sequence homology could not be ascertained with confidence, leaving only stems for this marker. Similarly, a very short segment of UPA, which aligned ambiguously on the scale of the entire data set, had to be removed. Concatenated sequences were analyzed with the rapid bootstrap (RBS) and maximum likelihood (ML) search algorithm of the RAxML Blackbox server http://phylobench.vital-it.ch/raxml-bb/. RAxML uses a hill-climbing algorithm that searches the tree space using Lazy Subtree Rearrangement  while allowing specification of partitioned models useful for concatenated data sets. An exploratory analysis of the sequences can be found in  which demonstrates the varying levels of substitution saturation and substitution rates for UPA, LSU, and all three codon positions of the COI marker. Following this five-partition scheme, the RAxML analysis was launched with a GTR+G4 model and branch length optimization per partition with the default 100 RBS. The returned ML tree was summarized by family and order in MEGA  and RBS values ≥70% were illustrated.
Phylogeographic patterns for Amansia glomerata, Asparagopsis taxiformis, Dichotomaria marginata and Dasya iridescens (representatives of widespread Hawaiian red algae with evident structuring at the molecular level) were estimated using COI marker sequences (which is commonly used for phylogeographic inference, e.g.  and the statistical parsimony method of Templeton et al.  using TCS v.1.21 , employing a 95% connection limit and treating gaps as a 5th state. Sequence divergences for the widespread taxa Acanthophora pacifica, Champia parvula, Coelothrix irregularis, Haliptilon subulatum, Laurencia majuscula, Lithophyllum insipidum, Polysiphonia howei and Tolypiocladia glomerulata were calculated based on p-distances using MEGA v.4  and reported as percent identities.
This research was funded by a grant from the National Science Foundation (DEB-0542608 to A.R.S. and G.G.P.) to support the Hawaiian Rhodophyta Biodiversity Project. T.S. was partially supported by The David and Lucile Packard Foundation (grant #2006-30569 to I.A. Abbott), and this grant also supported the collection and processing of many samples from the Northwestern Hawaiian Islands. Many thanks to Norman Wang for maintaining HADB and providing data downloads. We thank the many phycologists who helped by providing samples or collecting assistance. Napua Harbottle, Neal Evenhuis, Jack Fisher and Roy Tsuda of the Bernice P. Bishop Museum are acknowledged for their assistance with sampling from collections at the Herbarium Pacificum.
- Abbott IA: Marine Red Algae of the Hawaiian Islands. 1999, Honolulu, Bishop Museum PressGoogle Scholar
- Sherwood AR: Bibliographic checklist of the non-marine algae of the Hawaiian islands. Bish Mus Occas Pap Ser. 2004, 80: 1-23.Google Scholar
- Huisman JM, Sherwood AR, Abbott IA: Studies of Hawaiian Galaxauraceae (Nemaliales, Rhodophyta): Large subunit rDNA gene sequences support conspecificity of Galaxaura rugosa and G. subverticillata. Crypto Algol. 2004, 25: 337-352.Google Scholar
- Huisman JM, Abbott IA, Sherwood AR: Large subunit rDNA gene sequences and reproductive morphology reveal Stenopeltis as a member of the Liagoraceae (Nemaliales, Rhodophyta), with a description of Akalaphycus gen. nov. Eur J Phycol. 2004, 39: 257-272. 10.1080/09670260410001710123.View ArticleGoogle Scholar
- Huisman JM, Abbott IA, Sherwood AR: The Liagoraceae (Nemaliales, Rhodophyta) of the Hawaiian Islands III: The genus Ganonema, with a description of Ganonema yoshizakii sp. nov. Phycologia. 2004, 43: 296-310. 10.2216/i0031-8884-43-3-296.1.View ArticleGoogle Scholar
- Huisman JM, Sherwood AR, Abbott IA: Morphology, reproduction, and the 18 S rRNA gene sequence of Pihiella liagoraciphila gen. et sp. nov., (Rhodophyta), the so-called 'monosporangial discs' associated with members of the Liagoraceae (Rhodophyta), and proposal of the Pihiellales ord. nov. J Phycol. 2003, 39: 978-987. 10.1046/j.1529-8817.2003.03047.x.View ArticleGoogle Scholar
- Saunders GW, Hommersand MH: Assessing red algal supraordinal diversity and taxonomy in the context of contemporary systematic data. Am J Bot. 2004, 91: 1494-1507. 10.3732/ajb.91.10.1494.PubMedView ArticleGoogle Scholar
- Kurihara A, Abe T, Tani M, Sherwood AR: Molecular phylogeny and evolution of red algal parasites: a case study of Benzaitenia, Janczewskia, and Ululania (Ceramiales). J Phycol. 2010, 46: 580-590. 10.1111/j.1529-8817.2010.00834.x.View ArticleGoogle Scholar
- Saunders GW: Applying DNA barcoding to red macroalgae: a preliminary appraisal holds promise for future application. Phil Trans Roy Soc Ser B. 2005, 360: 1879-1888. 10.1098/rstb.2005.1719.View ArticleGoogle Scholar
- Sherwood A, Kurihara A, Conklin KY: Molecular diversity of Amansieae (Ceramiales, Rhodophyta) from the Hawaiian Islands: a multi-marker assessment reveals high diversity within Amansia glomerata. Phycol Res.
- Schmitz F: 6. Klasse Rhodophyceae. 2. Unterklasse Florideae. Syllabus der Vorlesungen über specielle und medicinisch-pharmaceutische Botanik...Grosse Ausgabe. Edited by: Engler A. 1892, Gebr. Borntraeger, Berlin, 16-23.Google Scholar
- Kylin H: Die Gattungen der Rhodophyceen. 1956, Lund, Gleerups FörlagGoogle Scholar
- Pueschel CM, Cole KM: Rhodophycean pit plugs: an ultrastructural survey with taxonomic implications. Am J Bot. 1982, 69: 703-720. 10.2307/2442960.View ArticleGoogle Scholar
- Saunders GW, Bailey JC: Phylogenesis of pit-plug associated features in the Rhodophyta: inferences from molecular systematic data. Can J Bot. 1997, 75: 1436-1447. 10.1139/b97-858.View ArticleGoogle Scholar
- Ragan MA, Bird CJ, Rice EL, Gutell RR, Murphy CA, Singh RK: A molecular phylogeny of the marine red algae (Rhodophyta) based on the nuclear small-subunit rRNA gene. Proc Natl Acad Sci USA. 1994, 91: 7276-7280. 10.1073/pnas.91.15.7276.PubMedPubMed CentralView ArticleGoogle Scholar
- Kamiya M, Tanaka J, King RJ, West JA, Zuccarello GC, Kawai H: Reproductive and genetic distinction between broad and narrow entities of Caloglossa continua (Delesseriaceae, Rhodophyta). Phycologia. 1999, 38: 356-367. 10.2216/i0031-8884-38-5-356.1.View ArticleGoogle Scholar
- Zuccarello GC, West JA: Multiple cryptic species: Molecular diversity and reproductive isolation in the Bostrychia radicans/B. moritziana complex (Rhodomelaceae, Rhodophyta) with a focus on North American isolates. J Phycol. 2003, 39: 948-959. 10.1046/j.1529-8817.2003.02171.x.View ArticleGoogle Scholar
- Kumano S: Freshwater Red Algae of the World. 2002, Bristol, Biopress LtdGoogle Scholar
- Sherwood AR: The first report of a marine representative of the genus Hildenbrandia in Hawai'i. Haw Bot Soc Newsl. 2002, 41: 7-Google Scholar
- Coomans RJ, Hommersand MH: Vegetative growth and organization. Biology of the Red Algae. Edited by: Cole KM, Sheath RG. 1990, Cambridge, Cambridge University Press, 275-304.Google Scholar
- Guiry MD, West JA: Life history and hybridization studies on Gigartina stellata and Petrocelis cruenta (Rhodophyta) in the North Atlantic. J Phycol. 1983, 19: 474-494. 10.1111/j.0022-3646.1983.00474.x.View ArticleGoogle Scholar
- Guiry MD, West JA, Kim D-H, Masuda M: Reinstatement of the genus Mastocarpus Kützing (Rhodophyta). Taxon. 1984, 33: 53-63. 10.2307/1222029.View ArticleGoogle Scholar
- Garbary DJ, Gabrielson PW: Taxonomy and evolution. Biology of the Red Algae. Edited by: Cole KM, Sheath RG. 1990, Cambridge, Cambridge University Press, 477-498.Google Scholar
- Hebert PDN, Cywinska A, Ball SL, deWaard JR: Biological identifications through DNA barcodes. Proc Roy Soc Lond Ser B. 2003, 270: 313-322. 10.1098/rspb.2002.2218.View ArticleGoogle Scholar
- Stoeckle M: Taxonomy, DNA and the bar code of life. Bioscience. 2003, 53: 2-3. 10.1641/0006-3568(2003)053[0796:TDATBC]2.0.CO;2.View ArticleGoogle Scholar
- Moritz C, Cicero C: DNA barcoding: Promise and pitfalls. PLoS Biol. 2004, 2: 1529-1531. 10.1371/journal.pbio.0020354.View ArticleGoogle Scholar
- Conklin KY, Kurihara A, Sherwood AR: A molecular method for identification of the morphologically plastic invasive algal species Eucheuma denticulatum and Kappaphycus spp. (Rhodophyta, Gigartinales) in Hawaii. J Appl Phycol. 2009, 21: 691-699. 10.1007/s10811-009-9404-2.View ArticleGoogle Scholar
- Presting GG: Identification of conserved regions in the plastid genome - implications for DNA barcoding and biological function. Can J Bot. 2006, 84: 1434-1443. 10.1139/B06-117.View ArticleGoogle Scholar
- Sherwood AR, Presting GG: Universal primers amplify a 23 S rDNA plastid marker in eukaryotic algae and cyanobacteria. J Phycol. 2007, 43: 605-608. 10.1111/j.1529-8817.2007.00341.x.View ArticleGoogle Scholar
- Agosti D: Encyclopedia of life: Should species description equal gene sequence?. Trends Ecol Evol. 2003, 18: 273-10.1016/S0169-5347(03)00099-5.View ArticleGoogle Scholar
- Juvik SP, Juvik JO: Atlas of Hawaii. 1998, Honolulu, University of Hawaii Press, 3Google Scholar
- Zeigler AC: Hawaiian Natural History, Ecology, and Evolution. 2002, Honolulu, University of Hawaii PressGoogle Scholar
- Smith JE, Hunter CL, Smith CM: Distribution and reproductive characteristics of nonindigenous and invasive marine algae in the Hawaiian islands. Pac Sci. 2002, 56: 299-315. 10.1353/psc.2002.0030.View ArticleGoogle Scholar
- Webster JM, Clague DA, Riker-Coleman K, Gallup C, Braga JC, Potts D, Moore JG, Winterer EL, Paull CK: Drowning of the -150 m reef off Hawaii: A casualty of global meltwater pulse 1A?. Geology. 2004, 32: 249-252. 10.1130/G20170.1.View ArticleGoogle Scholar
- Vis ML, Sheath RG, Hambrook JA, Cole KM: Stream macroalgae of the Hawaiian Islands: A Preliminary Survey. Pac Sci. 1994, 48: 175-187.Google Scholar
- Filkin NR, Sherwood AR, Vis M: Stream macroalgae of the Hawaiian Islands: 23 Additional Collection Sites. Pac Sci. 2003, 57: 421-431. 10.1353/psc.2003.0030.View ArticleGoogle Scholar
- Sherwood AR, Vis ML, Sheath RG: Phenology and phylogenetic positioning of the Hawaiian endemic freshwater alga Batrachospermum spermatiophorum (Rhodophyta, Batrachospermales). Phycol Res. 2004, 52: 193-203. 10.1111/j.1440-1835.2004.tb00329.x.View ArticleGoogle Scholar
- Sherwood AR: Stream macroalgae of the Hawaiian islands: a floristic survey. Pac Sci. 2006, 60: 191-205. 10.1353/psc.2006.0011.View ArticleGoogle Scholar
- Chiasson WB, Johansen KG, Sherwood AR, Vis ML: Phylogenetic affinities of Chantransia pygmaea (Rhodophyta) specimens from the Hawaiian islands. Phycologia. 2007, 46: 257-262. 10.2216/06-79.1.View ArticleGoogle Scholar
- Vroom PS: Dasya atropurpurea sp. nov. (Ceramiales, Rhodophyta), a deep-water species from the Hawaiian Archipelago. Phycologia. 2005, 44: 572-580. 10.2216/0031-8884(2005)44[572:DASNCR]2.0.CO;2.View ArticleGoogle Scholar
- Saunders GW: Routine DNA barcoding of Canadian Gracilariales (Rhodophyta) reveals the invasive species Gracilaria vermiculophylla in British Columbia. Mol Ecol Res. 2009, 9: 140-150. 10.1111/j.1755-0998.2009.02639.x.View ArticleGoogle Scholar
- Sherwood AR, Vis ML, Entwisle TJ, Necchi O, Presting GG: Contrasting intra- versus inter-species DNA sequence variation for representatives of the Batrachospermales (Rhodophyta). Phycol Res. 2008, 56: 269-279. 10.1111/j.1440-1835.2008.00508.x.View ArticleGoogle Scholar
- Sherwood AR: Phylogeography of Asparagopsis taxiformis (Bonnemaisoniales, Rhodophyta) in the Hawaiian Islands: two mtDNA markers support three separate introductions. Phycologia. 2008, 47: 79-88. 10.2216/07-39.1.View ArticleGoogle Scholar
- Kurihara A, Sherwood AR: Checklist of the Hawaiian (marine, freshwater and terrestrial) red algae. Edition 2010.ix. Worldwide electronic publication, University of Hawaii, Honolulu, Hawaii, [http://algae.manoa.hawaii.edu]
- Andreakis N, Procaccini G, Maggs C, Kooistra WHCF: Phylogeography of the invasive seaweed Asparagopsis (Bonnemaisoniales, Rhodophyta) reveals cryptic diversity. Mol Ecol. 2007, 16: 2285-2299. 10.1111/j.1365-294X.2007.03306.x.PubMedView ArticleGoogle Scholar
- Sheath RG: The Biology of Freshwater Red Algae. Progress in Phycological Research. 1984, 3: 89-157.Google Scholar
- Sherwood AR, Sauvage T, Kurihara A, Conklin KY, Presting GG: A comparative analysis of COI, LSU and UPA marker data for the Hawaiian florideophyte Rhodophyta: implications for DNA barcoding of red algae. Crypto Algol.
- Sherwood AR, Chan Y, Presting GG: Application of universal plastid primers to environmental sampling of a Hawaiian stream periphyton community. Mol Ecol Res. 2008, 8: 1011-1014. 10.1111/j.1755-0998.2008.02138.x.View ArticleGoogle Scholar
- Clarkston BE, Saunders GW: A comparison of two DNA barcode markers for species discrimination in the red algal family Kallymeniaceae (Gigartinales, Florideophyceae), with a description of Euthora timburtonii sp. nov. Botany. 2010, 88: 119-131. 10.1139/B09-101.View ArticleGoogle Scholar
- Ratnasingham S, Hebert PDN: BOLD: The Barcode of Life Data System. Mol Ecol Notes. 2007, 7: 355-364. 10.1111/j.1471-8286.2007.01678.x. [http://www.barcodinglife.org]PubMedPubMed CentralView ArticleGoogle Scholar
- Verbruggen H, Maggs CA, Saunders GW, LeGall L, Yoon HS, DeClerck O: Data mining approach identifies research priorities and data requirements for resolving the red algal tree of life. BMC Evol Biol. 2010, 10: 16-10.1186/1471-2148-10-16.PubMedPubMed CentralView ArticleGoogle Scholar
- Hart MW, Sunday J: Things fall apart: biological species form unconnected parsimony networks. Biol Lett. 2007, 3: 509-512. 10.1098/rsbl.2007.0307.PubMedPubMed CentralView ArticleGoogle Scholar
- LeGall L, Saunders GW: DNA barcoding is a powerful tool to uncover algal diversity: a case study of the Phyllophoraceae (Gigartinales, Rhodophyta) in the Canadian flora. J Phycol. 2010, 46: 374-389. 10.1111/j.1529-8817.2010.00807.x.View ArticleGoogle Scholar
- Ekrem T, Willassen E, Stur E: A comprehensive DNA sequence library is essential for identification with DNA barcodes. Mol Phyl Evol. 2007, 43: 530-542. 10.1016/j.ympev.2006.11.021.View ArticleGoogle Scholar
- Zhang AB, He LJ, Crozier RH, Muster C, Zhu C-D: Estimating sample sizes for DNA barcoding. Mol Phylo Evol. 2010, 54: 1035-1039. 10.1016/j.ympev.2009.09.014.View ArticleGoogle Scholar
- The Hawaiian Algal Database. [http://algae.manoa.hawaii.edu]
- Global Biodiversity Information Facility. [http://www.gbif.org]
- National Center for Biotechnology Information GenBank. [http://www.ncbi.nlm.nih.gov/]
- Wang N, Sherwood AR, Kurihara A, Conklin KY, Sauvage T, Presting GG: The Hawaiian Algal Database: a laboratory LIMS and online resource for biodiversity data. BMC Plant Biol. 2009, 9: 117-10.1186/1471-2229-9-117.PubMedPubMed CentralView ArticleGoogle Scholar
- GPS Visualizer. [http://www.gpsvisualizer.com/]
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignments through sequence weighting, position specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed CentralView ArticleGoogle Scholar
- Tamura K, Nei Dudley, Kumar MS: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Paradis E, Claude J, Strimmer K: APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004, 20: 289-290. 10.1093/bioinformatics/btg412.PubMedView ArticleGoogle Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria, [http://www.R-project.org]
- Stamatakis A, Hoover P, Rougemont JR: A Rapid Bootstrap Algorithm for the RAxML Web Servers. Syst Biol. 2008, 57: 758-771. 10.1080/10635150802429642.PubMedView ArticleGoogle Scholar
- Stamatakis A, Ludwig T, Meier H: RAxML-III: A fast program for maximum likelihood-based inference of large phylogenetic trees. Bioinformatics. 2005, 21: 456-463. 10.1093/bioinformatics/bti191.PubMedView ArticleGoogle Scholar
- Templeton AR, Crandall KA, Singh CF: A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992, 132: 619-633.PubMedPubMed CentralGoogle Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1660. 10.1046/j.1365-294x.2000.01020.x.PubMedView ArticleGoogle Scholar