Different complex regulatory phenotypes underlie hybrid male sterility in divergent rodent crosses.

Hybrid incompatibilities are a critical component of species barriers and may arise due to negative interactions between divergent regulatory elements in parental species. We used a comparative approach to identify common themes in the regulatory phenotypes associated with hybrid male sterility in two divergent rodent crosses, dwarf hamsters and house mice. We investigated three potential characteristic regulatory phenotypes in hybrids including the propensity towards over or underexpression relative to parental species, the influence of developmental stage on the extent of misexpression, and the role of the sex chromosomes on misexpression phenotypes. In contrast to near pervasive overexpression in hybrid house mice, we found that misexpression in hybrid dwarf hamsters was dependent on developmental stage. In both house mouse and dwarf hamster hybrids, however, misexpression increased with the progression of spermatogenesis, although to varying extents and with potentially different consequences. In both systems, we detected sex-chromosome specific overexpression in stages of spermatogenesis where inactivated X chromosome expression was expected, but the hybrid overexpression phenotypes were fundamentally different. Importantly, misexpression phenotypes support the presence of multiple histological blocks to spermatogenesis in dwarf hamster hybrids, including a potential role of meiotic stalling early in spermatogenesis. Collectively, we demonstrate that while there are some similarities in hybrid regulatory phenotypes of house mice and dwarf hamsters, there are also clear differences that point towards unique mechanisms underlying hybrid male sterility in each system. Our results highlight the potential of comparative approaches in helping to understand the importance of disrupted gene regulation in speciation.

The Evolution of Widespread Recombination Suppression on the Dwarf Hamster (Phodopus) X Chromosome.

The X chromosome of therian mammals shows strong conservation among distantly related species, limiting insights into the distinct selective processes that have shaped sex chromosome evolution. We constructed a chromosome-scale de novo genome assembly for the Siberian dwarf hamster (Phodopus sungorus), a species reported to show extensive recombination suppression across an entire arm of the X chromosome. Combining a physical genome assembly based on shotgun and long-range proximity ligation sequencing with a dense genetic map, we detected widespread suppression of female recombination across ∼65% of the Phodopus X chromosome. This region of suppressed recombination likely corresponds to the Xp arm, which has previously been shown to be highly heterochromatic. Using additional sequencing data from two closely related species (P. campbelli and P. roborovskii), we show that recombination suppression on Xp appears to be independent of major structural rearrangements. The suppressed Xp arm was enriched for several transposable element families and de-enriched for genes primarily expressed in placenta, but otherwise showed similar gene densities, expression patterns, and rates of molecular evolution when compared to the recombinant Xq arm. Phodopus Xp gene content and order was also broadly conserved relative to the more distantly related rat X chromosome. These data suggest that widespread suppression of recombination has likely evolved through the transient induction of facultative heterochromatin on the Phodopus Xp arm without major changes in chromosome structure or genetic content. Thus, substantial changes in the recombination landscape have so far had relatively subtle influences on patterns of X-linked molecular evolution in these species.

Genetic divergence among threespine stickleback that differ in nuptial coloration.

Sexual signals are shaped by their intended and unintended receivers as well as the signalling environment. This interplay between sexual and natural selection can lead to divergence in signals in heterogeneous environments. Yet, the extent to which gene flow is restricted when signalling phenotypes vary across environments and over what spatial scales remains an outstanding question. In this study, we quantify gene flow between two colour morphs, red and black, of freshwater threespine stickleback fish (Gasterosteus aculeatus). We capitalize on the very recent divergence of signalling phenotypes in this system to characterize within-species and among-morph genetic variation and to test for levels of gene flow between colour morphs in Oregon and Washington. Despite limited evidence for assortative mating between allopatric red and black populations, we found that black populations are genetically distinct from nearby red populations and that the black morph appears to have evolved independently at least twice in Oregon and Washington. Surprisingly, we uncovered a group of stickleback in one small coastal stream, Connor Creek, which is genetically and morphologically distinct from the red and black colour morphs and from marine stickleback. Historically, both colour morphs have coexisted in this location and sometimes hybridized, raising new questions about the origins and history of these fish, which were first described as anadromous-black hybrids >50 years ago. Understanding how genetic variation is currently partitioned within and among populations and colour morphs in this system should prompt future studies to assess the relative roles of habitat, ecological and pre- and post-reproductive barriers in the genetic divergence and phenotypic patterns we observe in nature.

