Article Text

PDF

Identification of a novel mutation disrupting the DNA binding activity of GCM2 in autosomal recessive familial isolated hypoparathyroidism
  1. L Baumber1,
  2. C Tufarelli1,
  3. S Patel2,
  4. P King3,
  5. C A Johnson4,
  6. E R Maher4,
  7. R C Trembath1
  1. 1Division of Medical Genetics, Departments of Genetics and Cardiovascular Sciences, University of Leicester, Leicester, UK
  2. 2Leicestershire Genetics Centre, University Hospitals of Leicester NHS Trust, Leicester, UK
  3. 3Derby Royal Infirmary, Derby, UK
  4. 4Section of Medical and Molecular Genetics, Department of Paediatrics and Child Health, University of Birmingham Medical School, Birmingham, UK
  1. Correspondence to:
 Richard C Trembath
 Professor of Medical Genetics, Division of Medical Genetics, Departments of Genetics and Cardiovascular Sciences, University of Leicester, University Road, Leicester LE1 7RH, UK; rtrembathgmp.mrc.ac.uk

Statistics from Altmetric.com

Hypoparathyroidism is a heterogeneous group of disorders with both acquired and inherited causes, each presenting clinically with hypocalcaemia. Familial cases of hypoparathyroidism may be due to an isolated defect of the parathyroid glands or be a component of a syndrome disorder, examples of which include DiGeorge, hypoparathyroidism-retardation-dysmorphism, and Kenny-Caffey syndrome.1 Familial isolated hypoparathyroidism (FIH) is characterised by hypocalcaemia and hyperphosphataemia and may be due to an inherited deficiency or the abnormal activity of parathyroid hormone (PTH). FIH is heterogeneous with X linked, autosomal dominant, and autosomal recessive modes of inheritance reported.2 Mutations in the PTH gene on chromosome 11p have been described in both autosomal dominant and autosomal recessive forms of the disorder.3,4 Mature PTH is generated by cleavage of the signal peptide from the prepro-PTH, followed by successive proteolysis of pro-PTH within the Golgi apparatus.5 A requirement for normal PTH processing is emphasised by the detection of FIH disease causing mutations, which include two missense amino acid substitutions within the signal region of the prepro-PTH gene, leading to inefficient cleavage by signal peptidases.3,5

Heterozygous defects in the calcium sensing receptor (CASR) have been detected in autosomal dominant FIH and sporadic cases of hypoparathyroidism.6–8 All known CASR mutations in FIH lead to receptor activation.9,10 The CASR normally functions to provide a negative feedback loop, whereby binding of extracellular Ca2+ activates the receptor, resulting in decreased PTH secretion by parathyroid cells, but also through inhibition of calcium reabsorption in the renal distal tubule. Constitutive activation of this receptor therefore results in hypocalcaemia and hypercalciuria.

Notwithstanding these reports, the molecular genetic basis of most cases of autosomal recessive FIH remains unknown.1 Recently, a homozygous deletion within the human GCM2 (glial cells missing, Drosophila homologue B, previously referred to as GCMB) gene has been shown to underlie hypoparathyroidism in one patient from a kindred with no consanguinity.11 GCM2 belongs to a family of novel transcription factors characterised by an N-terminal DNA binding domain, known as the gcm motif.12,13 This motif, unique to the GCM family, spans approximately 150 amino acids, exhibiting approximately 92% identity between human and mouse homologues.14 Human GCM1 and GCM2, the two identified isoforms of GCM, display additional homology to their mouse counterparts outside this region with GCM2 showing 68% and 47% similarity overall to Gcm2 and GCM1, respectively.15 Despite limited characterisation of human GCM, studies on the expression pattern of mouse Gcm2 protein have shown it to be confined to the developing parathyroid gland16 and homozygous mutant mice lacking the Gcm2 gene demonstrate absence of parathyroid tissue and consequent hypoparathyroidism.17 Taken together, these studies provide compelling evidence for a critical role of GCM2 in normal parathyroid gland function; however the molecular basis and interactions required in this process remain unclear.

Key points

  • Hypoparathyroidism is a heterogeneous group of disorders with both acquired and inherited causes, each presenting clinically with hypocalcaemia.

  • Familial isolated hypoparathyroidism is characterised by hypocalcaemia and hyperphosphataemia and may be due to an inherited deficiency or the abnormal activity of parathyroid hormone.

  • Recently, a homozygous deletion within the human GCM2 gene has been shown to underlie hypoparathyroidism in one patient, while other compelling evidence indicates a critical role for GCM2 in normal parathyroid gland function.

  • Mutation studies revealed a novel homozygous missense mutation of the DNA binding domain of GCM2, functional analysis of which now provides an important tool for the further delineation of GCM2 bioactivity and molecular mechanisms necessary for parathyroid gland development.

