Extreme genetic differences among populations of Gazella granti, Grant's gazelle in Kenya

Estimating admixture pedigrees of recent hybrids without a contiguous reference genome

The genome of recently admixed individuals or hybrids has characteristic genetic patterns that can be used to learn about their recent admixture history. One of these are patterns of interancestry heterozygosity, which can be inferred from SNP data from either called genotypes or genotype likelihoods, without the need for information on genomic location. This makes them applicable to a wide range of data that are often used in evolutionary and conservation genomic studies, such as low-depth sequencing mapped to scaffolds and reduced representation sequencing. Here we implement maximum likelihood estimation of interancestry heterozygosity patterns using two complementary models. We furthermore develop apoh (Admixture Pedigrees of Hybrids), a software that uses estimates of paired ancestry proportions to detect recently admixed individuals or hybrids, and to suggest possible admixture pedigrees. It furthermore calculates several hybrid indices that make it easier to identify and rank possible admixture pedigrees that could give rise to the estimated patterns. We implemented apoh both as a command line tool and as a Graphical User Interface that allows the user to automatically and interactively explore, rank and visualize compatible recent admixture pedigrees, and calculate the different summary indices. We validate the performance of the method using admixed family trios from the 1000 Genomes Project. In addition, we show its applicability on identifying recent hybrids from RAD-seq data of Grant's gazelle (Nanger granti and Nanger petersii) and whole genome low-depth data of waterbuck (Kobus ellipsiprymnus) which shows complex admixture of up to four populations.

Delineating origins of cheetah cubs in the illegal wildlife trade: Improvements based on the use of hair δ18O measurements

All African felids are listed as vulnerable or endangered according to the IUCN (International Union for Conservation of Nature) Red List of Threatened Species. Cheetahs (Acinonyx jubatus) in particular have declined rapidly as a result of human impacts so that development of effective strategies and tools for conservation of this highly vulnerable species, as well as African felids in general, are essential for their survival in the wild. Here we use the oxygen stable isotopic compositions of cheetah hair to determine origins of cheetah cubs destined for the illegal exotic pet trade by associating individual cubs with predicted δ18O isoscape locations. We found that cheetah cubs most likely originated in East Africa, close to the corridors responsible for this aspect of the illegal wildlife trade to the Middle East. Further refinement of these assignments using a two isotope analysis (δ18O and δ13C values) indicate that these cubs were likely sourced in Southern Ethiopia or possibly as far as Tanzania. We also demonstrate that δ18O values in tissues can provide provenance information in cases where results of δ2H analyses may be obscured by the effects of metabolic routing of nutrients during nursing, starvation, or dehydration. This study demonstrates the utility of stable isotopic tools for conservation and forensic uses for endangered mammalian species.

Vicariance followed by secondary gene flow in a young gazelle species complex

Grant's gazelles have recently been proposed to be a species complex comprising three highly divergent mtDNA lineages (Nanger granti, N. notata and N. petersii). The three lineages have nonoverlapping distributions in East Africa, but without any obvious geographical divisions, making them an interesting model for studying the early‐stage evolutionary dynamics of allopatric speciation in detail. Here, we use genomic data obtained by restriction site‐associated (RAD) sequencing of 106 gazelle individuals to shed light on the evolutionary processes underlying Grant's gazelle divergence, to characterize their genetic structure and to assess the presence of gene flow between the main lineages in the species complex. We date the species divergence to 134,000 years ago, which is recent in evolutionary terms. We find population subdivision within N. granti, which coincides with the previously suggested two subspecies, N. g. granti and N. g. robertsii. Moreover, these two lineages seem to have hybridized in Masai Mara. Perhaps more surprisingly given their extreme genetic differentiation, N. granti and N. petersii also show signs of prolonged admixture in Mkomazi, which we identified as a hybrid population most likely founded by allopatric lineages coming into secondary contact. Despite the admixed composition of this population, elevated X chromosomal differentiation suggests that selection may be shaping the outcome of hybridization in this population. Our results therefore provide detailed insights into the processes of allopatric speciation and secondary contact in a recently radiated species complex.

An ancient hybridization event reconciles mito-nuclear discordance among spiral-horned antelopes

