Population genetics of the yellow fever mosquito in Trinidad: comparisons of amplified fragment length polymorphism (AFLP) and restriction fragment length polymorphism (RFLP) markers

Recent development of DNA markers provides powerful tools for population genetic analyses. Amplified fragment length polymorphism (AFLP) markers result from a polymerase chain reaction (PCR)-based DNA fingerprinting technique that can detect multiple restriction fragments in a single polyacrylamide gel, and thus are potentially useful for population genetic studies. Because AFLP markers have to be analysed as dominant loci in order to estimate population genetic diversity and genetic structure parameters, one must assume that dominant (amplified) alleles are identical in state, recessive (unamplified) alleles are identical in state, AFLP fragments segregate according to Mendelian expectations and that the genotypes of an AFLP locus are in Hardy-Weinberg equilibrium (HWE). The HWE assumption is untestable for natural populations using dominant markers. Restriction fragment length polymorphism (RFLP) markers segregate as codominant alleles, and can therefore be used to test the HWE assumption that is critical for analysing AFLP data. This study examined whether the dominant AFLP markers could provide accurate estimates of genetic variability for the Aedes aegypti mosquito populations of Trinidad, West Indies, by comparing genetic structure parameters using AFLP and RFLP markers. For AFLP markers, we tested a total of five primer combinations and scored 137 putative loci. For RFLP, we examined a total of eight mapped markers that provide a broad coverage of mosquito genome. The estimated average heterozygosity with AFLP markers was similar among the populations (0.39), and the observed average heterozygosity with RFLP markers varied from 0.44 to 0.58. The average FST (standardized among-population genetic variance) estimates were 0.033 for AFLP and 0.063 for RFLP markers. The genotypes at several RFLP loci were not in HWE, suggesting that the assumption critical for analysing AFLP data was invalid for some loci of the mosquito populations in Trinidad. Therefore, the results suggest that, compared with dominant molecular markers, codominant DNA markers provide better estimates of population genetic variability, and offer more statistical power for detecting population genetic structure.     

Evaluation of factors affecting individual assignment precision using microsatellite data from horse breeds and simulated breed crosses

Assignment tests have been utilized to investigate population classification, measure genetic diversity and to solve forensic questions. Using microsatellite data from 26 loci genotyped in eight horse breeds we examined how population differentiation, number of scored loci, number of scored animals per breed and loci variability affected individual assignment precision applying log likelihood methods. We found that both genetic differentiation and number of scored loci were highly important for recognizing the breed of origin. When comparing two and two breeds, a proportion of 95% of the most differentiated breeds (0.200 < or = FST < or = 0.259) could be identified scoring only three loci, while the corresponding number was six for the least differentiated breeds (0.080 < or = FST < or = 0.139). An identical proportion of simulated breed crosses, differentiated from their parental breeds by FST estimates in the range 0.050-0.069, was identified when scoring 12 loci. This level of source identification was not obtained for the less differentiated breed crosses. The current data further suggested that population sample size and locus variability were not critical for the assignment precision as long as moderately large sample sizes (> or = 20 animals per population) and fairly variable loci were used.     

Genetic structure of the forest herb Primula elatior in a changing landscape

To investigate whether changes in land use and associated forest patch turnover affected genetic diversity and structure of the forest herb Primula elatior, historical data on landscape changes were combined with a population genetic analysis using dominant amplified fragment length polymorphism markers. Based on nine topographic maps, landscape history was reconstructed and forest patches were assigned to two age classes: young (less than 35 years) and old (more than 35 years). The level of differentiation among Primula populations in recently established patches was compared with the level of differentiation among populations in older patches. Genetic diversity was independent of population size (P > 0.05). Most genetic variation was present within populations. Within-population diversity levels tended to be higher for populations located in older forests compared with those for populations located in young forests (Hj = 0.297 and 0.285, respectively). Total gene diversity was also higher for old than for young populations (Ht = 0.2987 and 0.2828, respectively). The global fixation index FST averaged over loci was low, but significant. Populations in older patches were significantly more differentiated from each other than were populations in recently established patches and they showed significant isolation by distance. In contrast, no significant correlations between pairwise geographical distance and FST were found for populations in recently established patches. The location of young and old populations in the studied system and altered gene flow because of increased population density and decreased inter-patch distances between extant populations may explain the observed lower genetic differentiation in the younger populations. This study exemplifies the importance of incorporating data on historical landscape changes in population genetic research at the landscape scale.     

