Petr Simecek

Defining the consequences of genetic variation on a proteome-wide scale

Genetic variation modulates protein expression through both transcriptional and post-transcriptional mechanisms. To characterize the consequences of natural genetic diversity on the proteome, here we combine a multiplexed, mass spectrometry-based method for protein quantification with an emerging outbred mouse model containing extensive genetic variation from eight inbred founder strains. By measuring genome-wide transcript and protein expression in livers from 192 Diversity outbred mice, we identify 2,866 protein quantitative trait loci (pQTL) with twice as many local as distant genetic variants. These data support distinct transcriptional and post-transcriptional models underlying the observed pQTL effects. Using a sensitive approach to mediation analysis, we often identified a second protein or transcript as the causal mediator of distant pQTL. Our analysis reveals an extensive network of direct protein–protein interactions. Finally, we show that local genotype can provide accurate predictions of protein abundance in an independent cohort of collaborative cross mice.

Mechanistic basis of infertility of mouse intersubspecific hybrids

Significance Hybrid sterility contributes to speciation by restricting gene flow between related taxa. Although four hybrid sterility genes have been identified in Drosophila and mouse so far, the underlying molecular mechanisms are largely unknown. We describe extensive asynapsis of chromosomes in male and female meiosis of F1 hybrids between two closely related mouse subspecies. Using the intersubspecific chromosome-substitution strains, we demonstrate that the heterospecific pairing of homologous chromosomes is a preexisting condition of asynapsis and may represent a universal mechanism of pachytene arrest in interspecific hybrids. Sex-specific manifestation of asynapsis can explain the mechanism of Haldane’s rule. According to the Dobzhansky–Muller model, hybrid sterility is a consequence of the independent evolution of related taxa resulting in incompatible genomic interactions of their hybrids. The model implies that the incompatibilities evolve randomly, unless a particular gene or nongenic sequence diverges much faster than the rest of the genome. Here we propose that asynapsis of heterospecific chromosomes in meiotic prophase provides a recurrently evolving trigger for the meiotic arrest of interspecific F1 hybrids. We observed extensive asynapsis of chromosomes and disturbance of the sex body in >95% of pachynemas of Mus m. musculus × Mus m. domesticus sterile F1 males. Asynapsis was not preceded by a failure of double-strand break induction, and the rate of meiotic crossing over was not affected in synapsed chromosomes. DNA double-strand break repair was delayed or failed in unsynapsed autosomes, and misexpression of chromosome X and chromosome Y genes was detected in single pachynemas and by genome-wide expression profiling. Oocytes of F1 hybrid females showed the same kind of synaptic problems but with the incidence reduced to half. Most of the oocytes with pachytene asynapsis were eliminated before birth. We propose the heterospecific pairing of homologous chromosomes as a preexisting condition of asynapsis in interspecific hybrids. The asynapsis may represent a universal mechanistic basis of F1 hybrid sterility manifested by pachytene arrest. It is tempting to speculate that a fast-evolving subset of the noncoding genomic sequence important for chromosome pairing and synapsis may be the culprit.

Quantitative Trait Locus Mapping Methods for Diversity Outbred Mice

Genetic mapping studies in the mouse and other model organisms are used to search for genes underlying complex phenotypes. Traditional genetic mapping studies that employ single-generation crosses have poor mapping resolution and limit discovery to loci that are polymorphic between the two parental strains. Multiparent outbreeding populations address these shortcomings by increasing the density of recombination events and introducing allelic variants from multiple founder strains. However, multiparent crosses present new analytical challenges and require specialized software to take full advantage of these benefits. Each animal in an outbreeding population is genetically unique and must be genotyped using a high-density marker set; regression models for mapping must accommodate multiple founder alleles, and complex breeding designs give rise to polygenic covariance among related animals that must be accounted for in mapping analysis. The Diversity Outbred (DO) mice combine the genetic diversity of eight founder strains in a multigenerational breeding design that has been maintained for >16 generations. The large population size and randomized mating ensure the long-term genetic stability of this population. We present a complete analytical pipeline for genetic mapping in DO mice, including algorithms for probabilistic reconstruction of founder haplotypes from genotyping array intensity data, and mapping methods that accommodate multiple founder haplotypes and account for relatedness among animals. Power analysis suggests that studies with as few as 200 DO mice can detect loci with large effects, but loci that account for <5% of trait variance may require a sample size of up to 1000 animals. The methods described here are implemented in the freely available R package DOQTL.

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 F1 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)F1. 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)F1 hybrids, whereas removing one of the two Prdm9 copies in PWD or B6 background had no effect on male reproduction. Incompatibility(ies) not involving Prdm9B6 also acts in the (PWD×B6)F1 hybrids, since the correction of hybrid sterility by Prdm9B6 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)F1 males. Moreover, an increased dosage of Prdm9 and reciprocal cross enhanced fertility of other sperm-carrying male hybrids, (PWD×B6-C3H.Prdm9)F1, 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 F1 hybrids of PWD with B6 and B6-C3H.Prdm9 despite their different prospective fertility, but decreased to 53% after removal of Prdm9B6. Therefore, the Prdm9B6 allele probably takes part in posttranscriptional dominant-negative hybrid interaction(s) absent in the parental strains.