The spiral-horned antelopes (genus Tragelaphus) are among the most phenotypically diverse of all large mammals, and evolved in Africa during an adaptive radiation that began in the late Miocene, around 6 million years ago. Tragelaphus was able to exploit the habitat heterogeneity created by Plio-Pleistocene paleoclimatic fluctuations and tectonic processes to eventually occupy almost every habitat type in present day sub-Saharan Africa. The smallest of the spiral-horned antelopes, the bushbuck (T. scriptus), is also widely distributed across Africa, but is genetically divided into polyphyletic Scriptus and Sylvaticus mitochondrial (mt)DNA superlineages that inhabit opposite halves of the continent, suggesting the convergent evolution of independent bushbuck species. In this study, we provide a species tree reconstruction for the genus Tragelaphus and show that Scriptus and Sylvaticus are reciprocally monophyletic at nuclear DNA loci, comprising a single species across its African range. Given that mtDNA will sort into species-specific lineages more quickly than nuclear DNA, only an ancient interspecific hybridization event between a female from a now-extinct Tragelaphus species and a proto-Scriptus bushbuck male can reconcile the mito-nuclear incongruence. This extinct species diverged from the nyala (T. angasii) in the Pliocene about 4.1 million years ago. This study adds to an increasing body of evidence that suggests interspecific hybridization may be more common than previously thought.

A near-chromosome-scale genome assembly of the gemsbok (Oryx gazella): an iconic antelope of the Kalahari desert

Background

The gemsbok (Oryx gazella) is one of the largest antelopes in Africa. Gemsbok are heterothermic and thus highly adapted to live in the desert, changing their feeding behavior when faced with extreme drought and heat. A high-quality genome sequence of this species will assist efforts to elucidate these and other important traits of gemsbok and facilitate research on conservation efforts.

Findings

Using 180 Gbp of Illumina paired-end and mate-pair reads, a 2.9 Gbp assembly with scaffold N50 of 1.48 Mbp was generated using SOAPdenovo. Scaffolds were extended using Chicago library sequencing, which yielded an additional 114.7 Gbp of DNA sequence. The HiRise assembly using SOAPdenovo + Chicago library sequencing produced a scaffold N50 of 47 Mbp and a final genome size of 2.9 Gbp, representing 90.6% of the estimated genome size and including 93.2% of expected genes according to Benchmarking Universal Single-Copy Orthologs analysis. The Reference-Assisted Chromosome Assembly tool was used to generate a final set of 47 predicted chromosome fragments with N50 of 86.25 Mbp and containing 93.8% of expected genes. A total of 23,125 protein-coding genes and 1.14 Gbp of repetitive sequences were annotated using de novo and homology-based predictions.

Conclusions

Our results provide the first high-quality, chromosome-scale genome sequence assembly for gemsbok, which will be a valuable resource for studying adaptive evolution of this species and other ruminants.

Pan-African genetic structure in the African buffalo (Syncerus caffer): investigating intraspecific divergence

The African buffalo (Syncerus caffer) exhibits extreme morphological variability, which has led to controversies about the validity and taxonomic status of the various recognized subspecies. The present study aims to clarify these by inferring the pan-African spatial distribution of genetic diversity, using a comprehensive set of mitochondrial D-loop sequences from across the entire range of the species. All analyses converged on the existence of two distinct lineages, corresponding to a group encompassing West and Central African populations and a group encompassing East and Southern African populations. The former is currently assigned to two to three subspecies (S. c. nanus, S. c. brachyceros, S. c. aequinoctialis) and the latter to a separate subspecies (S. c. caffer). Forty-two per cent of the total amount of genetic diversity is explained by the between-lineage component, with one to seventeen female migrants per generation inferred as consistent with the isolation-with-migration model. The two lineages diverged between 145 000 to 449 000 years ago, with strong indications for a population expansion in both lineages, as revealed by coalescent-based analyses, summary statistics and a star-like topology of the haplotype network for the S. c. caffer lineage. A Bayesian analysis identified the most probable historical migration routes, with the Cape buffalo undertaking successive colonization events from Eastern toward Southern Africa. Furthermore, our analyses indicate that, in the West-Central African lineage, the forest ecophenotype may be a derived form of the savanna ecophenotype and not vice versa, as has previously been proposed. The African buffalo most likely expanded and diverged in the late to middle Pleistocene from an ancestral population located around the current-day Central African Republic, adapting morphologically to colonize new habitats, hence developing the variety of ecophenotypes observed today.

Population genetics of the malaria vector Anopheles aconitus in China and Southeast Asia

Anopheles aconitus is a well-known vector of malaria and is broadly distributed in the Oriental Region, yet there is no information on its population genetic characteristics. In this study, the genetic differentiation among populations was examined using 140 mtDNA COII sequences from 21 sites throughout Southern China, Myanmar, Vietnam, Thailand, Laos and Sri Lanka. The population in Sri Lanka has characteristic rDNA D3 and ITS2, mtDNA COII and ND5 haplotypes, and may be considered a distinct subspecies. Clear genetic structure was observed with highly significant genetic variation present among population groups in Southeast Asia. The greatest genetic diversity exists in Yunnan and Myanmar population groups. All population groups are significantly different from one another in pairwise Fst values, except Northern Thailand with Central Thailand. Mismatch distributions and extremely significant F(s) values suggest that the populations passed through a recent demographic expansion. These patterns are discussed in relation to the likely biogeographic history of the region and compared to other Anopheles species.

