Effective population size dynamics reveal impacts of historic climatic events and recent anthropogenic pressure in African elephants.

Two hundred years of elephant hunting for ivory, peaking in 1970-1980s, caused local extirpations and massive population declines across Africa. The resulting genetic impacts on surviving populations have not been studied, despite the importance of understanding the evolutionary repercussions of such human-mediated events on this keystone species. Using Bayesian coalescent-based genetic methods to evaluate time-specific changes in effective population size, we analysed genetic variation in 20 highly polymorphic microsatellite loci from 400 elephants inhabiting the greater Samburu-Laikipia region of northern Kenya. This area experienced a decline of between 80% and 90% in the last few decades when ivory harvesting was rampant. The most significant change in effective population size, however, occurred approximately 2500 years ago during a mid-Holocene period of climatic drying in tropical Africa. Contrary to expectations, detailed analyses of four contemporary age-based cohorts showed that the peak poaching epidemic in the 1970s caused detectable temporary genetic impacts, with genetic diversity rebounding as juveniles surviving the poaching era became reproductively mature. This study demonstrates the importance of climatic history in shaping the distribution and genetic history of a keystone species and highlights the utility of coalescent-based demographic approaches in unravelling ancestral demographic events despite a lack of ancient samples. Unique insights into the genetic signature of mid-Holocene climatic change in Africa and effects of recent poaching pressure on elephants are discussed.

Noninvasive genotyping and Mendelian analysis of microsatellites in African savannah elephants.

We obtained fresh dung samples from 202 (133 mother-offspring pairs) savannah elephants (Loxodonta africana) in Samburu, Kenya, and genotyped them at 20 microsatellite loci to assess genotyping success and errors. A total of 98.6% consensus genotypes was successfully obtained, with allelic dropout and false allele rates at 1.6% (n = 46) and 0.9% (n = 37) of heterozygous and total consensus genotypes, respectively, and an overall genotyping error rate of 2.5% based on repeat typing. Mendelian analysis revealed consistent inheritance in all but 38 allelic pairs from mother-offspring, giving an average mismatch error rate of 2.06%, a possible result of null alleles, mutations, genotyping errors, or inaccuracy in maternity assignment. We detected no evidence for large allele dropout, stuttering, or scoring error in the dataset and significant Hardy-Weinberg deviations at only two loci due to heterozygosity deficiency. Across loci, null allele frequencies were low (range: 0.000-0.042) and below the 0.20 threshold that would significantly bias individual-based studies. The high genotyping success and low errors observed in this study demonstrate reliability of the method employed and underscore the application of simple pedigrees in noninvasive studies. Since none of the sires were included in this study, the error rates presented are just estimates.

Mitochondrial DNA variation of the common hippopotamus: evidence for a recent population expansion.

Mitochondrial DNA control region sequence variation was obtained and the population history of the common hippopotamus was inferred from 109 individuals from 13 localities covering six populations in sub-Saharan Africa. In all, 100 haplotypes were defined, of which 98 were locality specific. A relatively low overall nucleotide diversity was observed (pi = 1.9%), as compared to other large mammals so far studied from the same region. Within populations, nucleotide diversity varied from 1.52% in Zambia to 1.92% in Queen Elizabeth and Masai Mara. Overall, low but significant genetic differentiation was observed in the total data set (F(ST) = 0.138; P = 0.001), and at the population level, patterns of differentiation support previously suggested hippopotamus subspecies designations (F(CT) = 0.103; P = 0.015). Evidence that the common hippopotamus recently expanded were revealed by: (i) lack of clear geographical structure among haplotypes, (ii) mismatch distributions of pairwise differences (r = 0.0053; P = 0.012) and site-frequency spectra, (iii) Fu's neutrality statistics (F(S) = -155.409; P < 0.00001) and (iv) Fu and Li's statistical tests (D* = -3.191; P < 0.01, F* = -2.668; P = 0.01). Mismatch distributions, site-frequency spectra and neutrality statistics performed at subspecies level also supported expansion of Hippopotamus amphibius across Africa. We interpret observed common hippopotamus population history in terms of Pleistocene drainage overflow and suggest recognising the three subspecies that were sampled in this study as separate management units in future conservation planning.

Population genetics of the roan antelope (Hippotragus equinus) with suggestions for conservation.

The roan antelope (Hippotragus equinus) is the second largest African antelope, distributed throughout the continent in sub-Saharan savannah habitat. Mitochondrial DNA (mtDNA) control region sequencing (401 bp, n = 137) and microsatellite genotyping (eight loci, n = 137) were used to quantify the genetic variability within and among 18 populations of this species. The within-population diversity was low to moderate with an average mtDNA nucleotide diversity of 1.9% and average expected heterozygosity with the microsatellites of 46%, but significant differences were found among populations with both the mtDNA and microsatellite data. Different levels of genetic resolution were found using the two marker sets, but both lent strong support for the separation of West African populations (samples from Benin, Senegal and Ghana) from the remainder of the populations studied across the African continent. Mismatch distribution analyses revealed possible past refugia for roan in the west and east of Africa. The West African populations could be recognized together as an evolutionarily significant unit (ESU), referable to the subspecies H. e. koba. Samples from the rest of the continent constituted a geographically more diverse assemblage with genetic associations not strictly corresponding to the other recognized subspecies.

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.

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.

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.