Here we report a consanguineous Pakistani family with multiple individuals presenting in early life with autosomal recessive FIH. Mutation studies revealed a novel homozygous missense mutation of the DNA binding domain of GCM2, functional analysis of which now provides an important tool for the further delineation of GCM2 bioactivity and molecular mechanisms necessary for parathyroid gland development.

METHODS

Clinical data

Subject IV:2 (fig 1), the second child of first cousin parents, was diagnosed as a neonate with primary hypoparathyroidism after presenting with convulsions. Her older sister had died in infancy following seizures and retrospectively was assumed to be similarly affected. The proband had no evidence of immunodeficiency, renal abnormalities, autoimmune disease, or dysmorphism. Physical examination was unremarkable; in particular, she had no evidence of neuromuscular irritability. Renal, liver, and thyroid functions were all normal. Biochemical investigations confirmed hypocalcaemia with a calcium level of 1.98 mmol/l (normal range 2.13–2.58), phosphate 2.02 mmol/l (normal range, 1.30–2.20), and a low PTH level of 5 ng/l (normal range, 12–72), findings compatible with primary hypoparathyroidism. She was treated with alphacalcidol and subsequent growth and development were reported to be normal.

Figure 1

 (A) Pedigree structure of the family with genotyping showing linkage to the GCM2 locus. Transmission of the phenotype is compatible with an autosomal recessive mode of inheritance. Genotype data are shown for the region encompassing markers D6S940 to D6S1653, spanning a 4.6 Mb interval. (B) Sequence chromatogram illustrating the GCM2 R47L mutation. Affected individuals carry a homozygous G→T substitution at nucleotide 140 (indicated by an asterisk). Nucleotides are numbered according to the cDNA sequence, with the adenosine of the initiation codon designated +1. (C) Segregation of the mutation was tracked through the family by HhaI digestion of PCR products. Mutation carriers are detected by the presence of a 246 bp band due to loss of an HhaI site. Gel lanes correspond to the individuals numbered in the pedigree (lane 1, normal control).

The unaffected sister of subject IV:2 married a first cousin and subsequently had two children (V:1 and V:2), both diagnosed clinically with primary hypoparathyroidism. Individual V:1 was a full term normal delivery, weighing 3.2 kg. She presented at 6 weeks of age with focal seizures involving the left arm and leg associated with left sided facial twitching. Laboratory investigation revealed hypocalcaemia at 1.35 mmol/l and hyperphosphataemia at 3.23 mmol/l. Alkaline phosphatase was markedly elevated at 1805 U/l (normal range, 76–196). Serum calcium levels were maintained by calcium and vitamin D supplementation. At presentation, magnesium was 0.69 mmol/l (normal range, 0.70–1.00), PTH level was less than 10 ng/l, and 25-hydroxy vitamin D level measured 11 nmol/l (normal range, 15.0–100.0). Chromosome analysis showed a normal female karyotype and FISH studies for the DiGeorge/velo-cardio-facial syndrome region on chromosome 22 showed no deletion.

Her younger sister (subject V:2) was a full term normal delivery weighing 2.75 kg. Because of the family history, she was screened for hypocalcaemia shortly after birth and initial serum calcium was 1.97 mmol/l. PTH level was at the lower end of the normal range (10 ng/l). She was also treated with vitamin D and calcium supplementation.

Informed consent was obtained from all subjects or their parents and the study was approved by the relevant local research ethics committees.

Molecular genetic analysis

Genomic DNA was extracted from 10 ml of peripheral blood by standard techniques. Autozygosity mapping was performed using fluorescently labelled microsatellite markers with an average heterozygosity of 75%. Markers for each locus are listed in table 1. PCR amplifications were achieved using MJ Research (Waltham, MA) DNA Engines (10 μl reactions; 50 ng DNA, 1.5 mM MgCl2, 0.2 mM dNTPs (Pharmacia, St Albans, Hertfordshire), 0.3 μM primer, 0.5 U Taq polymerase (ABgene, Epsom, UK) in a 10×KCl PCR buffer). PCR products were electrophoresed through 6% polyacrylamide gels (FMC Bioproducts, Rockland, ME) on an ABI 377 DNA sequencer (Applied Biosystems, Foster City, CA) and analysed using GeneScan v3.0 and Genotyper v2.1 software (PE-ABI, Foster City, CA). Linkage to candidate loci was examined under the assumption that FIH was segregating as an autosomal recessive trait due to homozygous mutation inherited common by descent.

