Jackie E. Lim

Latent TGF-β–binding protein 4 modifies muscular dystrophy in mice

Most single-gene diseases, including muscular dystrophy, display a nonuniform phenotype. Phenotypic variability arises, in part, due to the presence of genetic modifiers that enhance or suppress the disease process. We employed an unbiased mapping approach to search for genes that modify muscular dystrophy in mice. In a genome-wide scan, we identified a single strong locus on chromosome 7 that influenced two pathological features of muscular dystrophy, muscle membrane permeability and muscle fibrosis. Within this genomic interval, an insertion/deletion polymorphism of 36 bp in the coding region of the latent TGF-beta-binding protein 4 gene (Ltbp4) was found. Ltbp4 encodes a latent TGF-beta-binding protein that sequesters TGF-beta and regulates its availability for binding to the TGF-beta receptor. Insertion of 12 amino acids into the proline-rich region of LTBP4 reduced proteolytic cleavage and was associated with reduced TGF-beta signaling, decreased fibrosis, and improved muscle pathology in a mouse model of muscular dystrophy. In contrast, a 12-amino-acid deletion in LTBP4 was associated with increased proteolysis, SMAD signaling, and fibrosis. These data identify Ltbp4 as a target gene to regulate TGF-beta signaling and modify outcomes in muscular dystrophy.

A Common and Unstable Copy Number Variant Is Associated with Differences in Glo1 Expression and Anxiety-Like Behavior

Glyoxalase 1 (Glo1) has been implicated in anxiety-like behavior in mice and in multiple psychiatric diseases in humans. We used mouse Affymetrix exon arrays to detect copy number variants (CNV) among inbred mouse strains and thereby identified a ∼475 kb tandem duplication on chromosome 17 that includes Glo1 (30,174,390–30,651,226 Mb; mouse genome build 36). We developed a PCR-based strategy and used it to detect this duplication in 23 of 71 inbred strains tested, and in various outbred and wild-caught mice. Presence of the duplication is associated with a cis-acting expression QTL for Glo1 (LOD>30) in BXD recombinant inbred strains. However, evidence for an eQTL for Glo1 was not obtained when we analyzed single SNPs or 3-SNP haplotypes in a panel of 27 inbred strains. We conclude that association analysis in the inbred strain panel failed to detect an eQTL because the duplication was present on multiple highly divergent haplotypes. Furthermore, we suggest that non-allelic homologous recombination has led to multiple reversions to the non-duplicated state among inbred strains. We show associations between multiple duplication-containing haplotypes, Glo1 expression and anxiety-like behavior in both inbred strain panels and outbred CD-1 mice. Our findings provide a molecular basis for differential expression of Glo1 and further implicate Glo1 in anxiety-like behavior. More broadly, these results identify problems with commonly employed tests for association in inbred strains when CNVs are present. Finally, these data provide an example of biologically significant phenotypic variability in model organisms that can be attributed to CNVs.

A mutation in Sec15l1 causes anemia in hemoglobin deficit ( hbd ) mice

Hemoglobin deficit (hbd) mice carry a spontaneous mutation that impairs erythroid iron assimilation but does not cause other defects. Normal delivery of iron to developing erythroid precursors is highly dependent on the transferrin cycle. Through genetic mapping and complementation experiments, we show that the hbd mutation is an in-frame deletion of a conserved exon of the mouse gene Sec15l1, encoding one of two Sec15 proteins implicated in the mammalian exocyst complex. Sec15l1 is linked to the transferrin cycle through its interaction with Rab11, a GTPase involved in vesicular trafficking. We propose that inactivation of Sec15l1 alters recycling of transferrin cycle endosomes and increases the release of transferrin receptor exocytic vesicles. This in turn decreases erythroid iron uptake. Determining the molecular basis of the hbd phenotype provides new insight into the intricate mechanisms necessary for normal erythroid iron uptake and the function of a mammalian exocyst protein.

Genome-Wide Association Studies and the Problem of Relatedness Among Advanced Intercross Lines and Other Highly Recombinant Populations