Recombination among multiple mitochondrial pseudogenes from a Passerine genus.

PCR products of a fragment of the mitchondrial protein coding subunit 5 of NADH-dehydrogenase (ND5) from eight individuals representing five species of the South American bird genus Conirostrum were cloned. The 130 clones, which were subsequently sequenced, constituted 55 different sequences. Due to the observed differences in substitution patterns 58% of the cloned sequences were identified as pseudogenes. Recombination could be traced in 19% of the inferred nuclear pseudogenes, but this figure probably represents a significant underestimation of the factual recombination events. The nonrecombined pseudogenes consisted of multiple haplotypes found to diverge from 1 to 16% from the mitochondrial gene. The number of mitochondrial nuclear copies and their apparent frequent recombination suggest that pseudogenes constitute a serious potential risk in confounding phylogenetic studies and population genetic analysis.

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.

Hidden morphological support for the phylogenetic placement of Pseudoryx nghetinhensis with bovine bovids: a combined analysis of gross anatomical evidence and DNA sequences from five genes.

The saola, Pseudoryx nghetinhensis, was unknown to science until its formal description in 1993. This endangered species is a member of the ruminant artiodactyl family Bovidae (cattle, sheep, goats, and antelopes). However, given its puzzling combination of morphological traits, the specific affinities of Pseudoryx within Bovidae are controversial. A preliminary genetic investigation suggested that Pseudoryx should be placed in the subfamily Bovinae (cattle, buffaloes, spiral-horned antelopes, and nilgai), but a recent cladistic analysis of skeletal and dental characters allied Pseudoryx with caprine bovids (sheep, goats, musk oxen, goat antelopes, and Pantholops). The morphological and molecular hypotheses differ in assigning the saola to either of the two most divergent clades of Bovidae. In this report, phylogenetic analyses of DNA sequences from five genes are used to test these alternatives. Protein coding regions, introns, and ribosomal DNAs from the nuclear and mitochondrial genomes discount the hypothesis that Pseudoryx is a close relative of Caprinae. Instead, combined analyses of the DNA data and published morphological evidence place Pseudoryx with Bovini (cattle and buffaloes), a subclade of Bovinae. In a separate analysis, the matrix of morphological characters links Pseudoryx with caprine bovids, but in the context of the molecular data, the gross anatomical evidence strongly supports a grouping of Pseudoryx with Bovinae. Surprisingly, the morphological partition provides the most character support in the combined analysis. This striking result is obscured by separate analyses of the individual data sets and the taxonomic congruence approach.

Phylogeography of three closely related African bovids (tribe Alcelaphini).

The phylogeography of three species of African bovids, the hartebeest (Alcelaphus buselaphus), the topi (Damaliscus lunatus), and the wildebeest (Connochaetes taurinus), is inferred from sequence variation of 345 sequences at the control region (d-loop) of the mtDNA. The three species are closely related (tribe Alcelaphini) and share similar habitat requirements. Moreover, their former distribution extended over Africa, as a probable result of the expansion of open grassland on the continent during the last 2.5 Myr. A combination of population genetics (diversity and structure) and intraspecific phylogeny (tree topology and relative branch length) methods is used to substantiate scenarios of the species history. Population dynamics are inferred from the distribution of sequence pairwise differences within populations. In the three species, there is a significant structuring of the populations, as shown by analysis of molecular variance (AMOVA) pairwise and hierarchical differentiation estimations. In the wildebeest, a pattern of colonization from southern Africa toward east Africa is consistent with the asymmetric topology of the gene tree, showing a paraphyletic position of southern lineages, as well as their relatively longer branch lengths, and is supported by a progressive decline in population nucleotide diversity toward east Africa. The phylogenetic pattern found in the topi and the hartebeest differs from that of the wildebeest: lineages split into monophyletic clades, and no geographical trend is detected in population diversity. We suggest a scenario where these antelopes, previously with wide pan-African distributions, became extinct except in a few refugia. The hartebeest, and probably also the topi, survived in refugia north of the equator, in the east and the west, respectively, as well as one in the south. The southern refugium furthermore seems to have been the only place where the wildebeest has survived.

Population genetic structure of the African elephant in Uganda based on variation at mitochondrial and nuclear loci: evidence for male-biased gene flow.

