Impact of historical founder effects and a recent bottleneck on MHC variability in Commander Arctic foxes (Vulpes lagopus)

Genetic and Anatomical Determinants of Olfaction in Dogs and Wild Canids

Understanding the anatomical and genetic basis of complex phenotypic traits has long been a challenge for biological research. Domestic dogs offer a compelling model as they demonstrate more phenotypic variation than any other vertebrate species. Dogs have been intensely selected for specific traits and abilities, directly or indirectly, over the past 15,000 years since their initial domestication from the gray wolf. Because olfaction plays a central role in critical tasks, such as the detection of drugs, diseases, and explosives, as well as human rescue, we compared relative olfactory capacity across dog breeds and assessed changes to the canine olfactory system to their direct ancestors, wolves, and coyotes. We conducted a cross-disciplinary survey of olfactory anatomy, olfactory receptor (OR) gene variation, and OR gene expression in domestic dogs. Through comparisons to their closest wild canid relatives, the gray wolf and coyote, we show that domestic dogs might have lost functional OR genes commensurate with a documented reduction in nasal morphology as an outcome of the domestication process prior to breed formation. Critically, within domestic dogs alone, we found no genetic or morphological profile shared among functional or genealogical breed groupings, such as scent hounds, that might indicate evidence of any human-directed selection for enhanced olfaction. Instead, our results suggest that superior scent detection dogs likely owe their success to advantageous behavioral traits and training rather than an "olfactory edge" provided by morphology or genes.

The Amphibian Major Histocompatibility Complex-A Review and Future Outlook

The major histocompatibility complex (MHC) is a cluster of functionally related genes encoding proteins which, among other functions, mediate immune system activation. While the MHC of many vertebrates has been extensively studied, less is known about the amphibian MHC. This represents an important knowledge gap because amphibians mark the evolutionary transition from an aquatic to a terrestrial lifestyle and often maintain a biphasic lifestyle. Hence, they tend to be exposed to both aquatic and terrestrial pathogen communities, providing opportunities to gain fundamental insights into how the immune system responds to different environmental challenges. Moreover, amphibians are globally threatened by invasive pathogens and the MHC may play a role in combating population decline. In this review, we summarize the current state of knowledge regarding the amphibian MHC and identify the major differences with other vertebrates. We also review how the number of MHC gene copies varies across amphibian groups and how MHC-based variation relates to amphibian ontogeny, behaviour, disease, and phylogeography. We conclude by identifying knowledge gaps and proposing priorities for future research.

Diversity of the MHC class II DRB gene in the wolverine (Carnivora: Mustelidae: Gulo gulo) in Finland

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck.

Positive selection on the MHC class II DLA-DQA1 gene in golden jackals (Canis aureus) from their recent expansion range in Europe and its effect on their body mass index

Background

In Europe, golden jackals (Canis aureus) have been expanding their range out of the southern and southeastern Balkans towards central Europe continually since the 1960s. Here, we investigated the level of functional diversity at the MHC class II DLA-DQA1 exon 2 in golden jackal populations from Bulgaria, Serbia, and Hungary. Specifically, we tested for positive selection on and geographic variation at that locus due to adaptation to supposedly regionally varying pathogenic landscapes. To test for potential fitness effects of different protein variants on individual body condition, we used linear modeling of individual body mass indexes (bmi) and accounted for possible age, sex, geographical, and climatic effects. The latter approach was performed, however, only on Serbian individuals with appropriate data.

Results

Only three different DLA-DQA1 alleles were detected, all coding for different amino-acid sequences. The neutrality tests revealed no significant but positive values; there was no signal of spatial structuring and no deviation from the Hardy–Weinberg equilibrium across the studied range of expansion. However, we found a signal of trans-species polymorphism and significant test results for positive selection on three codons. Our information-theory based linear modeling results indicated an effect of ambient temperature on the occurrence of individual DLA-DQA1 genotypes in individuals from across the studied expansion range, independent from geographical position. Our linear modeling results of individual bmi values indicated that yearlings homozygous for DLA-DQA1*03001 reached values typical for adults contrary to yearlings carrying other genotypes (protein combinations). This suggested better growth rates and thus a possible fitness advantage of yearlings homozygous for DLA-DQA1*03001.

