Abstract
Free full text
GNB5 Mutations Cause an Autosomal-Recessive Multisystem Syndrome with Sinus Bradycardia and Cognitive Disability
Abstract
GNB5 encodes the G protein β subunit 5 and is involved in inhibitory G protein signaling. Here, we report mutations in GNB5 that are associated with heart-rate disturbance, eye disease, intellectual disability, gastric problems, hypotonia, and seizures in nine individuals from six families. We observed an association between the nature of the variants and clinical severity; individuals with loss-of-function alleles had more severe symptoms, including substantial developmental delay, speech defects, severe hypotonia, pathological gastro-esophageal reflux, retinal disease, and sinus-node dysfunction, whereas related heterozygotes harboring missense variants presented with a clinically milder phenotype. Zebrafish gnb5 knockouts recapitulated the phenotypic spectrum of affected individuals, including cardiac, neurological, and ophthalmological abnormalities, supporting a direct role of GNB5 in the control of heart rate, hypotonia, and vision.
Main Text
Heterotrimeric G proteins trigger a signal transduction cascade composed of α, β, and γ subunits. They are associated with G protein-coupled receptors (GPCRs) in modulating an array of cellular functions, including release of a multitude of hormones and growth factors, regulation of cell contraction and migration, and cell growth and differentiation during development.1, 2, 3, 4 G protein-coupled signaling plays a crucial role in neuronal communication, including regulation of the antagonistic effects of the parasympathetic and sympathetic branches of the autonomic nervous system throughout the body. We report a genetic disorder caused by mutations affecting GNB5 (MIM: 604447), encoding guanine nucleotide-binding protein subunit beta-5, and with disease manifestation in multiple systems.
We identified nine affected individuals (six females and three males) from six unrelated families presenting with a clinical overlap of neurological and cardiac conduction defects; all subjects were found to have variation in the same gene, GNB5, and share a similar rare phenotype. This work results from exome and phenotype data aggregation among independent groups engaged in studying the molecular basis of yet unsolved human genetic rare disease traits. Shared phenotypic features representing the cardinal characteristics of the syndrome include global developmental delay, seizures, generalized hypotonia, retinal disease, and the uncommon feature of early-onset sinus-node dysfunction (Table 1). Additional clinical investigations and diagnostic studies did not show any evidence of structural CNS, ocular, or cardiac anomalies. Affected individuals from four of the six families (families A–D) demonstrated the severe end of the disease spectrum, including substantial cognitive deficits, delayed motor development, severe hypotonia, retinal disease, pathological gastro-esophageal reflux, and sinus-node dysfunction. Affected individuals in families E and F presented with a milder phenotype, including mild intellectual impairment, language delay, and bradycardia (Figure 1, Table 1, Supplemental Note).
Table 1
Family A | Family B | Family C | Family D | Family E | Family F | ||||
---|---|---|---|---|---|---|---|---|---|
II.1 | II.2 | II.1 | II.2 | II.3 | II.2 | II.1 | II.2 | II.1 | |
Gender, age (years) | F, 22 | F, 20 | F, 6 | F, 11 | M, 9 | F, 12 | F, 13 | M, 8 | M, 23 |
Birth weight | 3,580 g (50th percentile) | NA | 2,980 g (15th percentile) | 2,751 g (15th percentile) | NA | 2,845 g (15th percentile) | NA | NA | NA |
Ethnicity | Italy | Italy | Jordan | Puerto Rico | Puerto Rico | India | Morocco | Morocco | Brazil |
Consanguinity | − | − | + | + | + | − | − | − | + |
Altered speech development | + | + | + | + | + | + | + | + | NA |
Verbal understanding | NA | NA | nonverbal | unremarkable | unremarkable | NA | NA | NA | NA |
Lexical production | NA | NA | nonverbal | delayed | delayed | nonverbal | delayed | delayed | NA |
Intellectual disability | + | + | + | + | + | + | mild | mild | mild |
Epilepsy | + | + | + | − | − | + | − | − | − |
Sinus sick syndrome | + | + | + | + | + | increased PR interval (intermittent Weckenbach) | + | + | + |
Minimum heart rate (bpm) | 24 | 39 | NA | paced | paced | NA | 20 | 16 | |
Maximum heart rate (bpm) | 163 | 192 | NA | paced (27% heartbeats on Holter) | paced (20% heartbeats on Holter) | NA | 176 | 180 | NA |
Chronotropic response | NA | NA | NA | + | + | NA | unremarkable | unremarkable | NA |
Escape beats | + | + | NA | paced | paced | NA | + | + | NA |
Pacemaker implantation | − | − | − | + | + | − | − | + | NA |
Heart structural abnormalities | − | PFO | NA | − | − | − | − | − | NA |
Hypotonia | + | + | + | + | + | + | − | impaired fine motor skills | − |
Pathological gastric reflux | + | + | − | + | + | + | − | − | NA |
Nystagmus | + | + | + | + | + | + | NA | − | NA |
Plasma amino acids chromatography | 938 μm/L (restored) | + (restored) | unremarkable | unremarkable | unremarkable | unremarkable | 444 μm/L | unremarkable | NA |
Urine organic acids | unremarkable | unremarkable | increased excretion of 3-methyl-glutaconic acid | unremarkable | unremarkable | unremarkable | NA | NA | NA |
Affected individual numbers refer to those in the pedigree in Figure 1. Complete pedigree charts, consanguinity status, variants, and related homozygous and/or compound heterozygous alleles are reported in Figure 1 and Table S1. Abbreviations are as follows: M, male; f, female; NA, not available; +, clinical trait present; −, clinical trait not present; PFO, patent foramen ovale; bpm, beats per minute.
Given that no potentially pathogenic genomic structural abnormalities were identified by array comparative genomic hybridization and karyotyping of the affected subjects, we applied whole-exome sequencing to all the affected individuals and their healthy parents. Families were recruited in Italy (family A), Brazil (B and F), the United States (C and D), and the Netherlands (E). The institutional review boards of the IRCCS Casa Sollievo Della Sofferenza Hospital, the Hospital das Clínicas da Universidade de São Paulo, the Baylor College of Medicine, the Amsterdam Academic Medical Center, and the University of Lausanne approved this study. Participants were enrolled after written informed consent was obtained from parents or legal guardians. The clinical evaluation included medical history interviews, a physical examination, and review of medical records. To uncover genetic variants associated with the complex phenotype shown by the nine affected subjects, we sequenced their exomes and those of their parents. DNA libraries were prepared from blood-derived genomic DNAs according to standard procedures. Exomes were captured and sequenced with different platforms to reach 50- to 120-fold coverage on average. Variants were called as previously described.5, 6, 7 Variants were filtered on the basis of inheritance patterns, including autosomal recessive, X-linked, and de novo and/or autosomal dominant. Variants with MAF < 0.05% in control cohorts (dbSNP, the 1000 Genome Project, NHLBI GO Exome Sequencing Project, the Exome Aggregation Consortium database, and our in-house databases) and predicted to be deleterious by SIFT,8 PolyPhen-2,9 and/or UMD-Predictor10 were prioritized.
Given a potential history of consanguinity reported in some families (families B, C, and F [Figure 1, Table 1]), we filtered variants by using Mendelian expectations for the assumption of a rare autosomal-recessive trait. We found GNB5 to be compliant with Mendelian expectations and bearing bi-allelic putative deleterious variants in all affected individuals (Figure 1, Table S1). Sanger sequencing in each family confirmed the anticipated segregation of the GNB5 variants. Strikingly, the variants found in the severely affected individuals (families A–D) were predicted to be loss-of-function (LoF) alleles, whereas the more mildly affected individuals from families E and F were homozygous for the same missense variant, c.242C>T (p.Ser81Leu [GenBank: NM_006578.3]) (Figures 1 and S1A). In families B, C, and D the affected individuals were homozygous for splice variants (c.249+1G>T [p.Asp84Leufs31] and c.249G+3A>G [p.Asp84Valfs31]) and a nonsense variant (c.906C>G [p.Tyr302]), respectively (Figures 1 and S1A, Table S1). In family A, the affected siblings were compound heterozygous for a maternally inherited nonsense variant (c.994C>T [p.Arg332]) and a paternally inherited splice-site change (c.249G>A [p.(=)]), which gives rise to an aberrantly spliced isoform containing an additional 25 nucleotides of the intervening intron 2 (Figure S2A). We experimentally show that the transcripts from both alleles are targeted by nonsense-mediated mRNA-decay (Figure S2B).
The five GNB5 LoF variants identified in families A–D are either not present or present with MAF ≤ 8.25 × 10−6 in ExAC (Exome Aggregation Consortium, v.0.3.1) (Table S1). Correspondingly, LoF variants in GNB5 are underrepresented in comparison to expectation in this database; specifically, ExAC reports 8 LoF variants whereas 19 were expected. The c.242C>T (p.Ser81Leu) missense variant identified in family E, of Moroccan ancestry, and family F, of Brazilian ancestry, has a MAF < 5 × 10−5 (6/121,000) in the human population and 4.3 × 10−4 in Latinos (5/11,574). A sample of individuals from Morocco identified a prevalence of 1 out of 1,260 (7.94 × 10−4) for this allele. We estimated the prevalence of the c.242C>T (p.Ser81Leu) variant in the Moroccan population by genotyping a total of 630 Moroccan individuals, including 394 Moroccans and 235 Dutch citizens of Moroccan descent by real-time PCR. Pathogenicity of this variant is further supported by three-dimensional representation of the encoded protein complexed with RGS9, a member of the R7 subfamily of regulators of G-protein signaling (RGS) proteins and common binding partner of GNB5. GNB5 is folded into essentially identical seven-bladed β-propellers (WD40 repeated domains) with equivalent N-terminal helical extensions.11 Replacement of the evolutionarily conserved serine 81 (Figure S1B) by leucine will induce localized structural changes in the immediate vicinity of this residue, which could impair both the central pore of the β-propeller and the binding kinetics of RGS proteins (Figures S3– S5).
In line with the clinical presentation of affected individuals, Gnb5 ablation in mice resulted in marked neurobehavioral abnormalities, including learning deficiencies, hyperactivity, impaired gross motor coordination, abnormal gait,12 defective visual adaptation,13 and perturbed development and functioning of retinal bipolar cells.14 Correspondingly, mice lacking Rgs6, the GNB5-dependent RGS protein enriched in heart tissue, exhibit bradycardia and hypersensitivity to parasympathomimetics.15, 16 To independently investigate the functional effects of variation of GNB5 in the full phenotypic spectrum of subjects reported herein, we engineered a zebrafish model knocked out for gnb5.
CRISPR/Cas9 genome editing was used to generate zebrafish with LoF mutations in gnb5a and gnb5b. This teleost has two GNB5 paralogs as a result of an ancient genome duplication event17 (Figure S6). We identified stable lines with a 7 bp insertion in gnb5a and a 8 bp deletion and 15 bp insertion in gnb5b, causing a frameshift and premature truncation of the encoded proteins, respectively (Figure S7). It was anticipated that gnb5a and gnb5b might have redundant functions, which was confirmed by the absence of overt phenotypes in embryos homozygous for either LoF mutations. As a consequence, a double knockout was generated to ensure complete loss of functional Gnb5. In-crosses of gnb5a and gnb5b double heterozygotes resulted in clutches of embryos containing the expected 6.25% of gnb5a−/−/gnb5b−/− double mutants (henceforth referred to as gnb5 mutants). Consistent with syndrome manifestations of affected individuals, zebrafish mutant embryos had no striking dysmorphologic features (Figure S7D). However, the larvae showed impaired swimming activity, remained small, and generally died 7–14 days post fertilization (dpf), most likely as a result of their inability to feed.
To assess the putative involvement of GNB5 in autonomic nervous system functions, we investigated the GNB5-RGS-GIRK channel pathway. As GNB5 recruits RGS proteins to G protein-coupled inward rectifier potassium (GIRK) channels involved in the hyperpolarization of cell membranes,16, 18 we investigated whether LoF of GNB5 could delay GIRK channel deactivation kinetics, increase hyperpolarization time of cell membranes, and impair cell responsiveness to new stimuli. Carbachol (PubChem CID: 5831) is a parasympathomimetic compound that activates acetylcholine receptors of the heart and the GNB5-RGS-GIRK channel pathway. Treatment of gnb5 mutant larvae with carbachol resulted in a strong decrease of the heart rate, whereas it had little effect on wild-type and sibling larvae (Figure 2), consistent with loss of negative regulation of the cardiac GIRK channel by GNB5-RGS. In contrast, treatment with the sympathetic agonist isoproterenol resulted in an increased heart rate that was similar in wild-type, sibling, and gnb5 mutant larvae (Figure 2). These results indicate that GNB5 is crucial for parasympathetic control of heart rate, but not for sympathetic control, suggesting that lack of GNB5 is associated with extreme bradycardia at rest. Correspondingly, affected individuals present with severe bradycardia at rest (minimal observed heart rates of <25 bpm [beats per minute]) combined with a normal chronotropic response (maximum heart rates >150 bpm).
The severe muscle hypotonia reported in affected individuals could result from GIRK-mediated hyperpolarization of neurons controlling skeletal muscle tone. gnb5 mutant embryos hatched normally from their chorion, a process that requires muscle contraction, but their swimming behavior appeared abnormal at 3 dpf. To investigate whether this abnormal behavior was linked to neurologic dysfunction and hypotonia, we examined the touch-evoked escape response. We anticipated that neurons would only become fully hyperpolarized after an initial stimulus and thus presented the embryos with three consecutive tactile stimuli. Whereas wild-type larvae rapidly swam away in response to repeated tactile stimuli, gnb5 mutants showed a significant decrease in swimming distance and swimming speed at stimuli two (p ≤ 0.0001) and three (p ≤ 0.01), but not after the first stimulus (Figures 3A–3C). Accordingly, gnb5 mutant larvae were predominantly unresponsive to repeated tactile stimuli (Movies S1 and S2). To test whether this abnormal escape response is the consequence of neurologic dysfunction rather than reduced muscle function, we performed a tail movement assay. 5 dpf larvae were given a strong tactile stimulus while we recorded the movement of the tail (Figures 3D and 3E). No significant difference in the maximum tail angle was detected between wild-type and gnb5 mutant larvae (Figure 3E). These results indicate that the tail muscles of gnb5 mutants are fully functional and that the abnormal escape response is associated with neurological dysfunction and possibly muscle hypotonia.
Given that affected individuals have visual problems, including nystagmus, we investigated the visual system by measuring the optokinetic response (OKR) of gnb5 mutant larvae. When wild-type larvae were placed in a drum with a rotating light stimulus (Figure S8A), the OKR consisted of smooth pursuit eye movements followed by rapid rest saccades in the opposite direction (Figure S8B, Movie S3). In contrast, OKR was completely absent in gnb5 mutant larvae although their eyes showed no morphological abnormalities and could make eye movements (Figure S8C, Movie S4). This indicates that the eye muscles are functional in gnb5 mutants but that proper eye-movement control depends on GNB5. Overall these data show that gnb5 mutants faithfully recapitulate the phenotypic spectrum of affected individuals, including cardiac, neurologic, and ophthalmologic abnormalities.
These results provide evidence for a direct role of GNB5 in the control of heart rate, motor capacity, and vision. Whereas GNB1 (MIM: 139380), GNB2 (MIM: 139390), GNB3 (MIM: 139130), and GNB4 (MIM: 610863) are widely expressed and encode highly homologous proteins,20 GNB5 is preferentially expressed in the brain and nervous system and encodes a peptide with less homology with its four paralogs.21, 22
Germline de novo GNB1 variants cause severe neurodevelopmental disability,23 hypotonia, and seizures. GNB3 bi-allelic LoF has been linked to congenital stationary night blindness (MIM: 610445, 163500, 610444, 613830, 616389, 310500, 257270, 613216, 614565, 615058, 300071, and 610427) and recessive retinopathy in humans,24, 25 retinal degeneration in chickens,26 and reduced cone sensitivity and mild bradycardia in mice.27, 28 A SNP in GNB3 was associated with postural tachycardia syndrome29 and incidence of cardiovascular disease and stroke.30 Similarly, GNB2 and GNB4 map to loci governing heart rate on chromosomes 7 and 3, respectively.31, 32 We hereby demonstrate that bi-allelic LoF and missense variants in GNB5 cause a multisystem syndrome with features that include global developmental delay, sinus-node dysfunction, seizures, eye abnormalities, gastric problems, and generalized hypotonia. We highlight the importance of GNB5 for neuronal signaling, including the regulation of the antagonistic effects of the parasympathetic and sympathetic nervous system.
Conflicts of Interest
J.R.L. has stock ownership in 23andMe, is a paid consultant for Regeneron Pharmaceuticals, has stock options in Lasergen, is a member of the Scientific Advisory Board of Baylor Miraca Genetics Laboratories, and is a co-inventor on multiple United States and European patents related to molecular diagnostics for inherited neuropathies, eye diseases, and bacterial genomic fingerprinting. Baylor College of Medicine (BCM) and Miraca Holdings have formed a joint venture with shared ownership and governance of the Baylor Miraca Genetics Laboratories (BMGL), which performs clinical exome sequencing. The Department of Molecular and Human Genetics at BCM derives revenue from the chromosomal microarray analysis and clinical exome sequencing offered in the BMGL (http://www.bmgl.com/BMGL/Default.aspx/website). G.M. is a paid consultant for Takeda Pharmaceutical Company.
Footnotes
Supplemental Data includes Supplemental Acknowledgments, a Supplemental Note, eight figures, one table, and four movies and can be found with this article online at http://dx.doi.org/10.1016/j.ajhg.2016.06.025.
Web Resources
1000 Genomes, http://www.1000genomes.org
Berkeley Drosophila Genome Project NNSplice 0.9, http://www.fruitfly.org/seq_tools/splice.html
Burrows-Wheeler Aligner, http://bio-bwa.sourceforge.net/
Ensembl Genome Browser, human genome, GRCh37, http://grch37.ensembl.org/Homo_sapiens/Info/Index
ExAC Browser, http://exac.broadinstitute.org/
ExomeDepth, https://cran.r-project.org/web/packages/ExomeDepth/index.html
GenBank, http://www.ncbi.nlm.nih.gov/genbank/
GraphPad, http://graphpad.com/
NetGene2, http://www.cbs.dtu.dk/services/NetGene2
NHLBI Exome Sequencing Project (ESP) Exome Variant Server, http://evs.gs.washington.edu/EVS/
OMIM, http://www.omim.org/
PolyPhen-2, http://genetics.bwh.harvard.edu/pph2/
SIFT, http://sift.jcvi.org/
SnpEff, http://snpeff.sourceforge.net/
Swiss PDB Viewer, http://www.expasy.org/spdbv/
UMD-Predictor, http://umd-predictor.eu/
Supplemental Data
Video demonstrating normal touch-evoked escape responses in a 3 dpf wild-type zebrafish.
Video demonstrating abnormal touch-evoked escape responses in a 3 dpf gnb5 mutant.
Video demonstration of a normal optokinetic response in a 7 dpf wild-type larvae, consisting of smooth pursuit eye movements followed by rapid rest saccades in the opposite direction.
Video demonstration of an abnormal optokinetic response in a 7 dpf gnb5 mutant.
References
Articles from American Journal of Human Genetics are provided here courtesy of American Society of Human Genetics
Full text links
Read article at publisher's site: https://doi.org/10.1016/j.ajhg.2016.06.025
Read article for free, from open access legal sources, via Unpaywall: http://www.cell.com/article/S0002929716302294/pdf
Citations & impact
Impact metrics
Citations of article over time
Alternative metrics
Article citations
The association of GNB5 with Alzheimer disease revealed by genomic analysis restricted to variants impacting gene function.
Am J Hum Genet, 111(3):473-486, 13 Feb 2024
Cited by: 1 article | PMID: 38354736
Prioritizing cardiovascular disease-associated variants altering NKX2-5 and TBX5 binding through an integrative computational approach.
J Biol Chem, 299(12):105423, 04 Nov 2023
Cited by: 0 articles | PMID: 37926287 | PMCID: PMC10750078
Visualizing the chaperone-mediated folding trajectory of the G protein β5 β-propeller.
Mol Cell, 83(21):3852-3868.e6, 17 Oct 2023
Cited by: 0 articles | PMID: 37852256
How genomic insights into the evolutionary history of clouded leopards inform their conservation.
Sci Adv, 9(40):eadh9143, 06 Oct 2023
Cited by: 0 articles | PMID: 37801506 | PMCID: PMC10558132
Genetics of sinoatrial node function and heart rate disorders.
Dis Model Mech, 16(5):dmm050101, 17 May 2023
Cited by: 4 articles | PMID: 37194974 | PMCID: PMC10214849
Review Free full text in Europe PMC
Go to all (31) article citations
Data
Data behind the article
This data has been text mined from the article, or deposited into data resources.
BioStudies: supplemental material and supporting data
Diseases (Showing 17 of 17)
- (1 citation) OMIM - 257270
- (1 citation) OMIM - 616389
- (1 citation) OMIM - 615058
- (1 citation) OMIM - 163500
- (1 citation) OMIM - 610427
- (1 citation) OMIM - 614565
- (1 citation) OMIM - 613830
- (1 citation) OMIM - 610863
- (1 citation) OMIM - 610445
- (1 citation) OMIM - 613216
- (1 citation) OMIM - 139130
- (1 citation) OMIM - 610444
- (1 citation) OMIM - 139380
- (1 citation) OMIM - 310500
- (1 citation) OMIM - 139390
- (1 citation) OMIM - 300071
- (1 citation) OMIM - 604447
Show less
RefSeq - NCBI Reference Sequence Database
- (2 citations) RefSeq - NM_006578.3
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Inhibition of G-protein signalling in cardiac dysfunction of intellectual developmental disorder with cardiac arrhythmia (IDDCA) syndrome.
J Med Genet, 58(12):815-831, 10 Nov 2020
Cited by: 2 articles | PMID: 33172956 | PMCID: PMC8639930
[Intellectual developmental disorder with cardiac arrhythmia syndrome in a family caused by GNB5 variation and literature review].
Zhonghua Er Ke Za Zhi, 58(10):833-837, 01 Oct 2020
Cited by: 2 articles | PMID: 32987464
Review
Genotype-phenotype correlation of congenital anomalies in multiple congenital anomalies hypotonia seizures syndrome (MCAHS1)/PIGN-related epilepsy.
Am J Med Genet A, 170A(1):77-86, 23 Sep 2015
Cited by: 24 articles | PMID: 26394714 | PMCID: PMC4886552
Unique retinal signaling defect in GNB5-related disease.
Doc Ophthalmol, 140(3):273-277, 12 Nov 2019
Cited by: 6 articles | PMID: 31720979
Funding
Funders who supported this work.
NINDS NIH HHS (1)
Grant ID: K23 NS078056
Telethon (1)
Telethon Network of Genetic Biobanks
Dott.ssa Chiara Baldo, E. O. OSPEDALI GALLIERA, GENOVA
Grant ID: GTB12001