The use of amplified fragment length polymorphism in determining species trees at fine taxonomic levels: analysis of a medically important snake, Trimeresurus albolabris

There are a huge number of phylogenetic studies based on mitochondrial DNA (mtDNA); however, these may represent gene trees that may not be congruent with the species tree. A solution to this problem is to include additional, independent, loci from the nuclear genome. At fine taxonomic levels, i.e. between populations and closely related species, previously suggested nuclear markers such as intron sequence data may not be appropriate. In this study we investigate the use of amplified fragment length polymorphisms (AFLP) to aid determination of the species tree for 24 specimens of a medically important snake, Trimeresurus albolabris. This is of particular importance for many venomous snakes as venom often varies intraspecifically. Five different primer combinations produced 434 bands and were analysed by constructing a phylogenetic tree using neighbour joining and principal component analysis. Results were similar across all methods and found distinct groupings. The results were compared with mtDNA data and a reconciled tree was constructed in order to determine the species tree for T. albolabris. We found that T. albolabris (sensu lato) is not monophyletic. Specimens from the Indonesian islands (except West Java) form a distinct clade and we propose elevation to species level. A specimen from Nepal is also distinct and suggests that this population also deserves specific status. We suggest that AFLPs may prove a valuable aid in determining species trees as opposed to gene trees at fine taxonomic levels and this should facilitate the incorporation of molecular data into such activities as antivenom production and conservation management.     

Genetic diversity in an endangered alpine plant, Eryngium alpinum L. (Apiaceae), inferred from amplified fragment length polymorphism markers

