Integration of QTL and bioinformatic tools to identify candidate genes for triglycerides in mice

To identify genetic loci influencing lipid levels, we performed quantitative trait loci (QTL) analysis between inbred mouse strains MRL/MpJ and SM/J, measuring triglyceride levels at 8 weeks of age in F2 mice fed a chow diet. We identified one significant QTL on chromosome (Chr) 15 and three suggestive QTL on Chrs 2, 7, and 17. We also carried out microarray analysis on the livers of parental strains of 282 F2 mice and used these data to find cis-regulated expression QTL. We then narrowed the list of candidate genes under significant QTL using a "toolbox" of bioinformatic resources, including haplotype analysis; parental strain comparison for gene expression differences and nonsynonymous coding single nucleotide polymorphisms (SNP); cis-regulated eQTL in livers of F2 mice; correlation between gene expression and phenotype; and conditioning of expression on the phenotype. We suggest Slc25a7 as a candidate gene for the Chr 7 QTL and, based on expression differences, five genes (Polr3 h, Cyp2d22, Cyp2d26, Tspo, and Ttll12) as candidate genes for Chr 15 QTL. This study shows how bioinformatics can be used effectively to reduce candidate gene lists for QTL related to complex traits.     

A novel candidate gene for mouse and human preaxial polydactyly with altered expression in limbs of Hemimelic extra-toes mutant mice

Polydactyly is a common malformation of vertebrate limbs. In humans a major locus for nonsyndromic pre-axial polydactyly (PPD) has been mapped previously to 7q36. The mouse Hemimelic extra-toes (Hx) mutation maps to a homologous chromosome segment and has been proposed to affect a homologous gene. To understand the molecular changes underlying PPD, we used a positional cloning approach to identify the gene or genes disrupted by the Hx mutation and a closely linked limb mutation, Hammertoe (Hm). High resolution genetic mapping identified a small candidate interval for the mouse mutations located 1.2 cM distal to the Shh locus. The nonrecombinant interval was completely cloned in bacterial artificial chromosomes and searched for genes using a combination of exon trapping, sample sequencing, and mapping of known genes. Two novel genes, Lmbr1 and Lmbr2, are entirely within the candidate interval we defined genetically. The open reading frame of both genes is intact in mutant mice, but the expression of the Lmbr1 gene is dramatically altered in developing limbs of Hx mutant mice. The correspondence between the spatial and temporal changes in Lmbr1 expression and the embryonic onset of the Hx mutant phenotype suggests that the mouse Hx mutation may be a regulatory allele of Lmbr1. The human ortholog of Lmbr1 maps within the recently described interval for human PPD, strengthening the possibility that both mouse and human limb abnormalities are due to defects in the same highly conserved gene.     

Mapping of a Chromosome 12 Region Associated with Airway Hyperresponsiveness in a Recombinant Congenic Mouse Strain and Selection of Potential Candidate Genes by Expression and Sequence Variation Analyses

In a previous study we determined that BcA86 mice, a strain belonging to a panel of AcB/BcA recombinant congenic strains, have an airway responsiveness phenotype resembling mice from the airway hyperresponsive A/J strain. The majority of the BcA86 genome is however from the hyporesponsive C57BL/6J strain. The aim of this study was to identify candidate regions and genes associated with airway hyperresponsiveness (AHR) by quantitative trait locus (QTL) analysis using the BcA86 strain. Airway responsiveness of 205 F2 mice generated from backcrossing BcA86 strain to C57BL/6J strain was measured and used for QTL analysis to identify genomic regions in linkage with AHR. Consomic mice for the QTL containing chromosomes were phenotyped to study the contribution of each chromosome to lung responsiveness. Candidate genes within the QTL were selected based on expression differences in mRNA from whole lungs, and the presence of coding non-synonymous mutations that were predicted to have a functional effect by amino acid substitution prediction tools. One QTL for AHR was identified on Chromosome 12 with its 95% confidence interval ranging from 54.6 to 82.6 Mbp and a maximum LOD score of 5.11 (p = 3.68×10⁻³). We confirmed that the genotype of mouse Chromosome 12 is an important determinant of lung responsiveness using a Chromosome 12 substitution strain. Mice with an A/J Chromosome 12 on a C57BL/6J background have an AHR phenotype similar to hyperresponsive strains A/J and BcA86. Within the QTL, genes with deleterious coding variants, such as Foxa1, and genes with expression differences, such as Mettl21d and Snapc1, were selected as possible candidates for the AHR phenotype. Overall, through QTL analysis of a recombinant congenic strain, microarray analysis and coding variant analysis we identified Chromosome 12 and three potential candidate genes to be in linkage with airway responsiveness.     

