Genetic structure of fragmented southern populations of African Cape buffalo (Syncerus caffer caffer)

Continent-wide genomic analysis of the African buffalo (Syncerus caffer)

The African buffalo (Syncerus caffer) is a wild bovid with a historical distribution across much of sub-Saharan Africa. Genomic analysis can provide insights into the evolutionary history of the species, and the key selective pressures shaping populations, including assessment of population level differentiation, population fragmentation, and population genetic structure. In this study we generated the highest quality de novo genome assembly (2.65 Gb, scaffold N50 69.17 Mb) of African buffalo to date, and sequenced a further 195 genomes from across the species distribution. Principal component and admixture analyses provided little support for the currently described four subspecies. Estimating Effective Migration Surfaces analysis suggested that geographical barriers have played a significant role in shaping gene flow and the population structure. Estimated effective population sizes indicated a substantial drop occurring in all populations 5-10,000 years ago, coinciding with the increase in human populations. Finally, signatures of selection were enriched for key genes associated with the immune response, suggesting infectious disease exert a substantial selective pressure upon the African buffalo. These findings have important implications for understanding bovid evolution, buffalo conservation and population management.

Genome-wide diversity, population structure and signatures of inbreeding in the African buffalo in Mozambique

The African buffalo, Syncerus caffer, is a key species in African ecosystems. Like other large herbivores, it plays a fundamental role in its habitat acting as an ecosystem engineer. Over the last few centuries, African buffalo populations have declined because of range contraction and demographic decline caused by direct or indirect human activities. In Mozambique, historically home to large buffalo herds, the combined effect of colonialism and subsequent civil wars has created a critical situation that urgently needs to be addressed. In this study, we focused on the analysis of genetic diversity of Syncerus caffer caffer populations from six areas of Mozambique. Using genome-wide SNPs obtained from ddRAD sequencing, we examined the population structure across the country, estimated gene flow between areas under conservation management, including national reserves, and assessed the inbreeding coefficients. Our results indicate that all studied populations of Syncerus caffer caffer are genetically depauperate, with a high level of inbreeding. Moreover, buffaloes in Mozambique present a significant population differentiation between southern and central areas. We found an unexpected genotype in the Gorongosa National Park, where buffaloes experienced a dramatic population size reduction, that shares a common ancestry with southern populations of Catuane and Namaacha. This could suggest the past occurrence of a connection between southern and central Mozambique and that the observed population structuring could reflect recent events of anthropogenic origin. All the populations analysed showed high levels of homozygosity, likely due to extensive inbreeding over the last few decades, which could have increased the frequency of recessive deleterious alleles. Improving the resilience of Syncerus caffer caffer in Mozambique is essential for preserving the ecosystem integrity. The most viable approach appears to be facilitating translocations and re-establishing connectivity between isolated herds. However, our results also highlight the importance of assessing intraspecific genetic diversity when considering interventions aimed at enhancing population viability such as selecting suitable source populations.

A continent-wide high genetic load in African buffalo revealed by clines in the frequency of deleterious alleles, genetic hitchhiking and linkage disequilibrium

A high genetic load can negatively affect population viability and increase susceptibility to diseases and other environmental stressors. Prior microsatellite studies of two African buffalo (Syncerus caffer) populations in South Africa indicated substantial genome-wide genetic load due to high-frequency occurrence of deleterious alleles. The occurrence of these alleles, which negatively affect male body condition and bovine tuberculosis resistance, throughout most of the buffalo's range were evaluated in this study. Using available microsatellite data (2-17 microsatellite loci) for 1676 animals from 34 localities (from 25°S to 5°N), we uncovered continent-wide frequency clines of microsatellite alleles associated with the aforementioned male traits. Frequencies decreased over a south-to-north latitude range (average per-locus Pearson r = -0.22). The frequency clines coincided with a multilocus-heterozygosity cline (adjusted R2 = 0.84), showing up to a 16% decrease in southern Africa compared to East Africa. Furthermore, continent-wide linkage disequilibrium (LD) at five linked locus pairs was detected, characterized by a high fraction of positive interlocus associations (0.66, 95% CI: 0.53, 0.77) between male-deleterious-trait-associated alleles. Our findings suggest continent-wide and genome-wide selection of male-deleterious alleles driven by an earlier observed sex-chromosomal meiotic drive system, resulting in frequency clines, reduced heterozygosity due to hitchhiking effects and extensive LD due to male-deleterious alleles co-occurring in haplotypes. The selection pressures involved must be high to prevent destruction of allele-frequency clines and haplotypes by LD decay. Since most buffalo populations are stable, these results indicate that natural mammal populations, depending on their genetic background, can withstand a high genetic load.

