Full length MHC IIβ exon 2 primers for salmonids: a new resource for next generation sequencing

Genomic and environmental influences on resilience in a cold-water fish near the edge of its range

Small, isolated populations present a challenge for conservation. The dueling effects of selection and drift in a limited pool of genetic diversity make the responses of small populations to environmental perturbations erratic and difficult to predict. This is particularly true at the edge of a species range, where populations often persist at the limits of their environmental tolerances. Populations of cisco, Coregonus artedi, in inland lakes have experienced numerous extirpations along the southern edge of their range in recent decades, which are thought to result from environmental degradation and loss of cold, well-oxygenated habitat as lakes warm. Yet, cisco extirpations do not show a clear latitudinal pattern, suggesting that local environmental factors and potentially local adaptation may influence resilience. Here, we used genomic tools to investigate the nature of this pattern of resilience. We used restriction site-associated DNA capture (Rapture) sequencing to survey genomic diversity and differentiation in southern inland lake cisco populations and compared the frequency of deleterious mutations that potentially influence fitness across lakes. We also examined haplotype diversity in a region of the major histocompatibility complex involved in stress and immune system response. We correlated these metrics to spatial and environmental factors including latitude, lake size, and measures of oxythermal habitat and found significant relationships between genetic metrics and broad and local factors. High levels of genetic differentiation among populations were punctuated by a phylogeographic break and residual patterns of isolation-by-distance. Although the prevalence of deleterious mutations and inbreeding coefficients was significantly correlated with latitude, neutral and non-neutral genetic diversity were most strongly correlated with lake surface area. Notably, differences among lakes in the availability of estimated oxythermal habitat left no clear population genomic signature. Our results shed light on the complex dynamics influencing these isolated populations and provide valuable information for their conservation.

Neutral variation does not predict immunogenetic variation in the European grayling (Thymallus thymallus)-implications for management

Preservation of genetic diversity is critical to successful conservation, and there is increasing demand for the inclusion of ecologically meaningful genetic information in management decisions. Supportive breeding programmes are increasingly implemented to combat declines in many species, yet their effect on adaptive genetic variation is understudied. This is despite the fact that supportive breeding may interfere with natural evolutionary processes. Here, we assessed the performance of neutral and adaptive markers (major histocompatibility complex; MHC) to inform management of European grayling (Thymallus thymallus), which routinely involves supplementation of natural populations with hatchery-reared fish (stocking). This study is the first to characterize MH II DAA and DAB loci in grayling and to investigate immune genetic variation in relation to management practice in this species. High-throughput Illumina sequencing of "introduced," "stocked native" and "non-stocked native" populations revealed significantly higher levels of allelic richness and heterozygosity for MH markers than microsatellites exclusively in non-stocked native populations. Likewise, significantly lower differentiation at the MH II than for microsatellites was apparent when considering non-stocked native populations, but not stocked populations. We developed a simulation model to test the effects of relaxation of selection during the early life stage within captivity. Dependent on the census population size and stocking intensity, there may be long-term effects of stocking on MH II, but not neutral genetic diversity. This is consistent with our empirical results. This study highlights the necessity for considering adaptive genetic variation in conservation decisions and raises concerns about the efficiency of stocking as a management practice.

Small-scale intraspecific patterns of adaptive immunogenetic polymorphisms and neutral variation in Lake Superior lake trout

Many fishes express high levels of intraspecific variability, often linked to resource partitioning. Several studies show that a species' evolutionary trajectory of adaptive divergence can undergo reversals caused by changes in its environment. Such a reversal in neutral genetic and morphological variation among lake trout Salvelinus namaycush ecomorphs appears to be underway in Lake Superior. However, a water depth gradient in neutral genetic divergence was found to be associated with intraspecific diversity in the lake. To investigate patterns of adaptive immunogenetic variation among lake trout ecomorphs, we used Illumina high-throughput sequencing. The population's genetic structure of the major histocompatibility complex (MHC Class IIβ exon 2) and 18 microsatellite loci were compared to disentangle neutral and selective processes at a small geographic scale. Both MHC and microsatellite variation were partitioned more by water depth stratum than by ecomorph. Several metrics showed strong clustering by water depth in MHC alleles, but not microsatellites. We report a 75% increase in the number of MHC alleles shared between the predominant shallow and deep water ecomorphs since a previous lake trout MHC study at the same locale (c. 1990s data). This result is consistent with the reverse speciation hypothesis, although adaptive MHC polymorphisms persist along an ecological gradient. Finally, results suggested that the lake trout have multiple copies of the MHC II locus consistent with a historic genomic duplication event. Our findings indicated that conservation approaches for this species could focus on managing various ecological habitats by depth, in addition to regulating the fisheries specific to ecomorphs.

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake, the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long-term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.

Neutral and selective processes shape MHC gene diversity and expression in stocked brook charr populations (Salvelinus fontinalis)

The capacity of an individual to battle infection is an important fitness determinant in wild vertebrate populations. The major histocompatibility complex (MHC) genes are crucial for a host's adaptive immune system to detect pathogens. However, anthropogenic activities may disrupt natural cycles of co-evolution between hosts and pathogens. In this study, we investigated the dynamic sequence and expression variation of host parasite interactions in brook charr (Salvelinus fontinalis) in a context of past human disturbance via population supplementation from domestic individuals. To do so, we developed a new method to examine selection shaping MHC diversity within and between populations and found a complex interplay between neutral and selective processes that varied between lakes that were investigated. We provided evidence for a lower introgression rate of domestic alleles and found that parasite infection increased with domestic genomic background of individuals. We also documented an association between individual MHC alleles and parasite taxa. Finally, longer cis-regulatory minisatellites were positively correlated with MHC II down-regulation and domestic admixture, suggesting that inadvertent selection during domestication resulted in a lower immune response capacity, through a trade-off between growth and immunity, which explained the negative selection of domestic alleles at least under certain circumstances.