Mitochondrial DNA diversity and PCR-based sex determination of Irrawaddy dolphin (<Emphasis Type="Italic">Orcaella</Emphasis><Emphasis Type="Italic">brevirostris</Emphasis>) from Chilika Lagoon,India

Of the only known two Lagoon populations of Irrawaddy dolphins (Orcaella) in the world, one is residing in the Chilika Lagoon in Orissa state, India. In addition to accidental deaths in gill net fishery and mechanized boat operations, there has been exploitation of the species for their oil. Extreme patchy distribution and vulnerability to becoming entangled in fishing gear has made it a focus of conservation concern. Information on genetic diversity of populations has considerable potential for informing conservation plans. The present paper reports the first genetic study of O. brevirostris from Chilika Lagoon based on mtDNA sequencing and PCR-based sex identification from 11 individuals. Control region sequence comparison showed two haplotypes and cytochrome b a single haplotype in the Chilika population of the species. Phylogenetic analysis indicated distinct clades within the Asian samples, with the Indian population showing closest genetic proximity to the haplotypes from Thailand. Sex of the animal was determined by PCR-based method. It is important to continue to examine the population discreteness and genetic variation of Irrawaddy dolphin in Chilika Lagoon vis-à-vis its global geographic distribution for formulating the conservation plans of the species.     

Geographic differentiation of <Emphasis Type="Italic">Cobitis shikokuensis</Emphasis> inferred from mtDNA RFLP analysis

Mitochondrial DNA divergence among populations of the Japanese spinous loach Cobitis shikokuensis, endemic to Shikoku Island, was investigated by restricted fragment length polymorphism analysis. A total of 68 restriction sites on DNA fragments from the cytochrome b to D-loop regions and from the 12S rRNA to 16S rRNA regions, amplified by PCR, were analyzed. A total of 12 haplotypes (plus 6 in outgroups) were detected in 268 specimens collected from 19 localities in seven rivers (and 41 specimens from four localities in three rivers in outgroups). Three of the seven river populations of C. shikokuensis were shown to have unique haplotypes, and four of the seven river populations were monomorphic. The nested structure of the haplotype network for populations of C. shikokuensis exhibited two large clades corresponding to (1) populations from the Shimanto River and its neighbors and (2) two genetically divergent populations in the Shigenobu and Iwamatsu Rivers. The population from the Shimanto River, the largest river inhabited by C. shikokuensis, maintains great haplotype diversity as well as the allozyme diversity previously reported. On the other hand, populations from the Hiji River, the second largest river, which exhibited the highest allozyme diversity, were monomorphic in their mtDNA. The nested clade analysis (NCA) revealed that past fragmentation between the above two clades could occur in the initial distribution process of C. shikokuensis. The large genetic divergence of two river populations from the Shigenobu and Iwamatsu Rivers was inferred to be caused by a process of long distance colonization and fragmentation. MtDNA introgression into the Hiji River population from southern river populations was suggested. Taking genetic divergence into consideration, each river population of C. shikokuensis should be conserved separately as like a distinct species, and conservation programs for the small populations showing less genetic variability should be invoked as soon as possible.     

Diversification and subspecies patterning of the goitered gazelle (Gazella subgutturosa) in Iran

Goitered gazelles, Gazella subgutturosa, exist in arid and semiarid regions of Asia from the Middle to the Far East. Although large populations were present over a vast area until recently, a decline of the population as a result of hunting, poaching, and habitat loss led to the IUCN classification of G. subgutturosa as “vulnerable." We examined genetic diversity, structure, and phylogeny of G. subgutturosa using mitochondrial cytochrome b sequences from 18 geographically distant populations in Iran. The median‐joining network of cyt b haplotypes indicated that three clades of goitered gazelles can be distinguished: a Middle Eastern clade west of the Zagros Mountains (and connected to populations in Turkey and Iraq), a Central Iranian clade (with connection to Azerbaijan), and an Asiatic clade in northeastern Iran (with connection to Turkmenistan, Uzbekistan, and other Asian countries as far as northeastern China and Mongolia). Based on our results, we argue that Iran is the center of diversification of goitered gazelles, due to the presence of large mountain ranges and deserts that lead to the separation of populations. In accordance with previous morphological studies, we identified the Asiatic clade as the subspecies G. s. yarkandensis, and the other two clades as the nominate form G. s. subgutturosa. The new genetic information for goitered gazelles in Iran provides the basis for future national conservation programs of this species.     

