Limited Maternal Gene Flow Amongst Elk in the Greater Yellowstone Ecosystem Revealed by Mitochondrial DNA

We quantified maternal patterns of population genetic structure to help understand gene flow among elk (Cervus elaphus) populations across the Greater Yellowstone Ecosystem. We sequenced 596 base pairs of the mitochondrial (mt)DNA control region of 407 elk from nine populations. Our analysis revealed high mtDNA diversity within populations including 12 haplotypes per population on average, and a mean haplotype diversity (i.e., gene diversity) of 0.84. The FST from mtDNA was high (mean FST = 0.162; P = 0.0001) compared to FST for nuclear microsatellites data (FST = 0.006, P = 0.125), which suggested relatively low female movement among populations, perhaps due to female philopatry. Genetic distance (mtDNA pair-wise FST) was not significantly correlated with geographic (Euclidean) distance between populations (Mantel’s r = 0.274, P = 0.168). The lack of isolation-by-distance and large genetic distance between geographically close populations (< 65 km) suggest that maternal gene flow is reduced by certain landscape features (e.g., large, non-forested valleys with roads), which is important for understanding and modeling landscape connectivity and related processes

Hidden Consequences of Living in a Wormy World: Nematode-Induced Immune Suppression Facilitates Tuberculosis Invasion in African Buffalo

Most hosts are infected with multiple parasites, and responses of the immune system to co occurring parasites may influence disease spread. Helminth infection can bias the host immune response toward a T-helper type 2 Th2) over a type 1 Th1) response, impairing the host\u27s ability to control concurrent intracellular microparasite infections and potentially modifying disease dynamics. In humans, immune-mediated interactions between helminths and microparasites can alter host susceptibility to diseases such as HIV, tuberculosis TB), and malaria. However, the extent to which similar processes operate in natural animal populations and influence disease spread remains unknown. We used cross-sectional, experimental, and genetic studies to show that gastrointestinal nematode infection alters immunity to intracellular microparasites in free-ranging African buffalo Syncerus caffer). Buffalo that were more resistant to nematode infection had weaker Th1 responses, there was significant genotypic variation in nematode resistance, and anthelminthic treatment enhanced Th1 immunity. Using a disease dynamic model parameterized with empirical data, we found that nematode-induced immune suppression can facilitate the invasion of bovine TB in buffalo. In the absence of nematodes, TB failed to invade the system, illustrating the critical role nematodes may play in disease establishment. Our results suggest that helminths, by influencing the likelihood of microparasite invasion, may influence patterns of disease emergence in the wild

Evidence for adaptive introgression of exons across a hybrid swarm in deer

Background: Secondary contact between closely related lineages can result in a variety of outcomes, including hybridization, depending upon the strength of reproductive barriers. By examining the extent to which different parts of the genome introgress, it is possible to infer the strength of selection and gain insight into the evolutionary trajectory of lineages. Following secondary contact approximately 8000 years ago in the Pacific Northwest, mule deer (Odocoileus hemionus hemionus) and black-tailed deer (O. h. columbianus) formed a hybrid swarm along the Cascade mountain range despite substantial differences in body size (up to two times) and habitat preference. In this study, we examined genetic population structure, extent of introgression, and selection pressures in freely interbreeding populations of mule deer and black-tailed deer using mitochondrial DNA sequences, 9 microsatellite loci, and 95 SNPs from protein-coding genes. Results: We observed bi-directional hybridization and classified approximately one third of the 172 individuals as hybrids, almost all of which were beyond the F1 generation. High genetic differentiation between black-tailed deer and mule deer at protein-coding genes suggests that there is positive divergent selection, though selection on these loci is relatively weak. Contrary to predictions, there was not greater selection on protein-coding genes thought to be associated with immune function and mate choice. Geographic cline analyses were consistent across genetic markers, suggesting long-term stability (over hundreds of generations), and indicated that the center of the hybrid swarm is 20-30 km to the east of the Cascades ridgeline, where there is a steep ecological transition from wet, forested habitat to dry, scrub habitat. Conclusions: Our data are consistent with a genetic boundary between mule deer and black-tailed deer that is porous but maintained by many loci under weak selection having a substantial cumulative effect. The absence of clear reproductive barriers and the consistent centering of geographic clines at a sharp ecotone suggests that ecology is a driver of hybrid swarm dynamics. Adaptive introgression in this study (and others) promotes gene flow and provides valuable insight into selection strength on specific genes and the evolutionary trajectory of hybridizing taxa