Stage-specific disruption of X chromosome expression during spermatogenesis in sterile house mouse hybrids.

Hybrid sterility is a complex phenotype that can result from the breakdown of spermatogenesis at multiple developmental stages. Here, we disentangle two proposed hybrid male sterility mechanisms in the house mice, Mus musculus domesticus and M. m. musculus, by comparing patterns of gene expression in sterile F1 hybrids from a reciprocal cross. We found that hybrid males from both cross directions showed disrupted X chromosome expression during prophase of meiosis I consistent with a loss of meiotic sex chromosome inactivation (MSCI) and Prdm9-associated sterility, but that the degree of disruption was greater in mice with an M. m. musculus X chromosome consistent with previous studies. During postmeiotic development, gene expression on the X chromosome was only disrupted in one cross direction, suggesting that misexpression at this later stage was genotype-specific and not a simple downstream consequence of MSCI disruption which was observed in both reciprocal crosses. Instead, disrupted postmeiotic expression may depend on the magnitude of earlier disrupted MSCI, or the disruption of particular X-linked genes or gene networks. Alternatively, only hybrids with a potential deficit of Sly copies, a Y-linked ampliconic gene family, showed overexpression in postmeiotic cells, consistent with a previously proposed model of antagonistic coevolution between the X- and Y-linked ampliconic genes contributing to disrupted expression late in spermatogenesis. The relative contributions of these two regulatory mechanisms and their impact on sterility phenotypes await further study. Our results further support the hypothesis that X-linked hybrid sterility in house mice has a variable genetic basis, and that genotype-specific disruption of gene regulation contributes to overexpression of the X chromosome at different stages of development. Overall, these findings underscore the critical role of epigenetic regulation of the X chromosome during spermatogenesis and suggest that these processes are prone to disruption in hybrids.

Unraveling patterns of disrupted gene expression across a complex tissue.

Whole tissue RNASeq is the standard approach for studying gene expression divergence in evolutionary biology and provides a snapshot of the comprehensive transcriptome for a given tissue. However, whole tissues consist of diverse cell types differing in expression profiles, and the cellular composition of these tissues can evolve across species. Here, we investigate the effects of different cellular composition on whole tissue expression profiles. We compared gene expression from whole testes and enriched spermatogenesis populations in two species of house mice, Mus musculus musculus and M. m. domesticus, and their sterile and fertile F1 hybrids, which differ in both cellular composition and regulatory dynamics. We found that cellular composition differences skewed expression profiles and differential gene expression in whole testes samples. Importantly, both approaches were able to detect large-scale patterns such as disrupted X chromosome expression, although whole testes sampling resulted in decreased power to detect differentially expressed genes. We encourage researchers to account for histology in RNASeq and consider methods that reduce sample complexity whenever feasible. Ultimately, we show that differences in cellular composition between tissues can modify expression profiles, potentially altering inferred gene ontological processes, insights into gene network evolution, and processes governing gene expression evolution.

Streetlights positively affect the presence of an invasive grass species.

Anthropogenic disturbances associated with urban ecosystems can create favorable conditions for populations of some invasive plant species. Light pollution is one of these disturbances, but how it affects the growth and establishment of invasive plant populations is unknown. Cheatgrass (Bromus tectorum) is a problematic invasive species where it has displaced native grassland communities in the United States, but to our knowledge, there have been no studies of the ecological factors that affect cheatgrass presence in urban ecosystems. We conducted field surveys in urban alleys in Denver, Colorado, to compare the presence of cheatgrass at sites with and without artificial light at night (hereafter artificial light) from streetlights. These streetlights are mounted on utility poles, which cause ground disturbance when installed in alleys; we were able to test the independent effect of poles on cheatgrass establishment because not all poles have streetlights on them. We found that cheatgrass was positively associated with the presence of streetlights and to a lesser extent poles. In addition to cheatgrass, we also found that other plants were positively associated with the presence of both poles and streetlights. Our results suggest that artificial light may benefit the occurrence of cheatgrass and other plant species in urban settings. While invasive populations of cheatgrass in wild habitats attract the most attention from managers, we suggest more consideration for this grass in urban environments where its growth and establishment benefit from anthropogenic changes.

Pedigree-based and phylogenetic methods support surprising patterns of mutation rate and spectrum in the gray mouse lemur.

