pmc logo imageJournal ListSearchpmc logo image
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS ONE. 2008; 3(3): e1835.
Published online 2008 March 19. doi: 10.1371/journal.pone.0001835.
PMCID: PMC2266797
Ancient DNA Elucidates the Controversy about the Flightless Island Hens (Gallinula sp.) of Tristan da Cunha
Dick S. J. Groenenberg,1* Albert J. Beintema,2 René W. R. J. Dekker,1 and Edmund Gittenberger1,3
1National Museum of Natural History Naturalis, Leiden, The Netherlands
2Biowrite, Gorssel, The Netherlands
3Institute of Biology, Leiden University, Leiden, The Netherlands
Jean-Nicolas Volff, Editor
Ecole Normale Supérieure de Lyon, France
* E-mail: groenenberg/at/naturalis.nl
Conceived and designed the experiments: DG. Performed the experiments: DG. Analyzed the data: DG. Wrote the paper: DG. Other: Provided a hypothesis and some of the specimens: AB. Reviewed the manuscript: EG RD AB.
Received October 3, 2007; Accepted February 17, 2008.
Abstract
A persistent controversy surrounds the flightless island hen of Tristan da Cunha, Gallinula nesiotis. Some believe that it became extinct by the end of the 19th century. Others suppose that it still inhabits Tristan. There is no consensus about Gallinula comeri, the name introduced for the flightless moorhen from the nearby island of Gough. On the basis of DNA sequencing of both recently collected and historical material, we conclude that G. nesiotis and G. comeri are different taxa, that G. nesiotis indeed became extinct, and that G. comeri now inhabits both islands. This study confirms that among gallinules seemingly radical adaptations (such as the loss of flight) can readily evolve in parallel on different islands, while conspicuous changes in other morphological characters fail to occur.
Introduction

Until recently it was assumed that the flightless moorhen of remote Tristan da Cunha in the southern Atlantic (Fig. 1), Gallinula nesiotis (Sclater, 1861) [1], became extinct by the end of the 19th century [2]. A few decades after its description, a very similar moorhen that was also flightless namely G. comeri (Allen, 1892), was described [3] from the island of Gough, ca. 400 km SE of Tristan. In the period between these descriptions G. nesiotis became rare [4], [5] and by the turn of the century it had probably gone extinct [7], [6]. Authentic remnants are two skins and a skeleton in the Natural History Museum, Tring [8]. Since unequivocal G. nesiotis had been collected only once from Tristan, and because of the presence of a healthy population of similar moorhens on the nearby island of Gough, some authors doubted whether an endemic moorhen had ever existed on Tristan [9]. Eber [10] compared Sclater's description of G. nesiotis from Tristan with her series of G. comeri from Gough and concluded that the differences fall within the range of variation of the latter. In her opinion it was very unlikely that moorhens from two islands in the same region would have independently lost the ability of flight, without differentiating in other characters. She suggested that Sclater's material might have been labelled inaccurately and that his specimens in fact also came from Gough. Consequently, Eber considered G. comeri a junior synonym of G. nesiotis [10] and controversy surrounded future illustrations of both taxa (Fig. 2).

Figure 1Figure 1
Map showing the location of Tristan da Cunha and Gough in the Mid Atlantic.
Figure 2Figure 2
Illustrative stamps, issued in 1987 and 2005.

Beintema [2] mentioned that there are old records of moorhens for Tristan, demonstrating that such birds were truly indigenous there. In his view, skeletal measurements differ slightly between G. nesiotis and G. comeri. Furthermore, rails are known to rapidly lose the ability of flight as soon as they arrive on remote islands [11]. When in 1972 live moorhens were discovered on Tristan [12], these birds were regarded as descendants of a small number of individuals brought from Gough [13]. Alternatively, Beintema suggested that G. nesiotis might have been temporarily rare on Tristan, but not extinct, and that the moorhens found there today are descendents of the original island population. Here we address this question, making use of DNA analyses of authentic material of G. nesiotis, recent specimens of the moorhens from both Tristan and Gough and some geographically and taxonomically close other taxa of moorhens.

Results

