Snowshoe hares (
) maintain seasonal camouflage by molting to a white winter coat, but some hares remain brown during the winter in regions with low snow cover. We show that cis-regulatory variation ...controlling seasonal expression of the
gene underlies this adaptive winter camouflage polymorphism. Genetic variation at
clustered by winter coat color across multiple hare and jackrabbit species, revealing a history of recurrent interspecific gene flow. Brown winter coats in snowshoe hares likely originated from an introgressed black-tailed jackrabbit allele that has swept to high frequency in mild winter environments. These discoveries show that introgression of genetic variants that underlie key ecological traits can seed past and ongoing adaptation to rapidly changing environments.
Hybridization may often be an important source of adaptive variation, but the extent and long-term impacts of introgression have seldom been evaluated in the phylogenetic context of a radiation. ...Hares (Lepus) represent a widespread mammalian radiation of 32 extant species characterized by striking ecological adaptations and recurrent admixture. To understand the relevance of introgressive hybridization during the diversification of Lepus, we analyzed whole exome sequences (61.7 Mb) from 15 species of hares (1-4 individuals per species), spanning the global distribution of the genus, and two outgroups. We used a coalescent framework to infer species relationships and divergence times, despite extensive genealogical discordance. We found high levels of allele sharing among species and show that this reflects extensive incomplete lineage sorting and temporally layered hybridization. Our results revealed recurrent introgression at all stages along the Lepus radiation, including recent gene flow between extant species since the last glacial maximum but also pervasive ancient introgression occurring since near the origin of the hare lineages. We show that ancient hybridization between northern hemisphere species has resulted in shared variation of potential adaptive relevance to highly seasonal environments, including genes involved in circadian rhythm regulation, pigmentation, and thermoregulation. Our results illustrate how the genetic legacy of ancestral hybridization may persist across a radiation, leaving a long-lasting signature of shared genetic variation that may contribute to adaptation. Adaptation; ancient introgression; hybridization; Lepus; phylogenomics..
Epigenetic modifications confer stable transcriptional patterns in the brain, and both normal and abnormal brain function involve specialized brain regions. We examined DNA methylation by ...whole-genome bisulfite sequencing in neuronal and non-neuronal populations from four brain regions (anterior cingulate gyrus, hippocampus, prefrontal cortex, and nucleus accumbens) as well as chromatin accessibility in the latter two. We find pronounced differences in both CpG and non-CpG methylation (CG-DMRs and CH-DMRs) only in neuronal cells across brain regions. Neuronal CH-DMRs were highly associated with differential gene expression, whereas CG-DMRs were consistent with chromatin accessibility and enriched for regulatory regions. These CG-DMRs comprise ~12 Mb of the genome that is highly enriched for genomic regions associated with heritability of neuropsychiatric traits including addictive behavior, schizophrenia, and neuroticism, thus suggesting a mechanistic link between pathology and differential neuron-specific epigenetic regulation in distinct brain regions.
DNA methylation dynamics in the brain are associated with normal development and neuropsychiatric disease and differ across functionally distinct brain regions. Previous studies of genome-wide ...methylation differences among human brain regions focus on limited numbers of individuals and one to two brain regions.
Using GTEx samples, we generate a resource of DNA methylation in purified neuronal nuclei from 8 brain regions as well as lung and thyroid tissues from 12 to 23 donors. We identify differentially methylated regions between brain regions among neuronal nuclei in both CpG (181,146) and non-CpG (264,868) contexts, few of which were unique to a single pairwise comparison. This significantly expands the knowledge of differential methylation across the brain by 10-fold. In addition, we present the first differential methylation analysis among neuronal nuclei from basal ganglia tissues and identify unique CpG differentially methylated regions, many associated with ion transport. We also identify 81,130 regions of variably CpG methylated regions, i.e., variable methylation among individuals in the same brain region, which are enriched in regulatory regions and in CpG differentially methylated regions. Many variably methylated regions are unique to a specific brain region, with only 202 common across all brain regions, as well as lung and thyroid. Variably methylated regions identified in the amygdala, anterior cingulate cortex, and hippocampus are enriched for heritability of schizophrenia.
These data suggest that epigenetic variation in these particular human brain regions could be associated with the risk for this neuropsychiatric disorder.
Aim Quaking aspen (Populus tremuloides) has the largest natural distribution of any tree native to North America. The primary objectives of this study were to characterize range-wide genetic ...diversity and genetic structuring in quaking aspen, and to assess the influence of glacial history and rear-edge dynamics. Location North America. Methods Using a sample set representing the full longitudinal and latitudinal extent of the species' distribution, we examined geographical patterns of genetic diversity and structuring using 8 nuclear microsatellite loci in 794 individuals from 30 sampling sites. Results Two major genetic clusters were identified across the range: a south-western cluster and a northern cluster. The south-western cluster, which included two subclusters, was bounded approximately by the Continental Divide to the east and the southern extent of the ice sheet at the Last Glacial Maximum to the north. Subclusters were not detected in the northern cluster, despite its continent-wide distribution. Genetic distance was significantly correlated with geographical distance in the south-western but not the northern cluster, and allelic richness was significantly lower in south-western sampling sites compared with northern sampling sites. Population structuring was low overall, but elevated in the south-western cluster. Main conclusions Aspen populations in the south-western portion of the range are consistent with expectations for a historically stable edge, with low within-population diversity, significant geographical population structuring, and little evidence of northward expansion. Structuring within the south-western cluster may result from distinct gene pools separated during the Pleistocene and reunited following glacial retreat, similar to patterns found in other forest tree species in the western USA. In aspen, populations in the south-western portion of the species range are thought to be at particularly high risk of mortality with climate change. Our findings suggest that these same populations may be disproportionately valuable in terms of both evolutionary potential and conservation value.
We document high rates of triploidy in aspen (Populus tremuloides) across the western USA (up to 69% of genets), and ask whether the incidence of triploidy across the species range corresponds with ...latitude, glacial history (as has been documented in other species), climate, or regional variance in clone size. Using a combination of microsatellite genotyping, flow cytometry, and cytology, we demonstrate that triploidy is highest in unglaciated, drought-prone regions of North America, where the largest clone sizes have been reported for this species. While we cannot completely rule out a low incidence of undetected aneuploidy, tetraploidy or duplicated loci, our evidence suggests that these phenomena are unlikely to be significant contributors to our observed patterns. We suggest that the distribution of triploid aspen is due to a positive synergy between triploidy and ecological factors driving clonality. Although triploids are expected to have low fertility, they are hypothesized to be an evolutionary link to sexual tetraploidy. Thus, interactions between clonality and polyploidy may be a broadly important component of geographic speciation patterns in perennial plants. Further, cytotypes are expected to show physiological and structural differences which may influence susceptibility to ecological factors such as drought, and we suggest that cytotype may be a significant and previously overlooked factor in recent patterns of high aspen mortality in the southwestern portion of the species range. Finally, triploidy should be carefully considered as a source of variance in genomic and ecological studies of aspen, particularly in western U.S. landscapes.
Abstract
Embryonic development in mammals is highly sensitive to changes in gene expression within the placenta. The placenta is also highly enriched for genes showing parent-of-origin or imprinted ...expression, which is predicted to evolve rapidly in response to parental conflict. However, little is known about the evolution of placental gene expression, or if divergence of placental gene expression plays an important role in mammalian speciation. We used crosses between two species of dwarf hamsters (Phodopus sungorus and Phodopus campbelli) to examine the genetic and regulatory underpinnings of severe placental overgrowth in their hybrids. Using quantitative genetic mapping and mitochondrial substitution lines, we show that overgrowth of hybrid placentas was primarily caused by genetic differences on the maternally inherited P. sungorus X chromosome. Mitochondrial interactions did not contribute to abnormal hybrid placental development, and there was only weak correspondence between placental disruption and embryonic growth. Genome-wide analyses of placental transcriptomes from the parental species and first- and second-generation hybrids revealed a central group of co-expressed X-linked and autosomal genes that were highly enriched for maternally biased expression. Expression of this gene network was strongly correlated with placental size and showed widespread misexpression dependent on epistatic interactions with X-linked hybrid incompatibilities. Collectively, our results indicate that the X chromosome is likely to play a prominent role in the evolution of placental gene expression and the accumulation of hybrid developmental barriers between mammalian species.
Discovery of cryptic species is essential to understand the process of speciation and assessing the impacts of anthropogenic stressors. Here, we used genomic data to test for cryptic species ...diversity within an ecologically well‐known radiation of North American rodents, western chipmunks (Tamias). We assembled a de novo reference genome for a single species (Tamias minimus) combined with new and published targeted sequence‐capture data for 21,551 autosomal and 493 X‐linked loci sampled from 121 individuals spanning 22 species. We identified at least two cryptic lineages corresponding with an isolated subspecies of least chipmunk (T. minimus grisescens) and with a restricted subspecies of the yellow‐pine chipmunk (Tamias amoenus cratericus) known only from around the extensive Craters of the Moon lava flow. Additional population‐level sequence data revealed that the so‐called Crater chipmunk is a distinct species that is abundant throughout the coniferous forests of southern Idaho. This cryptic lineage does not appear to be most closely related to the ecologically and phenotypically similar yellow‐pine chipmunk but does show evidence for recurrent hybridization with this and other species.
Abstract
Adaptive radiations are characterized by the diversification and ecological differentiation of species, and replicated cases of this process provide natural experiments for understanding the ...repeatability and pace of molecular evolution. During adaptive radiation, genes related to ecological specialization may be subject to recurrent positive directional selection. However, it is not clear to what extent patterns of lineage-specific ecological specialization (including phenotypic convergence) are correlated with shared signatures of molecular evolution. To test this, we sequenced whole exomes from a phylogenetically dispersed sample of 38 murine rodent species, a group characterized by multiple, nested adaptive radiations comprising extensive ecological and phenotypic diversity. We found that genes associated with immunity, reproduction, diet, digestion, and taste have been subject to pervasive positive selection during the diversification of murine rodents. We also found a significant correlation between genome-wide positive selection and dietary specialization, with a higher proportion of positively selected codon sites in derived dietary forms (i.e., carnivores and herbivores) than in ancestral forms (i.e., omnivores). Despite striking convergent evolution of skull morphology and dentition in two distantly related worm-eating specialists, we did not detect more genes with shared signatures of positive or relaxed selection than in a nonconvergent species comparison. Although a small number of the genes we detected can be incidentally linked to craniofacial morphology or diet, protein-coding regions are unlikely to be the primary genetic basis of this complex convergent phenotype. Our results suggest a link between positive selection and derived ecological phenotypes, and highlight specific genes and general functional categories that may have played an integral role in the extensive and rapid diversification of murine rodents.
Color molts from summer brown to winter white coats have evolved in several species to maintain camouflage year‐round in environments with seasonal snow. Despite the eco‐evolutionary relevance of ...this key phenological adaptation, its molecular regulation has only recently begun to be addressed. Here, we analyze skin transcription changes during the autumn molt of the mountain hare (Lepus timidus) and integrate the results with an established model of gene regulation across the spring molt of the closely related snowshoe hare (L. americanus). We quantified differences in gene expression among three stages of molt progression—“brown” (early molt), “intermediate,” and “white” (late molt). We found 632 differentially expressed genes, with a major pulse of expression early in the molt, followed by a milder one in late molt. The functional makeup of differentially expressed genes anchored the sampled molt stages to the developmental timeline of the hair growth cycle, associating anagen to early molt and the transition to catagen to late molt. The progression of color change was characterized by differential expression of genes involved in pigmentation, circadian, and behavioral regulation. We found significant overlap between differentially expressed genes across the seasonal molts of mountain and snowshoe hares, particularly at molt onset, suggesting conservatism of gene regulation across species and seasons. However, some discrepancies suggest seasonal differences in melanocyte differentiation and the integration of nutritional cues. Our established regulatory model of seasonal coat color molt provides an important mechanistic context to study the functional architecture and evolution of this crucial seasonal adaptation.
We analyze skin transcription changes during the autumn color molt of the mountain hare and integrate the results with an established model of gene regulation across the spring molt of the closely related snowshoe hare. Our work provides an important mechanistic context to study the functional architecture and evolution of this crucial seasonal adaptation.