Skip to main content

Advertisement

A novel NAA10 p.(R83H) variant with impaired acetyltransferase activity identified in two boys with ID and microcephaly

Abstract

Background

N-terminal acetylation is a common protein modification in human cells and is catalysed by N-terminal acetyltransferases (NATs), mostly cotranslationally. The NAA10-NAA15 (NatA) protein complex is the major NAT, responsible for acetylating ~ 40% of human proteins. Recently, NAA10 germline variants were found in patients with the X-linked lethal Ogden syndrome, and in other familial or de novo cases with variable degrees of developmental delay, intellectual disability (ID) and cardiac anomalies.

Methods

Here we report a novel NAA10 (NM_003491.3) c.248G > A, p.(R83H) missense variant in NAA10 which was detected by whole exome sequencing in two unrelated boys with intellectual disability, developmental delay, ADHD like behaviour, very limited speech and cardiac abnormalities. We employ in vitro acetylation assays to functionally test the impact of this variant on NAA10 enzyme activity.

Results

Functional characterization of NAA10-R83H by in vitro acetylation assays revealed a reduced enzymatic activity of monomeric NAA10-R83H. This variant is modelled to have an altered charge density in the acetyl-coenzyme A (Ac-CoA) binding region of NAA10.

Conclusions

We show that NAA10-R83H has a reduced monomeric catalytic activity, likely due to impaired enzyme-Ac-CoA binding. Our data support a model where reduced NAA10 and/or NatA activity cause the phenotypes observed in the two patients.

Background

N-terminal acetylation (Nt-acetylation) is an exceedingly prevalent protein modification, affecting an estimated 80% of human proteins [1]. An acetyl group donated by acetyl coenzyme A (Ac-CoA) is attached to the N-terminal amine group of a substrate protein through the catalytic action of N-terminal acetyltransferases (NATs), which are members of the GCN5-related N-acetyltransfererase (GNAT) superfamily [2]. Nt-acetylation may affect substrates in several ways, although for most substrates the importance of the Nt-acetyl group remains enigmatic. Nt-acetylation is reported to affect the half-life of substrates, destabilizing some substrates [3,4,5] and stabilizing others [6, 7]. Nt-acetylation also regulates folding [8,9,10,11], protein complex formation [12,13,14,15,16,17], subcellular localization and membrane interaction [16,17,18,19]. While several post-translational NATs have recently been described [20,21,22,23,24,25], most substrates appear to be Nt-acetylated at the ribosome by one of at least five ribosome-associated NATs – NatA, NatB, NatC, NatD, or NatE. Chief among these in terms of substrate diversity is NatA [1]. The human NatA complex is composed of the catalytic subunit NAA10 and the auxiliary subunits NAA15, NAA50 and HYPK [11, 26,27,28,29]. The ribosome-associated NatA complex Nt-acetylates an estimated 40% of the human proteome [30], preferring small, polar amino acids which are exposed after initiator methionine excision [26]. NAA10 is mainly found in complex with NAA15, both ribosome-and non-ribosome associated [31, 32]; however, there also exists a pool of non-NAA15 bound NAA10 [31]. Due to a conformational change in NAA10 upon NAA15 association, NAA10 alone has an altered substrate specificity, preferring acidic N-termini [31, 33]. Further, non-NAT functions of NAA10 as a lysine acetyltransferase (KAT) or as a non-catalytic regulator have been proposed [34,35,36,37,38]. Loss of NAA10 leads to developmental defects or embryonic death in several model organisms [34, 39,40,41]. Hereditary or de novo germline variants in the X-chromosomal NAA10 gene is associated with developmental syndromes and non-syndromic developmental delay in humans. A NAA10 S37P missense mutation is the cause of the Ogden syndrome, an extremely rare disease in which affected boys have an aged appearance, craniofacial anomalies, cardiac problems including arrhythmia, and where all affected boys have died by age 16 months [42]. This mutation affects NatA complex formation and leads to lowered cell proliferation, larger cell size and reduced Nt-acetylation of some NatA substrates [6, 42]. One splice-donor mutation was found to lead to Lenz microphthalmia syndrome, causing small or missing eyes, intellectual disability and skeletal, cardiac, and renal problems [43]. Various other mutations lead to non-syndromic developmental delay and seizures in males and females [44, 45], a novel intellectual disability syndrome in two brothers carrying the same mutation [46], intellectual disability, developmental delay and cardiac abnormalities in three brothers from two families [47], and non-syndromic intellectual disability with delayed language and motor development in a female proband [48]. The NAA10 c.247C > T p.R83C missense mutation is recurrent, previously having appeared de novo in one male and seven female patients, generally manifesting with moderate to severe intellectual disability and developmental delay, though only the boy had EEG anomalies [49]. While NAA10 mutations have a heterogenous clinical picture, with no clear genotype-catalytic activity-phenotype correlation [47], some features are seen in many or most patients; intellectual disability, developmental delay, growth failure, and cardiac anomalies. Here, we present the c.248G > A p.R83H variant, found in two boys, aged 15 and 12 with hyperactivity, limited language development, developmental delay, intellectual disability and hypertrophic cardiomyopathy. The NAA10 R83H mutation leads to a substantial decrease in NAA10 catalytic activity, supporting the hypothesis that this variant causes a loss of NAA10-mediated acetylation and is the cause of the observed phenotypes. Based on structural models of the variant, we predict that this reduced catalytic activity is due to impaired Ac-CoA binding.