Molecular diversity and phylogeography of Andricus curtisii (Hymenoptera,Cynipidae) in Turkey

A 433-bp fragment of the mitochondrial cytochrome b gene and the entire nuclear ITS2 region were sequenced in an oak gall wasp species, Andricus curtisii, which represents 21 populations across its range in Turkey. Forty cytochrome b haplotypes and 63 ITS2 alleles were discovered among the 161 individuals studied. Remarkable genetic variation was detected in the species, with an average haplotype diversity (h) of 0.72 and 0.76 and an average nucleotide diversity (π) of 0.06 and 0.01 for the cyt b gene and the ITS region, respectively. Most populations were genetically distinct, possessing unique cytochrome b haplotypes. Maximum likelihood and Bayesian phylogenetic analyses and application of the haplotype network to both datasets revealed that the most basally located haplotypes/alleles were from the eastern part of Turkey. The detection of older sequences in the phylogenetic trees of the easterly located populations implied an eastern origin of the species. The application of a conventional insect mitochondrial DNA clock to the dataset suggested a splitting of ingroup haplotypes from the outgroup lineages predating the Pleistocene epoch and the formation of two major haplogroups that coincided with the oscillation of the early Pleistocene glaciations.     

Genetic diversity and conservation units in wild and captive populations of endangered freshwater fishes: a case of <Emphasis Type="Italic">Hemigrammocypris rasborella</Emphasis> in Shizuoka,Japan

Population structure and genetic diversity were examined using partial mitochondrial cytochrome b gene sequences of four wild, one reintroduced, and five captive populations of the endangered cyprinid Hemigrammocypris rasborella from three river systems in the easternmost region of the species’ range in Shizuoka Prefecture, central Honshu, Japan. We detected loss of genetic diversity from portions of the wild and captive populations, as well as suspected nonindigenous haplotypes in some captive, reintroduced, and even wild populations. Given the population structure revealed, we suggest that the populations should be managed with consideration for both the endemism and viability (avoidance of inbreeding depression) of the local populations.     

Gorges partition diversity within New Zealand flathead Galaxias populations

Understanding the landscape factors governing population connectivity in riverine ecosystems represents an ongoing challenge for freshwater biologists. We used DNA sequence analysis to test the hypothesis that major geomorphological features underpin freshwater-limited fish diversity in a tectonically dynamic region of New Zealand. Phylogeographic analysis of 101 Galaxias depressiceps cytochrome b sequences, incorporating 55 localities from southern New Zealand, revealed 26 haplotypes, with only one shared among rivers. We detect strong hierarchical genetic differentiation both among and within river systems. Genetic structuring is particularly pronounced across the Taieri River system (63 individuals from 35 sites, 18 haplotypes), with 92% of variation partitioned among locations. Distinctive within-river genetic clusters are invariably associated with major subcatchment units, typically isolated by substantial gorges. The anomalous distribution of a single lineage across a major drainage divide is consistent with local, tectonically driven headwater capture. We conclude that major landscape features such as gorges can strongly partition riverine fish diversity and constrain freshwater biodiversity.     

Intraspecific molecular phylogeny and phylogeography of the Meriones meridianus (Rodentia: Cricetidae) complex in northern China reflect the processes of desertification and the Tianshan Mountains uplift

The phylogeographical patterns and demographic history of mitochondrial DNA (cytochrome b, N = 327; D‐loop, N = 252) and nuclear DNA (IRBP gene, N = 235) haplotypes were studied for the Meriones meridianus complex in northern China, a desert‐dwelling gerbil species complex. The phylogenetic analyses, which were performed on the separate and combined (mitochondrial + nuclear) datasets, revealed two divergent clades (Clade A and Clade B) corresponding to distinct geographical regions. Clade A contained the haplotypes found mostly in individuals from the Tianshan Mountains area. Clade B contained haplotypes from populations located in other deserts in northern China. The divergence times indicated that the history of the M. meridianus complex was influenced by the uplift of the Tianshan Mountains and climate‐induced habitat fluctuations. In the Pleistocene, the expansion of forests and grasslands during interglacial period led to the isolation of the M. meridianus complex, which preferred to inhabit deserts. Hence, long geological isolation and the M. meridianus complex adaptation to local ecological conditions led to its genetic divergence. Clade A had long‐lasting demographic stability, most likely because the populations of this clade remained in a stable desert environment for a long time. However, the extension of other deserts and disappearance of palaeolakes during the last glacial period resulted in demographic expansion of Clade B. Furthermore, our genetic data indicated that two subspecies may exist within the M. meridianus complex. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 362–383.     

