Kenneth Paigen

Mouse population-guided resequencing reveals that variants in CD44 contribute to acetaminophen-induced liver injury in humans

Interindividual variability in response to chemicals and drugs is a common regulatory concern. It is assumed that xenobiotic-induced adverse reactions have a strong genetic basis, but many mechanism-based investigations have not been successful in identifying susceptible individuals. While recent advances in pharmacogenetics of adverse drug reactions show promise, the small size of the populations susceptible to important adverse events limits the utility of whole-genome association studies conducted entirely in humans. We present a strategy to identify genetic polymorphisms that may underlie susceptibility to adverse drug reactions. First, in a cohort of healthy adults who received the maximum recommended dose of acetaminophen (4 g/d x 7 d), we confirm that about one third of subjects develop elevations in serum alanine aminotransferase, indicative of liver injury. To identify the genetic basis for this susceptibility, a panel of 36 inbred mouse strains was used to model genetic diversity. Mice were treated with 300 mg/kg or a range of additional acetaminophen doses, and the extent of liver injury was quantified. We then employed whole-genome association analysis and targeted sequencing to determine that polymorphisms in Ly86, Cd44, Cd59a, and Capn8 correlate strongly with liver injury and demonstrated that dose-curves vary with background. Finally, we demonstrated that variation in the orthologous human gene, CD44, is associated with susceptibility to acetaminophen in two independent cohorts. Our results indicate a role for CD44 in modulation of susceptibility to acetaminophen hepatotoxicity. These studies demonstrate that a diverse mouse population can be used to understand and predict adverse toxicity in heterogeneous human populations through guided resequencing.

The recombinational anatomy of a mouse chromosome.

Among mammals, genetic recombination occurs at highly delimited sites known as recombination hotspots. They are typically 1–2 kb long and vary as much as a 1,000-fold or more in recombination activity. Although much is known about the molecular details of the recombination process itself, the factors determining the location and relative activity of hotspots are poorly understood. To further our understanding, we have collected and mapped the locations of 5,472 crossover events along mouse Chromosome 1 arising in 6,028 meioses of male and female reciprocal F1 hybrids of C57BL/6J and CAST/EiJ mice. Crossovers were mapped to a minimum resolution of 225 kb, and those in the telomere-proximal 24.7 Mb were further mapped to resolve individual hotspots. Recombination rates were evolutionarily conserved on a regional scale, but not at the local level. There was a clear negative-exponential relationship between the relative activity and abundance of hotspot activity classes, such that a small number of the most active hotspots account for the majority of recombination. Females had 1.2× higher overall recombination than males did, although the sex ratio showed considerable regional variation. Locally, entirely sex-specific hotspots were rare. The initiation of recombination at the most active hotspot was regulated independently on the two parental chromatids, and analysis of reciprocal crosses indicated that parental imprinting has subtle effects on recombination rates. It appears that the regulation of mammalian recombination is a complex, dynamic process involving multiple factors reflecting species, sex, individual variation within species, and the properties of individual hotspots.

Health and population effects of rare gene knockouts in adult humans with related parents

Rare gene knockouts in adult humans On average, most peoples genomes contain approximately 100 completely nonfunctional genes. These loss-of-function (LOF) mutations tend to be rare and/or occur only as a single copy within individuals. Narasimhan et al. investigated LOF in a Pakistani population with high levels of consanguinity. Examining LOF alleles that were identical by descent, they found, as expected, an absence of homozygote LOF for certain protein-coding genes. However, they also identified many homozygote LOF alleles with no apparent deleterious phenotype, including some that were expected to confer genetic disease. Indeed, one family had lost the recombination-associated gene PRDM9. Science, this issue p. 474 The total loss of protein-coding genes, even those with the potential to confer genetic diseases, can be tolerated. Examining complete gene knockouts within a viable organism can inform on gene function. We sequenced the exomes of 3222 British adults of Pakistani heritage with high parental relatedness, discovering 1111 rare-variant homozygous genotypes with predicted loss of function (knockouts) in 781 genes. We observed 13.7% fewer homozygous knockout genotypes than we expected, implying an average load of 1.6 recessive-lethal-equivalent loss-of-function (LOF) variants per adult. When genetic data were linked to the individuals’ lifelong health records, we observed no significant relationship between gene knockouts and clinical consultation or prescription rate. In this data set, we identified a healthy PRDM9-knockout mother and performed phased genome sequencing on her, her child, and control individuals. Our results show that meiotic recombination sites are localized away from PRDM9-dependent hotspots. Thus, natural LOF variants inform on essential genetic loci and demonstrate PRDM9 redundancy in humans.

Evidence of a Large-Scale Functional Organization of Mammalian Chromosomes

Evidence from inbred strains of mice indicates that a quarter or more of the mammalian genome consists of chromosome regions containing clusters of functionally related genes. The intense selection pressures during inbreeding favor the coinheritance of optimal sets of alleles among these genetically linked, functionally related genes, resulting in extensive domains of linkage disequilibrium (LD) among a set of 60 genetically diverse inbred strains. Recombination that disrupts the preferred combinations of alleles reduces the ability of offspring to survive further inbreeding. LD is also seen between markers on separate chromosomes, forming networks with scale-free architecture. Combining LD data with pathway and genome annotation databases, we have been able to identify the biological functions underlying several domains and networks. Given the strong conservation of gene order among mammals, the domains and networks we find in mice probably characterize all mammals, including humans.

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.

Affinity-seq detects genome-wide PRDM9 binding sites and reveals the impact of prior chromatin modifications on mammalian recombination hotspot usage

