Positional cloning of the Hybrid sterility 1 gene: fine genetic mapping and evaluation of two candidate genes

Genetic Dissection of Hybrid Male Sterility Across Stages of Spermatogenesis

Hybrid sterility is a common form of reproductive isolation between nascent species. Although hybrid sterility is routinely documented and genetically dissected in speciation studies, its developmental basis is rarely examined, especially in generations beyond the F1 generation. To identify phenotypic and genetic determinants of hybrid male sterility from a developmental perspective, we characterized testis histology in 312 F2 hybrids generated by intercrossing inbred strains of Mus musculus domesticus and M. m. musculus, two subspecies of house mice. Hybrids display a range of histologic abnormalities that indicate defective spermatogenesis. Among these abnormalities, we quantified decreased testis size, reductions in spermatocyte and spermatid number, increased apoptosis of meiosis I spermatocytes, and more multinucleated syncytia. Collectively, our phenotypic data point to defects in meiosis I as a primary barrier to reproduction. We identified seven quantitative trait loci (QTL) controlling five histologic traits. A region of chromosome 17 that contains Prdm9, a gene known to confer F1 hybrid male sterility, affects multinucleated syncytia and round spermatids, potentially extending the phenotypic outcomes of this incompatibility. The X chromosome also plays a key role, with loci affecting multinucleated syncytia, apoptosis of round spermatids, and round spermatid numbers. We detected an epistatic interaction between QTL on chromosomes 17 and X for multinucleated syncytia. Our results refine the developmental basis of a key reproductive barrier in a classic model system for speciation genetics.

Gene Regulation and Speciation

Understanding the genetic architecture of speciation is a major goal in evolutionary biology. Hybrid dysfunction is thought to arise most commonly through negative interactions between alleles at two or more loci. Divergence between interacting regulatory elements that affect gene expression (i.e., regulatory divergence) may be a common route for these negative interactions to arise. We review here how regulatory divergence between species can result in hybrid dysfunction, including recent theoretical support for this model. We then discuss the empirical evidence for regulatory divergence between species and evaluate evidence for misregulation as a source of hybrid dysfunction. Finally, we review unresolved questions in gene regulation as it pertains to speciation and point to areas that could benefit from future research.

Interallelic and intergenic incompatibilities of the Prdm9 (Hst1) gene in mouse hybrid sterility

The Dobzhansky-Muller model of incompatibilities explains reproductive isolation between species by incorrect epistatic interactions. Although the mechanisms of speciation are of great interest, no incompatibility has been characterized at the gene level in mammals. The Hybrid sterility 1 gene (Hst1) participates in the arrest of meiosis in F₁ males of certain strains from two Mus musculus subspecies, e.g., PWD from M. m. musculus and C57BL/6J (henceforth B6) from M. m. domesticus. Hst1 has been identified as a meiotic PR-domain gene (Prdm9) encoding histone 3 methyltransferase in the male offspring of PWD females and B6 males, (PWD×B6)F₁. To characterize the incompatibilities underlying hybrid sterility, we phenotyped reproductive and meiotic markers in males with altered copy numbers of Prdm9. A partial rescue of fertility was observed upon removal of the B6 allele of Prdm9 from the azoospermic (PWD×B6)F₁ hybrids, whereas removing one of the two Prdm9 copies in PWD or B6 background had no effect on male reproduction. Incompatibility(ies) not involving Prdm9^B6 also acts in the (PWD×B6)F₁ hybrids, since the correction of hybrid sterility by Prdm9^B6 deletion was not complete. Additions and subtractions of Prdm9 copies, as well as allelic replacements, improved meiotic progression and fecundity also in the progeny-producing reciprocal (B6×PWD)F₁ males. Moreover, an increased dosage of Prdm9 and reciprocal cross enhanced fertility of other sperm-carrying male hybrids, (PWD×B6-C3H.Prdm9)F₁, harboring another Prdm9 allele of M. m. domesticus origin. The levels of Prdm9 mRNA isoforms were similar in the prepubertal testes of all types of F₁ hybrids of PWD with B6 and B6-C3H.Prdm9 despite their different prospective fertility, but decreased to 53% after removal of Prdm9^B6. Therefore, the Prdm9^B6 allele probably takes part in posttranscriptional dominant-negative hybrid interaction(s) absent in the parental strains.

