Conservation management strategy impacts inbreeding and mutation load in scimitar-horned oryx.

In an age of habitat loss and overexploitation, small populations, both captive and wild, are increasingly facing the effects of isolation and inbreeding. Genetic management has therefore become a vital tool for ensuring population viability. However, little is known about how the type and intensity of intervention shape the genomic landscape of inbreeding and mutation load. We address this using whole-genome sequence data of the scimitar-horned oryx (Oryx dammah), an iconic antelope that has been subject to contrasting management strategies since it was declared extinct in the wild. We show that unmanaged populations are enriched for long runs of homozygosity (ROH) and have significantly higher inbreeding coefficients than managed populations. Additionally, despite the total number of deleterious alleles being similar across management strategies, the burden of homozygous deleterious genotypes was consistently higher in unmanaged groups. These findings emphasize the risks associated with deleterious mutations through multiple generations of inbreeding. As wildlife management strategies continue to diversify, our study reinforces the importance of maintaining genome-wide variation in vulnerable populations and has direct implications for one of the largest reintroduction attempts in the world.

Chromosome-level genome provides insight into the evolution and conservation of the threatened goral (Naemorhedus goral).

BACKGROUND: Gorals Naemorhedus resemble both goats and antelopes, which prompts much debate about the intragenus species delimitation and phylogenetic status of the genus Naemorhedus within the subfamily Caprinae. Their evolution is believed to be linked to the uplift of the Qinghai-Tibet Plateau (QTP). To better understand its phylogenetics, the genetic information is worth being resolved. RESULTS: Based on a sample from the eastern margin of QTP, we constructed the first reference genome for Himalayan goral Naemorhedus goral, using PacBio long-read sequencing and Hi-C technology. The 2.59 Gb assembled genome had a contig N50 of 3.70 Mb and scaffold N50 of 106.66 Mb, which anchored onto 28 pseudo chromosomes. A total of 20,145 protein-coding genes were predicted in the assembled genome, of which 99.93% were functionally annotated. Phylogenetically, the goral was closely related to muskox on the mitochondrial genome level and nested into the takin-muskox clade on the genome tree, rather than other so-called goat-antelopes. The cladogenetic event among muskox, takin and goral occurred sequentially during the late Miocene (~ 11 - 5 Mya), when the QTP experienced a third dramatic uplift with consequent profound changes in climate and environment. Several chromosome fusions and translocations were observed between goral and takin/muskox. The expanded gene families in the goral genome were mainly related to the metabolism of drugs and diseases, so as the positive selected genes. The Ne of goral continued to decrease since ~ 1 Mya during the Pleistocene with active glaciations. CONCLUSION: The high-quality goral genome provides insights into the evolution and valuable information for the conservation of this threatened group.

Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini).

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

Genetic diversity and structure of mongolian gazelle (Procapra gutturosa) populations in fragmented habitats.

BACKGROUND: The Mongolian gazelle (Procapra gutturosa) population has shown a considerable range of contractions and local extinctions over the last century, owing to habitat fragmentation and poaching. A thorough understanding of the genetic diversity and structure of Mongolian gazelle populations in fragmented habitats is critical for planning effective conservation strategies. RESULT: In this study, we used eight microsatellite loci and mitochondrial cytochrome b (Cytb) to compare the levels of genetic diversity and genetic structure of Mongolian gazelle populations in the Hulun Lake National Nature Reserve (HLH) with those in the China-Mongolia border area (BJ). The results showed that the nucleotide diversity and observed heterozygosity of the HLH population were lower than those of the BJ population. Moreover, the HLH and BJ populations showed genetic differentiation. We concluded that the HLH population had lower genetic diversity and a distinct genetic structure compared with the BJ population. CONCLUSION: The genetic diversity of fragmented Mongolian gazelle populations, can be improved by protecting these populations while reinforcing their gene exchange with other populations. For example, attempts can be made to introduce new individuals with higher genetic diversity from other populations to reduce inbreeding.

Blue Turns to Gray: Paleogenomic Insights into the Evolutionary History and Extinction of the Blue Antelope (Hippotragus leucophaeus).