In a shake of a lamb's tail: using genomics to unravel a cause of chondrodysplasia in Texel sheep

Chondrodysplasia in Texel sheep is a recessively inherited disorder characterized by dwarfism and angular deformities of the forelimbs. A genome-wide association study using the Illumina OvineSNP50 BeadChip on 15 sheep diagnosed as affected and eight carriers descended from three affected rams was conducted to uncover the genetic cause. A homozygous region of 25 consecutive single nucleotide polymorphism (SNP) loci was identified in all affected sheep, covering a region of 1 Mbp on ovine chromosome 4. Seven positional candidate genes - including the solute carrier family 13 (sodium/sulphate symporters), member 1 (SLC13A1) - were identified and used to search for new SNPs for fine mapping of the causal locus. The SLC13A1 gene, encoding a sodium/sulphate transporter, was the primary candidate gene attributable to similar phenotypes observed in the Slc13a1 knockout mouse model. We discovered a 1-bp deletion of T (g.25513delT) at the 107 bp position of exon 3 in the SLC13A1 gene. Genotyping by direct sequencing and restriction fragment length polymorphism analysis for this mutation showed that all 15 affected sheep were g.25513delT/g.25513delT; the eight carriers were g.25513delT/T and 54 normal controls were T/T. The mutation g.25513delT shifts the open reading frame of SLC13A1 to introduce a stop codon and truncate C-terminal amino acids. It was concluded that the g.25513delT mutation in the SLC13A1 gene was responsible for the chondrodysplasia seen in these Texel sheep. This knowledge can be used to identify carriers with the defective g.[25513delT] allele to avoid at-risk matings to improve animal welfare and decrease economic losses.     

Characterization of Bglu3, a mouse fasting glucose locus, and identification of Apcs as an underlying candidate gene

Bglu3 is a quantitative trait locus for fasting glucose on distal chromosome 1 identified in an intercross between C57BL/6 (B6) and C3H/HeJ (C3H) apolipoprotein E-deficient (apoE(-/-)) mice. This locus was subsequently replicated in two separate mouse intercrosses. The objective of this study was to characterize Bglu3 through construction and analysis of a congenic strain and identify underlying candidate genes. Congenic mice were constructed by introgressing a genomic region harboring Bglu3 from C3H.apoE(-/-) into B6.apoE(-/-) mice. Mice were started with a Western diet at 6 wk of age and maintained on the diet for 12 wk. Gene expression in the liver was analyzed by microarrays. Congenic mice had significantly higher fasting glucose levels and developed more significant glucose intolerance compared with B6.apoE(-/-) mice on the Western diet. Microarray analysis revealed 336 genes to be differentially expressed in the liver of congenic mice. Further pathway analysis suggested a role for acute phase response signaling in regulating glucose intolerance. Apcs, encoding an acute phase response protein serum amyloid P (SAP), is located underneath the linkage peak of Bglu3. Multiple single nucleotide polymorphisms between B6 and C3H mice were detected within and surrounding Apcs. Apcs expression in the liver was significantly higher in congenic and C3H mice compared with B6 mice. The Western diet consumption led to a gradual rise in plasma SAP levels, which was accompanied by rising fasting glucose in both B6 and C3H apoE(-/-) mice. Expression of C3H Apcs in B6.apoE(-/-) mice aggravated glucose intolerance. Bglu3 is confirmed to be a locus affecting diabetes susceptibility, and Apcs is a probable candidate gene.     

A Spontaneous Fatp4/Scl27a4 Splice Site Mutation in a New Murine Model for Congenital Ichthyosis

Congenital ichthyoses are life-threatening conditions in humans. We describe here the identification and molecular characterization of a novel recessive mutation in mice that results in newborn lethality with severe congenital lamellar ichthyosis. Mutant newborns have a taut, shiny, non-expandable epidermis that resembles cornified manifestations of autosomal-recessive congenital ichthyosis in humans. The skin is stretched so tightly that the newborn mice are immobilized. The genetic defect was mapped to a region near the proximal end of chromosome 2 by SNP analysis, suggesting Fatp4/Slc27a4 as a candidate gene. FATP4 mutations in humans cause ichthyosis prematurity syndrome (IPS), and mutations of Fatp4 in mice have previously been found to cause a phenotype that resembles human congenital ichthyoses. Characterization of the Fatp4 cDNA revealed a fusion of exon 8 to exon 10, with deletion of exon 9. Genomic sequencing identified an A to T mutation in the splice donor sequence at the 3′-end of exon 9. Loss of exon 9 results in a frame shift mutation upstream from the conserved very long-chain acyl-CoA synthase (VLACS) domain. Histological studies revealed that the mutant mice have defects in keratinocyte differentiation, along with hyperproliferation of the stratum basale of the epidermis, a hyperkeratotic stratum corneum, and reduced numbers of secondary hair follicles. Since Fatp4 protein is present primarily at the stratum granulosum and the stratum spinosum, the hyperproliferation and the alterations in hair follicle induction suggest that very long chain fatty acids, in addition to being required for normal cornification, may influence signals from the stratum corneum to the basal cells that help to orchestrate normal skin differentiation.     