Model organisms offer many advantages for the genetic analysis of complex traits. However, identification of specific genes is often hampered by a lack of recombination between the genomes of inbred progenitors. Recently, genome-wide association studies (GWAS) in humans have offered gene-level mapping resolution that is possible because of the large number of accumulated recombinations among unrelated human subjects. To obtain analogous improvements in mapping resolution in mice, we used a 34th generation advanced intercross line (AIL) derived from two inbred strains (SM/J and LG/J). We used simulations to show that familial relationships among subjects must be accounted for when analyzing these data; we then used a mixed model that included polygenic effects to address this problem in our own analysis. Using a combination of F2 and AIL mice derived from the same inbred progenitors, we identified genome-wide significant, subcentimorgan loci that were associated with methamphetamine sensitivity, (e.g., chromosome 18; LOD = 10.5) and non-drug-induced locomotor activity (e.g., chromosome 8; LOD = 18.9). The 2-LOD support interval for the former locus contains no known genes while the latter contains only one gene (Csmd1). This approach is broadly applicable in terms of phenotypes and model organisms and allows GWAS to be performed in multigenerational crosses between and among inbred strains where familial relatedness is often unavoidable.

Fine mapping of QTL for prepulse inhibition in LG/J and SM/J mice using F2 and advanced intercross lines

Prepulse inhibition (PPI) of the startle response is a measure of sensorimotor gating, a process that filters out extraneous sensory, motor and cognitive information. Humans with neurological and psychiatric disorders, including schizophrenia, obsessive‐compulsive disorder and Huntingtons disease, exhibit a reduction in PPI. Habituation of the startle response is also disrupted in schizophrenic patients. In order to elucidate the genes involved in sensorimotor gating, we phenotyped 472 mice from an F2 cross between LG/J × SM/J for PPI and genotyped these mice genome‐wide using 162 single nucleotide polymorphism (SNP) markers. We used prepulse intensity levels that were 3, 6 and 12 dB above background (PPI3, PPI6 and PPI12, respectively). We identified a significant quantitative trait locus (QTL) on chromosome 12 for all three prepulse intensities as well as a significant QTL for both PPI6 and PPI12 on chromosome 11. We identified QTLs on chromosomes 7 and 17 for the startle response when sex was included as an interactive covariate and found a QTL for habituation of the startle response on chromosome 4. We also phenotyped 135 mice from an F34 advanced intercross line (AIL) between LG/J × SM/J for PPI and genotyped them at more than 3000 SNP markers. Inclusions of data from the AIL mice reduced the size of several of these QTLs to less than 5 cM. These results will be useful for identifying genes that influence sensorimotor gaiting and show the power of AIL for fine mapping of QTLs.

Fine-mapping of muscle weight QTL in LG/J and SM/J intercrosses.

Genetic variation plays a substantial role in variation in strength, but the underlying mechanisms remain poorly understood. The objective of the present study was to examine the mechanisms underlying variation in muscle mass, a predictor of strength, between LG/J and SM/J strains, which are the inbred progeny of mice selected, respectively, for high and low body weight. We measured weight of five hindlimb muscles in LG/J and SM/J males and females, in F(1) and F(2) intercrosses, and in an advanced intercross (AI), F(34), between the two. F(2) mice were genotyped using 162 SNPs throughout the genome; F(34) mice were genotyped at 3,015 SNPs. A twofold difference in muscle mass between the LG/J and SM/J mouse strains was observed. Integrated genome-wide association analysis in the combined population of F(2) and AI identified 22 quantitative trait loci (QTL; genome-wide P < 0.05) affecting muscle weight on Chr 2 (2 QTL), 4, 5, 6 (7 QTL), 7 (4 QTL), 8 (4 QTL), and 11 (3 QTL). The LG/J allele conferred greater muscle weight in all cases. The 1.5-LOD QTL support intervals ranged between 0.3 and 13.4 Mb (median 3.7 Mb) restricting the list of candidates to between 5 and 97 genes. Selection for body weight segregated the alleles affecting skeletal muscle, the most abundant tissue in the body. Combination of analyses in an F(2) and AI was an effective strategy to detect and refine the QTL in a genome-wide manner. The achieved resolution facilitates further elucidation of the underlying genetic mechanisms affecting muscle mass.

Fine-mapping alleles for body weight in LG/J × SM/J F2 and F34 advanced intercross lines