The blue antelope (Hippotragus leucophaeus) is the only large African mammal species to have become extinct in historical times, yet no nuclear genomic information is available for this species. A recent study showed that many alleged blue antelope museum specimens are either roan (Hippotragus equinus) or sable (Hippotragus niger) antelopes, further reducing the possibilities for obtaining genomic information for this extinct species. While the blue antelope has a rich fossil record from South Africa, climatic conditions in the region are generally unfavorable to the preservation of ancient DNA. Nevertheless, we recovered two blue antelope draft genomes, one at 3.4× mean coverage from a historical specimen (∼200 years old) and one at 2.1× mean coverage from a fossil specimen dating to 9,800-9,300 cal years BP, making it currently the oldest paleogenome from Africa. Phylogenomic analyses show that blue and sable antelope are sister species, confirming previous mitogenomic results, and demonstrate ancient gene flow from roan into blue antelope. We show that blue antelope genomic diversity was much lower than in roan and sable antelope, indicative of a low population size since at least the early Holocene. This supports observations from the fossil record documenting major decreases in the abundance of blue antelope after the Pleistocene-Holocene transition. Finally, the persistence of this species throughout the Holocene despite low population size suggests that colonial-era human impact was likely the decisive factor in the blue antelope's extinction.

Convergent Genomic Signatures of Cashmere Traits: Evidence for Natural and Artificial Selection.

Convergent evolution provides powerful opportunities to investigate the genetic basis of complex traits. The Tibetan antelope (Pantholops hodgsonii) and Siberian ibex (Capra sibirica) belong to different subfamilies in Bovidae, but both have evolved similar superfine cashmere characteristics to meet the cold temperature in plateau environments. The cashmere traits of cashmere goats underwent strong artificial selection, and some traces of domestication also remained in the genome. Hence, we investigated the convergent genomic signatures of cashmere traits between natural and artificial selection. We compared the patterns of convergent molecular evolution between Tibetan antelope and Siberian ibex by testing positively selected genes, rapidly evolving genes and convergent amino acid substitutions. In addition, we analyzed the selected genomic features of cashmere goats under artificial selection using whole-genome resequencing data, and skin transcriptome data of cashmere goats were also used to focus on the genes involved in regulating cashmere traits. We found that molecular convergent events were very rare, but natural and artificial selection genes were convergent enriched in similar functional pathways (e.g., ECM-receptor interaction, focal adhesion, PI3K-Akt signaling pathway) in a variety of gene sets. Type IV collagen family genes (COL4A2, COL4A4, COL4A5, COL6A5, COL6A6) and integrin family genes (ITGA2, ITGA4, ITGA9, ITGB8) may be important candidate genes for cashmere formation and development. Our results provide a comprehensive approach and perspective for exploring cashmere traits and offer a valuable reference for subsequent in-depth research on the molecular mechanisms regulating cashmere development and fineness.

African climate and geomorphology drive evolution and ghost introgression in sable antelope.

The evolutionary history of African ungulates has been explained largely in the light of Pleistocene climatic oscillations and the way these influenced the distribution of vegetation types, leading to range expansions and/or isolation in refugia. In contrast, comparatively fewer studies have addressed the continent's environmental heterogeneity and the role played by its geomorphological barriers. In this study, we performed a range-wide analysis of complete mitogenomes of sable antelope (Hippotragus niger) to explore how these different factors may have contributed as drivers of evolution in southcentral Africa. Our results supported two sympatric and deeply divergent mitochondrial lineages in west Tanzanian sables, which can be explained as the result of introgressive hybridization of a mitochondrial ghost lineage from an archaic, as-yet-undefined, congener. Phylogeographical subdivisions into three main lineages suggest that sable diversification may not have been driven solely by climatic events affecting populations differently across a continental scale. Often in interplay with climate, geomorphological features have also clearly shaped the species' patterns of vicariance, where the East Africa Rift System and the Eastern Arc Mountains acted as geological barriers. Subsequent splits among southern populations may be linked to rearrangements in the Zambezi system, possibly framing the most recent time when the river attained its current drainage profile. This work underlines how the use of comprehensive mitogenomic data sets on a model species with a wide geographical distribution can contribute to a much-enhanced understanding of environmental, geomorphological and evolutionary patterns in Africa throughout the Quaternary.

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.