Limitations to estimating bacterial cross-species transmission using genetic and genomic markers: inferences from simulation modeling

Cross-species transmission (CST) of bacterial pathogens has major implications for human health, livestock, and wildlife management because it determines whether control actions in one species may have subsequent effects on other potential host species. The study of bacterial transmission has benefitted from methods measuring two types of genetic variation: variable number of tandem repeats (VNTRs) and single nucleotide polymorphisms (SNPs). However, it is unclear whether these data can distinguish between different epidemiological scenarios. We used a simulation model with two host species and known transmission rates (within and between species) to evaluate the utility of these markers for inferring CST. We found that CST estimates are biased for a wide range of parameters when based on VNTRs and a most parsimonious reconstructed phylogeny. However, estimations of CST rates lower than 5% can be achieved with relatively low bias using as low as 250 SNPs. CST estimates are sensitive to several parameters, including the number of mutations accumulated since introduction, stochasticity, the genetic difference of strains introduced, and the sampling effort. Our results suggest that, even with whole-genome sequences, unbiased estimates of CST will be difficult when sampling is limited, mutation rates are low, or for pathogens that were recently introduced

High Polymorphism in the K-Casein (CSN3) Gene from Wild and Domestic Caprine Species Revealed by DNA Sequencing

We assessed polymorphisms in exon IV of the k-casein gene (CSN3) in ten different breeds of domestic goat (Capra hircus) from three continents and in three related wild caprine taxa (Capra ibex, Capra sibirica and Capra aegagrus). Thirty-five DNA samples were sequenced within a 558 bp fragment of exon IV. Nine polymorphic sites were identified in domestic goat, including four new polymorphisms. In addition to four previously described polymorphic positions, a total of 13 polymorphisms allowed the identification of 13 DNA variants, corresponding to 10 protein variants. Because of conflicting nomenclature of these variants, we propose a standardized allele designation. CSN3*A, CSN3*B, and CSN3*D were found as widely distributed alleles in European goat breeds. Within Capra ibex we identified three variants and showed that the sequence of Capra aegagrus is identical to the most common Capra hircus variant, consistent with Capra aegagrus being the wild progenitor of domestic goats. A dendrogram was drawn to represent the molecular network between the caprine CSN3 variants

Hidden Consequences of Living in a Wormy World: Nematode‐Induced Immune Suppression Facilitates Tuberculosis Invasion in African Buffalo

Most hosts are infected with multiple parasites, and responses of the immune system to co-occurring parasites may influence disease spread. Helminth infection can bias the host immune response toward a T-helper type 2 (Th2) over a type 1 (Th1) response, impairing the host’s ability to control concurrent intracellular microparasite infections and potentially modifying disease dynamics. In humans, immune-mediated interactions between helminths and microparasites can alter host susceptibility to diseases such as HIV, tuberculosis (TB), and malaria. However, the extent to which similar processes operate in natural animal populations and influence disease spread remains unknown.We used cross-sectional, experimental, and genetic studies to show that gastrointestinal nematode infection alters immunity to intracellular microparasites in free-ranging African buffalo (Syncerus caffer). Buffalo that were more resistant to nematode infection had weaker Th1 responses, there was significant genotypic variation in nematode resistance, and anthelminthic treatment enhanced Th1 immunity. Using a disease dynamic model parameterized with empirical data, we found that nematode-induced immune suppression can facilitate the invasion of bovine TB in buffalo. In the absence of nematodes, TB failed to invade the system, illustrating the critical role nematodes may play in disease establishment. Our results suggest that helminths, by influencing the likelihood of microparasite invasion, may influence patterns of disease emergence in the wild

Estimating Bighorn Sheep (\u3ci\u3eOvis Canadensis\u3c/i\u3e) Abundance Using Noninvasive Sampling at a Mineral Lick within a National Park Wilderness Area