Low diversity and high levels of population genetic structuring in introduced eastern mosquitofish (<Emphasis Type="Italic">Gambusia holbrooki</Emphasis>) in the greater Melbourne area,Australia

Eastern mosquitofish (Gambusia holbrooki) were introduced into Australia in 1925 and released to control mosquitoes. Gambusia holbrooki rapidly became invasive in recipient environments and now threaten native fauna. In this study, we used five polymorphic microsatellite loci and sequence from two mitochondrial genes, cytochrome b and cytochrome oxidase I, to evaluate genetic variation, colonisation and movement patterns of introduced G. holbrooki in the greater Melbourne area, and to assist in identifying the feasibility of local eradication. Microsatellite variation was consistently low within populations and there was evidence of bottleneck events for several populations. Populations displayed significant structuring associated with river basins rather than geographic distance, suggesting that habitat connectivity is important for dispersal. However, a few populations within river basins were more closely related to populations in other river basins than within their own basin, most likely reflecting a role of human-assisted dispersal in population establishment. Mitochondrial sequencing revealed only a single haplotype and suggested all populations were founded by individuals from a common source. These genetic data help delineate boundaries for local management strategies.     

Genetic Variability in the Yellow-Necked Field Mouse (<Emphasis Type="Italic">Sylvaemus flavicollis</Emphasis> Melch., 1834, Muridae,Rodentia) at the Eastern Border of the Range

The first characterization of the genetic variability in populations of S. flavicollis on the eastern border of the range is presented. Seven individuals from the most northeastern habitat of the species (Middle Urals) were karotyped. No deviations from the standard chromosome set, either by the chromosome number or morphology, were revealed. Analysis of the mtDNA cytochrome b gene sequences (1133 bp) in 44 individuals from five populations on the eastern border of the species range (Middle and Southern Urals) resulted in identification of 17 haplotypes. All haplotypes were new and not found earlier in other parts of the species range. The genetic diversity indices and analysis of the demographic and genetic structure indicate a relatively recent origin of the populations under study as a result of rapid expansion. Phylogenetic analysis (97 haplotypes, including the GenBank data) showed that all haplotypes described at the eastern border of the range belonged to the same phylogroup distributed in the Balkan region, Northern and Eastern Europe, and Russia (Samara oblast). Close relationships between the examined populations and the populations from the northern part of the species range in Eastern Europe was demonstrated.     

Mitochondrial DNA divergence in the critically imperilled pygmy madtom,Noturus stanauli (Siluriformes: Ictaluridae)

The mitochondrial cytochrome b gene was sequenced for six individuals of the pygmy madtom, Noturus stanauli, a globally imperilled catfish, from both known localities in Tennessee (U.S.A.) separated by over 1055 river km. Phylogenetic and population genetic analyses revealed little divergence between these populations for this locus.     

基于ISSR分子标记的褐马鸡亲缘关系分析

褐马鸡(Crossoptilon mantchuricum)是中国特有濒危鸟类,国家一级保护动物。为了保护褐马鸡种质资源,保障驯养繁育种群和再引入种群的遗传基因结构优化,采用ISSR分子标记技术,对山西庞泉沟国家自然保护区和太原动物园两个种群35个褐马鸡进行了亲缘关系分析。从20条ISSR引物中筛选出10条扩增条带清晰稳定、重复性好的引物,共扩增出65条DNA条带,其中77%呈多态性。两个种群各个体之间的Nei's无偏差遗传距离和遗传相似度分析表明:个体间遗传相似性的平均值分别为0.5061与0.7591,遗传距离的平均值分别为0.4939与0.2409。用组平均法对35个褐马鸡个体进行聚类分析,结果显示:大多数同一种群的个体首先聚在一起,显示出生态地理条件相似的种群亲缘关系较近。     

Mitochondrial DNA variation in natural populations of endangered Indian Feather-Back Fish, Chitala chitala