Conclusions

Our results indicate a demographic (stochastic) signal of reduced DLA-DQA1 exon 2 variation, in line with the documented historical demographic bottleneck. At the same time, however, allelic variation was also affected by positive selection and adaptation to varying ambient temperature, supposedly reflecting geographic variation in the pathogenic landscape. Moreover, an allele effect on body mass index values of yearlings suggested differential fitness associated with growth rates. Overall, a combination of a stochastic effect and positive selection has shaped and is still shaping the variation at the studied MHC locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01856-z.

Genetic origins and diversity of bushpigs from Madagascar (Potamochoerus larvatus, family Suidae)

The island of Madagascar, situated off the southeast coast of Africa, shows the first evidence of human presence ~ 10,000 years ago; however, other archaeological data indicates a settlement of the modern peoples of the island distinctly more recent, perhaps > 1500 years ago. Bushpigs of the genus Potamochoerus (family Suidae), are today widely distributed in Madagascar and presumed to have been introduced from Africa at some stage by human immigrants to the island. However, disparities about their origins in Madagascar have been presented in the literature, including the possibility of endemic subspecies, and few empirical data are available. Furthermore, the separation of bushpigs in Madagascar from their mainland relatives may have favoured the evolution of a different repertoire of immune genes first due to a founder effect and then as a response to distinct pathogens compared to their ancestors. Molecular analysis confirmed the species status of the bushpig in Madagascar as P. larvatus, likely introduced from the central region of southern Africa, with no genetic evidence for the recognition of eastern and western subspecies as suggested from previous cranial morphology examination. Investigation of the immunologically important SLA-DQB1 peptide-binding region showed a different immune repertoire of bushpigs in Madagascar compared to those on the African mainland, with seventeen exon-2 haplotypes unique to bushpigs in Madagascar (2/28 haplotypes shared). This suggests that the MHC diversity of the Madagascar populations may have enabled Malagasy bushpigs to adapt to new environments.

Severe population bottleneck and cranial morphology change in the Mednyi Island subspecies of Arctic fox Vulpes lagopus (Carnivora: Canidae)

Arctic foxes, Vulpes lagopus living on Mednyi Island suffered a drastic decline in population size in the late 1970s due to an outbreak of mange epizootic. This dramatic fall in numbers rendered the subspecies endangered, and the concomitant loss of variability resulted in a population bottleneck. Here, we investigate whether differences in cranial morphology between Mednyi Island Arctic foxes and Bering Island Arctic foxes could be attributed to the severe population bottleneck suffered by the Mednyi population in the 1970s. We used morphometric traits as proxies for genetic data to provide estimates of FST. Results show higher FST estimates for the Mednyi population than for the Bering population, which we interpret as a bottleneck signature. FST results also indicate a pattern of divergence between the two populations consistent with random genetic drift. Bottleneck detection is critical for the interpretation of the demographic history of the endangered Mednyi Island Arctic fox, with consequences for conservation management.

Molecular evolution of MHC class II DRB exon 2 in Japanese and Russian raccoon dogs, Nyctereutes procyonoides (Carnivora: Canidae)

Raccoon dogs, Nyctereutes procyonoides, are native to East Asia, but have been introduced into western Russia and eastern Europe. To determine allelic diversity and elucidate the evolution of major histocompatibility complex (MHC) class II genes in the raccoon dog, we analysed a 237-bp region of DRB exon 2 from 36 individuals of native and introduced populations from Japan and Russia. We detected 23 DRB alleles (Nypr-DRBs), 22 of which were novel. Some alleles were found across the species’ range, while others were geographically restricted. For both native and introduced populations, the ratio of non-synonymous to synonymous substitution rates for codons at predicted antigen-binding sites was significantly greater than 2, indicating that Nypr-DRBs have evolved under positive selection. Mixed effect model evolution analysis and an algorithm to detect recombination showed five positively selected codons and one recombination breakpoint, respectively. Overall, our results suggest that the diversity of MHC class II DRB in N. procyonoides was influenced and maintained by recombination, pathogen-driven positive selection, geographical barriers and the founder effect. A Bayesian phylogenetic tree revealed no evidence of trans-species polymorphism (TSP), but instead showed monophyly for the Nypr-DRB alleles within a larger clade of canid sequences. The lack of TSP may have been due to long-term divergence of raccoon dogs from other canids, or to their having encountered different sets of pathogens due to occupying a different ecological niche.