Disturbed gametogenesis in the progeny of two fertile parental forms is called hybrid sterility; it is an important part of reproductive barriers between species. The Dobzhansky-Muller model of incompatibilities explains reproductive isolation between species by incorrect interactions between genes. Hybrid sterility 1 (Hst1) is one of the genes causing meiotic arrest in F₁ male hybrids between certain Mus musculus musculus (e.g., the PWD strain) and M. m. domesticus (C57BL/6J etc.) mice. Hst1, the first mammalian candidate for a speciation gene, was identified as a meiotic PR/SET-domain gene, Prdm9, but the mechanism causing sterility has remained unknown. While the F₁ male offspring of C57BL/6J males and PWD females produce no sperm, the males from the reciprocal cross using PWD males and C57BL/6J females yield progeny. Here we show that the meiotic progress and fertility of hybrid males from both F₁ crosses improved by removal as well as overexpression of the C57BL/6J allele of Prdm9, suggesting that Prdm9 interactions not present in the parental species (incompatibilities) play a role in hybrid sterility. Furthermore, the Prdm9 dosage also controlled fecundity in other F₁ hybrids, indicating that this gene is an important regulator of mouse hybrid fertility.

The microevolution of V1r vomeronasal receptor genes in mice

Vomeronasal sensitivity is important for detecting intraspecific pheromonal cues as well as environmental odorants and is involved in mating, social interaction, and other daily activities of many vertebrates. Two large families of seven-transmembrane G-protein-coupled receptors, V1rs and V2rs, bind to various ligands to initiate vomeronasal signal transduction. Although the macroevolution of V1r and V2r genes has been well characterized throughout vertebrates, especially mammals, little is known about their microevolutionary patterns, which hampers a clear understanding of the evolutionary forces behind the rapid evolutionary turnover of V1r and V2r genes and the great diversity in receptor repertoire across species. Furthermore, the role of divergent vomeronasal perception in enhancing premating isolation and maintaining species identity has not been evaluated. Here we sequenced 44 V1r genes and 25 presumably neutral noncoding regions in 14 wild-caught mice belonging to Mus musculus and M. domesticus, two closely related species with strong yet incomplete reproductive isolation. We found that nucleotide changes in V1rs are generally under weak purifying selection and that only ∼5% of V1rs may have been subject to positive selection that promotes nonsynonymous substitutions. Consistent with the low functional constraints on V1rs, 18 of the 44 V1rs have null alleles segregating in one or both species. Together, our results demonstrate that, despite occasional actions of positive selection, the evolution of V1rs is in a large part shaped by purifying selection and random drift. These findings have broad implications for understanding the driving forces of rapid gene turnovers that are often observed in the evolution of large gene families.

A mouse speciation gene encodes a meiotic histone H3 methyltransferase

Speciation genes restrict gene flow between the incipient species and related taxa. Three decades ago, we mapped a mammalian speciation gene, hybrid sterility 1 (Hst1), in the intersubspecific hybrids of house mouse. Here, we identify this gene as Prdm9, encoding a histone H3 lysine 4 trimethyltransferase. We rescued infertility in male hybrids with bacterial artificial chromosomes carrying Prdm9 from a strain with the "fertility" Hst1(f) allele. Sterile hybrids display down-regulated microrchidia 2B (Morc2b) and fail to compartmentalize gammaH2AX into the pachynema sex (XY) body. These defects, seen also in Prdm9-null mutants, are rescued by the Prdm9 transgene. Identification of a vertebrate hybrid sterility gene reveals a role for epigenetics in speciation and opens a window to a hybrid sterility gene network.