DISSECTING THE GENETIC ARCHITECTURE OF F1 HYBRID STERILITY IN HOUSE MICE

Hybrid sterility as a postzygotic reproductive isolation mechanism has been studied for over 80 years, yet the first identifications of hybrid sterility genes in Drosophila and mouse are quite recent. To study the genetic architecture of F1 hybrid sterility between young subspecies of house mouse Mus m. domesticus and M. m. musculus, we conducted QTL analysis of a backcross between inbred strains representing these two subspecies and probed the role of individual chromosomes in hybrid sterility using the intersubspecific chromosome substitution strains. We provide direct evidence that the asymmetry in male infertility between reciprocal crosses is conferred by the middle region of M. m. musculus Chr X, thus excluding other potential candidates such as Y, imprinted genes, and mitochondrial DNA. QTL analysis identified strong hybrid sterility loci on Chr 17 and Chr X and predicted a set of interchangeable autosomal loci, a subset of which is sufficient to activate the Dobzhansky–Muller incompatibility of the strong loci. Overall, our results indicate the oligogenic nature of F1 hybrid sterility, which should be amenable to reconstruction by proper combination of chromosome substitution strains. Such a prefabricated model system should help to uncover the gene networks and molecular mechanisms underlying hybrid sterility.

X Chromosome Control of Meiotic Chromosome Synapsis in Mouse Inter-Subspecific Hybrids

Hybrid sterility (HS) belongs to reproductive isolation barriers that safeguard the integrity of species in statu nascendi. Although hybrid sterility occurs almost universally among animal and plant species, most of our current knowledge comes from the classical genetic studies on Drosophila interspecific crosses or introgressions. With the house mouse subspecies Mus m. musculus and Mus m. domesticus as a model, new research tools have become available for studies of the molecular mechanisms and genetic networks underlying HS. Here we used QTL analysis and intersubspecific chromosome substitution strains to identify a 4.7 Mb critical region on Chromosome X (Chr X) harboring the Hstx2 HS locus, which causes asymmetrical spermatogenic arrest in reciprocal intersubspecific F1 hybrids. Subsequently, we mapped autosomal loci on Chrs 3, 9 and 13 that can abolish this asymmetry. Combination of immunofluorescent visualization of the proteins of synaptonemal complexes with whole-chromosome DNA FISH on pachytene spreads revealed that heterosubspecific, unlike consubspecific, homologous chromosomes are predisposed to asynapsis in F1 hybrid male and female meiosis. The asynapsis is under the trans- control of Hstx2 and Hst1/Prdm9 hybrid sterility genes in pachynemas of male but not female hybrids. The finding concurred with the fertility of intersubpecific F1 hybrid females homozygous for the Hstx2Mmm allele and resolved the apparent conflict with the dominance theory of Haldanes rule. We propose that meiotic asynapsis in intersubspecific hybrids is a consequence of cis-acting mismatch between homologous chromosomes modulated by the trans-acting Hstx2 and Prdm9 hybrid male sterility genes.

A Multi-Megabase Copy Number Gain Causes Maternal Transmission Ratio Distortion on Mouse Chromosome 2