Spatiotemporal spread of sarcoptic mange in the red fox (Vulpes vulpes) in Switzerland over more than 60 years: lessons learnt from comparative analysis of multiple surveillance tools

BACKGROUND

Sarcoptic mange is a contagious skin disease of wild and domestic mammals caused by the mite Sarcoptes scabiei. Reports of sarcoptic mange in wildlife increased worldwide in the second half of the 20th century, especially since the 1990s. The aim of this study was to provide new insights into the epidemiology of mange by (i) documenting the emergence of sarcoptic mange in the red fox (Vulpes vulpes) in the last decades in Switzerland; and (ii) describing its spatiotemporal spread combining data obtained through different surveillance methods.

METHODS

Retrospective analysis of archived material together with prospective data collection delivered a large dataset from the 19th century to 2018. Methods included: (i) a review of historical literature; (ii) screening of necropsy reports from general health surveillance (1958-2018); (iii) screening of data on mange (1968-1992) collected during the sylvatic rabies eradication campaign; (iv) a questionnaire survey (<1980-2017) and (v) evaluation of camera-trap bycatch data (2005-2018).

RESULTS

Sarcoptic mange in red foxes was reported as early as 1835 in Switzerland. The first case diagnosed in the framework of the general health surveillance was in 1959. Prior to 1980, sarcoptic mange occurred in non-adjacent surveillance districts scattered all over the country. During the period of the rabies epidemic (1970s-early 1990s), the percentage of foxes tested for rabies with sarcoptic mange significantly decreased in subregions with rabies, whereas it remained high in the few rabies-free subregions. Sarcoptic mange re-emerged in the mid-1990s and continuously spread during the 2000-2010s, to finally extend to the whole country in 2017. The yearly prevalence of mange in foxes estimated by camera-trapping ranged from 0.1-12%.

CONCLUSIONS

Sarcoptic mange has likely been endemic in Switzerland as well as in other European countries at least since the mid-19th century. The rabies epidemics seem to have influenced the pattern of spread of mange in several locations, revealing an interesting example of disease interaction in free-ranging wildlife populations. The combination of multiple surveillance tools to study the long-term dynamics of sarcoptic mange in red foxes in Switzerland proved to be a successful strategy, which underlined the usefulness of questionnaire surveys.

Phenotypic integration and modularity drives skull shape divergence in the Arctic fox (Vulpes lagopus) from the Commander Islands

Phenotypic integration and modularity influence morphological disparity and evolvability. However, studies addressing how morphological integration and modularity change for long periods of genetic isolation are scarce. Here, we investigate patterns of phenotypic integration and modularity in the skull of phenotypically and genetically distinct populations of the Artic fox (Vulpes lagopus) from the Commander Islands of the Aleutian belt (i.e. Bering and Mednyi) that were isolated ca 10 000 years by ice-free waters of the Bering sea. We use three-dimensional geometric morphometrics to quantify the strength of modularity and integration from inter-individual variation (static) and from fluctuating asymmetry (random developmental variation) in both island populations compared to the mainland population (i.e. Chukotka) and we investigated how changes in morphological integration and modularity affect disparity and the directionality of trait divergence. Our results indicate a decrease in morphological integration concomitant to an increase in disparity at a developmental level, from mainland to the smallest and farthest population of Mednyi. However, phenotypic integration is higher in both island populations accompanied by a reduction in disparity compared to the population of mainland at a static level. This higher integration may have favoured morphological adaptive changes towards specific feeding behaviours related to the extreme environmental settings of islands. Our study demonstrates how shifts in phenotypic integration and modularity can facilitate phenotypic evolvability at the intraspecific level that may lead to lineage divergence at macroevolutioanry scales.