Discovering of Genomic Variations Associated to Growth Traits by GWAS in Braunvieh Cattle

A genome-wide association study (GWAS) was performed to elucidate genetic architecture of growth traits in Braunvieh cattle. Methods: The study included 300 genotyped animals by the GeneSeek^® Genomic Profiler Bovine LDv.4 panel; after quality control, 22,734 SNP and 276 animals were maintained in the analysis. The examined phenotypic data considered birth (BW), weaning (WW), and yearling weights. The association analysis was performed using the principal components method via the egscore function of the GenABEL version 1.8-0 package in the R environment. The marker rs133262280 located in BTA 22 was associated with BW, and two SNPs were associated with WW, rs43668789 (BTA 11) and rs136155567 (BTA 27). New QTL associated with these liveweight traits and four positional and functional candidate genes potentially involved in variations of the analyzed traits were identified. The most important genes in these genomic regions were MCM2 (minichromosome maintenance complex component 2), TPRA1 (transmembrane protein adipocyte associated 1), GALM (galactose mutarotase), and NRG1 (neuregulin 1), related to embryonic cleavage, bone and tissue growth, cell adhesion, and organic development. This study is the first to present a GWAS conducted in Braunvieh cattle in Mexico providing evidence for genetic architecture of assessed growth traits. Further specific analysis of found associated genes and regions will clarify its contribution to the genetic basis of growth-related traits.

Host Genetic Diversity and Infectious Diseases. Focus on Wild Boar, Red Deer and Tuberculosis

Host genetic diversity tends to limit disease spread in nature and buffers populations against epidemics. Genetic diversity in wildlife is expected to receive increasing attention in contexts related to disease transmission and human health. Ungulates such as wild boar (Sus scrofa) and red deer (Cervus elaphus) are important zoonotic hosts that can be precursors to disease emergence and spread in humans. Tuberculosis is a zoonotic disease with relevant consequences and can present high prevalence in wild boar and red deer populations. Here, we review studies on the genetic diversity of ungulates and determine to what extent these studies consider its importance on the spread of disease. This assessment also focused on wild boar, red deer, and tuberculosis. We found a disconnection between studies treating genetic diversity and those dealing with infectious diseases. Contrarily, genetic diversity studies in ungulates are mainly concerned with conservation. Despite the existing disconnection between studies on genetic diversity and studies on disease emergence and spread, the knowledge gathered in each discipline can be applied to the other. The bidirectional applications are illustrated in wild boar and red deer populations from Spain, where TB is an important threat for wildlife, livestock, and humans.

High diversity, inbreeding and a dynamic Pleistocene demographic history revealed by African buffalo genomes

Genomes retain records of demographic changes and evolutionary forces that shape species and populations. Remnant populations of African buffalo (Syncerus caffer) in South Africa, with varied histories, provide an opportunity to investigate signatures left in their genomes by past events, both recent and ancient. Here, we produce 40 low coverage (7.14×) genome sequences of Cape buffalo (S. c. caffer) from four protected areas in South Africa. Genome-wide heterozygosity was the highest for any mammal for which these data are available, while differences in individual inbreeding coefficients reflected the severity of historical bottlenecks and current census sizes in each population. PSMC analysis revealed multiple changes in N_e between approximately one million and 20 thousand years ago, corresponding to paleoclimatic changes and Cape buffalo colonisation of southern Africa. The results of this study have implications for buffalo management and conservation, particularly in the context of the predicted increase in aridity and temperature in southern Africa over the next century as a result of climate change.

Genetic diversity, relatedness and inbreeding of ranched and fragmented Cape buffalo populations in southern Africa

Wildlife ranching, although not considered a conventional conservation system, provides a sustainable model for wildlife utilization and could be a source of valuable genetic material. However, increased fragmentation and intensive management may threaten the evolutionary potential and conservation value of species. Disease-free Cape buffalo (Syncerus caffer caffer) in southern Africa exist in populations with a variety of histories and management practices. We compared the genetic diversity of buffalo in national parks to private ranches and found that, except for Addo Elephant National Park, genetic diversity was high and statistically equivalent. We found that relatedness and inbreeding levels were not substantially different between ranched populations and those in national parks, indicating that breeding practices likely did not yet influence genetic diversity of buffalo on private ranches in this study. High genetic differentiation between South African protected areas highlighted their fragmented nature. Structure analysis revealed private ranches comprised three gene pools, with origins from Addo Elephant National Park, Kruger National Park and a third, unsampled gene pool. Based on these results, we recommend the Addo population be supplemented with disease-free Graspan and Mokala buffalo (of Kruger origin). We highlight the need for more research to characterize the genetic diversity and composition of ranched wildlife species, in conjunction with wildlife ranchers and conservation authorities, in order to evaluate the implications for management and conservation of these species across different systems.