Cross-amplification of ungulate microsatellite markers in the endemic Indian antelope or blackbuck (Antilope cervicapra) for population monitoring and conservation genetics studies in south Asia.

The Indian antelope or blackbuck (Antilope cervicapra) is endemic to the Indian subcontinent, inhabiting scrublands and dry grasslands. Most of the blackbuck populations are small, isolated, and threatened by habitat fragmentation and degradation. Management of such disjunct populations requires genetic characterization, which is critical for assessing hazards of stochastic events and inbreeding. Addressing the scarcity of such information on the blackbuck, we describe a novel panel of microsatellite markers that could be used to monitor blackbuck demography and population genetic parameters using non-invasive faecal sampling. We screened microsatellites (n = 40) that had been reported to amplify in bovid and cervid species using faecal samples of the blackbuck collected from Kaimoor Wildlife Sanctuary, Uttar Pradesh, India and its vicinities. We selected 12 markers for amplification using faecal DNA extracts (n = 140) in three multiplex reactions. We observed a mean amplification success rate of 72.4% across loci (92.1-25.7%) with high allele diversity (mean number of alleles/locus = 8.67 ± 1.03). Mean genotyping error rates across the markers were low to moderate (allelic drop-out rate = 0.09; false allele rate = 0.11). The proportions of first- and second-order relatives in the study population were 0.69% and 6.21%, respectively. Based on amplification success, genotyping error rates and the probability of identity (PID), we suggest (i) a panel of five microsatellite markers (cumulative PID = 1.24 × 10-5) for individual identification and population monitoring and (ii) seven additional markers for conservation genetics studies. This study provides essential tools capable of augmenting blackbuck conservation strategies at the landscape level, integral to protecting the scrubland-grassland ecosystem.

Genetic diversity and phylogenetic analysis of blackbuck (Antilope cervicapra) in southern India.

Blackbuck (Antilope cervicapra) is a threatened species endemic to the Indian subcontinent. Many populations of blackbuck are found in southern India. Populations of blackbuck are negatively affected in many places for various reasons, such as habitat destruction and poaching. Their range decreased sharply during the 20th century. There is very limited information available on the population dynamics of blackbuck in southern India. For the phylogenetic and genetic diversity analyses of blackbuck populations among different distribution ranges in southern India, we sequenced mt DNA of cytochrome b (Cyt b) for 120, cytochrome c oxidase subunit-1 (COI) for 137 and the control region (CR) for 137 fecal pellets from eleven different locations in southern India. We analyzed the genetic structure of three mitochondrial markers, the CR, Cyt b and the COI region, separately and in a combined dataset. The haplotype diversity and nucleotide diversity of CR were 0.969 and 0.047, respectively, and were higher than those of Cyt b and COI. A Bayesian phylogeny and an MJ network based on the CR and combined dataset (105 sequences) signified several distinct haplotype clusters within blackbuck, whereas no clusters were identified with the Cyt b and COI phylogenetic analyses. The analysis of molecular variance of the combined data set revealed 52.46% genetic variation within the population. Mismatch distribution analysis revealed that blackbuck populations underwent complex changes with analysis of the combined dataset in each population and analysis of each marker separately in the overall population. The results provide evidence that blackbuck in different geographic locations has a distinct population structure due to habitat fragmentation after the formation of the Western and Eastern Ghats.

Shift of Maternal Gut Microbiota of Tibetan Antelope (Pantholops hodgsonii) During the Periparturition Period.