Species identification of Tanzanian antelopes using DNA barcoding

Efficient tools for consistent species identification are important in wildlife conservation as it can provide information on the levels of species exploitation and assist in solving forensic-related problems. In this study, we evaluated the effectiveness of the mitochondrial cytochrome c oxidase subunit I (COI) barcode in species identification of Tanzanian antelope species. A 470 base-pair region of the COI gene was examined in 95 specimens representing 20 species of antelopes, buffalo and domestic Bovidae. All the Tanzanian species showed unique clades, and sequence divergence within species was <1%, whereas divergence between species ranged from 6.3% to 22%. Lowest interspecific divergence was noted within the Tragelaphus genus. Neighbour-joining phylogenetic analyses demonstrated that the examined COI region provided correct and highly supported species clustering using short fragments down to 100 base-pair lengths. This study demonstrates that even short COI fragments can efficiently identify antelope species, thus demonstrating its high potential for use in wildlife conservation activities.

Microsatellite analysis of genetic variation and population subdivision for the black muntjac, Muntiacus crinifrons

The black muntjac (Muntiacus crinifrons) is a rare deer found only in a restricted region in east China. Recent studies of mitochondrial DNA diversity have shown a markedly low level of nucleotide diversity for the species, and the Suichang population was genetically differentiated from the two other populations, in Huangshan and Tianmushan mountains. In this study, we extended the analysis of genetic diversity and population subdivision for the black muntjac using data from 11 highly polymorphic nuclear DNA microsatellite loci. Contrary to the results based on mtDNA data, the microsatellite loci revealed that the black muntjac retained a rather high nuclear genetic diversity (overall average H (E) = 0.78). Nevertheless, both types of markers supported the idea that the extant black muntjac population is genetically disrupted (overall phi (ST) = 0.16 for mtDNA and overall F (ST) = 0.053 for microsatellite, both P < 0.001). The correlation between genetic differentiation and geographic distance was not significant (Mantel test; P > 0.05), implying that the patterns of genetic differentiation observed in this study might result from recent habitat fragmentation or loss. Based on the results from the mtDNA and nuclear DNA data sets, two management units were defined for the species, Huangshan/Tianmushan and Suichang. We also recommend that a new captive population be established with individuals from the Suichang region as a founder source.

Mitochondrial DNA diversity and population structure among southern right whales (Eubalaena australis)

The population structure and mitochondrial (mt) DNA diversity of southern right whales (Eubalaena australis) are described from 146 individuals sampled on 4 winter calving grounds (Argentina, South Africa, Western Australia, and the New Zealand sub-Antarctic) and 2 summer feeding grounds (South Georgia and south of Western Australia). Based on a consensus region of 275 base pairs of the mtDNA control region, 37 variable sites defined 37 unique haplotypes, of which only one was shared between regional samples of the Indo-Pacific and South Atlantic Oceans. Phylogenetic reconstruction of the southern right whale haplotypes revealed 2 distinct clades that differed significantly in frequencies between oceans. An analysis of molecular variance confirmed significant overall differentiation among the 4 calving grounds at both the haplotype and the nucleotype levels (F(ST) = 0.159; Phi(ST) = 0.238; P < 0.001). Haplotype diversity was significantly lower in the Indo-Pacific (h = 0.701 +/- 0.037) compared with the South Atlantic (h = 0.948 +/- 0.013), despite a longer history of exploitation and larger catches in the South Atlantic. In fact, the haplotype diversity in the Indo-Pacific basin was similar to that of the North Atlantic right whale that currently numbers about 300 animals. Multidimensional scaling of genetic differentiation suggests that gene flow occurred primarily between adjacent calving grounds within an ocean basin, with mixing of lineages from different calving grounds occurring on feeding grounds.

Mitochondrial phylogeography and genetic diversity of Tibetan gazelle (Procapra picticaudata): implications for conservation