Methods

Trio exome sequencing

A trio-based whole-exome sequencing approach was undertaken. For patient 1, whole exome sequence was performed as described [50]. The NAA10 variant was verified by targeted Sanger sequencing. DNA from patient 2 and parents were subjected to exome capture using NimbleGen SeqCap EZ MedExome (Roche), followed by sequencing on an Illumina NextSeq550 to a mean coverage of 91x, with 94% of targeted bases covered with minimum 20x coverage. Raw reads were aligned using the Burrows-Wheeler Alignment tool (BWA-MEM) v. 0.7.15 [51] and the GATK Best Practice pipeline v. 3.8–0 was used for variant calling [52]. Annotation and filtering of variants was performed using VarSeq 1.5.0 (Golden Helix). The NAA10 variant was verified by targeted Sanger sequencing. Informed consents were obtained from patient indexes and family members.

Multiple sequence alignment, conservation scores and structural model

Multiple sequence alignments were created using ClustalX [53] and the illustration and conservation scores using Jalview [54]. The input sequences were: human NAA10 (Uniprot ID: P41227), Mus musculus Naa10 (Uniprot ID: Q9QY36), Rattus norvegicus Naa10 (Uniprot ID: D3ZUQ2), Xenopus laevis naa10 (Uniprot ID: Q7ZXW3), Danio rerio naa10 (Uniprot ID: Q7T3B8) and Saccharomyces cerevisia ARD1 (Uniprot ID: P07347). The NAA10 structural model is from the human NatA complex (PDB ID: 6C9M) [55]. The image was obtained by aligning this structure to the structure of Schizosaccharomyces pombe (4KVM) [33], which was solved with CoA and the substrate peptide SASE. The structure alignment and figure generation was performed using PyMol [56].

In silico mutagenesis

NAA10-WT and NAA10-R83H from the human NatA complex (PDB ID: 6C9M [55]) were aligned in PyMOL to NAA10 from the S. pombe NatA complex (PDB ID: 4KVM [33]). In order to calculate the electrostatic potential of the two variants, these chains were uploaded to the Adaptive Poisson-Boltzmann Solver (APBS) PDB2PQR webserver [57,58,59] PDB2PQR makes preparatory changes to the PDB file by optimizing hydrogen bond, repairing heavy atoms, assessing pKa values and assigning charge and radius parameters [59]. The resulting PQR file was then solved for continuum electrostatic potential by the APBS software. Output files from the APBS were visualized in PyMOL and the chains were aligned. Both PDB files were prepared and calculated with the same parameters, optimized for cellular conditions. The molecular mechanical forcefield, Amber [60] and the heuristic pKa calculation software PROPKA [61] were used to generate the PQR at pH 7.4. Following the PQR preparations, APBS was set to utilize a manually-configured multigrid calculation, where the automatically suggested dimensions were applied. “Ionic strength of the solvent” was set to that of 0.15 M NaCl with ion radius of 1.5 Å. No electrostatic energies were calculated, and the output was set to “write out the electrostatic potential in units of kbT/ec (multigrid and finite element)”. The APBS was then launched with the remaining parameters unchanged. The output PQR file from the APBS was then opened in PyMOL and visualized in surface view.

Preparation of plasmids

In order to study the NM_003491.3 c.248G > A, p.(R83H) missense variant, a bacterial expression vector pETM-41/His-MBP-NAA10 was modified by site-directed mutagenesis (Q5® Site Directed Mutagenesis Kit, New England Biolabs) according to the manufacturer’s protocol. The primers used for mutagenesis were NAA10 G248A p.R83H fwd (5′-CAACCTCCTCGGCCTGCCCGA) and NAA10 G248A p.R83H rev (CAGTGCTGCATGTTCATTAGGTC), with an annealing temperature of 69 °C. The plasmid was verified by sequencing.

Protein expression in E.coli BL21 cells

To study the impact of the novel variant on in vitro catalytic activity, a two-step purification of His/MBP-NAA10-WT and His/MBP-NAA10-R83H expressed in BL21 Star DE3 E. coli was carried out essentially as described [44].

In vitro colorimetric acetylation assays

