Both isoforms of ketohexokinase are dispensable for normal growth and development.

Dietary fructose intake has dramatically increased over recent decades and is implicated in the high rates of obesity, hypertension, and type 2 diabetes (metabolic syndrome) in Western societies. The molecular determinants of this epidemiologic correlation are incompletely defined, but high-flux fructose catabolism initiated by ketohexokinase (Khk, fructokinase) is believed to be important. The Khk gene encodes two enzyme isoforms with distinctive substrate preferences, the independent physiological roles of which are unclear. To investigate this question, and for testing the importance of Khk in metabolic syndrome, isoform-selective genetic lesions would be valuable. Two deficiency alleles of the mouse Khk gene were designed. The first, Khk(3a), uses targeted "knock-in" of a premature termination codon to induce a selective deficiency of the minor Khk-A isoform, preserving the major Khk-C isoform. The second, the Khk(Δ) allele, ablates both isoforms. Mice carrying each of these Khk-deficiency alleles were generated and validated at the DNA, RNA, and protein levels. Comparison between normal and knockout animals confirmed the specificity of the genetic lesions and allowed accurate analysis of the cellular distribution of Khk within tissues such as gut and liver. Both Khk(3a/3a) and Khk(Δ/Δ) homozygous mice were healthy and fertile and displayed minimal biochemical abnormalities under basal dietary conditions. These studies are the first demonstration that neither Khk isoform is required for normal growth and development. The new mouse models will allow direct testing of various hypotheses concerning the role of this enzyme in metabolic syndrome in humans and the value of Khk as a pharmacological target.

Inferring relative proportions of DNA variants from sequencing electropherograms.

MOTIVATION: Determination of the relative copy number of single-nucleotide sequence variants (SNVs) within a DNA sample is a frequent experimental goal. Various methods can be applied to this problem, although hybridization-based approaches tend to suffer from high-setup cost and poor adaptability, while others (such as pyrosequencing) may not be accessible to all laboratories. The potential to extract relative copy number information from standard dye-terminator electropherograms has been little explored, yet this technology is cheap and widely accessible. Since several biologically important loci have paralogous copies that interfere with genotyping, and which may also display copy number variation (CNV), there are many situations in which determination of the relative copy number of SNVs is desirable. RESULTS: We have developed a desktop application, QSVanalyzer, which allows high-throughput quantification of the proportions of DNA sequences containing SNVs. In reconstruction experiments, QSVanalyzer accurately estimated the known relative proportions of SNVs. By analyzing a large panel of genomic DNA samples, we demonstrate the ability of the software to analyze not only common biallelic SNVs, but also SNVs within a locus at which gene conversion between four genomic paralogs operates, and within another that is subject to CNV. AVAILABILITY AND IMPLEMENTATION: QSVanalyzer is freely available at http://dna.leeds.ac.uk/qsv/. It requires the Microsoft .NET framework version 2.0, which can be installed on all Microsoft operating systems from Windows 98 onwards. CONTACT: msjimc@leeds.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Phenotype associated with recessively inherited mutations in DNA mismatch repair (MMR) genes.

The MMR (DNA mismatch repair) system helps to maintain the integrity of the genome. This involves eliminating base-base mismatches and insertion/deletion loops, which can lead to microsatellite instability, as seen in tumour cells. Hereditary non-polyposis colon cancer is the result of dominant mutations in MMR genes, such as MLH1, MSH2 and MSH6. More recently there have been case reports of biallelic mutations in the MMR genes MLH1, MSH2 and PMS2. These result in a distinct autosomal recessive cancer predisposition syndrome. The syndrome is characterized by childhood haematological malignancies, brain tumours and the presence of café au lait patches. Second primaries occur frequently in this condition, and survival into adulthood is rare.

Mutations in the glucokinase regulatory protein gene in 2p23 in obese French caucasians.