The Tibetan gazelle (Procapra picticaudata) is a threatened species and distributed on the Qinghai-Tibet Plateau of China (Qinghai Province, Tibet Autonomous Region and the adjacent Gansu Province, Sichuan Province, and Xinjiang Uigur Autonomous Region). Small peripheral populations of Tibetan gazelle were once found in northern Sikkim and Ladakh, but now these are close to extinction. To describe the evolutionary history and to assess the genetic diversity within this monotypic species and population structure among different geographic locations in China, we sequenced mitochondrial DNA from the control region (CR) and cytochrome (cyt) b gene for 46 individuals from 12 geographic localities in Qinghai, Tibet, Xinjiang, Gansu, and Sichuan. A total of 25 CR haplotypes and 16 cyt b haplotypes were identified from these gazelle samples. CR haplotype diversity (0.98+/-0.01) and nucleotide diversity (0.08+/-0.009) were both high. Phylogenetic trees indicate that the Tibetan gazelle in China can be divided into three main clades: Tibet, Sichuan (SCH) and Qinghai-Arjin Shan-Kekexili (QH-ARJ-KKXL). Analysis of molecular variance (AMOVA) and network analysis consistently support this geographic structure in both datasets. Significant differentiation between populations argues for the presence of management units (MUs). Such differentiation may reflect a geographic separation resulting from the uplift of the Qinghai-Tibet Plateau during the Late Pliocene and Pleistocene. Mismatch distribution analysis implies that Tibetan gazelle has undergone complex population changes. We suggest that the present population structure has resulted from habitat fragmentation during the recent glacial period on the Qinghai-Tibet Plateau and population expansion from glacial refugia after the glacial period. It is likely that the present populations of Tibetan gazelle exhibit a pattern reminiscent of several bottlenecks and expansions in the recent past.

No suggestion of hybridization between the vulnerable black-faced impala (Aepyceros melampus petersi) and the common impala (A. m. melampus) in Etosha National Park, Namibia

There are two recognized subspecies of impala in sub-Saharan Africa: the common impala (Aepyceros melampus melampus) -- widespread in southern and east Africa -- and the vulnerable black-faced impala (A. m. petersi) -- found naturally in only a small enclave in southwest Africa. The Etosha National Park (NP) in Namibia harbours the largest and only protected-area population of black-faced impala, numbering some 1500 individuals. Due to translocations of the exotic common impala to commercial farms in Namibia during the past decades, the black-faced impala in Etosha is faced with the potentially serious threat of hybridization posed by secondary contact with the common impala inhabiting bordering farms. Using eight microsatellite DNA markers, we analysed 127 black-faced impala individuals from the five subpopulations in Etosha NP, to determine the degree, if any, of hybridization within the park. We found that (a) the black-faced impala were highly genetically differentiated from the common impala (pairwise theta-values ranged from 0.18 to 0.39 between subspecies; overall value = 0.27) and (b) black-faced samples showed high levels of genetic variability [average expected heterozygosity (H(E)) = 0.61 +/- 0.01 SE], although not as high as that observed in the common impala (average H(E) = 0.69 +/- 0.02 SE). (c) No hybridization between the subspecies in Etosha was suggested. A Bayesian Markov Chain Monte Carlo approach revealed clear distinction of individuals into groups according to their subspecies of origin, with a zero level of 'genetic admixture' among subspecies.

Molecular architecture of Pipistrellus pipistrellus/Pipistrellus pygmaeus complex (Chiroptera: Vespertilionidae): further cryptic species and Mediterranean origin of the divergence

Previous genetic analyses have demonstrated that two phonic types of one of the most common European bats, the Common pipistrelle, belong to distinct species, although they are almost identical morphologically (45 kHz Pipistrellus pipistrellus and 55 kHz Pipistrellus pygmaeus). To reconstruct the history of the species complex and explain the codistribution of both forms in Europe and the Mediterranean, we performed phylogenetic analysis based on a 402-bp portion of the cytochrome b gene. Particular attention was paid to the eastern and southern parts of the range where no data were available. We found further distinctive allopatric haplotypes from Libya and Morocco. The difference of about 6-7% described in the Libyan population suggests the occurrence of a new species in the southern Mediterranean. The species status of Moroccan population is also discussed. The phylogeographic patterns obtained and analysis of fossil records support the hypothesis of expansion of both species into Europe from the Mediterranean region during the Holocene. The allopatric speciation model fits our data best. The paleobiographic scenario envisaged is corroborated also by molecular clock estimations and correlations with Late Neogene environmental changes in the Mediterranean region which ended with the Messinian salinity crisis.

Phylogeography and population structure of the common warthog (Phacochoerus africanus) inferred from variation in mitochondrial DNA sequences and microsatellite loci

Global climate fluctuated considerably throughout the Pliocene and Pleistocene, influencing the evolutionary history of a wide range of species. Using both mitochondrial sequences and microsatellites, we have investigated the evolutionary consequences of such environmental fluctuation for the patterns of genetic variation in the common warthog, sampled from 24 localities in Africa. In the sample of 181 individuals, 70 mitochondrial DNA haplotypes were identified and an overall nucleotide diversity of 4.0% was observed. The haplotypes cluster in three well-differentiated clades (estimated net sequence divergence of 3.1-6.6%) corresponding to the geographical origins of individuals (i.e. eastern, western and southern African clades). At the microsatellite loci, high polymorphism was observed both in the number of alleles per locus (6-21), and in the gene diversity (in each population 0.59-0.80). Analysis of population differentiation indicates greater subdivision at the mitochondrial loci (FST=0.85) than at nuclear loci (FST=0.20), but both mitochondrial and nuclear loci support the existence of the three warthog lineages. We interpret our results in terms of the large-scale climatic fluctuations of the Pleistocene.