The catalytic activity of the novel NAA10-R83H variant was compared with NAA10-WT by performing both time dependent and substrate dependent DTNB-based acetylation assays. In the substrate dependent assay, 50 μL reactions containing 100 nM His-MBP-NAA10-WT or NAA10-R83H, 300 μM peptide substrate (either EEEIA, DDDIA, MLGPE or SESSS, short for EEEIAALRWGRPVGRRRRPVRVYP,DDDIAALRWGRPVGRRRRPVRVYP, MLGPEGGRWGRPVGRRRRPVRVYP, and SESSSKSRWGRPVGRRRRPVRVYP, where the part in bold is the variable N-terminus and the rest of the amino acid sequence is identical between the peptides), 300 μM Ac-CoA, 1 x acetylation buffer (50 mM Tris-HCl, 1 mM EDTA and 10% glycerol, pH 8.5), incubated for 20 min at 37 °C. Blank reactions incubated in the absence of enzyme. The reactions were stopped using 100 μL quenching buffer (3.2 M guanidine-HCl, 100 mM Na2HPO4, pH 6.8), at which time the blanks were added enzyme in equal proportion to the reactions. To indicate the degree of acetylation for the individual reactions, 25 μL saturated solution of DTNB in DTNB buffer (100 mM Na2HPO4, 10 mM EDTA, pH 6.8) was added to each reaction. The time dependent assay was carried out using only EEEIA as a substrate, and reactions were quenched at different times (after 10, 20 and 30 min); the procedures for the assays were otherwise the same. Absorbance was measured at 412 nm using TECAN Infinite® 200 PRO plate reader.

Results

Description of patients