An overview of the specimens that were used in this study, with taxon names, locality data, and year of acquisition, is given in Table 1. Alignments of all cloned sequences of G. nesiotis (two independent amplifications per target region) are shown in Dataset S1, S2, S3. On the one hand, none of the sequences of genuine, historical G. nesiotis was identical to those of G. comeri and, on the other hand, all sequences of the moorhens collected in Tristan da Cunha in 1993, are identical to those found for specimens of G. comeri from Gough, dated 1960. The genetic distances between G. nesiotis and the other gallinules, are of the same magnitude as the distances between G. comeri and the other gallinules (Table 2). A pairwise relative rates test did not reveal significantly different substitution rates for any of the lineages (P>0.18). Of the selected markers, most variation was detected in the control region (D-loop) sequences, viz. 9.6% TSH for G. nesiotis compared with G. comeri ( = Total Sequence Heterogeneity [15], [14]) versus 2.1% and 0.3% TSH for tRNA-Lys/ATP synthase subunit 8 (ATP8) and cytochrome b, respectively.

Table 1Table 1
Taxa and collection information.
Table 2Table 2
Genetic distances.

The results of phylogenetic analyses (Neighbour-Joining, Maximum Likelihood and Bayes) based on a combined dataset (all taxa, all regions, Dataset S4) are shown in Fig. 35. In these cladograms Gallinula nesiotis and G. comeri form a clade with the moorhens of Africa/Eurasia, whereas the other taxa that were investigated are less closely related.

Figure 3Figure 3
Bootstrap 50% majority rule consesus NJ tree.
Figure 4Figure 4
Bootstrap 50% majority rule consensus ML tree.
Figure 5Figure 5
Bayes 50% majority rule consensus (‘halfcompat’) phylogram (branch lengths are proportional to the expected number of substitutions).
Discussion

Our results show that genuine G. nesiotis, identified on the basis of historical material from the island of Tristan da Cunha, differs genetically from G. comeri, which has been described from the island of Gough. For each marker the sequence of G. nesiotis differs from that of G. comeri, as well as from all other Gallinula taxa that were analysed, but the amount of variation differs strongly between the regions studied. The position of G. nesiotis on the cladograms (Fig. 35) makes sense biologically. Apparently, G. nesiotis became extinct on Tristan and G. comeri from Gough was introduced there, resulting in the current situation with G. comeri occurring on both islands. This implies that modern illustrations of so-called G. nesiotis from Tristan (Fig. 2) probably show introduced G. comeri from Gough.

The difference between G. nesiotis and G. comeri is most conspicuous in the D-loop sequence. This is a marker from a non-coding region, which makes it more difficult to exclude it as a potential pseudogene [16]. In some cases preferential amplification of numt (nuclear mitochondrial insertion) sequences has been observed [17], [18], but in most ancient DNA studies only occasional co-amplification of numts has been reported [19], [20]. This is not surprising since (particularly) in ancient DNA samples, mitochondrial DNA will be in excess over nuclear DNA. Consequently the incidence of numts should be reduced in such samples. To minimize the chance of amplifying numts, we did not use blood as a source for DNA-extractions in our recently collected material [17], [21] (in the older specimens it was not possible anyway). Because ‘universal’ primers may also be particularly prone to amplification of numts [17], the primers for the D-loop were made ‘gallinule-specific’. They did not even work for the closely related coot, Fulica atra. Products that were sequenced directly (both strands) showed only one signal, whereas multiple signals can be expected if both the target product and a numt would have been amplified. Interclone variation was low (Dataset S1, S2, S3). Only one sequence (G. nesiotis, marker ATP8) out of 96 clones (Table 3) could clearly be identified as a numt (Dataset S2). No stop-codons or frame-shift mutations were observed for the coding-region datasets (ATP8 and cytochrome b). No obvious deviations in either substitution rate (pairwise relative rate test) or base composition, like a decrease in GC content [22], [23], were observed.

Table 3Table 3
Number of colonies sequenced per target region per taxon.

All sequences of recently collected moorhens from Tristan were identical to those of G. comeri from Gough and should be considered conspecific therefore. Cross contamination is very unlikely, since specimens from Tristan and Gough were amplified in different PCR-batches and contamination was not detected in other, partly much older specimens. Most probably the sequences are identical because G. comeri was introduced only recently on Tristan. Genetic variation within island populations is generally small compared to mainland populations [24][26]. For example, the giant tortoises (Aldabrachelys) of Mahé (Seychelles) and Mauritius (Mascare islands) still have identical sequences compared to those of Aldabra, from where they were shipped since the 1820s or earlier [27]. G. comeri may have been introduced on Tristan somewhere in the mid 1950s [12]. Assuming a generation time of two years, as known for Gallinula chloropus, G. comeri would only have had about only 20 generations (40 years) to differentiate on Tristan.