Spatial and seasonal patterns of FMD primary outbreaks in cattle in Zimbabwe between 1931 and 2016

Foot and mouth disease (FMD) is an important livestock disease impacting mainly intensive production systems. In southern Africa, the FMD virus is maintained in wildlife and its control is therefore complicated. However, FMD control is an important task to allow countries access to lucrative foreign meat market and veterinary services implement drastic control measures on livestock populations living in the periphery of protected areas, negatively impacting local small-scale livestock producers. This study investigated FMD primary outbreak data in Zimbabwe from 1931 to 2016 to describe the spatio-temporal distribution of FMD outbreaks and their potential drivers. The results suggest that: (i) FMD outbreaks were not randomly distributed in space across Zimbabwe but are clustered in the Southeast Lowveld (SEL); (ii) the proximity of protected areas with African buffalos was potentially responsible for primary FMD outbreaks in cattle; (iii) rainfall per se was not associated with FMD outbreaks, but seasons impacted the temporal occurrence of FMD outbreaks across regions; (iv) the frequency of FMD outbreaks increased during periods of major socio-economic and political crisis. The differences between the spatial clusters and other areas in Zimbabwe presenting similar buffalo/cattle interfaces but with fewer FMD outbreaks can be interpreted in light of the recent better understanding of wildlife/livestock interactions in these areas. The types of wildlife/livestock interfaces are hypothesized to be the key drivers of contacts between wildlife and livestock, triggering a risk of FMD inter-species spillover. The management of wildlife/livestock interfaces is therefore crucial for the control of FMD in southern Africa.

A natural gene drive system influences bovine tuberculosis susceptibility in African buffalo: Possible implications for disease management

Bovine tuberculosis (BTB) is endemic to the African buffalo (Syncerus caffer) of Hluhluwe-iMfolozi Park (HiP) and Kruger National Park, South Africa. In HiP, the disease has been actively managed since 1999 through a test-and-cull procedure targeting BTB-positive buffalo. Prior studies in Kruger showed associations between microsatellite alleles, BTB and body condition. A sex chromosomal meiotic drive, a form of natural gene drive, was hypothesized to be ultimately responsible. These associations indicate high-frequency occurrence of two types of male-deleterious alleles (or multiple-allele haplotypes). One type negatively affects body condition and BTB resistance in both sexes. The other type has sexually antagonistic effects: negative in males but positive in females. Here, we investigate whether a similar gene drive system is present in HiP buffalo, using 17 autosomal microsatellites and microsatellite-derived Y-chromosomal haplotypes from 401 individuals, culled in 2002–2004. We show that the association between autosomal microsatellite alleles and BTB susceptibility detected in Kruger, is also present in HiP. Further, Y-haplotype frequency dynamics indicated that a sex chromosomal meiotic drive also occurred in HiP. BTB was associated with negative selection of male-deleterious alleles in HiP, unlike positive selection in Kruger. Birth sex ratios were female-biased. We attribute negative selection and female-biased sex ratios in HiP to the absence of a Y-chromosomal sex-ratio distorter. This distorter has been hypothesized to contribute to positive selection of male-deleterious alleles and male-biased birth sex ratios in Kruger. As previously shown in Kruger, microsatellite alleles were only associated with male-deleterious effects in individuals born after wet pre-birth years; a phenomenon attributed to epigenetic modification. We identified two additional allele types: male-specific deleterious and beneficial alleles, with no discernible effect on females. Finally, we discuss how our findings may be used for breeding disease-free buffalo and implementing BTB test-and-cull programs.

Genetic responsiveness of African buffalo to environmental stressors: A role for epigenetics in balancing autosomal and sex chromosome interactions?