BackgroundGenetic recombination plays an important role in evolution, facilitating the creation of new, favorable combinations of alleles and the removal of deleterious mutations by unlinking them from surrounding sequences. In most mammals, the placement of genetic crossovers is determined by the binding of PRDM9, a highly polymorphic protein with a long zinc finger array, to its cognate binding sites. It is one of over 800 genes encoding proteins with zinc finger domains in the human genome.ResultsWe report a novel technique, Affinity-seq, that for the first time identifies both the genome-wide binding sites of DNA-binding proteins and quantitates their relative affinities. We have applied this in vitro technique to PRDM9, the zinc-finger protein that activates genetic recombination, obtaining new information on the regulation of hotspots, whose locations and activities determine the recombination landscape. We identified 31,770 binding sites in the mouse genome for the PRDM9Dom2 variant. Comparing these results with hotspot usage in vivo, we find that less than half of potential PRDM9 binding sites are utilized in vivo. We show that hotspot usage is increased in actively transcribed genes and decreased in genomic regions containing H3K9me2/3 histone marks or bound to the nuclear lamina.ConclusionsThese results show that a major factor determining whether a binding site will become an active hotspot and what its activity will be are constraints imposed by prior chromatin modifications on the ability of PRDM9 to bind to DNA in vivo. These constraints lead to the presence of long genomic regions depleted of recombination.

Nuclear localization of PRDM9 and its role in meiotic chromatin modifications and homologous synapsis

Developmental progress of germ cells through meiotic phases is closely tied to ongoing meiotic recombination. In mammals, recombination preferentially occurs in genomic regions known as hotspots; the protein that activates these hotspots is PRDM9, containing a genetically variable zinc finger (ZNF) domain and a PR-SET domain with histone H3K4 trimethyltransferase activity. PRDM9 is required for fertility in mice, but little is known about its localization and developmental dynamics. Application of spermatogenic stage-specific markers demonstrates that PRDM9 accumulates in male germ cell nuclei at pre-leptonema to early leptonema but is no longer detectable in nuclei by late zygonema. By the pachytene stage, PRDM9-dependent histone H3K4 trimethyl marks on hotspots also disappear. PRDM9 localizes to nuclei concurrently with the deposition of meiotic cohesin complexes, but is not required for incorporation of cohesin complex proteins into chromosomal axial elements, or accumulation of normal numbers of RAD51 foci on meiotic chromatin by late zygonema. Germ cells lacking PRDM9 exhibit inefficient homology recognition and synapsis, with aberrant repair of meiotic DNA double-strand breaks and transcriptional abnormalities characteristic of meiotic silencing of unsynapsed chromatin. Together, these results on the developmental time course for nuclear localization of PRDM9 establish its direct window of function and demonstrate the independence of chromosome axial element formation from the concurrent PRDM9-mediated activation of recombination hotspots.

PRDM9 interactions with other proteins provide a link between recombination hotspots and the chromosomal axis in meiosis

Meiotic recombination hotspots activated by PRDM9 are associated with the chromosomal axis and synaptonemal complex via their interaction with other proteins, including CDYL, EHMT2, EWSR1, and CXXC1.

Clusters of Ancestrally Related Genes That Show Paralogy in Whole or in Part Are a Major Feature of the Genomes of Humans and Other Species

Arrangements of genes along chromosomes are a product of evolutionary processes, and we can expect that preferable arrangements will prevail over the span of evolutionary time, often being reflected in the non-random clustering of structurally and/or functionally related genes. Such non-random arrangements can arise by two distinct evolutionary processes: duplications of DNA sequences that give rise to clusters of genes sharing both sequence similarity and common sequence features and the migration together of genes related by function, but not by common descent [1], [2], [3]. To provide a background for distinguishing between the two, which is important for future efforts to unravel the evolutionary processes involved, we here provide a description of the extent to which ancestrally related genes are found in proximity. Towards this purpose, we combined information from five genomic datasets, InterPro, SCOP, PANTHER, Ensembl protein families, and Ensembl gene paralogs. The results are provided in publicly available datasets (http://cgd.jax.org/datasets/clustering/paraclustering.shtml) describing the extent to which ancestrally related genes are in proximity beyond what is expected by chance (i.e. form paraclusters) in the human and nine other vertebrate genomes, as well as the D. melanogaster, C. elegans, A. thaliana, and S. cerevisiae genomes. With the exception of Saccharomyces, paraclusters are a common feature of the genomes we examined. In the human genome they are estimated to include at least 22% of all protein coding genes. Paraclusters are far more prevalent among some gene families than others, are highly species or clade specific and can evolve rapidly, sometimes in response to environmental cues. Altogether, they account for a large portion of the functional clustering previously reported in several genomes.

Colocalization of somatic and meiotic double strand breaks near the Myc oncogene on mouse chromosome 15

Both somatic and meiotic recombinations involve the repair of DNA double strand breaks (DSBs) that occur at preferred locations in the genome. Improper repair of DSBs during either mitosis or meiosis can lead to mutations, chromosomal aberration such as translocations, cancer, and/or cell death. Currently, no model exists that explains the locations of either spontaneous somatic DSBs or programmed meiotic DSBs or relates them to each other. One common class of tumorigenic translocations arising from DSBs is chromosomal rearrangements near the Myc oncogene. Myc translocations have been associated with Burkitt lymphoma in humans, plasmacytoma in mice, and immunocytoma in rats. Comparing the locations of somatic and meiotic DSBs near the mouse Myc oncogene, we demonstrated that the placement of these DSBs is not random and that both events clustered in the same short discrete region of the genome. Our work shows that both somatic and meiotic DSBs tend to occur in proximity to each other within the Myc region, suggesting that they share common originating features. It is likely that some regions of the genome are more susceptible to both somatic and meiotic DSBs, and the locations of meiotic hotspots may be an indicator of genomic regions more susceptible to DNA damage.