The maternal gut microbiota can influence and be affected by the substantial physiological changes taking place during the periparturition period. However, little information is known about the changes in the maternal gut microbiota and hormonal variations during this period in nonmodel organisms. Tibetan antelope (Pantholops hodgsonii) provide a unique system to address this issue because their summer migration cycle is synchronized with the periparturition period. Here, we used fecal microbiota as a proxy of gut microbiota. We characterized fecal microbial community of female migratory Tibetan antelope in the late pregnancy and postpartum periods using 16S rRNA gene sequencing and quantified fecal glucocorticoids (GCs) and triiodothyronine (T3) metabolite concentrations through enzyme immunoassays to identify the associations between maternal gut microbiota and physiological changes related with reproduction. We found that the fecal microbiota of Tibetan antelope was dominated by Firmicutes and Bacteroidetes. The microbial composition was significantly altered during the transition from late pregnancy to the postpartum period. Fecal T3 concentration was significantly higher in the postpartum period compared to late pregnancy, whereas GC metabolite concentration did not significantly differ between two reproductive states. We identified six genera (Anaerofustis, Bacteroides, Coprococcus_2, Ruminiclostridium_5, Ruminococcaceae_UCG-007, and Tyzzerella) that were significantly associated with reproductive states. We also found two genera (Christensenellaceae_R-7_group and Rikenellaceae_RC9_gut_group) significantly associated with GC metabolite concentration and two genera (Agathobacter and Papillibacter) significantly associated with T3 metabolite concentration, though these correlations were weak with coefficient values ranging from - 0.007 to 0.03. Our results indicate that many members of the gut microbiota are associated with the physiological changes in the transition from late pregnancy to the postpartum period, likely reflecting the metabolic and immune system dynamics during the periparturition period. This study highlights the importance of integrating microbiota, hormones and migration pattern to study the reproductive health of wildlife. By establishing a baseline of the physiological changes during the migration/periparturition period, we can have a better understanding of the impacts of increasing human activities on the Tibetan Plateau on the reproductive health of Tibetan antelope.

The development of genome-wide single nucleotide polymorphisms in blue wildebeest using the DArTseq platform.

Blue wildebeest (Connochaetes taurinus taurinus) are economically important antelope that are widely utilised in the South African wildlife industry. However, very few genomic resources are available for blue wildebeest that can assist in breeding management and facilitate research. This study aimed to develop a set of genome-wide single nucleotide polymorphism (SNP) markers for blue wildebeest. The DArTseq genotyping platform, commonly used in polyploid plant species, was selected for SNP discovery. A limited number of published articles have described the use of the DArTseq platform in animals and, therefore, this study also provided a unique opportunity to assess the performance of the DArTseq platform in an animal species. A total of 20,563 SNPs, each located within a 69 bp sequence, were generated. The developed SNP markers had a high average scoring reproducibility (>99%) and a low percentage missing data (~9.21%) compared to other reduced representation sequencing approaches that have been used in animal studies. Furthermore, the number of candidate SNPs per nucleotide position decreased towards the 3' end of sequence reads, and the ratio of transitions (Ts) to transversions (Tv) remained similar for each read position. These observations indicate that there was no read position bias, such as the identification of false SNPs due to low sequencing quality, towards the tail-end of sequencing reads. The DArTseq platform was also successful in identifying a large number of informative SNPs with desirable polymorphism parameters such as a high minor allele frequency (MAF). The Bos taurus genome was used for the in silico mapping of the marker sequences and a total of 6020 (29.28%) sequences were successfully mapped against the bovine genome. The marker sequences mapped to all of the bovine chromosomes establishing the genome-wide distribution of the SNPs. Moreover, the high observed Ts:Tv ratio (2.84:1) indicate that the DArTseq platform targeted gene-rich regions of the blue wildebeest genome. Finally, functional annotation of the marker sequences revealed a wide range of different putative functions indicating that these SNP markers can be useful in functional gene studies. The DArTseq platform, therefore, represents a high-throughput, robust and cost-effective genotyping platform, which may find adoption in several other African antelope and animal species.

Subspecific identity and a comparison of genetic diversity between wild and ex situ wildebeest.