Table 1

 Microsatellite markers spanning the known FIH loci

Mutation analysis of GCM2

The complete protein coding region and intron/exon boundaries of the GCM2 gene were amplified from all available family members and one normal control individual by PCR in 40 μl reactions. Primer sequences for PCR products spanning the five coding exons are listed in table 2. PCR products were separated by electrophoresis through 2% LE agarose (Cambrex, East Rutherford, NJ) to ensure the presence of sufficient quantities for sequence analysis, and were purified using the QIAquick PCR purification kit (Qiagen, Valencia, CA). Purified PCR products were sequenced, using the same forward primers as for the PCR reactions, with the Applied Biosystems BigDye terminator kit on an ABI 377 sequencer.

Table 2

 Primer sequences for GCM2 gene amplification

The mutation identified in family 1, at nucleotide position 140 in exon 2, abolished an HhaI restriction endonuclease recognition site. We therefore independently confirmed segregation of the variant within the family by restriction analysis with HhaI (New England Biolabs, Beverly, MA), and tested a panel of chromosomes from demographically matched control individuals. An 1198 bp fragment containing exons 2 and 3 was amplified by PCR under identical conditions as those described above, using primers GCM2ex2F and GCM2ex3R. Each PCR product was digested with 10 U enzyme at 37°C for 90 min. After digestion, the DNA fragments were separated on a 3% LE agarose gel and visualised with ethidium bromide.

Preparation of GCM constructs

Full length human GCM2 cDNA in pBluescript was kindly donated by A Giangrande.14 Using this plasmid as template, the arginine to leucine (R47L) point mutation was introduced by QuickChange (Stratagene, La Jolla, CA) site directed mutagenesis, according to the manufacturer’s protocol. The entire GCM2 and R47L coding sequences were amplified by PCR using primers GCM2-5′BamHI and GCM2-3′XhoI (primer sequences available on request). Products were digested with BamHI and XhoI and inserted into the same sites of the eukaryotic expression vector pCMV-Tag3B (Stratagene) to facilitate the incorporation of a 5′ c-Myc epitope. All constructs were verified by automated DNA sequencing.

Transfection and immunoblotting

COS cells were maintained at 37°C and 5% CO2 in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal calf serum and antibiotics. Cells were seeded onto 10 cm dishes and transiently transfected with 8 μg of expression plasmid using FuGene 6 reagent (Roche Molecular Biochemicals, Indianapolis, IN) according to the manufacturer’s recommendation. COS cells were harvested after 48 h and total cell extracts were prepared as described.18 Protein concentrations were quantified using the Bio-Rad (Hercules, CA) protein assay. Samples were separated by SDS-PAGE (10% acrylamide), transferred to nitrocellulose membrane, and stained with Ponceau S to visualise protein bands. After immunoblot analysis, the membrane was washed with PBS-Tween wash buffer (80 mM Na2HPO4, 20 mM KH2PO4, 100 mM NaCl, 0.1% Tween-20), supplemented with 5% non-fat dry milk (Marvel) for both the blocking solution and antibody diluent. Primary c-Myc antibody (Zymed, San Francisco, CA) was used at a 1:500 dilution. Horseradish peroxidase conjugated secondary antibody (1:2000 dilution; Sigma, St Louis, MO) was used to visualise bound primary antibody with ECL Plus Western Blotting Detection Reagent (Amersham Biosciences, Little Chalfont, UK).

Electrophoretic mobility shift assays

Sequences of DNA probes containing consensus (wt) and mutant (m3) GCM recognition sites were as described in Schreiber et al.18 Double stranded probes were generated by annealing of complementary oligonucleotides and labelled with γ-dATP 32P using T4 polynucleotide kinase (New England Biolabs).

Electrophoretic mobility shift assays (EMSA) were performed as described18 using 3 μg whole cell extracts from transiently transfected COS cells. Protein samples were incubated with 0.5 ng of 32P-labelled double stranded probe for 30 min at room temperature. Samples were loaded onto native 4% polyacrylamide gels (Flowgen, Nottingham, UK) and electrophoresed in 0.5×TBE (45 mM Tris, 45 mM boric acid, 1 mM EDTA, pH 8.3) at 180 V, 4°C for 75 min. Competition assays were performed using unlabelled double stranded wildtype and competitor (m3) probes in 1-, 10-, 20-, and 50-fold molar excess.

RESULTS

Autozygosity mapping

All available family members were typed with polymorphic microsatellite markers spanning the PTH locus over a 4.3 Mb interval. There was no evidence of linkage to PTH, with affected individuals displaying heterozygous genotypes for all markers across this region (data not shown).

