Abstract
NAA15 is a member of the NatA N-terminal acetyltransferase complex which acts by interacting with the NAA10 enzymatic sub-unit. Individuals with variants in the NAA15 coding region develop NAA15-related neurodevelopmental syndrome which presents with a wide array of manifestations that affect the heart, brain, musculoskeletal system, and behavioral and cognitive development. We tracked a cohort of 26 participants (8 females and 18 males) over time, each with a pathogenic NAA15 variant, and administered the Vineland-3 assessment to them to assess their adaptive functioning. We found that the cohort performed significantly worse compared to the normalized Vineland values. On average, females performed better than males across all domains. They performed significantly better on the Motor Domain and Fine Motor Sub-Domain portions of the assessment. Over time, females showed a decrease in adaptive functioning with the decline being especially strongly correlated at the Coping sub-domain, Domestic sub-domain, and Fine motor sub-domains. It is difficult to determine the strength of these correlations due to limited power. Males (after excluding one outlier) showed a moderate positive correlation between age and ABC standard score. Ultimately, additional longitudinal data should be collected to determine the validity of the between sex-differences and to better understand the change in adaptive behavioral outcomes of individuals with NAA15-neurodevelopmental disorder as they age.
Introduction
N-terminal acetylation is a prevalent protein modification that involves the addition of an acetyl group to the alpha-amino group of the N-terminal amino acid of a protein (Liszczak et al., 2013). This enzymatic process is catalyzed by a group of N-terminal acetyltransferases (NATs) and plays a pivotal role in various cellular processes, including protein stability, protein-protein interactions, and cellular localization (Aksnes et al., 2015; Bienvenut et al., 2012; Dikiy & Eliezer, 2014). There are seven known NAT types, with NatA being the most prominent, believed to be responsible for modifying 40-50% of the human proteome(Deng et al., 2019; Feng & Ma, 2016; Van Damme, 2021). The NatA complex is composed of the catalytic sub-unit NAA10, and the auxiliary sub-units NAA15 and HYPK (Deng et al., 2019; Dörfel et al., 2017; Feng et al., 2016; Knorr et al., 2019; Weidenhausen et al., 2021).
Dysregulation of NAA10 and NAA15 expression or activity has been associated with the pathogenesis of various diseases, including cancer (Kuhns et al., 2018; Le et al., 2023; Zhang et al., 2020; Zhu et al., 2024), neurodevelopmental disorders (Cheng et al., 2018; Lyon et al., 2023; Makwana et al., 2024; Sidhu et al., 2017; Wu & Lyon, 2018), and cardiovascular diseases (Belbachir et al., 2023; Huth et al., 2023; Ritter et al., 2021; Støve et al., 2018; Ward et al., 2021). More notably, mutations in the X-linked gene NAA10 lead to Ogden Syndrome, which has numerous effects on individual’s development and overall health. NAA15-related neurodevelopmental syndrome has overlapping symptoms with Ogden Syndrome, including intellectual disability, atypical or delayed development, heart anomalies, and more features previously described (Cheng et al., 2020; Dörfel & Lyon, 2015; Lyon et al., 2023; McTiernan et al., 2018; Myklebust et al., 2015; Patel et al., 2024; Pesz et al., 2018; Stessman et al., 2017; Støve et al., 2018; Umlai et al., 2021; Wong et al., 2019; Zhao et al., 2018). Additionally, NAA15 variants are associated with musculoskeletal and neuromuscular abnormalities that can appear in childhood or may have a delayed appearance until adulthood (Danti et al., 2024; Monestier et al., 2019; Pesz et al., 2018; Tian et al., 2022; Yubero et al., 2022).
The aim of the current study is to further describe adaptive functioning in individuals with NAA15-related neurodevelopmental syndrome by analyzing prospective data and understand how this disorder progresses over time.
Methods
Participants
All families that participated in the investigation signed Institutional Review Board (IRB)-approved consents and HIPAA forms. The scope of the project was explained in lay language. The participant population included children with NAA15 pathogenic variants that the principal investigator has previously worked with for other research projects. Participants were not compensated for their time. Overall, there were 26 participants who were administered the assessment.
Cognitive Assessment
The participants were administered Vineland-3 which is composed of three core domains: Communication Daily Living Skills, and Socialization. These domains are further split into sub-domains that are used to assess an individual’s competency in tasks of personal and social sufficiency (Perry et al., 2009). The scores for each core domain are norm referenced and summed to generate an Adaptive Behavior composite (ABC) score that provides a holistic picture of an individual’s adaptive behavior across domains. The Vineland 3 tools were administered by three trained assessors at various timepoints in the participants’ lives ranging from 1-23 years old. Individual assessor scores were cross validated with one another to ensure equivalency in how assessments were administered and scored. Participant caregivers were those interviewed due to their proximity and knowledge of the patients. Participants were given access to the Vineland assessment scores after they had been scored.
Analysis
The Principal Investigator and associated staff collected at least one Vineland 3 score for each participant over multiple years. Some participants were able to take the assessment more than once providing a longitudinal outlook. Natural history analysis and visualization was then performed using the GraphPad Prism software. All pathogenic variants were analyzed together and graphed by comparing their ABC standard score against the age at which the assessment was administered. Further analysis was performed by filtering graphs by sex. All pathogenic variants, male participants, and female participants were also graphed according to their Communication (com), Daily Living Skills (dls), Social (soc), and Motor (mot) domain standard scores. The domain scores were further broken down by sub-domain standard score to allow for more granular identification of participant strength and weakness over time. The Communication sub-domain scores were Receptive (rec), Expressive (exp), and Written (wrn). The Daily Living Skills sub-domains analyzed were Personal (per), Domestic (dom), and Community (cmm). The Socialization sub-domains of interest were Interpersonal Relationships (ipr), Play and Leisure (pla), and Coping Skills (cop). The Motor sub-domains of Interest were Fine Motor (fmo) and Gross Motor (gmo). Scores were only collected and included for this motor domain if patients were under the age of nine to stay within the validated age range of the exam (Perry et al., 2009). Lastly, the Internalizing (int) and Externalizing (ext) sub-domain components of the Maladaptive Trait domain were also included in analysis, however the Critical Items component of the domain were not.
The scores were compared in Microsoft Excel with two-tailed equal variance t-tests to compare the adaptive behavior outcomes between males and females to determine if there were sex differences. The p-value was set at .05 for all calculations. Correlation coefficients were also calculated for each of the graphs made using Excel. A second round of calculations for the t-test and the correlation coefficients in the male cohort was performed to exclude NAA15-011 due to being three standard deviations above the mean in age compared to the other participants and having performed two standard deviations below the mean in Vineland ABC, Motor, Social, and Daily Living Skills standard scores. Data presented in tables includes calculations performed including participant NAA15-011 unless otherwise specified.
Results
There were 26 participants who were administered the Vineland and included in this study. A breakdown of the different pathogenic variants in each sex can be seen in Table 1. Males and females each had unique pathogenic variants. Twenty-four participants developed their variant de novo and 2 participants have variants of unknown origin (due to the parents not being tested). Of the 26 participants, 8 are female and 18 are male with the age range being from ∼0 to ∼25 years old (mean = 8.8 years, standard deviation (SD) = 4.6).
Compared to the general population, individuals with NAA15 variants scored below average on the Vineland 3 Assessment (where average is mean=100, sd=15). This was true for ABC standard scores as well as the main domain scores. The average ABC score among NAA15 variants was 66.7 (SD = 21.1). Females tended to outperform males across adaptive behavior domains. However, except for the Motor domain score these differences were not significant. A summary of the Vineland 3 scores can be seen in Table 2.
Graphical representations of the ABC standard scores over time can be seen in Figure 1. Before stratifying by sex, there appears to be an overall downward trend in score with age. Figure 2 shows the ABC standard scores separated by male and female. Separation by sex shows a moderate linear increase in adaptive function in males over time, when excluding NAA15-011, whereas females appear to show a moderately strong linear decay. There is also a strong linear decay in the female Daily Living Skills domain over times whereas the males showed a moderate positive linear growth when excluding NAA15-011. Table 3 showcases the remaining calculated correlation coefficient values for the Vineland ABC and Domain standard scores for the females, males, and males excluding NAA15-011.
The Communication, Daily Living Skills, Social, and Motor standard score graphs can be seen in Figure 3 and follow similar trends for females and males, when excluding NAA15-011, as the ABC scores over time. It should be noted that excluding NAA15-011 did not significantly change the t-test statistics of the standard scores when comparing the males and the females. The magnitude at which the females outperformed the males decreased. However, except for the Motor domain score, their differences were still not significant.
Dividing the Motor domain into its Fine Motor and Gross Motor sub-domains elucidated the difference between sexes. Males, with and without the inclusion of NAA15-011, on average scored 7.0 (SD = 2.5) which is significantly less than the females who had an average score of 14.0 (SD = 3.6). On average for Gross Motor subdomain scores, Males (mean = 8.5, SD = 3.5) also performed worse than their female counterparts (mean = 10.8, SD = 1.8). This difference was not significant. There were no other significant differences between how the sexes performed in the various sub-domains. A summary of the breakdown of each domain into its sub-domain components can be seen in Table 4.
Each domain sub-domain score can be seen in Figure 4. There are strong negative linear associations between the Coping, Fine Motor, and Domestic sub-domains over time in the NAA15 females. There are also moderately strong negative linear associations in the Receptive and Personal sub domains. The correlation coefficients in the males were not as strong as those in the females. However, when excluding NAA15-011, the Personal sub-domain had a moderate positive linear relationship over time. Calculated correlation coefficients for each sub-domain can be seen in Table 5.
Discussion
Individuals with NAA15-related neurodevelopmental syndrome performed significantly worse than the mean on the Vineland 3. This difference was expected due to the expression of the NatA complex being ubiquitous across tissue types and playing a prominent role in development (Cheng et al., 2020; Liszczak et al., 2013; Lyon et al., 2023). When comparing performance by sex, the females on average performed better than the males. Except for two participants with an unknown inheritance pattern, the NAA15 pathogenic variants in the study cohort occurred de-novo. Without evidence that the disease is maternally inherited, the difference in sex is unexpected. A possible explanation could be that the NAA10 catalytic subunit of the NatA complex is present on the X-chromosome (Bader et al., 2020; Lyon et al., 2023; Rope et al., 2011; Wu & Lyon, 2018). An increase in NAA10 dosage due to partial inactivation or a different X-chromosome escape mechanism could partially offset the loss of active site conformational change of NAA10 by NAA15 (Carrel & Willard, 2005; Cotton et al., 2013; Schultz et al., 2015; Tukiainen et al., 2017). The mechanism of NAA15-related neurodevelopmental syndrome appearing to be a milder form of NAA10-related neurodevelopmental syndrome (Ogden Syndrome) also provides some support for this hypothesis as the NAA10 catalytic unit has function outside of the NatA complex but is unable to acetylate NatA specific sites without NAA15 (Cheng et al., 2020; Lyon et al., 2023; Makwana et al., 2024; Patel et al., 2024; Wu & Lyon, 2018). However, the differences seen between sexes could also be due to small sample sizes. Also, because of a lack of overlap between pathogenic variants in males and females, it was not possible to do a direct comparison in function while controlling for sex. Further genetic analysis and cognitive testing needs to be performed to determine if the differences between sexes are significant.
While females performed on average better than males across all adaptive domains, they performed significantly better on the motor domain and fine motor sub-domain portions of the Vineland. Pathogenic variants in NAA15 have been associated with infant onset dystonia in a few individuals (Cheng et al., 2018; Danti et al., 2024; Lyon et al., 2023; Yubero et al., 2022) which could be one reason why the cohort performed worse on motor tasks compared to the norm population. However, there did not appear to be a difference between the sexes in the referenced literature in how musculoskeletal manifestations presented. Furthermore, there is a case where a NAA15 variant led to the adult-onset presentation of parkinsonism in a man presenting with other manifestations of NAA15-related neurodevelopmental syndrome such as speech delay and intellectual disability (Straka et al., 2022). Continued genetic testing of individuals with adult-onset dystonia and longitudinal monitoring of the current cohort for motor disorders will be necessary to help determine if the difference in motor function between male and females is due to chance, sex differences, or pathogenic variant type.
In aggregate, it is difficult to determine if there is an overall increase or decrease in function over time in the participants. However, when comparing performance by sex, there is a strong negative correlation between adaptive behavior in females and age. The decline with age has the strongest correlation between the coping, domestic, and fine motor sub-domains. With a decline in fine motor capacity, it would stand to reason that an individual’s ability to take care of oneself (domestic) would also decline. However, both the domestic and fine motor sub-domains have age restrictions leading to unequal administration of those portions of the exam creating an even smaller population size to compare from (Farmer et al., 2020; Perry et al., 2009, p. 3).
Furthermore, this decline in function contradicts the mild recovery in gross motor, personal-social, and language development function previously reported in the literature (Cheng et al., 2018, 2020; Tian et al., 2022). The decline in the coping sub-domain in the females could also be partially explained by the decrease in motor skills as there is evidence that similar motor declines present in those with Autism Spectrum Disorder (ASD) have led to decreased coping and socialization skills overall (Peyre et al., 2024; Raditha et al., 2023). The overall decline could also be partially explained by the upper half of the age range being from 11-17, a time period that is traditionally tumultuous for adolescents due to various physical and hormonal changes, and is associated with increased stress for those with ASD (Cridland et al., 2014; Steward et al., 2018; Tager-Flusberg & Kasari, 2013). Ultimately, follow up studies need to be performed on this cohort to determine if they are continuing to show decline over time, if there is a bottoming out effect, or if they begin to improve after navigating adolescence.
In contrast to the female decline over time, males when analyzed without NAA15-011 are correlated to have a moderate increase in adaptive function with age. However, there are not strong correlations present in any subsequent sub-domain analysis performed. Additionally, when including NAA15-011, there were no or small correlations between Vineland score and age.
Patient NAA15-011 was removed from the analysis due to their having performed much more poorly than the rest of their cohort. They had a de novo p.Asn864Ser pathogenic variant which is located in part of the region (residues 500-866 (Arnesen et al., 2010)) that interacts with the HYPK chaperone like protein which is important for NatA mediated acetylation and Huntingtin protein aggregation (Arnesen et al., 2010; Ghosh & Ranjan, 2022; Miklánková et al., 2022; Weyer et al., 2017). In addition to dysfunctional N-terminal acetylation being associated with neurodevelopmental delay, abnormal Huntingtin protein aggregation also can lead to cognitive decline (Alberti & Hyman, 2021; Liu et al., 2023; Soltani Khaboushan et al., 2023).
The stark difference in function between NAA15-011 and the rest of the male cohort could be explained by a substantial decrease in HYPK interaction. However, additional molecular testing would need to be conducted to be certain due to others in the cohort having mutations in this cohort without having similar decreases in function. If abnormally aggregated Huntingtin protein is associated with the steeper decline in function in patient NAA15-011, then future potential therapies that are effective in decreasing Huntingtin aggregation could possible rescue partial function in these individuals (Im et al., 2015; Jain & Roy, 2023; Posey et al., 2018; Roth et al., 2023). The increase in age could also be the reason for the sharp decline in function as NAA15-011 is 10 years older than the second oldest male. However, additional longitudinal data needs to be collected before definitive conclusions can be drawn.
The overall decrease in female adaptive behavior scores and the overall significantly decreased male adaptive scores suggest that there is an urgent need for additional study of and the development of interventions for NAA15-related neurodevelopmental disorder. Future research should aim to complete the clinical timeline as best as possible. The present cohort was only able to include individuals up to the age of ∼20-24 years due to the rarity of the disease.
Without insight into the later stages of the disorder, it will be more difficult for clinicians to help provide expectations of what to expect at different life stages to caregivers. NAA15-related neurodevelopmental disorder is thought to be caused by haploinsufficiency (Cheng et al., 2018; Tian et al., 2022) of which there are several therapies in development that show promise in cell lines and animal models that have been used in similar disorders such as Rett Syndrome and Angelman Syndrome to help restore protein function and dosage to adequate levels (Albadri et al., 2017; Burbano et al., 2022; Hill & Meisler, 2021; Milazzo et al., 2021; Palmieri et al., 2023; Protic et al., 2019). If interventions are implemented early enough, it may be possible to restore protein function and improve adaptive behavior skills.
Conclusion
NAA15-related neurodevelopmental disorder is closely related to Ogden Syndrome (NAA10-related neurodevelopmental disorder) due to both the NAA15 and NAA10 proteins being part of the NatA complex. NAA15 plays a supporting role in complex with the NAA10 enzymatic subunit, with dysfunction leading to clinical manifestations in the heart, musculoskeletal system, and brain. The pathogenic variants in NAA15-related neurodevelopmental disorder are heterozygous, so there is still one functioning copy of NAA15, meaning that the NatA complex can still form, but likely with haploinsufficiency. As such, NAA15 pathogenic variants usually present less severely than those persons with pathogenic NAA10 variants. Repeated Vineland-3 administration to a cohort of 26 patients with unique de novo NAA15 pathogenic variants has revealed the natural history of the disease, showcasing an overall decrease in adaptive functioning compared to the norm from infancy till adolescence.
Females tended to perform on average better than males. However, over time, females tended to show a decline in function whereas males showed an increase. Additional natural history studies need to be performed to understand disease outcomes outside of adolescence, difference between sexes, and to identify outliers who perform significantly better or worse than their peers so genetic or environmental factors influencing disease can be identified.
Data Availability
All data produced in the present study are available upon reasonable request to the authors.
Ethical Approval
Both oral and written patient consent were obtained for research and publication, with approval of protocol #7659 for the Jervis Clinic by the New York State Psychiatric Institute - Columbia University Department of Psychiatry Institutional Review Board.
Funding
This work is supported by New York State Office for People with Developmental Disabilities (OPWDD) and NIH NIGMS R35-GM-133408.
Competing Interests
The authors declare that they have no competing interests or personal relationships that could have influenced the work reported in this paper.