- Research article
- Open Access
Phylogenetic relationships between Hapalemur species and subspecies based on mitochondrial DNA sequences
BMC Evolutionary Biology volume 2, Article number: 4 (2002)
Phylogenetic relationships of the genus Hapalemur remains controversial, particularly within the Hapalemur griseus species group. In order to obtain more information on the taxonomic status within this genus, and particularly in the cytogenetic distinct subspecies group of Hapalemur griseus, 357 bp sequence of cytochrome b and 438 bp of 12S mitochondrial DNAs were analyzed on a sample of animals captured in areas extending from the north to the south-east of Madagascar. This sample covers all cytogenetically defined types recognized of the genus Hapalemur.
Phylogenetic trees and distances analyses demonstrate a first emergence of Hapalemur simus followed by H. aureus which is the sister clade of the H. griseus subspecies. Hapalemur griseus is composed of 4 subspecies separated into two clades. The first contains H. g. griseus, H. g. alaotrensis and H. g. occidentalis. The second consists of H. g. meridionalis. A new chromosomal polymorphic variant from the region of Ranomafana, H. griseus ssp, has been analysed and was found in both clades.
Our results support the raising of H. g. meridionalis to the specific rank H. meridionalis, while neither cytogenetic nor molecular evidences support the raising of H. g. alaotrensis to a species rank despite its morphological characteristics. The new cytotype H. g. ssp which has been previously characterized by cytogenetic studies contains animals clustering either with the group of Hapalemur griseus griseus or with that of Hapalemur meridionalis. This suggests the existence of an ancestral polymorphism or an introgression of mitochondrial DNA between subspecies.
The genus Hapalemur of the Malagasy Prosimian family Lemuridae comprises three clearly different species: Hapalemur griseus, H. simus and H. aureus. Classification of different types within the H. griseus group has been controversial, particularly in the light of the variations existing in coat colour . Systematic cytogenetic studies led to the description of four subspecies of H. griseus characterised by different karyotypes: H. g. griseus (HGG) (2N = 54) [5, 6], initially considered as the most widespread and supposed to live in the Eastern forest of Madagascar from the north to the south ; H. g. occidentalis (HGO) (2N = 58)  restricted in the north-west and the eastern part of Madagascar as east as Maroantsetra [4, 8]; H. g. meridionalis (HGM) (2N = 54)  found in the south, near Fort Dauphin, and differing from H. g. griseus in one metacentric chromosome ; H. g. alaotrensis (HGA) (2N = 54) , restricted to the bed of the Alaotra lake, with a karyotype similar to HGG containing more juxtacentromeric heterochromatin, and characterized by a larger body size . In addition, a new variant of HG, H. griseus ssp (HGssp), has recently been described. This type is characterized by a chromosomal polymorphism (2N = 54, 55 and 56) and is living in the area of Ranomafana-Kianjavato .
The comparison of the karyotypes of the different Hapalemur led to a phylogenetic tree characterised by an earlier emergence of H. simus (HSI), followed by H. aureus (HAU) which represents the sister clade of the different groups of the H. griseus subspecies [8, 10]. As molecular studies were fruitfully used for taxonomic investigations of several taxa, including lemurs [11–15], and as no large morphological differences could be found between the cytogenetically close subspecies of H. griseus, mitochondrial cytochrome b and 12S DNA analyses were performed in order to determinate more accurately their species and/or subspecies status.
A total of 115 animals covering all the genus Hapalemur were captured in bamboo forests extending from the north to the south of Madagascar. The different capture areas of each species and subspecies as the number of animals are indicated in Fig. 1.
The PCR amplification of the cytochrome b resulted in a 357 bp fragment. Alignments using blast search demonstrated that this fragment corresponds to the beginning of the mitochondrial cytochrome b gene. All the sequences were aligned and the substitutions were scored. Animals sharing the same sequences were regrouped into a same haplotype. These different haplotypes are summarised in Fig. 1 and table 2. Consensus sequences constructed with the most common sequence in each species and subspecies (not shown) demonstrated that transition substitutions occur more frequently than transversions (an average of eight transitions for 1 transversion) and gaps are absent. As expected , most of the substitutions occur on the third position of the codon (54 out of 64 substitutions), while the first position changed nine times and the second once. The portion amplified comprises the four invariant codons considered necessary for cytochrome b function . Changes were sometimes noticed in the third base of the codon in these positions, but the amino acid remains the same so that the functionality of the cytochrome b is not affected.
Kimura two parameter distances  and absolute distances calculated with the different haplotypes (Table 1) as phylogenetic trees constructed with the computer program Pylogenetic Analysis Using Parsimony (PAUP) 4.0b.4a  (Figs 2, 3, 4, 5) were used to compare the different Hapalemur species and subspecies. Neighbor-joining analyses of cytochrome b haplotypes demonstrated that H. simus emerged first followed by H. aureus which is the sister clade of the H. griseus subspecies. Inside H. griseus subspecies, two subclades are present and the H. griseus ssp (HGssp) are distributed in both the subclades (Figs 2, 3). For this reason, the H. griseus ssp were separated in HGsspa and Hgsspb. The first subclade (S1) contains HGA, HGG, HGsspa and HGO. The second subclade (S2) is composed of HGM and HGsspb. Maximum parsimony analyses (data not shown) demonstrated the same topology, unless the bootstrap support for S1 dropped from 94 to 66.
The use of the Kimura two parameter method (Table 1) shows inside S1 a very close distance between HGA and HGG (ranging from 0.00281 to 0.00847), while a larger distance occurs between HGG and HGsspa (0.00847 to 0.02875) as well as between HGG and HGO (0.01412 to 0.03180). In S2 (Table 1), short distances are observed between HGM and HGsspb (0–0.02855), (in this subclade, the sequence of HGM 02 and HGssp 01 are identical). Larger distances ranging from 0.10473 (HGsspa/HGM) to 0.06240 (HGM-HGsspb/HGG-HGsspa) separates S1 from S2 (Table 1).
In order to confirm the position of H. griseus ssp inside the H. griseus subclades, a portion of 12S mitochondrial DNA was amplified and sequenced. The resulted 438 bp fragments were aligned, and the sequences were analyzed with the program PAUP. The trees obtained with HGO, HGG, HGssp and HGM 12S sequences (Figs 4, 5) also demonstrated a separation of HGssp into two subclades. The first (HGsspa), which is composed of animals clustering with HGO and HGG when cytochrome b sequences are analyzed (Figs 2, 3), also clusters with HGO and HGG when 12S sequences are employed (Figs 4, 5). The second (HGsspb), constructed with animals clustering with HGM when cytochrome b sequences are analyzed (Figs 2, 3), also clusters with HGM when 12S sequences are employed (Figs 4, 5).
During the last decade, comparisons of mitochondrial DNA sequences have been very useful for the analyses of the phylogenetic species relationships including lemurs [11–15, 20]. In our study, cytochrome b and 12S mtDNA sequences were used, in order to clarify the position of each species and subspecies in the genus Hapalemur. Despite the fact that we analysed only short sequences, the genetic distances are considered to be relevant because of the large number of animals involved, since large sample size reduces errors in the estimation of evolutionary relatedness [11, 21].
Comparisons of cytochrome b sequences demonstrated that pairwise genetic distances between H. simus and H. aureus and the other Hapalemur (distances ranging from 0.09098 to 0.15338) are the highest. The corresponding phylogram also shows that these taxa are well separated. Our molecular data thus strongly support the species status for H. simus, H. aureus and H. griseus.
Among the H. griseus, the genetic divergence observed between H. g. meridionalis and the HGG, HGA, HGsspa and HGO (ranging from 0.06240 to 0.10473) is markedly higher than all other values of the same level of intra H. griseus comparison (ranging from 0.00281 to 0.03781). The phylogram also clearly separates the H. g. meridionalis from other H. griseus. This high genetic divergence in the ranging of those observed between species of lemurs, would support the classification of H. g. meridionalis in a separate species "H. meridionalis" despite the relatively small chromosomal differences observed .
Lower distances (ranging from 0.01133 to 0.03781) are observed between H. griseus and H. g. occidentalis. Moreover, these two forms cluster in the same clade. So, despite the cytogenetic differences existing between these two forms, which differ by two chromosomal rearrangements , we propose the maintenance of the subspecies status for H. g. occidentalis.
On the basis of our analyses, H. g. alaotrensis is indistinguishable from H. g. griseus. The very short genetic distances found between H. g. griseus and H. g. alaotrensis (0.00281 to 0.00847) suggests a combination of these two groups into a single subspecies. However, the differences of the body sizes  and the differences in the content of heterochromatin found in both karyotypes  support the separation of H. g. griseus and H. g. alaotrensis in separate subspecies , but in no case in species apart as it has been previously proposed .
The systematic position of HGssp is more difficult to establish when cytogenetic and molecular data are compared. Cytogenetic data demonstrated the existence of an unique chromosomal polymorphic subspecies characterized by two karyotypes, 2N = 56 and 2N = 54, and their hybrids 2N = 55 . In the view of molecular data, HGssp contains two groups, HGsspa and HGsspb, separated by an important genetic distance (0.06240–0.09719) in the range of those observed between species. Each of these two groups contains both karyotypes, 2N = 56 and 2N = 54, as well as their hybrids. The HGsspa are separated from HGG, HGA and HGO by a genetic distance in the range of those observed between subspecies (0.00563–0.03781), while HGsspb appears similar to HGM. Moreover, the HGsspb haplotypes are mixed with those of HGM. Taking into account these results, HGsspa should be considered as belonging to the group of HGG, HGA and HGO, and HGsspb as belonging to the group of HGM.
The existence of two well-separated clades among the Hapalemur griseus ssp originating from the Ranomafana region could possibly resulted from either recent mitochondrial DNA introgression or ancestral polymorphism. As the boundaries between HGM and HGssp are still unknown, we could hypothesize that the limits could be close to Ranomafana, allowing hybridization between these two forms. A transfer of mitochrondrial DNA from HGM into the HGssp population could have occurred through a matrilinear process resulting in the HGsspb haplotype. New investigations in areas located in the south of Ranomafana should thus allow the finding of HGM populations containing introgressed HGssp mitochondrial DNAs, unless this transfer occurs only from HGM to HGssp. A second hypothesis is that an ancestral polymorphic population containing both HGM and HGssp haplotypes and a chromosomal polymorphism have diverged in two separated populations. In the population of Ranomafana, the HGM and HGssp haplotypes as well as the chromosomal polymorphism were maintained. In the second population, only the HGM haplotype remained present, and a gain of a large block of heterochromatin gave rise to the karyotype characteristic of the HGM .
The comparisons of the phylogenetic trees based on mitochondrial DNA sequences with those previously obtained from cytogenetic data are only partially concordant. In both trees, H. simus emerges first, followed by H. aureus and then by the different H. griseus forms. The cytogenetic data allow to propose an evolutionary tree in which H. g. occidentalis emerges first, followed by H. g. griseus, and H. g. meridionalis, whereas on the tree based on cytochrome b and 12S sequences, H. g. meridionalis appears as a sister clade of the other H. griseus. This difference may be related to the short distances observed inside each clade which allowed no branching order.
Our molecular studies of the Hapalemur genus raise the question of the classification of H. g. meridionalis in the species status H. meridionalis. They also confirm the subspecies status of H. g. occidentalis and the absence of arguments in favour of the classification of H. g. alaotrensis as a separate species. The sequencing of the Hapalemur griseus ssp originating from Ranomafana reveals animals clustering either with the H. meridionalis or with the group of Hapalemur griseus griseus. No monophyletic clade could be determined in this new cytotype, so that the taxonomic status of the Hapalemur griseus ssp remains undefined. As our molecular data did not match the branching sequence within the H. griseus group based on cytogenetic data, further investigations including nuclear DNA will be necessary in order to resolve this issue.
Materials and methods
Six survey were organised from 1997 to 2001 allowing the capture of 115 animals in bamboo forests extending from the north to the south of Madagascar. Animals were captured using blowpipe projections and then sexed, weighted and measured. Skin samples were cut off under general anaesthesia with a 2 mg/kg injection of ketamine solution (Ketalar® Parke-Davis) and conserved deep frozen in liquid nitrogen. The different capture areas and the number of animals are indicated in Figure 1. As the H. griseus taxonomy is essentially based on cytogenetic criteria, karyotypes were made in order to confirm the species and subspecies rank. From the north to the south, the following subspecies were caught: a total of nine H. g. occidentalis (four at Analamera (a), two at Ambato (b), two at Ambakoany (c), and one at Maroantsetra (d)). Three H. g. alaotrensis at the Alaotra lake (e); 18 H. g. griseus at Maromiza (f); 50 H. griseus ssp. (37 at Ranomafana (h), five at Ambolomavo (i), eight at Kianjavato (j)); 16 H. g. meridionalis (14 in Andohahela (k) and two in Mandena (l)). In addition, 12 H. simus (four from Ranomafana (h), and eight from Ambolomavo (i)), as well as two H. simus and five H. aureus from the Zoological Park of Tsimbazaza (g) were studied. Animals were released at their capture location, immediately after recovery from anaesthesia, except two H. g. occidentalis captured in Ambato as well as three H. g. alaotrensis which are kept in the Zoological Park of Mulhouse and one H. g. occidentalis from Maroan setra which is kept in the private Zoological Park of Mandraka. As outgroup, we used a sample of Eulemur macaco flavifrons, from one of the individuals captured in the Sahamalaza forest for an earlier study .
DNA samples were extracted from the skin biopsies using the standard proteinase K digestion followed by a phenol chloroform extraction as described by Sambrook et al. with minor modifications. Small pieces of skin (~9mm2) were suspended in 200μl of extraction buffer (Tris 0.2M pH 8.4; KCl 0.5M; proteinase K: 5 mg/ml), and incubated at 37°C overnight. The samples were than heated at 95°C and mixed with an equal volume (200μl) of phenol/chloroform (1/1). After centrifugation (6 mn at 8500 g) the aqueous part was precipitated with 2 volumes of absolute ethanol at -20°C in presence of 1/10 (V/V) of ammonium acetate 5M. After centrifugation (5 mn at 8500 g), pellets were rinsed with 200μl of 70% ethanol and dried at 37°C. Pellets were resuspended in sterile double-distilled water and the concentration of the DNAs were measured by absorption at 260 nm. DNA samples were then stored at -30°C.
The polymerase chain reaction (PCR) was employed to generate a double-stranded fragment of 357 bp corresponding to a part of the mitochondrial cytochrome b gene. Each amplification was performed in the presence of 20 mM Tris-HCl pH 8.4, 50 mM KCl, 1.5mM MgCl2 (Gibco BRL), 200 μg BSA, 0.5 μM of each primer, 5U of Taq DNA polymerase (Perkin Elmer Cetus), 200 μM of each dNTP (Boehringer Mannheim) and 360ng of template DNA in a volume of 200 μl. The following primers derived from those described by Kocher et al. were employed: Pr181 (5'-CCATCCAACATGTCAGCATGATGAAA-3') and Pr182 (5'-CCCTCAGAATGATATTTGTCCTCA-3'). Reactions were done in a Perkin Elmer Cetus DNA thermocycler 480 as follows: predenaturation (7 mn at 93 °C) and 35 cycles of denaturation (30 s at 94°C), annealing (45 s at 42°C), extension (1 mn 30 s at 72°C), followed by a final extension step (10 mn at 72°C).
12S mitochondrial DNA
The same PCR conditions than those described for the cytochrome b amplifications were applied, with a concentration of 2 mM MgCl2 and the following set of primers: Pr179 (5'-AAACTAGGATTAGATACCCTATTAT-3') and Pr180 (5'-AAGAGCGACGGGCGATGTGT-3'). Amplifications were performed in a Perkin Elmer Cetus DNA thermocycler 480 as follows: predenaturation (10 mn at 94°C) and 40 cycles of denaturation (30 s at 94°C), annealing (45 s at 45°C), extension (1 mn 30 s at 72°C), followed by a final extension step (10 mn at 72°C).
Amplification products were electrophoresed on 1.3% agarose gels in TBE buffer (Tris base 87 mM, boric acid 89 mM, EDTA 2 mM, pH 8.0) in the presence of 0.5 μg/ml of ethidium bromide. Electrophoresis were performed at 200 mA for 2 h in TBE. Gels were than examined and photographed under UV light with a Polaroid system. The major bands on the gels were cut and centrifuged (8500 g 10 mn) in order to recover the DNA fragments . The fragments were then precipitated with 0.1 volumes of ammonium acetate 5 M and 2 volumes of absolute ethanol at -20°C. After centrifugation (8500 g 10-mn), the pellets were air dried, resuspended in 50 μl of sterile bidistilled water and sequenced on an automatic ABI PRISM sequencer with the Taq dye deoxy terminator cycle sequencing kit. Each sample was sequenced from 5' to 3' and 3' to 5'.
Cytochrome b sequences were aligned and for each taxa, similar sequences were grouped under one haplotype. These haplotypes and the EMBL GenBank accession numbers are listed in table 2. EMBL GenBank accession numbers for the 12S sequences are listed in table 3.
The aligned sequences were analysed using neighbor-joining and maximum parsimony methods with the computer program Pylogenetic Analysis Using Parsimony (PAUP) 4.0b.4a . Genetic distances (d) measured in neighbor-joining analyses were calculated using the Kimura two parameter method  with the following formula: d=-1/21n [(1-2P-Q)Xsqr(1-2Q)], with P=transitions/positions scored and Q=transversions/positions scored. For neighbor-joining and maximum parsimony analyses, bootstraps of 10000 replicates were performed to examine the relative support of each relationship in the resultant topologies. Maximum parsimony trees were calculated via fast stepwise addition with random addition sequence.
Link JA:Lemur griseus.Beytr Naturg. 1785, 65.
Gray JE: Notes on Hapalemur (Prolemur) simus, a new species lately living in the gardens of the society. Proc Zool Soc Lond. 1870, 828-831.
Meier B, Albignac R, Peyrieras A, Rumpler Y, Wright P: A new species of Hapalemur (primates) from south east Madagascar. Folia Primatol. 1987, 48: 211-215.
Petter JJ, Albignac R, Rumpler Y: Mammifères lémuriens (primates prosimiens). In: Faune de Madagascar. Edited by: ORSTOM, CNRS. 1977, 44:
Chu EHY, Swomley BA: Chromosomes of lemurine lemurs. Science. 1961, 133: 1925-1926.
Rumpler Y, Albignac R: Cytogenetic study of the endemic Malagasy lemur: Hapalemur, I. Geoffroy, 1851. J Hum Evol. 1973, 2: 267-270.
Tattersall I: The primates of Madagascar. New York, Columbia University Press. 1982
Rumpler Y, Prosper , Hauwy M, Rabarivola C, Rakotoarisoa G, Dutrillaux B: Chromosomal evolution of the Hapalemur griseus subspecies (Malagasy Prosimian), including a new chromosomal polymorphic variant. Chromosome Research (accepted).
Warter S, Randrianasolo G, Dutrillaux B, Rumpler Y: Cytogenetic study of a new subspecies of Hapalemur griseus. Folia Primatol. 1987, 48: 50-55.
Rumpler Y, Warter S, Ishak B: Chromosomal evolution in prosimians. Hum Evol. 1989, 4: 157-170.
Montagnon D, Ravaoarimanana B, Rumpler Y: Taxonomic relationships and sampling effects among Lepilemuridae and Lemuridae using partial cytochrome b gene. CR Acad Sci Paris/Life Sci Serie III. 2001, 334 (7): 647-656. 10.1016/S0764-4469(01)01331-2.
Montagnon D, Ravaoarimanana B, Rakotosamimana B, Rumpler Y: Ancient DNA from Megaladapis edwardsi (Malagasy Subfossil): Preliminary results using partial cytochrome b sequence. Folia Primatol. 2001, 72: 30-32. 10.1159/000049916.
Pastorini J, Forstner MRJ, Martin RD: Relationships among Brown lemurs (Eulemur fulvus) based on mitochondrial DNA sequences. Mol Phylogenet Evol. 2000, 16 (3): 418-429. 10.1006/mpev.2000.0796.
Pastorini J, Martin RD, Ehresmann P, Zimmerman E, Forstner MRJ: Molecular phylogeny of the lemur family Cheirogaleidae (primates) based on mitochondrial DNA sequences. Mol Phylogenet Evol. 2001, 19 (1): 45-56. 10.1006/mpev.2000.0904.
Wyner Y, DeSalle R, Absher R: Phylogeny and character behaviour in the family lemuridae. Mol Phylogenet Evol. 2000, 15 (1): 124-134. 10.1006/mpev.1999.0723.
Yoder AD, Vilgalys R, Ruvolo M: Molecular evolutionary dynamics of cytochrome b in strepsirrhine primates: the phylogenetic significance of third-position transvertion. Mol Biol Evol. 1996, 13 (10): 1339-1350.
Howell N, Gilbert K: Mutational analysis of the mouse mitochondrial cytochrome b gene. J Mol Biol. 1988, 203: 607-618.
Kimura MA: Simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 11-120.
Swofford D: PAUP*4.0b.4a. Phylogenetic analysis using parsimony* (and other methods),. Sunderland, MA, Sinauer Associates. 2000
Hayasaka K, Fujii K, Horai S: Molecular phylogeny of macaques: Implication of nucleotide sequences from a 896-base pair region of mitochondrial DNA. Mol Biol Evol. 1996, 13 (7): 1044-1053.
Martin AP, Bailey D, Kessing D, Palumbi SR: Accuracy of estimating genetic distances between species from short sequences of mitochondrial DNA. Mol Biol Evol. 1990, 7 (5): 485-488.
Rumpler Y, Dutrillaux B: Chromosomal evolution in Malagasy lemurs. III. Chromosome banding studies in the genus Hapalemur and the species Lemur catta. Cytogenet Cell Genet. 1978, 21: 201-211.
Rumpler Y: The significance of chromosomal studies in the systematics of the Malagasy lemurs. In: Lemur Biology . Edited by: Tattersall I, Sussman RW. 1975, New York, Plenum Press, 25-40.
Groves CP: A theory of human and primate evolution. Oxford, Clarendon press. 1989
Fausser JL, Rabarivola C, Meir B, Hahn T, Rumpler Y: Genetic comparison between different populations of Eulemur macaco flavifrons in Northwest Madagascar using RAPD markers. Am J Primatol. 2000, 51: 249-255. 10.1002/1098-2345(200008)51:4<249::AID-AJP4>3.0.CO;2-4.
Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. Cold Spring Harbor, Cold Spring Harbor press. 1989, 2
Kocher TD, Thomas WK, Meyer A, Edwards SV, Pääbo S, Villablanca FX, Wilson AC: Dynamics of mitochondrial DNA evolution in animals: Amplification and sequencing with conserved primers. Proc Natl Acad Sci USA. 1989, 86: 6196-6200.
Tautz D, Rentz M: An optimized free-squeeze method for the recovery of DNA fragments from agarose gels. Anal Bio chem. 1983, 132: 14-19.
This work was partially supported by the Association Européenne pour l'Etude et la Conservation des Lémuriens (A.E.E.C.L.), the DREIF-MENESR (EA3428), and the Institut de Génétique et de Biologie Moléculaire et Cellulaire (Pr. P. Chambon). The authors thank the Commission Tripartite of the Malagasy Government and the Ministère pour la Production Animale et des Eaux et Forêts for their permission to capture the animals, to take samples and for delivery CITES. Thanks are also due to Dr J. Ganzhorn for the correction of the manuscript, Mrs M. Lavaux for her secretarial assistance, Dr B. Ravaoarimanana and M. Hauwy for their helpful comments, and to S. Vicaire (IGBMC, Strasbourg) for sequencing the cytochrome b and 12S fragments. Special thanks for insight comments on the cladograms are due to Dr. D. Montagnon.