Individuals were also typed for polymorphic markers flanking the GCM2 gene. The finding of homozygosity in affected individuals for all or some of these markers was considered compatible with linkage to GCM2 (fig 1A).

Mutation analysis of GCM2

We next sequenced all five coding exons and intron/exon boundaries of GCM2 for all available family members and identified a homozygous G to T transition at nucleotide position 140 (c.140G→T) in exon 2, present in all affected individuals. This substitution results in a change in the amino acid at codon 47, within the DNA binding domain of GCM2, from the normal arginine to leucine. All unaffected parents were heterozygous for this R47L missense mutation.

The G to T transition eliminates an HhaI restriction site, hence we independently confirmed segregation of the mutation within the family by restriction analysis. An 1198 bp fragment containing exons 2 and 3 was digested with HhaI. Normal individuals (homozygous wildtype) produce five fragments of size 352, 324, 276, 186, and 60 bp. Affected individuals IV:2, V:1, and V:2, homozygous for the R47L mutation, lose the cut site at nucleotide position 140 and produce a 246 bp fragment in place of the 186 and 60 bp fragments. Heterozygous parents of affected individuals, carrying one normal and one mutant allele, display all six fragments. Figure 1B and C illustrates the R47L mutation identified and segregation of this mutation within the family. A cohort of 108 chromosomes from healthy unrelated Pakistani control subjects was screened for the presence of the mutation. All controls demonstrated a wildtype genotype, implying that the DNA variant is highly unlikely to represent a sequence polymorphism.

Effect of R47L mutation on DNA binding ability of GCM2

To assess whether the R47L amino acid substitution altered the DNA binding ability of the mutant GCM2 protein, we performed EMSA. Significant levels of expression of both wildtype and mutant GCM proteins following transfection were determined by means of western blotting. Normal levels of both proteins were detected in cell extracts subsequently used for EMSA experiments (fig 2C).

Figure 2

 (A) Sequence alignment of human (GCM1, GCM2), mouse (mGcm1, mGcm2), rat (rGcm1), and fly (dGCM, dGCM2) proteins in the DNA binding domain. Conserved amino acids are boxed; numbers indicate amino acid positions of human GCM2 protein. The arrow identifies the position of the R47L substitution. (B) Predicted structure of normal GCM2 protein backbone, with arginine residues highlighted in ball-and-stick format. The arrow identifies the predicted DNA binding region. Position of the native arginine residue at position 47, substituted for leucine in affected individuals, is indicated. (C) Western blot demonstrating expression of wildtype and mutant GCM2, indicated by the horizontal arrow. Size ladder is shown in kDa. (D) Electrophoretic mobility shift assay showing binding ability of normal and mutant proteins to radio labelled probe containing the GCM DNA binding site. Contents of lanes are indicated above the gel image. Competition assays were performed with 1-, 10-, 20-, and 50-fold molar excess of self and non-self competitors. The R47L mutation abolishes DNA binding ability of the GCM2 protein. Figure indicates DNA-protein complexes (C) and free probes (F).

Using a 32P-labelled oligonucleotide containing the known DNA binding recognition site for GCM proteins, EMSA reactions were performed with extracts from COS cells transfected with full length wildtype GCM2 cDNA or full length GCM2 containing the R47L substitution. Cells transfected with empty vector (pCMV-Tag3B) were used as a negative control. Binding of GCM2 to the recognition site was readily detected for wildtype protein and was competed out with increasing amounts of unlabelled self-competitor probe, indicating the binding to be specific. Protein containing the R47L mutation was unable to bind to the GCM binding site, implying that substitution of leucine for arginine at this residue abolishes normal DNA binding ability of the GCM2 protein (fig 2D).

DISCUSSION

FIH is known to be a heterogeneous disorder with mutations in the PTH gene, encoding parathyroid hormone, being described in both autosomal dominant and autosomal recessive cases.3–5 PTH secretion is normally regulated by the G protein-coupled CASR, functioning as a negative feedback loop; heterozygous defects in the CASR gene resulting in constitutive receptor activation have also been detected in autosomal dominant FIH.7,9 In addition to these two genes described, there has been a single report of a mutation in the Drosophila glial cells missing (gcm) human homologue, GCM2, encoding a transcription factor implicated in parathyroid gland development.