A drastic decline has occurred in the size of the Uganda elephant population in the last 40 years, exacerbated by two main factors; an increase in the size of the human population and poaching for ivory. One of the attendant consequences of such a decline is a reduction in the amount of genetic diversity in the surviving populations due to increased effects of random genetic drift. Information about the amount of genetic variation within and between the remaining populations is vital for their future conservation and management. The genetic structure of the African elephant in Uganda was examined using nucleotide variation of mitochondrial control region sequences and four nuclear microsatellite loci in 72 individuals from three localities. Eleven mitochondrial DNA (mtDNA) haplotypes were observed, nine of which were geographically localized. We found significant genetic differentiation between the three populations at the mitochondrial locus while three out of the four microsatellite loci differentiated KV and QE, one locus differentiated KV and MF and no loci differentiated MF and QE. Expected heterozygosity at the four loci varied between 0.51 and 0.84 while nucleotide diversity at the mitochondrial locus was 1.4%. Incongruent patterns of genetic variation within and between populations were revealed by the two genetic systems, and we have explained these in terms of the differences in the effective population sizes of the two genomes and male-biased gene flow between populations.

A case of rapid diversification in the neotropics: phylogenetic relationships among Cranioleuca spinetails (Aves, Furnariidae).

Relationships among the 18-19 species of spinetails of the genus Cranioleuca are difficult to establish. Attempts based on traditional taxonomic characters have failed because of a high degree of homoplasy. Most morphological characters vary independently, producing leap-frog patterns of variation along the Eastern Brazilian Andean track, and behavior and vocalizations vary little. We use mtDNA sequence data from the cyt b and ND2 genes in an attempt to clarify relationships within the genus. We show (i) that Cranioleuca represents a recent burst of speciation and (ii) that a set of species thought by Maijer and Fjeldså (1997) to form a natural group is in fact a paraphyletic assemblage which also includes humid forest species with different pigmentations and vocalizations. However, synapomorphic variation in the sequences is not sufficient to unambiguously resolve the relationships within the genus. Several species (C. baroni, C. antisiensis, C. pyrrhophia, C. albiceps) show more than one haplotype, without any obvious correlation between genetic and geographic or morphological variation, and the different species do not always show reciprocal monophyly in haplotype diversity. Nevertheless, low genetic differentiation characterizes not only allopatric taxa but also some forms which are essentially sympatric, supporting species rank for the former. Our data suggest a recent diversification and proliferation possibly linked to Pleistocene climatic variation and its consequent vegetational shifts, at least in the Andean species.

Diversity of Holocene life forms in fossil glacier ice.

Studies of biotic remains of polar ice caps have been limited to morphological identification of plant pollen and spores. By using sensitive molecular techniques, we now demonstrate a much greater range of detectable organisms; from 2000- and 4000-year-old ice-core samples, we obtained and characterized 120 clones that represent at least 57 distinct taxa and reveal a diversity of fungi, plants, algae, and protists. The organisms derive from distant sources as well as from the local arctic environment. Our results suggest that additional taxa may soon be readily identified, providing a plank for future studies of deep ice cores and yielding valuable information about ancient communities and their change over time.

Serological and genetic characterisation of bovine respiratory syncytial virus (BRSV) indicates that Danish isolates belong to the intermediate subgroup: no evidence of a selective effect on the variability of G protein nucleotide sequence by prior cell culture adaption and passages in cell culture or calves.

Danish isolates of bovine respiratory syncytial virus (BRSV) were characterised by nucleotide sequencing of the G glycoprotein and by their reactivity with a panel of monoclonal antibodies (MAbs). Among the six Danish isolates, the overall sequence divergence ranged between 0 and 3% at the nucleotide level and between 0 and 5% at the amino acid level. Sequence divergences of 7-8%, 8-9% and 2-3% (nucleotide) and 9-11%, 12-16% and 4-6% (amino acid) were obtained in the comparison made between the group of Danish isolates and the previously sequenced 391-2USA, 127UK and 220-69Bel isolates, respectively. Phylogenetic analysis showed that the Danish isolates formed three lineages within a separate branch of the phylogenetic tree. Nevertheless, the Danish isolates were closely related to the 220-69Bel isolate, the prototype of the intermediate antigenic subgroup. The sequencing of the extracellular part of the G gene of additional 11 field BRSV viruses, processed directly from lung samples without prior adaption to cell culture growth, revealed sequence variabilities in the range obtained with the propagated virus. In addition, several passages in cell culture and in calves had no major impact on the nucleotide sequence of the G protein. These findings indicated that the previously established variabilities of the G protein of RS virus isolates were not attributable to mutations induced during the propagation of the virus. The reactivity of the Danish isolates with G protein-specific MAbs were similar to that of the 220-69Bel isolate. Furthermore, the sequence of the immunodominant region was completely conserved among the Danish isolates on one side and the 220-69Bel isolate on the other. When combined, these data strongly suggested that the Danish isolates belong to the intermediate subgroup.