AIMS/HYPOTHESIS: Glucokinase regulatory protein (GKRP) controls the activity of glucokinase in liver but possibly also in some areas of the central nervous system, suggesting that it could play a role in body mass control. Its gene is located in a region (2p21-23) linked to serum leptin levels. Our goal was to investigate whether mutations in the GKRP gene were associated with obesity. METHODS: Mutations were sought in the GKRP gene of 57 patients from the families of the French genome-wide scan for obesity that contributed most to the positive LOD score with 2p21-23. The identified mutations were further sought in 720 unrelated obese individuals and 384 individuals of normal weight and their effect on the properties of recombinant GKRP were investigated. RESULTS: The most frequent mutation (Pro446Leu) had a similar allele frequency in the obese (0.63) and normal weight (0.64) subjects and did not affect the properties of GKRP. Similarly, no effect on the properties of GKRP was observed with Arg590Tyr, found in 10 out of 720 obese subjects and in 2 out of 384 control subjects (p=0.18). Mutation Arg227Stop was found in one obese family and in 1 out of 384 control subjects and led to an insoluble protein. Mutation Arg518Gln, replacing a conserved residue, led to a marked decrease in the affinity of GKRP for both fructose 6-phosphate and fructose 1-phosphate and to a destabilization of GKRP. However, this mutation did not co-segregate with obesity in the single family in which it was found. CONCLUSIONS/INTERPRETATION: Mutations that affect the properties of GKRP are found in the French population, but they do not seem to account for the linkage between the 2p23 locus and quantitative markers of obesity.

Lack of involvement of known DNA methyltransferases in familial hydatidiform mole implies the involvement of other factors in establishment of imprinting in the human female germline.

BACKGROUND: Differential methylation of the two alleles is a hallmark of imprinted genes. Correspondingly, loss of DNA methyltransferase function results in aberrant imprinting and abnormal post-fertilization development. In the mouse, mutations of the oocyte-specific isoform of the DNA methyltransferase Dnmt1 (Dnmt1o) and of the methyltransferase-like Dnmt3L gene result in specific failures of imprint establishment or maintenance, at multiple loci. We have previously shown in humans that an analogous inherited failure to establish imprinting at multiple loci in the female germline underlies a rare phenotype of recurrent hydatidiform mole. RESULTS: We have identified a human homologue of the murine Dnmt1o and assessed its pattern of expression. Human DNMT1o mRNA is detectable in mature oocytes and early fertilized embryos but not in any somatic tissues analysed. The somatic isoform of DNMT1 mRNA, in contrast, is not detectable in human oocytes. In the previously-described family with multi-locus imprinting failure, mutation of DNMT1o and of the other known members of this gene family has been excluded. CONCLUSIONS: Mutation of the known DNMT genes does not underlie familial hydatidiform mole, at least in the family under study. This suggests that trans-acting factors other than the known methyltransferases are required for imprint establishment in humans, a concept that has indirect support from recent biochemical studies of DNMT3L.

Imprinting of the G(s)alpha gene GNAS1 in the pathogenesis of acromegaly.

Approximately 40% of growth hormone-secreting pituitary adenomas have somatic mutations in the GNAS1 gene (the so-called gsp oncogene). These mutations at codon 201 or codon 227 constitutively activate the alpha subunit of the adenylate cyclase-stimulating G protein G(s). GNAS1 is subject to a complex pattern of genomic imprinting, its various promoters directing the production of maternally, paternally, and biallelically derived gene products. Transcripts encoding G(s)alpha are biallelically derived in most human tissues. Despite this, we show here that in 21 out of 22 gsp-positive somatotroph adenomas, the mutation had occurred on the maternal allele. To investigate the reason for this allelic bias, we also analyzed GNAS1 imprinting in the normal adult pituitary and found that G(s)alpha is monoallelically expressed from the maternal allele in this tissue. We further show that this monoallelic expression of G(s)alpha is frequently relaxed in somatotroph tumors, both in those that have gsp mutations and in those that do not. These findings imply a possible role for loss of G(s)alpha imprinting during pituitary somatotroph tumorigenesis and also suggest that G(s)alpha imprinting is regulated separately from that of the other GNAS1 products, NESP55 and XLalphas, imprinting of which is retained in these tumors.