In the African buffalo (Syncerus caffer) population of the Kruger National Park (South Africa) a primary sex-ratio distorter and a primary sex-ratio suppressor have been shown to occur on the Y chromosome. A subsequent autosomal microsatellite study indicated that two types of deleterious alleles with a negative effect on male body condition, but a positive effect on relative fitness when averaged across sexes and generations, occur genome-wide and at high frequencies in the same population. One type negatively affects body condition of both sexes, while the other acts antagonistically: it negatively affects male but positively affects female body condition. Here we show that high frequencies of male-deleterious alleles are attributable to Y-chromosomal distorter-suppressor pair activity and that these alleles are suppressed in individuals born after three dry pre-birth years, likely through epigenetic modification. Epigenetic suppression was indicated by statistical interactions between pre-birth rainfall, a proxy for parental body condition, and the phenotypic effect of homozygosity/heterozygosity status of microsatellites linked to male-deleterious alleles, while a role for the Y-chromosomal distorter-suppressor pair was indicated by between-sex genetic differences among pre-dispersal calves. We argue that suppression of male-deleterious alleles results in negative frequency-dependent selection of the Y distorter and suppressor; a prerequisite for a stable polymorphism of the Y distorter-suppressor pair. The Y distorter seems to be responsible for positive selection of male-deleterious alleles during resource-rich periods and the Y suppressor for positive selection of these alleles during resource-poor periods. Male-deleterious alleles were also associated with susceptibility to bovine tuberculosis, indicating that Kruger buffalo are sensitive to stressors such as diseases and droughts. We anticipate that future genetic studies on African buffalo will provide important new insights into gene fitness and epigenetic modification in the context of sex-ratio distortion and infectious disease dynamics.

Ancient diversity and geographical sub-structuring in African buffalo Theileria parva populations revealed through metagenetic analysis of antigen-encoding loci

An infection and treatment protocol involving infection with a mixture of three parasite isolates and simultaneous treatment with oxytetracycline is currently used to vaccinate cattle against Theileria parva. While vaccination results in high levels of protection in some regions, little or no protection is observed in areas where animals are challenged predominantly by parasites of buffalo origin. A previous study involving sequencing of two antigen-encoding genes from a series of parasite isolates indicated that this is associated with greater antigenic diversity in buffalo-derived T. parva. The current study set out to extend these analyses by applying high-throughput sequencing to ex vivo samples from naturally infected buffalo to determine the extent of diversity in a set of antigen-encoding genes. Samples from two populations of buffalo, one in Kenya and the other in South Africa, were examined to investigate the effect of geographical distance on the nature of sequence diversity. The results revealed a number of significant findings. First, there was a variable degree of nucleotide sequence diversity in all gene segments examined, with the percentage of polymorphic nucleotides ranging from 10% to 69%. Second, large numbers of allelic variants of each gene were found in individual animals, indicating multiple infection events. Third, despite the observed diversity in nucleotide sequences, several of the gene products had highly conserved amino acid sequences, and thus represent potential candidates for vaccine development. Fourth, although compelling evidence for population differentiation between the Kenyan and South African T. parva parasites was identified, analysis of molecular variance for each gene revealed that the majority of the underlying nucleotide sequence polymorphism was common to both areas, indicating that much of this aspect of genetic variation in the parasite population arose prior to geographic separation.

Transmission of Foot-and-Mouth Disease SAT2 Viruses at the Wildlife-Livestock Interface of Two Major Transfrontier Conservation Areas in Southern Africa

Over a decade ago, foot-and-mouth disease (FMD) re-emerged in Southern Africa specifically in beef exporting countries that had successfully maintained disease-free areas in the past. FMD virus (FMDV) serotype SAT2 has been responsible for a majority of these outbreaks. Epidemiological studies have revealed the importance of the African buffalo as the major wildlife FMD reservoir in the region. We used phylogeographic analysis to study dynamics of FMD transmission between buffalo and domestic cattle at the interface of the major wildlife protected areas in the region currently encompassing two largest Transfrontier conservation areas: Kavango-Zambezi (KAZA) and Great Limpopo (GL). Results of this study showed restricted local occurrence of each FMDV SAT2 topotypes I, II, and III, with occasional virus migration from KAZA to GL. Origins of outbreaks in livestock are frequently attributed to wild buffalo, but our results suggest that transmission from cattle to buffalo also occurs. We used coalescent Bayesian skyline analysis to study the genetic variation of the virus in cattle and buffalo, and discussed the association of these genetic changes in the virus and relevant epidemiological events that occurred in this area. Our results show that the genetic diversity of FMDV SAT2 has decreased in buffalo and cattle population during the last decade. This study contributes to understand the major dynamics of transmission and genetic variation of FMDV SAT2 in Southern Africa, which will could ultimately help in designing efficient strategies for the control of FMD at a local and regional level.