Mouse consomic strains: exploiting genetic divergence between Mus m. musculus and Mus m. domesticus subspecies

Consomic (chromosome substitution) strains (CSs) represent the most recent addition to the mouse genetic resources aimed to genetically analyze complex trait loci (QTLs). In this study, we report the development of a set of 28 mouse intersubspecific CSs. In each CS, we replaced a single chromosome of the C57BL/6J (B6) inbred strain (mostly Mus m. domesticus) with its homolog from the PWD/Ph inbred strain of the Mus m. musculus subspecies. These two progenitor subspecies diverged less than 1 million years ago and accumulated a large number of genetic differences that constitute a rich resource of genetic variation for QTL analyses. Altogether, the 18 consomic, nine subconsomic, and one conplastic strain covered all 19 autosomes, X and Y sex chromosomes, and mitochondrial DNA. Most CSs had significantly lower reproductive fitness compared with the progenitor strains. CSs homosomic for chromosomes 10 and 11, and the C57BL/6J-Chr X males, failed to reproduce and were substituted by less affected subconsomics carrying either a proximal, central, or distal part of the respective chromosome. A genome-wide scan of 965 DNA markers revealed 99.87% purity of the B6 genetic background. Thirty-three nonsynonymous substitutions were uncovered in the protein-coding regions of the mitochondrial DNA of the B6.PWD-mt conplastic strain. A pilot-phenotyping experiment project revealed a high number of variations among B6.PWD consomics.

Fine haplotype structure of a chromosome 17 region in the laboratory and wild mouse

Extensive linkage disequilibrium among classical laboratory strains represents an obstacle in the high-resolution haplotype mapping of mouse quantitative trait loci (QTL). To determine the potential of wild-derived mouse strains for fine QTL mapping, we constructed a haplotype map of a 250-kb region of the t-complex on chromosome 17 containing the Hybrid sterility 1 (Hst1) gene. We resequenced 33 loci from up to 80 chromosomes of five mouse (sub)species. Trans-species single-nucleotide polymorphisms (SNPs) were rare between Mus m. musculus (Mmmu) and Mus m. domesticus (Mmd). The haplotypes in Mmmu and Mmd differed and therefore strains from these subspecies should not be combined for haplotype-associated mapping. The haplotypes of t-chromosomes differed from all non-t Mmmu and Mmd haplotypes. Half of the SNPs and SN indels but only one of seven longer rearrangements found in classical laboratory strains were useful for haplotype mapping in the wild-derived M. m. domesticus. The largest Mmd haplotype block contained three genes of a highly conserved synteny. The lengths of the haplotype blocks deduced from 36 domesticus chromosomes were in tens of kilobases, suggesting that the wild-derived Mmd strains are suitable for fine interval-specific mapping.

Asymmetry and polymorphism of hybrid male sterility during the early stages of speciation in house mice

House mice offer a powerful system for dissecting the genetic basis of phenotypes that isolate species in the early stages of speciation. We used a series of reciprocal crosses between wild-derived strains of Mus musculus and M. domesticus to examine F(1) hybrid male sterility, one of the primary phenotypes thought to isolate these species. We report four main results. First, we found significantly smaller testes and fewer sperm in hybrid male progeny of most crosses. Second, in some crosses hybrid male sterility was asymmetric and depended on the species origin of the X chromosome. These observations confirm and extend previous findings, underscoring the central role that the M. musculus X chromosome plays in reproductive isolation. Third, comparisons among reciprocal crosses revealed polymorphism at one or more hybrid incompatibilities within M. musculus. Fourth, the spermatogenic phenotype of this polymorphic interaction appears distinct from previously described hybrid incompatibilities between these species. These data build on previous studies of speciation in house mice and show that the genetic basis of hybrid male sterility is fairly complex, even at this early stage of divergence.

