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.

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.

Targeted conservation genetics of the endangered chimpanzee.

Populations of the common chimpanzee (Pan troglodytes) are in an impending risk of going extinct in the wild as a consequence of damaging anthropogenic impact on their natural habitat and illegal pet and bushmeat trade. Conservation management programmes for the chimpanzee have been established outside their natural range (ex situ), and chimpanzees from these programmes could potentially be used to supplement future conservation initiatives in the wild (in situ). However, these programmes have often suffered from inadequate information about the geographical origin and subspecies ancestry of the founders. Here, we present a newly designed capture array with ~60,000 ancestry informative markers used to infer ancestry of individual chimpanzees in ex situ populations and determine geographical origin of confiscated sanctuary individuals. From a test panel of 167 chimpanzees with unknown origins or subspecies labels, we identify 90 suitable non-admixed individuals in the European Association of Zoos and Aquaria (EAZA) Ex situ Programme (EEP). Equally important, another 46 individuals have been identified with admixed subspecies ancestries, which therefore over time, should be naturally phased out of the breeding populations. With potential for future re-introduction to the wild, we determine the geographical origin of 31 individuals that were confiscated from the illegal trade and demonstrate the promises of using non-invasive sampling in future conservation action plans. Collectively, our genomic approach provides an exemplar for ex situ management of endangered species and offers an efficient tool in future in situ efforts to combat the illegal wildlife trade.

Unrecognized circulation of SAT 1 foot-and-mouth disease virus in cattle herds around Queen Elizabeth National Park in Uganda.

BACKGROUND: Foot-and-mouth disease (FMD) is endemic in Uganda in spite of the control measures used. Various aspects of the maintenance and circulation of FMD viruses (FMDV) in Uganda are not well understood; these include the role of the African buffalo (Syncerus caffer) as a reservoir for FMDV. To better understand the epidemiology of FMD at the livestock-wildlife-interface, samples were collected from young, unvaccinated cattle from 24 pastoral herds that closely interact with wildlife around Queen Elizabeth National Park in Uganda, and analysed for evidence of FMDV infection. RESULTS: In total, 37 (15%) of 247 serum samples had detectable antibodies against FMDV non-structural proteins (NSPs) using a pan-serotypic assay. Within these 37 sera, antibody titres ≥ 80 against the structural proteins of serotypes O, SAT 1, SAT 2 and SAT 3 were detected by ELISA in 5, 7, 4 and 3 samples, respectively, while neutralizing antibodies were only detected against serotype O in 3 samples. Two FMDV isolates, with identical VP1 coding sequences, were obtained from probang samples from clinically healthy calves from the same herd and are serotype SAT 1 (topotype IV (EA-I)). Based on the VP1 coding sequences, these viruses are distinct from previous cattle and buffalo SAT 1 FMDV isolates obtained from the same area (19-30% nucleotide difference) and from the vaccine strain (TAN/155/71) used within Uganda (26% nucleotide difference). Eight herds had only one or a few animals with antibodies against FMDV NSPs while six herds had more substantial evidence of prior infection with FMDV. There was no evidence for exposure to FMDV in the other ten herds. CONCLUSIONS: The two identical SAT 1 FMDV VP1 sequences are distinct from former buffalo and cattle isolates from the same area, thus, transmission between buffalo and cattle was not demonstrated. These new SAT 1 FMDV isolates differed significantly from the vaccine strain used to control Ugandan FMD outbreaks, indicating a need for vaccine matching studies. Only six herds had clear serological evidence for exposure to O and SAT 1 FMDV. Scattered presence of antibodies against FMDV in other herds may be due to the occasional introduction of animals to the area or maternal antibodies from past infection and/or vaccination. The evidence for asymptomatic FMDV infection has implications for disease control strategies in the area since this obstructs early disease detection that is based on clinical signs in FMDV infected animals.

Foot-and-mouth disease virus serotype SAT 3 in long-horned Ankole calf, Uganda.

After a 16-year interval, foot-and-mouth disease virus serotype SAT 3 was isolated in 2013 from an apparently healthy long-horned Ankole calf that grazed close to buffalo in Uganda. The emergent virus strain is ≈20% different in nucleotide sequence (encoding VP1 [viral protein 1]) from its closest relatives isolated previously from buffalo in Uganda.