The original North American ex situ wildebeest population was believed to originate from the white-bearded wildebeest (Connochaetes taurinus albojubatus), which is both morphologically distinct and geographically separated from the brindled wildebeest (C. t. taurinus). However, after an import of wildebeest into North America in 2001, managers have suspected that white-bearded and brindled wildebeest were mixed in herds at multiple institutions. We sequenced the mitochondrial control region (d-loop) from a portion of the managed North American population and compared our sequences with previously published sequences from wild individuals to determine the subspecific identity and genetic diversity of our ex situ population. We were able to confidently identify C. t. albojubatus as the subspecies identity of the sampled portion of our population. Within our population, haplotype and nucleotide diversity were low (0.169 and 0.001, respectively) with a single common haplotype (H1) containing 41 of the 45 individuals sequenced, while two rare haplotypes (H2 and H3) were derived from three individuals and a single individual, respectively. Nucleotide and haplotype diversity were greater overall in the wild populations compared with our managed population. However, C. t. albojubatus was found to exhibit lower nucleotide diversity in both wild and ex situ populations when compared to other wild subspecies. Though the overall goal of the North American wildebeest population is for public education and not reintroduction, maintaining genetic diversity is vital for the long-term viability of this managed population, which may benefit from periodic supplementation of wild animals.

A high-density BAC physical map covering the entire MHC region of addax antelope genome.

BACKGROUND: The mammalian major histocompatibility complex (MHC) harbours clusters of genes associated with the immunological defence of animals against infectious pathogens. At present, no complete MHC physical map is available for any of the wild ruminant species in the world. RESULTS: The high-density physical map is composed of two contigs of 47 overlapping bacterial artificial chromosome (BAC) clones, with an average of 115 Kb for each BAC, covering the entire addax MHC genome. The first contig has 40 overlapping BAC clones covering an approximately 2.9 Mb region of MHC class I, class III, and class IIa, and the second contig has 7 BAC clones covering an approximately 500 Kb genomic region that harbours MHC class IIb. The relative position of each BAC corresponding to the MHC sequence was determined by comparative mapping using PCR screening of the BAC library of 192,000 clones, and the order of BACs was determined by DNA fingerprinting. The overlaps of neighboring BACs were cross-verified by both BAC-end sequencing and co-amplification of identical PCR fragments within the overlapped region, with their identities further confirmed by DNA sequencing. CONCLUSIONS: We report here the successful construction of a high-quality physical map for the addax MHC region using BACs and comparative mapping. The addax MHC physical map we constructed showed one gap of approximately 18 Mb formed by an ancient autosomal inversion that divided the MHC class II into IIa and IIb. The autosomal inversion provides compelling evidence that the MHC organizations in all of the ruminant species are relatively conserved.

Drivers of sex ratio bias in the eastern bongo: lower inbreeding increases the probability of being born male.

Parent sex ratio allocation has consequences for individual fitness, population dynamics, and conservation. Theory predicts that parents should adjust offspring sex ratio when the fitness returns of producing male or female offspring varies. Previous studies have assumed that only mothers are capable of biasing offspring sex ratios, but have neglected fathers, given the expectation of an equal proportion of X- and Y-chromosome-bearing (CBS) sperm in ejaculates due to sex chromosome segregation at meiosis. This assumption has been recently refuted and both paternal fertility and paternal genetic quality have been shown to bias sex ratios. Here, we simultaneously test the relative contribution of paternal, maternal, and individual genetic quality, as measured by inbreeding, on the probability of being born a son or a daughter, using pedigree and lifelong offspring sex ratio data for the eastern bongo ( Tragelaphus eurycerus isaaci). Our models showed first, that surprisingly, as individual inbreeding decreases the probability of being born male increases, second, that paternal genetic effects on sex ratio were stronger than maternal genetic effects (which were absent). Furthermore, paternal effects were opposite in sign to those predicted; father inbreeding increases the probability of having sons. Previous paternal effects have been interpreted as adaptive due to sex-specific inbreeding depression for reproductive traits. We argue that in the eastern bongo, the opposite sign of the paternal effect on sex ratios results from a reversed sex-specific inbreeding depression pattern (present for female but not male reproductive traits). We anticipate that this research will help stimulate research on evolutionary constraints to sex ratios. Finally, the results open a new avenue of research to predict sex ratio allocation in an applied conservation context. Future models of sex ratio allocation should also include the predicted inbreeding level of the offspring and paternal inbreeding levels.