Genome-wide patterns of gene flow across a house mouse hybrid zone

Hybrid zones between closely related species or subspecies provide useful settings for studying the genetic architecture of speciation. Using markers distributed throughout the mouse genome, we use a hybrid zone between two recently diverged species of house mice (Mus musculus and Mus domesticus) as a natural mapping experiment to identify genomic regions that may be involved in reproductive isolation. Using cline analysis we document a nearly 50-fold variation in level of introgression among markers. Some markers have extremely narrow cline widths; these genomic regions may contribute to reproductive isolation. Biological processes associated with these narrow clines include physiological and immune responses to the environment as well as physiological and behavioral aspects of reproduction. Other autosomal markers exhibit asymmetrically broad clines, usually with high frequencies of M. domesticus alleles on the M. musculus side of the hybrid zone. These markers identify genome regions likely housing genes with alleles that are spreading from one species to the other. Biological processes associated with these wide clines include cell signaling, olfaction, and pheromone response. These processes play important roles in survival and reproduction, and associated genes are likely targets of selection. Patterns of linkage disequilibrium in the center of the hybrid zone suggest that isolation may be caused by multiple epistatic interactions between sets of genes. These data highlight the complex genetic architecture underlying speciation even at early stages of divergence and point to some of the biological processes that may govern this architecture.

On the subspecific origin of the laboratory mouse

The genome of the laboratory mouse is thought to be a mosaic of regions with distinct subspecific origins. We have developed a high-resolution map of the origin of the laboratory mouse by generating 25,400 phylogenetic trees at 100-kb intervals spanning the genome. On average, 92% of the genome is of Mus musculus domesticus origin, and the distribution of diversity is markedly nonrandom among the chromosomes. There are large regions of extremely low diversity, which represent blind spots for studies of natural variation and complex traits, and hot spots of diversity. In contrast with the mosaic model, we found that most of the genome has intermediate levels of variation of intrasubspecific origin. Finally, mouse strains derived from the wild that are supposed to represent different mouse subspecies show substantial intersubspecific introgression, which has strong implications for evolutionary studies that assume these are pure representatives of a given subspecies.

Conserved alternative and antisense transcripts at the programmed cell death 2 locus

Background

The programmed cell death 2 (Pdcd2) gene on mouse chromosome 17 was evaluated as a member of a highly conserved synteny, a candidate for an imprinted locus, and a candidate for the Hybrid sterility 1 (Hst1) gene.

Results

New mouse transcripts were identified at this locus: an alternative Pdcd2 mRNA skipping the last two coding exons and two classes of antisense RNAs. One class of the antisense RNA overlaps the alternative exon and the other the entire Pdcd2 gene. The antisense RNAs are alternative transcripts of the neighboring TATA-binding protein gene (Tbp) that are located mainly in the cell nucleus. Analogous alternative PDCD2 forms truncating the C-terminal domain were also detected in human and chicken. Alternative transcripts of the chicken PDCD2 and TBP genes also overlap. No correlation in the transcription of the alternative and overlapping mRNAs was detected. Allelic sequencing and transcription studies did not reveal any support for the candidacy of Pdcd2 for Hst1. No correlated expression of Pdcd2 with the other two genes of the highly conserved synteny was observed. Pdcd2, Chd1, and four other genes from this region were not imprinted in the embryo.

Conclusion

The conservation of alternative transcription of the Pdcd2 gene in mouse, human and chicken suggests the biological importance of such truncated protein. The biological function of the alternative PDCD2 is likely to be opposite to that of the constitutive form. The ratio of the constitutive and alternative Pdcd2 mRNAs differs in the tissues, suggesting a developmental role.