Patterns of linkage disequilibrium in mitochondrial DNA of 16 ruminant populations

Mitochondrial DNA (mtDNA) is a widely employed molecular tool in phylogeography, in the inference of human evolutionary history, in dating the domestication of livestock and in forensic science. In humans and other vertebrates the popularity of mtDNA can be partially attributed to an assumption of strict maternal inheritance, such that there is no recombination between mitochondrial lineages. The recent demonstration that linkage disequilibrium (LD) declines as a function of distance between polymorphic sites in hominid mitochondrial genomes has been interpreted as evidence of recombination between mtDNA haplotypes, and hence nonclonal inheritance. However, critics of mtDNA recombination have suggested that this association is an artefact of an inappropriate measure of LD or of sequencing error, and subsequent studies of other populations have failed to replicate the initial finding. Here we report the analysis of 16 ruminant populations and present evidence that LD significantly declines with distance in five of them. A meta-analysis of the data indicates a nonsignificant trend of LD declining with distance. Most of the earlier criticisms of patterns between LD and distance in hominid mtDNA are not applicable to this data set. Our results suggest that either ruminant mtDNA is not strictly clonal or that compensatory selection has influenced patterns of variation at closely linked sites within the mitochondrial control region. The potential impact of these processes should be considered when using mtDNA as a tool in vertebrate population genetic, phylogenetic and forensic studies.

The genetic impact of demographic decline and reintroduction in the wild boar (Sus scrofa): a microsatellite analysis

The reintroduction of wild boar from central Europe after World War II has contributed substantially to the range expansion of this species in Italy, where indiscriminate hunting in earlier times resulted in extreme demographic reduction. However, the genetic impact of such processes is not well-understood. In this study, 105 individuals from Italian and Hungarian wild boar populations were characterized for nine autosomal microsatellite loci. The Hungarian samples, and two central Italian samples from protected areas (parks) where reintroduction is not documented, were assumed to be representative of the genetic composition of the source and the target populations in the reintroduction process, respectively. Animals hunted in the wild in the Florence area of Tuscany (Italy) were then studied to identify the effects of reintroduction. The results we obtained can be summarized as follows: (i) none of the populations analysed shows genetic evidence of demographic decline; (ii) the three parental populations from Italy and Hungary are genetically distinct; however, the low level of divergence appears in conflict with the naming of the Italian and the European subspecies (Sus scrofa majori and Sus scrofa scrofa, respectively); in addition, the Italian groups appear to be as divergent from each other as they are from the Hungarian population; (iii) most of the individuals hunted near Florence are genetically intermediate between the parental groups, suggesting that hybridization has occurred in this area, the average introgression of Hungarian genotypes is 13%, but approximately 45% of the genetic pool of these individuals can not be directly attributed to any of the parental populations we analysed; (iv) analysis of microsatellite loci, though in a limited number, is an important tool for estimating the genetic effect of reintroduction in the wild boar, and therefore for the development of conservation and management strategies for this species.

Genetic structure and phylogeography of the takin (Budorcas taxicolor) as inferred from mitochondrial DNA sequences

To investigate the genetic structure within the takin (Budorcas taxicolor), a little-studied and threatened species, portions of the mitochondrial DNA control region and cytochrome b gene were sequenced from 40 takins from three subspecies and 10 separate geographical locations in China. Overall variation among individuals was low. Both neighbor-joining and parsimony analyses did separate animals into three distinct clades reflecting the subspecific definitions, with only one locality grouping with a clade distinct from its purported subspecies. Animals within each clade were genetically similar to each other, which we hypothesize is due to the mobility of the animals leading to significant gene flow within each region until recent habitat alteration limited movements. We also suggest that genetic bottlenecks survived by the species since the end of the Pleistocene may have contributed to overall low levels of variability within the species.

Inferring the phylogeny of disjunct populations of the azure-winged magpie Cyanopica cyanus from mitochondrial control region sequences

The azure-winged magpie (AWM), Cyanopica cyanus, is found in Asia and Iberia. This remarkable disjunct distribution has been variously explained by either the sixteenth-century introduction of birds into Iberia from the Far East, or by the loss of individuals from the central part of their range as a result of Pleistocene glaciations. We have used the mitochondrial control region to undertake a molecular phylogenetic analysis of the AWM, with sequences examined from individuals collected from across the current distribution range and incorporating representatives of all currently defined subspecies. The Western birds are genetically distinct from their Asian congeners and their divergence is basal in the phylogenetic tree. This indicates that the AWM is native to Iberia and not the result of a recent introduction from Asia. In Asia, two major mitochondrial DNA lineages were identified. These correspond to an Inland Asia group and a Pacific Seaboard group, and are separated topographically by the Da Hingan Ling mountains and the Yellow Sea. Molecular clock estimates suggest that these divergences are associated with Pleistocene glaciations. Furthermore, our data do not support the current classification of the AWM into 10 subspecies, as defined based on morphology and geographical distribution.