Pla2g12b and Hpn Are Genes Identified by Mouse ENU Mutagenesis That Affect HDL Cholesterol

Despite considerable progress understanding genes that affect the HDL particle, its function, and cholesterol content, genes identified to date explain only a small percentage of the genetic variation. We used N-ethyl-N-nitrosourea mutagenesis in mice to discover novel genes that affect HDL cholesterol levels. Two mutant lines (Hlb218 and Hlb320) with low HDL cholesterol levels were established. Causal mutations in these lines were mapped using linkage analysis: for line Hlb218 within a 12 Mbp region on Chr 10; and for line Hlb320 within a 21 Mbp region on Chr 7. High-throughput sequencing of Hlb218 liver RNA identified a mutation in Pla2g12b. The transition of G to A leads to a cysteine to tyrosine change and most likely causes a loss of a disulfide bridge. Microarray analysis of Hlb320 liver RNA showed a 7-fold downregulation of Hpn; sequencing identified a mutation in the 3' splice site of exon 8. Northern blot confirmed lower mRNA expression level in Hlb320 and did not show a difference in splicing, suggesting that the mutation only affects the splicing rate. In addition to affecting HDL cholesterol, the mutated genes also lead to reduction in serum non-HDL cholesterol and triglyceride levels. Despite low HDL cholesterol levels, the mice from both mutant lines show similar atherosclerotic lesion sizes compared to control mice. These new mutant mouse models are valuable tools to further study the role of these genes, their affect on HDL cholesterol levels, and metabolism.     

Molecular genetic analysis of two loci (Ity2 and Ity3) involved in the host response to infection with Salmonella typhimurium using congenic mice and expression profiling

Numerous genes have been identified to date that contribute to the host response to systemic Salmonella Typhimurium infection in mice. We have previously identified two loci, Ity2 and Ity3, that control survival to Salmonella infection in the wild-derived inbred MOLF/Ei mouse using a (C57BL/6J x MOLF/Ei)F(2)cross. We validated the existence of these two loci by creating congenic mice carrying each quantitative trait locus (QTL) in isolation. Subcongenic mice generated for each locus allowed us to define the critical intervals underlying Ity2 and Ity3. Furthermore, expression profiling was carried out with the aim of identifying differentially expressed genes within the critical intervals as potential candidate genes. Genomewide expression arrays were used to interrogate expression differences in the Ity2 congenics, leading to the identification of a new candidate gene (Havcr2, hepatitis A virus cellular receptor 2). Interval-specific oligonucleotide arrays were created for Ity3, identifying one potential candidate gene (Chi3l1, chitinase 3-like 1) to be pursued further. The combination of the use of congenics in QTL confirmation and fine mapping and in the identification of candidate genes by expression profiling has been successful and represents a step toward quantitative gene(s) identification.     

Lith6: a new QTL for cholesterol gallstones from an intercross of CAST/Ei and DBA/2J inbred mouse strains

A complex genetic basis determines the individual predisposition to develop cholesterol gallstones in response to environmental factors. We employed quantitative trait locus/loci (QTL) analyses of an intercross between inbred strains CAST/Ei (susceptible) and DBA/2J (resistant) to determine the subset of gallstone susceptibility (Lith) genes these strains possess. Parental and first filial generation mice of both genders and male intercross offspring were evaluated for gallstone formation after feeding a lithogenic diet. Linkage analysis was performed using a form of multiple interval mapping. One significant QTL colocalized with Lith1 [chromosome (chr) 2, 50 cM], a locus identified previously. Significantly, new QTL were detected and named Lith10 (chr 6, 4 cM), Lith6 (chr 6, 54 cM), and Lith11 (chr 8, 58 cM). Statistical and genetic analyses suggest that Lith6 comprises two QTL in close proximity. Our molecular and genetic data support the candidacy of peroxisome proliferator-activated receptor gamma (Pparg) and Slc21a1, encoding Pparg, and the basolateral bile acid transporter SLC21A1 (Slc21a1/Oatp1), respectively, as genes underlying Lith6.     