Significant departures from expected Mendelian inheritance ratios (transmission ratio distortion, TRD) are frequently observed in both experimental crosses and natural populations. TRD on mouse Chromosome (Chr) 2 has been reported in multiple experimental crosses, including the Collaborative Cross (CC). Among the eight CC founder inbred strains, we found that Chr 2 TRD was exclusive to females that were heterozygous for the WSB/EiJ allele within a 9.3 Mb region (Chr 2 76.9 – 86.2 Mb). A copy number gain of a 127 kb-long DNA segment (designated as responder to drive, R2d) emerged as the strongest candidate for the causative allele. We mapped R2d sequences to two loci within the candidate interval. R2d1 is located near the proximal boundary, and contains a single copy of R2d in all strains tested. R2d2 maps to a 900 kb interval, and the number of R2d copies varies from zero in classical strains (including the mouse reference genome) to more than 30 in wild-derived strains. Using real-time PCR assays for the copy number, we identified a mutation (R2d2WSBdel1) that eliminates the majority of the R2d2WSB copies without apparent alterations of the surrounding WSB/EiJ haplotype. In a three-generation pedigree segregating for R2d2WSBdel1, the mutation is transmitted to the progeny and Mendelian segregation is restored in females heterozygous for R2d2WSBdel1, thus providing direct evidence that the copy number gain is causal for maternal TRD. We found that transmission ratios in R2d2WSB heterozygous females vary between Mendelian segregation and complete distortion depending on the genetic background, and that TRD is under genetic control of unlinked distorter loci. Although the R2d2WSB transmission ratio was inversely correlated with average litter size, several independent lines of evidence support the contention that female meiotic drive is the cause of the distortion. We discuss the implications and potential applications of this novel meiotic drive system.

Genetic Architectures of Quantitative Variation in RNA Editing Pathways

RNA editing refers to post-transcriptional processes that alter the base sequence of RNA. Recently, hundreds of new RNA editing targets have been reported. However, the mechanisms that determine the specificity and degree of editing are not well understood. We examined quantitative variation of site-specific editing in a genetically diverse multiparent population, Diversity Outbred mice, and mapped polymorphic loci that alter editing ratios globally for C-to-U editing and at specific sites for A-to-I editing. An allelic series in the C-to-U editing enzyme Apobec1 influences the editing efficiency of Apob and 58 additional C-to-U editing targets. We identified 49 A-to-I editing sites with polymorphisms in the edited transcript that alter editing efficiency. In contrast to the shared genetic control of C-to-U editing, most of the variable A-to-I editing sites were determined by local nucleotide polymorphisms in proximity to the editing site in the RNA secondary structure. Our results indicate that RNA editing is a quantitative trait subject to genetic variation and that evolutionary constraints have given rise to distinct genetic architectures in the two canonical types of RNA editing.

High-Density Genotypes of Inbred Mouse Strains: Improved Power and Precision of Association Mapping

Human genome-wide association studies have identified thousands of loci associated with disease phenotypes. Genome-wide association studies also have become feasible using rodent models and these have some important advantages over human studies, including controlled environment, access to tissues for molecular profiling, reproducible genotypes, and a wide array of techniques for experimental validation. Association mapping with common mouse inbred strains generally requires 100 or more strains to achieve sufficient power and mapping resolution; in contrast, sample sizes for human studies typically are one or more orders of magnitude greater than this. To enable well-powered studies in mice, we have generated high-density genotypes for ∼175 inbred strains of mice using the Mouse Diversity Array. These new data increase marker density by 1.9-fold, have reduced missing data rates, and provide more accurate identification of heterozygous regions compared with previous genotype data. We report the discovery of new loci from previously reported association mapping studies using the new genotype data. The data are freely available for download, and Web-based tools provide easy access for association mapping and viewing of the underlying intensity data for individual loci.

Genetic Analysis of Substrain Divergence in Non-Obese Diabetic (NOD) Mice

The non-obese diabetic (NOD) mouse is a polygenic model for type 1 diabetes that is characterized by insulitis, a leukocytic infiltration of the pancreatic islets. During ~35 years since the original inbred strain was developed in Japan, NOD substrains have been established at different laboratories around the world. Although environmental differences among NOD colonies capable of impacting diabetes incidence have been recognized, differences arising from genetic divergence have not been analyzed previously. We use both mouse diversity array and whole-exome capture sequencing platforms to identify genetic differences distinguishing five NOD substrains. We describe 64 single-nucleotide polymorphisms, and two short indels that differ in coding regions of the five NOD substrains. A 100-kb deletion on Chromosome 3 distinguishes NOD/ShiLtJ and NOD/ShiLtDvs from three other substrains, whereas a 111-kb deletion in the Icam2 gene on Chromosome 11 is unique to the NOD/ShiLtDvs genome. The extent of genetic divergence for NOD substrains is compared with similar studies for C57BL6 and BALB/c substrains. As mutations are fixed to homozygosity by continued inbreeding, significant differences in substrain phenotypes are to be expected. These results emphasize the importance of using embryo freezing methods to minimize genetic drift within substrains and of applying appropriate genetic nomenclature to permit substrain recognition when one is used.