Population structure of the African savannah elephant inferred from mitochondrial control region sequences and nuclear microsatellite loci

Two hundred and thirty-six mitochondrial DNA nucleotide sequences were used in combination with polymorphism at four nuclear microsatellite loci to assess the amount and distribution of genetic variation within and between African savannah elephants. They were sampled from 11 localities in eastern, western and southern Africa. In the total sample, 43 haplotypes were identified and an overall nucleotide diversity of 2.0% was observed. High levels of polymorphism were also observed at the microsatellite loci both at the level of number of alleles and gene diversity. Nine to 14 alleles per locus across populations and 44 alleles in the total sample were found. The gene diversity ranged from 0.51 to 0.72 in the localities studied. An analysis of molecular variance showed significant genetic differentiation between populations within regions and also between regions. The extent of subdivision between populations at the mtDNA control region was approximately twice as high as shown by the microsatellite loci (mtDNA F(ST) = 0.59; microsatellite R(ST) = 0.31). We discuss our results in the light of Pleistocene refugia and attribute the observed pattern to population divergence in allopatry accompanied by a recent population admixture following a recent population expansion.

An exceptional case of historical outbreeding in African sable antelope populations

Empirical investigations of intraspecific outbreeding and subsequent introgressive hybridization in natural populations are rare, particularly among conspecific populations of large mammals. Using mitochondrial DNA data [partial control region (496 basepairs - bp) and cytochrome b gene (343 bp) sequences analysed from 95 individuals representing 17 sampling locations scattered through the African miombo (Brachystegia) woodland ecosystem] and phylogeographical statistical procedures (gene genealogy, nested cladistic and admixture proportion analyses), we (i) give a detailed dissection of the geographical genetic structure of Hippotragus niger; (ii) infer the processes and events potentially involved in the population history; and (iii) trace extensive introgressive hybridization in the species. The present-day sable antelope population shows a tripartite pattern of genetic subdivision representing West Tanzanian, Kenya/East Tanzanian and Southern Africa locations. Nested clade analysis revealed that past allopatric fragmentation, caused probably by habitat discontinuities associated with the East African Rift Valley system, together with intermediary episodic long-distance colonization and restricted, recurrent gene flow have played an predominant role in shaping the extent of maternal genetic diversity (10.4%) and population structure. An extensive (average rate of admixture = 20.0%), but geographically circumscribed and unidirectional hybridization event in the past was inferred, resulting in an extreme (the highest discovered so far in mammals) intraspecific difference of 18.2% among morphologically monotypic sable antelopes from West Tanzania. The results are used to provide an evolutionary framework within which taxonomic implications and conservation decisions can be evaluated.

Patterns of human diversity, within and among continents, inferred from biallelic DNA polymorphisms

Previous studies have reported that about 85% of human diversity at Short Tandem Repeat (STR) and Restriction Fragment Length Polymorphism (RFLP) autosomal loci is due to differences between individuals of the same population, whereas differences among continental groups account for only 10% of the overall genetic variance. These findings conflict with popular notions of distinct and relatively homogeneous human races, and may also call into question the apparent usefulness of ethnic classification in, for example, medical diagnostics. Here, we present new data on 21 Alu insertions in 32 populations. We analyze these data along with three other large, globally dispersed data sets consisting of apparently neutral biallelic nuclear markers, as well as with a beta-globin data set possibly subject to selection. We confirm the previous results for the autosomal data, and find a higher diversity among continents for Y-chromosome loci. We also extend the analyses to address two questions: (1) whether differences between continental groups, although small, are nevertheless large enough to confidently assign individuals to their continent on the basis of their genotypes; (2) whether the observed genotypes naturally cluster into continental or population groups when the sample source location is ignored. Using a range of statistical methods, we show that classification errors are at best around 30% for autosomal biallelic polymorphisms and 27% for the Y chromosome. Two data sets suggest the existence of three and four major groups of genotypes worldwide, respectively, and the two groupings are inconsistent. These results suggest that, at random biallelic loci, there is little evidence, if any, of a clear subdivision of humans into biologically defined groups.

Restricted effective queen dispersal at a microgeographic scale in polygynous populations of the ant Formica exsecta