Multilocus nuclear markers provide new insights into the origin and evolution of the blackbuck (Antilope cervicapra, Bovidae).

Evolutionary relationships between members of the Antilopina taxon have been much debated in recent years. The 'true antelope' clade is currently comprised of 4 genera viz., Gazella, Nanger, Eudorcas and the monotypic genus Antilope, that includes A. cervicapra. Most studies have focused on the mitochondrial genome or morphological data to study their relationships. However, signals from mitochondrial data can often be misleading when compared with nuclear markers, as has been shown in multiple taxonomic groups. In this study, we revisit the phylogenetic relationships among members of Antilopina, particularly the phylogenetic position of A. cervicapra, using 12 nuclear markers and compare it with the mitochondrial tree. Furthermore, we explore the implications of the results of this study on the taxonomy and biogeography of Indian antelopes. The nuclear phylogenetic trees built using multiple coalescent and concatenated methods all supported a paraphyletic genus Gazella. Antilope was nested within Gazella as opposed to being sister to it, which was suggested by previous studies and our results based on mitochondrial markers. Our fossil-calibrated larger bovid phylogeny, based on nuclear markers, suggested that the Antilope lineage diverged from its sister species more recently in the Pleistocene, rather than in late Miocene as per previous studies. Our biogeographic analyses suggest that the lineage leading to genus Antilope dispersed into India from the Saharo-Arabian realm around 2 mya, post the expansion of grasslands. We speculate that the adaptations of this savanna-grassland specialist did not allow them to extend their range beyond the Indian subcontinent. Whereas, the only other true antelope in India, G. bennetti, extended its range into India more recently, probably after the establishment of the Thar desert in northwest India.

Investigating the genetic diversity and presence of forensically informative nucleotide sequences in Indian antelope (Antilope cervicapra) using multiple genes of the mitochondrial genome.

Indian antelope or Blackbuck (Antilope cervicapra) is one of the widely distributed endemic species in India among wild bovids and a majority of preferred habitats are in human-dominated landscapes. Poaching threats and habitat degradation are major factors for the decline in Blackbuck population from its distribution range. Till date, there is no detailed study using molecular techniques in India on Blackbuck, except a few studies entailing phylogenetic scenario based on inadequate sampling and DNA sequences restricted over limited geographic areas. In view of this, the present study is aimed to screen the Blackbuck samples from a large part of its distribution range and to investigate the genetic diversity as well as to identify the forensically informative nucleotide sequences (FINS) for species identification. We relied on multi-genes approach using three genes of mtDNA genome viz. Cytochrome Oxidase I, Cytochrome b and 16S rRNA and identified the FINS in the Blackbuck population along with conspecific sequences divergence and genetic diversity indices. In all three genes, we observed 8 to 17 haplotypes with the intra-species sequence divergence of 0.004-0.016. Inter-species sequence divergence with the other closely related species of the Blackbuck was 0.0225-0.033. We report the presence of FINS across three genes from 12 to 18 and found more informative nucleotide sites using Cytochrome Oxidase I genes compared to Cytochrome b and 16S rRNA gene. We did not observe the presence of geographic-specific FINS amongst Blackbuck population that can be used to assign individuals to geographic origin. Besides, in the phylogenetic tree, samples from different locations did not cluster into geographic-specific clade and exhibited mixed homology for these sequences. We suggest exploring the feasibility of using nuclear markers for population assignment.

Multi-locus phylogeny of the tribe Tragelaphini (Mammalia, Bovidae) and species delimitation in bushbuck: Evidence for chromosomal speciation mediated by interspecific hybridization.