Mutations are the raw material on which evolution acts, and knowledge of their frequency and genomic distribution is crucial for understanding how evolution operates at both long and short timescales. At present, the rate and spectrum of de novo mutations have been directly characterized in relatively few lineages. Our study provides the first direct mutation-rate estimate for a strepsirrhine (i.e., the lemurs and lorises), which comprises nearly half of the primate clade. Using high-coverage linked-read sequencing for a focal quartet of gray mouse lemurs (Microcebus murinus), we estimated the mutation rate to be among the highest calculated for a mammal at 1.52 × 10-8 (95% credible interval: 1.28 × 10-8-1.78 × 10-8) mutations/site/generation. Further, we found an unexpectedly low count of paternal mutations, and only a modest overrepresentation of mutations at CpG sites. Despite the surprising nature of these results, we found both the rate and spectrum to be robust to the manipulation of a wide range of computational filtering criteria. We also sequenced a technical replicate to estimate a false-negative and false-positive rate for our data and show that any point estimate of a de novo mutation rate should be considered with a large degree of uncertainty. For validation, we conducted an independent analysis of context-dependent substitution types for gray mouse lemur and five additional primate species for which de novo mutation rates have also been estimated. These comparisons revealed general consistency of the mutation spectrum between the pedigree-based and the substitution-rate analyses for all species compared.

Cryptic Patterns of Speciation in Cryptic Primates: Microendemic Mouse Lemurs and the Multispecies Coalescent.

Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic revision has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we investigate mouse lemur diversity in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for two pairs of sister lineages that include three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among the two pairs of lineages. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species, M. mittermeieri and M. lehilahytsara is not supported-a result that is particularly striking when using the genealogical discordance index (gdi). Nonsister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated relaxed-clock estimates. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case. [Cryptic species; effective population size; microendemism; multispecies coalescent; speciation; species delimitation.].

Comparative Genomic Analysis of the Pheromone Receptor Class 1 Family (V1R) Reveals Extreme Complexity in Mouse Lemurs (Genus, Microcebus) and a Chromosomal Hotspot across Mammals.

Sensory gene families are of special interest for both what they can tell us about molecular evolution and what they imply as mediators of social communication. The vomeronasal type-1 receptors (V1Rs) have often been hypothesized as playing a fundamental role in driving or maintaining species boundaries given their likely function as mediators of intraspecific mate choice, particularly in nocturnal mammals. Here, we employ a comparative genomic approach for revealing patterns of V1R evolution within primates, with a special focus on the small-bodied nocturnal mouse and dwarf lemurs of Madagascar (genera Microcebus and Cheirogaleus, respectively). By doubling the existing genomic resources for strepsirrhine primates (i.e. the lemurs and lorises), we find that the highly speciose and morphologically cryptic mouse lemurs have experienced an elaborate proliferation of V1Rs that we argue is functionally related to their capacity for rapid lineage diversification. Contrary to a previous study that found equivalent degrees of V1R diversity in diurnal and nocturnal lemurs, our study finds a strong correlation between nocturnality and V1R elaboration, with nocturnal lemurs showing elaborate V1R repertoires and diurnal lemurs showing less diverse repertoires. Recognized subfamilies among V1Rs show unique signatures of diversifying positive selection, as might be expected if they have each evolved to respond to specific stimuli. Furthermore, a detailed syntenic comparison of mouse lemurs with mouse (genus Mus) and other mammalian outgroups shows that orthologous mammalian subfamilies, predicted to be of ancient origin, tend to cluster in a densely populated region across syntenic chromosomes that we refer to as a V1R "hotspot."

Conservation genomic analysis reveals ancient introgression and declining levels of genetic diversity in Madagascar's hibernating dwarf lemurs.

Madagascar's biodiversity is notoriously threatened by deforestation and climate change. Many of these organisms are rare, cryptic, and severely threatened, making population-level sampling unrealistic. Such is the case with Madagascar's dwarf lemurs (genus Cheirogaleus), the only obligate hibernating primate. We here apply comparative genomic approaches to generate the first genome-wide estimates of genetic diversity within dwarf lemurs. We generate a reference genome for the fat-tailed dwarf lemur, Cheirogaleus medius, and use this resource to facilitate analyses of high-coverage (~30×) genome sequences for wild-caught individuals representing species: C. sp. cf. medius, C. major, C. crossleyi, and C. sibreei. This study represents the largest contribution to date of novel genomic resources for Madagascar's lemurs. We find concordant phylogenetic relationships among the four lineages of Cheirogaleus across most of the genome, and yet detect a number of discordant genomic regions consistent with ancient admixture. We hypothesized that these regions could have resulted from adaptive introgression related to hibernation, indeed finding that genes associated with hibernation are present, though most significantly, that gene ontology categories relating to transcription are over-represented. We estimate levels of heterozygosity and find particularly low levels in an individual sampled from an isolated population of C. medius that we refer to as C. sp. cf. medius. Results are consistent with a recent decline in effective population size, which is evident across species. Our study highlights the power of comparative genomic analysis for identifying species and populations of conservation concern, as well as for illuminating possible mechanisms of adaptive phenotypic evolution.