Genetic Diversity and Differentiation at Structurally Varying MHC Haplotypes and Microsatellites in Bottlenecked Populations of Endangered Crested Ibis

Investigating adaptive potential and understanding the relative roles of selection and genetic drift in populations of endangered species are essential in conservation. Major histocompatibility complex (MHC) genes characterized by spectacular polymorphism and fitness association have become valuable adaptive markers. Herein we investigate the variation of all MHC class I and II genes across seven populations of an endangered bird, the crested ibis, of which all current individuals are offspring of only two pairs. We inferred seven multilocus haplotypes from linked alleles in the Core Region and revealed structural variation of the class II region that probably evolved through unequal crossing over. Based on the low polymorphism, structural variation, strong linkage, and extensive shared alleles, we applied the MHC haplotypes in population analysis. The genetic variation and population structure at MHC haplotypes are generally concordant with those expected from microsatellites, underlining the predominant role of genetic drift in shaping MHC variation in the bottlenecked populations. Nonetheless, some populations showed elevated differentiation at MHC, probably due to limited gene flow. The seven populations were significantly differentiated into three groups and some groups exhibited genetic monomorphism, which can be attributed to founder effects. We therefore propose various strategies for future conservation and management.

Patterns of MHC-DRB1 polymorphism in a post-glacial island canid, the Newfoundland red fox (Vulpes vulpes deletrix), suggest balancing selection at species and population timescales

As the only native insular Newfoundland canid between the extinction of the wolf in the 1930s and the recent arrival of coyotes, the red fox (Vulpes vulpes deletrix Bangs 1898) poses interesting questions about genetic distinctiveness and the post-glacial colonization history of the island’s depauperate mammalian fauna. Here, we characterized genetic variability at the major histocompatibility complex (MHC) class II DR β1 domain (DRB1) locus in 28 red foxes from six sampling localities island-wide and compared it with mitochondrial control region (CR) diversity and DRB1 diversity in other canids. Our goals were to describe novel DRB1 alleles in a new canid population and to make inferences about the role of selection in maintaining their diversity. As in numerous studies of vertebrates, we found an order-of-magnitude higher nucleotide diversity at the DRB1 locus compared with the CR and significantly positive nonsynonymous-to-synonymous substitution ratios, indicative of selection in the distant past. Although the evidence is weaker, the Ewens-Watterson test of neutrality and the geographical distribution of variation compared with the CR suggest a role for selection over the evolutionary timescale of populations. We report the first genetic data from the DRB1 locus in the red fox and establish baseline information regarding immunogenetic variation in this island canid population which should inform continued investigations of population demography, adaptive genetic diversity, and wildlife disease in red foxes and related species.

Infectious disease and red wolf conservation: assessment of disease occurrence and associated risks

Infectious diseases pose a significant threat to global biodiversity and may contribute to extinction. As such, establishing baseline disease prevalence in vulnerable species where disease could affect persistence is important to conservation. We assessed potential disease threats to endangered red wolves (Canis rufus) by evaluating regional (southeastern United States) disease occurrences in mammals and parasite prevalence in red wolves and sympatric coyotes (Canis latrans) in North Carolina. Common viral pathogens in the southeast region, such as canine distemper and canine parvovirus, and numerous widespread endoparasites could pose a threat to the red wolf population. The most prevalent parasites in red wolves and sympatric coyotes were heartworm (Dirofilaria immitis), hookworm (Ancylostoma caninum), and Ehrlichia spp.; several red wolves and coyotes were also positive for bacteria causing Lyme disease (Borrelia burgdorferi). Coyotes had a more species-rich parasite community than red wolves, suggesting they could harbor more parasites and act as a disease reservoir. Species identity and sex did not significantly affect parasite loads, but young canids were less likely to have heartworm and more likely to have high levels of endoparasites. Continued disease monitoring is important for red wolf recovery because low levels of genetic variability may compromise the wolves' abilities to combat novel pathogens from closely related species, such as domestic dogs and coyotes.