Fine localization of a new cataract locus, Kec, on mouse chromosome 14 and exclusion of candidate genes as the gene that causes cataract in the Kec mouse

A mouse with cataract, Kec, was generated from N-ethyl-N-nitrosourea (ENU) mutagenesis. Cataract in the Kec mouse was observable at about 5 weeks after birth and this gradually progressed to become completely opaque by 12 weeks. Dissection microscopy revealed that vacuoles with a radial or irregular shape were located primarily in the cortex of the posterior and equatorial regions of the lens. At the late stage, the lens structure was distorted, but not ruptured. This cataract phenotype was inherited in an autosomal recessive manner. We performed a genetic linkage analysis using 133 mutant and 67 normal mice produced by mating Kec mutant (BALB/c) and F1 (C57BL/6 x Kec) mice. The Kec locus was mapped to the 3 cM region encompassed by D14Mit34 and D14Mit69. In addition we excluded coding sequences of 9 genes including Rcbtb2, P2ry5, Itm2b, Med4, Nudt15, Esd, Lcp1, Slc25a30, and 2810032E02Rik as the candidate gene that causes cataract in the Kec mouse.     

Characterization of the quantitative trait locus for haloperidol-induced catalepsy on distal mouse chromosome 1

We report here the confirmation of the quantitative trait locus for haloperidol-induced catalepsy on distal chromosome (Chr) 1. We determined that this quantitative trait locus was captured in the B6.D2- Mtv7a /Ty congenic mouse strain, whose introgressed genomic interval extends from approximately 169.1 to 191.3 Mb. We then constructed a group of overlapping interval-specific congenic strains to further break up the interval and remapped the locus between 177.5 and 183.4 Mb. We next queried single nucleotide polymorphism (SNP) data sets and identified three genes with nonsynonymous coding SNPs in the quantitative trait locus. We also queried two brain gene expression data sets and found five known genes in this 5.9-Mb interval that are differentially expressed in both whole brain and striatum. Three of the candidate quantitative trait genes were differentially expressed using quantitative real-time polymerase chain reaction analyses. Overall, the current study illustrates how multiple approaches, including congenic fine mapping, SNP analysis and microarray gene expression screens, can be integrated both to reduce the quantitative trait locus interval significantly and to detect promising candidate quantitative trait genes.     

Genome-wide microarray expression analysis of CD4+ T Cells from nonobese diabetic congenic mice identifies Cd55 (Daf1) and Acadl as candidate genes for type 1 diabetes

NOD.Idd3/5 congenic mice have insulin-dependent diabetes (Idd) regions on chromosomes 1 (Idd5) and 3 (Idd3) derived from the nondiabetic strains B10 and B6, respectively. NOD.Idd3/5 mice are almost completely protected from type 1 diabetes (T1D) but the genes within Idd3 and Idd5 responsible for the disease-altering phenotype have been only partially characterized. To test the hypothesis that candidate Idd genes can be identified by differential gene expression between activated CD4+ T cells from the diabetes-susceptible NOD strain and the diabetes-resistant NOD.Idd3/5 congenic strain, genome-wide microarray expression analysis was performed using an empirical Bayes method. Remarkably, 16 of the 20 most differentially expressed genes were located in the introgressed regions on chromosomes 1 and 3, validating our initial hypothesis. The two genes with the greatest differential RNA expression on chromosome 1 were those encoding decay-accelerating factor (DAF, also known as CD55) and acyl-coenzyme A dehydrogenase, long chain, which are located in the Idd5.4 and Idd5.3 regions, respectively. Neither gene has been implicated previously in the pathogenesis of T1D. In the case of DAF, differential expression of mRNA was extended to the protein level; NOD CD4+ T cells expressed higher levels of cell surface DAF compared with NOD.Idd3/5 CD4+ T cells following activation with anti-CD3 and -CD28. DAF up-regulation was IL-4 dependent and blocked under Th1 conditions. These results validate the approach of using congenic mice together with genome-wide analysis of tissue-specific gene expression to identify novel candidate genes in T1D.     

Microarray analysis of broad-spectrum resistance derived from an <Emphasis Type="Italic">indica</Emphasis> cultivar Rathu Heenati