Genetic variation at mitochondrial cytochrome b (cyt b) and D-loop region reveals the evidence of population sub-structuring in Indian populations of highly endangered primitive feather-back fish Chitala chitala. Samples collected through commercial catches from eight riverine populations from different geographical locations of India were analyzed for cyt b region (307 bp) and D-loop region (636–716 bp). The sequences of the both the mitochondrial regions revealed high haplotype diversity and low nucleotide diversity. The patterns of genetic diversity, haplotypes networks clearly indicated two distinct mitochondrial lineages and mismatch distribution strongly suggest a historical influence on the genetic structure of C. chitala populations. The baseline information on genetic variation and the evidence of population sub-structuring generated from this study would be useful for planning effective strategies for conservation and rehabilitation of this highly endangered species.     

Biogeography and molecular phylogeny of the genus Schizothorax (Teleostei: Cyprinidae) in China inferred from cytochrome b sequences

Aim To test a vicariant speciation hypothesis derived from geological evidence of large‐scale changes in drainage patterns in the late Miocene that affected the drainages in the south‐eastern Tibetan Plateau. Location The Tibetan Plateau and adjacent areas. Methods The cytochrome b DNA sequences of 30 species of the genus Schizothorax from nine different river systems were analysed. These DNA sequences were analysed using parsimony, maximum likelihood and Bayesian methods. The approximately unbiased and Shimodaira–Hasegawa tests were applied to evaluate the statistical significance of the shortest trees relative to alternative hypotheses. Dates of divergences between lineages were estimated using the nonparametric rate smoothing method, and confidence intervals of dates were obtained by parametric bootstrapping. Results The phylogenetic relationships recovered from molecular data were inconsistent with traditional taxonomy, but apparently reflected geographical associations with rivers. Within the genus Schizothorax, we observed a divergence between the lineages from the Irrawaddy–Lhuit and Tsangpo–Parlung rivers, and tentatively dated this vicariant event back to the late Miocene (7.3–6.8 Ma). We also observed approximately simultaneous geographical splits within drainages of the south‐eastern Tibetan Plateau, the Irrawaddy, the Yangtze and the Mekong–Salween rivers in the late Miocene (7.1–6.2 Ma). Main conclusions Our molecular evidence tentatively highlights the importance of palaeoriver connections and the uplift of the Tibetan Plateau in understanding the evolution of the genus Schizothorax. Molecular estimates of divergence times allowed us to date these vicariant scenarios back to the late Miocene, which agrees with geological suggestions for the separation of these drainages caused by tectonic uplift in south‐eastern Tibet. Our results indicated the substantial role of vicariant‐based speciation in shaping the current distribution pattern of the genus Schizothorax.     

祁连山区特有物种黑虎耳草的居群遗传多样性研究

利用叶绿体基因(rbcL和trnS-G)及核糖体DNA内转录间隔区(ITS),对祁连山区狭域分布的特有物种黑虎耳草8个居群(115个个体)进行遗传多样性研究,旨在揭示黑虎耳草的居群遗传结构及其历史进化过程。结果表明:(1)所有取样个体共检测到4个cpDNA单倍型和9个ITS单倍型,其中祁连山东南部的居群固定较多的单倍型和特有单倍型,而西北部居群只固定少数几个广泛分布的单倍型,且遗传多样性普遍较低。(2)基于cpDNA数据和ITS数据的分子变异分析(AMOVA)显示,遗传变异主要来源于居群内。(3)基于cpDNA数据的中性检验表明,Tajima’s D(-1.012 30,P> 0.05)和Fu&Li’s D*(-2.066 77,P> 0.05)均为负值,均不显著;歧点分布分析结果显示,黑虎耳草居群经历过近期的扩张事件。根据物种现有遗传分布格局推测,黑虎耳草在第四纪冰期时退缩到祁连山东南部的边缘避难所,间冰期或冰期后回迁到祁连山西北部地区,在回迁过程中由于奠基者效应导致祁连山西北部的居群仅固定少数广泛分布的单倍型,并呈现出较低的遗传多样性;由于居群较小且相互隔离,该...     

Geoclimatic factors influence the population genetic connectivity of Incarvillea arguta (Bignoniaceae) in the Himalaya–Hengduan Mountains biodiversity hotspot