Expression, purification and preliminary crystallographic studies of human ketohexokinase.

Ketohexokinase (KHK; E.C. 2.7.1.3) catalyses the (reversible) phosphorylation of fructose to fructose-1-phosphate. KHK is the first enzyme in a specialized catabolic pathway metabolizing dietary fructose to the glycolytic intermediate glyceraldehyde-3-phosphate. Mutations inactivating KHK underlie the metabolic disorder essential fructosuria. The primary structure of KHK shows no significant homology to other mammalian hexokinases. It is most similar to prokaryotic ribokinases, but catalyses a distinct phosphorylation reaction. Recombinant human KHK has been crystallized in the orthorhombic form (space group P2(1)2(1)2 or P2(1)2(1)2(1)). Single crystals of this polymorph suitable for X-ray diffraction have been obtained by vapour diffusion using 2-propanol and MPD as precipitants (pH 7.5). The crystals have unit-cell parameters a = 93.4, b = 121.5, c = 108.4 A. Diffraction data were collected to 4.3 A resolution. The asymmetric unit contains four protein molecules.

Tissue-specific expression of antisense and sense transcripts at the imprinted Gnas locus.

The mouse Gnas gene encodes an important signal transduction protein, the alpha subunit of the stimulatory G protein, G(s). In humans, partial deficiency of G(s)alpha, the alpha subunit of G(s), results in the hormone-resistance syndrome pseudohypoparathyroidism type 1a. The mouse Gnas (and the human GNAS1) locus is transcribed from three promoter regions. Transcripts from P1, which encode Nesp55, are derived from the maternal allele only. Transcripts from P2 encode Xlalphas and are derived only from the paternal allele, while transcripts from P3 encode the alpha subunit and are from both parental alleles. The close proximity of reciprocal imprinting suggests the presence of important putative imprinting elements in this region. In this report, we demonstrate that the reciprocal imprinting occurs in normal tissues of interspecific (Mus spretus x C57BL/6) mice. Transcripts from P1 are most abundant in CNS (pons and medulla) in contrast to the more ubiquitous expression from P2 and P3. In the P1-P2 genomic region, we have identified an antisense transcript that starts 2.2 kb upstream of the P2 exon and spans the P1 region. While the P1 transcript is derived from the maternal allele, the P1-antisense (Gnas-as) is derived only from the paternal allele in most but not all tissues. Although both the Nesp55 region and the Gnas-as transcripts are present in cerebral cortex, adrenal, and spleen, Gnas-as is abundant in some tissues in which transcription from the Nesp55 region is negligible. Furthermore, the Nesp55 region transcripts remain strictly imprinted in tissues that lack Gnas-as. Our results suggest that multiple imprinting elements, including the unique Gnas-as, regulate the allelic expression of the Nesp55 region sense transcript.

Characterization of TH1 and CTSZ, two non-imprinted genes downstream of GNAS1 in chromosome 20q13.

The clustering and coordinate regulation of many imprinted genes justifies positional searches for imprinted genes adjacent to known ones. We recently characterized a locus on 20q13, containing GNAS1, which has a highly complex imprinted expression pattern. In a search for neighbouring genes, we have now characterized a new gene, TH1, downstream of GNAS1. TH1 and GNAS1 are separated by more than 70 kb consisting largely of interspersed repetitive DNA. TH1 is the homologue of a gene that, in Drosophila, lies adjacent to the DNA repair gene mei-41. We have determined the full-length structures of human, mouse and Drosophila TH1. Though of unknown function, TH1 is highly conserved and widely expressed. Nonetheless, there is no similar Caenorhabditis elegans protein. We have also determined the complete genomic structures of human and Drosophila TH1. The Drosophila gene has five exons spanning 2.6 kb. The last three introns have precise equivalents in the human gene, which has 15 exons spanning 14 kb and is transcribed away from GNAS1. Using a single-nucleotide polymorphism in the 3' untranslated region, we have demonstrated biallelic TH1 expression in human fetal tissues, suggesting that, unlike GNAS1, TH1 is probably not imprinted. Immediately downstream of TH1 lies CTSZ, encoding the recently described cysteine protease, cathepsin Z. We have also elucidated the genomic structure of this gene; it has six exons spanning 12 kb and is oriented tail-to-tail with TH1, only 70 bp separating their polyadenylation sites. A polymorphism was again identified within the CTSZ 3' untranslated region and used to demonstrate biallelic expression in fetal tissues.