Both the genetic distances (Table 2) and the fact that G. nesiotis and G. comeri form a clade with the investigated moorhens of Africa/Eurasia (Fig. 35), suggest that the ancestor(s) of these island gallinules originated from Africa and not America, as suggested by Eber [10]. Our data do not allow us to distinguish between a single dispersal event to the archipelago, followed by allopatric differentiation, or two separate introductions from the continent to both Tristan and Gough.

As is inevitably the case with isolated island populations, the question of whether G. nesiotis and G. comeri were reproductively isolated under natural circumstances cannot be answered. Our limited data from a small number of specimens and sequences of only the mitochondrial lineage are insufficient to demonstrate hybridisation. Even though island populations generally show lower genetic variation than related mainland populations [24], the genetic distances between G. nesiotis and G. comeri are of at least the same magnitude as those found between taxa that figure as subspecies of G. chloropus in the literature (Fig. 5, Table 2). Therefore, we propose that the extinct moorhen of Tristan and the moorhens that live on Gough and Tristan today be regarded as subspecies, viz. G. n. nesiotis and G. n. comeri, respectively. This is in conformity with two recent, general checklists of the birds of the world [28], [29] and a detailed monograph of the rails of the world [13], but is different from a morphological study by Eber [10] in which both taxa are considered synonyms.

Materials and Methods

Taxa
Tissues from fourteen gallinules and a coot were put at our disposal by various institutes (Table 1). These include tissues of (I) ‘recently’ collected moorhens from Tristan da Cunha, (II) moorhens from Gough from the collection ZMA, (III) the 1864 specimen of G. nesiotis from the Natural History Museum, Tring, and (IV) a number of subspecies of Gallinula chloropus from South-America, Africa, Europe, Taiwan and Java from the National Museum of Natural History Naturalis, Leiden. The coot, Fulica atra, was used as outgroup.

DNA extraction
DNA extractions on specimens from 1968 (Table 1) and older were carried out in a dedicated aDNA facility (LAF, Leiden, the Netherlands), which is physically isolated from the main laboratories. Before extractions took place, the extraction room was cleaned with a 0.05% bleach solution and the extraction-cabinet was decontaminated by turning on the UV lights at least 1 hour prior to the start of the extractions. No more than four extractions were done at once and negative controls were included with each set of extractions. Pippetes were cleaned with bleach and subsequently decontaminated (together with the dispossables) by UV irradiation (UV linker). Tissues were cut into small pieces to enlarge the contact surface between tissue and buffer. Total genomic DNA was extracted with a DNeasy Tissue Kit (Qiagen) using a prolonged incubation (24 hours). Proteinase K was added twice, once at the start and after 6 hours of incubation. To concentrate the extract, elution volume was decreased to 40 µl. Extractions on recently collected specimens (1993–2005) were done in a common lab, also using a DNeasy Tissue Kit (Qiagen).

PCR and sequencing
PCRs were never performed in the aDNA facility and amplicons were never stored in this building. PCRs on the extract of the 1864 specimen were duplicated in different laboratories that are physically separated from each other as well as from the aDNA facility. In none of these labs had ever been worked on any species of gallinule before. Fragments from three non-adjacent mitochondrial gene regions (679 basepairs in total; primersites excluded) were amplified by PCR: the D-loop, tRNA-Lysine/ATP8 and cytochrome b. The length of these fragments (primer sites included) was 234, 236 and 375 bp, respectively. Primer sequences and references are described in Table 4. For G. nesiotis and (most) specimens of 1960 and older, cytochrome b could not be amplified directly using primers L14841 and H15149 [30]. Presumably because the DNA within these specimens got too degraded over time. Therefore, internal primers were designed (Table 4) to amplify this fragment in two overlapping parts: L14841- Rev219 (219 bp) and Fwd141- H15149 (249 bp). The primers for the D-loop (CR-OUD-F and CR-OUD-R) are ‘gallinule-specific’. For the coot, Fulica atra, the same region had to be amplified with other primers: CR-175-F and 12S-29-R (Table 4). PCRs were done using a standard Taq DNA Polymerase kit (Qiagen). Reaction volume was 25 µl and PCR conditions were 0.4 µM of each primer, 0.2 mM dNTP's and 5 units of Taq DNA Polymerase. For amplification of the cytochrome b and ATP8 regions, the final concentration of MgCl2 was 2.5 mM. For amplification of the D-loop fragment no MgCl2 was added (1.5 mM was already in the Qiagen PCR-buffer). Thermocycling conditions were 3 min. at 94°C (initial denaturation), followed by 40 cycles (15 sec. at 94°C, 30 sec. at AT°C and 40 sec. at 72°C) and final extension 5 min. at 72°C. Where AT is the anealling temperature for each primerset; 50°C for both cytochrome b and ATP8, 55°C for the D-loop fragment and 57°C for reamplification of cloned products (see below).
Table 4Table 4
PCR and sequencing primers.