Geoclimatic factors related to the uplift of the Himalaya and the Quaternary climatic oscillations influence the population genetic connectivity in the Himalaya–Hengduan Mountains (HHM) biodiversity hotspot. Therefore, to explore the relative roles played by these two factors, we examined the population dynamics and dispersal corridors of Incarvillea arguta (Royle) Royle incorporating ensemble species distribution modelling (SDM). Thirty‐seven populations were genotyped using plastid chloroplast DNA and low copy nuclear gene (ncpGS) sequences. Phylogeographic analysis was carried out to reveal the genetic structure and lineage differentiation. Ensemble SDMs were carried out for distributional change in the last glacial maximum, present, and future. Finally, the least cost path method was used to trace out possible dispersal corridors. The haplotypes were divided into four clades with strong geographical structure. The late Miocene origin of I. arguta in the western Himalaya ca. 7.92 Ma indicates lineage diversification related to the uplift of the HHM. The variability in habitat connectivity revealed by SDM is due to change in suitability since the Pleistocene. A putative dispersal corridor was detected along the drainage systems and river valleys, with strong support in the eastern Hengduan Mountains group. Our results support the signature of geoclimatic influence on population genetic connectivity of I. arguta in the HHM. We proposed that the major drainage systems might have assisted the rapid dispersal of isolated riverine plant species I. arguta in the HHM. The population genetic connectivity, using the fine‐tuned ensemble SDMs, enables scientists and policymakers to develop conservation strategies for the species gene pool in the HHM biodiversity hotspots.     

Phylogenetic Relationships within GenusLeuciscus(Pisces, Cyprinidae) in Portuguese Fresh Waters, Based on Mitochondrial DNA CytochromebSequences

To investigate phylogenetic relationships amongLeuciscusspecies occurring in Portuguese inland waters, the cytochromebgene was sequenced from representatives of the main rivers. This study supports the recognition of the species level forL. pyrenaicus,including populations from the southern Portuguese drainages (Tejo, Sado, and Guadiana drainages), and forL. carolitertii,including populations from the northern Portuguese drainages. The existence of two new species occurring in the extreme southwestern drainages of Mira and Arade is also suggested. The present results support the monophyly of the Mira and the Arade populations, as well as an early divergence of these two lineages. The present-day distribution ofLeuciscusspecies is seen as a consequence of Pliocene and Pleistocene events, such as river disjunctions and posterior confluence in epicontinental seas and river captures. A mixture of haplotypes was observed in the Mondego and the Tejo drainages, which could be a consequence of ancient river captures, with a possible mitochondrial DNA introgression in the Tejo drainage and a recent introduction by man in the Mondego drainage. The pattern of differentiation among mtDNA haplotypes and their geographic distribution is discussed in terms of evolutionary aspects.     

Genetic bottleneck in invasive species: the potato tuber moth adds to the list

The level of genetic diversity within populations of introduced species has received increasing attention as an important factor influencing their survival and adaptive potential. We examined this issue with the Guatemalan potato tuber moth Tecia solanivora, an agricultural pest which has successfully invaded South America and the Canary Islands within the last 20 years. To analyse changes in T. solanivora genetic diversity, the mitochondrial marker cytochrome b was sequenced from individuals collected across its known distribution area. High haplotypic diversity was observed in Guatemala, whereas only three haplotypes have been found in Venezuela and a single one in the remaining invaded South American countries and the Canary Islands. Invasive haplotypes were not observed in our samples from Guatemala but are closely related to Guatemalan haplotypes. These results are consistent with the hypotheses that (i) either a few individuals were introduced into Venezuela leading to a strong initial genetic bottleneck, or selection pressure may have lead to the disparition of all but a reduced number of introduced haplotypes, (ii) a second bottleneck occured between Venezuela and Colombia, and (iii) the invasion of the Canary Islands originated from South America. We further reviewed the recent literature to compare this change in genetic diversity with those reported for other invasive species. We quantified the changes in genetic diversity between native and introduced ranges for 57 biological invasions. We found that the genetic homogenization in T. solanivora was among the strongest reported and discuss factors that can explain the success of invasive populations with low genetic diversity.     

The role of waterholes as 'refugia' in sustaining genetic diversity and variation of two freshwater species in dryland river systems (Western Queensland, Australia)