The present study measured variation in body weight using a combined analysis in an F2 intercross and an F34 advanced intercross line (AIL). Both crosses were derived from inbred LG/J and SM/J mice, which were selected for large and small body size prior to inbreeding. Body weight was measured at 62 (±5) days of age. Using an integrated GWAS and forward model selection approach, we identified 11 significant QTLs that affected body weight on ten different chromosomes. With these results we developed a full model that explained over 18% of the phenotypic variance. The median 1.5-LOD support interval was 5.55 Mb, which is a significant improvement over most prior body weight QTLs. We identified nonsynonymous coding SNPs between LG/J and SM/J mice in order to further narrow the list of candidate genes. Three of the genes with nonsynonymous coding SNPs (Rad23b, Stk33, and Anks1b) have been associated with adiposity, waist circumference, and body mass index in human GWAS, thus providing evidence that these genes may underlie our QTLs. Our results demonstrate that a relatively small number of loci contribute significantly to the phenotypic variance in body weight, which is in marked contrast to the situation in humans. This difference is likely to be the result of strong selective pressure and the simplified genetic architecture, both of which are important advantages of our system.

Anxiety and fear in a cross of C57BL/6J and DBA/2J mice: mapping overlapping and independent QTL for related traits.

Anxiety, like other psychiatric disorders, is a complex neurobehavioral trait, making identification of causal genes difficult. In this study, we examined anxiety‐like behavior and fear conditioning (FC) in an F2 intercross of C57BL/6J and DBA/2J mice. We identified numerous quantitative trait loci (QTL) influencing anxiety‐like behavior in both open field (OF) and FC tests. Many of these QTL were mapped back to the same chromosomal regions, regardless of behavior or test. For example, highly significant overlapping QTL on chromosome 1 were found in all FC measures as well as in center time measures in the OF. Other QTL exhibited strong temporal profiles over testing, highlighting dynamic relationship between genotype, test and changes in behavior. Next, we implemented a factor analysis design to account for the correlated nature of the behaviors measured. OF and FC behaviors loaded onto four main factors representing both anxiety and fear behaviors. Using multiple QTL modeling, we calculated the percentage variance in anxiety and fear explained by multiple QTL using both additive and interactive terms. Quantitative trait loci modeling resulted in a broad description of the genetic architecture underlying anxiety and fear accounting for 14–37% of trait variance. Factor analysis and multiple QTL modeling showed both unique and shared QTL for anxiety and fear; suggesting a partially overlapping genetic architecture for these two different models of anxiety.

Hnrnph1 Is A Quantitative Trait Gene for Methamphetamine Sensitivity

Psychostimulant addiction is a heritable substance use disorder; however its genetic basis is almost entirely unknown. Quantitative trait locus (QTL) mapping in mice offers a complementary approach to human genome-wide association studies and can facilitate environment control, statistical power, novel gene discovery, and neurobiological mechanisms. We used interval-specific congenic mouse lines carrying various segments of chromosome 11 from the DBA/2J strain on an isogenic C57BL/6J background to positionally clone a 206 kb QTL (50,185,512–50,391,845 bp) that was causally associated with a reduction in the locomotor stimulant response to methamphetamine (2 mg/kg, i.p.; DBA/2J < C57BL/6J)—a non-contingent, drug-induced behavior that is associated with stimulation of the dopaminergic reward circuitry. This chromosomal region contained only two protein coding genes—heterogeneous nuclear ribonucleoprotein, H1 (Hnrnph1) and RUN and FYVE domain-containing 1 (Rufy1). Transcriptome analysis via mRNA sequencing in the striatum implicated a neurobiological mechanism involving a reduction in mesolimbic innervation and striatal neurotransmission. For instance, Nr4a2 (nuclear receptor subfamily 4, group A, member 2), a transcription factor crucial for midbrain dopaminergic neuron development, exhibited a 2.1-fold decrease in expression (DBA/2J < C57BL/6J; p 4.2 x 10−15). Transcription activator-like effector nucleases (TALENs)-mediated introduction of frameshift deletions in the first coding exon of Hnrnph1, but not Rufy1, recapitulated the reduced methamphetamine behavioral response, thus identifying Hnrnph1 as a quantitative trait gene for methamphetamine sensitivity. These results define a novel contribution of Hnrnph1 to neurobehavioral dysfunction associated with dopaminergic neurotransmission. These findings could have implications for understanding the genetic basis of methamphetamine addiction in humans and the development of novel therapeutics for prevention and treatment of substance abuse and possibly other psychiatric disorders.