Population genomics of the muskox' resilience in the near absence of genetic variation.

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Characterization of 29 polymorphic microsatellite markers developed by genomic screening of Sumatran rhinoceros (Dicerorhinus sumatrensis).

OBJECTIVE: The Sumatran rhinoceros is critically endangered, with fewer than 100 individuals surviving across its current range. Accurate census estimates of the remaining populations are essential for development and implementation of conservation plans. In order to enable molecular censusing, we here develop microsatellite markers with amplicon sizes of short length, appropriate for non-invasive fecal sampling. RESULTS: Due to limited sample quantity and potential lack of genome-wide diversity, Illumina sequence reads were generated from two Sumatran rhinoceros samples. Genomic screening identified reads with short tandem repeats and loci that were polymorphic within the dataset. Twenty-nine novel polymorphic microsatellite markers were characterized (A = 2.4; HO = 0.30). These were sufficient to distinguish among individuals (PID < 0.0001), and to distinguish among siblings (PID(sib) < 0.0001). Among rhinos in Indonesia, almost all markers were established as polymorphic and effective for genotyping DNA from fecal samples. Notably, the markers amplified and displayed microsatellite polymorphisms using DNA extracted from 11 fecal samples collected non-invasively from wild Sumatran rhinoceros. These microsatellite markers provide an important resource for a census and genetic studies of wild Sumatran rhinos.

The Muskox Lost a Substantial Part of Its Genetic Diversity on Its Long Road to Greenland.

The muskox (Ovibos moschatus) is the largest terrestrial herbivore in the Arctic and plays a vital role in the tundra ecosystem [1-4]. Its range, abundance, and genetic diversity have declined dramatically over the past 30,000 years [5]. Two subspecies are recognized, but little is known about the genetic structure and how this relates to the species history. One unresolved question is how and when the species dispersed into its present range, notably the present strongholds in the Canadian archipelago and Greenland. We used genotyping by sequencing (GBS) data from 116 muskox individuals and genotype likelihood-based methods to infer the genetic diversity and distribution of genetic variation in the species. We identified a basal split separating the two recognized subspecies, in agreement with isolation of the muskox into several refugia in the Nearctic around 21,000 years ago [6], near the last glacial maximum (LGM). In addition, we found evidence of strong, successive founder effects inflicting a progressive loss of genetic diversity as the muskox colonized the insular High Arctic from an unknown Nearctic origin. These have resulted in exceptionally low genetic diversity in the Greenlandic populations, as well as extremely high genetic differentiation among regional populations. Our results highlight the need for further investigations of genetic erosion in Nearctic terrestrial mammals, of which several show similar colonization histories in the High Artic.

Genetic Structure and Diversity Among Historic and Modern Populations of the Sumatran Rhinoceros (Dicerorhinus sumatrensis).

The Sumatran rhinoceros (Dicerorhinus sumatrensis), once widespread across Southeast Asia, now consists of as few as 30 individuals within Sumatra and Borneo. To aid in conservation planning, we sequenced 218 bp of control region mitochondrial (mt) DNA, identifying 17 distinct mitochondrial haplotypes across modern (N = 13) and museum (N = 26) samples. Museum specimens from Laos and Myanmar had divergent mtDNA, consistent with the placement of western mainland rhinos into the distinct subspecies D. s. lasiotis (presumed extinct). Haplotypes from Bornean rhinos were highly diverse, but dissimilar from those of other regions, supporting the distinctiveness of the subspecies D. s. harrissoni. Rhinos from Sumatra and Peninsular Malaysia shared mtDNA haplotypes, consistent with their traditional placement into a single subspecies D. s sumatrensis. Modern samples of D. s. sumatrensis were genotyped at 18 microsatellite loci. Rhinos within Sumatra formed 2 sub-populations, likely separated by the Barisan Mountains, though with only modest genetic differentiation between them. There are so few remaining Sumatran rhinoceros that separate management strategies for subspecies or subpopulations may not be viable, while each surviving rhino pedigree is likely to retain alleles found in no other individuals. Given the low population size and low reproductive potential of Sumatran rhinos, rapid genetic erosion is inevitable, though an under-appreciated concern is the potential for fixation of harmful genetic variants. Both concerns underscore 2 overriding priorities for the species: 1) translocation of wild rhinos to ex situ facilities, and 2) collection and storage of gametes and cell lines from every surviving captive and wild individual.

Genetic connectivity across marginal habitats: the elephants of the Namib Desert.

Locally isolated populations in marginal habitats may be genetically distinctive and of heightened conservation concern. Elephants inhabiting the Namib Desert have been reported to show distinctive behavioral and phenotypic adaptations in that severely arid environment. The genetic distinctiveness of Namibian desert elephants relative to other African savanna elephant (Loxodonta africana) populations has not been established. To investigate the genetic structure of elephants in Namibia, we determined the mitochondrial (mt) DNA control region sequences and genotyped 17 microsatellite loci in desert elephants (n = 8) from the Hoanib River catchment and the Hoarusib River catchment. We compared these to the genotypes of elephants (n = 77) from other localities in Namibia. The mtDNA haplotype sequences and frequencies among desert elephants were similar to those of elephants in Etosha National Park, the Huab River catchment, the Ugab River catchment, and central Kunene, although the geographically distant Caprivi Strip had different mtDNA haplotypes. Likewise, analysis of the microsatellite genotypes of desert-dwelling elephants revealed that they were not genetically distinctive from Etosha elephants, and there was no evidence for isolation by distance across the Etosha region. These results, and a review of the historical record, suggest that a high learning capacity and long-distance migrations allowed Namibian elephants to regularly shift their ranges to survive in the face of high variability in climate and in hunting pressure.

