A highly polymorphic insertion in the Y-chromosome amelogenin gene can be used for evolutionary biology, population genetics and sexing in Cetacea and Artiodactyla
© Macé and Crouau-Roy; licensee BioMed Central Ltd. 2008
Received: 28 September 2007
Accepted: 16 October 2008
Published: 16 October 2008
The early radiation of the Cetartiodactyla is complex, and unambiguous molecular characters are needed to clarify the positions of hippotamuses, camels and pigs relative to the remaining taxa (Cetacea and Ruminantia). There is also a need for informative genealogic markers for Y-chromosome population genetics as well as a sexing method applicable to all species from this group. We therefore studied the sequence variation of a partial sequence of the evolutionary conserved amelogenin gene to assess its potential use in each of these fields.
Results and discussion
We report a large interstitial insertion in the Y amelogenin locus in most of the Cetartiodactyla lineages (cetaceans and ruminants). This sex-linked size polymorphism is the result of a 460–465 bp inserted element in intron 4 of the amelogenin gene of Ruminants and Cetaceans. Therefore, this polymorphism can easily be used in a sexing assay for these species.
When taking into account this shared character in addition to nucleotide sequence, gene genealogy follows sex-chromosome divergence in Cetartiodactyla whereas it is more congruent with zoological history when ignoring these characters. This could be related to a loss of homology between chromosomal copies given the old age of the insertion.
The 1 kbp Amel-Y amplified fragment is also characterized by high nucleotide diversity (64 polymorphic sites spanning over 1 kbp in seven haplotypes) which is greater than for other Y-chromosome sequence markers studied so far but less than the mitochondrial control region.
The gender-dependent polymorphism we have identified is relevant not only for phylogenic inference within the Cetartiodactyla but also for Y-chromosome based population genetics and gender determination in cetaceans and ruminants. One single protocol can therefore be used for studies in population and evolutionary genetics, reproductive biotechnologies, and forensic science.
About 240 to 320 million years ago, shortly after the divergence of mammalian and avian lineages, progressive X-Y differentiation began, following chromosomal interstitial rearrangements . This resulted in a partial loss of homology between both chromosomes which reached its maximal extent in primates . Amelogenin is the enamel matrix protein that combines with hydroxyapatite crystals to form enamel prisms in teeth . The gene encoding the amelogenin protein (Amel) is among the few genes expressed from both X and Y chromosomes in placental mammals (Eutheria) .
Evolutionary uncertainties about the basal diversification of Cetartiodactyla
The Cetartiodactyla (even-toed ungulates, whales and dolphins) radiated approximately 70–80 Myrs ago. The relative positions of the camelid, suiform (pigs), hippopotamus, ruminant and cetacean (whales and dolphins) groups remain unclear, whether morphological or molecular characters are used for attribution [5–9]. Moreover, polytomies (unresolved tree nodes) within some Cetartiodactyla taxa  highlight areas for further data collection (both species and markers) and phylogenetic research. This is a particularly delicate problem within cetaceans, probably due to adaptative radiations within a short period of time [10, 11].
Y-chromosome sequence markers are needed for population genetics
Males are the heterogametic sex in mammals, and usually, unequal numbers of males and females transmit genes from one generation to the next. Y-specific polymorphisms should allow the inference of sex-specific population parameters and decryption of breeding system patterns and dispersal strategies. Overall, the use of Y-specific markers has been restricted to evolutionary studies of human history and some scarce studies in population genetics, perhaps because of the low diversity of these markers . Within Cetartiodactyla, genetic structure or admixture, e.g. in sheep  or cattle [14, 15], has made use of a few Y-specific markers including microsatellites, SNPs and indels.
The matrilineally transmitted mitochondrial control region is commonly used as an informative sequence for population genetics. An equivalent had not been found to date on the Y chromosome. We considered the well-known amelogenin gene to be of particular interest because parts of it do not recombine between X and Y chromosomes.
Sex-chromosome recombination discrepancies have been exploited to develop many molecular sexing techniques. Although it can be ambiguous in some small populations, the amelogenin locus is the most commonly used for gender determination in humans . Accurate gender determination in mammals is crucial to various applications in reproductive technologies, forensic investigations and population management. Some techniques rely on specific amplification of loci localized on the Y chromosome (such as Sry ) while others are based on amplification of homologous fragments from both X and Y chromosomes (e.g. ZFX/ZFY ) or use both markers . Each of these has limitations, such as the need for multiplexing with other markers or additional steps such as digestion, labelling or sequencing. Several amelogenin-based techniques have expanded the taxonomic coverage of molecular sexing for Artiodactyla [20–22] but they have not yet been extended to Cetacea. Therefore, there is a need both for new methods that apply to a greater number of species and to increase the number of cross-checking sexing methods, especially in conservation biology .
In this study, we found out that sequence variations in the amelogenin locus can be applied in evolutionary and population genetics as well as for molecular sexing in the highly diversified Cetartiodactyla group. We therefore carried out an evolutionary study of orthologous Amel-Y and Amel-X sequences (exons 4 to 5) in Cetartiodactyla. We studied four Cetacea (the striped dolphin, the bottlenosed dolphin, Risso's dolphin and the fin whale) and three Artiodactyla (cow, pig and sheep) species.
Results and Discussion
Amelogenin can be used for molecular sexing and evolutionary genetics in Cetartiodactyla
List of Amel-X and Amel-Y haplotype names in Cetaceans and their EMBL accession numbers
Stenella cœruleoalba YA1
Stenella cœruleoalba YA2
Stenella cœruleoalba YB1
Stenella cœruleoalba YB2
Stenella cœruleoalba YB3
Stenella cœruleoalba YB4
Stenella cœruleoalba YB5
Delphinus delphis Y1
Delphinus delphis Y2
Stenella cœruleoalba X
Grampus griseus Y
Balænoptera physalus Y
Balænoptera physalus X
Eschrichtius robustus X
Tursiops truncatus X
It would be interesting to study this region at the whole clade level by combining sequence and indel characters in the same analysis. This could give clues to test the many hypotheses about basal radiation of Cetartiodactyla (e.g. [5, 6, 8]). Given the presumably basal position of the Suioidea and Tylopoda in the Cetartiodactyla phylogeny ( and Fig. 3), we hypothesize that the major evolutionary event represented by the insertion (illustrated by an arrow Figure 4a) occurred once in the Cetacea-Ruminantia clade and not in the remaining Cetartiodactyla.
The presence of this large insertion in the Amel-Y copy can be useful for sex determination.
Evolutionary history also indicates that our sexing technique is applicable, in addition to cetaceans, to over a wide range of Cetartiodactyla species including domestic and wild species, in particular the widespread Ruminantia (Bovidae, Capridae and most likely Cervidae). It is however not suitable to Suiformes and further studies are required to confirm that the technique is also not applicable to Camelidae, given their even more basal position in the Cetartiodactyla phylogeny.
Use in pedigree assessment and population genetics
In dolphins, the Amel-Y fragments amplified with the SC1-SC2 primer pair were easily sequenced without the need of cloning since amplification was Y chromosome-specific. From the ten Striped dolphin samples sequenced, nine were males, and we could deduce seven distinct Y-haplotypes (one haplotype represented by three individuals and four individual haplotypes) bearing 64 polymorphic sites (nucleotide diversity π = 0.004 ± 0.0007). Half of these were in the ~460 bp insertion. An alignment of polymorphic sites is presented in Figure 5 (a). Strikingly, these sequences showed two highly divergent haplogroups, diverging by a mean of 49 substitutions. This concords with our results that support the probable existence of two subspecies within the Mediterranean sea (unpublished data). Moreover, one of these haplogroups displayed a high degree of polymorphism, with 24 informative sites, whereas the others showed only eight. These values are sufficient for use in pedigree analysis and population genetics, as the Y chromosome counterpart of the mitochondrial d-loop in this species. Indeed, in striped dolphin the intra-specific (inter-group) divergence is greater than inter-specific divergence with a mean of 45 nucleotide substitutions between the striped dolphin and fin whale sequences. There is an average of 0.048 ± 0.01 substitutions per site when comparing the two striped dolphin populations. This is comparable to the divergence observed between each population and the Common dolphin (0.058 ± 0.03) and confirms that nucleotide diversity is one order of magnitude higher than the range observed (10-4) for Y chromosome markers in mammals . As for the mitochondrial d-loop, the size of the amplified fragment slightly limits the use of the technique. Some degraded samples do not amplify; even so, a particularly degraded sperm whale sample was still amplifiable (data not shown).
Since the Y chromosome population is expected to have a small effective size, it is more likely to be affected by genetic drift. Thus, it reflects more recent demographic events such as bottlenecks, expansions or founder effects . To study this sort of event, one needs a marker whose diversity is high enough to allow the reconstruction of gene genealogies with the least ambiguities and in regions where recombination does not interfere with the uniqueness of the trees. For this purpose, highly variable microsatellites represent valuable markers but they require intensive computing methods to take into account uncertainties in the trees arising from alleles that are identical by state and not by descent (homoplasies) [29, 30]. Adding a new sequence marker is therefore of interest for Y-chromosome population genetics in Cetartiodactyla. Moreover, the Bayesian estimate of mutation rate on each edge of both trees in Fig. 4, jointly computed with phylogenetic inference, shows high values for a marker of nuclear DNA: between 10-8 and 10-10 substitutions per site and per year in all Cetartiodactyla branches. This value is intermediate between those of mitochondrial d-loop and nuclear DNA in mammals [31, 32].
Functional perspectives in amelogenin evolution
We found two stop codons at amino acid positions 98 and 99 of exon 5 in all Y chromosome copies of amelogenin in the four studied cetacean species (positions 988–993 of sequence EMBL:AM744959). The Amel-Y gene product may therefore be truncated in these species or represent a pseudogene as already observed in species from most of the other eutherian clades 
The 460 bp insertion studied represents a single-event synapomorphy among most Cetartiodactyla. Together with the presence of other numerous indels informative at the order-level, it could help resolve the phylogenic discrepancies between hippopotamuses, pigs, camels and other Cetartiodactyla observed by many authors [7–9]. In addition, we demonstrate higher diversity within a single sequence than has yet been observed in multi-sequence assays . This high diversity should allow the use of this sequence as the male counterpart of the mitochondrial control region. The applications would include inference of male-driven evolution in population genetics as applied to breeding management; domestication studies in archaeogenetics ; conservation biology (population history, sex-biased dispersal, admixture); or for testing sex-biased selection . Amelogenin intron 4 amplification will also be an efficient tool for sexing Ruminantia and Cetacea. This will be useful for many fields of veterinary and forensic science (embryo technologies, in vitro fertilization, meat products). Finally, amelogenin amplification could also be a helpful tool for conservation biology through sampling of dead animals, faecal remains and biopsies of free-ranging animals like whales and dolphins. Amelogenin amplification in Cetartiodactyla therefore is a simple, single-step procedure with a wide range of applications.
Biological material was isolated from soft tissues sampled from dead stranded cetaceans and extracted using a classical phenol-chloroform protocol . We used heterologous primers, X5 (5'-GTGCTTACCCCTTTGAAGTG-3') and X6 (5'-CTTCCTCCCGCTTGGTCTTG-3'), designed from the amelogenin intron 4 and exon 5 of Homo sapiens X chromosome (GenBank:NC 000023, chrX:11221454–11228802) reference assembly Build 36.3, to amplify the homologous region in cetacean amelogenin. This region of Amel-X is 92% identical to Amel-Y.
We subsequently cloned and sequenced these PCR fragments and designed oligonucleotide primers specific to the Artiodactyla. They are anchored in exons 4 and 5 of Amel-X and Amel-Y, which allows complete amplification of the 4th intron (SC1: 5'-CAAGCATGCATTTCAATTCCC-3' and SC2: 5'-CTGCATGGGGAACATCGGAG-3'). Optimal PCR conditions were adjusted by using temperature (49–62°C) and MgCl2 (1–2.5 mM) gradients. Following this optimization step, the PCR amplifications were conducted in a reaction mix consisting of 477 mM KCl, 1 mM Tris/HCl pH 8.3, 1.5 mM MgC12, 250 μM each dNTP, 2 pmol/μl of each primer, 3 units of Taq polymerase in a final volume of 25 μl. Cycling was conducted as follows: 95°C for 2 min, followed by 30 cycles of denaturation at 94°C for 1 min, annealing at 55°C for 45 s, extension at 72°C for 1 min, and a final extension at 70°C for 10 min. PCR products were run on 1.2% agarose gel ethidium bromide stained, alongside a 1 kbp ladder (New England Biolabs, County Road, MA). In order to validate the assay, the gender, if identified during examination of the stranded carcasses, was recorded (22 out of 38 samples). These amplifications were conducted in eight Cetacean species: 20 Striped dolphin (Stenella cœruleoalba), five Fin whales (Balaenoptera physalus), four Bottlenosed dolphins (Tursiops truncatus), three Common dolphins (Delphinus delphis), three Gray whales (Eschrichtius robustus), two Sperm whales (Physeter macrocephalus), one Minke whale (Balænoptera acutorostrata) and one Risso's dolphin (Grampus griseus).
Of these, we sequenced striped dolphins (9 males and one female), fin whales (two males and two females), common dolphin (2 males), gray whale (one female) and Risso's dolphin (one female). Sequencing was performed on an ABI Prism sequencer (Applied Biosystems, Foster City, CA) with the dye terminator protocol directly for fragments amplified using the same (SC1-SC2) primer pair or after cloning in pGEM-T (Promega, Madison, WI) for fragments amplified using the X5-X6 pair. Each sample was sequenced from at least two independent amplifications for greater reliability.
The obtained sequences were aligned using ClustalW, a clustering-based sequence alignment algorithm , with bottlenosed dolphin (GenBank:AY787743S2), cow (GenBank:AB091789–AB091790), pig (GenBank:AF328419 and AB091792), sheep (GenBank:AY604731) and human (GenBank:NT_011757 from 9098117 to 9098612 and GenBank:NC_000024 from 6796200 to 6796719) sequences. Striped dolphin haplotypes were named according to membership in one of the two Mediterranean populations inferred by population genetic data on mitochondrial and microsatellite DNA (data not shown). Stenella cœruleoalba YAn thus represents the nth haplotype from population 1 while Stenella cœruleoalba YBn represents the nth haplotype from population 2.
Phylogenic inference was performed using a Bayesian Monte Carlo Markov Chain (MCMC) approach, implemented in Mrbayes , that allows determining posterior probability of each node of the tree.
Given the high indel density in this region, data were divided into two partitions: a sequence matrix (1120 characters corresponding to nucleotides) and a binary character matrix (presence/absence of the indels). Gap coding was achieved using Modified Complex Indel Coding  in which parsimony-informative gap characters were scored from pairwise unambiguously aligned sequences. Hence, a total of 47 gap characters were included (summarized in Figure 5b). We used a General Time Reversible (GTR) model, allowing for a proportion of invariant sites for sequence partition, and a binary F81-like model with a variable ascertainment bias for the insert's absence/presence (state 0/state 1) for the binary partition. In this model, π0 and π1 are the stationary probabilities of the two states and π0/π1 the rate of transition from state 1 to 0.
Many thanks go to H. Etchevers, N. Lambert and P. Monget for their much valuable comments on the manuscript as well as to P. Sudour and J. Butterworth.
We also wish to thank B. Jacobsen from the UAF Museum and F. Dhermain from the GECEM who kindly provided a significant part of the samples as well as A.-M. Gasc for technical help.
We are indebted to three anonymous reviewers for valuable comments that significantly helped to improve manuscript.
- Lahn BT, Page DC: Four evolutionary strata on the human X chromosome. Science. 1999, 286 (5441): 964-967. 10.1126/science.286.5441.964.View ArticlePubMedGoogle Scholar
- Lahn BT, Pearson NM, Jegalian K: The human Y chromosome, in the light of evolution. Nat Rev Genet. 2001, 2 (3): 207-216. 10.1038/35056058.View ArticlePubMedGoogle Scholar
- Satchell PG, Anderton X, Ryu OH, Luan XH, Ortega AJ, Opamen R, Berman BJ, Witherspoon DE, Gutmann JL, Yamane A, et al: Conservation and variation in enamel protein distribution during vertebrate tooth development. J Exp Zool (Mol Dev Evol). 2002, 294 (2): 91-106. 10.1002/jez.10148.View ArticleGoogle Scholar
- Delgado S, Girondot M, Sire JY: Molecular evolution of amelogenin in mammals. J Mol Evol. 2005, 60 (1): 12-30. 10.1007/s00239-003-0070-8.View ArticlePubMedGoogle Scholar
- Boisserie JR, Lihoreau F, Brunet M: The position of hippopotamidae within cetartiodactyla. Proc Natl Acad Sci USA. 2005, 102 (5): 1537-1541. 10.1073/pnas.0409518102.PubMed CentralView ArticlePubMedGoogle Scholar
- Gatesy J, O' Leary M: Deciphering whale origins with molecules and fossils. Trends Ecol Evol. 2001, 16 (10): 562-570. 10.1016/S0169-5347(01)02236-4.View ArticleGoogle Scholar
- Nikaido M, Rooney AP, Okada N: Phylogenetic relationships among cetartiodactyls based on insertions of short and long interpersed elements: hippopotamuses are the closest extant relatives of whales. Proc Natl Acad Sci USA. 1999, 96 (18): 10261-10266. 10.1073/pnas.96.18.10261.PubMed CentralView ArticlePubMedGoogle Scholar
- Price SA, Bininda Emonds ORP, Gittleman AL: A complete phylogeny of the whales, dolphins and even-toed hoofed mammals (Cetartiodactyla). Biol Rev. 2005, 80 (3): 445-473. 10.1017/S1464793105006743.View ArticlePubMedGoogle Scholar
- Thewissen JG, Williams EM, Roe LJ, Hussain ST: Skeletons of terrestrial cetaceans and the relationship of whales to artiodactyls. Nature. 2001, 413 (6853): 277-281. 10.1038/35095005.View ArticlePubMedGoogle Scholar
- Leduc RG, Perrin WF, Dizon AE: Phylogenetic relationships among the delphinid cetaceans based on full cytochrome B sequences. Mar Mamm Sci. 1999, 15 (3): 619-648. 10.1111/j.1748-7692.1999.tb00833.x.View ArticleGoogle Scholar
- Milinkovitch MC: Molecular phylogeny of cetaceans prompts revision of morphological transformations. Trends Ecol Evol. 1995, 10 (8): 328-334. 10.1016/S0169-5347(00)89120-X.View ArticlePubMedGoogle Scholar
- Hellborg L, Ellegren H: Low levels of nucleotide diversity in mammalian Y chromosomes. Mol Biol Evol. 2004, 21 (1): 158-163. 10.1093/molbev/msh008.View ArticlePubMedGoogle Scholar
- Meadows JR, Hanotte O, Drogemuller C, Calvo J, Godfrey R, Coltman D, Maddox JF, Marzanov N, Kantanen J, Kijas JW: Globally dispersed Y chromosomal haplotypes in wild and domestic sheep. Anim Genet. 2006, 37 (5): 444-453. 10.1111/j.1365-2052.2006.01496.x.View ArticlePubMedGoogle Scholar
- Anderung C, Hellborg L, Seddon J, Hanotte O, Gotherstrom A: Investigation of X- and Y-specific single nucleotide polymorphisms in taurine (Bos taurus) and indicine (Bos indicus) cattle. Anim Genet. 2007, 38 (6): 595-600.View ArticlePubMedGoogle Scholar
- Edwards CJ, Gaillard C, Bradley DG, MacHugh DE: Y-specific microsatellite polymorphisms in a range of bovid species. Anim Genet. 2000, 31 (2): 127-130. 10.1046/j.1365-2052.2000.00602.x.View ArticlePubMedGoogle Scholar
- Brinkmann B: Is the amelogenin sex test valid?. Int J Legal Med. 2002, 116 (2): 63-10.1007/s00414-001-0263-x.View ArticlePubMedGoogle Scholar
- Bryja J, Konečný A: Fast sex identification in wild mammals using PCR amplification of the Sry gene. Folia Zool. 2003, 52 (3): 269-274. [http://www.ivb.cz/folia/52/3/269-274.pdf]Google Scholar
- Shaw CN, Wilson PJ, White BN: A reliable molecular method of gender determination for mammals. J Mammal. 2003, 84 (1): 123-128. 10.1644/1545-1542(2003)084<0123:ARMMOG>2.0.CO;2.View ArticleGoogle Scholar
- Pomp D, Good BA, Geisert RD, Corbin CJ, Conley AJ: Sex identification in mammals with polymerase chain reaction and its use to examine sex effects on diameter of day-10 or -11 pig embryos. J Anim Sci. 1995, 73 (5): 1408-1415.PubMedGoogle Scholar
- Ennis S, Gallagher TF: A PCR-based sex-determination assay in cattle based on the bovine amelogenin locus. Anim Genet. 1994, 25 (6): 425-427.View ArticlePubMedGoogle Scholar
- Pfeiffer I, Brenig B: X- and Y-chromosome specific variants of the amelogenin gene allow sex determination in sheep (Ovis aries) and European red deer (Cervus elaphus). BMC Genet. 2005, 6 (1): 16-10.1186/1471-2156-6-16.PubMed CentralView ArticlePubMedGoogle Scholar
- Weikard R, Pitra C, Kuhn C: Amelogenin cross-amplification in the family Bovidae and its application for sex determination. Mol Reprod Dev. 2006, 73 (10): 1333-1337. 10.1002/mrd.20554.View ArticlePubMedGoogle Scholar
- Robertson BC, Gemmell NJ: PCR-based sexing in conservation biology: Wrong answers from an accurate methodology?. Conserv Genet. 2006, 7 (2): 267-271. 10.1007/s10592-005-9105-6.View ArticleGoogle Scholar
- Findlay I, Quirke P: Fluorescent polymerase chain reaction: Part I. A new method allowing genetic diagnosis and DNA fingerprinting of single cells. Hum Reprod Update. 1996, 2 (2): 137-152. 10.1093/humupd/2.2.137.View ArticlePubMedGoogle Scholar
- Walsh PS, Erlich HA, Higuchi R: Preferential PCR amplification of alleles: mechanisms and solutions. PCR Methods Appl. 1992, 1 (4): 241-250.View ArticlePubMedGoogle Scholar
- Findlay I, Ray P, Quirke P, Rutherford A, Lilford R: Allelic drop-out and preferential amplification in single cells and human blastomeres: implications for preimplantation diagnosis of sex and cystic fibrosis. Hum Reprod. 1995, 10 (6): 1609-1618.View ArticlePubMedGoogle Scholar
- Bininda-Emonds OR, Cardillo M, Jones KE, MacPhee RD, Beck RM, Grenyer R, Price SA, Vos RA, Gittleman JL, Purvis A: The delayed rise of present-day mammals. Nature. 2007, 446 (7135): 507-512. 10.1038/nature05634.View ArticlePubMedGoogle Scholar
- Petit E, Balloux F, Excoffier L: Mammalian population genetics: why not Y?. Trends Ecol Evol. 2002, 17 (1): 28-33. 10.1016/S0169-5347(01)02356-4.View ArticleGoogle Scholar
- Beerli P, Felsenstein J: Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach. Proc Natl Acad Sci USA. 2001, 98 (8): 4563-4568. 10.1073/pnas.081068098.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilson IJ, Balding DJ: Genealogical inference from microsatellite data. Genetics. 1998, 150 (1): 499-510.PubMed CentralPubMedGoogle Scholar
- Ayala FJ: Vagaries of the molecular clock. Proc Natl Acad Sci USA. 1997, 94 (15): 7776-7783. 10.1073/pnas.94.15.7776.PubMed CentralView ArticlePubMedGoogle Scholar
- Moritz C, Dowling TE, Brown WM: Evolution of Animal Mitochondrial DNA: Relevance for Population Biology and Systematics. Annu Rev Ecol Syst. 1987, 18: 269-292. 10.1146/annurev.es.18.110187.001413.View ArticleGoogle Scholar
- Sire JY, Delgado S, Fromentin D, Girondot M: Amelogenin: lessons from evolution. Arch Oral Biol. 2005, 50 (2): 205-212. 10.1016/j.archoralbio.2004.09.004.View ArticlePubMedGoogle Scholar
- Hatch LT, Dopman EB, Harrison RG: Phylogenetic relationships among the baleen whales based on maternally and paternally inherited characters. Mol Phylogenet Evol. 2006, 41 (1): 12-27. 10.1016/j.ympev.2006.05.023.View ArticlePubMedGoogle Scholar
- Zeder MA, Emshwiller E, Smith BD, Bradley DG: Documenting domestication: the intersection of genetics and archaeology. Trends Genet. 2006, 22 (3): 139-155. 10.1016/j.tig.2006.01.007.View ArticlePubMedGoogle Scholar
- Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 2001, Cold Spring Harbor: Cold Spring Harbor Laboratory Press, 3Google Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Müller K: Incorporating information from length-mutational events into phylogenetic analysis. Mol Phylogenet Evol. 2006, 38 (3): 667-676. 10.1016/j.ympev.2005.07.011.View ArticlePubMedGoogle Scholar
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19 (18): 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Nei M: Molecular evolutionary genetics. 1987, United States: Columbia University Press, New York, NYGoogle Scholar
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.