Ecological constraints on effective dispersal have been suggested to be a key factor influencing social evolution in animal societies as well as the shift from single queen colonies (monogyny) to multiple queen colonies (polygyny) in ants. However, little is known about the effective dispersal patterns of ant queens. Here we investigate the microgeographic genetic structure of mitochondrial haplotypes in polygynous populations of the ant Formica exsecta, both between pastures and among nests within pastures. An analysis of molecular variance revealed a very high genetic differentiation (phiST = 0.72) between pastures, indicating that queens rarely disperse successfully between pastures, despite the fact that pastures were sometimes as close as 1 km. Most of the pastures contained only a single haplotype, and haplotypes were frequently distinct between nearby pastures and even between groups of nests within the same pasture. In the three pastures that contained several haplotypes, haplotypes were not randomly distributed, the genetic differentiation between nests being phiST = 0.17, 0.52, and 0.69. This indicates that most queens are recruited within their parental colonies. However, a large proportion of nests contained more than one haplotype, demonstrating that colonies will sometimes accept foreign queens. The relatedness of mitochondrial genes among nestmates varied between 0.62 and 0.75 when relatedness was measured within each pasture and ranged between 0.72 and 1.0 when relatedness was assessed with all pastures as a reference population. Neighboring nests were more genetically similar than distant ones, and there was significant isolation by distance. This pattern may be due to new nests being formed by budding or by limited effective queen dispersal, probably on foot between neighboring nests. These results show that effective queen dispersal is extremely restricted even at a small geographical scale, a pattern consistent with the idea that ecological constraints are an important selective force leading to the evolution and maintenance of polygyny.

Phylogeography of the asian elephant (Elephas maximus) based on mitochondrial DNA

Populations of the Asian elephant (Elephas maximus) have been reduced in size and become highly fragmented during the past 3,000 to 4,000 years. Historical records reveal elephant dispersal by humans via trade and war. How have these anthropogenic impacts affected genetic variation and structure of Asian elephant populations? We sequenced mitochondrial DNA (mtDNA) to assay genetic variation and phylogeography across much of the Asian elephant's range. Initially we compare cytochrome b sequences (cyt b) between nine Asian and five African elephants and use the fossil-based age of their separation (approximately 5 million years ago) to obtain a rate of about 0.013 (95% CI = 0.011-0.018) corrected sequence divergence per million years. We also assess variation in part of the mtDNA control region (CR) and adjacent tRNA genes in 57 Asian elephants from seven countries (Sri Lanka, India, Nepal, Myanmar, Thailand, Malaysia, and Indonesia). Asian elephants have typical levels of mtDNA variation, and coalescence analyses suggest their populations were growing in the late Pleistocene. Reconstructed phylogenies reveal two major clades (A and B) differing on average by HKY85/gamma-corrected distances of 0.020 for cyt b and 0.050 for the CR segment (corresponding to a coalescence time based on our cyt b rate of approximately 1.2 million years). Individuals of both major clades exist in all locations but Indonesia and Malaysia. Most elephants from Malaysia and all from Indonesia are in well-supported, basal clades within clade A. thus supporting their status as evolutionarily significant units (ESUs). The proportion of clade A individuals decreases to the north, which could result from retention and subsequent loss of ancient lineages in long-term stable populations or, perhaps more likely, via recent mixing of two expanding populations that were isolated in the mid-Pleistocene. The distribution of clade A individuals appears to have been impacted by human trade in elephants among Myanmar, Sri Lanka, and India, and the subspecies and ESU statuses of Sri Lankan elephants are not supported by molecular data.

Phylogeography and conservation of impala and greater kudu

The phylogeography of the bush habituated African bovid species impala (Aepyceros melampus) and greater kudu (Tragelaphus strepsiceros) is investigated using mitochondrial DNA (mtDNA) markers. Combined analysis of individual lineages, relationships and population genetics suggest a colonization process from Southern Africa toward Eastern regions in the greater kudu. Results are less clear for the impala, although remaining consistent with a similar pattern of historical dispersion. The study reveals a similar pattern, that is a marked divergence of lineages from South-western Africa relative to other regions. This pattern is opposed to previously published findings in other African bovid species. In the impala, the genetically isolated region is consistent with morphology because it is recognized as the subspecies A. m. petersi, the black-faced impala. In contrast, the similar split of South-western mitochondrial lineages was not expected in the greater kudu on the basis of morphology. Both species show a significant population genetic differentiation. Beyond their phylogeographical value, our results should raise conservation concerns about South-western populations of both species. The black-faced impala is categorized as vulnerable and our data show indications of hybridization with common impala A. m. melampus. The previously unrecognized genetic status of the South-western kudus could also imply conservation regulations.

Large sequence divergence of mitochondrial DNA genotypes of the control region within populations of the african antelope, kob (Kobus kob)