Extensive sampling of polar bears (Ursus maritimus) in the Northwest Passage (Canadian Arctic Archipelago) reveals population differentiation across multiple spatial and temporal scales.

As global warming accelerates the melting of Arctic sea ice, polar bears (Ursus maritimus) must adapt to a rapidly changing landscape. This process will necessarily alter the species distribution together with population dynamics and structure. Detailed knowledge of these changes is crucial to delineating conservation priorities. Here, we sampled 361 polar bears from across the center of the Canadian Arctic Archipelago spanning the Gulf of Boothia (GB) and M'Clintock Channel (MC). We use DNA microsatellites and mitochondrial control region sequences to quantify genetic differentiation, estimate gene flow, and infer population history. Two populations, roughly coincident with GB and MC, are significantly differentiated at both nuclear (F ST = 0.01) and mitochondrial (ΦST = 0.47; F ST = 0.29) loci, allowing Bayesian clustering analyses to assign individuals to either group. Our data imply that the causes of the mitochondrial and nuclear genetic patterns differ. Analysis of mtDNA reveals the matrilineal structure dates at least to the Holocene, and is common to individuals throughout the species' range. These mtDNA differences probably reflect both genetic drift and historical colonization dynamics. In contrast, the differentiation inferred from microsatellites is only on the scale of hundreds of years, possibly reflecting contemporary impediments to gene flow. Taken together, our data suggest that gene flow is insufficient to homogenize the GB and MC populations and support the designation of GB and MC as separate polar bear conservation units. Our study also provide a striking example of how nuclear DNA and mtDNA capture different aspects of a species demographic history.

Distinguishing forest and savanna African elephants using short nuclear DNA sequences.

A more complete description of African elephant phylogeography would require a method that distinguishes forest and savanna elephants using DNA from low-quality samples. Although mitochondrial DNA is often the marker of choice for species identification, the unusual cytonuclear patterns in African elephants make nuclear markers more reliable. We therefore designed and utilized genetic markers for short nuclear DNA regions that contain fixed nucleotide differences between forest and savanna elephants. We used M13 forward and reverse sequences to increase the total length of PCR amplicons and to improve the quality of sequences for the target DNA. We successfully sequenced fragments of nuclear genes from dung samples of known savanna and forest elephants in the Democratic Republic of Congo, Ethiopia, and Namibia. Elephants at previously unexamined locations were found to have nucleotide character states consistent with their status as savanna or forest elephants. Using these and results from previous studies, we estimated that the short-amplicon nuclear markers could distinguish forest from savanna African elephants with more than 99% accuracy. Nuclear genotyping of museum, dung, or ivory samples will provide better-informed conservation management of Africa's elephants.

Maternal influences on population dynamics: evidence from an exploited freshwater fish.

We used a field experiment, population modeling, and an analysis of 30 years of data from walleye (Sander vitreus; a freshwater fish) in Lake Erie to show that maternal influences on offspring survival can affect population dynamics. We first demonstrate experimentally that the survival of juvenile walleye increases with egg size (and, to a lesser degree, female energy reserves). Because egg size in this species tends to increase with maternal age, we then model these maternal influences on offspring survival as a function of maternal age to show that adult age structure can affect the maximum rate at which a population can produce new adults. Consistent with this hypothesis, we present empirical evidence that the maximum reproductive rate of an exploited population of walleye was approximately twice as high when older females were abundant as compared to when they were relatively scarce. Taken together, these results indicate that age- or size-based maternal influences on offspring survival can be an important mechanism driving population dynamics and that exploited populations could benefit from management strategies that protect, rather than target, reproductively valuable individuals.

Appraisal of the consequences of the DDT-induced bottleneck on the level and geographic distribution of neutral genetic variation in Canadian peregrine falcons, Falco peregrinus.

Peregrine falcon populations underwent devastating declines in the mid-20th century due to the bioaccumulation of organochlorine contaminants, becoming essentially extirpated east of the Great Plains and significantly reduced elsewhere in North America. Extensive re-introduction programs and restrictions on pesticide use in Canada and the United States have returned many populations to predecline sizes. A proper population genetic appraisal of the consequences of this decline requires an appropriate context defined by (i) meaningful demographic entities; and (ii) suitable reference populations. Here we explore the validity of currently recognized subspecies designations using data from the mitochondrial control region and 11 polymorphic microsatellite loci taken from 184 contemporary individuals from across the breeding range, and compare patterns of population genetic structure with historical patterns inferred from 95 museum specimens. Of the three North American subspecies, the west coast marine subspecies Falco peregrinus pealei is well differentiated genetically in both time periods using nuclear loci. In contrast, the partitioning of continental Falco peregrinus anatum and arctic Falco peregrinus tundrius subspecies is not substantiated, as individuals from these subspecies are historically indistinguishable genetically. Bayesian clustering analyses demonstrate that contemporary genetic differentiation between these two subspecies is mainly due to changes within F. p. anatum (specifically the southern F. p. anatum populations). Despite expectations and a variety of tests, no genetic bottleneck signature is found in the identified populations; in fact, many contemporary indices of diversity are higher than historical values. These results are rationalized by the promptness of the recovery and the possible introduction of new genetic material.