The functional gene diversity in natural populations over postglacial areas: the shaping mechanisms behind genetic composition of longnose dace (Rhinichthys cataractae) in northeastern North America

Expression and Variations in EPAS1 Associated with Oxygen Metabolism in Sheep

Endothelial PAS domain protein 1 gene (EPAS1) is a member of the HIF gene family. This gene encodes a transcription factor subunit that is involved in the induction of oxygen-regulated genes. Several studies have demonstrated that a mutation in EPAS1 could affect oxygen sensing, polycythemia, and hemoglobin level. However, whether EPAS1 mutation affects sheep oxygen metabolism is still unknown. Therefore, we explored the relationship between the variation of EPAS1 and oxygen metabolism in sheep. In this study, variations in ovine EPAS1 exon 15 were investigated in 332 Tibetan sheep and 339 Hu sheep by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) analysis. In addition, we studied the effect of these variations on blood gas in 176 Tibetan sheep and 231 Hu sheep. Finally, the mRNA expression of EPAS1 in six tissues of Hu sheep and Tibetan sheep living at different altitudes (2500 m, 3500 m, and 4500 m) was analyzed by real-time quantitative PCR (RT-qPCR). Four alleles (A, B, C, and D) were detected, and their distributions highly differed between Tibetan sheep and Hu sheep. In Tibetan sheep, B was the dominant allele, and C and D alleles were rare, whereas all four alleles were common in Hu sheep. Six single nucleotide polymorphisms (SNPs) were identified between the four alleles and one of them was non-synonymous (p.F606L). While studying the blood gas levels in Tibetan sheep and Hu sheep, one variant region was found to be associated with an elevated pO₂ and sO₂, which suggested that variations in EPAS1 are associated with oxygen metabolism in sheep. RT-qPCR results showed that EPAS1 was expressed in the six tissues of Hu sheep and Tibetan sheep at different altitudes. In addition, the expression of EPAS1 in four tissues (heart, liver, spleen, and longissimus dorsi muscle) of Hu sheep was lower than that in Tibetan sheep from three different altitudes, and the expression of EPAS1 was positively correlated with the altitude. These results indicate that the variations and expression of EPAS1 is closely related to oxygen metabolism.

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.

Spatio-temporal variation in parasite communities maintains diversity at the major histocompatibility complex class IIβ in the endangered Rio Grande silvery minnow

Climate change will strongly impact aquatic ecosystems particularly in arid and semi-arid regions. Fish-parasite interactions will also be affected by predicted altered flow and temperature regimes, and other environmental stressors. Hence, identifying environmental and genetic factors associated with maintaining diversity at immune genes is critical for understanding species' adaptive capacity. Here, we combine genetic (MHC class IIβ and microsatellites), parasitological and ecological data to explore the relationship between these factors in the remnant wild Rio Grande silvery minnow (Hybognathus amarus) population, an endangered species found in the southwestern United States. Infections with multiple parasites on the gills were observed and there was spatio-temporal variation in parasite communities and patterns of infection among individuals. Despite its highly endangered status and chronically low genetic effective size, Rio Grande silvery minnow had high allelic diversity at MHC class IIβ with more alleles recognized at the presumptive DAB1 locus compared to the DAB3 locus. We identified significant associations between specific parasites and MHC alleles against a backdrop of generalist parasite prevalence. We also found that individuals with higher individual neutral heterozygosity and higher amino acid divergence between MHC alleles had lower parasite abundance and diversity. Taken together, these results suggest a role for fluctuating selection imposed by spatio-temporal variation in pathogen communities and divergent allele advantage in maintenance of high MHC polymorphism. Understanding the complex interaction of habitat, pathogens and immunity in protected species will require integrated experimental, genetic and field studies.

Drift, not selection, shapes toll-like receptor variation among oceanic island populations

Understanding the relative role of different evolutionary forces in shaping the level and distribution of functional genetic diversity among natural populations is a key issue in evolutionary and conservation biology. To do so accurately genetic data must be analysed in conjunction with an unambiguous understanding of the historical processes that have acted upon the populations. Here, we focused on diversity at toll‐like receptor (TLR) loci, which play a key role in the vertebrate innate immune system and, therefore, are expected to be under pathogen‐mediated selection. We assessed TLR variation within and among 13 island populations (grouped into three archipelagos) of Berthelot's pipit, Anthus berthelotii, for which detailed population history has previously been ascertained. We also compared the variation observed with that found in its widespread sister species, the tawny pipit, Anthus campestris. We found strong evidence for positive selection at specific codons in TLR1LA, TLR3 and TLR4. Despite this, we found that at the allele frequency level, demographic history has played the major role in shaping patterns of TLR variation in Berthelot's pipit. Levels of diversity and differentiation within and across archipelagos at all TLR loci corresponded very closely with neutral microsatellite variation and with the severity of the bottlenecks that occurred during colonization. Our study shows that despite the importance of TLRs in combating pathogens, demography can be the main driver of immune gene variation within and across populations, resulting in patterns of functional variation that can persist over evolutionary timescales.

Phylogeography of the longnose dace (Rhinichthys cataractae) species group in northwestern North America — the Nooksack dace problem

The longnose dace (Rhinichthys cataractae (Valenciennes, 1842); Cyprinidae) is one of the most widespread freshwater fishes in North America, and across its range there have been several divergent forms described that are of uncertain taxonomic status. One of these forms, the Nooksack dace, is found in southwestern British Columbia and adjacent portions of western Washington, and is distinguished from longnose dace by a lower number of lateral-line scales. We sequenced a total of approximately 1400 base pairs (bp) of mitochondrial DNA (mtDNA) and noted that the longnose dace found west of the Continental Divide and Nooksack dace constituted reciprocally monophyletic clades that differed from each other by between 2% and 3% sequence divergence. Sequence analysis at two nuclear loci (the S7 ribosomal protein intron 1 (S7) and recombination activation gene 1 (RAG1)), however, showed no consistent difference between longnose dace and Nooksack dace and several alleles were shared between them. By contrast, consistent differences at both mtDNA and nuclear DNA loci were resolved between R. cataractae samples from east and west of the Continental Divide. The Nooksack dace does not appear to warrant separate taxonomic status from the longnose dace, but the mtDNA differences support its recognition as an important component of the evolutionary and biogeographic legacy of R. cataractae.