Eryngium alpinum L. is an endangered species found across the European Alps. In order to obtain base-line data for the conservation of this species, we investigated levels of genetic diversity within and among 14 populations from the French Alps. We used the amplified fragment length polymorphism (AFLP) technique with three primer pairs and scored a total of 62 unambiguous, polymorphic markers in 327 individuals. Because AFLP markers are dominant, within-population genetic structure (e.g. FIS) could not be assessed. Analyses based either on the assumption of random-mating or on complete selfing lead to very similar results. Diversity levels within populations were relatively high (mean Nei's expected heterozygosity = 0.198; mean Shannon index = 0.283), and a positive correlation was detected between both genetic diversity measurements and population size (Spearman rank correlation: P = 0. 005 and P = 0.002, respectively). Moreover, FST values and exact tests of differentiation revealed high differentiation among populations (mean pairwise FST = 0.40), which appeared to be independent of geographical distance (nonsignificant Mantel test). Founder events during postglacial colonizations and/or bottlenecks are proposed to explain this high but random genetic differentiation. By contrast, we detected a pattern of isolation by distance within populations and valleys. Predominant local gene flow by pollen or seed is probably responsible for this pattern. Concerning the management of E. alpinum, the high genetic differentiation leads us to recommend the conservation of a maximum number of populations. This study demonstrates that AFLP markers enable a quick and reliable assessment of intraspecific genetic variability in conservation genetics.     

A Bayesian approach to inferring population structure from dominant markers

Molecular markers derived from polymerase chain reaction (PCR) amplification of genomic DNA are an important part of the toolkit of evolutionary geneticists. Random amplified polymorphic DNA markers (RAPDs), amplified fragment length polymorphisms (AFLPs) and intersimple sequence repeat (ISSR) polymorphisms allow analysis of species for which previous DNA sequence information is lacking, but dominance makes it impossible to apply standard techniques to calculate F-statistics. We describe a Bayesian method that allows direct estimates of FST from dominant markers. In contrast to existing alternatives, we do not assume previous knowledge of the degree of within-population inbreeding. In particular, we do not assume that genotypes within populations are in Hardy-Weinberg proportions. Our estimate of FST incorporates uncertainty about the magnitude of within-population inbreeding. Simulations show that samples from even a relatively small number of loci and populations produce reliable estimates of FST. Moreover, some information about the degree of within-population inbreeding (FIS) is available from data sets with a large number of loci and populations. We illustrate the method with a reanalysis of RAPD data from 14 populations of a North American orchid, Platanthera leucophaea.     

Genetic diversity in Tetrachaetum elegans, a mitosporic aquatic fungus

Tetrachaetum elegans Ingold is a saprobic aquatic hyphomycete for which no sexual stage has yet been described. It occurs most commonly during the initial decay of tree leaves in temperate freshwater habitats and typically sporulates under water. Dispersal of the aquatic fungus takes place primarily in the water column and has a large passive component. Differences in substrate composition (e.g. quality of leaf litter) may also play a role in the distribution of different species or genotypes. The population genetic structure of T. elegans was studied using amplified fragment length polymorphism (AFLP) multilocus fingerprints. The populations were isolated from the leaf litter of three different tree genera, sampled in nine streams distributed throughout a mixed deciduous forest. Molecular markers were developed for 97 monosporic isolates using four selective primer pairs. A total of 247 fragments were scored, of which only 32 were polymorphic. Significant stream differentiation was detected for the isolates considered in this study. Analysis of molecular variance revealed that 20% of the genetic variation observed was the result of differences between streams. No correlation between genetic and geographical distances was found but a few multilocus genotypes were observed in different locations. Altogether these results suggest that environmental barriers play a role in the population structure of this aquatic fungus. No clear-cut effect of leaf litter composition on genetic variation could be demonstrated. Finally, tests of linkage disequilibrium between the 32 polymorphic AFLP loci as well as simulations did not provide a final answer regarding clonality in T. elegans. Indeed, it was possible to reject linkage equilibrium at different sampling levels and show that full linkage was unlikely.     

Genetic analysis reveals population structure and recent migration within the highly fragmented range of the Cross River gorilla (Gorilla gorilla diehli)

Recently developed methods of individual-based analysis of genetic data allow an unprecedented opportunity to understand the relationships among fragmented populations. By defining population structure and identifying migrant individuals, such analyses can provide a framework to aid in evaluating the threats posed by inbreeding and reduced genetic variability as a consequence of limited gene flow among fragments. Here we investigate population structure in the critically endangered Cross River gorilla (Gorilla gorilla diehli) by applying a suite of individual-based analyses to data obtained from between one-quarter and one-third of the estimated total population through the use of noninvasively collected DNA samples. The population structure inferred using data from 11 autosomal microsatellite loci was broadly consistent with geography and habitat fragmentation, but showed no simple isolation-by-distance effects. In contrast to previous field surveys, which suggested that all gorilla localities were isolated from one another, we infer low levels of gene flow and identify migrants between habitat fragments as well as individuals of admixed ancestry, suggesting persistent recent reproductive contact between many of the localities. These results are encouraging for the conservation of the Cross River gorilla population. Conservation efforts should strive to maintain connectivity between subpopulations that are still in migratory contact and attempt to restore connectivity where it has been lost.     

Do polymorphic loci require large sample sizes to estimate genetic distances?

The coefficient of variation of estimates of three genetic distances (standard genetic distance of Nei, chord distance, FST) was examined with computer simulation to determine if large samples (per population) are necessary to precisely estimate genetic distances at loci with high levels of polymorphism. These simulations showed that loci with high mutation rates produce estimates of genetic distance with lower coefficients of variation than loci with lower mutation rates--without requiring larger sample sizes from each population. In addition, the rate at which increasing sample sizes decreases the coefficient of variation of estimates of genetic distances was shown to be approximately determined by the value of FST between the populations being sampled. When FST was greater than 0.05, sampling fewer than 20 individuals (per population) should be sufficient. When FST was less than 0.01, sampling 100 individuals (per population) or more will be useful.     

Global Population Genetic Structure and Male-Mediated Gene Flow in the Green Turtle (Chelonia Mydas): RFLP Analyses of Anonymous Nuclear Loci

We introduce an approach for the analysis of Mendelian polymorphisms in nuclear DNA (nDNA), using restriction fragment patterns from anonymous single-copy regions amplified by the polymerase chain reaction, and apply this method to the elucidation of population structure and gene flow in the endangered green turtle, Chelonia mydas. Seven anonymous clones isolated from a total cell DNA library were sequenced to generate primers for the amplification of nDNA fragments. Nine individuals were screened for restriction site polymorphisms at these seven loci, using 40 endonucleases. Two loci were monomorphic, while the remainder exhibited a total of nine polymorphic restriction sites and three size variants (reflecting 600-base pair (bp) and 20-bp deletions and a 20-bp insertion). A total of 256 turtle specimens from 15 nesting populations worldwide were then scored for these polymorphisms. Genotypic proportions within populations were in accord with Hardy-Weinberg expectations. Strong linkage disequilibrium observed among polymorphic sites within loci enabled multisite haplotype assignments. Estimates of the standardized variance in haplotype frequency among global collections (FST = 0.17), within the Atlantic-Mediterranean (FST = 0.13), and within the Indian-Pacific (FST = 0.13), revealed a moderate degree of population substructure. Although a previous study concluded that nesting populations appear to be highly structured with respect to female (mitochondrial DNA) lineages, estimates of Nm based on nDNA data from this study indicate moderate rates of male-mediated gene flow. A positive relationship between genetic similarity and geographic proximity suggests historical connections and/or contemporary gene flow between particular rookery populations, likely via matings on overlapping feeding grounds, migration corridors or nonnatal rookeries.     

Character compatibility of molecular markers to distinguish asexual and sexual reproduction

Particularly in polyploids, the potential of the high variability of dominant markers such as random amplified polymorphic DNA fragments (RAPDs) and amplified fragment length polymorphisms (AFLPs) in population genetic studies and analysis of breeding systems is reduced due to their dominant nature. In contrast, the criterion of character compatibility is hindered neither by dominance nor by polyploidy as allelic interpretation is not necessary. Character compatibility, which can be used to detect events of genetic exchange (or recombination), is particularly informative if these events are expected to be rare such as in taxa with extensive vegetative reproduction or apomixis. Binary unordered characters such as presence and absence of anonymous DNA markers are incompatible if all four pairwise combinations of character states are present among the individuals studied. Because incompatible character state distributions defy any progenitor–derivative relationship among individuals, they provide strong evidence for genetic exchange. Both the absolute number of incompatible character combinations and the probability of compatibility can be used as a measure of incompatibility. Although these measures may not directly relate to the frequency of genetic exchange, they provide a useful tool to heuristically explore data sets. The most commonly used input for multivariate analyses and analysis of molecular variance in population genetic studies of (dis)similarity of marker distributions are amalgamates of mutation and recombination. Character compatibility can be used to complement these traditional methods of analysis. Advantages and disadvantages of character incompatibility relative to multilocus analysis of modes of reproduction and population genetics are demonstrated with data from RAPDs, isozymes, and restriction fragment length polymorphisms (RFLPs) of the nuclear ribosomal and chloroplast genome.     

Accurate gene diversity estimates from amplified fragment length polymorphism (AFLP) markers

Three procedures for the estimation of null allele frequencies and gene diversity from dominant multilocus data were empirically tested in natural populations of the outcrossing angiosperm Persoonia mollis (Proteaceae). The three procedures were the square root transform of the null homozygote frequency, the Lynch & Milligan procedure, and the Bayesian method. Genotypes for each of 116 polymorphic loci generated by amplified fragment length polymorphism (AFLP) were inferred from segregation patterns in progeny arrays. Therefore, for the plus phenotype (band present), heterozygotes were distinguished from homozygotes. In contrast to previous studies, all three procedures produced very similar mean estimates of heterozygosity, which were in turn accurate estimators of the direct value (HO = 0.28). A second population of P. mollis displayed markedly lower levels of heterozygosity (HO = 0.20) but approximately twice as many polymorphic loci (284). These AFLP results show that biases in estimates of average null allele frequency and heterozygosity are largely eliminated in highly polymorphic dominant marker data sets displaying a J-shaped beta distribution with a high percentage of loci containing more than three null homozygotes and relatively few loci with no null homozygotes. This distribution may be typical of outcrossing angiosperms.     

Population genetic structure of brown trout (Salmo trutta L.) within a northern boreal forest stream

Brown trout (Salmo trutta L.) population genetic structure and its temporal stability were studied within a small forest stream in central Sweden using five microsatellite loci. Both resident and migrating brown trout are present in the watershed. Tissue samples were collected from seven sections of the stream during two consecutive years. No differences were found in multilocus FST estimates between years within sections except in one case. Moreover, differences between age cohorts within sample sections were rare. The low interannual variation and the low heterogeny between cohorts is interpreted as indications of temporal stability. Pairwise multilocus FST estimates increased with increasing geographic distance, indicating isolation by distance. It is argued that the brown trout of the F?rs?n stream represent a population complex. The structure is probably maintained by precise natal homing and a limited, but important, amount of gene flow between closely situated sections within the stream supporting a stepping-stone model of gene flow.     

Genetic diversity and population structure of Tasmanian devils, the largest marsupial carnivore

Genetic diversity and population structure were investigated across the core range of Tasmanian devils (Sarcophilus laniarius; Dasyuridae), a wide-ranging marsupial carnivore restricted to the island of Tasmania. Heterozygosity (0.386-0.467) and allelic diversity (2.7-3.3) were low in all subpopulations and allelic size ranges were small and almost continuous, consistent with a founder effect. Island effects and repeated periods of low population density may also have contributed to the low variation. Within continuous habitat, gene flow appears extensive up to 50 km (high assignment rates to source or close neighbour populations; nonsignificant values of pairwise FST), in agreement with movement data. At larger scales (150-250 km), gene flow is reduced (significant pairwise FST) but there is no evidence for isolation by distance. The most substantial genetic structuring was observed for comparisons spanning unsuitable habitat, implying limited dispersal of devils between the well-connected, eastern populations and a smaller northwestern population. The genetic distinctiveness of the northwestern population was reflected in all analyses: unique alleles; multivariate analyses of gene frequency (multidimensional scaling, minimum spanning tree, nearest neighbour); high self-assignment (95%); two distinct populations for Tasmania were detected in isolation by distance and in Bayesian model-based clustering analyses. Marsupial carnivores appear to have stronger population subdivisions than their placental counterparts.     

Estimating species trees using multiple-allele DNA sequence data

Several techniques, such as concatenation and consensus methods, are available for combining data from multiple loci to produce a single statement of phylogenetic relationships. However, when multiple alleles are sampled from individual species, it becomes more challenging to estimate relationships at the level of species, either because concatenation becomes inappropriate due to conflicts among individual gene trees, or because the species from which multiple alleles have been sampled may not form monophyletic groups in the estimated tree. We propose a Bayesian hierarchical model to reconstruct species trees from multiple-allele, multilocus sequence data, building on a recently proposed method for estimating species trees from single allele multilocus data. A two-step Markov Chain Monte Carlo (MCMC) algorithm is adopted to estimate the posterior distribution of the species tree. The model is applied to estimate the posterior distribution of species trees for two multiple-allele datasets--yeast (Saccharomyces) and birds (Manacus-manakins). The estimates of the species trees using our method are consistent with those inferred from other methods and genetic markers, but in contrast to other species tree methods, it provides credible regions for the species tree. The Bayesian approach described here provides a powerful framework for statistical testing and integration of population genetics and phylogenetics.     

Genetic structure of the Iberian pig breed using microsatellites

An analysis of 25 microsatellite loci in 210 animals has been used to define the genetic structure of the Iberian pig, traditionally classified into several varieties. In addition, a sample of 20 Duroc pigs was used as an outgroup for topology trees. Inter-variety genetic variation was estimated by unbiased average heterozygosity and the number of alleles observed. Significant deviations from the Hardy-Weinberg equilibrium (HWE) were shown for 19 loci across the whole population. By contrast, equilibrium deviation within varieties was much lower. Genetic variation measures, genetic distance values and a neighbour-joining tree were used to estimate subdivision. In addition, an individual tree was constructed to contrast the assignation of animals into varieties. Despite the low bootstrap values obtained in the varieties neighbour-joining tree, the degree of genetic variation found was sufficient to support the division of the Iberian pig into varieties, although in some cases the traditional classification cannot be accepted. These results have shown the value of this marker panel in the study of intra-breed genetic structures.     

Measures of divergence between populations and the effect of forces that reduce variability

Wright's FST and related statistics are often used to measure the extent of divergence among populations of the same species relative to the net genetic diversity within the species. This paper compares several definitions of FST which are relevant to DNA sequence data, and shows that these must be used with care when estimating migration parameters. It is also pointed out that FST is strongly influenced by the level of within-population diversity. In situations where factors such as selection on closely linked sites are expected to have stronger effects on within-population diversity at some loci than at others, differences among loci can result entirely from differences in within- population diversities. It is shown that several published cases of differences in FST among regions of high and low recombination in Drosophila may be caused in this way. For the purpose of comparisons of levels of between-population differences among loci or species which are subject to different intensities of forces that reduce variability within local populations, absolute measures of divergence between populations should be used in preference to relative measures such as FST.      

Species trees from gene trees: reconstructing Bayesian posterior distributions of a species phylogeny using estimated gene tree distributions

The desire to infer the evolutionary history of a group of species should be more viable now that a considerable amount of multilocus molecular data is available. However, the current molecular phylogenetic paradigm still reconstructs gene trees to represent the species tree. Further, commonly used methods of combining data, such as the concatenation method, are known to be inconsistent in some circumstances. In this paper, we propose a Bayesian hierarchical model to estimate the phylogeny of a group of species using multiple estimated gene tree distributions, such as those that arise in a Bayesian analysis of DNA sequence data. Our model employs substitution models used in traditional phylogenetics but also uses coalescent theory to explain genealogical signals from species trees to gene trees and from gene trees to sequence data, thereby forming a complete stochastic model to estimate gene trees, species trees, ancestral population sizes, and species divergence times simultaneously. Our model is founded on the assumption that gene trees, even of unlinked loci, are correlated due to being derived from a single species tree and therefore should be estimated jointly. We apply the method to two multilocus data sets of DNA sequences. The estimates of the species tree topology and divergence times appear to be robust to the prior of the population size, whereas the estimates of effective population sizes are sensitive to the prior used in the analysis. These analyses also suggest that the model is superior to the concatenation method in fitting these data sets and thus provides a more realistic assessment of the variability in the distribution of the species tree that may have produced the molecular information at hand. Future improvements of our model and algorithm should include consideration of other factors that can cause discordance of gene trees and species trees, such as horizontal transfer or gene duplication.     

Genetic diversity of 15 autosomal short tandem repeats loci using the AmpFLSTR? Identifiler? kit in a Bhil Tribe Population from Gujarat state,India

MATERIALS AND METHODS:

The genetic diversity and forensic parameters based on 15 autosomal short tandem repeats (STR) loci; D8S1179,D21S11, D7S820, CSF1PO, D3S1358, TH01, D13S317,D16S539, D2S1338, D19S433, vWA, TPOX, D18S51,D5S818, and FGA in AmpFLSTR® Identifiler™ kit from Applied Biosystems, Foster City, CA, USA were evaluated in saliva samples of 297 unrelated individuals from the Bhil Tribe population of Gujarat state, India to study genetic diversities and relatedness of this population with other national and international populations.

RESULTS:

Statistical analysis of the data revealed all loci were within Hardy-Weinberg Equilibrium expectations with the exception of the locus vWA (0.019) and locus D18S51 (0.016). The neighbour joining phylogeny tree and Principal Co-ordinate Analysis plot constructed based on Fst distances from autosomal STRs allele frequencies of the present study and other national as well as international populations show clustering of all the South Asian populations in one branch of the tree, while Middle Eastern and African populations cluster in a separate branch.

CONCLUSION:

Our findings reveal strong genetic affinities seen between the Indo-European (IE) speaking Bhil Tribe of Gujarat and Dravidian groups of South India.     

Effects of Asymmetric Nuclear Introgression,Introgressive Mitochondrial Sweep,and Purifying Selection on Phylogenetic Reconstruction and Divergence Estimates in the Pacific Clade of Locustella Warblers

When isolated but reproductively compatible populations expand geographically and meet, simulations predict asymmetric introgression of neutral loci from a local to invading taxon. Genetic introgression may affect phylogenetic reconstruction by obscuring topology and divergence estimates. We combined phylogenetic analysis of sequences from one mtDNA and 12 nuDNA loci with analysis of gene flow among 5 species of Pacific Locustella warblers to test for presence of genetic introgression and its effects on tree topology and divergence estimates. Our data showed that nuDNA introgression was substantial and asymmetrical among all members of superspecies groups whereas mtDNA showed no introgression except a single species pair where the invader''s mtDNA was swept by mtDNA of the local species. This introgressive sweep of mtDNA had the opposite direction of the nuDNA introgression and resulted in the paraphyly of the local species'' mtDNA haplotypes with respect to those of the invader. Тhe multilocus nuDNA species tree resolved all inter- and intraspecific relationships despite substantial introgression. However, the node ages on the species tree may be underestimated as suggested by the differences in node age estimates based on non-introgressing mtDNA and introgressing nuDNA. In turn, the introgressive sweep and strong purifying selection appear to elongate internal branches in the mtDNA gene tree.