Conservation of species requires accurate population estimates. We used genetic markers from feces to determine bighorn sheep abundance for a herd that was hypothesized to be declining and in need of population status monitoring. We sampled from a small but accessible portion of the population\u27s range where animals natural congregate at a natural mineral lick to test whether we could accurately estimate population size by sampling from an area where animals concentrate. We used mark-recapture analysis to derive population estimates, and compared estimates from this smaller spatial sampling to estimates from sampling of the entire bighorn sheep range. We found that estimates were somewhat comparable; in 2009, the mineral lick sample and entire range sample differed by 20 individuals, and in 2010 they differed by only one individual. However, we captured 13 individuals in the entire range sample that were not captured at the mineral lick, and thus violated a model assumption that all individuals had an equal opportunity of being captured. This eliminated the possibility of inferring a total population estimate from just animals visiting the mineral lick, but because estimates were relatively similar, monitoring at the mineral lick can provide a useful index for management and conservation. We compared our results to a radio-collar study conducted in 2003-2004 and confirmed that the population remained stable since 2004. Our population estimates were 78 (CI 62-114) in 2009 and 95 (CI 77-131) in 2010. Between 7 and 11 sampling dates were needed to achieve a CV of 20% for population estimates, assuming a capture probability between 0.09 and 0.13. We relied on citizen science volunteers to maximize data collection and reduce costs; 71% of all fecal samples were collected by volunteers, compared to 29% collected by paid staff. We conclude that our technique provides a useful monitoring tool for managers. The technique could be tested and applied in similar populations where animals congregate with high fidelity at a mineral lick or other area

Lessons for conservation management: Monitoring temporal changes in genetic diversity of Cape mountain zebra (Equus zebra zebra)

The Cape mountain zebra (Equus zebra zebra) is a subspecies of mountain zebra endemic to South Africa. The Cape mountain zebra experienced near extinction in the early 1900's and their numbers have since recovered to more than 4,800 individuals. However, there are still threats to their long-term persistence. A previous study reported that Cape mountain zebra had low genetic diversity in three relict populations and that urgent conservation management actions were needed to mitigate the risk of further loss. As these suggestions went largely unheeded, we undertook the present study, fifteen years later to determine the impact of management on genetic diversity in three key populations. Our results show a substantial loss of heterozygosity across the Cape mountain zebra populations studied. The most severe losses occurred at De Hoop Nature Reserve where expected heterozygosity reduced by 22.85% from 0.385 to 0.297. This is alarming, as the De Hoop Nature Reserve was previously identified as the most genetically diverse population owing to its founders originating from two of the three remaining relict stocks. Furthermore, we observed a complete loss of multiple private alleles from all populations, and a related reduction in genetic structure across the subspecies. These losses could lead to inbreeding depression and reduce the evolutionary potential of the Cape mountain zebra. We recommend immediate implementation of evidence-based genetic management and monitoring to prevent further losses, which could jeopardise the long term survival of Cape mountain zebra, especially in the face of habitat and climate change and emerging diseases

Sex-Biased Gene Flow Among Elk in the Greater Yellowstone Ecosystem

We quantified patterns of population genetic structure to help understand gene flow among elk populations across the Greater Yellowstone Ecosystem. We sequenced 596 base pairs of the mitochondrial control region of 380 elk from eight populations. Analysis revealed high mitochondrial DNA variation within populations, averaging 13.0 haplotypes with high mean gene diversity (0.85). The genetic differentiation among populations for mitochondrial DNA was relatively high (FST = 0.161; P = 0.001) compared to genetic differentiation for nuclear microsatellite data (FST = 0.002; P = 0.332), which suggested relatively low female gene flow among populations. The estimated ratio of male to female gene flow (mm/mf = 46) was among the highest we have seen reported for large mammals. Genetic distance (for mitochondrial DNA pairwise FST) was not significantly correlated with geographic (Euclidean) distance between populations (Mantel’s r = 0.274, P = 0.168). Large mitochondrial DNA genetic distances (e.g., FST . 0.2) between some of the geographically closest populations (,65 km) suggested behavioral factors and/or landscape features might shape female gene flow patterns. Given the strong sex-biased gene flow, future research and conservation efforts should consider the sexes separately when modeling corridors of gene flow or predicting spread of maternally transmitted diseases. The growing availability of genetic data to compare male vs. female gene flow provides many exciting opportunities to explore the magnitude, causes, and implications of sex-biased gene flow likely to occur in many species