The genetic differentiation among kob populations (Kobus kob) representing two recognized subspecies was examined using mitochondrial control region sequences. Two distinct lineages (estimated sequence divergence of 9.8%) exhibited different geographical distributions and do not coincide with previously recognized ranges of subspecies. The presence of the two lineages was further supported with sequences of mitochondrial cytochrome b gene. One lineage was predominant in the west and southern ranges of the populations studied and the other was commonly found in a more northern distribution (Murchison populations) in Uganda. Murchison and the geographically intermediate Toro populations (Uganda) represented the area of overlap. The existence of the two lineages in the area of overlap is hypothesized to have resulted from a range expansion and secondary contact of the two lineages of kob that evolved in allopatry. The existence of the kob during the Pleistocene offers a plausible explanation for the observed biogeographic pattern. Our mitochondrial data reveal two examples of discordance between a gene tree and presumed species tree as: (i) the two lineages co-occur in the kob subspecies, Kobus kob thomasi (Uganda kob); and (ii) the puku, which was included in the analysis because of its controversial taxonomic status (currently recognized as a distinct species from the kob), is paraphyletic with respect to the kob. Significant degrees of heterogeneity were detected between populations. Relatively high genetic variation was observed in the populations, however, the inclusion of distinct lineages influences the population structure and nucleotide diversity of the kob populations.

Microsatellite analysis of genetic diversity in African buffalo (Syncerus caffer) populations throughout Africa

Genetic diversity in nine African buffalo (Syncerus caffer) populations throughout Africa was analysed with 14 microsatellites to study the effects of rinderpest epidemics and habitat fragmentation during the 20th century. A gradient of declining expected heterozygosity was observed among populations in Save Valley Conservancy (Zimbabwe), and northern and southern Kruger National Park (South Africa). This was explained by a high mortality in northern Kruger National Park during the rinderpest pandemic at the end of the 19th century followed by recolonization from neighbouring populations, resulting in intermediate heterozygosity levels in northern Kruger National Park. In other populations expected heterozygosity was very high, indicating that rinderpest and recent habitat fragmentation had a limited effect on genetic diversity. From expected heterozygosity, estimates of long-term effective population size were derived. Migration rates among populations in eastern and southern Africa were very high, as shown by a weak isolation by distance and significant correlation in allele frequencies between populations. However, there were indications that dry habitats could limit migration. Genetic distances within buffalo in central Africa were relatively large, supporting their status as distinct subspecies. Finally, it was observed that the higher polymorphic microsatellites were less sensitive at detecting isolation by distance and differences in Ne, which may be a result of the high mutation pressure at these loci.

Consequences of recombination on traditional phylogenetic analysis

We investigate the shape of a phylogenetic tree reconstructed from sequences evolving under the coalescent with recombination. The motivation is that evolutionary inferences are often made from phylogenetic trees reconstructed from population data even though recombination may well occur (mtDNA or viral sequences) or does occur (nuclear sequences). We investigate the size and direction of biases when a single tree is reconstructed ignoring recombination. Standard software (PHYLIP) was used to construct the best phylogenetic tree from sequences simulated under the coalescent with recombination. With recombination present, the length of terminal branches and the total branch length are larger, and the time to the most recent common ancestor smaller, than for a tree reconstructed from sequences evolving with no recombination. The effects are pronounced even for small levels of recombination that may not be immediately detectable in a data set. The phylogenies when recombination is present superficially resemble phylogenies for sequences from an exponentially growing population. However, exponential growth has a different effect on statistics such as Tajima's D. Furthermore, ignoring recombination leads to a large overestimation of the substitution rate heterogeneity and the loss of the molecular clock. These results are discussed in relation to viral and mtDNA data sets.

Two genetically distinct lineages of the sika deer, Cervus nippon, in Japanese islands: comparison of mitochondrial D-loop region sequences

To investigate genetic diversity among populations of the sika deer, Cervus nippon, nucleotide sequences (705-824 bases) of the mitochondrial D-loop regions were determined in animals from 13 localities in the Japanese islands. Phylogenetic trees constructed by the sequences indicated that the Japanese sika deer is separated into two distinct lineages: the northern Japan group (the Hokkaido island and most of the Honshu mainland) and the southern Japan group (a part of the southern Honshu mainland, the Kyushu island, and small islands around the Kyushu island). All sika deer examined in this study shared four to seven units of repetitive sequences (37 to 40 bases each) within the D-loop sequences. The number of tandem repeats was different among the populations, and it was specific to each population. Six or seven repeats occurred in populations of the northern Japan group, while four or five repeats occurred in populations of the southern Japan group. Each repeat unit included several nucleotide substitutions, compared with others, and 26 types were identified from 31 animals. Sequences of the first, second, and third units in arrays were clearly different between the northern and the southern groups. Based on these D-loop data, colonization and separation of the sika deer populations in the Japanese islands were estimated to have occurred less than 0.5 million years before present. Our results provide an invaluable insight into better understanding the evolutionary history, phylogeny, taxonomy, and population genetics of the sika deer.