1. Episodic floods and extended low or no flow periods characterise dryland river systems in Western Queensland, Australia. During protracted intervals between floods, rivers consist of a series of isolated waterholes, which serve as ‘refugia’ for aquatic species and much of the channel is dry. We categorised these waterholes into ‘main waterholes’, which are located in the main part of the river channel and ‘satellite waterholes’, which are located in distributary river channels. 2. We used mitochondrial sequences and allozymes to investigate levels of genetic diversity and patterns of connectivity among waterholes for two obligate freshwater species: Macrobrachium australiense (Decapoda: Palaemonidae) and Notopala sublineata (Gastropoda: Viviparidae). 3. We sampled 31 waterholes for M. australiense and 12 for N. sublineata. Based on a 505‐bp fragment of cytochrome oxidase subunit I, we identified 54 haplotypes in a sample of 232 individuals for M. australiense and based on a 457‐bp fragment of the same gene, 36 haplotypes in a sample of 145 individuals for N. sublineata. 4. Both nuclear and mitochondrial genetic data sets indicated that estimates of genetic diversity were not different in populations inhabiting main and satellite waterholes for either species. Also, there was generally very limited genetic differentiation among populations at any site. 5. We suggest that levels of connectivity among populations inhabiting waterholes at most sites are higher than expected. High levels of connectivity may help to maintain overall high levels of genetic diversity as well as low levels of genetic differentiation among waterholes within sites.     

Mitochondrial DNA analysis reveals hidden genetic diversity in captive populations of the threatened American crocodile (Crocodylus acutus) in Colombia

Identification of units within species worthy of separate management consideration is an important area within conservation. Mitochondrial DNA (mtDNA) surveys can potentially contribute to this by identifying phylogenetic and population structure below the species level. The American crocodile (Crocodylus acutus) is broadly distributed throughout the Neotropics. Its numbers have been reduced severely with the species threatened throughout much of its distribution. In Colombia, the release of individuals from commercial captive populations has emerged as a possible conservation strategy that could contribute to species recovery. However, no studies have addressed levels of genetic differentiation or diversity within C. acutus in Colombia, thus complicating conservation and management decisions. Here, sequence variation was studied in mtDNA cytochrome b and cytochrome oxidase I gene sequences in three Colombian captive populations of C. acutus. Two distinct lineages were identified: C. acutus‐I, corresponding to haplotypes from Colombia and closely related Central American haplotypes; and C. acutus‐II, corresponding to all remaining haplotypes from Colombia. Comparison with findings from other studies indicates the presence of a single “northern” lineage (corresponding to C. acutus‐I) distributed from North America (southern Florida), through Central America and into northern South America. The absence of C. acutus‐II haplotypes from North and Central America indicates that the C. acutus‐II lineage probably represents a separate South American lineage. There appears to be sufficient divergence between lineages to suggest that they could represent two distinct evolutionary units. We suggest that this differentiation needs to be recognized for conservation purposes because it clearly contributes to the overall genetic diversity of the species. All Colombian captive populations included in this study contained a mixture of representatives of both lineages. As such, we recommend against the use of captive‐bred individuals for conservation strategies until further genetic information is available.     

Species boundaries in Malagasy snakes of the genus Madagascarophis (Serpentes: Colubridae sensu lato) assessed by nuclear and mitochondrial markers

We studied mitochondrial divergence in 27 individuals of colubrid snakes of the genus Madagascarophis Mertens from most of its distribution area in Madagascar. Combined analyses of 16S rRNA and cytochrome b sequences identified six major clades which only partly agreed with previously proposed classifications. Analysis of nuclear DNA sequences of the c-mos gene as well as of ISSR fingerprints revealed consistent differences only among three clades which we consider as distinct species: a widespread Madagascarophis colubrinus (Schlegel), with M. citrinus (Boettger) as a junior synonym, a southern M. meridionalis Domergue, and a presumably undescribed species from the extreme north of Madagascar. The species M. ocellatus Domergue was not available for our study. Within M. colubrinus there are two populations from the north-west, each showing two divergent haplotypes with pairwise divergences of up to 5.2% in the cytochrome b gene. Maximum divergence in this gene within M. colubrinus was 7.1%. These high values emphasise that caution needs to be applied before genetic distance values are used for species delimitation. Phylogeographically, most of the genetic variation in M. colubrinus is found in northern Madagascar, indicating that the species might have originated in this region. Later one haplotype clade colonised western and eastern Madagascar, with a putative secondary introgression into north-western populations.