Structure and mutation analysis of the gene encoding DNA fragmentation factor 40 (caspase-activated nuclease), a candidate neuroblastoma tumour suppressor gene.

We have characterised the DFFB gene, encoding the active subunit of the apoptotic nuclease DNA fragmentation factor (DFF40). DFFB maps to 1p36, near the imprinted putative tumour suppressor gene TP73. The DFFA gene (encoding the inhibitory DFF45 subunit) also maps to 1p36.2-36.3, and we show by FISH that DFFB lies distal to DFFA. We have also mapped a processed DFFB pseudogene to chromosome 9. DFFB itself has seven coding exons spanning 10 kb. Exhaustive mutation screening of 41 neuroblastomas and other tumours in which a 1p36 tumour suppressor gene is implicated showed no tumour-specific mutations. A coding region polymorphism was used to demonstrate uniformly biallelic expression in human fetal DFFB transcripts. Since the putative neuroblastoma tumour suppressor gene in distal 1p36 is predicted to be maternally expressed, the lack of imprinting and absence of somatic mutations in DFFB indicate that it is probably not the neuroblastoma tumour suppressor gene.

An imprinted antisense transcript at the human GNAS1 locus.

Recent studies of the GNAS1 gene have shown a highly complex imprinted expression pattern, with paternally, maternally and biallelically derived protein products, raising questions regarding how such transcriptional complexity is established and maintained. GNAS1 was originally identified as the gene encoding an important and widely expressed signal transduction protein, the alpha subunit of the stimulatory G protein G(s). Partial G(s)alpha deficiency results in the hormone resistance syndrome, pseudohypoparathyroidism type 1a. G(s)alpha is encoded by exons 1-13 of GNAS1 and, in most tissues at least, expression of this transcript is biallelic. Two large upstream exons, however, have monoallelic expression patterns, and in each case their transcripts splice onto GNAS1 exon 2. The most 5' of these is maternally expressed, and encodes neuroendocrine secretory protein 55 (NESP55), whose coding region does not overlap with that of G(s)alpha. The other exon, 14 kb further 3', is paternally expressed, and encodes XL(alpha)s (extra large alphas-like protein), translated in-frame with G(s)alpha exons 2-13. This close proximity of two oppositely imprinted promoters suggested the likelihood of important regulatory interactions between them, and to investigate this possibility we have performed a search for other transcripts in the region. Here we show that the maternally methylated region upstream of the XL(alpha)s exon gives rise to a spliced polyadenylated antisense transcript, which spans the upstream NESP55 region. This antisense transcript is imprinted, and expressed only from the paternal allele, suggesting that it may have a specific role in suppressing in cis the activity of the paternal NESP55 allele.

The cell cycle control gene ZAC/PLAGL1 is imprinted--a strong candidate gene for transient neonatal diabetes.