All PCR products from G. nesiotis and a number of PCR product from the D-loop of selected taxa (Table 3) were cloned using either pGEM®-T Easy Vector Sytems from Promega or Topo TA Cloning® from Invitrogen. At least three colonies were picked per plate and used to initiate reamplifications with primers 21M13_F and 21M13_R (Table 4). Reamplified products were cleaned using a Nucleospin® kit (Macherey-Nagel). Subsequently these products were sequenced either in-house on a Megabace™ 1000 DNA Analysis System (Amersham), or on a 3730xl DNA analyzer (Applied Biosystems) at Macrogen Inc. (Korea) using only primer 21M13_F. All other PCR products were cleaned (same procedure) and sequenced directly (both directions) with their respective PCR primers (Table 4). A summary of the specimens and the number of colonies sequenced per target region is given in Table 3. Sequences were assembled using Sequencher version 4.2 (Gene Codes Corporation) and aligned manualy using MacClade version 4.08 [31]. The sequences were deposited in GenBank (accession numbers EF681971-EF682015).

Phylogenetic analysis
For phylogenetic analyses, all sequences (all regions; consensus sequences when products were cloned) were put in a single datamatrix (Dataset S4; an ILD-test showed no incongruence between the regions, p = 0.971) and Fulica atra was designated as outgroup. To get branch support values, we performed phylogeneticanalyses with three methods: Neighbour-Joining (PAUP ver. 4.0b2a [32]), Maximum Likelihood (PAUP ver. 4.0b2a [32]) and Bayesian analysis (Mr.Bayes ver. 3.1.2 [33]). For the NJ analysis, we performed a bootstrap analysis (1000 replicates, optimality criterion set to distance) and calculated a 50% majority rule consensus cladogram (Fig. 3). For the ML analysis, the HKY+G model was selected by Modeltest ver. 3.7 [34] with the following parameters: Tratio = 10.660, gamma shape parameter = 0.0941, base frequencies A = 0.3473, C = 0.3004, G = 0.1341, T = 0,2182 and proportion of invariable sites (pinvar) = 0. A bootstrap analysis (1000 replicates, 5 random additions per bootstrap replicate and TBR branch swapping) was performed and a 50% majority rule consensus tree was calculated (Fig. 4). For the MrBayes analysis, the best-fit model for each partition (four partitions: tRNA-Lysine and ATP8 were considered as two partitions) was selected by hLRT in MrModeltest ver. 2.2 [35]: D-loop (HKY+Y), tRNA-lysine (HKY), ATP8 (GTR+I) and cytochrome b (GTR+I). A dirichlet (1,1,1,1) prior was specified on the state frequencies for all partitions, except for the tRNA-Lysine partition, where the frequencies were equal. All partitions had different rates for transition and transversions (nst = 2), except for tRNA-Lysine (nst = 1). Among-site rate variation was equal for tRNA-Lysine, gamma-distributed for both D-loop and ATP8 and for cytochrome b a proportion of the sites was invariant. Two runs (set up for 10 000 000 generations) were performed simultaneously (4 chains per run) in MrBayes ver. 3.1.2 [34] and the convergence diagnostic was set to 0.009. A Markov chain Monte Carlo (MCMC) analysis was done with swapfreq = 2, temp = 0.002 and samplefreq = 100; convergence was reached after 165 000 generations. The trees of both runs (3302 in total) were combined (2702: burnin was set to 300) and a 50% majoritiy rule consensus tree (contype = halfcompat) was calculated (Fig. 5).

Genetic distances
Genetic distances (absolute number of changes and uncorrected “p” distances) were calculated with Paup ver. 4.0b2a [32] based on the combined dataset (Dataset S4).