Patient 1 is the only child born to non-consanguineous parents. His birth weight was 3.25 kg (25th centile). There were no antenatal or neonatal complications. His father has microcephaly with mild learning difficulties; his mother attended mainstream school but left school early. Patient 1 presented with developmental delay at age 2. Early milestones: he sat at 9 month, walked age 2 years 2 months, he had very poor speech development, only has ~ 300 words with 2–3 worded sentences aged 15 years. He was diagnosed with coeliac disease aged 5 years 5 months (his father also has coeliac disease). He had short stature but this resolved once treated for coeliac disease. He has transient neutropenia of infancy which resolved. He developed epilepsy aged 13 years 10 months. He wears glasses for astigmatism. He has a history of chronic constipation. He has behavioural issues including poor concentration, and can be quite volatile and aggressive. He has been assessed as having a moderate learning disability & attends a special school. His cardiac assessment revealed a normal cardiac examination aged 2 years 9 months and again 15 years with a normal long QT interval. However, he had been noted to have mild concentric left ventricular hypertrophy on an assessment in between. He is hypermobile. He takes sodium valproate for epilepsy and movicol for his constipation. At 2 years 1 month of age his weight was 11.4 kg (9th centile)/height 80.5 cm (> 0.4th centile), and 24.6 kg (>3rd centile)/height 124.4 (> 0.4th centile) cm at 10 years 7 months. Occipitofrontal circumference (OFC) aged 14 years 7 months 52.5 cm (<10th centile) (His father’s OFC is 53.5 cm (< 0.4th centile). He has a tented upper lip but no other dysmorphism. Trio exome revealed a maternally inherited NAA10 (NM_003491.3) c.248G > A, p.R83H variant which was further confirmed by Sanger sequencing. The mother’s parents are both dead precluding further segregation studies; however her two healthy brothers tested negative for this variant.

Patient 2 is the second child to non-consanguineous parents. He was born after 39 weeks of gestation. Birth weight 3344 g (median), length 51 cm (median), and head circumference 36 cm (+1SD). He had neonatal-onset hypotonia and poor feeding. He was breastfed to some extent until he was five years old. He presented with developmental delay at 12 months. Early milestones: he sat at 6 month, walked at age 2 years 3 months, from 18 months he rolled around indoors. He had very poor speech development and only uses a limited number of words with 2–3 worded sentences aged 12 years. He has behavioural issues including poor concentration, but is a quite happy and very active boy. He had eczema and very fine and sparse scalp hair until 3 years of age. He has a tendency to develop mild fevers of unknown origin. He does not have epilepsy. His height has been constantly at -2SD since 2 years of age; at 12 years his height was 142 cm. Dysmorphic features at 12 years of age was rather mild in the form of large ears, and rather closely spaced eyes. Cardiac examination revealed a structurally normal heart, i.e. no congenital malformations. Measurements of ventricular wall thickness were taken at end-diastole, indexed to body surface area and z-score measurements calculated as described by the Pediatric Heart Network [62]. This revealed a predominantly septal hypertrophy (z-score + 8.3) and reduced end-diastolic diameter of the left ventricle (Fig. 1). Although the septal hypertrophy caused narrowing of the left ventricular outflow tract, it did not cause obstruction to flow (Fig. 1 A and B, Additional file 1: Video S1). The ECG was normal except for borderline prolongation of the corrected QT-interval (420 milliseconds, Fig. 1 C). To assess for potential arrhythmias a 48-h Holter monitoring was performed. The heart rate varied between 66 to 165 bpm, and on average was 106 bpm. No arrhythmias were detected. Trio exome sequencing revealed a mosaic de novo NAA10 ((NM_003491.3) c.248G > A, p.(Arg83His)) hemizygote missense variant, with a mosaic degree of 75%. The variant was confirmed by Sanger sequencing.

Fig. 1
figure1

Cardiac hypertrophy in patient 2. Transthoracic echocardiography showing the left ventricle in a parasternal long axi and b apical 3-chamber view. Please note how the hypertrophy predominantly affects the interventricular septum with relative sparring of the apical and posterior regions. Red marker on EKG trace shows temporal relation to heart cycle. LV, left ventricle; Ao, aorta; LA, left atrium; IVS, interventricular septum. c Standard electrocardiogram (EKG) in patient 2. Note the borderline prolonged QTc

Functional testing

To assess the effect of the R83H mutation on NAA10 catalytic activity, we performed site-directed mutagenesis to obtain the R83H mutant NAA10 protein. We expressed wild-type NAA10 (NAA10-WT) and NAA10-R83H as 6xHis-tagged maltose-binding protein (MBP)-fusion proteins in E. coli, purified them, and performed in vitro Nt-acetylation assays against model substrate peptides. These peptides have identical C-termini but differ in the 7 N-terminal amino acids. Such setups are commonly used to determine the substrate preferences of NATs, and to measure catalytic activities of NAT mutant enzymes [35, 48, 49]. We were able to purify NAA10-R83H in monomeric form under the same conditions as NAA10-WT. Both WT and R83H eluted at around 80 mL, and the R83H mutant was at the expected size. In vitro, monomeric NAA10, not bound to the NatA complex, prefers substrates with acidic N-termini [31], represented in our panel by EEEIA and DDDIA. The NatA canonical substrate SESSS and the NatE substrate MLGPE were also included. We found no difference in the preference for different peptides between NAA10-WT and NAA10-R83H, but the catalytic efficiency was severely diminished in the R83H mutant (Fig. 2a). Timecourse experiments with the EEEIA peptide confirmed that the catalytic activity is significantly reduced (Fig. 2c).

Fig. 2
figure2

Functional testing of NAA10-R83H. MBP-NAA10 fusion proteins (wild-type (WT) or R83H) were tested for in vitro NAT activity against substrate peptides with indicated N-termini a or against EEEIA-peptide for the indicated reaction times c. Sequences of the peptide substrates are EEEIAALRWGRPVGRRRRPVRVYP, DDDIAALRWGRPVGRRRRPVRVYP, MLGPEGGRWGRPVGRRRRPVRVYP, and SESSSKSRWGRPVGRRRRPVRVYP (bold indicating the variable N-termini of the peptides, while the rest of the peptide sequence is identical between them). b Structure of the human NAA10 protein (PDB ID: 6C9M) superimposed on the substrate from the Schizosaccharomyces pombe NAA10 structure (4KVM) Substrate peptide (SASE) is shown in green, while the mutation site Arg83 is shown in red. CoA is shown as a licorice model colored according to element. d Multiple sequence alignment of NAA10 from Homo sapiens, Mus musculus, Rattus norvegicus, Xenopus laevis, Danio rerio, and Saccharomyces cerevisiae. Conservation score is calculated by Jalview (http://www.jalview.org/), with * signifying perfect conservation

Structural conservation and surface charge of NAA10-R83H

A structural model was made by aligning the human NAA10 chain from the structure of the human NatA complex (PDB ID: 6C9M, [55]) to the corresponding chain of the S. pombe (PDB ID: 4KVM, [33]), the latter of which was solved with CoA and substrate peptide. This showed that the mutation site R83 is in proximity to the negatively charged phosphate group on the ribose ring of CoA (Fig. 2b). A defining feature of GNAT acetyltransferases is a core fold which includes an Ac-CoA binding region, four alpha-helices and six or seven beta-sheets [1]. R83 is part of the Ac-CoA binding region of NAA10, though R82 is even more conserved (Fig. 2d). We find an arginine in the same position in crystal structures of NAA20 [63], NAA40 [64] NAA50 [65], and NAA60 [66]. At physiological pH, arginine is expected to be protonated to a greater degree than histidine. We performed in silico mutagenesis in NAA10 from the hNatA structure [55] to determine whether the Ac-CoA binding region had altered electrostatic properties (Fig. 3), and found that the surface near Ac-CoA appears less basic when R83 is mutated into a histidine.

Fig. 3
figure3

NAA10-R83H has altered charge distribution at the surface facing Ac-CoA WT R83 a or H83 mutant NAA10 b is shown in a surface model, with residue 83 marked with an arrow. Charge was visualized in PyMOL. Blue denotes positive charge and red denotes negative charge. Ac-CoA is shown as a licorice model

Discussion

NAA10 is predicted to Nt-acetylate around 40% of the human proteome [1] as part of the NatA complex on the ribosome [29], but it may also be a monomeric NAT enzyme not associated with the ribosome [31]. In addition, both lysine acetyltransferase activities and acetyltransferase-independent functions have been proposed for NAA10 [34,35,36,37,38]. NAA10 is essential for normal development in all surveyed model organisms [34, 39,40,41]. Several mutations in human NAA10 are known, with substantial heterogeneity in presentation and severity of symptoms, yet with the common features of intellectual disability, delayed developmental and growth failure, and occasionally cardiac anomalies [42,43,44, 46,47,48,49, 67, 68]. In many cases, NAA10 mutations result in a decrease in catalytic activity (S37P, Y43S, I72T, R83C, V107F, V111G, R116W, and F128 L), while others destabilize NAA10 (Y43S, I72T, V111G, F128I, F128 L). The mutation causative of Ogden syndrome, S37P, leads to impaired NAA10 interaction with other NatA complex components [6], in addition to having decreased catalytic activity [42]. A new potential mechanism for how NAA10 mutation may cause developmental phenotypes is acetyltransferase-independent, however. Three previously described variants (S37P, V107F and R116W) were found to have decreased binding to imprinting control regions, potentially leading to dysregulation of genomic imprinting [34]. NAA15 mutations have also been described, with an overlapping phenotype to NAA10 mutations. Patients had several different indels or splice-site mutations, and presented with neurodevelopmental problems including intellectual disability, autism, motor function impairment and developmental delay [69]. Taken together, these findings suggest that NatA may have a critical role in nervous system development.

In this study, we present two male patients with a novel NAA10 variant, R83H. One patient presented with developmental delay, very limited language development at age 15, epilepsy, behavioural issues, but normal cardiac function. The other patient with a similar phenotype with developmental delay, very limited language development at age 12, ADHD like behaviour and hypertrophic cardiomyopathy. This is a similar clinical picture to what is known in some of the other described NAA10 patients with developmental delay and a cardiac phenotype.

Functional testing of in vitro catalytic activity was performed to determine whether the mutant enzyme was active. We found that NAA10-R83H had an unchanged substrate preference profile (Fig. 2a), but a greatly reduced catalytic activity. A time-course assay further bolstered this conclusion (Fig. 2c). The recurring R83C variant of NAA10 likewise had a greatly reduced catalytic activity [49], supporting the importance of R83 for NAA10 catalytic function.

We performed in silico mutagenesis to model the impact of this mutation on the surface charge of NAA10 (Fig. 3). While an arginine in this position contributes to a positively charged pocket, which can interact favorably with the phosphate groups on Ac-CoA, mutation to histidine decreases the positive charge density in this region. Several patients with a R83C variant at the same site have been reported. This variant likewise leads to a sharp decrease in catalytic activity [49]. The previously described R83C mutation and the loss of activity in R83H described here point to R83 being important for NAA10 enzymatic activity. The altered charge density in our in silico model suggests that perturbed interaction between NAA10-R83H and Ac-CoA may be the reason for the impaired catalytic activity. This is consistent with the hypothesis that the R83H variant leads to decreased NAA10 catalytic activity and that this is causing the symptoms observed in these two cases. Whether this is due to a loss of NAT- or KAT-activity, and if so which subset of substrates are mechanistically responsible for the phenotypes, is not clear.

Conclusions

We identify a novel missense variant of NAA10, and present evidence to support that this variant is causative of the symptoms seen in the two cases.

Abbreviations

Ac-CoA:

Acetyl coenzyme A

DD:

Development delay

DNMT1:

DNA methyltransferase 1

ID:

Intellectual disability

IP:

Immunoprecipitation

NAA10:

N-alpha acetyltransferase 10

NAA15:

N-alpha acetyltransferase 15

NAT:

N-terminal acetyltransferase

NatA:

N-terminal acetyltransferase A

PIX:

p21-activated kinase-interacting exchange factor

WES:

Whole exome sequencing

References

  1. 1.

    Aksnes H, Drazic A, Marie M, Arnesen T. First things first: vital protein Marks by N-terminal acetyltransferases. Trends Biochem Sci. 2016;41(9):746–60.

  2. 2.

    Vetting MW, et al. Structure and functions of the GNAT superfamily of acetyltransferases. Arch Biochem Biophys. 2005;433(1):212–26.

  3. 3.

    Hwang C-S, Shemorry A, Varshavsky A. N-terminal acetylation of cellular proteins creates specific degradation signals. Science. 2010;327(5968):973–7.

  4. 4.

    Shemorry A, Hwang C, Varshavsky A. Article control of protein quality and Stoichiometries by N-terminal acetylation and the N-end rule pathway. Mol Cell. 2013;50:1–12.

  5. 5.

    Park S-E, et al. Control of mammalian G protein signaling by N-terminal acetylation and the N-end rule pathway. Science. 2015;347(6227):1249–52.

  6. 6.

    Myklebust LM, et al. Biochemical and cellular analysis of Ogden syndrome reveals downstream Nt-acetylation defects. Hum Mol Genet. 2015;24(7):1956–76.

  7. 7.

    Xu F, et al. Two N-terminal acetyltransferases antagonistically regulate the stability of a nod-like receptor in Arabidopsis. Plant Cell. 2015;27:1547–62.

  8. 8.

    Holmes WM, Mannakee BK, Gutenkunst RN, Serio TR. Loss of amino-terminal acetylation suppresses a prion phenotype by modulating global protein folding. Nat Commun. 2014;5:4383.

  9. 9.

    Trexler AJ, Rhoades E. N-terminal acetylation is critical for forming alpha-helical oligomer of alpha-synuclein. Protein Sci. 2012;21(5):601–5.

  10. 10.

    Bartels T, Kim NC, Luth ES, Selkoe DJ. N-alpha-acetylation of α-synuclein increases its helical folding propensity, GM1 binding specificity and resistance to aggregation. PLoS One. 2014;9(7):e103727.

  11. 11.

    Arnesen T, et al. The chaperone-like protein HYPK acts together with NatA in Cotranslational N-terminal acetylation and prevention of huntingtin aggregation †. Mol Cell Biol. 2010;30(8):1898–909.

  12. 12.

    Scott DC, Monda JK, Bennett EJ, Harper JW, Schulman BA. N-terminal acetylation acts as an avidity enhancer within an interconnected multiprotein complex. Science. 2011;334(6056):674–8.

  13. 13.

    Monda JK, et al. Structural conservation of distinctive N-terminal acetylation-dependent interactions across a family of mammalian NEDD8 ligation enzymes. Structure. 2013;21(1):42–53.

  14. 14.

    Yang D, et al. Nα-acetylated Sir3 stabilizes the conformation of a nucleosome-binding loop in the BAH domain. Nat Struct Mol Biol. 2013;20(9):1116–8.

  15. 15.

    Arnaudo N, Fernández IS, Mclaughlin SH, Peak-chew SY, Rhodes D, Martino F. The N-terminal acetylation of Sir3 stabilizes its binding to the nucleosome core particle. Nat Struct Mol Biol. 2013;20(9):1119–21.

  16. 16.

    Setty SRG, Strochlic TI, Tong AHY, Boone C, Burd CG. Golgi targeting of ARF-like GTPase Arl3p requires its Nalpha-acetylation and the integral membrane protein Sys1p. Nat Cell Biol. 2004;6(5):414–9.

  17. 17.

    Behnia R, Panic B, Whyte JRC, Munro S. Targeting of the Arf-like GTPase Arl3p to the Golgi requires N-terminal acetylation and the membrane protein Sys1p. Nat Cell Biol. 2004;6(5):405–13.

  18. 18.

    Dikiy I, Eliezer D. N-terminal acetylation stabilizes N-terminal helicity in lipid- and micelle-bound α-Synuclein and increases its affinity for physiological membranes. J Biol Chem. 2014;289(6):3652–65.

  19. 19.

    Forte GMA, Pool MR, Stirling CJ. N-terminal acetylation inhibits protein targeting to the endoplasmic reticulum. PLoS Biol. 2011;9(5):e1001073.

  20. 20.

    Aksnes H, et al. An organellar nα-acetyltransferase, naa60, acetylates cytosolic n termini of transmembrane proteins and maintains golgi integrity. Cell Rep. 2015;10(8):1362–74.

  21. 21.

    Van Damme P, et al. NatF contributes to an evolutionary shift in protein N-terminal acetylation and is important for Normal chromosome segregation. PLoS Genet. Jul. 2011;7(7):e1002169.

  22. 22.

    Dinh TV, et al. Molecular identification and functional characterization of the first N-alpha-acetyltransferase in plastids by global acetylome profiling. Proteomics. 2015;15(14):2426–35.

  23. 23.

    Drazic A, et al. NAA80 is actin’s N-terminal acetyltransferase and regulates cytoskeleton assembly and cell motility. Proc Natl Acad Sci. 2018;115(17):4399–404.

  24. 24.

    Goris M, et al. Structural determinants and cellular environment define processed actin as the sole substrate of the N-terminal acetyltransferase NAA80. Proc Natl Acad Sci. 2018;115(17):4405–10.

  25. 25.

    Wiame E, et al. NAT6 acetylates the N-terminus of different forms of actin. FEBS J. 2018;285(17):3299–316.

  26. 26.

    Mullen JR, et al. Identification and characterization of genes and mutants for an N-terminal acetyltransferase from yeast. EMBO J. 1989;8(7):2067–75.

  27. 27.

    Park E, Szostak JW. ARD1 and NAT1 proteins form a complex that has N-terminal acetyltransferase activity. EMBO J. 1992;11(6):2087–93.

  28. 28.

    Arnesen T, Anderson D, Torsvik J, Halseth HB, Varhaug JE, Lillehaug JR. Cloning and characterization of hNAT5/hSAN: an evolutionarily conserved component of the NatA protein N-α-acetyltransferase complex. Gene. 2006;371(2):291–5.

  29. 29.

    Arnesen T, Anderson D, Baldersheim C, Lanotte M, Varhaug JE, Lillehaug JR. Identification and characterization of the human ARD1–NATH protein acetyltransferase complex. Biochem J. 2005;386(3):433–43.

  30. 30.

    Arnesen T, Van Damme P, Polevoda B, Helsens K, Evjenth R, Colaert N. Proteomics analyses reveal the evolutionary conservation and divergence of N-terminal acetyltransferases from yeast and humans. Proc Natl Acad Sci. 2009;106(20):8157–62.

  31. 31.

    Van Damme P, et al. Proteome-derived peptide libraries allow detailed analysis of the substrate specificities of Nα-acetyltransferases and point to hNaa10p as the post-translational actin Nα-acetyltransferase. Mol Cell Proteomics. 2011;10(5):M110.004580.

  32. 32.

    Gautschi M, et al. The yeast Nα-acetyltransferase NatA is quantitatively anchored to the ribosome and interacts with nascent polypeptides. Mol Cell Biol. 2003;23(20):7403–14.

  33. 33.

    Liszczak G, Goldberg JM, Foyn H, Petersson EJ, Arnesen T, Marmorstein R. Molecular basis for N-terminal acetylation by the heterodimeric NatA complex. Nat Struct Mol Biol. 2013;20(9):1098–105.

  34. 34.

    Lee CC, et al. The role of N-α-acetyltransferase 10 protein in DNA methylation and genomic imprinting. Mol Cell. 2017;68(1):89–103.e7.

  35. 35.

    Ree R, Varland S, Arnesen T. Spotlight on protein N-terminal acetylation. Exp Mol Med. 2018;50:90.

  36. 36.

    Kang J, Chun Y, Huh J, Park J. FIH permits NAA10 to catalyze the oxygen-dependent lysyl-acetylation of HIF-1α. Redox Biol. 2018;19(September):364–74.

  37. 37.

    Qian X, et al. Phosphoglycerate kinase 1 phosphorylates Beclin1 to induce autophagy. Mol Cell. 2017;65:917–31.

  38. 38.

    Lim JH, Chun YS, Park JW. Hypoxia-inducible factor-1α obstructs a wnt signaling pathway by inhibiting the hARD1-mediated activation of β-catenin. Cancer Res. 2008;68(13):5177–84.

  39. 39.

    Ree R, Myklebust LM, Thiel P, Foyn H, Fladmark KE, Arnesen T. The N-terminal acetyltransferase Naa10 is essential for zebrafish development. Biosci Rep. 2015;35:e00249.

  40. 40.

    Ingram AK, Cross GAM, Horn D. Genetic manipulation indicates that ARD1 is an essential Nα-acetyltransferase in Trypanosoma brucei. Mol Biochem Parasitol. 2000;111:309–17.

  41. 41.

    Wang Y, et al. Drosophila variable nurse cells encodes arrest defective 1 (ARD1), the catalytic subunit of the major N-terminal acetyltransferase complex. Dev Dyn. 2010;239:2813–27.

  42. 42.

    Rope AF, et al. Using VAAST to identify an X-linked disorder resulting in lethality in male infants due to N-terminal acetyltransferase deficiency. Am J Hum Genet. 2011;89:1–16.

  43. 43.

    Esmailpour T, et al. A splice donor mutation in NAA10 results in the dysregulation of the retinoic acid signalling pathway and causes Lenz microphthalmia syndrome. J Med Genet. 2014;51(3):185–96.

  44. 44.

    Popp B, et al. De novo missense mutations in the NAA10 gene cause severe non-syndromic developmental delay in males and females. Eur J Hum Genet. 2015;23(5):602–9.

  45. 45.

    Valentine V, Sogawa Y, Rajan D, Ortiz D. A case of de novo NAA10 mutation presenting with eyelid myoclonias (AKA Jeavons syndrome). Seizure. 2018;60:120–2.

  46. 46.

    Casey JP, et al. NAA10 mutation causing a novel intellectual disability syndrome with long QT due to N-terminal acetyltransferase impairment. Sci Rep. 2015;5:16022.

  47. 47.

    Støve SI, et al. A novel NAA10 variant with impaired acetyltransferase activity causes developmental delay, intellectual disability, and hypertrophic cardiomyopathy. Eur J Hum Genet. 2018;26(9):1294–305.

  48. 48.

    McTiernan N, et al. NAA10 dysfunction with normal NatA-complex activity in a girl with non-syndromic ID and a de novo NAA10 p.(V111G) variant - a case report. BMC Med Genet. 2018;19:47.

  49. 49.

    Saunier C, et al. Expanding the phenotype associated with NAA10-related N-terminal acetylation deficiency. Hum Mutat. 2016;37(8):755–64.

  50. 50.

    The Deciphering Developmental Disorders Study. Large-scale discovery of novel genetic causes of developmental disorders. Nature. 2015;519(7542):223–8.

  51. 51.

    Li H. “Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM,” arXiv; 2013. p. 1303.3997v2 [q-bio.GN].

  52. 52.

    G. A. Van der Auwera et al., From fastQ data to high-confidence variant calls: The genome analysis toolkit best practices pipeline, no. SUPL.43. 2013.

  53. 53.

    Larkin MA, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.

  54. 54.

    Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ. Jalview version 2-a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25(9):1189–91.

  55. 55.

    Gottlieb L, Marmorstein R. Structure of human NatA and its regulation by the huntingtin interacting protein HYPK article. Structure. 2018;26:1–11.

  56. 56.

    Schrödinger LLC. The PyMOL molecular graphics system, version 1.8. 2015.

  57. 57.

    Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA. PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res. 2004;32:665–7.

  58. 58.

    Dolinsky TJ, et al. PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations. Nucleic Acids Res. May 2007;35:W522–5.

  59. 59.

    Jurrus E, et al. Improvements to the APBS biomolecular solvation software suite. Protein Sci. 2018;27:112–28.

  60. 60.

    Weiner SJ, et al. A new force field for molecular mechanical simulation of nucleic acids and proteins. J Am Chem Soc. 1984;106(3):765–84.

  61. 61.

    Søndergaard CR, Olsson MHM, Rostkowski M, Jensen JH. Improved treatment of ligands and coupling effects in empirical calculation and rationalization of p Kavalues. J Chem Theory Comput. 2011;7(7):2284–95.

  62. 62.

    P. H. Network, “How do I calculate Echo Z-scores?” [Online]. Available: http://www.pediatricheartnetwork.org/ForResearchers/EchoZ-Scores.aspx.

  63. 63.

    Hong H, Cai Y, Zhang S, Ding H, Wang H, Han A. Molecular basis of substrate specific acetylation by N-terminal acetyltransferase NatB. Structure. 2017;25(4):641–9.

  64. 64.

    Magin RS, Liszczak GP, Marmorstein R. The molecular basis for histone H4- and H2A-specific amino-terminal acetylation by NatD. Structure. 2015;23(2):332–41.

  65. 65.

    Liszczak G, Arnesen T, Marmorsteins R. Structure of a ternary Naa50p (NAT5/SAN) N-terminal acetyltransferase complex reveals the molecular basis for substrate-specific acetylation. J Biol Chem. 2011;286(42):37002–10.

  66. 66.

    Støve SI, Magin RS, Foyn H, Haug BE, Marmorstein R, Arnesen T. Crystal structure of the Golgi-associated human Nα-acetyltransferase 60 reveals the molecular determinants for substrate-specific acetylation. Structure. 2016;24:1–13.

  67. 67.

    Rauch A, et al. Range of genetic mutations associated with severe non-syndromic sporadic intellectual disability: an exome sequencing study. Lancet. 2012;380(9854):1674–82.

  68. 68.

    Thevenon J, et al. Diagnostic odyssey in severe neurodevelopmental disorders: toward clinical whole-exome sequencing as a first-line diagnostic test. Clin Genet. 2016;89(6):700–7.

  69. 69.

    Cheng H, et al. Truncating variants in NAA15 are associated with variable levels of intellectual disability, autism Spectrum disorder, and congenital anomalies. Am J Hum Genet. 2018;102(5):985–94.

Download references

Acknowledgements

We are very grateful for the participation of the families involved in this study. The DDD study presents independent research commissioned by the Health Innovation Challenge Fund (grant number HICF-1009-003). See [50] or www.ddduk.org/access.html for full acknowledgement.

Funding

The work has been supported by the Research Council of Norway (project 249,843), the Norwegian health authorities of Western Norway (project 912,176), the Norwegian Cancer Society, and Bergen Research Foundation (BFS). The funding bodies had no role in the design of the study, collection, analysis, or interpretation of data or in writing the manuscript

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files. Generated plasmids are available from the corresponding author on request. Information on the NAA10 missense variant c.248G > A, p.(R83H) has been submitted to Decipher (DDD-286692).

Author information

RR and ASG designed and performed experiments, analysed results, made figures and wrote the manuscript. SAL, PMT, KPS, MD consulted the patients, designed and analysed WES data, led the study. TA designed and analysed experiments, led the study. All authors read, approved and commented on the final manuscript.

Correspondence to Sally Ann Lynch or Thomas Arnesen.

Ethics declarations

Ethics approval and consent to participate

All parents have consented in writing to publication of these findings. Since the findings were a consequence of routine clinical evaluation and diagnostics, and further research did not require patient investigations that would not otherwise have been performed, ethical review board approvals are not required according to national rules.

Consent for publication

Parents of both children gave written consent to publish this case report, including the medical data.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Video S1. Transthoracic echocardiography showing the left ventricle from an apical 3-chamber view (identical to Fig. 1) during two heart cycles in color compared mode. Please note that although there is no obstruction to flow in the left ventricular outflow tract (LVOT), a slight turbulence does occur along where the hypertrophic septum protrudes into the LVOT. (MOV 793 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • NAA10
  • X-linked
  • XLID
  • Microcephaly
  • Intellectual disability
  • N-alpha-acetyltransferase
  • Acetylation
  • NatA