The major histocompatibility complex of tassel-eared squirrels. II. Genetic diversity associated with Abert squirrels

Genome-wide markers reveal a complex evolutionary history involving divergence and introgression in the Abert's squirrel (Sciurus aberti) species group

Background

Genetic introgression between divergent lineages is now considered more common than previously appreciated, with potentially important consequences for adaptation and speciation. Introgression is often asymmetric between populations and patterns can vary for different types of loci (nuclear vs. organellar), complicating phylogeographic reconstruction. The taxonomy of the ecologically specialized Abert’s squirrel species group has been controversial, and previous studies based on mitochondrial data have not fully resolved the evolutionary relationships among populations. Moreover, while these studies identified potential areas of secondary contact between divergent lineages, the possibility for introgression has not been tested.

Results

We used RAD-seq to unravel the complex evolutionary history of the Abert’s squirrel species group. Although some of our findings reinforce inferences based on mitochondrial data, we also find significant areas of discordance. Discordant signals generally arise from previously undetected introgression between divergent populations that differentially affected variation at mitochondrial and nuclear loci. Most notably, our results support earlier claims (disputed by mitochondrial data) that S. aberti kaibabensis, found only on the north rim of the Grand Canyon, is highly divergent from other populations. However, we also detected introgression of S. aberti kaibabensis DNA into other S. aberti populations, which likely accounts for the previously inferred close genetic relationship between this population and those south of the Grand Canyon.

Conclusions

Overall, the evolutionary history of Abert’s squirrels appears to be shaped largely by divergence during periods of habitat isolation. However, we also found evidence for interbreeding during periods of secondary contact resulting in introgression, with variable effects on mitochondrial and nuclear markers. Our results support the emerging view that populations often diversify under scenarios involving both divergence in isolation and gene flow during secondary contact, and highlight the value of genome-wide datasets for resolving such complex evolutionary histories.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1248-4) contains supplementary material, which is available to authorized users.

Phylogeny of mitochondrial DNA clones in tassel-eared squirrels Sciurus aberti

The tassel-eared squirrel, Sciurus aberti, includes six subspecies which occupy restrictive and apparently identical habitats in Ponderosa pine forests in the south-western United States and Mexico; the strict habitat requirement of this species is based on dietary requirements which are only fulfilled in these forests. To examine evolutionary relationships among certain subspecies of S. aberti, we obtained estimates of nucleotide diversity within subspecies as well as nucleotide divergence between subspecies using mitochondrial DNA (mtDNA) analysis. Restriction site polymorphisms were identified in samples of the four US subspecies: S. a. aberti (Abert), S. a. kaibabensis (Kaibab), S. a. ferreus (Ferreus), and S. a. chuscensis (Chuska) Fourteen mtDNA clones were resolved that were, with one exception, uniquely subspecific. Dendrograms constructed by neighbour-joining and maximum parsimony methods revealed two major assemblages: (1) an Abert/Kaibab group; and (2) a Ferreus/Chuska group. The Abert vs. Ferreus clones exhibited the greatest net nucleotide divergence, with a lineage separation estimate approximating 572,000 years ago assuming a nucleotide substitution rate of 7.15 x 10(-9)/year/site. Five out of ten Chuska squirrels shared a clone with one Abert sample; the relative sizes of these two populations and their respective ranges as well as their close proximity support the proposal for relatively recent intermixing of Abert and Chuska populations resulting in what appears to be Abert-->Chuska migration. Nucleotide diversity within subspecies ranked as Kaibab < Ferreus < Abert < Chuska; the relatively high diversity for the Chuska sample is based on the apparent introgression of Abert mtDNA. The relative diversity exhibited by Kaibab, Ferreus and Aberti samples corresponds to the range size of the respective subspecies.

T-cell receptor genes in tassel-eared squirrels (Sciurus aberti). I. Genetic polymorphism and divergence in the Abert and Kaibab subspecies

The role of environmental factors in the evolution and maintenance of diversity of antigen receptor gene families which participate in the immune response in mammals is inadequately understood. In order to elucidate the impact of these factors, we have undertaken the analysis of these gene families in the tassel-eared squirrel (Sciurus aberti) which has been separated into discrete subspecies by geographic barriers and whose food resources can be quantitated for estimating environmental quality. In this communication we describe the initial analysis of the complexity and polymorphism of sequences related to T-cell receptor (Tcr) alpha and beta chain genes in two subspecies, Sciurus aberti aberti (Abert) and Sciurus aberti kaibabensis (Kaibab) which have identical habitats and are separated by the Grand Canyon in Arizona, USA. Genomic blot analysis of 60 Abert and 62 Kaibab individuals collected over a 3-year period was performed with mouse Tcrb and Tcra cDNA probes. Sequences homologous to Tcrb-C, Tcrb-J1, and Tcrb-J2 genes were observed in all individuals from both subspecies; although Tcrb-J1 fragments were monomorphic. Tcrb-C and Tcrb-J2 fragments were polymorphic with both species- and subspecies-specific sequences. A single, monomorphic Tcra-C fragment was observed in addition to multiple Tcra-V fragments homologous to the mouse Tcra-V1 subfamily. Abert samples exhibited greater numbers of Tcra-V1 fragments as well as greater polymorphism than Kaibab samples. Heterozygosity estimates of Tcrb-C and Tcra-V1 sequences were determined for annually collected samples and compared with the yearly estimates of availability of hypogeous fungi, one of the major diet items of tassel-eared squirrels. In the Kaibab annual collections, Tcra-V1 heterozygosity declined with the decline in food resource, whereas heterozygosity of Tcrb-C sequences was inversely related to food resource. Similarly, a reduction in food resource for Abert squirrels in 1985 coincided with an increase in Tcrb-C heterozygosity in the same year. These results suggest that the diversity of gene families which participate in the immune response in mammals may be affected by environmental factors.

The major histocompatibility complex of tassel-eared squirrels. I. Genetic diversity associated with Kaibab squirrels

The complexity and polymorphism of sequences related to the class I and class II genes of mammalian major histocompatibility complexes (MHCs) were investigated in the tassel-eared squirrel subspecies Sciurus aberti kaibabensis or Kaibab squirrel. Kaibab squirrels are geographically isolated on the Kaibab plateau north of the Grand Canyon in Arizona. Genomic DNA from 22 individuals was digested with Eco RI and Bam HI, electrophoresed, blotted, and hybridized with a panel of human class I and class II probes. Sequences homologous to DR alpha, DR beta, DQ beta probes were observed. A single, nonpolymorphic DR alpha-related sequence and multiple, polymorphic DQ alpha-related sequences were observed. Hybridization with DR beta and DQ beta probes revealed multiple, polymorphic sequences with such specificity that no bands were observed to hybridize with both probes. The level of polymorphism of beta sequences exceeded that observed with alpha sequences. Further, three Eco RI bands apparently included at least parts of both alpha and beta sequences. Hybridization of genomic blots with the HLA-B7 class I probe revealed a number of bands comparable in size range and number to other mammalian species. However, only a minor percentage of bands were observed to segregate. The inheritance of these five families of sequences appeared to be neither concordant nor random in the sample population. Based on prior conclusions in other species, these class I and class II sequences are presumed to map to the Kabib MHC, TLSA. Although DQ alpha- and DQ beta-related sequences were concordantly inherited, segregating sequences in the other families could not be assigned to identifiable, segregating haplotypes. These observations suggest that the present-day TSLA haplotypes have been derived from a limited number of progenitor haplotypes through repeated, intra-TSLA recombination.