Article Text
Statistics from Altmetric.com
The 100 000 Genomes Project (100KGP) is a UK-wide initiative that has a goal of using whole genome sequencing (WGS) to identify genetic causes of rare inherited diseases and embed the use of this technology within the NHS.1 Using data from this resource alongside international gene-matching efforts, four individuals from two independent families were identified harbouring homozygous frameshift or stop-gain variants in PRKG2, a recently described skeletal dysplasia gene.2 Detailed clinical and radiological assessments helped extend the phenotypic range associated with this autosomal recessive condition while functional studies indicated that both variants had a similar impact on FGF-induced MAPK signalling.
PRKG2 encodes the cyclic guanosine monophosphate dependent protein kinase II (cGKII), which acts downstream of the natriuretic peptide receptor-B/C- natriuretic peptide (NPR-B/CNP). NPR-B is encoded by NPR2, biallelic variants in which are responsible for acromesomelic dysplasia, Maroteaux type (AMDM; MIM 602875). Rodent models further implicate PRKG2 in skeletal development3 4 and cGKII deficiency was shown to be the cause of the dwarfism phenotype observed in Angus cattle.5 Building on support from pathway analysis and model organisms, a recent study showed that biallelic PRKG2 variants can result in acromesomelic dysplasia, PRKG2-type (AMDP) in humans,2 adding PRKG2 to a list of >400 genes associated with genetic skeletal disorders.6 As only two affected individuals were reported, it is important that the full clinical range of this condition is described.
In this study, we searched for rare biallelic PRKG2 variants using data from the 100KGP via the LabKey application available within Genomic England’s research environment. Researchers can apply for access online (www.genomicsengland.co.uk/join-a-gecip-domain). Initial filtering employed a 1% population allele frequency threshold based on data from the 1000 Genomes Project as well as in-house frequency information. An additional family was identified via a network of collaborators and variants were classified using ACMG criteria (online supplemental table 1).
Supplemental material
In family 1, WGS and subsequent Sanger sequencing uncovered a homozygous pathogenic PRKG2 variant, NM_006259.3:c.2282dup (p.Asp761Glufs*34; online supplemental figure 1) in three brothers referred with spondylometaphyseal dysplasia (figure 1A). Interestingly, the middle-affected brother (F1-IV-6) also has type I osteogenesis imperfecta (OI). An early clinical exome sequencing study found that for 4.6% of cases with a molecular diagnosis, more than one gene was contributing to a blended phenotype.7 Complex cases such as these are expected to be more common in highly consanguineous families where large regions of homozygosity (ROHs) make up a significant proportion of the genome; however, for F1-IV-6 the secondary diagnosis of OI was due to a COL1A1 frameshift, which had arisen de novo. OI was suspected in this child because of multiple fractures in childhood (arm as an infant, wrist aged 8 and thoracic T6 wedge fracture) combined with blue sclerae. It is certainly possible that the coexistent OI may have had an impact on the severity of the phenotype in this individual, not least because his height was more significantly reduced than for his two brothers and OI (type 1) is a known cause of reduced stature in its own right.
Supplemental material
In family 2, exome sequencing for a girl with acromesomelic dysplasia revealed a homozygous pathogenic PRKG2 variant c.1705C>T; p.(Arg569*) (online supplemental figure 2), observed previously in a patient with similar clinical and radiological features.2 Comparison of the available genomic data for F2-V-3 and the previously published case was not able to detect a shared haplotype across the PRKG2 locus. However, exome sequencing has limited resolution to detect small regions of identity by descent and so a founder mutation cannot be ruled out. Given the differing ethnicities and the fact that c.1705C>T lies at a CpG dinucleotide, the recurrence of c.1705C>T being due to separate mutational events seems a more likely scenario.
Both variants described here are extremely rare; p.(Asp761Glufs*34) is absent from gnomAD (https://gnomad.broadinstitute.org), while p.(Arg569*) is present as a singleton allele. In both families, the disease-causing variants lay within large ROHs (online supplemental table 2). Pathogenic variants are overrepresented in the largest ROHs and it has been proposed that lying in one of the top 10 such regions can be used as evidence supporting pathogenicity.8 While the p.(Arg569*) variant has already been demonstrated to affect the downstream MAPK pathway,2 the p.(Asp761Glufs*34) in family 1 is likely to be disruptive given the switch of the final Asp-Phe residues for 33 alternative amino acids at the C terminus. In silico modelling highlights the structural importance of this region, in particular the final Phe762 residue (figure 1B, supplementary methods; interactive version at https://michelanglo.sgc.ox.ac.uk/r/prkg2).
To functionally confirm the pathogenicity of the newly identified p.(Asp761Glufs*34) variant, we first analysed cGKII expression by western blot analysis. Plasmid construction for p.(Asp761Glufs*34) involved a sequential PCR strategy (supplemental methods), with the previously characterised p.(Arg569*) variant employed as a positive control. For both variants, cGKII was detected at the predicted size (figure 1C), although at dramatically reduced levels (≥80%) compared with the wild type (figure 1D). Next, we evaluated whether the p.Asp761Glufs*34 mutant was able to inhibit FGF2-induced MAPK pathway by analysing its ability to induce phosphorylation of Raf-1 at Ser-43 and ERK1/2, as described previously.2 Wild-type cGKII downregulated MAPK signalling by reducing ERK1/2 activation through the upstream phosphorylation of Raf-1 at Ser-43 in a cGMP-dependent manner. However, the p.Asp761Glufs*34 mutant failed to phosphorylate Raf-1 at Ser-43 and thus, reduced FGF2-induced ERK1/2 phosphorylation (figure 1E–F), similar to results for the p.Arg569* variant.2
Detailed phenotypic information is provided for both families and compared with the two published cases (online supplemental table 2, figure 3). Radiological findings for F2-V-3 were very similar to those observed for ‘Proband 1’ described previously,2 which is unsurprising given that both individuals harbour the same homozygous p.(Arg569*). In contrast, for PRKG2 family 1 there was a consistent radiological phenotype distinct from previously reported AMDP and AMDM. The three brothers reported here (F1-IV-3, IV-6 and IV-7) had no evidence of acromesomelic shortening, except for mild shortening of toes observed for individual F1-IV-7. The main findings were platyspondyly with anterior vertebral body projections, long slender femoral necks and some metaphyseal irregularity (most evident in the radius and ulna) and striations (figure 2). The metaphyses of the distal phalanges were somewhat cone-shaped in one child, but not pronounced, generalised or associated with shortening, as seen in AMDM.9 In summary, family 1 exhibited a skeletal phenotype characterised by spondylometaphyseal dysplasia, rather than acromesomelic dysplasia as expected in AMDP and AMDM.
Interestingly, the PRKG2 locus has been identified in several genome-wide association studies on height (www.ebi.ac.uk/gwas/genes/PRKG2). Therefore, the description of this now confirmed Mendelian condition constitutes an additional example of rare variants in a gene causing a severe condition, where common variants in the same gene are associated with a related trait.10 In summary, analysis of 100KGP data combined with gene-matching efforts identified four affected individuals with biallelic loss of function variants in PRKG2, extending the phenotypic range of this condition to include spondylometaphyseal dysplasia. The patients described here were the only individuals harbouring severe biallelic PRKG2 variants across all rare disease areas within the 100KGP. These data include 295 patients recruited due to an unexplained skeletal dysplasia and therefore our results are consistent with this condition being extremely rare in humans.
Ethics statements
Patient consent for publication
Ethics approval
The 100K Genomes Project has approval from the HRA Committee East of England (Cambridge South REC: 14/EE/1112).
Acknowledgments
We thank the two families for their involvement in our study. This research was made possible through access to the data and findings generated by the 100,000 Genomes Project. The 100KGP is managed by Genomics England Limited (a wholly owned company of the Department of Health and Social Care). The 100KGP uses data provided by patients and collected by the National Health Service as part of their care and support.
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Collaborators The Genomics England Research Consortium (27th May 2021): John C. Ambrose (Genomics England, London, UK); Prabhu Arumugam (Genomics England, London, UK); Roel Bevers (Genomics England, London, UK); Marta Bleda (Genomics England, London, UK); Freya Boardman-Pretty (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Christopher R Boustred (Genomics England, London, UK); Helen Brittain (Genomics England, London, UK); Mark J Caulfield (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Georgia C Chan (Genomics England, London, UK); Greg Elgar (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Tom Fowler (Genomics England, London, UK); Adam Giess (Genomics England, London, UK); Angela Hamblin (Genomics England, London, UK); Shirley Henderson (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Tim J P Hubbard (Genomics England, London, UK); Rob Jackson (Genomics England, London, UK); Louise J Jones (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Dalia Kasperaviciute (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Melis Kayikci (Genomics England, London, UK); Athanasios Kousathanas (Genomics England, London, UK); Lea Lahnstein (Genomics England, London, UK); Sarah E A Leigh (Genomics England, London, UK); Ivonne U S Leong (Genomics England, London, UK); Javier F Lopez (Genomics England, London, UK); Fiona Maleady-Crowe (Genomics England, London, UK); Meriel McEntagart (Genomics England, London, UK); Federico Minneci (Genomics England, London, UK); Loukas Moutsianas (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Michael Mueller (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Nirupa Murugaesu (Genomics England, London, UK); Anna C Need (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Peter O’Donovan (Genomics England, London, UK); Chris A Odhams (Genomics England, London, UK); Christine Patch (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Mariana Buongermino Pereira (Genomics England, London, UK); Daniel Perez-Gil (Genomics England, London, UK); John Pullinger (Genomics England, London, UK); Tahrima Rahim (Genomics England, London, UK); Augusto Rendon (Genomics England, London, UK); Tim Rogers (Genomics England, London, UK); Kevin Savage (Genomics England, London, UK); Kushmita Sawant (Genomics England, London, UK); Richard H Scott (Genomics England, London, UK); Afshan Siddiq (Genomics England, London, UK); Alexander Sieghart (Genomics England, London, UK); Samuel C Smith (Genomics England, London, UK); Alona Sosinsky (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Alexander Stuckey (Genomics England, London, UK); Mélanie Tanguy (Genomics England, London, UK); Ana Lisa Taylor Tavares (Genomics England, London, UK); Ellen R A Thomas (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Simon R Thompson (Genomics England, London, UK); Arianna Tucci (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Matthew J Welland (Genomics England, London, UK); Eleanor Williams (Genomics England, London, UK); Katarzyna Witkowska (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK); Suzanne M Wood (Genomics England, London, UK, William Harvey Research Institute, Queen Mary University of London, London, EC1M 6BQ, UK).
Contributors ATP and JCT conceived the project. ATP, BB-P, MPF, FD-G, TL, RM and KEH performed data analysis. BB-P and FD-G generated experimental data. The 57 members of GERC performed WGS and provided access data via a secure research environment. MBT, ADC, EGK, MD, AW, PW, KB, KE, RM, HT and DS recruited the patients and collected/interpreted clinical information. ATP drafted the manuscript, which was revised and approved by all.
Funding This research was supported by the National Institute for Health Research (NIHR) Oxford Biomedical Research Centre Programme, the Wellcome Trust (203141/Z/16/Z) and the European Union’s Seventh Framework Programme for research, technological development and demonstration under grant agreement no. 608473 (to HT). Additional support was from MINECO grant no. SAF2017-84646-R from (to KEH) and via an FPU studentship (to FD-G) from the Spanish Ministry of Education. The 100KGP is funded by the National Institute for Health Research and NHS England. The Wellcome Trust, Cancer Research UK and the Medical Research Council have also funded research infrastructure.
Competing interests None declared.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.