In a highly consanguineous Pakistani family within which autosomal recessive FIH is segregating, we report a novel missense mutation of the GCM2 gene. This is only the second report of a mutation of the GCM2 gene, and the first to identify a critical domain of the GCM2 protein necessary for control of parathyroid tissue development. We used autozygosity mapping, initially to exclude linkage to a known gene causative of autosomal recessive FIH, namely the PTH locus. We next identified a region of homozygosity extending over an interval of 180 kb, proximal to the GCM2 gene (fig 1A). Direct sequence analysis of all codingGCM2 exons showed all three affected individuals to be homozygous for an Arg47Leu amino acid substitution, the first missense mutation to be reported in this gene which encodes a transcription factor. Unaffected parents were heterozygous for the mutation, which co-segregated with the disease phenotype and was absent from normal controls, producing compelling evidence that this missense mutation represents the inherited basis of hypoparathyoidism in this kindred (fig 1B and C). The only previously reported mutation of the GCM2 gene in FIH is a large homozygous deletion of exons 1–4.11 The present study demonstrates that a missense mutation within the conserved DNA binding domain of GCM2 is sufficient to cause an FIH phenotype.

The substituted arginine residue is conserved amongst its mouse and Drosophila homologues, indicating a potentially important and functional role in GCM2 activity (fig 2A). Arginine residues are positively charged, polar amino acids, which play a significant role in structure. They are frequently involved in salt bridges, where they pair with negatively charged amino acids to form stable hydrogen bonds, important for protein stability.19 The replacement of arginine with the similarly positively charged lysine in the GCM1 isoforms provides further evidence of an important role for this amino acid residue within the GCM DNA binding domain. The affected individuals in our FIH kindred have a homozygous substitution of leucine for arginine at position 47. The amino acid leucine, in contrast to arginine residues, is aliphatic and hydrophobic, preferring to be situated within protein hydrophobic cores. In addition, this highly non-reactive amino acid favours a position within alpha helices, rather than in β strands.19 Recently, Cohen et al20 solved the crystal structure of the GCM domain-DNA complex, defining the amino acid residues that make contact with the DNA. The GCM domain exists as two β sheet domains, which form a clamp, seizing the DNA from two sides of the major groove. The upper jaw of the clamp is formed by the β sheet of the larger domain, with its strands orientated orthogonally to the DNA axis. Analysis of the predicted protein structure of GCM2 indicates that the arginine residue substituted in this hypoparathyroidism family is not involved directly in DNA binding (fig 2B). However, EMSA analysis has demonstrated abrogation of DNA binding activity when the arginine to leucine mutation is introduced at amino acid position 47 (fig 2D). This residue lies in the large β sheet domain, but is located on the outer surface of the protein, and does not contact the DNA substrate. Taken together, these data suggest DNA binding activity of GCM2 may be dependent on complex formation with one or more additional proteins, and we hypothesise that disruption of this protein-protein interaction represents the molecular mechanism leading to FIH in this family.

Although the function of human GCM2 remains unknown, studies using mouse Gcm2 and Drosophila GCM2 provide some insight into the expression patterns and transcriptional activity of this protein. In Drosophila, GCM proteins are transiently expressed in glial precursors and act as a binary switch between neuronal and glial cell determination.12,21 From sequence homology, indicating a potential conservation of function, mammalian GCM proteins may be expected to be involved in gliogenesis.22 However, expression studies have shown murine Gcm1 and Gcm2 to be expressed in trophoblast cells during placental development and in developing parathyroid glands, respectively.16,23 Despite containing two separate transactivation domains, Gcm2 is a poor transcriptional activator when compared to Gcm1, and is unable to initiate the neuron-glia transformation when expressed as a transgene in the developing nervous system of Drosophila.16,22,24 Targeted deletion of the mouse Gcm2 gene results in absence of parathyroid glands and undetectable PTH expression in the thyroid and neck of homozygous mice.17 These mice develop hypoparathyroidism but with normal serum PTH levels, insufficient to correct the hypocalcaemia. This compensatory PTH is synthesised in the thymus, where PTH positive cells are shown to express Gcm1 and Casr.17 Of interest, affected members in the family we investigated each had detectable but sub-normal levels of circulating PTH.

In summary, the present study reinforces the critical role of GCM2 activity in human parathyroid gland development through molecular analysis of an autosomal recessive FIH family. Detailed characterisation of GCM2 mutations associated with FIH is warranted and will likely provide insight into the molecular pathways necessary for PTH secreting cells. These findings have more immediate relevance for the diagnosis and management of families presenting with isolated hypoparathyroidism.

Acknowledgments

We thank the family members for their participation.

REFERENCES

View Abstract

Footnotes

  • This work was supported by a Wellcome Trust grant (programme grant 062346/Z/00/Z, Functional Genomics National Autozygosity Mapping Resource). CT is a Royal Society Dorothy Hodgkin Research Fellow.

  • Competing interests: none declared

  • Ethics approval: approval was obtained from the Leicestershire research ethics committee.

Request permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.