Pairwise relative rates test
With Fulica atra specified as outgroup, a Pairwise Relative Rate Test [36] as implemented in HyPhy [37] using the HKY model (as specified by Modeltest ver. 3.7 [34]) was performed on the combined dataset (Dataset S4).

Supporting Information
Dataset S1

Cloning results of G. nesiotis for ATP8.

(0.00 MB TXT)

Dataset S2

Cloning results of G. nesiotis for D-loop.

(0.00 MB TXT)

Dataset S3

Cloning results of G. nesiotis for cytochrome b.

(0.01 MB TXT)

Dataset S4

Combined dataset showing all taxa and all markers (D-loop, ATP8 and cytochrome b) used in this study.

(0.01 MB TXT)

Acknowledgments

We thank H. van Grouw (National Museum of Natural History Naturalis, Leiden, The Netherlands), Dr R. Prys-Jones (The Natural History Museum, Tring, U.K.) and T. Prins (Zoological Museum, University of Amsterdam, the Netherlands) who provided most valuable material for DNA analyses. We also thank H. Cross (National Herbarium the Netherlands, Leiden) for technical assistance, K. Vrieling (Leiden University, the Netherlands) for putting at our disposal the Topo TA Cloning® System (Invitrogen) and P. Brakefield (Leiden University, the Netherlands) for reviewing our manuscript. Finally we acknowledge the Tristan da Cunha Post Office, who gave permission to illustrate this manuscript with some of their stamps (Fig. 1) and we acknowledge the Crown Agents Stamp Bureau for facilitating correspondence.