We describe a screen for new imprinted human genes, and the identification in this way of ZAC (zinc finger protein which regulates apoptosis and cell cycle arrest)/ PLAGL1 (pleomorphicadenoma of the salivary gland gene like 1) as a strong candidate gene for transient neonatal diabetes mellitus (TNDM). To screen for imprinted genes, we compared parthenogenetic DNA from the chimeric patient FD and androgenetic DNA from hydatidiform mole, using restriction landmark genome scanning for methylation. This resulted in identification of two novel imprinted loci, one of which (NV149) we mapped to the TNDM region of 6q24. From analysis of the corresponding genomic region, it was determined that NV149 lies approximately 60 kb upstream of the ZAC / PLAGL1 gene. RT-PCR analysis was used to confirm that this ZAC / PLAGL1 is expressed only from the paternal allele in a variety of tissues. TNDM is known to result from upregulation of a paternally expressed gene on chromosome 6q24. The paternal expression, map position and known biological properties of ZAC / PLAGL1 make it highly likely that it is the TNDM gene. In particular, ZAC / PLAGL1 is a transcriptional regulator of the type 1 receptor for pituitary adenylate cyclase-activating polypeptide, which is the most potent known insulin secretagog and an important mediator of autocrine control of insulin secretion in the pancreatic islet.

Complex patterns of intragenic polymorphism at the PDGFA locus.

The human platelet-derived growth factor A chain gene (PDGFA) on chromosome 7p22 encodes an important mitogen. Within PDGFA lies a complex minisatellite structure that results in partial duplications of exon 4 and the IVS4 splice donor site. Here, we show that the PDGFA genes of four ape species and an Old-World monkey all have similar complex minisatellites at this position. Comparison of their structures suggests evolutionary constraints resulting from the protein-coding function of the minisatellite. Nonetheless, the IVS4 minisatellite seems to have undergone independent expansion events in different primate lineages. Within the human IVS4 minisatellite, an embedded pentanucleotide repeat, based on the sequence (CCTCC)n, shows frequent subunit sequence variation but only rare length polymorphism. In contrast, within IVS3 of human PDGFA, we have discovered a second minisatellite which, unlike the IVS4 minisatellite, is highly polymorphic. The subunit sequences of these two minisatellites, which lie less than 0.5 kb apart, are non-identical, but share a CnT-rich core. Two new single nucleotide polymorphisms (SNPs), in exon 3 and IVS4, are in linkage disequilibrium, despite flanking the two minisatellite regions. Reverse transcription-polymerase chain reaction analysis of the exon 3 SNP in human foetal tissues demonstrated biallelic expression of PDGFA in all tissues examined. The unusual location of PDGFA exon 4 between two minisatellite sequences, together with its partial duplication, may have functional implications, particularly for the splicing of the gene. The high level of polymorphism demonstrated in this region will also be valuable for disease-association and linkage studies of the PDGFA locus.

A locus for isolated cleft palate, located on human chromosome 2q32.

We present evidence for the existence of a novel chromosome 2q32 locus involved in the pathogenesis of isolated cleft palate. We have studied two unrelated patients with strikingly similar clinical features, in whom there are apparently balanced, de novo cytogenetic rearrangements involving the same region of chromosome 2q. Both children have cleft palate, facial dysmorphism, and mild learning disability. Their karyotypes were originally reported as 46, XX, t(2;7)(q33;p21) and 46, XX, t(2;11)(q33;p14). However, our molecular cytogenetic analyses localize both translocation breakpoints to a small region between markers D2S311 and D2S116. This suggests that the true location of these breakpoints is 2q32 rather than 2q33. To obtain independent support for the existence of a cleft-palate locus in 2q32, we performed a detailed statistical analysis for all cases in the human cytogenetics database of nonmosaic, single, contiguous autosomal deletions associated with orofacial clefting. This revealed 2q32 to be one of only three chromosomal regions in which haploinsufficiency is significantly associated with isolated cleft palate. In combination, our data provide strong evidence for the location at 2q32 of a gene that is critical to the development of the secondary palate. The close proximity of these two translocation breakpoints should also allow rapid progress toward the positional cloning of this cleft-palate gene.

Bidirectional imprinting of a single gene: GNAS1 encodes maternally, paternally, and biallelically derived proteins.