The bushbuck is the most widespread bovid species in Africa. Previous mitochondrial studies have revealed a polyphyletic pattern suggesting the possible existence of two distinct species. To assess this issue, we have sequenced 16 nuclear genes and one mitochondrial fragment (cytochrome b gene + control region) for most species of the tribe Tragelaphini, including seven bushbuck individuals belonging to the two divergent mtDNA haplogroups, Scriptus and Sylvaticus. Our phylogenetic analyses show that the Scriptus lineage is a sister-group of Sylvaticus in the nuclear tree, whereas it is related to Tragelaphus angasii in the mitochondrial tree. This mito-nuclear discordance indicates that the mitochondrial genome of Scriptus was acquired by introgression after one or several past events of hybridization between bushbuck and an extinct species closely related to T. angasii. The division into two bushbuck species is supported by the analyses of nuclear markers and by the karyotype here described for T. scriptus (2n = 57 M/58F), which is strikingly distinct from the one previously found for T. sylvaticus (2n = 33 M/34F). Molecular dating estimates suggest that the two species separated during the Early Pleistocene after an event of interspecific hybridization, which may have mediated massive chromosomal rearrangements in the common ancestor of T. scriptus.

Evaluating the role of Pleistocene refugia, rivers and environmental variation in the diversification of central African duikers (genera Cephalophus and Philantomba).

BACKGROUND: This study aims to assess the role that Pleistocene refugia, rivers and local habitat conditions may have played in the evolutionary diversification of three central African duiker species (Cephalophus dorsalis, C. callipygus and Philantomba monticola). Genetic data from geo-referenced feces were collected from a wide range of sites across Central Africa. Historical patterns of population genetic structure were assessed using a ~ 650 bp fragment of the mitochondrial control region and contemporary patterns of genetic differentiation were evaluated using 12 polymorphic microsatellite loci. RESULTS: Mitochondrial analyses revealed that populations of C. callipygus and P. monticola in the Gulf of Guinea refugium are distinct from other populations in west central Africa. All three species exhibit signatures of past population expansion across much of the study area consistent with a history of postglacial expansion. There was no strong evidence for a riverine barrier effect in any of the three species, suggesting that duikers can readily cross major rivers. Generalized dissimilarity models (GDM) showed that environmental variation explains most of the nuclear genetic differentiation in both C. callipygus and P. monticola. The forest-savanna transition across central Cameroon and the Plateaux Batéké region in southeastern Gabon show the highest environmentally-associated turnover in genetic variability. A pattern of genetic differentiation was also evident between the coast and forest interior that may reflect differences in precipitation and/or vegetation. CONCLUSIONS: Findings from this study highlight the historical impact of Pleistocene fragmentation and current influence of environmental variation on genetic structure in duikers. Conservation efforts should therefore target areas that harbor as much environmentally-associated genetic variation as possible in order to maximize species' capacity to adapt to environmental change.

Effect of species-specific differences in chromosome morphology on chromatin compaction and the frequency and distribution of RAD51 and MLH1 foci in two bovid species: cattle (Bos taurus) and the common eland (Taurotragus oryx).

Meiotic recombination between homologous chromosomes is crucial for their correct segregation into gametes and for generating diversity. We compared the frequency and distribution of MLH1 foci and RAD51 foci, synaptonemal complex (SC) length and DNA loop size in two related Bovidae species that share chromosome arm homology but show an extreme difference in their diploid chromosome number: cattle (Bos taurus, 2n = 60) and the common eland (Taurotragus oryx, 2nmale = 31). Compared to cattle, significantly fewer MLH1 foci per cell were observed in the common eland, which can be attributed to the lower number of initial double-strand breaks (DSBs) detected as RAD51 foci in leptonema. Despite the significantly shorter total autosomal SC length and longer DNA loop size of the common eland bi-armed chromosomes compared to those of bovine acrocentrics, the overall crossover density in the common eland was still lower than in cattle, probably due to the reduction in the number of MLH1 foci in the proximal regions of the bi-armed chromosomes. The formation of centric fusions during karyotype evolution of the common eland accompanied by meiotic chromatin compaction has greater implications in the reduction in the number of DSBs in leptonema than in the decrease of MLH1 foci number in pachynema.