Drift rather than selection dominates MHC class II allelic diversity patterns at the biogeographical range scale in natterjack toads Bufo calamita

Study of major histocompatibility complex (MHC) loci has gained great popularity in recent years, partly due to their function in protecting vertebrates from infections. This is of particular interest in amphibians on account of major threats many species face from emergent diseases such as chytridiomycosis. In this study we compare levels of diversity in an expressed MHC class II locus with neutral genetic diversity at microsatellite loci in natterjack toad (Bufo (Epidalea) calamita) populations across the whole of the species' biogeographical range. Variation at both classes of loci was high in the glacial refugium areas (REF) and much lower in postglacial expansion areas (PGE), especially in range edge populations. Although there was clear evidence that the MHC locus was influenced by positive selection in the past, congruence with the neutral markers suggested that historical demographic events were the main force shaping MHC variation in the PGE area. Both neutral and adaptive genetic variation declined with distance from glacial refugia. Nevertheless, there were also some indications from differential isolation by distance and allele abundance patterns that weak effects of selection have been superimposed on the main drift effect in the PGE zone.

Towards the non-invasive assessment of MHC genotype in wild primates: analysis of wild Assamese macaque MHC-DRB from fecal samples

The major histocompatibility complex (MHC) plays an important role in the immune response and may thus crucially affect an individual's fitness, relevant also for studies on evolutionary ecology and wildlife conservation. Detailed knowledge on the genomic organization, polymorphism and diversity of the MHC has a narrow taxonomic focus though and among macaques is only available for rhesus and long-tailed macaques-the species most commonly kept for biomedical research. The lack of data on wild populations is largely due to the difficulty of obtaining blood or tissue samples necessary for genotyping approaches. Here, we aimed at analyzing MHC-DRB from non-invasively collected fecal samples in wild Assamese macaques (Macaca assamensis), utilizing the MHC-DRB-STR (D6S2878) microsatellite marker. Due to the fecal DNA source incomplete genotypes occurred, which may be improved in the future by method refinement. We detected 28 distinct DRB-STR lengths in 43 individuals with individual genotypes containing 1-9 MHC-DRB-STRs and defined four haplotypes segregating between families in Mendelian fashion. Our results indicate that variability and diversity of MHC-DRB in Assamese macaques is comparable to that of other macaque species and importantly, that fecal samples can be used for non-invasive analysis of MHC genes after refinement of the applied methods, opening a number of opportunities for MHC research on natural populations.

Comparison of 454 pyrosequencing methods for characterizing the major histocompatibility complex of nonmodel species and the advantages of ultra deep coverage

Characterization and population genetic analysis of multilocus genes, such as those found in the major histocompatibility complex (MHC) is challenging in nonmodel vertebrates. The traditional method of extensive cloning and Sanger sequencing is costly and time-intensive and indirect methods of assessment often underestimate total variation. Here, we explored the suitability of 454 pyrosequencing for characterizing multilocus genes for use in population genetic studies. We compared two sample tagging protocols and two bioinformatic procedures for 454 sequencing through characterization of a 185-bp fragment of MHC DRB exon 2 in wolverines (Gulo gulo) and further compared the results with those from cloning and Sanger sequencing. We found 10 putative DRB alleles in the 88 individuals screened with between two and four alleles per individual, suggesting amplification of a duplicated DRB gene. In addition to the putative alleles, all individuals possessed an easily identifiable pseudogene. In our system, sequence variants with a frequency below 6% in an individual sample were usually artefacts. However, we found that sample preparation and data processing procedures can greatly affect variant frequencies in addition to the complexity of the multilocus system. Therefore, we recommend determining a per-amplicon-variant frequency threshold for each unique system. The extremely deep coverage obtained in our study (approximately 5000×) coupled with the semi-quantitative nature of pyrosequencing enabled us to assign all putative alleles to the two DRB loci, which is generally not possible using traditional methods. Our method of obtaining locus-specific MHC genotypes will enhance population genetic analyses and studies on disease susceptibility in nonmodel wildlife species.