Footnotes
Competing Interests: The authors have declared that no competing interests exist.
Funding: The authors have no support or funding to report.
References
1.
Sclater, PL. On the Island-Hen of Tristand d'Acunha. Proceedings of the Zoological Society Londen. 1861:260–263.
2.
Beintema, AJ. The history of the Island Hen (Gallinula nesiotis), the extinct flightless gallinule of Tristan da Cunha. Bulletin of the British Ornithologists' Club. 1972;92:106–113.
3.
Allen, JA. Description of a new Gallinule from Gough island. Bulletin of the American Museum of Natural History. 1892;4:57–58.
4.
Milner, J; Brierly, OW. London: W.H. Allen; 1869. The Cruise of H.M.S. ‘Glatea’, Captain H.R.H. the Duke of Edinburgh, K.G., in 1867–1868.
5.
Sperling, RM. Letter on Tristan da Cunha. Ibis. 1872;3:74–79.
6.
Nicoll, MJ. On the Birds collected and observed during the Voyage of the ‘Valhalla’. 1906;8:666–712.
7.
Nicoll, MJ. London: R.Y.S. Witherby & Co; 1908. Three Voyages of a Naturalist; being an Account to Many Little-Known Islands in the Three Oceans visited by the ‘Valhalla’.
8.
Knox, AG; Walters, MP. Extinct and endangered birds in the collections of the Natural History Museum. British Ornithologist's Club Occasional Publications. 1994;1
9.
Broekhuysen, GJ; Macnae, B. Obersvations on the birds of Tristan da Cunha Islands and Gough Island in February and early March, 1948. Ardea. 1949;37:97–113.
10.
Eber, G. Vergleichende untersuchungen am flugfähigen Teichhuhn Gallinula chl. chloropus und an der flugunfähigen Inselralle Gallinula nesiotis. Bonner Zoologische Beitrage. 1961;12:247–315.
11.
Slikas, B; Olson, SL; Fleischer, RC. Rapid, independent evolution of flightlessness in four species of Pacific Island rails (Rallidae): an analysis based on mitochondrial sequence data. Journal of Avian Biology. 2002;33:5–14.
12.
Richardson, ME. Aspects of the ornithology of the Tristan da Cunha group and Gough Island, 1972–1974. Cormorant. 1984;12:123–201.
13.
Taylor, B; van Perlo, B. Sussex: Pica Press; 1998. Rails: A guide to the Rails, Coots, Crakes and Gallinules of the world.
14.
Binladen, J; Wiuf, C; Gilbert, MTP; Bunce, M; Barnett, R, et al. Assessing the fidelity of ancient DNA sequences amplified from nuclear genes. Genetics Genetics. February 2006;172:733–741.
15.
Gilbert, MTP; Willerslev, E; Hansen, AJ; Barnes, I; Rudbeck, L, et al. Distribution patterns of postmortem damage in human mitochondrial DNA. Am J Hum Genet. January 2003;72:32–47. [PubMed]
16.
Perna, NT; Kocher, TD. Mitochondrial DNA: molecular fossils in the nucleus. Curr Biol. 1996;6:128–129. [PubMed]
17.
Sorenson, MD; Quinn, TW. Numts: A Challenge for Avian Systematics and Population Biology. The Auk. 1998;115:214–221.
18.
Collura, RV; Stewart, CB. Insertions and duplications of mtDNA in the nuclear genomes of Old World monkeys and hominoids. Nature. 1995;378:485–489. [PubMed]
19.
Kolokotronis, SO; MacPhee, RDE; Greenwood, AD. Detection of mitochondrial insertions in the nucleus(NuMts) of Pleistocene and modern muskoxen. BMC Evolutionary Biology. 2007;7:67. [PubMed]
20.
Orlando, L; Leonard, JA; Thenot, A; Laudet, V; Guerin, C, et al. Ancient DNA analysis reveals woolly rhino evolutionary relationships. Molecular Phylogenetics and Evolution. 2003;28:485–499. [PubMed]
21.
Quinn, TW. The genetic legacy of Mother Goose- phylogeographic patterns of lesser snow goose Chen caerulescens caerulescens maternal lineages. Molecular Ecology. 1992;1:105–117. [PubMed]
22.
Casane, D. Mutation Pattern Variation Among Regions of the Primate Genome. Journal of Molecular Evolution. 1997;45:216–226. [PubMed]
23.
Bensasson, D; Zhang, DX; Hartl, DL; Hewitt, GM. Mitochondrial pseudogenes: evolution's misplaced witnesses. Trends in Ecology & Evolution. 2001;16:314–321. [PubMed]
24.
Frankham, R. Do island populations have less genetic variation than mainland populations? Heredity. 1997;78:311–327. [PubMed]
25.
Frankham, R. Inbreeding and Extinction: Island Populations. Conservation Biology. 1998;12:665–675.
26.
Eldridge, MDB; King, JM; Loupis, AK; Spencer, PBS; Taylor, AC, et al. Unprecedented Low Levels of Genetic Variation and Inbreeding Depression in an Island Population of the Black-Footed Rock-Wallaby. Conservation Biology. 1999;13:531–541. doi:10.1046/j.1523-1739.1999.98115.x.
27.
Austin, JJ; Arnold, EN; Bour, R. Was there a second adaptive radiation of giant tortoises in the Indian Ocean? Using mitochondrial DNA to investigate speciation and biogeography of Aldabrachelys (Reptilia, Testudinidae). Molecular Ecology. 2003;12:1415–1424. doi:10.1046/j.1365-294X.2003.01842.x. [PubMed]
28.
del Hoyo, J; Elliott, A; Sargatal, J. Barcelona: Lynx Edicions; 1996. Handbook of the Birds of the World. 3. Hoatzin to Auks.
29.
Dickinson, EC. The Howard and Moore Complete Checklist of the Birds of the World. Christopher Helm. 2003
30.
Kocher, TD; Thomas, WK; Meyer, A; Edwards, SV; Paabo, S, et al. Dynamics of Mitochondrial-DNA Evolution in Animals-Amplification and Sequencing with Conserved Primers. Proc. Natl. Acad. Sci. U. S. A. 1989;86:6196–6200. [PubMed]
31.
Maddison, WP; Maddison, DR. MacClade: analysis of phylogeny and character evolution, version 4.8. Sinauer, Sunderland. 2005
32.
Swofford, DL. Sunderland, MA: Sinauer Associates; 1998. PAUP*. Phylogenetic Analysis Using Parsimony (* and Other Methods). Version 4.
33.
Ronquist, F; Huelsenbeck, JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. August 12, 2003;19:1572–1574. doi:10.1093/bioinformatics/btg180. [PubMed]
34.
Posada, D; Crandall, KA. MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998;14:817–818. doi:10.1093/bioinformatics/14.9.817. [PubMed]
35.
Nylander, JAA. MrModeltest 2.2. Program distributed by the author. Uppsala University: Evolutionary Biology Centre. 2004
36.
Muse, SV; Weir, BS. Testing for Equality of Evolutionary Rates. Genetics. 1992;132:269–276. [PubMed]
37.
Pond, SLK; Frost, SDW; Muse, SV. Oxford Univ Press; 2005. HyPhy: hypothesis testing using phylogenies. pp. 676–679.