Rathu Heenati (RHT) is a Sri Lankan rice cultivar that carries a brown planthopper (BPH) resistance gene, Bph3, and shows broad-spectrum resistance to all four biotypes of BPH. The BPH-resistance loci in RHT has been studied extensively and assigned to four different rice chromosomes (3, 4, 6, and 10) by different research groups, but the gene has not been cloned previously. An Affymetrix rice genome array containing 48,564 japonica and 1,260 indica sequences was used to analyze the potential resistance-related genes on the four chromosomes by comparative analysis of the differentially expressed genes between resistant and susceptible rice cultivars exposed to BPH attack. The microarray results showed that at least 17 genes related to induced resistance and at least 193 genes related to constitutive resistance in RHT. On chromosome 3, the AOC4 was hypothesized to be the most important candidate gene. On chromosome 6, no valuable candidate resistance gene was identified in the Bph3 localization region. In the three Quantitative trait locus regions of chromosomes 3, 4, and 10, the numbers of constitutive and induced resistance-related genes found were 17, 26, and 12, respectively. The major probe on chromosome 10 represents a constitutive expression gene with a very high absolute fold-change of 2,588.82. The microarray analysis indicated that BPH resistance in RHT is probably controlled by a series of resistance-related genes. This study provides valuable information for cloning, functional analysis and marker-assisted breeding of these BPH resistance genes.     

Genetic regulation of endotoxin-induced airway disease

To identify novel genes regulating the biologic response to lipopolysaccharide (LPS), we used a combination of quantitative trait locus (QTL) analysis and microarray-based gene expression studies of C57BL/6J x DBA/2J(BXD) F2 and recombinant inbred (RI) mice. A QTL affecting pulmonary TNF-alpha production was identified on chromosome 2, and a region affecting both polymorphonuclear leukocyte recruitment and TNF-alpha levels was identified on chromosome 11. Microarray analyses of unchallenged and LPS-challenged BXD RI strains identified approximately 500 genes whose expression was significantly changed by inhalation of LPS. Of these genes, 28 reside within the chromosomal regions identified by the QTL analyses, implicating these genes as high priority candidates for functional studies. Additional high priority candidate genes were identified based on their differential expression in mice having high and low responses to LPS. Functional studies of these genes are expected to reveal important molecular mechanisms regulating the magnitude of biologic responses to LPS.     

Characterization of Ath29, a major mouse atherosclerosis susceptibility locus, and identification of Rcn2 as a novel regulator of cytokine expression

Ath29 is an atherosclerosis susceptibility locus on chromosome 9 identified in an intercross between C57BL/6 (B6) and C3H/HeJ (C3H) apolipoprotein E-deficient (apoE(-/-)) mice. This locus was subsequently replicated in two separate intercrosses that developed early or advanced atherosclerotic lesions. The objective of this study was to characterize Ath29 through construction and analysis of a congenic strain and identify underlying candidate genes. A congenic line was constructed by introgressing the chromosomal segment harboring Ath29 from C3H.apoE(-/-) into B6.apoE(-/-) mice. Congenic mice developed significantly smaller early and advance atherosclerotic lesions than B6.apoE(-/-) mice. Microarray analysis revealed 317 genes to be differentially expressed in the aorta of congenic mice compared with B6.apoE(-/-) mice. Pathway analysis of these genes suggested the Ca(2+) signaling pathway to be implicated in regulating atherosclerosis susceptibility. Rcn2 is located underneath the linkage peak of Ath29 and involved in Ca(2+) signaling. Multiple single-nucleotide polymorphisms between B6 and C3H mice were detected within and surrounding Rcn2 with one single-nucleotide polymorphism falling within an upstream cAMP response element. Immunostaining demonstrated its expression in atherosclerotic lesions. Knockdown of Rcn2 with small interfering RNAs resulted in significant reductions in both baseline and oxidized phospholipid-induced VCAM-1 and monocyte chemoattractant protein-1 expression by endothelial cells. Ath29 is confirmed to be a major atherosclerosis susceptibility locus affecting both early and advanced lesion formation in mice, and Rcn2 is identified as a novel regulator of cytokine expression.     

Early transposable element insertion in intron 9 of the Hsf4 gene results in autosomal recessive cataracts in lop11 and ldis1 mice

Lens opacity 11 (lop11) is an autosomal recessive mouse cataract mutation that arose spontaneously in the RIIIS/J strain. At 3 weeks of age mice exhibit total cataracts with vacuoles. The lop11 locus was mapped to mouse chromosome 8. Analysis of the mouse genome for the lop11 critical region identified Hsf4 as a candidate gene. Molecular evaluation of Hsf4 revealed an early transposable element (ETn) in intron 9 inserted 61 bp upstream of the intron/exon junction. The same mutation was also identified in a previously mapped cataract mutant, ldis1. The ETn insertion altered splicing and expression of the Hsf4 gene, resulting in the truncated Hsf4 protein. In humans, mutations in HSF4 have been associated with both autosomal dominant and recessive cataracts. The lop11 mouse is an excellent resource for evaluating the role of Hsf4 in transparency of the lens.