Characterisation of recent foot-and-mouth disease viruses from African buffalo (Syncerus caffer) and cattle in Kenya is consistent with independent virus populations.

BACKGROUND: Understanding the epidemiology of foot-and-mouth disease (FMD), including roles played by different hosts, is essential for improving disease control. The African buffalo (Syncerus caffer) is a reservoir for the SAT serotypes of FMD virus (FMDV). Large buffalo populations commonly intermingle with livestock in Kenya, yet earlier studies have focused on FMD in the domestic livestock, hence the contribution of buffalo to disease in livestock is largely unknown. This study analysed 47 epithelia collected from FMD outbreaks in Kenyan cattle between 2008 and 2012, and 102 probang and serum samples collected from buffalo in three different Kenyan ecosystems; Maasai-Mara (MME) (n = 40), Tsavo (TSE) (n = 33), and Meru (ME) (n = 29). RESULTS: Antibodies against FMDV non-structural proteins were found in 65 of 102 (64%) sera from buffalo with 44/102 and 53/102 also having neutralising antibodies directed against FMDV SAT 1 and SAT 2, respectively. FMDV RNA was detected in 42% of the buffalo probang samples by RT-qPCR (Cycle Threshold (Ct) ≤32). Two buffalo probang samples were positive by VI and were identified as FMDV SAT 1 and SAT 2 by Ag-ELISA, while the latter assay detected serotypes O (1), A (20), SAT 1 (7) and SAT 2 (19) in the 47 cattle epithelia. VP1 coding sequences were generated for two buffalo and 21 cattle samples. Phylogenetic analyses revealed SAT 1 and SAT 2 virus lineages within buffalo that were distinct from those detected in cattle. CONCLUSIONS: We found that FMDV serotypes O, A, SAT 1 and SAT 2 were circulating among cattle in Kenya and cause disease, but only SAT 1 and SAT 2 viruses were successfully isolated from clinically normal buffalo. The buffalo isolates were genetically distinct from isolates obtained from cattle. Control efforts should focus primarily on reducing FMDV circulation among livestock and limiting interaction with buffalo. Comprehensive studies incorporating additional buffalo viruses are recommended.

Large-scale ruminant genome sequencing provides insights into their evolution and distinct traits.

The ruminants are one of the most successful mammalian lineages, exhibiting morphological and habitat diversity and containing several key livestock species. To better understand their evolution, we generated and analyzed de novo assembled genomes of 44 ruminant species, representing all six Ruminantia families. We used these genomes to create a time-calibrated phylogeny to resolve topological controversies, overcoming the challenges of incomplete lineage sorting. Population dynamic analyses show that population declines commenced between 100,000 and 50,000 years ago, which is concomitant with expansion in human populations. We also reveal genes and regulatory elements that possibly contribute to the evolution of the digestive system, cranial appendages, immune system, metabolism, body size, cursorial locomotion, and dentition of the ruminants.

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.

The confounding effect of population structure on Bayesian skyline plot inferences of demographic history.

Many coalescent-based methods aiming to infer the demographic history of populations assume a single, isolated and panmictic population (i.e. a Wright-Fisher model). While this assumption may be reasonable under many conditions, several recent studies have shown that the results can be misleading when it is violated. Among the most widely applied demographic inference methods are Bayesian skyline plots (BSPs), which are used across a range of biological fields. Violations of the panmixia assumption are to be expected in many biological systems, but the consequences for skyline plot inferences have so far not been addressed and quantified. We simulated DNA sequence data under a variety of scenarios involving structured populations with variable levels of gene flow and analysed them using BSPs as implemented in the software package BEAST. Results revealed that BSPs can show false signals of population decline under biologically plausible combinations of population structure and sampling strategy, suggesting that the interpretation of several previous studies may need to be re-evaluated. We found that a balanced sampling strategy whereby samples are distributed on several populations provides the best scheme for inferring demographic change over a typical time scale. Analyses of data from a structured African buffalo population demonstrate how BSP results can be strengthened by simulations. We recommend that sample selection should be carefully considered in relation to population structure previous to BSP analyses, and that alternative scenarios should be evaluated when interpreting signals of population size change.