Warning SINEs: Alu elements, evolution of the human brain, and the spectrum of neurological disease.

Alu elements are a highly successful family of primate-specific retrotransposons that have fundamentally shaped primate evolution, including the evolution of our own species. Alus play critical roles in the formation of neurological networks and the epigenetic regulation of biochemical processes throughout the central nervous system (CNS), and thus are hypothesized to have contributed to the origin of human cognition. Despite the benefits that Alus provide, deleterious Alu activity is associated with a number of neurological and neurodegenerative disorders. In particular, neurological networks are potentially vulnerable to the epigenetic dysregulation of Alu elements operating across the suite of nuclear-encoded mitochondrial genes that are critical for both mitochondrial and CNS function. Here, we highlight the beneficial neurological aspects of Alu elements as well as their potential to cause disease by disrupting key cellular processes across the CNS. We identify at least 37 neurological and neurodegenerative disorders wherein deleterious Alu activity has been implicated as a contributing factor for the manifestation of disease, and for many of these disorders, this activity is operating on genes that are essential for proper mitochondrial function. We conclude that the epigenetic dysregulation of Alu elements can ultimately disrupt mitochondrial homeostasis within the CNS. This mechanism is a plausible source for the incipient neuronal stress that is consistently observed across a spectrum of sporadic neurological and neurodegenerative disorders.

The Alu neurodegeneration hypothesis: A primate-specific mechanism for neuronal transcription noise, mitochondrial dysfunction, and manifestation of neurodegenerative disease.

It is hypothesized that retrotransposons have played a fundamental role in primate evolution and that enhanced neurologic retrotransposon activity in humans may underlie the origin of higher cognitive function. As a potential consequence of this enhanced activity, it is likely that neurons are susceptible to deleterious retrotransposon pathways that can disrupt mitochondrial function. An example is observed in the TOMM40 gene, encoding a ß-barrel protein critical for mitochondrial preprotein transport. Primate-specific Alu retrotransposons have repeatedly inserted into TOMM40 introns, and at least one variant associated with late-onset Alzheimer's disease originated from an Alu insertion event. We provide evidence of enriched Alu content in mitochondrial genes and postulate that Alus can disrupt mitochondrial populations in neurons, thereby setting the stage for progressive neurologic dysfunction. This Alu neurodegeneration hypothesis is compatible with decades of research and offers a plausible mechanism for the disruption of neuronal mitochondrial homeostasis, ultimately cascading into neurodegenerative disease.

Geogenetic patterns in mouse lemurs (genus Microcebus) reveal the ghosts of Madagascar's forests past.

Phylogeographic analysis can be described as the study of the geological and climatological processes that have produced contemporary geographic distributions of populations and species. Here, we attempt to understand how the dynamic process of landscape change on Madagascar has shaped the distribution of a targeted clade of mouse lemurs (genus Microcebus) and, conversely, how phylogenetic and population genetic patterns in these small primates can reciprocally advance our understanding of Madagascar's prehuman environment. The degree to which human activity has impacted the natural plant communities of Madagascar is of critical and enduring interest. Today, the eastern rainforests are separated from the dry deciduous forests of the west by a large expanse of presumed anthropogenic grassland savanna, dominated by the Family Poaceae, that blankets most of the Central Highlands. Although there is firm consensus that anthropogenic activities have transformed the original vegetation through agricultural and pastoral practices, the degree to which closed-canopy forest extended from the east to the west remains debated. Phylogenetic and population genetic patterns in a five-species clade of mouse lemurs suggest that longitudinal dispersal across the island was readily achieved throughout the Pleistocene, apparently ending at ∼55 ka. By examining patterns of both inter- and intraspecific genetic diversity in mouse lemur species found in the eastern, western, and Central Highland zones, we conclude that the natural environment of the Central Highlands would have been mosaic, consisting of a matrix of wooded savanna that formed a transitional zone between the extremes of humid eastern and dry western forest types.