<title>Introduction</title> <p>Type 2 diabetes (T2D) afflicts ∼246 million people worldwide, including ∼21 million in the United States (7% of the population); another 54 million Americans are pre-diabetic. If the incidence of T2D continues to increase at the present rate, one in three Americans, and one in two minorities born in 2000 will develop diabetes in their lifetimes ##REF##12966357##[1]##. Direct medical costs associated with diabetes in the United States exceed $132 billion a year ##REF##11239409##[2]##, and consume ∼10% of health care costs in industrialized nations.</p> <p>Peripheral hyporesponsiveness to insulin increases metabolic demands on the insulin-producing β-cells of the pancreatic islets. Many obese individuals are insulin-resistant, but do not become overtly diabetic provided that the increased demand for insulin is effectively met ##REF##16931493##[3]##,##REF##17229948##[4]##. However, if β-cell mass and/or function are insufficient to meet this requirement, overt hyperglycemia and T2D ensue ##REF##1532777##[5]##. In autopsy series of subjects with T2D, total β-cell mass is decreased ##REF##3901180##[6]##,##REF##12502499##[7]##. Primary reductions of β-cell mass predispose to diabetes in rodent models ##REF##11272209##[8]##,##REF##2673897##[9]##,##REF##5012744##[10]## and in autosomal dominant forms of diabetes (e.g., MODY; maturity onset diabetes of youth) ##REF##11272211##[11]##. Such primary reductions might predispose to some instances of T2D.</p> <p>Susceptibility to T2D is strongly inherited as evidenced by the &gt;80% concordance rates in monozygotic twins ##REF##7193616##[12]##,##REF##1813278##[13]##,##REF##8712800##[14]##,##REF##10064093##[15]##, familial aggregation, and ethnic predispositions ##REF##10837696##[16]##. Heritability of sub-phenotypes related to T2D, e.g. insulin resistance and β-cell hypofunction is even higher ##REF##15931378##[17]##. Environmental factors are also important ##REF##15931378##[17]##,##REF##14527304##[18]##. Although several genes for relatively rare monogenic forms of diabetes, including MODY, syndromic (Wolfram syndrome), lipoatrophic, and mitochondrial-inherited diabetes have been identified ##REF##11239409##[2]##,##REF##11317350##[19]##, the underlying genetic bases for the genetically complex T2D, accounting for &gt;95% of diabetes patients, have remained elusive. The identification of susceptibility genes is made difficult by the polygenic nature of the phenotype ##REF##1607066##[20]##, its reflection of convergent, distinct metabolic processes producing identical phenotypes (phenocopies), and the potent gene-gene and gene-environment (e.g. obesity) interactions that characterize the disease. Clear genetic influences on the endophenotypes (intermediate phenotypes) of β-cell mass/function and insulin resistance have been shown, and vary among racial groups. ##REF##7776502##[21]##,##REF##7758230##[22]##,##REF##8422779##[23]##,##REF##12540622##[24]##. Some notable earlier successes (e.g. <italic>PPARG</italic>, <italic>CAPN10</italic>), and a recent series of genome-wide association studies of large numbers of well-phenotyped subjects ##REF##17293876##[25]##,##REF##17460697##[26]##,##REF##16415884##[27]##,##REF##17434869##[28]##,##REF##17463249##[29]##,##REF##17463248##[30]##,##REF##17463246##[31]## have identified T2D susceptibility loci/genes (e.g. <italic>TCF7L2</italic>) whose functions with regard to the implicated phenotypes are uncertain. As no single implicated gene or allele accounts for more than a small fraction of risk to develop T2D, there are still many genes/molecular mechanisms awaiting identification.</p> <p>In mice, there is striking strain-dependent susceptibility to T2D in the context of obesity ##REF##17018838##[32]##. We exploited the differential diabetes susceptibilities of the B6 and DBA strains segregating for the obesity mutation <italic>Lep<sup>ob</sup></italic>
##REF##17018838##[32]## to identify a diabetes susceptibility QTL in B6xDBA progeny and then used congenic lines derived from the implicated interval to clone a candidate gene accounting for the QTL. Similar strategies have been used to identify QTLs (and responsible genes) for other complex phenotypes in mice ##REF##15803197##[33]## such as type 1 diabetes ##REF##10193189##[34]##, diet-induced obesity ##REF##8703121##[35]##, tuberculosis susceptibility ##REF##11197687##[36]##, atherosclerosis ##REF##11438740##[37]##, epilepsy ##REF##10645948##[38]##, schizophrenia ##REF##12431851##[39]## and, also, T2D ##REF##16682971##[40]##,##REF##17426289##[41]##,##REF##16399503##[42]##,##REF##16804088##[43]##.</p> <p>We identified, “Lisch-like” (<italic>Ll</italic>), a novel gene, encoding multiple, tissue-specific transcripts in brain, liver and islets. The functional consequences of the hypomorphic DBA allele (diabetes-prone) in <italic>Lep<sup>ob/ob</sup></italic> mice appear to be late embryonic to early postnatal reductions in β-cell mass due to diminished rates of β-cell replication, some “catch-up” of β-cell mass by 2–3 months, followed by mild glucose intolerance at &gt;6 months of age. These phenotypes are recapitulated in mice with an ENU-induced null allele of <italic>Ll</italic>.</p>
<title>Methods</title> <title>Animal Husbandry</title> <p>Mice were housed in a barrier facility in ventilated Plexiglas cages under pathogen-free conditions at room temperature (22±1°C) with a 12 h light/dark cycle. Mice were weaned at 21 d and given <italic>ad libitum</italic> access to water and 9% Kcal fat Picolab Rodent Chow 20 (Purina Mills; <ext-link ext-link-type="uri" xlink:href="http://www.purinamills.com/">www.purinamills.com/</ext-link>).The high fat diet protocol used in some animals is described below. Columbia University's Institutional Animal Care and Use Committee (IACUC) approved all protocols. After a 4 h morning fast, mice were sacrificed by carbon dioxide asphyxiation and phenotyped for weight, naso-anal length, and glycosuria. Blood was collected by cardiac puncture and aliquoted into microfuge tubes containing an anticoagulant cocktail of 10 µl of 1 mM EDTA and 1.5 mg/ml aprotinin (Sigma A-6279). Plasma and red blood cell pellets were used to measure glucose, insulin, and glycosylated hemoglobin as previously described ##REF##9287054##[110]##. Tissues (skeletal muscle, pancreas/pancreatic islets, liver, brain, hypothalamus, kidney, spleen, heart, visceral fat, retroperitoneal fat) were collected and immediately frozen in liquid N<sub>2</sub>, and stored at −80°C for further studies. Pancreata were dissected under stereoscope, weighed, and fixed in Z-fix zinc-formalin fixative (Anatech; <ext-link ext-link-type="uri" xlink:href="http://www.anatechltdusa.com/">www.anatechltdusa.com/</ext-link>).</p> <title>Genotyping</title> <p>Liver tissue or tail tips were used for genomic DNA isolation according to standard procedures ##REF##3004966##[111]##. A mutation-specific assay was used to confirm that all phenotypically obese animals were <italic>Lep<sup>ob</sup>/Lep<sup>ob</sup></italic> and all lean animals either +/+ or heterozygous at the <italic>Lep</italic> locus ##REF##9192391##[112]## Animals were genotyped using MapPairs Microstaellite Markers (Invitrogen; <ext-link ext-link-type="uri" xlink:href="http://www.invitrogen.com/">www.invitrogen.com/</ext-link>) as previously described ##REF##12438443##[113]##.</p> <title>Mapping T2D-related Phenotypes in B6xDBA F2 Progeny</title> <p>Maps were created using MapMarkerQTL (<ext-link ext-link-type="uri" xlink:href="http://www.broad.mit.edu/genome_software/other/qtl.html">www.broad.mit.edu/genome_software/other/qtl.html</ext-link>) on a dataset representing 404 obese F2 progeny of a B6xDBA cross segregating for <italic>Lep<sup>ob</sup></italic> at 120–150 days of age. The QTL for T2D was most significantly associated with fasting blood glucose, glycosylated hemoglobin, and islet histology in male mice to a region of Chr1, with peak statistical significance at D1Mit110 at 169.6 Mb from the centromere (p&lt;10<sup>−8</sup>) (##FIG##11##Figure 12##). Other QTLs were identified on other chromosomes (for example Chr5 at 78cM), but none had as great an effect on the phenotype or demonstrated consistent effects on all aspects of the phenotype. We tested for interactions for QTLs and identified a modest interaction between the locus on Chr1 and a second locus at D4Mit286 (p = 0.008).</p> <title>B6.DBA Congenic Lines: Creation and Fine Mapping</title> <p>B6.DBA congenic mice were generated by intercrossing <italic>Lep<sup>ob</sup>/Lep<sup>+</sup></italic> B6 X DBA mice from Jackson Laboratory (<ext-link ext-link-type="uri" xlink:href="http://www.jax.org/">www.jax.org/</ext-link>) to generate F1 progeny, followed by backcrossing to the recurrent B6 strain using a “speed congenic” approach in subsequent generations ##REF##10505408##[114]##. At the eighth backcross, a genome scan was performed in all breeders using polymorphic markers at 20 cM intervals. In the mouse line that was continued, all non-contiguous markers outside the DBA interval were homozygous B6. Over the next two generations, there were two recombination events, one that eliminated a telomeric portion of the DBA interval (line 1jc) and one that preserved approximately half of the originally defined DBA interval (line 1jcd). The 1jcd mouse was bred repeatedly to B6 mice, giving rise, by meiotic recombination, to two additional subcongenic lines (1jcdt and 1jcdc) (see ##FIG##0##Figure 1##). Preservation of the phenotypes present in the original B6xDBA and DBAxB6 F2/F3 progeny was assessed by longitudinal and end-point measurements of fasting glucose, insulin, glycosylated hemoglobin and islet morphology. At N12, <italic>Lep<sup>ob/+</sup></italic> mice B6/DBA (B/D) for the respective congenic intervals were intercrossed to produce N12F1 progeny. Obese progeny were used for fine mapping and phenotyping experiments. <italic>Lep<sup>ob/+</sup></italic> animals D/D for the congenic interval were recurrently intercrossed or crossed to B6 <italic>Lep<sup>ob/+</sup></italic> animals to generate <italic>ob/ob Lep<sup>ob</sup></italic>/<italic>Lep<sup>ob</sup></italic> animals with D/D and B/D genotypes for the Chr1 interval, respectively.</p> <title>Studies of Glucose Homeostasis</title> <p>For longitudinal phenotyping studies, mice were fasted for 4 h and restrained for blood collection by a trained individual. Blood was collected from unanesthetized animals by capillary tail bleed into heparinized tubes and stored at −80°C. Glucose was measured with a FreeStyle Flash Blood Glucose Monitor (Abbott; <ext-link ext-link-type="uri" xlink:href="http://www.abbottdiabetescare.com/">www.abbottdiabetescare.com/</ext-link>). Insulin was measured by ultra-sensitive rat insulin ELISA (ALPCO; <ext-link ext-link-type="uri" xlink:href="http://www.alpco.com/">www.alpco.com/</ext-link>). HbA1c was measured by affinity chromatography (Mega Diagnostics; <ext-link ext-link-type="uri" xlink:href="http://www.mega-dx.com/">www.mega-dx.com/</ext-link>). Urine ketones were measured using Chemistrip Test Strips (Roche Diagnostics; <ext-link ext-link-type="uri" xlink:href="http://us.labsystems.roche.com/index.shtml">http://us.labsystems.roche.com/index.shtml</ext-link>). For ipGTT, mice were fasted overnight and 0.5 g/kg body weight of 50% dextrose was administered intra-peritoneally at time 0. Plasma glucose was measured at 15–30 min intervals for 3 h, as above. Terminal phenotypic characterization consisted of measurements of fasting glucose, insulin, glycosuria, and glycosylated hemoglobin as previously described ##REF##9287054##[110]##. To control for stress-induced hyperglycemia at the time of sacrifice, tail blood glucose was also measured by glucometer one day prior to sacrifice.</p> <title>High Fat and “Surwit” Diet Studies</title> <p>High fat chow pellets (#D12492i: 60% kcal from fat, 20% kcal from protein, 20% kcal from carbohydrate) and “Surwit” ##REF##3044882##[115]## (#D12331i; 58% kcal from fat, 16.4% kcal from protein, 25.5% kcal from carbohydrate) (Research Diets; <ext-link ext-link-type="uri" xlink:href="http://www.researchdiets.com/">www.researchdiets.com/</ext-link>) were used as described in the text.</p> <title>Morphometric and β-cell replication analysis of Pancreatic Islets</title> <p>Pancreatic tissues were dissected under stereoscope to avoid contamination with adipose tissue, and weighed.</p> <title>Islet Morphometry</title> <p>Non-overlapping images of longitudinal pancreatic sections were acquired and analyzed using ImageProPlus software version 5.0 (Media Cybernetics; <ext-link ext-link-type="uri" xlink:href="http://www.mediacy.com/">www.mediacy.com/</ext-link>) to calculate insulin-positive area, insulin-positive area as % total area, and number of islets (defined by an area containing a minimum of 8 contiguous insulin-positive cells). For β-cell replication studies, we recorded the number of Ki67-positive or negative, insulin-positive cells. Replication of β-cells was expressed as % of cells (Ki67-positive and insulin-positive)/ total insulin-positive. For replication studies, ∼100 islets were examined per animal from several different non-overlapping sections through the pancreas. ImageProPlus or Image J (1.37 V; NIH) were used to determine the relative area of each section occupied by β-cells or the actual of number of β-cells for each representative longitudinal pancreatic section (50 µm apart) that had been immunochemically stained for insulin as previously described ##REF##11334404##[116]##. We analyzed 5–7 sections from different regions of the pancreas. Apoptosis rates were assessed using the DeadEnd Fluormetric TUNEL System G3250 (Promega; <ext-link ext-link-type="uri" xlink:href="http://www.promega.com/">www.promega.com/</ext-link>) TUNEL assay and cleaved Caspase-3 (Asp175) Antibody 9661S (Cell Signaling Technology; <ext-link ext-link-type="uri" xlink:href="http://www.cellsignal.com/">www.cellsignal.com/</ext-link>).</p> <title>Pancreatic Islet Isolation</title> <p>Pancreatic perfusion and islet collection were performed as previously described ##REF##10969832##[117]##. Each pancreas was perfused via the bile duct with 1.5 mg/ml collagenase P (Roche Applied Science; <ext-link ext-link-type="uri" xlink:href="http://www.roche-applied-science.com/">www.roche-applied-science.com/</ext-link>) and incubated at 37°C for 17 min. Following disaggregation of pancreatic tissue, pancreata were rinsed with M199 medium containing 10% NCS. Islets were collected by density- gradient centrifugation in Histopaque (Sigma-Aldrich; <ext-link ext-link-type="uri" xlink:href="http://www.sigmaaldrich.com/">www.sigmaaldrich.com/</ext-link>) ##REF##10969832##[117]##, and washed several times with M199 medium. For glucose-stimulated insulin release studies ##REF##5333500##[118]##,##REF##2996187##[119]##, islets were incubated overnight in RPMI medium 1640 (Invitrogen).</p> <title>Glucose-Stimulated Insulin Secretion (GSIS)</title> <p>The GSIS procedure has been described previously ##REF##2666106##[120]##. Islets were hand-picked into tissue culture dishes containing cold Kreb's buffer (118.5 mM NaCl, 2.54 mM CaCl<sub>2</sub>, 1.19 mM KH<sub>2</sub>PO<sub>4</sub>, 1.19 mM MgSO<sub>4</sub>, 10 mM HEPES, pH 7.4), and 2% BSA (Sigma-Aldrich), 5.5 mM glucose, and incubated overnight at 37°C. Islets were hand-picked and incubated another 15 min. in Kreb's buffer+BSA, containing 11.2 mM glucose. Hand-picked islets are then resuspended in Kreb's buffer plus BSA, supplemented with 2.8 mM glucose, and shaken at 37°C for 15 min. The pellet was spun down gently and resuspended in triplicate (5–10 islets each) in 500 µl Kreb's buffer, supplemented with glucose at final concentrations of 2.8 mM, 5.6 mM, 11.2 mM or 16.8 mM, or supplemented with 10 mM arginine and incubated for 1 h in a water bath at 37°C with constant shaking (300 rpm). After 1 h incubation, islets were gently pelleted and the supernatant collected and assayed for insulin by ELISA. Islet pellets were dissolved in high salt buffer (2.15 M NaCl, 0.01 M NaH<sub>2</sub>PO<sub>4</sub>, 0.04 M Na<sub>2</sub>HPO<sub>4</sub>, EDTA 0.672 g/L, pH 7.4) and sonicated at 4–5 W for 30 s and DNA concentration was measured using a TKO100 fluorometer (Hoefer; <ext-link ext-link-type="uri" xlink:href="http://www.hoeferinc.com/">www.hoeferinc.com/</ext-link>) with Hoechst #33258 dye (Polysciences; <ext-link ext-link-type="uri" xlink:href="http://www.polysciences.com">www.polysciences.com</ext-link>). Results were expressed as concentration of secreted insulin/[DNA]/h.</p> <title>Testing for Predicted Transcripts in cDNA Pools</title> <p>Putative transcripts, identified from public annotation and local sequencing, were validated by PCR-amplification from tissue-specific cDNA pools prepared from male and female B6 mice. Two cDNA pools were used: 1. An <underline>inclusive cDNA pool</underline> was prepared from E7 and E20 fetuses and P1 pups, and included the following tissues of 60-day old mice: eyes, large intestine, skin, tongue, spinal cord, kidney, testes/ovaries, pancreatic islets, whole brain, hypothalamus, skeletal muscle, and liver. This pool was used for transcript validation. 2. <underline>A diabetes-relevant cDNA pool</underline>, from 90-day old mice, was comprised of only the following tissues and organs: pancreatic islets, whole brain, hypothalamus, skeletal muscle, liver, and adipose tissue. This pool was used to quantify transcripts identified by computational approaches and the microarrays. Nominal intron-spanning primers were generated using the Primer3 program (<ext-link ext-link-type="uri" xlink:href="http://www.genome.wi.mit.edu/cgi-bin/primer/primer3_www.cgi">www.genome.wi.mit.edu/cgi-bin/primer/primer3_www.cgi</ext-link>). Amplification was first performed on the diabetes-relevant pool at an annealing temperature of 60°C. If we detected no PCR-product, we performed gradient temperature PCR on the same pool using eight different annealing temperatures from 58–68°C. Gradient temperature PCR was then used to amplify the inclusive cDNA pool. If no product was detected in this pool, a 2nd set of intron-spanning primers was used before we interpreted negative amplification as failure to substantiate a predicted transcript. Positive amplification products of predicted sizes, and those that did not match the expected sizes, were gel-purified and sequenced for confirmation. The final set of primer-pairs is listed in Real-time qPCR.</p> <title>Microarray Gene Expression Analysis</title> <p>RNA extraction, purification, labeling, hybridization and analysis were performed as described ##REF##14679176##[121]##. 10 BB and 10 DD 21-day old <italic>Lep<sup>ob/ob</sup></italic> 1jc males were dissected and RNA was extracted from hypothalamus, liver, isolated islets, EDL muscle, and soleus muscle. Individually labeled RNA (by mouse and organ) was interrogated with Affymetrix MOE-430A expression arrays. For further details, see legends to ##TAB##0##Table 1## and ##FIG##6##Figure 7##. For all transcripts in the region of interest, where possible, only probes that spanned multiple exons and clearly represented each of the 14 genes in the interval were used. If &gt;1 probe met these conditions, we used only, the probe that gave the strongest signal. Organs were grouped into two groups by genotype and were compared using a two tailed T-test. The Affymetrix probe IDs selected for this analysis are shown in ##SUPPL##2##Table S3##.</p> <title>Real-Time qPCR</title> <p>Effects of the DBA congenic interval on the levels of confirmed transcripts expressed in diabetes-relevant organs were assessed on an organ-specific basis. We made separate pools from 90-day old <italic>Lep<sup>ob/ob</sup></italic> 1jc D/D and B/B mice for each of the diabetes-relevant organs (see above). Each individual organ pool was generated on 2 occasions from 5 mice. RNA was extracted from organs with TRIzol acid-phenol reagent (Invitrogen). 2 µg of RNA were reverse-transcribed using SuperScript III reverse transcriptase (cDNA First Synthesis Kit, Invitrogen) with random hexamer priming. The cDNA was diluted 4-fold using nuclease-free water (QIAGEN; <ext-link ext-link-type="uri" xlink:href="http://www.qiagen.com">www.qiagen.com</ext-link>). 2 µl of diluted cDNA were amplified by PCR in Roche LightCycler. A standard curve for each transcript was generated using cDNA diluted 1∶1, 1∶10, and 1∶100. We assessed the number of mRNA molecules in each sample using the slope and intercepts of PCR product appearance during the exponential phase of the PCR reactions optimized for transcript-specific product using specific primers. Each sample was run in triplicate in the same LightCycler run. Using LightCycler Software, we calculated the crossing point (CP) for each sample. The CP is the first maximum of the second derivative of the fluorescence curve, and is equivalent to the number of cycles at which the fluorescence first exceeds background. In the exponential phase, the relationship between CP and initial transcript concentration is linear. We calculated relative concentration ratios, normalized to actin, as follows:
</p> <p>In this expression, ΔCP<sub>gene</sub> is the CP of the gene in the sample minus the CP of the gene in the relevant reference; ΔCP<sub>hg</sub> is the CP of the housekeeping gene in the sample minus the CP of the housekeeping gene in the reference (“ref”) sample; and η is the efficiency (where 2 is perfectly efficient) as determined by the negative slope of the plot generated when CP is plotted as a function of the log of initial concentration determined in the standard curve. Each CP listed is the mean of CP values of the triplicates for each sample. Results are summarized in ##TAB##0##Table 1##. Primers used are listed in ##SUPPL##3##Table S4## (A).</p> <title>Cloning and Sequencing of <italic>Lisch-like</italic> Isoforms</title> <p>We amplified full-length <italic>Ll</italic> cDNAs from either B6 islets (isolated by us) or from Clontech MTC Panels 1 #636745 and 3 #636757, containing pooled multiple tissue cDNAs from 8–12 week old BALB/c mice and from Swiss Webster embryos. In a final volume of 50 µl, we added 0.5 µl LA Taq (TaKaRa; <ext-link ext-link-type="uri" xlink:href="http://www.takara-bio.com/">www.takara-bio.com/</ext-link>) to a cocktail containing TaKaRa GC Buffer II, 400 µm each dNTP, 1 µl cDNA and 1 µl each primer (300 ng/µl). Primers are listed in ##SUPPL##3##Table S4## (B). Samples were cycled in an MJ Tetrad Thermalcycler (BioRad; <ext-link ext-link-type="uri" xlink:href="http://www.bio-rad.com">www.bio-rad.com</ext-link>) using a Touchdown protocol of a 2 min. extension and decreasing annealing temperature from 60°C to 55°C for 10 cycles, followed by 25 cycles with an annealing temperature of 55°C. Each sample was TOPO TA cloned (Invitrogen) and plated. From all three libraries, a total of 140 colonies were picked and grown overnight in LB buffer. Inserts were amplified by colony PCR and sized by gel-fractionation. Inserts representing each unique size were then sequenced. The isoforms and the exons deleted (Δ): iso1 (intact 10 exons); iso2, Δ6; iso3, Δ4,5,6; iso4, Δ4; iso5, Δ5,6; iso6, Δ9; iso7, Δ5,6,7,8,9.</p> <title>Zebra Fish Analyses</title> <title>A. Zebra Fish Strains and Embryo Culture</title> <p>Zebra fish and embryos were raised, maintained and staged according to standard procedures ##UREF##4##[122]##. The AB* (Eugene, OR) line and Tg(gut GFP)s854 transgenic line (gutGFP; ##REF##12645931##[73]##) were used in natural matings to obtain embryos. The gutGFP line was provided by Didier Stainier. Embryos examined at stages later than 24 hpf were maintained in embryo medium containing 0.003% phenylthiourea to inhibit pigmentation.</p> <title>B. Morpholino Injections</title> <p>Morpholino antisense and control oligonucleotides, listed in ##SUPPL##3##Table S4## (C), were purchased from Gene Tools (<ext-link ext-link-type="uri" xlink:href="http://www.gene-tools.com/">www.gene-tools.com/</ext-link>) and injected into 1–2 cell stage embryos at concentrations from 7–20 ng/embryo as previously described ##REF##11017081##[67]##.</p> <title>C. RT-PCR</title> <p>Total RNA was extracted from morpholino-injected and uninjected sibling embryos at 29 hpf with TRIzol; cDNA was synthesized with SuperScript II Reverse Transcriptase (Invitrogen) using primer-pairs shown in ##SUPPL##3##Table S4## (D).</p> <title>D. Immunofluorescence and RNA <italic>in situ</italic> Hybridization</title> <p>Zebra fish gene sequences were amplified using the primer-pairs shown in ##SUPPL##3##Table S4## (E) and cloned into the PSTBlue-1 vector (Novagen) and used for antisense probe synthesis with T7 RNA polymerase after <italic>XhoI</italic> linearization (<italic>Lsr-like</italic>) and SP6 polymerase following <italic>Bam</italic>H1 linearization (<italic>Lisch-like</italic>). Whole-mount <italic>in situ</italic> hybridization was performed as described ##UREF##5##[123]##. For immunofluorescence, embryos were fixed at room temperature (rt) in 4% paraformaldehyde for 2 h. After fixation, yolks were manually removed and embryos were permeabilized in acetone at −20°C for 7 min. Embryos were washed briefly in PBS +0.1% Triton ×100 (PBSTx) and incubated for 1 h in antibody hybridization buffer (PBSTx with 2% DMSO, 2% BSA and 2% sheep serum). Guinea pig anti-insulin antibody (Biomeda V2024) was diluted 1∶1000 in antibody hybridization buffer and incubated with embryos for 2 h at rt. Following antibody hybridization, embryos were washed extensively with PBSTx and incubated with Cy3-labelled donkey anti-guinea pig secondary antibody diluted 1∶500 in antibody hybridization buffer for 2 h at rt. Embryos were washed extensively with PBSTx and cleared in 80% glycerol/20% PBS. Images of optical sections were captured using a confocal microscope and 2-D projections were generated from optical sections using MetaMorph software.</p> <title>Computational Methods for Evaluating Effect of nsSNPs</title> <p>We used five methods to compute the likelihood of a functional change due to single amino acid substitutions (see ##FIG##8##Figure 9##). SNAP, PolyPhen, and SIFT predict changes in protein function due to a single amino acid substitution. SNAP ##REF##17526529##[57]## is a neural-network based method that considers protein features predicted from sequence (<italic>e.g.</italic>, residue solvent accessibility and chain flexibility). Scores from −9 to +9 are estimates of accuracy of prediction, computed using a testing set of ∼80,000 mutants. A low negative score indicates confidence in prediction of neutrality (functional change absent), whereas a high positive score indicates confidence in prediction of non-neutrality (functional change present). Accuracy was computed for neutrals using the equation below:
</p> <p>PolyPhen considers structural and functional information and alignments. Predictions are sorted into 4 classes: benign, possibly damaging, probably damaging, and unknown.</p> <p>
<bold>SIFT predictions.</bold> SIFT ##REF##12824425##[59]## is a statistical method that only considers alignments. Scores range from 0 to 1. Scores &gt;0.05 indicate neutrality of a substitution.</p> <p>
<bold>PAM250 matrix substitutions.</bold> PAM matrix ##REF##202030##[124]## (Percent Accepted Mutations) reflects frequency of amino acid interchange throughout evolution (by evaluating alignments of proteins in a family). Scores range from a low of −8 for rare substitutions (<italic>e.g.</italic> W to C) to a high of 17 (same residue found in almost all proteins in alignment).</p> <p>
<bold>Percentage in alignment (PROFacc).</bold> The score is reported as the difference in observed percentages of wild-type and mutated residues in alignments against a non-redundant UniProt ##REF##15608167##[125]## and PDB ##REF##10592235##[126]## database (at 80% sequence identity).Scores range from −100 (if the mutant is observed in all instances) to +100 (if the wild type is observed in all instances); 0 if the mutant is observed as often as the wild type. Scores near 0 favor the likelihood of a mutation being neutral.</p> <title>DBA BAC Shotgun Sequencing</title> <p>BAC 95f9 DNA (5 µg) was fragmented to 1–5 kb using a nebulizer supplied with the TOPO Shotgun Subcloning kit (Invitrogen) and checked for size and quantity on an agarose gel. The shotgun library was constructed with 2 µg of sheared DNA. Blunt-end repair, dephosphorylation, ligation into PCR 4Blunt-TOPO vector, and transformation into TOP10 Electrocompetent <italic>E. coli</italic> were performed with the TOPO Shotgun Subcloning kit, following the manufacturer's protocol. Phenol∶chloroform extraction of the dephosphorylated DNA was replaced with Qiagen QIAquick PCR Purification spin columns (QIAGEN). Recombinant colonies were selected by blue/white screening and incubated in LB medium supplemented with 50 µg/ml ampicillin for 20 h at 37°C in 96-well deepwell plates. Plasmid miniprep was conducted in 96-well plates using QIAGEN Turbo Miniprep kits on a QIAGEN BioRobot 9600. DNA sequencing was performed on a 3730xl Genetic Analyzer (Applied Biosystems; <ext-link ext-link-type="uri" xlink:href="http://www.appliedbiosystems.com/">www.appliedbiosystems.com/</ext-link>) using BigDye® Terminator v3.1 Cycle Sequencing Kits with M13 forward and reverse sequencing primers.</p> <title>Statistical Analyses</title> <p>ANOVA and ANCOVA were used to assess effects of genotype in congenic interval. Comparisons at individual time points, or pairs of means were performed using Student's t-test. P values are 2-tailed. The Statistica package (StatSoft; <ext-link ext-link-type="uri" xlink:href="http://www.statsoft.com/">www.statsoft.com/</ext-link>) was used for ANOVAE; Excel (Microsoft, <ext-link ext-link-type="uri" xlink:href="http://office.microsoft.com/en-us/default.aspx">http://office.microsoft.com/en-us/default.aspx</ext-link>) for t-testing.</p> <title>Western Blot</title> <p>Hypothalamic extracts were prepared using M-PER Mammalian Protein Extraction Reagent (Pierce Biotechnology, <ext-link ext-link-type="uri" xlink:href="http://www.piercenet.com/">www.piercenet.com/</ext-link>). Hypothalamic extracts (85 mg for B/B and D/D congenics and 175 mg for wild-type and mutant ENU mice) were resolved by 8% SDS-PAGE, transferred to nitrocellulose membrane (Invitrogen). We generated a set of polyclonal rabbit antibodies (Covance Research Products; <ext-link ext-link-type="uri" xlink:href="http://www.covance.com">www.covance.com</ext-link>) against the predicted ICD, spanning residues 298–401 (exons 7,8) and verified that the α-ICD rabbit antibodies detected the appropriate fusion proteins, with only minor cross-reactivity in cultured cells. We hybridized the blot with anti-LL anti-sera at a dilution of 1∶5,000 in TBS/0.05%Tween/5% milk (TBSTM) or with blocked anti-LL anti-sera diluted 1∶10,000 in TBSTM. To prepare blocked anti-sera, liver sections from C3HeB/FeJ knock-out mice were fixed overnight in phosphate-buffered paraformaldehyde at 4°C and rinsed in PBS. Sections equivalent to one-third of a liver were fragmented and mixed with 1 ml anti-sera diluted 1/1000 in PBS/0.1% Triton. Liver fragments were spun out and the supernatant was used to probe filters from ENU mice. We detected bound antibody with horseradish peroxidase-coupled antibody against rabbit IgG (Amersham Biosciences; <ext-link ext-link-type="uri" xlink:href="http://www.amershambiosciences.com">www.amershambiosciences.com</ext-link>) at a dilution of 1∶5,000 using the SuperSignal West Pico Chemiluminescent Substrate kit (Pierce Biotechnology).</p> <title>Immunohistochemical and Immunofluorescnce Analysis of Pancreatic Islets</title> <title>For β-Cell Replication Studies</title> <p>Pancreata were fixed overnight in 10% formalin, embedded the specimens in paraffin, and consecutive 5 µm-thick sections were mounted on slides. For immunofluoresence and diaminobenzidine (DAB) staining of Ki67 and for insulin immunoreactivity, tissue sections were de-waxed in xylene, hydrated through a descending ethanol series and subjected to an antigen retrieval step using a heated citrate buffer solution. Several longitudinal sections &gt;100 µm apart were used to assess β-cell replication and double staining for the nuclear proliferation marker Ki67 and insulin. Sections were incubated with Novocastra rabbit polyclonal anti-Ki67 antibody (Leica Microsystems; <ext-link ext-link-type="uri" xlink:href="http://www.leica-microsystems.com">www.leica-microsystems.com</ext-link>) diluted 1∶200 and an insulin polyclonal guinea pig anti-swine antibody (Vector Lab; <ext-link ext-link-type="uri" xlink:href="http://www.vectorlabs.com/">www.vectorlabs.com/</ext-link>) diluted 1∶2000 overnight at 4°C.</p> <title>For Immunofluorescence Detection</title> <p>Sections were washed in PBS and incubated with secondary anti-guinea pig IgG (1∶200) and fluorescein isothiocyanate-conjugated rabbit secondary antibody (1∶200) (Vector Labs) for 1 hr and counterstained with DAPI before the addition of mounting medium. Non-overlapping images of longitudinal pancreatic sections were acquired using a Nikon Eclipse microscope and images imported into ImageJ (1.37 V, NIH) to count insulin-positive and Ki67-insulin-positive cells. β-cell replication is expressed as % Ki67-positive+insulin-positive/total insulin-positive cells. For diaminobenzidine staining, sections were incubated with secondary biotinylated rabbit and quinea pig IgG for 1 hr and then subjected to an avidin∶biotyinylated enzyme complex (ABC Kit; Vector Labs) with DAB as substrate. Sections were counterstained with hematoxylin. Images of pancreatic sections were acquired using SpotAdvanced version 5 software (Diagnostic Instruments; <ext-link ext-link-type="uri" xlink:href="http://www.diaginc.com/">www.diaginc.com/</ext-link>) and analyzed using Image Pro Plus software to calculate the % of β-cell area occupied by Ki67-positive cells. We examined 30–50 islets per animal from several non-overlapping sections through the pancreas.</p> <title>Accession Numbers</title> <p>Genbank (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/">www.ncbi.nlm.nih.gov/</ext-link>(Genbank) accession numbers for the <italic>M. musculus</italic> genes: <italic>Lisch-like</italic>, lipolysis-stimulated remnant receptor-related (XM_001473525); Lsr (NM_017405); <italic>Ildr1</italic> (NM_134109); <italic>Tada1l</italic> SPT3-associated factor 42 (NM_030245); <italic>Pogk</italic> pogo transposable element with KRAB domain (NM_175170); <italic>FMO13</italic>, flavin-containing monooxygenase family; FMO-like (XM_136366); <italic>FMO9</italic>, flavin-containing monooxygenase family; FMO-like (NM_172844) <italic>FMO12</italic>, flavin-containing monooxygenase family; FMO-like (XM_136368); <italic>C030014K22Rik</italic>, unknown (NM_175461); <italic>Uck2</italic>, uridine monophosphate kinase (NM_030724); <italic>Tmco1</italic>, membrane protein of unknown function (NM_001039483); <italic>Aldh9a1</italic>, aldehyde dehydrogenase 9, subfamily A1 (NM_019993); <italic>Mgst3</italic>, microsomal glutathione-S-transferase 3 (NM_025569); <italic>Lrrc52</italic>, leucine-rich repeat (LRR) protein of unknown function (NM_00103382); <italic>Rxrg</italic>, retinoid X receptor, gamma (NM_009107); <italic>Lmx1a</italic>, LIM homeobox transcription factor 1, α (NM_033652); <italic>Pbx1</italic> (NM_008783); <italic>H.sapiens C1orf32</italic> (NM_199351); <italic>LSR</italic> (NM_015925); <italic>ILDR1</italic> (NM_175924); <italic>D.rerio Ll</italic> ortholog, zgc:110016 (NM_001030192.1); <italic>D. rerio Lsr</italic> ortholog, zgc:114089 (NM_001025472.1); <italic>R. rattus Lsr</italic> (NM_032616).</p> <p>The Genbank accession numbers for protein sequences: <italic>M. musculus</italic> Lisch-like (amino acid residues 150–795, XP_001473575); (Lsr) (NP_059101); Ildr1 (NP_598870); <italic>H. sapiens</italic> C1orf32 (NP_955383); LSR (NP_057009); ILDR1 (NP_787120); <italic>D. rerio</italic> Lisch-like (NP_001025363); <italic>D. rerio</italic> Lsr (NP_001020643); <italic>R. rattus</italic> Lsr (NP_116005)</p>
<title>Results</title> <title>Genetic Map of Diabetes QTL and Related Congenic Lines</title> <p>We identified a QTL for diabetes-related phenotypes in obese F2 and F3 progeny of an intercross between diabetes-resistant C57BL/6J (B6) and diabetes-susceptible DBA/2J (DBA) mice segregating for <italic>Lep<sup>ob</sup></italic>. Phenotypes including fasting blood glucose, HbA1c and islet histology mapped with LOD &gt;8 around D1Mit110 on distal Chr 1 at 169.6 Mb (details in <xref ref-type="sec" rid="s4">Methods</xref>: Mapping T2D-related Phenotypes). By producing congenic and sub-congenic B6.DBA lines also segregating for <italic>Lep<sup>ob</sup></italic>, we refined the interval to 5.0 Mb between rs31968429 at 168.1 Mb and rs31547961 at 173.1 Mb where all four congenic lines overlap for DBA (##FIG##0##Figure 1##; details in <xref ref-type="sec" rid="s4">Methods</xref>: B6.DBA Congenic Lines: Creation and Fine Mapping).</p> <p>We further restricted the search (##FIG##0##Figure 1##) by identifying a haplotype block ##REF##12466852##[44]## conserved between B6 and DBA that extends 3.2 Mb from rs30708865 at 169.9 Mb to rs31547961 at 173.1 Mb. Only eleven unvalidated B6 vs. DBA single nucleotide polymorphisms (SNPs) in this interval are listed in the Mouse SNP database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/SNP/MouseSNP.cgi">www.ncbi.nlm.nih.gov/SNP/MouseSNP.cgi</ext-link>); however, among fragments we could amplify containing nine of these putative SNPs, we detected no sequence variants. Moreover, we found no coding sequence/expression difference between B6 and DBA among all genes and transcripts in the “conserved” interval by computation, direct sequencing, and quantitative mRNA expression analysis. Thus, it is unlikely that the variant(s) in the genetically-defined interval with peak at 169.6 Mb mediating differential diabetes susceptibility between these two strains is within the “conserved region.” We sequenced the 3 kb interval between rs31968429 and rs33860076 at the centromeric end of subcongenic line 1jcdt and detected no variants between the two strains. Therefore, we focused our efforts on the 1.8 Mb B6 vs. DBA “variable” interval, between rs33860076 at 168.1 Mb and rs30708865 at 169.9 Mb.</p> <title>Metabolic and Anatomic Phenotypes of Congenic Lines</title> <p>The congenic/sub-congenic lines shown in ##FIG##0##Figure 1## displayed phenotypes of hypoinsulinemic hyperglycemia in association with relative reductions in β-cell mass due to reduced β-cell proliferation (see Islet Morphology and β-cell Replication and Apoptosis). Phenotypes were generally more salient in male animals. Genotype in the congenic interval (B6 or DBA) <italic>per se</italic> did not affect their body weight or composition. Supporting experiments are described below.</p> <p>By 4 weeks of age, fasting plasma glucose was elevated in <italic>Lep<sup>ob/ob</sup></italic> males who were D/D (DBA/DBA) for the congenic interval 1jcd and fed standard (9% fat) chow; glucose concentrations were higher up to 120 days. After 120 days, there were no significant differences in fasting glucose between D/D (DBA/DBA) and B/B (B6/B6) mice (##FIG##1##Figure 2A##). The decline in pre-prandial blood glucose levels in <italic>Lep<sup>ob/ob</sup></italic> males between 90 and 200 days is probably attributable to a slight expansion of β-cell mass in response to transient insulin resistance occurring as a normal consequence of sexual maturation (∼60 days of age) ##REF##2673897##[9]##,##REF##2912706##[45]##. To examine diabetes susceptibility in D/D animals that were obese independent of leptin deficiency, we fed lean (<italic>Lep<sup>+/+</sup></italic>) 1jcd males a high-fat diet (60% kcal from fat) for 13 weeks, starting at 7 weeks of age. These mice became more hyperglycemic than B/B mice (##FIG##1##Figure 2B##), showing a persistence of this difference – similar to the animals in 2A – up to age ∼140 days when the study ended.</p> <p>To delineate differences in acute glucose handling in D/D vs. B/B animals, we used intraperitoneal glucose tolerance testing (ipGTT). At 60 days (##FIG##1##Figure 2C##), and even up to 200 days, when the study ended (##FIG##1##Figure 2D##), <italic>Lep<sup>ob/ob</sup></italic> 1jcdc males were less glucose tolerant than B/B. The relative reduction in glucose tolerance in D/D vs. B/B animals that are not overtly diabetic is likely related to reduction in the number of islets. The occurrence of the diabetes-related phenotype is independent of <italic>Lep<sup>ob</sup></italic>, since 100-day old <italic>Lep<sup>+/+</sup></italic> 1jc D/D males fed the Surwit (high fat, high sucrose) diet for 10 weeks were also less glucose tolerant than littermate B/B males (##FIG##1##Figure 2E##).</p> <p>Hyperglycemia due to relative hypoinsulinemia, was evident in 1jc <italic>Lep<sup>ob/ob</sup></italic> D/D animals fed a chow diet as early as 4 weeks (##FIG##2##Figure 3A##). At mean ages of 30- and 62-days, age-adjusted plasma insulin concentrations per mg blood glucose were lower in D/D than in B/B animals. This difference was due to lower plasma insulin in D/D (p = 0.0004) and not higher blood glucose in D/D (p = 0.916). Consistent with these ratios, D/D <italic>Lep<sup>+/+</sup></italic> males showed a 40% decrease in insulin secretion when clamped at a blood glucose level of 250 mg/dl for an hour (##FIG##2##Figure 3B##). No difference in insulin sensitivity was detected by euglycemic – hyperinsulinemic clamping (data not shown).</p> <p>Consistent with their hypoinsulinemic hyperglycemia, 21-day old 1jcd D/D males had smaller islets than their B/B counterparts (##FIG##3##Figure 4A##). A qualitative cell-autonomous β-cell defect in insulin secretion, however, is unlikely to be the primary functional defect in D/D animals, since islets isolated from 28-day old 1jcd D/D males responded to graded glucose concentrations (2.8 mM–16.8 mM) or 10 mM arginine by secreting amounts of insulin comparable to age- and sex-matched B/B littermates (##FIG##3##Figure 4B##). Also consistent with insulin/glucose ratios and hyperglycemic clamp results, isolated islets from 60-day old 1jc <italic>Lep<sup>ob/ob</sup></italic> males fed normal chow and 100-day old 1jc <italic>Lep</italic> +/+ on the Surwit diet showed reduced insulin secretion at 2.8 mM and 5.6 mM [glucose] in D/D vs. B/B littermates. For reasons indicated below, the early glucose intolerance of D/D mice is probably due, in part, to a deficiency of β-cell mass.</p> <title>Islet Morphology and β-cell Replication and Apoptosis</title> <p>The fractional area of the pancreas accounted for by β-cells ##REF##10642598##[46]## in <italic>Lep<sup>ob/ob</sup></italic> 1jcd males was examined in 20-, 60- and 150-day old mice. By 60 days a trend to reduced β-cell area was apparent in D/D, and by 150 days of age, β-cell mass of the 1jcd D/D sub-congenics was about half that of B/B littermate controls. B/D animals had β-cell masses that were about two-thirds of B/B littermate controls (##FIG##4##Figure 5A##). These findings are consistent with <italic>in vivo</italic> data showing onset of elevated blood glucose (see ##FIG##1##Figure 2A##) and lower circulating insulin concentrations (relative to glucose) in D/D sub-congenics at ∼60 days of age (see ##FIG##2##Figure 3A##), and persistence of decreased glucose tolerance at 200 days of age. The lower relative β-cell mass in D/D animals reflects fewer numbers of β-cells, rather than smaller sized β-cells. There were no differences in pancreatic weight between D/D and B/B male animals.</p> <p>To assess the basis for the difference in β-cell mass by 60 days, we measured rates of β-cell replication and apoptosis. We co-stained pancreatic sections in 1jcd congenic 1- and 21-day old <italic>Lep<sup>ob/ob</sup></italic> male mice with antibodies to insulin and Ki67 antigen, a nuclear marker of proliferation expressed during all stages of the cell cycle except G0 ##REF##14599612##[47]##. To estimate the proportion of dividing β-cells, we normalized the number of Ki67 positive β-cells to the total number of insulin positive cells. Groups consisted of 4 B/B and 4 D/D 1-day old mice or 4 B/B, and 8 D/D 21-day old mice. In 1-day old D/D males, the rate of β-cell replication was ∼1/3 that of B/B littermates, whereas there was no difference in 21-day old animals due to normally reduced β-cell replication by the time of weaning (##FIG##4##Figure 5B##) ##REF##11042468##[48]##,##REF##10830272##[49]##,##REF##11272192##[50]##.</p> <p>The proportion of small islets (250–2000 µm<sup>2</sup>) in 21-day old <italic>Lep<sup>ob/ob</sup></italic> males was greater in D/D (1jc and 1jcd) mice (73%) than in B/B (60%); whereas the proportion of large islets (10,000–50,000 µm<sup>2</sup>) was lower (9% in D/D and 14% in B/B). This finding is consistent with the β-cell replication studies in P1 mice (##FIG##4##Figure 5B##), and recently reported evidence that new β-cells are derived from replication of pre-existing β-cells ##REF##15129273##[51]##.</p> <p>In 13-day old 1jc mice, when β-cell apoptosis is active ##REF##9075738##[52]##, we did not detect significant differences between B/B and D/D islets in β-cell apoptosis using a TUNEL assay ##REF##1400587##[53]## and caspase-3 staining ##REF##7983002##[54]## (data not shown). Thus, the lower number of β-cells in D/D mice is primarily a result of lower rates of proliferation of β-cells in the perinatal period.</p> <title>Genes in the Minimal DBA Interval Conveying Diabetes Susceptibility</title> <p>To identify all genes in the minimal DBA variable interval, (see above for definition) we screened 277 genes and transcripts, computationally predicted by GenScan, TwinScan, FGeneSH, Otto, or SGP2 that map to the interval. We excluded 50 single-exon transcripts (probably pseudogenes ##REF##12951575##[55]##) that did not belong to a transcript cluster and were not homologous to transcripts in the syntenic human interval, and 16 ribosomal gene transcripts, unique to this interval, that could not be specifically amplified due to their genomic redundancy, and manually curated the remaining 211 predicted transcripts. We rejected 63 that did not amplify in RNA/cDNA pools from multiple organs/ages of B6 and DBA mice (see <xref ref-type="sec" rid="s4">Methods</xref>: Testing for Predicted Transcripts in cDNA Pools) and, using BLASTn, clustered the remaining 148 transcripts into 14 groups. These, correspond to 11 known genes and 3 predicted genes that we validated by amplification in cDNA pools.</p> <p>A map of the “variable” interval shows 14 genes, flanked by <italic>Mael</italic> and <italic>Pbx1</italic> (##FIG##5##Figure 6##). We analyzed all transcripts in the entire “variable” region.</p> <title>Analysis of Genes in the Variable Interval</title> <p>The genetic variation accounting for differential diabetes-susceptibility in mice segregating B/B vs. D/D in the congenic intervals could be due to: 1) coding sequence variant(s) that alter the amino acid sequence of a protein(s); 2) regulatory variants, including anti-sense transcripts that affect expression and stability, and 3′ untranslated region (UTR) variants; or 3) splicing variants. We investigated all hypotheses.</p> <title>Non-Synonymous Sequence Variants</title> <p>To identify all non-synonymous B6 vs. DBA sequence variants in the “variable” interval, we collected genomic sequence for B6 and DBA strains from databases at NCBI and Celera ##REF##11783003##[56]##, filled gaps using bi-directional sequencing to achieve 100% coverage of all coding sequences in both strains, and validated coding sequence variants by bi-directionally re-sequencing gene fragments encompassing each variant in both B6 and DBA strains. Consequently, we identified five non-synonymous single nucleotide variants: one in each of three FMO-like (flavin mono-oxygenase) genes, and two in chr1.1224.1 (##FIG##5##Figure 6##). The latter gene, we designated “<italic>Lisch-like</italic>” (<italic>Ll</italic>) because of its sequence similarity to a gene in mouse and rat, formerly known as <italic>Lisch</italic>7 (<ext-link ext-link-type="uri" xlink:href="http://rgd.mcw.edu/">http://rgd.mcw.edu/</ext-link>), but now known as <italic>Lsr</italic> (lipolysis stimulated receptor).</p> <p>Computational analysis of LL and the three FMO-like proteins using SNAP ##REF##17526529##[57]##, PolyPhen ##REF##12202775##[58]##, SIFT ##REF##12824425##[59]##, PAM250 matrix substitution weights ##UREF##0##[60]## and PROFacc ##REF##8066087##[61]## predicted that all of the amino acid substitutions were benign with respect to function. The SNAP scores obtained for our variant alleles, -1 (FMO13, K282E), -2 (FMO12, V239I), -3 (LL, A647V), and -6 (LL, T587A; FMO9, Q5R), indicate that there is a ∼60%, ∼69%, ∼79%, and ∼90% respective chance of the non-synonymous variants being neutral. Similarly, PolyPhen classified all variations as “benign” and SIFT scores were well above 0.05 (neutral). PAM weights of 0 and above suggest interchangeability of the respective amino acids throughout evolution. The % differences were low, suggesting that the DBA and B6 variants are equally likely to occur in related sequences (see <xref ref-type="sec" rid="s4">Methods</xref>: Computational Methods for Evaluating Effects of nsSNPs).</p> <title>Expression Differences</title> <p>We used Affymetrix microarrays to quantify those transcripts in the minimum congenic interval that we had validated by PCR-amplification (see <xref ref-type="sec" rid="s4">Methods</xref>: Testing for Predicted Transcripts in cDNA Pools). We examined hypothalamus, islets, liver, soleus and EDL (extensor digitorum longus) skeletal muscle from DD and BB <italic>Lep<sup>ob/ob</sup></italic> congenic animals (see <xref ref-type="sec" rid="s4">Methods</xref>: Microarray Gene Expression Analysis). These arrays did not contain elements for all of the 14 genes we confirmed in the interval: missing from the array were the 3 FMO genes. Therefore, we also used real-time qPCR, to quantify expression of each gene and confirmed transcript in tissues and organs central to diabetes (pancreatic islets, liver, skeletal muscle, adipose tissue and hypothalamus) in 90-day old male <italic>Lep<sup>ob/ob</sup></italic> 1jc D/D and B/B animals (see <xref ref-type="sec" rid="s4">Methods</xref>: real time qPCR). Results of the microarray and qPCR experiments are shown in ##TAB##0##Table 1## and summarized in ##FIG##6##Figure 7A##.</p> <p>Among genes in the region, including <italic>Lmx1a</italic>
##REF##7698771##[62]##, and <italic>Rxrg</italic>
##REF##16893912##[63]##, that constitute candidates for susceptibility to T2D, we identified no non-synonymous SNPs (nsSNPs) and no multi-organ differences in expression levels between B/B and D/D animals. The most prominent and consistent differences in expression we did observe were for chr1.1224.1 (<italic>Ll</italic>), which was two to four-fold lower in 21-day old <italic>Lep<sup>ob/ob</sup></italic> D/D mice than in B/B mice in the diabetes-relevant tissues/organs by microarray analysis and up to twenty-fold lower by qPCR (##FIG##6##Figure 7A##). (We later show that <italic>Ll</italic> protein in hypothalamus is strikingly reduced in 1jc D/D vs. B/B; see ##FIG##10##Figure 11A##). The difference in <italic>Ll</italic> gene expression in liver persists with age (##FIG##6##Figure 7B##) as does the difference in glucose tolerance in response to overt glucose challenge (see ##FIG##1##Figure 2D##). Whether the differences in hepatic <italic>Ll</italic> expression are mechanistically related to differences in glucose homeostasis are unknown at this point; LL may influence hepatic gluconeogenesis, or the hepatic differences could simply mirror parallel and more physiological relevant changes in β-cells.</p> <p>We also detected (by PCR) <italic>Ll</italic> transcripts in e7, e11, e15, and e17 whole mouse embryos, and in testis, kidney, heart, lung, uterus, eye, thymus and spleen. For the anti-sense interval between intron 9 and intron 7 (see below and ##FIG##0##Figures 1## and ##FIG##7##8##), we found higher expression levels in liver and hypothalamus of D/D v. B/B animals. This difference is consistent with a possible suppressive role for the D/D anti-sense transcript (see below). The <italic>Aldh9a1</italic> gene, known to be highly expressed in human embryonic brain and involved in glycolysis and fatty acid metabolism, showed qualitative changes comparable to those seen in <italic>Ll</italic>. The mapping experiment that identified the interval of mouse Chr1 containing statistical signals related to T2D phenotypes would be expected to enrich for regions in which several genes might contribute to the phenotypes. Although <italic>Aldh9a1</italic> may be such a gene, we chose to focus initially on <italic>Ll</italic>, since it showed the most striking quantitative differences in expression between D/D and B/B animals.</p> <title>“Lisch-like” (<italic>Ll</italic>) Gene Structure and Splice Variants</title> <title>Complete Gene Sequence</title> <p>To identify 3′ and 5′ UTRs flanking the isolated transcripts of Lisch-like (<italic>Ll</italic>), we mapped each transcript onto the UCSC Mouse Genome Browser (<ext-link ext-link-type="uri" xlink:href="http://genome.ucsc.edu/">http://genome.ucsc.edu/</ext-link>) and included contiguous 5′ and 3′ ESTs. The <italic>Ll</italic> gene spans 62,714 bp on mouse Chr. 1, from 168,090,795–168,153,508 (##FIG##7##Figure 8##). The full-length, 10-exon transcript, isoform 1 (iso1), is 8,279 nucleotides. It comprises a 301 nt 5′ non-coding sequence, a 1941 nt coding sequence (including stop codon), encoding a 646 amino acid polypeptide, and a 6,037 nt 3′ UTR. The predicted protein includes a cleavable, signal peptide (SP; exon 1), an extra-cellular domain (ECD; exons 2–4), a trans-membrane domain (TMD; the amino-half of exon 5), and a large intra-cellular domain (ICD; from the cysteine-rich, carboxy-half of exon 5–exon 10). Exons 2 and 3 of the ECD are immunoglobulin-like (Ig-like) V-type domains. Exon 6 is proline-rich and the ICD is overall serine/threonine-rich.</p> <title>Isoforms</title> <p>We isolated complete transcripts for 7 isoforms of <italic>Ll</italic> by PCR amplification of cDNAs using primer-pairs flanking the first and last predicted exons (see <xref ref-type="sec" rid="s4">Methods</xref>: Cloning and Sequencing of Lisch-like Isoforms). We identified 4 major isoforms shown in ##FIG##7##Figure 8## and 3 minor isoforms. Exons 5 and 6 are absent in iso5; exon 9 is absent in iso6; and exons 5–9 are absent in iso7.</p> <title>5′ Upstream Interval</title> <p>The 5′ upstream interval shown (##FIG##7##Figure 8A##) includes 569 nt upstream of the predicted first transcribed base of the 5′ UTR. A CpG island is predicted to overlap the 5′ UTR. By sequencing this interval in DBA BAC 95f9 (MM_DBA library, Clemson University Genomics Institute; <ext-link ext-link-type="uri" xlink:href="http://www.genome.clemson.edu/">www.genome.clemson.edu/</ext-link>), we discovered 8 DBA vs. B6 nucleotide variants not in the public database. Of these, only one variant, (a C to T substitution within a CpG island) is outside a repeat element.</p> <title>Anti-sense Interval</title> <p>An unspliced 2,845 nt anti-sense transcript (##FIG##7##Figure 8B##) of <italic>Ll</italic>, from adult male mouse B6 pituitary gland (5330438I03Rik; red bar in ##FIG##0##Figure 1##), starts 42 bp telomeric of exon 9, crosses exons 9 and 8, and terminates in the intron between exons 7 and 8. This transcript (see ##FIG##6##Figure 7A##) is expressed 2–3 fold higher in DBA vs. B6 in hypothalamus and liver. The centromeric end of the anti-sense transcript is just 506 bp from rs33860076 at the centromeric end of the region of DBA overlap among congenic lines 1jcd, 1jcdt and 1jcdc. An open reading frame (ORF) encodes a predicted polypeptide of 271 amino acids, but with no identifiable domain, and homologous only to ORF segments in anti-sense strands of <italic>Ll</italic> in other species. The interval contains 45 DBA vs. B6 variants, five of which, underlying exon 9, are listed in dbSNP. One newly discovered variant in the intron preceding exon 8, is an insertion in DBA of a 37 nt unique sequence that is homologous to a sequence in an intron of the mouse otoancorin gene on chromosome 7 and to an intronic sequence of an N-scan predicted gene on chromosome 11.</p> <title>3′ UTR</title> <p>Of 52 B/D sequence variants in the long (6 kb) 3′ UTR of the <italic>Ll</italic> transcript (##FIG##7##Figure 8C##), 20 were newly discovered by our “in-house” sequencing.</p> <title>Cross-Species Comparisons of <italic>Ll</italic> Sequence</title> <p>From the Ensembl database, we identified zebra fish orthologs of <italic>Ll and Lsr</italic>. The clustalW pair-wise similarity scores for the predicted protein coded for by the zebra fish gene zgc:114089 (<italic>Lsr</italic> ortholog) is 42 vs, the mouse LSR protein, and 29 vs. the mouse LL protein. The similarity scores for the predicted protein coded for by the zebra fish gene zgc:110016 (<italic>Lisch-like</italic> ortholog) are 36 vs. LL and 28 vs. LSR. We performed clustalW analysis (##FIG##8##Figure 9##) between the mouse LL-iso1 protein and three related proteins: 1) the human C1orf32 protein at 1q24.1 (chr.1 165,154,620–165,211,185; NCBI Build 36.1), which is the product of a gene highly expressed in the developing human retina and brain ##UREF##1##[64]##; 2) the predicted protein sequence for the zebra fish <italic>Lisch-like</italic> ortholog, zgc:110016 located on zebra fish chromosome 9 at 31.6 Mb; and 3) the mouse <italic>LSR</italic> protein, transcribed from a gene on chromosome 7 at 30.7 Mb. Pair-wise similarity scores for the intact proteins and major domains are shown in the legend. The human homolog is similar throughout, but diverges slightly in the putative ICD. The zebra fish Lisch-like ortholog and mouse LSR proteins are most alike in the TMD, less so in the Ig-like domain, and most dissimilar in the ICD. The <italic>Lsr</italic> protein has a short extension to exon 6, and no exon 8 equivalent. <italic>Ll</italic> and <italic>Lsr</italic> also have splicing patterns similar to the mouse <italic>Ildr1</italic> (Ig-like domain receptor 1) gene ##REF##15381095##[65]##, and the proteins they encode all belong to the Lisch7 family (IPR008664; <ext-link ext-link-type="uri" xlink:href="http://www.ebi.ac.uk/interpro">www.ebi.ac.uk/interpro</ext-link>).</p> <title>Knockdown of <italic>Ll</italic> and <italic>Lsr</italic> Orthologs in Zebra Fish</title> <p>To assess the function of <italic>Ll</italic> in islet/β-cell ontogenesis, we examined expression patterns and the effects of morpholino-mediated knockdown in zebra fish embryos. Morpholinos are modified anti-sense oligonucleotides that produce a strong hypomorphic “knockdown” phenotype ##REF##11477696##[66]## either by inhibiting proper splicing of the pre-RNA transcript ##REF##11477696##[66]## or by ATG-blocking of translation ##REF##11017081##[67]##. Morpholino knockdown has been used to demonstrate a role for the endocrine hormones GnRH, GHRH and PACAP during development ##REF##16099813##[68]##,##REF##12941629##[69]##,##REF##15862551##[70]##,##REF##16731775##[71]##. Many of the molecular mechanisms regulating pancreas development appear to be conserved among zebra fish and other vertebrates ##REF##15602891##[72]##, and the single zebra fish islet provides an excellent model of vertebrate development.</p> <p>Using whole mount <italic>in situ</italic> hybridization (##FIG##9##Figure 10A##), we observed that the <italic>Lisch-like</italic> ortholog zgc:110016 was expressed in the brain and otocyst by 48 hours post fertilization (hpf), and by 72 hpf expression was evident in the intestine. The <italic>Lsr</italic> ortholog zgc:114089, located on Chr 15 at 39.0 Mb, was expressed in pancreas at 48 and 72 hpf, (similar to our postnatal observations in mouse with <italic>Ll</italic>), intestine, liver, pharynx, pronehphros and otocyst for 48 hpf (72 hpf not shown), and, at 34 hpf, in both pancreatic buds. Since the anterior bud gives rise to exocrine tissue, pancreatic duct, and a small number of endocrine cells, while the posterior bud gives rise only to endocrine tissue ##REF##12941629##[69]##, expression of the <italic>Lsr-like</italic> paralog throughout this stage is consistent with a role in the ontogeny of pancreatic endocrine tissue.</p> <p>The close structural similarities among Lisch-related genes (see ##FIG##8##Figure 9##) suggested that functional data on both zebra fish genes could be physiologically relevant and, therefore, we studied the involvement in islet development of both orthologs. We injected (in separate experiments) morpholinos for both genes into embryos homozygous for the gut-GFP (green fluorescent protein) transgene to visualize developing endodermal organs (##FIG##9##Figure 10B##) ##REF##12645931##[73]##. We assessed β-cell development with an anti-insulin antibody at 48 hpf or by insulin <italic>in situ</italic> hybridization at 24 hpf (not shown). To assess morpholino specificity, we analyzed the effects of two separate, non-overlapping morpholinos for each gene. Both morpholinos for each ortholog independently produced similar phenotypes, providing evidence that the effects (described below) were the result of specific gene knockdown and not due to nonspecific morpholino-related effects.</p> <p>
##FIG##9##Figure 10B## shows that both <italic>Lsr-like</italic> and <italic>Ll</italic> morpholinos injected at 15 ng/embryo produced general developmental delay in the endodermal organs, evidenced by a smaller liver, a smaller, straighter intestine, and a smaller pancreas that does not extend as much as in wild-type. The <italic>Lsr-like</italic> morpholinos disrupt β-cells more severely (note ectopic insulin-positive cells in the cephalad region of the pancreas) than do the <italic>Ll</italic> morpholinos (note the milder local dispersion of insulin-positive cells); 48/72 and 25/144 embryos injected with morpholinos targeting <italic>Lsr-like</italic> and <italic>Ll</italic>, respectively, displayed a scattered β-cell phenotype. These effects were rarely observed in uninjected sibling embryos (0/25) or embryos injected with a control morpholino (1/35). Lower doses of <italic>Lsr-like</italic> and <italic>Ll</italic> morpholinos (∼7–10 ng) resulted in a lower frequency of β-cell scattering and higher doses (∼20–25 ng) resulted in embryonic toxicity, which is common with high doses of morpholinos. The efficacy of the splice-blocking <italic>Lsr-like</italic> and <italic>Ll</italic> morpholinos was assessed via RT-PCR and all were found to strongly and specifically inhibit proper splicing of their respective target transcripts at the 15 ng dose (not shown). In combination, the expression analyses and morpholino knockdown studies provide support for a role of <italic>Lisch</italic> gene family members in endodermal development, and suggest specific effects on the embryonic β-cell. The relevance of such zebra fish studies to mammalian pancreas development has been shown earlier for <italic>Ptf1a</italic>
##REF##15031114##[74]##,##REF##15570689##[75]## and for <italic>Pdx1</italic>
##REF##11477692##[76]##.</p> <title>W87* Stop Mutation of <italic>Ll</italic> in C3HeB/FeJ Mice</title> <p>To examine phenotypes of mice segregating for a null allele for <italic>Ll</italic>, we screened a repository of ENU-generated (N-ethyl-N-nitrosourea) mutant sperm DNAs from 18,000 C3HeB/FeJ G1 males (Ingenium; <ext-link ext-link-type="uri" xlink:href="http://www.ingenium-pharmaceuticals.com/">http://www.ingenium-pharmaceuticals.com/</ext-link>) for mutations in <italic>Lisch-like</italic>
##REF##16075367##[77]##. We detected a G/A substitution that encodes an amber stop mutation at threonine-87 [W87*] and also creates an EcoN1 cleavage site, which we used to genotype for the mutation. By <italic>in vitro</italic> fertilization, we generated W87* heterozygotes on the C3HeB/FeJ background, and bred these animals to generate progeny that were homozygous wild-type (+/+), homozygous mutant (−/−) or heterozygous (+/−) for the W87* mutation. Progeny were born at the anticipated Mendelian ratios, and the −/− animals did not appear grossly compromised.</p> <p>To verify that the W87* homozygous mutant was hypomorphic for LL protein, we compared a Western blot of hypothalamic extracts prepared from C3HeBFeJ wild-type (+/+) and mutant (−/−) mice, with a second blot of hypothalamic extracts prepared from B/B and 1jc-D/D congenic mice. We probed both sets of filters with a polyclonal rabbit antibody generated to a conjugated polypeptide, corresponding to exons 7 and 8 of isoform 1, in the predicted ICD of LL. As anticipated, LL protein was greatly reduced in the brains of D/D vs B/B congenics and in the ENU-treated W87* homozygotes vs. the wild-type animals (##FIG##10##Figure 11A##).</p> <p>In mice at 14 days of age we can detect reductions in β-cell replication rates that are similar to those seen in the DD congenic lines (##FIG##4##Figure 5B##) There is a &gt;2-fold difference in the proportion of Ki67-positive β-cells in 14-day old wild-type (3.75%) vs. homozygous W87* mice (1.75%), with heterozygotes intermediate (2.5%) (##FIG##10##Figure 11B##). Plasma insulin concentrations in <italic>Ll</italic> W87* homozygotes are reduced by the time of sexual maturation (##FIG##10##Figure 11C##) and, consistent with this difference, at 50 days of age, homozygous W87* males show an increased glucose AUC during iPGTT (##FIG##10##Figure 11D##). A significant decrease in β-cell mass is also detected in W87* homozygotes (1.05%±.117, n = 3, p = .0113) v. +/+ littermates (2.74±.364; n = 3) at 150 days of age.</p> <p>It is important to note that these phenotypes were detected despite the segregation of the mutation on a different background strain (C3HeB/FeJ) than our congenics (C57BL/6J), and in the absence of co-segregation of the <italic>Lep<sup>ob</sup></italic>. These preliminary data strongly support the candidacy of <italic>Ll</italic> as the gene accounting for the diabetes-related phenotypes of the DD congenic lines.</p>
<title>Discussion</title> <p>Based upon a QTL analysis of modifiers of T2D in B6xDBA F2 <italic>Lep<sup>ob/ob</sup></italic> mice, we identified a novel gene, Lisch-like (<italic>Ll</italic>), whose apparent effect on β-cell development, and possibly other aspects of β-cell/islet biology, qualify it as a strong candidate mediator of susceptibility to T2D. On the C57BL/6J strain background, the presence of the DBA/2J congenic interval(s) produced mild hypoinsulinemic hyperglycemia (in association with reduced β-cell replication and mass). Our preliminary data in ENU-mutagenized mice with a null <italic>Ll</italic> allele are consistent with a role for LL in β-cell development.</p> <p>Three of the <italic>Ll</italic> subcongenic lines (1jcd, 1jcdt and 1jcdc) contain only DBA DNA 3′ of exon 7, while line 1jc is DBA for the entire gene and extends DBA for another 3 Mb 5′ of <italic>Ll</italic>. We infer, therefore, that coding and/or non-coding DBA vs. B6 variant(s) in the region of DBA overlap accounts for the phenotypic differences between the DBA congenic lines and animals segregating for B6 alleles in this region. In the region of overlap that includes the DBA vs. B6 “variable region” (##FIG##5##Figure 6##), <italic>Ll</italic> is the only gene showing anticipated differences in coding sequence and gene expression. These findings strongly support, but do not prove, the putative role of <italic>Ll</italic> alleles in conveying the phenotypic differences seen between the various DD and BB congenic lines. The phenotypes of the <italic>Ll</italic> W87* C3H mice also support our inferences regarding the candidacy of <italic>Ll</italic> based upon the B.D congenics.</p> <p>There are two non-synonymous SNPs in <italic>Ll</italic> within the region of overlap among the congenic lines, in exon 9. However, their effects on protein function are predicted to be minor and it is unlikely that they determine the differences in either transcript abundance or protein level seen in the congenics. Variants in other regions of the gene are likely more relevant.</p> <p>In the 5′ UTR, all but one of the eight variants are in simple repeats, where they are likely less significant. The interval underlying the anti-sense transcript contains 45 D/B variants, including a long, unique insertion. A regulatory role for the <italic>Ll</italic> anti-sense transcript is suggested by the similar location of anti-sense transcripts at the 3′ ends of the human <italic>C1orf32</italic> (human ortholog of <italic>Ll</italic>) gene (<italic>e.g.</italic>, DA322725 from hippocampus), the human <italic>LSR</italic> gene (DA320945, also from hippocampus), the human <italic>ILDR1</italic> gene (AW851103), and the mouse <italic>Lsr</italic> gene (BY747866). Moreover, comparative inter-species transcriptomic analysis has identified the 3′ regions of transcripts as important in anti-sense regulation, and conserved overlap between species may be evidence of function ##REF##17250976##[78]##. For a recent review of anti-sense regulatory mechanisms, see ##REF##17139297##[79]##.</p> <p>We identified 52 B/D variants in the 3′ UTR, and it is estimated that the stability of 35% of yeast transcripts are regulated by motifs in the 3′ UTR ##REF##16207357##[80]##. Regulatory motifs, at a similar density, have been identified in the 3′ UTRs of several mammals, including mice ##REF##15735639##[81]##. A 3′ UTR polymorphism between two putative mRNA destabilizing motifs in <italic>PPPIR3</italic> (muscle-specific glycogen-targeting regulatory PP1 subunit) has been genetically ##REF##9726244##[82]## and functionally ##REF##10479482##[83]## related to T2D. Variants in the 3′ UTR may also affect regulation by microRNAs (miRNAs). The 3′ UTR is the target of mammalian microRNAs (miRNAs) ##REF##17612493##[84]## and their relevance to diabetes is underscored by the finding that mouse islet-specific miR-375 affects insulin secretion ##REF##15538371##[85]##.</p> <p>The physiological role of <italic>Ll</italic> is unknown. Based upon the effects of D alleles of <italic>Ll</italic> on β-cell proliferation rates, β-cell mass, <italic>in vivo</italic> insulin release and glucose tolerance, (##FIG##4##Figure 5##) it is likely that <italic>Ll</italic> influences early β-cell differentiation/turnover in a manner that predisposes obese animals to later failure of β-cells by effects on mass and possibly function ##REF##16823478##[86]##,##REF##17259975##[87]##. The fact that these phenotypes are substantially recapitulated in W87* <italic>Ll</italic> C3H mice supports this inference.</p> <p>In the neonatal rodent, extensive remodeling of β-cells occurs as a result of simultaneous activation of both apoptosis and β-cell replication ##REF##10830272##[49]##. Between 4 and 24 weeks, postnatally, β-cell mass is estimated to increase 10 fold, related in part to increased body mass ##REF##10830272##[49]##. Compensation for β-cell stress/loss in adult rodents is primarily by β-cell hypertrophy and β-cell proliferation ##REF##15129273##[51]##. In rats, β-cell proliferation rates decline from ∼20% per day in pups, to ∼10% per day at 6–8 weeks, and to ∼2% shortly thereafter ##REF##7883109##[88]##. However, even this low rate of turnover apparently does not persist in adulthood. Using continuous long term BrdU labeling in B6x129Sv and BALB/C one year-old mice, replacement rates as low as ∼1/1400 mature β-cells/day have been reported ##REF##16123343##[89]##. Consistent with this finding, pancreas mass in the mouse was recently shown to be irreversibly constrained by the size of a progenitor pool in the embryonic pancreatic bud ##REF##17259975##[87]##. These data suggest that β-cell mass established in the first 6–8 weeks of life may be critical to the ability to meet subsequent stresses on β-cell function imposed by e.g. obesity, hyperglycemia, and dyslipidemia. The molecular regulation of these processes is incompletely understood, but even transient interruptions may, based upon this formulation, result in permanent effects on cell mass, or function, or both ##REF##11809615##[90]##. Hypoactivity of the candidate T2D modifier gene (<italic>Ll</italic>) reported here could mediate such effects on establishment of initial β-cell mass, and/or later responses of cell hypertrophy/replication by β-cell-autonomous effects or in response to an exogenous ligand for this putative receptor.</p> <p>Observations that expression levels of <italic>Ll</italic> are most strikingly affected in liver, the effects of the zebra fish knockdowns on general endodermal development, and structure/function considerations raised by the homologous LSR molecule ##REF##10224102##[91]##, are consistent with the possibility that the mechanism(s) by which <italic>Ll</italic> conveys effects on cell mass/function might relate, in part, to consequences of putative effects on hepatic development/ function. IGF1 ##REF##2407519##[92]## and hepatic growth factor ##REF##10625667##[93]## are examples of such β-cell “hepatokines” affecting β-cell function.</p> <title>Similarities to Trans-Membrane Receptors LSR AND ILDR1</title> <p>Insight into the function(s) of the mouse <italic>Lisch-like</italic> protein may be gained from similarities in structure, expression, and cellular location with the human paralog, <italic>C1orf32</italic>, and with genes encoding related trans-membrane receptors, <italic>Ildr1</italic>
##REF##15381095##[65]## and <italic>Lsr</italic>
##REF##10224102##[91]##. Splicing patterns of these genes generate isoforms, similar to those of <italic>Ll</italic>. Each gene's largest isoform includes an extra-cellular Ig-like domain, a single TMD, and a similar set of ICDs in related order. In one isoform of each protein, the TMD and cysteine-rich domains are absent. An evolutionary, regulatory relationship is suggested by the observation that the <italic>Ll</italic>-paralog and <italic>lldr1</italic> are adjacent in the zebra fish genome (Zv6 assembly, UCSC Genome Browser). All three genes are abundantly expressed in the brain, liver and pancreas (and islets, where studied), and all are predicted to have 14-3-3 interacting domains (thus far experimentally verified for the human <italic>LSR</italic>) ##REF##15324660##[94]##. Although 14-3-3 interacting domains may be present on as many as 0.6% of human proteins, their occurrence on all of these Lisch-related proteins is notable, since among known 14-3-3-interacting proteins is phoshodiesterase-3B, which is implicated in diabetes and pancreatic β-cell physiology ##REF##12453887##[95]##,##REF##14693720##[96]##,##REF##14744259##[97]##, and others, such as the Cdc25 family members, important in regulating cell proliferation and survival ##REF##15161933##[98]##,##REF##16697662##[99]##.</p> <title>T2D Genetics for Region of <italic>C1orf32</italic>: Chr1q23</title> <p>The human ortholog of <italic>Ll</italic>, <italic>C1orf32</italic>, which is 90% identical to <italic>Ll</italic> at the amino acid level, maps to a region of Chr1q23 that has been implicated in T2D in seven ethnically diverse populations including Caucasians (Northern Europeans in Utah) ##REF##10331426##[100]##, Amish Family Study ##REF##12540634##[101]##,##UREF##2##[102]##, United Kingdom Warren 2 study ##REF##11484155##[103]##, French families ##REF##11067779##[104]##, and Framingham Offspring study ##REF##11872688##[105]##, Pima Indians ##REF##9758619##[106]##, and Chinese ##REF##14693720##[96]## with LOD scores as high as 4.3. The mouse congenic interval examined here is in the middle of, and physically ∼10× smaller than, the 30 Mb human interval. Recent analysis of the broad interval ascertained in Utah identified two peaks, one of which, at D1S2762 (at 163.6 Mb), is just 12 kb telomeric to the 5′ end of <italic>C1orf32</italic>
##REF##17425920##[107]##. The genes, and gene order, are generally conserved between mouse and human in the region syntenic to the congenic interval. The metabolic phenotypes documented in human subjects with T2D linked to 1q23 resemble diabetic phenotypes observed in congenic mice segregating for the DBA interval in B6.DBA congenics examined here ##UREF##3##[108]##, suggesting that the diabetes-susceptibility gene in congenic mice and human subjects may be the same gene, or among the genes, acting in the same genetic pathway(s). The syntenic interval in the Goto-Kakizaki (GK) rat also correlates with diabetes-susceptibility ##REF##9169130##[109]##.</p> <title>Summary</title> <p>We report the molecular cloning and preliminary characterization of a candidate gene for a mouse QTL modifying T2D phenotypes in mice. The gene, Lisch-like, is novel in structure among diabetes susceptibility genes, and appears to modify β-cell development. Amino acid sequence analysis is consistent with the possibility that hypomorphism for this gene could affect β-cell development by a number of possible molecular mechanisms. Proof of the role of this gene in the imputed phenotypes and molecular processes awaits its further analysis in transgenic animals and cell-based systems.</p>
<p>Conceived and designed the experiments: MDC WKC JCC EW NW CVW JO CEA KKB RI CL KS RLL. Performed the experiments: MDC WKC JCC EW NW CVW JO SG GMS CEA KKB RI CL KS RLL. Analyzed the data: JY BW YB SGF RLL. Wrote the paper: MDC WKC JCC KS SGF RLL.</p> <p>In 404 <italic>Lep<sup>ob/ob</sup></italic> F2 progeny of a C57BL/6J (B6) x DBA/2J (DBA) intercross, we mapped a DBA-related quantitative trait locus (QTL) to distal Chr1 at 169.6 Mb, centered about D1Mit110, for diabetes-related phenotypes that included blood glucose, HbA1c, and pancreatic islet histology. The interval was refined to 1.8 Mb in a series of B6.DBA congenic/subcongenic lines also segregating for <italic>Lep<sup>ob</sup></italic>. The phenotypes of B6.DBA congenic mice include reduced β-cell replication rates accompanied by reduced β-cell mass, reduced insulin/glucose ratio in blood, reduced glucose tolerance, and persistent mild hypoinsulinemic hyperglycemia. Nucleotide sequence and expression analysis of 14 genes in this interval identified a predicted gene that we have designated “Lisch-like” (<italic>Ll</italic>) as the most likely candidate. The gene spans 62.7 kb on Chr1qH2.3, encoding a 10-exon, 646–amino acid polypeptide, homologous to <italic>Lsr</italic> on Chr7qB1 and to <italic>Ildr1</italic> on Chr16qB3. The largest isoform of <italic>Ll</italic> is predicted to be a transmembrane molecule with an immunoglobulin-like extracellular domain and a serine/threonine-rich intracellular domain that contains a 14-3-3 binding domain. Morpholino knockdown of the zebrafish paralog of <italic>Ll</italic> resulted in a generalized delay in endodermal development in the gut region and dispersion of insulin-positive cells. Mice segregating for an ENU-induced null allele of <italic>Ll</italic> have phenotypes comparable to the B.D congenic lines. The human ortholog, <italic>C1orf32</italic>, is in the middle of a 30-Mb region of Chr1q23-25 that has been repeatedly associated with type 2 diabetes.</p> <title>Author Summary</title> <p>Type 2 diabetes (T2D) accounts for over 90% of instances of diabetes and is a leading cause of medical morbidity and mortality. Twin studies indicate a strong polygenic contribution to susceptibility within the context of obesity. Although approximately ten genes making important contributions to individual risk have been identified, it is clear that others remain to be identified. In this study, we intercrossed obese, diabetes-resistant and diabetes-prone mouse strains to implicate a genetic interval on mouse Chr1 associated with reduced β-cell numbers and elevated blood glucose. We narrowed the region using molecular genetics and computational approaches to identify a novel gene we designated “Lisch-like” (<italic>Ll</italic>). The orthologous human genetic interval has been repeatedly implicated in T2D. Mice with an induced mutation that reduces <italic>Ll</italic> expression are impaired in both β-cell development and glucose metabolism, and reduced expression of the homologous gene in zebrafish disrupts islet development. <italic>Ll</italic> is expressed in organs implicated in the pathophysiology of T2D (hypothalamus, islets, liver, and skeletal muscle) and is predicted to encode a transmembrane protein that could mediate cholesterol transport and/or convey signals related to cell division. Either mechanism could mediate effects on β-cell mass that would predispose to T2D.</p>
<title>Supporting Information</title>
<p>This manuscript is dedicated to our friend and associate, Dr. Marija Dokmanovic-Chouinard, who died on 25 July 2006. Her efforts were critical to the project reported here. We thank Rick Rausch, Brandi Davis, Rosa Rosario and Qiong Li for technical assistance, Michael Rosenbaum for help with statistical analyses, Streamson C. Chua for technical assistance and helpful discussions, Larry L. Shapiro for technical help with antibodies, and Domenico Accili for helpful advice. Derek J. Nunez (GlaxoSmithKline) provided helpful consultations.</p>
<fig id="pgen-1000137-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g001</object-id><label>Figure 1</label><caption><title>Genetic Map of Sub-Congenic Lines for Diabetes-Related Phenotypes in the Interval Chr 1 164–194 Mb.</title><p>Genetic map shows sub-congenic lines (1jc, 1jcdc, 1jcd, 1jcdt) in the interval Chr1:164–194 Mb that display hypoinsulinemic hyperglycemia in association with histological evidence of a relative reduction in β-cell mass in the first 21–28 days of life due to reduced β-cell proliferation. An expanded view of the <italic>Ll</italic> gene (chr1.1224.1) is shown at bottom. Above the map scale, in black type, are microsatellite markers that were used to genotype B6 and DBA alleles to establish general boundaries of these congenic intervals. D1mit110 is the peak of the F2/F3 QTL linkage map (see Mapping T2D-related Phenotypes in B6xDBA F2/F3 Progeny). Below map scale, RefSNP (rs) and D-markers in red type identify DBA sequence limits of the respective congenic lines. Markers in blue type identify the closest, confirmed non-DBA (B6) sequence. Sequences in intervals between markers in red and blue type are DBA vs. B6 invariant. Gray bars are DBA-derived sequences. Yellow box corresponds to a 3.2 Mb interval, conserved between DBA and B6. The red box identifies the N-scan predicted gene, chr1.1224.1, subsequently identified as <italic>Lisch-like</italic> (<italic>Ll</italic>), extending centromerically from line 1jcdt. In the expanded view of <italic>Ll</italic>, the B6 boundary (rs31968429) for lines 1jcdc, 1jcd, 1jcdt is 333 bp centromeric of exon 7; the DBA boundary, (rs33860076) is 2,700 bp telomeric of exon 7. 5330438I03Rik is an anti-sense transcript described in detail in the text. Marker positions are from the mouse genome annotation (NCBI Build 36, February 2006).</p></caption></fig> <fig id="pgen-1000137-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g002</object-id><label>Figure 2</label><caption><title>Fasting Blood Glucose and Glucose Tolerance Tests in Congenic Lines.</title><p>A) Blood glucose in <italic>Lep<sup>ob/ob</sup></italic> males congenic for the interval 1jcd (BB; homozygous for B6; DD; homozygous for DBA in the interval) and fed regular mouse chow diet (9% fat) <italic>ad libitum</italic>. Determinations made following a 4 h morning fast. 4–13 animals per genotype group. B) Blood glucose in <italic>Lep<sup>+/+</sup></italic> males congenic for the interval 1jcd and fed high fat diet (60% of calories as fat) <italic>ad libitum</italic> for 13 wks, starting at 7 wks of age. Determinations made following a 4 h morning fast. C) ipGTT in 60-day old <italic>Lep<sup>ob/ob</sup></italic> males congenic for the interval 1jcdc. D) ipGTT in 200-day old <italic>Lep<sup>ob/ob</sup></italic> males congenic for the interval 1jcdc. E) ipGTT in 14-wk old male <italic>Lep<sup>+/+</sup></italic> males congenic for the interval 1jc who had been fed the “Surwit” diet for 10 wks. In all panels, * indicates p&lt;0.05 (2 tailed t-test for genotype effect of congenic interval; mean is +/−SEM.</p></caption></fig> <fig id="pgen-1000137-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g003</object-id><label>Figure 3</label><caption><title>Plasma Insulin/Glucose Ratios and Hyperglycemic Clamps in Age-Grouped 1jc Congenic Males.</title><p>A) Fasting plasma insulin/glucose ratios in 30- and 62-day old 1jc <italic>Lep<sup>ob/ob</sup></italic> B/B and D/D male mice, chow-fed since weaning. Raw data are shown in ##SUPPL##0##Table S1##. Asterisk (*) indicates significant difference between B/B and D/D animals; p-value &lt;0.05 for 2-tailed t-test. B) Hyperglycemic clamping in 100-day old 1jc males on Surwit Diet for 18 weeks. 1jc DD male mice fed a Surwit diet for 18 wks were clamped at a blood glucose concentration of 250 mg/dl for 1 hr and serum insulin concentrations measured at 1 hr. Asterisk (*) indicates p-value &lt;0.05 for 2-tailed t-test.</p></caption></fig> <fig id="pgen-1000137-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g004</object-id><label>Figure 4</label><caption><title>Relationship between Islet Histology and Insulin Secretion.</title><p>A) Islet histology in 21-day old 1jcd <italic>Lep<sup>ob/ob</sup></italic> B/B and D/D male mice. 4 µm pancreatic sections from 21-day old 1jcd <italic>Lep<sup>ob/ob</sup></italic> B/B and D/D male mice were insulin-stained with anti-guinea pig IgG and visualized by light microscopy at 10× magnification. In D/D animals, islets were smaller and less numerous. By histomorphometry, the proportion of small islets (250–2000 µm<sup>2</sup>) in 21-day old <italic>Lep<sup>ob/ob</sup></italic> males was greater in D/D (1jc and 1jcd) mice (73%) than in B/B (60%); whereas the proportion of large islets (10,000–50,000 µm<sup>2</sup>) was lower (9% in D/D and 14% in B/B). B) <italic>In vitro</italic> glucose-stimulated insulin secretion in pancreatic islets in 28-day old 1jc <italic>Lep<sup>ob/ob</sup></italic> B/B and D/D males. Each congenic genotype group consisted of 3 male animals. Negative control consisted of 3 4-week old diabetes-prone <italic>Lepr<sup>db/db</sup></italic> KsJ male animals that are hypo-responsive to glucose stimulation ##REF##2673897##[9]##; positive control was 3 4-week old insulin-resistant animals segregating for a diabetes-susceptibility QTL on Chr5 at 78cM, characterized by hyperglycemia and hyperinsulinemia. B/B and D/D show dose response, but no B/B vs. D/D difference at any concentration of glucose. Response to 10 mM arginine in the same animals confirms that the β-cells of the B/B and D/D congenics are comparable with regard to insulin release to a non-glucose stimulus. The 0 mM arginine control in B/B is shown to establish baseline insulin levels.</p></caption></fig> <fig id="pgen-1000137-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g005</object-id><label>Figure 5</label><caption><title>β-Cell Mass and Replication Rates in 1jcd <italic>Lep<sup>ob/ob</sup></italic> Males.</title><p>A) Relative β-cell area in 20-, 60, and 150-day old <italic>Lep<sup>ob/ob</sup></italic> males segregating for B and D 1jcd congenic intervals. In 60 and 150-day old males segregating for the D/D sub-congenic interval, relative β-cell mass was approximately half that of B/B littermates; B/D animals were intermediate at 150 days. N = 10 for each of the 3 groups of animals. Mean+/−SEM. The asterisk (*) indicates that p&lt;0.05 for D/D vs. B/B at 60 days, and D/D and B/D vs B/B at 150 days. These findings are consistent with <italic>in vivo</italic> data (see ##FIG##1##Figure 2C##) showing onset of elevated blood glucose at rest and during ipGTT by 60 days. B) β-cell replication rates (Ki67) in 1- and 21-day old <italic>Lep<sup>ob/ob</sup></italic> B/B and D/D 1jcd males. To estimate the proportion of dividing cells, the number of Ki67-positive β-cells was normalized to the total number of insulin-positive cells. Replication of β-cells in 1-day old D/D males was ∼1/3 that of B/B littermates (p = 0.017). This difference, not present in 21-day old animals, was probably due to normally reduced β-cell replication by the time of weaning. Mean+/−SEM. The asterisk (*) indicates that p&lt;0.05 for D/D vs. B/B in 1-day old animals.</p></caption></fig> <fig id="pgen-1000137-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g006</object-id><label>Figure 6</label><caption><title>Genes in the Minimal Congenic Interval on Chr1:168.1–170.3 Mb.</title><p>Gray background corresponds to the minimal DBA “variable” interval from 168.1 Mb–169.9 Mb, between markers rs33860076 and rs30708865. Yellow background corresponds to the centromeric end of the DBA vs. B6 “conserved” interval (i.e. nominally invariant). Genes in blue are from RefSeq; genes in black are predicted and locally confirmed as described in the text. The N-scan predicted gene chr1.1224.1 is designated here as “Lisch-like”. Amino acid variants are shown in red to the right of the corresponding gene. Nucleotide substitutions were confirmed by bidirectional sequencing in both C57BL/6J and DBA/2J DNA.</p></caption></fig> <fig id="pgen-1000137-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g007</object-id><label>Figure 7</label><caption><title>Expression Analysis of Candidate Genes and Liver Expression of <italic>Lisch-like</italic>.</title><p>A) Tissue-specific expression analysis of genes in the “variable” portion of the minimum congenic interval. Data for relative expression (B/B to D/D) from ##TAB##0##Table 1## for hypothalamus, islets, liver and EDL-muscle are displayed graphically and numerically below the graph. 21-day old DD and BB <italic>Lep<sup>ob/ob</sup></italic> 1jc congenic males were analyzed using Affymetrix #430A microarrays. B) Liver expression of <italic>Lisch-like</italic> in 1jc B/B and D/D males from 21–120 days. Samples from <italic>Lep<sup>ob/ob</sup></italic> 1jc males were analyzed by qPCR.</p></caption></fig> <fig id="pgen-1000137-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g008</object-id><label>Figure 8</label><caption><title>Predicted Structure of <italic>Ll</italic> Gene with Expanded Views of Critical Regions.</title><p>
<italic>Lisch-like</italic> gene (middle of figure) is the full-length, 10-exon, splice variant (iso1) and includes 872 bp upstream of the transcriptional start site. <underline>Predicted</underline> domains are below exons. Exon 1 includes the 5′ UTR (narrow orange bar) and cleavable signal peptide (SP). Exons 2–4 are extra-cellular, within which exons 2–3 code for an Ig-like domain. Exon 5 includes the TMD with a very cysteine-rich cluster in the carboxyl half; exons 6–10 code for a serine- and proline-rich intracellular domain; exon 10 also includes a long 3′ UTR. The “xs” identify exons deleted in isoforms 2–4. A. 5′ upstream interval (expanded view); Black bars correspond to BLAT displays vs. the reference B6 genome. DBA variants are below the DBA bar. Annotations are composites of displays from the UCSC Genome Browser on Mouse February 2006 Assembly. Royal blue peaks correspond to sequences with predicted ESPERR regulatory potential (<ext-link ext-link-type="uri" xlink:href="http://www.bx.psu.edu/projects/esperr">http://www.bx.psu.edu/projects/esperr</ext-link>). The darker blue blocks correspond to evolutionarily conserved sequences, from the UCSC “Conservation” track. Lighter blue blocks show positions of simple sequence motifs, with the consensus motif shown below. The CpG island track, provided by the UCSC Genome Browser, was generated using the unpublished <italic>cpglh</italic> program from Washington University (St. Louis) Genome Sequencing Center. B. Anti-sense interval corresponds to the sequences overlapping the Riken transcript 5339438I03Rik. The green triangle identifies a 37 nt unique sequence insertion in DBA. The two non-synonymous sequence variants in exon 9 sense transcript are shown. The SNP rs33860076 is the DBA marker at the centromeric end of congenic lines 1jcdc, 1jcd, and 1jcdt. C. 3′ UTR interval; vertical black bars represent positions of 52 B6 vs. DBA nucleotide sequence variants.</p></caption></fig> <fig id="pgen-1000137-g009" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g009</object-id><label>Figure 9</label><caption><title>ClustalW Analysis of <italic>Lisch-like</italic> Homologs and the <italic>LSR</italic> Protein.</title><p>ClustalW analysis was performed on the EMBL-EBI server (<ext-link ext-link-type="uri" xlink:href="http://www.ebi.ac.uk/clustalw/">www.ebi.ac.uk/clustalw/</ext-link>) using their default settings. We modified the display to emphasize exonic alignments. Positions of the two non-synonymous variants in exon 9 of <italic>Ll</italic> are identified by blue shading. The non-homologous extension of mouse <italic>Lsr</italic> exon 6 (pink background) is shown beneath exon 7. Mm_Ll; <italic>Mus musculus Lisch-like</italic>; Hs_C1orf32, <italic>Homo sapiens C1orf32</italic>; Dr_Ll, <italic>Danio rerio</italic> (zebra fish) <italic>Lisch-like</italic> ortholog; Mm_LSR, <italic>Mus musculus LSR-like</italic> ortholog. Pair-wise similarity scores by isoform and domain are shown in ##SUPPL##1##Table S2##.</p></caption></fig> <fig id="pgen-1000137-g010" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g010</object-id><label>Figure 10</label><caption><title>Expression Patterns and Morpholino Knockdown in Zebra Fish Embryos.</title><p>A) Developmental expression of zebra fish <italic>Lisch-like</italic> and <italic>Lsr-like</italic> orthologs. <italic>Lisch-like</italic> RNA was hybridized <italic>in situ</italic> to whole-mount zebra fish embryos at 48 hours post-fertilization (hpf), dorsal view with anterior towards the top; and 72 hpf, lateral view with anterior towards the top, ventral towards the right and yolk removed. <italic>Lsr-like</italic> RNA was hybridized at 48 hpf and 34 hpf. <italic>Ll</italic> panels show ventral views of embryos with yolks removed and anterior towards the top. <italic>Lsr-like</italic> panels show the same image captured in the focal plane of the anterior (ap) and posterior (pp) pancreatic buds, respectively. i, intestine; ph, pharynx; pn, pronephric ducts; l, liver; ap, anterior pancreatic bud; pp, posterior pancreatic bud; p, pancreas (after anterior and posterior bud fusion); b, brain; o, otic vesicle. B) Morpholino knockdown of <italic>Lisch-like</italic> and <italic>Lsr-like</italic> orthologs at 48 hpf. Two dimensional ventral views (anterior towards top) of confocal stacks of 48 hpf embryos, uninjected or injected with 15 ng morpholino: control, <italic>Lsr-like</italic> sp1, and <italic>Lisch-like</italic> ATG. Gut-GFP transgene expression (green); insulin immunolabelling (red).</p></caption></fig> <fig id="pgen-1000137-g011" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g011</object-id><label>Figure 11</label><caption><title>Phenotypes of Mice Segregating for the W87* Allele of <italic>Lisch-like</italic>.</title><p>A) Western analysis of <italic>Lisch-like</italic> in hypothalamus of 1jc and homozygous W87* mice. The Western immunoblot shows differences in <italic>Ll</italic> expression in hypothalami of 1jc <italic>Lep<sup>ob/ob</sup></italic> B/B vs. D/D congenic males (left panel), and between wild-type C3HeB/FeJ and W87* C3HeB/FeJ males (right panel). The right panel immunoblot was incubated with rabbit anti-LL antiserum, prepared against a polypeptide corresponding to exons 7 and 8 of the ICD. The antiserum had been absorbed to fixed liver extracts from knock-out mice in order to block non-specific proteins from interacting with the antibody. The LL transcript isomers are visible as a 65 and 70 kD doublet in the B/B and C3HeB/FeJ wild-type lanes, but absent in the lanes of the 1jc-D/D congenic and C3HeB/FeJ W87* homozygous ENU mutants. B) Percent Replicating β-cells in 14-day old ENU-mutagenized mice. The percentage of Ki67-positive β-cells was estimated in 14-day old C3HeB/FeJ ENU-mutagenized mice, who were either homozygous wild-type (+/+), heterozygous (+/−), or homozygous for the W87* LL amber mutation (−/−). At 14 days there was a 2-fold difference in the % of Ki67<sup>+</sup> β-cells in +/+ (3.75%) vs. −/− (1.75%) ENU W87* mice; +/− were intermediate (2.5%). Non-overlapping images of longitudinal pancreatic sections (200 µm apart) were acquired and analyzed using ImageJ software version 1.37 (NIH) to count insulin-positive and Ki67<sup>+</sup> cells. Pancreatic weights of +/+ and −/− were not different. C) Fasting blood glucose (squares) and insulin/glucose ratios (diamonds) in W87* (−/−) and wild-type (+/+) littermates. P-value &lt;0.05 for 2-tailed t-test at 63 days of age. Data points at other ages show trends. D) ipGTT on 50-day old Surwit-fed B6.CH3. N3F1 W87* males. Glucose intolerance is seen in W87* mice. Mice were fasted overnight prior to dextrose injection (50% dextrose solution, 0.5 g/kg, ip). Capillary tail bleeds were performed at the specified time points to determine circulating glucose levels by glucometer (FreeStyle Flash, Abbott). Blood glucose concentrations that are marked with an asterisk are significantly different (t-test; p&lt;0.05; mean±SEM). Area under curve +/+ vs. −/− (p = 0.02).</p></caption></fig> <fig id="pgen-1000137-g012" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.g012</object-id><label>Figure 12</label><caption><title>LOD Scores for Markers along Chromosome 1.</title><p>LOD scores are shown for fasting blood glucose (black) and pancreatic grade (blue). Terminal phenotypes by genotype at D1Mit110 at 169.6 Mb are summarized in ##SUPPL##4##Table S5##.</p></caption></fig>
<table-wrap id="pgen-1000137-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000137.t001</object-id><label>Table 1</label><caption><title>BB/DD Transcript Ratios of Genes in the Variable Interval 168.1–169.9 Mb.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="10" align="left" rowspan="1">BB/DD Transcript ratio</td></tr><tr><td align="left" rowspan="1" colspan="1">Confirmed Genes<xref ref-type="table-fn" rid="nt101">a</xref>
</td><td colspan="2" align="left" rowspan="1">Liver</td><td colspan="2" align="left" rowspan="1">Brain</td><td colspan="2" align="left" rowspan="1">Islets</td><td align="left" rowspan="1" colspan="1">EDL</td><td align="left" rowspan="1" colspan="1">Soleus</td><td align="left" rowspan="1" colspan="1">Muscle</td><td align="left" rowspan="1" colspan="1">Adipose</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">µ-array<xref ref-type="table-fn" rid="nt102">b</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt103">c</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt104">d</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt105">e</xref>
</td><td align="left" rowspan="1" colspan="1">qPCR<xref ref-type="table-fn" rid="nt106">f</xref>
<sup>.</sup>
<xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">µ-array</td><td align="left" rowspan="1" colspan="1">qPCR</td><td align="left" rowspan="1" colspan="1">µ-array</td><td align="left" rowspan="1" colspan="1">qPCR</td><td align="left" rowspan="1" colspan="1">µ-array</td><td align="left" rowspan="1" colspan="1">µ-array</td><td align="left" rowspan="1" colspan="1">qPCR</td><td align="left" rowspan="1" colspan="1">qPCR</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>chr1.1224.1</italic>
</td><td align="left" rowspan="1" colspan="1">2.5</td><td align="left" rowspan="1" colspan="1">20<xref ref-type="table-fn" rid="nt108">h</xref>
</td><td align="left" rowspan="1" colspan="1">1.7</td><td align="left" rowspan="1" colspan="1">2.0</td><td align="left" rowspan="1" colspan="1">2.9</td><td align="left" rowspan="1" colspan="1">3.0</td><td align="left" rowspan="1" colspan="1">2.6</td><td align="left" rowspan="1" colspan="1">4.1</td><td align="left" rowspan="1" colspan="1">2.9</td><td align="left" rowspan="1" colspan="1">2.7</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>(Lisch-like)</italic>
</td><td align="left" rowspan="1" colspan="1">9×10<sup>−4</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">5×10<sup>−7</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2×10<sup>−3</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4×10<sup>−4</sup>
</td><td align="left" rowspan="1" colspan="1">7×10<sup>−4</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Lisch-like antisense</italic></bold>
</td><td align="left" rowspan="1" colspan="1">0.5</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.3</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NE<sup>i</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">6×10<sup>−3</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3×10<sup>−9</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Tada1l</italic></bold>
</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1.4</td><td align="left" rowspan="1" colspan="1">1.4</td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">1.05</td><td align="left" rowspan="1" colspan="1">1.1</td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">0.5</td><td align="left" rowspan="1" colspan="1">1.1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS<xref ref-type="table-fn" rid="nt110">j</xref>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1×10<sup>−8</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Pogk</italic></bold>
</td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">1.3</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0.73</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.2</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">0.7</td><td align="left" rowspan="1" colspan="1">1.05</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1×10<sup>−3</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>FMO13</italic>
</td><td colspan="10" align="left" rowspan="1">Not on array - <italic>inclusive-only</italic>
<xref ref-type="table-fn" rid="nt111">k</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt112">l</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>FMO12</italic>
</td><td colspan="10" align="left" rowspan="1">Not on array - <italic>inclusive-only</italic>
<xref ref-type="table-fn" rid="nt111">k</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt112">l</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>FMO9</italic></bold>
</td><td colspan="10" align="left" rowspan="1">Not on array - <italic>inclusive-only</italic>
<xref ref-type="table-fn" rid="nt111">k</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt112">l</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>C030014K22</italic></bold>
</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">1.1</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">1.4</td><td align="left" rowspan="1" colspan="1">1.1</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Uck2</italic></bold>
</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">1.4</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">1.4</td><td align="left" rowspan="1" colspan="1">1.2</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.3</td><td align="left" rowspan="1" colspan="1">0.7</td><td align="left" rowspan="1" colspan="1">0.8</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2×10<sup>−2</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Tmco1</italic></bold>
</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0.7</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">0.7</td><td align="left" rowspan="1" colspan="1">1.5</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2×10<sup>−5</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3×10<sup>−4</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2×10<sup>−4</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2×10<sup>−2</sup>
</td><td align="left" rowspan="1" colspan="1">3×10<sup>−2</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Aldh9a1</italic></bold>
</td><td align="left" rowspan="1" colspan="1">1.4</td><td align="left" rowspan="1" colspan="1">1.3</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">1.4</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">1.6</td><td align="left" rowspan="1" colspan="1">1.7</td><td align="left" rowspan="1" colspan="1">1.3</td><td align="left" rowspan="1" colspan="1">1.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1×10<sup>−6</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3×10<sup>−7</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4×10<sup>−5</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3×10<sup>−4</sup>
</td><td align="left" rowspan="1" colspan="1">3×10<sup>−2</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Mgst3</italic></bold>
</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">1.3</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0.7</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">0.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1×10<sup>−9</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2×10<sup>−8</sup>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">7×10<sup>−6</sup>
</td><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Lrrc52</italic></bold>
</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Rxrg</italic></bold>
</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">0.7</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">1.6</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.2</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">1.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">NS</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>Lmx1a</italic></bold>
</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">nd</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">2.0</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">NE</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.7</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula></disp-formula> <disp-formula></disp-formula>
<supplementary-material content-type="local-data" id="pgen.1000137.s001"><label>Table S1</label><caption><p>Data for ##FIG##2##Figure 3A##: Plasma Insulin/Glucose Ratios in Age-Grouped 1jc <italic>Lep<sup>ob/ob</sup></italic> Males.</p><p>(0.03 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000137.s002"><label>Table S2</label><caption><p>Pair-Wise Similarity Scores by Isoform and Domain for ##FIG##8##Figure 9##: ClustalW Analysis of <italic>Lisch-like</italic> Homologs and the <italic>LSR</italic> Protein.</p><p>(0.03 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000137.s003"><label>Table S3</label><caption><p>ID numbers for Affymetrix MOE-430A Probes used in <xref ref-type="sec" rid="s4">Methods</xref>: Microarray Gene Expression Analysis.</p><p>(0.03 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000137.s004"><label>Table S4</label><caption><p>Oligonucleotide Sequences.</p><p>(0.05 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000137.s005"><label>Table S5</label><caption><p>Data for ##FIG##11##Figure 12##: Terminal Phenotypes in 404 Obese F2 Progeny by Genotype at D1mit110.</p><p>(0.04 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>a</label><p>Boldface indicates RefSeq gene; regular type is locally-confirmed, predicted transcript.</p></fn><fn id="nt102"><label>b</label><p>Expression profiling was carried out using Affymetrix #430A mouse gene chips.</p></fn><fn id="nt103"><label>c</label><p>Ten 21-day old DD and BB <italic>Lep<sup>ob/ob</sup></italic> 1jc congenic males were analyzed.</p></fn><fn id="nt104"><label>d</label><p>Ratio is average B/B signal divided by average D/D signal in the organ.</p></fn><fn id="nt105"><label>e</label><p>Lower line in µ-array box is p-value, 2-sided t-test, comparing each set of 10 BB and 10 DD mice in the same organ.</p></fn><fn id="nt106"><label>f</label><p>Samples for qPCR prepared from 5 BB and 5 DD 90-day old <italic>Lep<sup>ob/ob</sup></italic> 1jc males on 2 occasions.</p></fn><fn id="nt107"><label>g</label><p>Ratios represent relative numbers of transcripts (BB/DD).</p></fn><fn id="nt108"><label>h</label><p>Primer-pairs amplify transcripts of <italic>Ll</italic> isoforms 1, 2, 4, and 5; these comprise &gt;90% of <italic>Ll</italic> transcripts.</p></fn><fn id="nt109"><label>i</label><p>Expression not detected.</p></fn><fn id="nt110"><label>j</label><p>Not significant (p&gt;0.05).</p></fn><fn id="nt111"><label>k</label><p>
<italic>Inclusive-only</italic> transcripts were detected in a cDNA pool that included whole embryos, 1-day old pups, and other tissues, but not in the cDNA pool prepared from diabetes-relevant organs.</p></fn><fn id="nt112"><label>l</label><p>Probes for these genes were neither on the Affymetrix #430A nor analyzed by qPCR.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This work supported by NIH grant RO1-DK66518 (RL), DERC 5P30DK063608, NYORC 5T32DK007599-18. The Berrie Foundation Program in Cellular Therapies of Diabetes, Albert Renold Fellowship of the European Foundation for the Study of Diabetes (EFSD) and grants from the Eli Lilly Company and GlaxoSmithKline.</p></fn></fn-group>
<graphic xlink:href="pgen.1000137.g001"/> <graphic xlink:href="pgen.1000137.g002"/> <graphic xlink:href="pgen.1000137.g003"/> <graphic xlink:href="pgen.1000137.g004"/> <graphic xlink:href="pgen.1000137.g005"/> <graphic xlink:href="pgen.1000137.g006"/> <graphic xlink:href="pgen.1000137.g007"/> <graphic id="pgen-1000137-t001-1" xlink:href="pgen.1000137.t001"/> <graphic xlink:href="pgen.1000137.g008"/> <graphic xlink:href="pgen.1000137.g009"/> <graphic xlink:href="pgen.1000137.g010"/> <graphic xlink:href="pgen.1000137.g011"/> <graphic xlink:href="pgen.1000137.g012"/> <graphic xlink:href="pgen.1000137.e001.jpg" mimetype="image" position="float"/> <graphic xlink:href="pgen.1000137.e002.jpg" mimetype="image" position="float"/>
<media xlink:href="pgen.1000137.s001.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000137.s002.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000137.s003.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000137.s004.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000137.s005.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["60"], "element-citation": ["\n"], "surname": ["Dayhoff", "Dayhoff"], "given-names": ["M", "M"], "year": ["1978"], "article-title": ["Atlas of Protein Sequence and Structure."], "publisher-loc": ["Washington, DC"], "publisher-name": ["National Biochemical Research Foundation"], "fpage": ["353"], "lpage": ["358"]}, {"label": ["64"], "element-citation": ["\n"], "surname": ["Schulz"], "given-names": ["H"], "year": ["2003"], "article-title": ["Towards a Comprehensive Description of the Human Retinal Transcriptome: Identification and Characterization of Differentially Expressed Genes [PhD dissertation]: University of Wurzberg"], "fpage": ["5"]}, {"label": ["102"], "element-citation": ["\n"], "surname": ["St. Jean", "Hsueh", "Mitchell", "Ehm", "Wagner", "Burns", "Shuldiner"], "given-names": ["P", "W-C", "B", "M", "M", "D", "AR"], "year": ["2000"], "article-title": ["Association between diabetes, obesity, glucose and insulin levels in the Old Amish and SNP's on 1q21-23 (Abstract)."], "source": ["American Journal of Human Genetics"], "volume": ["67"], "fpage": ["1848"]}, {"label": ["108"], "element-citation": ["\n"], "surname": ["McCarthy", "Shuldiner", "Bogardus", "Hanson", "Elbein"], "given-names": ["M", "AR", "C", "RL", "S"], "year": ["2004"], "article-title": ["Positional Cloning of a Type 2 Diabetes Susceptibility Gene on Chromosome 1q: A collaborative effort by the Chromosome 1q Diabetes Positional Cloning Consortium"], "fpage": ["1"], "lpage": ["39"]}, {"label": ["122"], "element-citation": ["\n"], "surname": ["Westerfield"], "given-names": ["M"], "year": ["2000"], "source": ["The Zebrafish Book: A Guide for the Laboratory Use of Zebrafish (Danio rerio)"], "publisher-loc": ["Eugene, OR"], "publisher-name": ["University of Oregon Press"]}, {"label": ["123"], "element-citation": ["\n"], "surname": ["Thisse", "Thisse"], "given-names": ["C", "B"], "year": ["1998"], "article-title": ["High resolution whole-mount in situ hybridization."], "source": ["Zebrafish Sci Monit"], "volume": ["5"], "fpage": ["8"], "lpage": ["9"]}]
PLoS Genet. 2008 Jul 25; 4(7):e1000137