The GNAS1 gene encodes the alpha subunit of the guanine nucleotide-binding protein Gs, which couples signaling through peptide hormone receptors to cAMP generation. GNAS1 mutations underlie the hormone resistance syndrome pseudohypoparathyroidism type Ia (PHP-Ia), so the maternal inheritance displayed by PHP-Ia has raised suspicions that GNAS1 is imprinted. Despite this suggestion, in most tissues Gsalpha is biallelically encoded. In contrast, the large G protein XLalphas, also encoded by GNAS1, is paternally derived. Because the inheritance of PHP-Ia predicts the existence of maternally, rather than paternally, expressed transcripts, we have investigated the allelic origin of other mRNAs derived from GNAS1. We find this gene to be remarkable in the complexity of its allele-specific regulation. Two upstream promoters, each associated with a large coding exon, lie only 11 kb apart, yet show opposite patterns of allele-specific methylation and monoallelic transcription. The more 5' of these exons encodes the neuroendocrine secretory protein NESP55, which is expressed exclusively from the maternal allele. The NESP55 exon is 11 kb 5' to the paternally expressed XLalphas exon. The transcripts from these two promoters both splice onto GNAS1 exon 2, yet share no coding sequences. Despite their structural unrelatedness, the encoded proteins, of opposite allelic origin, both have been implicated in regulated secretion in neuroendocrine tissues. Remarkably, maternally (NESP55), paternally (XLalphas), and biallelically (Gsalpha) derived proteins all are produced by different patterns of promoter use and alternative splicing of GNAS1, a gene showing simultaneous imprinting in both the paternal and maternal directions.

Structure and alternative splicing of the ketohexokinase gene.

Ketohexokinase (fructokinase, KHK) catalyses the phosphorylation of fructose to fructose-l-phosphate. It thereby initiates the intracellular catabolism of a large proportion of dietary carbohydrate. Although found at high level in liver, renal cortex and small intestine, fructokinase activity has also been known for many years to be present at lower levels in most other tissues. We previously found that there appeared to be two isoforms of human KHK, and have now investigated the molecular basis for this in human, rat and mouse. Cloning of the human KHK gene, on chromosome 2p23.2-2p23.3, shows that it has nine exons, spanning 14 kb. An intragenic duplication has resulted in two similar 135-bp exons (designated 3a and 3c), separated by a short intron. Exon 3a and exon 3c are mutually exclusively spliced into KHK mRNA. This exon-intron structure and the pattern of alternative splicing are conserved in both the rat and mouse, suggesting distinct conserved functions for the two KHK isoforms. The alternative splicing is also tissue specific, since in both rat and human, tissues expressing high levels of KHK (liver, kidney and duodenum) utilise exclusively the 3c exon, while other tissues use only 3a. Furthermore, comparison of human foetal and adult tissues indicates a developmental splicing shift from use of exon 3a to exon 3c.

The human GNAS1 gene is imprinted and encodes distinct paternally and biallelically expressed G proteins.

The GNAS1 gene encodes the alpha subunit of the G protein Gs, which couples receptor binding by several hormones to activation of adenylate cyclase. Null mutations of GNAS1 cause pseudohypoparathyroidism (PHP) type Ia, in which hormone resistance occurs in association with a characteristic osteodystrophy. The observation that PHP Ia almost always is inherited maternally has led to the suggestion that GNAS1 may be an imprinted gene. Here, we show that, although Gsalpha expression (directed by the promoter upstream of exon 1) is biallelic, GNAS1 is indeed imprinted in a promoter-specific fashion. We used parthenogenetic lymphocyte DNA to screen by restriction landmark genomic scanning for loci showing differential methylation between paternal and maternal alleles. This screen identified a region that was found to be methylated exclusively on a maternal allele and was located approximately 35 kb upstream of GNAS1 exon 1. This region contains three novel exons that are spliced into alternative GNAS1 mRNA species, including one exon that encodes the human homologue of the large G protein XLalphas. Transcription of these novel mRNAs is exclusively from the paternal allele in all tissues examined. The differential imprinting of separate protein products of GNAS1 therefore may contribute to the anomalous inheritance of PHP Ia.