The identified Tbp-alternative Pdcd2-antisense transcripts may interfere with the transcription of the Pdcd2 gene, as they are transcribed at a comparable level. The conservation of the Pdcd2/Tbp sense-antisense overlap in the mouse and chicken points out its biological relevance. Our results also suggest that some cDNAs in databases labeled as noncoding are incomplete alternative cDNAs of neighboring protein-coding genes.

Searching the genomes of inbred mouse strains for incompatibilities that reproductively isolate their wild relatives

Identification of the genes that underlie reproductive isolation provides important insights into the process of speciation. According to the Dobzhansky-Muller model, these genes suffer disrupted interactions in hybrids due to independent divergence in separate populations. In hybrid populations, natural selection acts to remove the deleterious heterospecific combinations that cause these functional disruptions. When selection is strong, this process can maintain multilocus associations, primarily between conspecific alleles, providing a signature that can be used to locate incompatibilities. We applied this logic to populations of house mice that were formed by hybridization involving two species that show partial reproductive isolation, Mus domesticus and Mus musculus. Using molecular markers likely to be informative about species ancestry, we scanned the genomes of 1) classical inbred strains and 2) recombinant inbred lines for pairs of loci that showed extreme linkage disequilibria. By using the same set of markers, we identified a list of locus pairs that displayed similar patterns in both scans. These genomic regions may contain genes that contribute to reproductive isolation between M. domesticus and M. musculus. This hypothesis can now be tested using laboratory crosses and surveys of introgression in the wild.

Genomic islands of differentiation between house mouse subspecies

Understanding the genes that contribute to reproductive isolation is essential to understanding speciation, but isolating such genes has proven very difficult. In this study I apply a multilocus test statistic to >10,000 SNP markers assayed in wild-derived inbred strains of house mice to identify genomic regions of elevated differentiation between two subspecies of house mice, Mus musculus musculus and M. m. domesticus. Differentiation was high through approximately 90% of the X chromosome. In addition, eight regions of high differentiation were identified on the autosomes, totaling 7.5% of the autosomal genome. Regions of high differentiation were confirmed by direct sequencing of samples collected from the wild. Some regions of elevated differentiation have an overrepresentation of genes with host-pathogen interactions and olfaction. The most strongly differentiated region on the X has previously been shown to fail to introgress across a hybrid zone between the two subspecies. This survey indicates autosomal regions that should also be examined for differential introgression across the hybrid zone, as containing potential genes causing hybrid unfitness.

Signatures of reproductive isolation in patterns of single nucleotide diversity across inbred strains of mice

Reproductive isolation is often caused by the disruption of genic interactions that evolve in geographically separate populations. Identifying the genomic regions and genes involved in these interactions, known as "Dobzhansky-Muller incompatibilities," can be challenging but is facilitated by the wealth of genetic markers now available in model systems. In recent years, the complete genome sequence and thousands of single nucleotide polymorphisms (SNPs) from laboratory mice, which are largely genetic hybrids between Mus musculus and M. domesticus, have become available. Here, we use these resources to locate genomic regions that may underlie reproductive isolation between these two species. Using genotypes from 332 SNPs that differ between wild-derived strains of M. musculus and M. domesticus, we identified several physically unlinked SNP pairs that show exceptional gametic disequilibrium across the lab strains. Conspecific alleles were associated in a disproportionate number of these cases, consistent with the action of natural selection against hybrid gene combinations. As predicted by the Dobzhansky-Muller model, this bias was differentially attributable to locus pairs for which one hybrid genotype was missing. We assembled a list of potential Dobzhansky-Muller incompatibilities from locus pairs that showed extreme associations (only three gametic types) among conspecific alleles. Two SNPs in this list map near known hybrid sterility loci on chromosome 17 and the X chromosome, allowing us to nominate partners for disrupted interactions involving these genomic regions for the first time. Together, these results indicate that patterns produced by speciation between M. musculus and M. domesticus are visible in the genomes of lab strains of mice, underscoring the potential of these genetic model organisms for addressing general questions in evolutionary biology.