A rare missense variant in the ATP2C2 gene is associated with language impairment and related measures
======================================================================================================

* Angela Martinelli
* Mabel Rice
* Joel B. Talcott
* Rebeca Diaz
* Shelley Smith
* Muhammad Hashim Raza
* Margaret J. Snowling
* Charles Hulme
* John Stein
* Marianna E. Hayiou-Thomas
* Ziarih Hawi
* Lindsey Kent
* Samantha J. Pitt
* Dianne F. Newbury
* Silvia Paracchini

## ABSTRACT

At least 5% of children present unexpected difficulties in expressing and understanding spoken language. This condition is highly heritable and often co-occurs with other neurodevelopmental disorders such as dyslexia and ADHD. Through an exome sequencing analysis, we identified a rare missense variant (chr16:84405221, GRCh38.p12) in the *ATP2C2* gene. *ATP2C2* was implicated in language disorders by linkage and association studies, and exactly the same variant was reported previously in a different exome sequencing study for language impairment (LI). We followed up this finding by genotyping the mutation in cohorts selected for LI and comorbid disorders. We found that the variant had a higher frequency in LI cases (1.8%, N=360) compared to cohorts selected for dyslexia (0.8%, N = 520) and ADHD (0.7%, N = 150), which presented frequencies comparable to reference databases (0.9%, N = 24,046 gnomAD controls). Additionally, we observed that carriers of the rare variant identified from a general population cohort (N=42, ALSPAC cohort) presented, as a group, lower scores on a range of reading and language-related measures compared to controls (N=1825) (minimum p = 0.002 for nonword reading). *ATP2C2* encodes for an ATPase (SPCA2) that transports calcium and manganese ions into the Golgi lumen. Our functional characterization suggested that the rare variant influences the ATPase activity of SPCA2. Thus, our results further support the role of *ATP2C2* locus in language-related phenotypes and pinpoint the possible effects of a specific rare variant at molecular level.

Keywords
*   ATP2C2
*   SPCA2
*   language impairment
*   DLD
*   ALSPAC
*   neurodevelopment

## INTRODUCTION

Speech and language disorders represent one of the most common childhood disabilities and account for a large proportion (23%) of paediatric referrals and statements of special educational need (School census, England, January 2019, [https://www.gov.uk/government/statistics/special-educational-needs-in-england-january-2019](https://www.gov.uk/government/statistics/special-educational-needs-in-england-january-2019)).

Speech and language deficits are often the secondary manifestations of a pre-existent neurological or medical condition (e.g. Down syndrome, intellectual disability, autism spectrum disorder, hearing loss, traumatic brain injury and mental retardation). Severe and persistent language difficulties can be observed also despite normal intelligence and access to adequate educational opportunities. In this case children were diagnosed with specific language impairment (SLI) (1). Recently, the use of this term has been replaced by the developmental language disorder (DLD) definition (2,3) which captures a wider group of children presenting co-occurring conditions such as intellectual disability, ADHD and developmental coordination disorder (DCD). Previous genetic research, and the samples described here, refer to participants recruited under the SLI criteria and therefore we will use the terms SLI and LI hereafter.

SLI affects up to 8% of pre-school children (4,5) and has a significant impact on social, academic and psychological outcomes (6–10). SLI often co-occurs with other neurodevelopmental disorders such as dyslexia (11,12) and attention deficit/hyperactivity disorder (ADHD) (13). The comorbidity observed between these conditions suggests that SLI, dyslexia and ADHD may be different manifestations of the same cognitive deficit or share common genetic factors (14). This hypothesis is supported by significant genetic overlap across studies as demonstrated by a recent genome-wide association study (GWAS) for dyslexia which showed significant shared contributions across disorders (15).

A strong genetic component, of ∼ 70%, was first proposed by twin and familial studies (16–20). Initially, linkage studies conducted across multiplex families led to the identification of 5 susceptibility loci on chromosomes 2q36 (SLI5, OMIM #615432) (21), 7q35-q36 (SLI4—OMIM #612514) (22), 13q21 (SLI3—OMIM #607134) (23,24), 16q (SLI1—OMIM #606711) and 19q (SLI2— OMIM #606712) (25–28). Fine mapping within these regions reported associations of common variants at the *ATP2C2, CMIP, CNTNAP2* genes with quantitative measures of language abilities (22,29), including the phonological memory nonword repetition task (NWR), which is considered as a strong endophenotype for SLI (30).

GWAS for language skills have been conducted in clinical and epidemiological samples, but have been limited to a few studies of relatively small sample sizes (31). In the first case, the only locus reaching genome-wide significance in SLI individuals was found in the *NOP9* gene, when including parent-of-origin effects (32,33). In the second case, the most significant result was the association between the locus 3p12.3, near *ROBO2*, and expressive vocabulary in the early phase of language acquisition (34). A GWAS conducted in two population cohorts reported an association close to statistical significance between *ABCC13* gene and the NWR test (35). Other studies combined analyses for language and reading abilities using both a case-control (36) and quantitative design (37) reporting only marginally significant associations.

Chromosomal rearrangements and copy number variants (CNV) have also been implicated in the aetiology of language impairment. Children with SLI show a higher frequency of sex chromosome aneuploidies (38), and an increased CNV burden compared to population controls (39,40). A balanced translocation interrupting the gene *SEMA6D* (41) has been identified in a child with SLI. Rare deletions have been reported in the *ZNF277* (42) and *TM4SF20* (SLI5) genes (21). Larger *de novo* deletions have been reported across the *ATP2C2* gene in a 10 year old boy with receptive and expressive language delay (43) and at chromosome 15q13.1–13.3 in a child with LI (44). Single locus mutations have been implicated in severe forms of language disorders (31,45). De novo mutations in the *FOXP2* gene cause childhood apraxia of speech (CAS), a severe oral-motor delay accompanied by deficits in expressive and receptive language (46).

Next-generation sequencing (NGS) technology is uncovering the role of rare variants in complex traits. NGS in individuals with CAS or SLI pinpointed to novel pathogenic variants in coding and non-coding regions (45,47,48). Whole exome sequencing (WES) analysis in 43 unrelated probands affected by severe SLI (49) reported several risk variants but suggested that a single mutation might not be sufficient to lead to a disorder unless in combination with other factors in line with previous studies (50–52).

Here, we followed up a specific finding from ongoing exome sequencing analyses on large pedigrees selected for having multiple members affected by SLI. We identified, in one family, a rare non-synonymous variant in the *ATP2C2* gene (Chr16:84405221 GRCh38.p12, rs78887288, NM_001286527:exon3:c.G304A:p.V102M), previously implicated in SLI through association studies (29,43)(29,43). Exactly the same variant, was also reported in the WES study by Chen and collaborators (49) bolstering the idea that this rare variant has a role in susceptibility to language disorders. Nonetheless, given its rarity, larger screens and functional studies are required to support this role.

We screened the chr16:84405221 variant in larger cohorts of children with language difficulties. Due to the high rate of comorbidity between SLI, dyslexia and ADHD, the screening included cohorts selected for these conditions.

We found that this variant had a higher frequency in the SLI and LI cases compared to the control group and cohorts selected for dyslexia or ADHD. Additionally, this variant was associated with low performance in a range of reading and language measures. By means of a colorimetric ATPase assay, we further investigated the effects of the variant on the function of SPCA2, the calcium and manganese transporter encoded by *ATP2C2*. The mutation displayed a lower, albeit not significant, ATPase activity compared to SPCA2 WT.

## MATERIAL AND METHODS

### Discovery pedigree

An exome sequencing analysis was conducted in six pedigrees of SLI probands with multiple affected members per family, selected from an ongoing longitudinal study in the Midwestern United States (53,54). The full results of this screening will be reported in a separate manuscript. SLI probands were identified based on language performance 1 standard deviation (SD) or more below the mean on an age appropriate omnibus language test (53), and met the following entrance screening criteria: (i) a standard score above 85 on an initial age appropriate test of non-verbal IQ, (ii) normal hearing acuity, (iii) no history of neurological disorders or diagnosis of autism, and (iv) intelligible speech sufficient for language transcription. Except for normal hearing acuity, the entrance screening criteria did not apply to siblings and parents of SLI probands, as their intellectual status was an outcome of the study (53).

Here we report the analysis and follow up of a specific variant of interest (Chr16: 84405221) in one family (Figure 1). A male SLI proband, his younger male sibling (Sib 1) and their parents were selected for WES. All family members including the other two siblings were genotyped for the *ATP2C2* variant. Longitudinal phenotype data on the proband were assessed from age 6 to 19 years and the siblings were assessed also at multiple time points. Both mother and father were tested at a single time point. All siblings and the father met the criteria above for a clinical diagnosis of SLI. The sibling selected for WES was also delayed on the Test of Early Grammatical Impairment (TEGI) until age 7 years, 6 months (55). The children were also delayed in receptive vocabulary and oral reading across age, as well as in mean length of utterance when young. The study was approved by the institutional review boards at the University of Kansas and at the University of Nebraska Medical Center. Appropriate informed consent was obtained from the subjects and/or their caregivers.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.12.21249315/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/01/15/2021.01.12.21249315/F1)

Figure 1. The rare *ATP2C2* missense variant (Chr16:84405221).
**A**) Discovery pedigree. Father (1), mother (2), proband (3) and siblings (4, 5, 6). Subjects with SLI are shaded black and unaffected are white. Individuals included in the WES screening are indicated by an asterisk. **B**) Allele-specific PCR. The presence of the missense variant was validated in follow up screening using PCR 3’-end mismatch primers. Depending on the allele present, the forward primer anneals perfectly to the template or forms a mismatch. The reaction gives an amplicon only in the first case. **C**) On top, a schematic of the *ATP2C2* gene. Exons are represented by numbered boxes. The rare variant is localised in the third exon (red arrow) and causes the substitution of the residue from valine to methionine. Exons and introns of *ATP2C2* gene are not to scale. At the bottom, predicted domains of the protein ([http://www.rcsb.org/pdb/protein/O75185](http://www.rcsb.org/pdb/protein/O75185)) (82). The variant changes the amino acid sequence which codes for the cation transporting ATPase domain at the N-terminal of the protein.

### Validation, genotyping and bioinformatics prediction for the variant

Sanger sequencing was used to validate the variant and test all available family members in the discovery pedigree.

The variant was then screened in larger cohorts with a TaqMan™ SNP Genotyping Assay run on a ViiA7 instrument (Life Technologies, Paisley, UK). The presence of the mutation was validated using PCR 3’-end mismatch primers (Figure 1B). All primer sequences are reported in Supplementary Table S1. For protocol details see Supplementary Table S2 and S3, and Supplementary Figure S1. The functional impact of the rare *ATP2C2* variant was estimated using the default settings of seven different *in silico* prediction algorithms (Table 1), SIFT (56), Mutation Assessor (57), PROVEAN (58), PolyPhen-2 (59), MutPred (60), INPS-MD (61), MutationTaster (62).

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/01/15/2021.01.12.21249315/T1)

Table 1. Prediction of the possible effects of the *ATP2C2* missense variant (chr16:84405221).

### Cohorts for follow-up analysis

The genotyping was conducted across different UK cohorts of children with language impairment, dyslexia or ADHD used in previous genetic studies (Table 2). We use the SLI term when participants were recruited under SLI criteria and the term LI to refer to individuals who did not receive a specific diagnosis but presented language difficulties on psychometric tests. Assignment to a case group was based either on existing clinical diagnoses (York and ADHD cohorts) or derived from poor performance on different psychometric tests as specified below. Exclusion criteria included non-European ethnicity, English as a second language, signs of other sensory or neurological conditions and missing information for the selected phenotypes.

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/01/15/2021.01.12.21249315/T2)

Table 2. Frequency of the rare allele A (Chr16:84405221) across different cohorts.

The SLI Consortium (SLIC) cohort consisted of 263 families (263 probands, 258 siblings, 551 parents) (age range 5-19 years) recruited across different UK clinics (29). Participants were assigned to a case category when they scored ≥ 80 on performance IQ (PIQ) and had impaired expressive or receptive language skills (ELS or RLS score ≥ 1.5 SD below the normative mean of the general population) as assessed through the Clinical Evaluation of Language Fundamentals test (CELF-R, revised version) (63). If the same family comprised several cases, the most severe sibling was included into the analysis. Severity was determined by comparing the scores for ELS and RLS composite measure. If the cases performed at the same level on these measures, inclusion was determined by the most severe score on the NWR task. Overall, 137 unrelated cases met these criteria. A subset of this sample was used for the WES analysis conducted by Chen and colleagues where the rare variant was also reported (49). The ethical approval was given by local ethics committees of the hospitals involved in the consortium, and all subjects provided informed consent.

The York cohort consisted of 115 families (115 probands, 27 siblings, 184 parents), collected for a longitudinal study designed to investigate language and reading development (64). Participants entered the study for either having a family history of dyslexia or showing preschool language difficulties. A third group was recruited as control or typically developing children (TD). At 8 years the children were given a clinical diagnosis for SLI (N = 23, including 11 individuals with dyslexia) and dyslexia (N = 18). Ethical approval was obtained from the Yorkshire & The Humber - Humber Bridge NHS Research Ethics and the Ethics Committee at the University of York - Department of Psychology. Informed consent was given by all the participants and/or their parents.

Participants were assigned to a category of dyslexia from two other separate cohorts of families and individual cases. The criteria for inclusion were PIQ ≥ 85 and a score of 1 SD below the population mean on the single-word reading (READ) test measured through the British Ability Scales (BAS) (65) or the Wide Range Achievement Test (WRAT-R, revised version) (66). The first cohort comprised 290 families (689 siblings, 580 parents) (age range 6–27 years) described before (67). All the subjects were recruited by the Dyslexia Clinic at the Royal Berkshire Hospital, Reading, United Kingdom. When multiple children from the same family met the above criteria, we referred to other psychometric measures and selected the individual with the most severe phenotype on the following tasks: BAS single-word spelling accuracy (SPELL) (65), phonological decoding ability (CCN) (68), phonemic awareness (SPOON) (69) and orthographic coding skill assessed through the irregular word (CCI) (68) or forced word choice test (OLSON) (70). The Dyslexia Cases cohort, N = 592 singletons (age range 6–18 years), was collected from either the Dyslexia Research Centre clinics in Oxford and Reading, or the Aston Dyslexia and Development Clinic in Birmingham (71). In total, N = 471 participants met the above criteria for dyslexia. The study was approved by the Oxfordshire Psychiatric Research Ethics Committee and the Aston University Ethics Committee and informed consent was given by all the participants and/or their caregivers.

The ADHD cohort included 159 trios (age range 5-16 years) recruited from child psychiatry clinics in United Kingdom and Ireland (72). Participants were interviewed employing the Child and Adolescent Psychiatric Assessment (CAPA) (73). All the probands fulfilled the DSM-IV (American Psychiatric Society, 1994) diagnostic criteria for ADHD. Ethical approvals were acquired from the Trinity College Dublin in Ireland, and the University of Birmingham in UK. Informed consent was obtained from the parents of all the participants.

In addition to the clinical cohorts, we analysed a subset of samples from the Avon Longitudinal Study of Parents and Children (ALSPAC, see Supplementary material) (74,75) for which whole-genome (WGS) sequence data were generated as part of the UK10K project (76). Participants were excluded if they reported a non-white European ethnicity, presented with overt hearing difficulties or were affected by neurological conditions which may explain their language difficulties, such as autism, Asperger syndrome or pervasive developmental disorder. Individuals with missing information for the phenotypes of interest were also excluded.

The criteria for an assignment to a LI category were a PIQ ≥ 80 at age of 8.5 and a score either 1 SD below the population mean on the Wechsler Objective Language Dimensions (WOLD) comprehension test (77) at 8 years, or 1 SD below the population mean on both tests of syntax, and intelligibility and verbal fluency, as assessed through the Children’s Communication Checklist (CCC) (78) at 7.5 years. Assignment to the dyslexia category was given based on a diagnosis reported by the mother through a questionnaire and a PIQ score higher or equal to 85 at the age of 8.5. Assignment to the ADHD category was given if the subject met the DAWBA DSM-IV diagnostic criteria and had a PIQ > 80 at the age of 8.5. Please note that the ALSPAC study website contains details of all the data that is available through a fully searchable data dictionary and variable search tool ([https://www.bristol.ac.uk/alspac/researchers/our-data/](https://www.bristol.ac.uk/alspac/researchers/our-data/)).

Participants were assigned to a ‘super-control’ group (N=681) when they performed above the mean on tests of receptive (WOLD comprehension) and expressive (CCC syntax and intelligibility and fluency scores) language and had both verbal and non-verbal IQ ≥ 80. From the controls we also excluded any child who self-reported a need for special educational support or speech and language therapy and children who self-reported a diagnosis of neurodevelopmental disorder (ADHD, conduct disorder, anxiety disorder, depressive disorder, dyslexia, dyspraxia, dyscalculia, learning difficulties, speech and language difficulties or developmental delay).

Overall, in the ALSPAC sample, we identified: N = 170 LI cases, N = 31 dyslexia cases, N = 10 ADHD cases and N = 681 controls. Individuals with comorbid LI and dyslexia were included in the LI category.

Ethical approval for the study was obtained from the ALSPAC Law and Ethics Committee and the Local Research Ethics Committees. Consent for biological samples has been collected in accordance with the Human Tissue Act (2004). Informed consent was signed by the parents after having received a complete description of the study at the time of enrolment into the ALSPAC project, with the option for them or their children to withdraw at any time.

The gnomAD v2.1.1 database ([https://gnomad.broadinstitute.org/](https://gnomad.broadinstitute.org/)) was used as a reference European population (79). Of note, this source includes the UK10K samples.

### Quantitative analysis

While in most cohorts the number of carriers of the rare variant was too small to conduct any meaningful quantitative analysis, the ALSPAC dataset was suitable for this option. All the individuals with available genotype (N = 1867) were included into the analysis. Carriers (N = 42) and non-carriers (N = 1825) were compared for the mean score calculated for reading and language-related measures including: single-word reading accuracy (READ, target age: 7.5 year), single-word spelling accuracy (SPELL, 7.5 year), phoneme awareness (PHONEME, 7.5 year), working memory (MEMSPAN, 10.5 year), listening and comprehension test (WOLD, 8.5 year), phonological short-term memory test (NWR, 8.5 year), single-non-word reading accuracy (NW-READ, 9.5 year), average of the first seven scales from the Children’s Communication Checklist (CCC-average7, 7.5 year) and performance IQ (PIQ, 8.5 year) (See Supplementary Table S4 for more details). Statistical significance was calculated with the Student’s t-test (equal variances, two-tailed distribution).

### ATPase assay

The human *ATP2C2* coding sequence was cloned into the expression vector pcDNA™5/TO (Invitrogen, kind gift from Professor Galione (80). *ATP2C2* cDNA was obtained from SK-BR-3 cells maintained under standard conditions and amplified by high fidelity PCR. The primers were designed to include a *HindIII* site prior the Kozak consensus sequence and a *NotI* site immediately before the stop codon (See Supplementary Table S1 for the primer sequences). The *ATP2C2* open reading frame was placed under the control of a constitutive CMV promoter and in frame with a downstream HA epitope tag. The SPCA2-HA tagged plasmid was modified to introduce the missense variant (allele A). The mutated construct (SPCA2 p.V102M) was built using the GeneArt™ Site-Directed Mutagenesis kit (Thermo Fisher) according to the manufacturer’s instructions. The sequences were verified by Sanger sequencing (See Supplementary Table S1 for primer sequences).

Microsomes were prepared as described by Vandecaetsbeek (81) with minor modifications (See Supplementary material for details). HEK293 cells were transfected with the *ATP2C2* expressing constructs. The empty vector and untransfected cells were used as controls.

ATPase activity was evaluated using the PiColorLock™ Gold kit (Expedeon, Abcam) according to the manufacturer’s specifications. The assay was conducted in a 96-well plate and each sample run in duplicates. The amount of free inorganic phosphate (Pi) released from each reaction was calculated by plotting the absorbance against a series of Pi standards of known concentration.

The data generated were presented as mean percentages ± SD of two independent experiments. The activity of each sample was expressed relative to that detected in the empty vector, normalized to 100%. Pair-wise comparison between samples was performed by the Welch’s t-test (unequal variances, two-tailed distribution).

For Western Blot analysis, 5 µg of microsomes were mixed with Laemmli buffer and incubated at room temperature (RT) for 30 minutes. Samples (5 µg per well) were separated on a precast 4-12% NuPAGE™ polyacrylamide gradient gel (Invitrogen). The anti-SPCA2 (Abcam) and anti-GAPDH (Abcam) antibodies were diluted 1:4000 and 1:5000 in 0.1% TBS-Tween™20 (Fisher BioReagents™), respectively.

## RESULTS

### Identification of the missense variant Chr16:84405221

WES analysis in six families selected for SLI led to the identification of a rare variant (Chr16:84405221) in the *ATP2C2* gene (Figure 1). This variant was found in two children and their father, all diagnosed with SLI. The same variant was previously reported in an exome sequencing study performed on 43 unrelated probands affected by severe language impairment (49). *ATP2C2* was previously associated with SLI (29). Sanger sequencing in the other family members showed imperfect segregation as the variant was not transmitted to the other two affected sons (Figure 1). The variant was inherited from the father who also scored below expected (i.e., < 15th %ile) on a standardized expressive language measure. The variant leads to a missense substitution from a valine to a methionine in exon three spanning the ATPase domain. Bioinformatics tools for pathogenicity suggest the change is likely to be benign (Table 1). The two algorithms that predicted a functional effect, albeit subtle, on protein activity and stability were the Mutation Assessor and the INPS-MD tools (Table 1).

### Genotyping of Chr16:84405221 in the SLI, dyslexia and ADHD cohorts

The variant was genotyped in independent cohorts with language impairment, dyslexia, or ADHD (Table 2). The frequency of the variant was 2.2% in the SLIC cohort selected for SLI and 1.8% across all the unrelated LI cases (combined LI samples, N = 360), including the discovery family and previously published data (49). The frequency was 0.8% in the combined dyslexia samples (N = 520), and 0.7% in the combined ADHD cases (N = 151).

In the ALSPAC super-controls (N = 681) the variant occurred at a frequency of 1%. The frequency in the global gnomAD database is of 0.5% (N = 59718) and, of 0.9% in the European population (N = 24046).

Although the frequency was higher in the LI cases (1.5%) than the controls (1%) derived from the ALSPAC cohort, the difference was not statistically significant when examined with a two-tailed Fisher’s exact test (p = 0.152).

### Mean score of language and reading-related phenotypes

We evaluated the effects of the variant on a range of language-related quantitative measures available for the ALSPAC cohort (Supplementary Table S4). The carriers (N = 42) had lower mean scores compared to the non-carriers (N = 1825) on all tested measures, reaching nominal statistical significance (t-test, two-tailed distribution) for six traits: READ (p = 0.016), SPELL (p = 0.006), PHONEME (p = 0.026), NW_READ (p = 0.002), VIQ (p = 0.006) and PIQ (p = 0.022) (Table 3). NW_READ survived a conservative Bonferroni correction for multiple comparison (10 tests, threshold for significance = 0.005; see also Supplementary Table S5 for analysis through the Wald test and Supplementary Figure S2 for box plots). It is worth noting that the 10 tests are not completely independent given the correlation across these measures (83>.

View this table:
[Table 3.](http://medrxiv.org/content/early/2021/01/15/2021.01.12.21249315/T3)

Table 3. Mean score of different phenotypes across carriers and non-carriers (ALSPAC cohort)

### ATPase assay

An effect of allele A on the ATPase activity of SPCA2 was tested with a colorimetric assay to quantify the inorganic phosphate released during the catalytic reaction (Figure 2A). Microsomes were isolated from HEK293 cells overexpressing the wildtype or mutated SPCA2-HA tagged protein. The activity of the overexpressed SPCA2 WT sample was higher compared to the negative controls (empty vector and untransfected cells) across all the conditions. SPCA2 p.V102M displayed on average a lower ATPase activity compared to that of the WT protein (Figure 2A). WB confirmed the overexpression of SPCA2 for the WT and mutated constructs while no protein was detected for the negative controls, suggesting a very minimal effect of baseline endogenous protein (Figure 2B).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.12.21249315/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/01/15/2021.01.12.21249315/F2)

Figure 2. Effect of the mutation on ATPase activity.
**A**) Concentration of free inorganic phosphate (Pi) released during the ATPase catalytic reaction. The activity is expressed relative to the empty vector, normalized to 100%. The results are shown as mean percentages ± SD of two independent experiments. The samples are: untransfected HEK293 cells (UT), microsomes isolated from HEK293 cells transiently transfected to overexpress the recombinant SPCA2-HA tagged (SPCA2 WT) or mutated protein (SPCA2 p.V102M) and microsomes isolated from HEK293 cells transiently transfected with the empty vector (E). **B**) Western blot. Each lane was loaded with 5 µg of microsomes. The membranes were probed with antibodies against SPCA2 (top, green) or GAPDH (36 kDa) (bottom, red). SPCA2 can be detected when overexpressed, while no signal is observed in the untransfected cells and the control samples. The band corresponds to the predicted SPCA2-HA tagged molecular weight (105.2 kDa).

## DISCUSSION

In this study we investigated the effects of a rare *ATP2C2* variant on language impairment and related disorders. Overall, we found that this variant had a higher frequency in LI samples (1.8%, N = 360) compared to controls (1%, N = 681, ALSPAC) and subjects with dyslexia (0.8%, N = 520) or ADHD (0.7%, N = 150). The frequency reached 2.2% in a clinical cohort ascertained for SLI (N = 161; Table 2). In a general population sample derived from the ALSPAC cohort, carriers of the rare variants presented lower scores compared to non-carries on a range of language-related measures, with the stronger effect observed for the nonword reading measure (p = 0.002; Table 3). Taken together, our results suggest that the *ATP2C2* missense variant is associated to LI and lower scores on reading and language measures. The presence of the *ATP2C2* variant in controls shows that it is not sufficient on its own to cause a disorder, and its association with a range of measures suggests a pleiotropic effect on cognitive traits. These observations are in line with the emerging findings in the field of complex traits genetics and, specifically, for neurodevelopmental and psychiatric disorders. For example, under the double-hit model, a risk factor might lead to a specific disorder or a range of conditions depending on the interactions with other genetic risk variants (51,52,84>. Under this model, the rare *ATP2C2* variant might increase the risk for LI or affect performance on reading- and language-related measures when combined with different susceptibility factors, as observed for other rare variants implicated in LI. These include a microdeletion in *ZNF277* (42), a rare stop-gain variant in *ERC1* (49) and nonsynonymous mutations in *NFXL1* (85), *GRIN2B* and *SRPX2* (Chen et al., 2017). Increasing evidence shows that the penetrance of rare risk factors is modulated by the cumulative effect of common risk alleles, even in severe neurodevelopmental disorders usually assumed to have a monogenic origin (86). Different manifestations determined by the interplay between common and rare variants have been reported for conditions such as schizophrenia (87,88), ADHD (89), autism spectrum disorder (90), intellectual disability (91) and outcomes such as educational attainment (92). Rare risk factors causing psychiatric disorders have also been reported to have an effect on cognitive abilities in the general population (93,94). This is in line with what we observe for the *ATP2C2* variant which showed associations with poor performance on cognitive traits in the ALSPAC sample regardless of any disorder group assignement (Table 3). The variant is associated to lower mean scores on a range of cognitive measures, with the strongest effect observed for reading abilities rather than language skills. However, although tempting, we should not speculate on specific genotype-phenotype relationship. Instead our results support the emerging idea that genetic risk factors have generalised effects across cognitive and neurodevelopmental traits. The same concept is illustrated also for common variants through polygenic risk scores (PRS) analysis. For example, PRS for schizophrenia showed associations for a number of traits including poor cognitive outcomes (95) but also enhanced creativity (96), language performance (97) and hearing (98). Similarly, polygenic risk scores (PRS) for general cognitive abilities show effects on reading and language measures (15,99–101>.

The *ATP2C2* gene encodes for an ATPase that transports Ca2+ and Mn2+ into the Golgi lumen (102>. ATPases are emerging as a new class of proteins involved in neurodevelopmental disorders (103– 106). Ca2+ is involved in the regulation of several neuronal processes such as neuronal migration, axon guidance, synaptic plasticity, memory, neuronal excitability, neurosecretion and aging (107– 110). Excessive accumulation of Mn2+ inside the cells is toxic, induces apoptosis and causes manganism, a neurological disorder characterised by tremors and spasm (111,112). Additionally, dysregulation of these cations has been linked to Alzheimer’s, Parkinson’s and Huntington’s disease (113,114). The transport of Ca2+ and Mn2+ ions inside the Golgi is essential for the activity of enzymes involved in mechanisms such as neurotransmitter synthesis, membrane trafficking, metabolism, protein modification, sorting and folding (110). Previous studies have shown that *ATP2C2* is expressed in postnatal rat hippocampal neurons (115) and in adult human brain (102,115,116), suggesting a possible role of *ATP2C2* in maintaining ionic homeostasis in neurons.

The missense variant described in the current study is localised in the cation transporting N-terminal domain and causes the amino acid substitution from valine to methionine (Figure 1). The latter residue is characterised by an unbranched side chain that provides increased flexibility. In addition to this, methionine has a sulphur group that can influence the binding preference for metal ions, and be reversibly oxidised converting the hydrophobic amino acid into a highly polar one (117,118>. Bioinformatic predictions do not attribute substantial detrimental effects to the substitution, and therefore any potential effects are expected to be subtle (Table 1). We tested experimentally whether the rare allele A had an effect on the SPCA2 activity, through an ATPase assay. This approach allowed us to measure the hydrolysis of ATP in a quick and sensitive way. Independently of the concentrations tested, the rare variant displayed on average a lower ATPase activity compared to SPCA2 WT (Figure 2). Although not statistically significant, the data suggested a possible impact of the variant on the velocity of the reaction and therefore the cation binding affinity. More specifically, the mutated transporter tended to present a lower ATP hydrolysis indicative of a lower turnover rate and consequently a higher ion affinity. In line with the bioinformatics prediction, the rare variant does not affect the functionality of the protein but might alter its performance which in turn, in combination with other factors, might contribute to neurodevelopmental phenotypes.

In summary, we identified a rare missense variant (chr16:84405221, GRCh38.p12) in *ATP2C2*, a previously reported candidate gene for LI. We showed that the rare allele A increases the risk for LI cohorts and is associated with more general cognitive abilities in a population-based cohort. Preliminary functional characterization of the rare allele suggested that the rare variant influences the ATPase activity of SPCA2 supporting follow up functional studies to fully understand the role of this variant and the ATP2C2 gene. WES studies for language disorders have been limited to samples, which makes the result interpretation particularly challenging. Our results demonstrate that combined analyses across disorders and phenotypes provide a valid route to validate WES findings. Taken together, our results further support a role of *ATP2C2* in language impairment and pinpoint to a specific rare variant which might impact ATPase activity of the SPCA2 protein.

## Supporting information

Supplementary Information [[supplements/249315_file02.pdf]](pending:yes)

## Data Availability

Data used in this study are available from UK10K Data Access Committee and the ALSPAC Executive Committee for researchers who meet the criteria for access to confidential data.

## Disclosure

The authors have declared that no competing interests existed at the time of publication.

## Acknowledgement

SP if funded by the Royal Society. This work was supported by an Action Medical Research Action/The Chief Scientist (CSO) Office, Scotland grant (GN2614) and a Cunningham Trust grant to SP and SJP. Support to the analysis was provided by the St Andrews Bioinformatics Unit funded by the Wellcome Trust [grant 105621/Z/14/Z]. Assessment of the Aston cohort was supported by funding from The Waterloo Foundation to JBT and SP [797–1720]. Analysis of the discovery pedigree was supported by the University of Kansas grant (NIH:NIDCD 5 R01 DC001803). Analysis of the York cohort was funded by Wellcome Trust Programme Grant 082036/B/07/Z. DFN is currently supported by Oxford Brookes University funds, the Leverhulme Trust and the British Academy. This work was, in part, completed while DFN was at the Wellcome Trust Centre for Human Genetics, Oxford as an MRC Career Development Fellow (G1000569/1). We would like to thank the members of all teams who collected the data and the families who participated.

We are extremely grateful to all the families who took part in the ALSPAC study, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists and nurses. The UK Medical Research Council and Wellcome (Grant ref: 217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. This publication is the work of the authors and DFN will serve as guarantors for the analysis of the ALSPAC data presented in this paper. A comprehensive list of grants funding is available on the ALSPAC website ([http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf](http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf)). Whole-genome sequencing of the ALSPAC samples was performed as part of the UK10K onsortium (a full list of investigators who contributed to the generation of the data is available from [www.UK10K.org.uk](http://www.UK10K.org.uk)). ALSPAC GWAS data was generated by Sample Logistics and Genotyping Facilities at Wellcome Sanger Institute and LabCorp (Laboratory Corporation of America) using support from 23andMe. We thank Michael Gill for sharing data for the ADHD cohort and Arathi Jeyaratnam for the technical support.

*   Received January 12, 2021.
*   Revision received January 12, 2021.
*   Accepted January 15, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## REFERENCES

1.  1.Bishop, D. V. M. (2014) Ten questions about terminology for children with unexplained language problems. Ten questions about terminology for children with unexplained language problems. Int. J. Lang. Commun. Disord. (2014), 49, 381–415.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/1460-6984.12101&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25142090&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

2.  2.Bishop, D. V. M., Snowling, M. J., Thompson, P. A., et al. (2016) CATALISE: A multinational and multidisciplinary Delphi consensus study. Identifying language impairments in children. PLoS One, 11, e0158753.
    
    

3.  3.Bishop, D. V. M. (2017) Why is it so hard to reach agreement on terminology? The case of developmental language disorder (DLD). Why is it so hard to reach agreement on terminology? The case of developmental language disorder (DLD). Int. J. Lang. Commun. Disord. (2017), 52, 671–680.
    
    

4.  4.Norbury, C. F., Gooch, D., Wray, C., et al. (2016) The impact of nonverbal ability on prevalence and clinical presentation of language disorder: evidence from a population study. J. Child Psychol. Psychiatry Allied Discip., 57, 1247–1257.
    
    

5.  5.Tomblin, J. B., Records, N. L., Buckwalter, P., et al. (1997) Prevalence of specific language impairment in kindergarten children. J. Speech. Lang. Hear. Res., 40, 1245–60.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1044/jslhr.4006.1245&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9430746&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000071027200003&link_type=ISI) 

6.  6.Conti-Ramsden, G., Durkin, K., Simkin, Z., et al. (2009) Specific language impairment and school outcomes. I: Identifying and explaining variability at the end of compulsory education. Int. J. Lang. Commun. Disord., 44, 15–35.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/13682820801921601&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18608604&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

7.  7.Conti-Ramsden, G., Mok, P. L. H., Pickles, A., et al. (2013) Adolescents with a history of specific language impairment (SLI): Strengths and difficulties in social, emotional and behavioral functioning. Res. Dev. Disabil., 34, 4161–4169.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ridd.2013.08.043&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24077068&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

8.  8. J. Snowling, John W. Adams, D. V. M, M., Adams, J. W., Bishop, D. V. M., et al. (2001) Educational attainments of school leavers with a preschool history of speech-language impairments. Int J Lang Commun Disord, 36, 173–183.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/13682820010019892&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11344593&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000168381500002&link_type=ISI) 

9.  9.Conti-Ramsden, G., Durkin, K., Toseeb, U., et al. (2018) Education and employment outcomes of young adults with a history of developmental language disorder. Int. J. Lang. Commun. Disord., 53, 237–255.
    
    

10. 10.Botting, N. (2020) Language, literacy and cognitive skills of young adults with developmental language disorder (DLD). Int. J. Lang. Commun. Disord., 55, 255–265.
    
    

11. 11.McArthur, G. M., Hogben, J. H., Edwards, V. T., et al. (2000) On the “specifics” of specific reading disability and specific language impairment. J Child Psychol Psychiatry, 41, 869–874.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0021963099006186&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11079429&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000165233400006&link_type=ISI) 

12. 12.Snowling, M. and Stothard, S. E. (2000) Is Preschool Language Impairment a Risk Factor for Dyslexia in Adolescence? 41, 587–600.
    
    

13. 13.Mueller, K. L. and Tomblin, J. B. (2012) Exmaining the comorbidity of language disorders and ADHD. Biol. Bull., 221, 18–34.
    
    

14. 14.Newbury, D. F., Paracchini, S., Scerri, T. S., et al. (2011) Investigation of dyslexia and SLI risk variants in reading- and language-impaired subjects. Behav. Genet., 41, 90–104.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10519-010-9424-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21165691&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

15. 15.Gialluisi, A., Andlauer, T. F. M., Mirza-Schreiber, N., et al. (2020) Genome-wide association study reveals new insights into the heritability and genetic correlates of developmental dyslexia. Mol. Psychiatry.
    
    

16. 16.Barry, J. G., Yasin, I. and Bishop, D. V (2007) Heritable risk factors associated with language impairments. Genes Brain Behav, 6, 66–76.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1601-183X.2006.00232.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17233642&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000243473900006&link_type=ISI) 

17. 17.Bishop, D. V. M., North, T. and Donlan, C. (1995) GENETIC BASIS OF SPECIFIC LANGUAGE IMPAIRMENT: EVIDENCE FROM A TWIN STUDY. Dev. Med. Child Neurol., 37, 56–71.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7828787&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995PX60500006&link_type=ISI) 

18. 18.Lewis, B. A. and Thompson, L. A. (1992) A study of developmental speech and language disorders in twins. J Speech Hear Res, 35, 1086–1094.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1447920&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

19. 19.Choudhury, N., & Benasich, A. A. (2003) A Family Aggregation Study: The Influence of Family History and Other Risk Factors on Language Development. J Speech Lang Hear Res., 46, 261–272.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1044/1092-4388(2003/021)&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14700370&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184906200002&link_type=ISI) 

20. 20.Tomblin, J. B. and Buckwalter, P. R. (1998) Heritability of poor language achievement among twins. J. Speech, Lang. Hear. Res., 41, 188–199.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9493744&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000071915700016&link_type=ISI) 

21. 21.Wiszniewski, W., Hunter, J. V, Hanchard, N. A., et al. (2013) TM4SF20 ancestral deletion and susceptibility to a pediatric disorder of early language delay and cerebral white matter hyperintensities. Am. J. Hum. Genet., 93, 197–210.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2013.05.027&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23810381&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

22. 22.Vernes, S. C., Newbury, D. F., Abrahams, B. S., et al. (2008) A functional genetic link between distinct developmental language disorders. N. Engl. J. Med., 359, 2337–2345.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa0802828&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18987363&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000261270300005&link_type=ISI) 

23. 23.Bartlett, C. W., Flax, J. F., Logue, M. W., et al. (2002) A major susceptibility locus for specific language impairment is located on 13q21. Am. J. Hum. Genet., 71, 45–55.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/341095&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12048648&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176307700005&link_type=ISI) 

24. 24.Bartlett, C. W., Flax, J. F., Logue, M. W., et al. (2004) Examination of potential overlap in autism and language loci on chromosomes 2, 7, and 13 in two independent samples ascertained for specific language impairment. Hum. Hered., 57, 10–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000077385&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15133308&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221270600002&link_type=ISI) 

25. 25.Newbury, D. F., Cleak, J. D., Ishikawa-Brush, Y., et al. (2002) A genomewide scan identifies two novel loci involved in specific language impairment. Am. J. Hum. Genet., 70, 384–398.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/338649&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11791209&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173254300011&link_type=ISI) 

26. 26.Monaco, A. P. (2004) Highly significant linkage to the SLI1 locus in an expanded sample of individuals affected by specific language impairment. Am. J. Hum. Genet., 74, 1225–1238.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/421529&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15133743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221651900015&link_type=ISI) 

27. 27.Monaco, A. P. (2007) Multivariate linkage analysis of specific language impairment (SLI). Ann. Hum. Genet., 71, 660–673.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1469-1809.2007.00361.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17388790&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248121700011&link_type=ISI) 

28. 28.Falcaro, M., Pickles, A., Newbury, D. F., et al. (2008) Genetic and phenotypic effects of phonological short-term memory and grammatical morphology in specific language impairment. Genes, Brain Behav., 7, 393–402.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1601-183X.2007.00364.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18005161&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255286200002&link_type=ISI) 

29. 29.Newbury, D. F., Winchester, L., Addis, L., et al. (2009) CMIP and ATP2C2 Modulate Phonological Short-Term Memory in Language Impairment. Am. J. Hum. Genet., 85, 264–272.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2009.07.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19646677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269332900015&link_type=ISI) 

30. 30.Bishop, D. V. M., North, T. and Donlan, C. (1996) Nonword repetition as a behavioural marker for inherited language impairment: Evidence from a twin study. J. Child Psychol. Psychiatry Allied Discip., 37, 391–403.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1469-7610.1996.tb01420.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8735439&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996UL11000005&link_type=ISI) 

31. 31.Graham, S. A. and Fisher, S. E. (2015) Understanding Language from a Genomic Perspective. Annu. Rev. Genet., 49, 131–160.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-genet-120213-092236&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26442845&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

32. 32.Nudel, R., Simpson, N. H., Baird, G., et al. (2014) Genome-wide association analyses of child genotype effects and parent-of-origin effects in specific language impairment. Genes, Brain Behav., 13, 418–429.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/gbb.12127&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24571439&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

33. 33.Pettigrew, K. A., Frinton, E., Nudel, R., et al. (2016) Further evidence for a parent-of-origin effect at the NOP9 locus on language-related phenotypes. J. Neurodev. Disord., 8.
    
    

34. 34.St Pourcain, B., Cents, R. A. M., Whitehouse, A. J. O., et al. (2014) ARTICLE Common variation near ROBO2 is associated with expressive vocabulary in infancy. Nat. Commun., 5.
    
    

35. 35.Luciano, M., Evans, D. M., Hansell, N. K., et al. (2013) A genome-wide association study for reading and language abilities in two population cohorts. Genes, Brain Behav., 12, 645–652.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/gbb.12053&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23738518&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

36. 36.Eicher, J. D., Powers, N. R., Miller, L. L., et al. (2013) Genome-wide association study of shared components of reading disability and language impairment. Genes, Brain Behav., 12, 792–801.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/gbb.12085&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24024963&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

37. 37.Gialluisi, A., Newbury, D. F., Wilcutt, E. G., et al. (2014) Genome-wide screening for DNA variants associated with reading and language traits. Genes, Brain Behav., 13, 686–701.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/gbb.12158&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25065397&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

38. 38.Simpson, N. H., Addis, L., Brandler, W. M., et al. (2014) Increased prevalence of sex chromosome aneuploidies in specific language impairment and dyslexia. Dev. Med. Child Neurol., 56, 346–353.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/dmcn.12294&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24117048&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

39. 39.Simpson, N. H., Ceroni, F., Reader, R. H., et al. (2015) Genome-wide analysis identifies a role for common copy number variants in specific language impairment. Eur. J. Hum. Genet., 23, 1370–1377.
    
    

40. 40.Kalnak, N., Stamouli, S., Peyrard-Janvid, M., et al. (2018) Enrichment of rare copy number variation in children with developmental language disorder. Clin. Genet., 94, 313–320.
    
    

41. 41.Gulhan Ercan-Sencicek, A., Davis Wright, N. R., Sanders, S. J., et al. (2011) A balanced t(10;15) translocation in a male patient with developmental language disorder.
    
    

42. 42.Ceroni, F., Simpson, N. H., Francks, C., et al. (2014) Homozygous microdeletion of exon 5 in ZNF277 in a girl with specific language impairment. Eur. J. Hum. Genet., 22, 1165–1171.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ejhg.2014.4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24518835&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

43. 43.Smith, A. W., Holden, K. R., Dwivedi, A., et al. (2014) Deletion of 16q24.1 Supports a Role for the ATP2C2 Gene in Specific Language Impairment. J. Child Neurol., 30, 517–521.
    
    

44. 44.Pettigrew, K. A., Reeves, E., Leavett, R., et al. (2015) Copy number variation screen identifies a rare de Novo deletion at chromosome 15q13.1-13.3 in a child with language impairment. PLoS One, 10.
    
    

45. 45.Hildebrand, M. S., Jackson, V. E., Scerri, T. S., et al. (2020) Severe childhood speech disorder: Gene discovery highlights transcriptional dysregulation. Neurology, 94, e2148–e2167.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0000000000009441&link_type=DOI) 

46. 46.Lai, C. S. L., Fisher, S. E., Hurst, J. A., et al. (2001) A forkhead-domain gene is mutated in a severe speech and language disorder. Nature, 413, 519–523.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/35097076&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11586359&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171340500046&link_type=ISI) 

47. 47.Devanna, P., Chen, X. S., Ho, J., et al. (2018) Next-gen sequencing identifies non-coding variation disrupting miRNA-binding sites in neurological disorders. Mol. Psychiatry, 23, 1375–1384.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2017.30&link_type=DOI) 

48. 48.Eising, E., Carrion-Castillo, A., Vino, A., et al. (2019) A set of regulatory genes co-expressed in embryonic human brain is implicated in disrupted speech development. Mol. Psychiatry, 24, 1065–1078.
    
    

49. 49.Chen, X. S., Reader, R. H., Hoischen, A., et al. (2017) Next-generation DNA sequencing identifies novel gene variants and pathways involved in specific language impairment. Sci. Rep., 7, 46105.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/srep46105&link_type=DOI) 

50. 50.Girirajan, S., Brkanac, Z., Coe, B. P., et al. (2011) Relative Burden of Large CNVs on a Range of Neurodevelopmental Phenotypes. PLoS Genet., 7, e1002334.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1002334&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22102821&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

51. 51.Girirajan, S., Rosenfeld, J. A., Cooper, G. M., et al. (2010) A recurrent 16p12.1 microdeletion supports a two-hit model for severe developmental delay. Nat. Genet., 42, 203–209.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.534&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20154674&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274912400007&link_type=ISI) 

52. 52.O ‘ Roak, B. J., Vives, L., Girirajan, S., et al. (2012) Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations. Nature, 485, 246–250.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature10989&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22495309&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000303799800043&link_type=ISI) 

53. 53.Rice, M. L., Smith, S. D. and Gayán, J. (2009) Convergent genetic linkage and associations to language, speech and reading measures in families of probands with Specific Language Impairment. J. Neurodev. Disord., 1, 264–282.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11689-009-9031-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19997522&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

54. 54.Andres, E. M., Earnest, K. K., Smith, S. D., et al. (2020) Pedigree-Based Gene Mapping Supports Previous Loci and Reveals Novel Suggestive Loci in Specific Language Impairment. J. Speech, Lang. Hear. Res., 1–16.
    
    

55. 55.Rice, Mabel L.; Wexler, K. (2001) Rice/Wexler Test of Early Grammatical Impairment. Rice/Wexler Test of Early Grammatical Impairment.; San Antoni.; (2001).
    
    

56. 56.Sim, N. L., Kumar, P., Hu, J., et al. (2012) SIFT web server: Predicting effects of amino acid substitutions on proteins. Nucleic Acids Res., 40.
    
    

57. 57.Reva, B., Antipin, Y. and Sander, C. (2011) Predicting the functional impact of protein mutations: Application to cancer genomics. Nucleic Acids Res., 39.
    
    

58. 58.Choi, Y. and Chan, A. P. (2015) PROVEAN web server: A tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics, 31, 2745–2747.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv195&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25851949&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

59. 59.Adzhubei, I. A., Schmidt, S., Peshkin, L., et al. (2010) A method and server for predicting damaging missense mutations. Nat Methods, 7, 248–249.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth0410-248&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20354512&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276150600004&link_type=ISI) 

60. 60.Pejaver, V., Urresti, J., Lugo-Martinez, J., et al. (2017) MutPred2: inferring the molecular and phenotypic impact of amino acid variants..
    
    

61. 61.Savojardo, C., Fariselli, P., Martelli, P. L., et al. (2016) INPS-MD: A web server to predict stability of protein variants from sequence and structure. Bioinformatics, 32, 2542–2544.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw192&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27153629&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

62. 62.Schwarz, J. M., Cooper, D. N., Schuelke, M., et al. (2014) Mutationtaster2: Mutation prediction for the deep-sequencing age. Mutationtaster2: Mutation prediction for the deep-sequencing age. Nat. Methods (2014), 11, 361–362.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.2890&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24681721&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000333749900007&link_type=ISI) 

63. 63.Semel, E. M., Wiig, E. H., Secord, W. (1992) Clinical Evaluation of Language Fundamental, Revised. Psychol. Corp. San Antonio.
    
    

64. 64.Nash, H. M., Hulme, C., Gooch, D., et al. (2013) Preschool language profiles of children at family risk of dyslexia: Continuities with specific language impairment. J. Child Psychol. Psychiatry Allied Discip.
    
    

65. 65.Elliot, C. D., Murray, D. J., and Pearson, L. S. (1983) British Abilities Scales. NFER-Nelson, Wind. UK.
    
    

66. 66.Jastak, J., Wilkinson, G. S. (1984) Wide Range Achievement Test – Revised (WRAT-R). Psychol. Corp. San Antonio, TX.
    
    

67. 67.Scerri, T. S., Paracchini, S., Morris, A., et al. (2010) Identification of candidate genes for dyslexia susceptibility on chromosome 18. PLoS One, 5, e13712.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0013712&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21060895&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

68. 68.Castles, A. and Coltheart, M. (1993) Varieties of developmental dyslexia. Cognition, 47, 149–180.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0010-0277(93)90003-E&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8324999&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1993LD02000003&link_type=ISI) 

69. 69.Frederickson, N., Frith, U., Reason, R. (1997) Phonological assessment battery (PhAB). Wind. NFER- Nelson.
    
    

70. 70.Olson, R., Forsberg, H., Wise, B., & Rack, J. (1994) Measurement of word recognition, orthographic, and phonological skills. In G. R. Lyon (Ed.), Frames of reference for the assessment of learning disabilities: New views on measurement issues. Paul H Brookes Publ. Co.
    
    

71. 71.Scerri, T. S., Macpherson, E., Martinelli, A., et al. (2017) The DCDC2 deletion is not a risk factor for dyslexia. Transl. Psychiatry, 7, e1182.
    
    

72. 72.Kent, L., Doerry, U., Hardy, E., et al. (2002) Evidence that variation at the serotonin transporter gene influences susceptibility to attention deficit hyperactivity disorder (ADHD): Analysis and pooled analysis. Mol. Psychiatry, 7, 908–912.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.mp.4001100&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12232786&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000178018000017&link_type=ISI) 

73. 73.Angold, A., Prendergast, M., Cox, A., et al. (1995) The Child and Adolescent Psychiatric Assessment (CAPA). Psychol. Med., 25, 739–753.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S003329170003498X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7480451&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995RT07200007&link_type=ISI) 

74. 74.Boyd, A., Golding, J., Macleod, J., et al. (2013) Cohort profile: The ‘Children of the 90s’-The index offspring of the avon longitudinal study of parents and children. Int. J. Epidemiol., 42, 111–127.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dys064&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22507743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316699300012&link_type=ISI) 

75. 75.Fraser, A., Macdonald-wallis, C., Tilling, K., et al. (2013) Cohort profile: The avon longitudinal study of parents and children: ALSPAC mothers cohort. Int. J. Epidemiol., 42, 97–110.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dys066&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22507742&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316699300011&link_type=ISI) 

76. 76.Walter, K., Min, J. L., Huang, J., et al. (2015) The UK10K project identifies rare variants in health and disease. Nature, 526, 82–89.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature14962&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26367797&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

77. 77.Rust, J. (1996) WOLD Wechsler Objective Language Dimensions Manual. WOLD Wechsler Objective Language Dimensions Manual ; The Psychological Corporation, London, UK, (1996).
    
    

78. 78.Bishop, D. V. M. (1998) Development of the Children’s Communication Checklist (CCC): a method for assessing qualitative aspects of communicative impairment in children. J Child Psychol Psychiatry, 39, 879–891.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0021963098002832&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9758196&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000076168300009&link_type=ISI) 

79. 79.Karczewski, K. J., Francioli, L. C., Tiao, G., et al. (2020) The mutational constraint spectrum quantified from variation in 141,456 humans. Nature, 581, 434–443.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2308-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32461654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

80. 80.Rietdorf, K., Funnell, T. M., Ruas, M., et al. (2011) Two-pore channels form homo- and heterodimers. J. Biol. Chem., 286, 37058–37062.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODYvNDMvMzcwNTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMS8xNS8yMDIxLjAxLjEyLjIxMjQ5MzE1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

81. 81.Vandecaetsbeek, I., Holemans, T., Wuytack, F., et al. (2014) High-throughput measurement of the Ca2+-dependent ATPase activity in COS microsomes. Cold Spring Harb. Protoc., 2014, 865–875.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25086012&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

82. 82.Berman, H. M., Battistuz, T., Bhat, T. N., et al. (2002) The protein data bank. Acta Crystallogr. Sect. D Biol. Crystallogr., 58, 899–907.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1107/S0907444902003451&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12037327&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176271000003&link_type=ISI) 

83. 83.Scerri, T. S., Morris, A. P., Buckingham, L. L., et al. (2011) DCDC2, KIAA0319 and CMIP are associated with reading-related traits. Biol. Psychiatry, 70, 237–245.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2011.02.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21457949&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000292633600007&link_type=ISI) 

84. 84.Parenti, I., Rabaneda, L. G., Schoen, H., et al. (2020) Neurodevelopmental Disorders: From Genetics to Functional Pathways. Neurodevelopmental Disorders: From Genetics to Functional Pathways. Trends Neurosci. (2020).
    
    

85. 85.Villanueva, P., Nudel, R., Hoischen, A., et al. (2015) Exome Sequencing in an Admixed Isolated Population Indicates NFXL1 Variants Confer a Risk for Specific Language Impairment. PLoS Genet., 11, e1004925.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004925&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25781923&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

86. 86.Niemi, M. E. K., Martin, H. C., Rice, D. L., et al. (2018) Common genetic variants contribute to risk of rare severe neurodevelopmental disorders. Nature, 562, 268–271.
    
    

87. 87.Tansey, K. E., Rees, E., Linden, D. E., et al. (2016) Common alleles contribute to schizophrenia in CNV carriers. Mol. Psychiatry, 21, 1085–1089.
    
    

88. 88.Bergen, S. E., Ploner, A., Howrigan, D., et al. (2019) Joint contributions of rare copy number variants and common SNPs to risk for schizophrenia. Am. J. Psychiatry, 176, 29–35.
    
    

89. 89.Martin, J., O’Donovan, M. C., Thapar, A., et al. (2015) The relative contribution of common and rare genetic variants to ADHD. Transl. Psychiatry, 5.
    
    

90. 90.Weiner, D. J., Wigdor, E. M., Ripke, S., et al. (2017) Polygenic transmission disequilibrium confirms that common and rare variation act additively to create risk for autism spectrum disorders. Nat. Genet., 49, 978–985.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3863&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28504703&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

91. 91.Kurki, M. I., Saarentaus, E., Pietiläinen, O., et al. (2019) Contribution of rare and common variants to intellectual disability in a sub-isolate of Northern Finland. Nat. Commun., 10, 1–15.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-09078-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30602773&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

92. 92.Ganna, A., Genovese, G., Howrigan, D. P., et al. (2016) Ultra-rare disruptive and damaging mutations influence educational attainment in the general population HHS Public Access. Nat Neurosci, 19, 1563–1565.
    
    

93. 93.Kendall, K. M., Bracher-Smith, M., Fitzpatrick, H., et al. (2019) Cognitive performance and functional outcomes of carriers of pathogenic copy number variants: Analysis of the UK Biobank. Br. J. Psychiatry, 214, 297–304.
    
    

94. 94.Stefansson, H., Meyer-Lindenberg, A., Steinberg, S., et al. (2014) CNVs conferring risk of autism or schizophrenia affect cognition in controls. Nature, 505, 361–366.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12818&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24352232&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000329621800034&link_type=ISI) 

95. 95.Hubbard, L., Tansey, K. E., Rai, D., et al. (2016) Evidence of common genetic overlap between schizophrenia and cognition. Schizophr. Bull., 42, 832–842.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/schbul/sbv168&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26678674&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

96. 96.Power, R. A., Steinberg, S., Bjornsdottir, G., et al. (2015) Polygenic risk scores for schizophrenia and bipolar disorder predict creativity. Nat. Neurosci., 18, 953–955.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.4040&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26053403&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

97. 97.Rajagopal, V. M., Ganna, A., Coleman, J. R. I., et al. (2020) Genome-wide association study of school grades identifies a genetic overlap between language ability, psychopathology and creativity. BioRxiv, 16, 1–29.
    
    

98. 98.Schmitz, S. J. J., Abbondanza, A., et al. (2020) Polygenic risk scores implicate genetic pathways involved in neurodevelopmental disorders in hearing thresholds and hearing asymmetry in children. bioRxiv Genet., 2020.07.22.215376.
    
    

99. 99.Gialluisi, A., Andlauer, T. F. M., Mirza-Schreiber, N., et al. (2019) Genome-wide association scan identifies new variants associated with a cognitive predictor of dyslexia. Transl. Psychiatry, 9, 15.
    
    

100.100.Newbury, D. F., Gibson, J. L., Conti-Ramsden, G., et al. (2019) Using polygenic profiles to predict variation in language and psychosocial outcomes in early and middle childhood. J. Speech, Lang. Hear. Res., 62, 3381–3396.
    
    

101.101.Verhoef, E., Demontis, D., Burgess, S., et al. (2019) Disentangling polygenic associations between attention-deficit/hyperactivity disorder, educational attainment, literacy and language. Transl. Psychiatry, 9.
    
    

102.102.Vanoevelen, J., Dode, L., Van Baelen, K., et al. (2005) The secretory pathway Ca2+/Mn2+-ATPase 2 is a Golgi-localized pump with high affinity for Ca2+ ions. J. Biol. Chem., 280, 22800–22808.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODAvMjQvMjI4MDAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMS8xNS8yMDIxLjAxLjEyLjIxMjQ5MzE1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

103.103.Guillen Sacoto, M. J., Tchasovnikarova, I. A., Torti, E., et al. (2020) De Novo Variants in the ATPase Module of MORC2 Cause a Neurodevelopmental Disorder with Growth Retardation and Variable Craniofacial Dysmorphism. Am. J. Hum. Genet., 107, 352–363.
    
    

104.104.Kaplanis, J., Samocha, K., Wiel, L., et al. (2018) Integrating healthcare and research genetic data empowers the discovery of 28 novel developmental disorders..
    
    

105.105.Reuter, M. S., Tawamie, H., Buchert, R., et al. (2017) Diagnostic yield and novel candidate genes by exome sequencing in 152 consanguineous families with neurodevelopmental disorders. JAMA Psychiatry, 74, 293–299.
    
    

106.106.Snijders Blok, L., Rousseau, J., Twist, J., et al. (2018) CHD3 helicase domain mutations cause a neurodevelopmental syndrome with macrocephaly and impaired speech and language. Nat. Commun., 9, 1–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-02974-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29317637&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

107.107.Zheng, J. Q. and Poo, M.-M. (2007) Calcium Signaling in Neuronal Motility. Annu. Rev. Cell Dev. Biol, 23, 375–404.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.cellbio.23.090506.123221&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17944572&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250896200015&link_type=ISI) 

108.108.Dash, P. K., Moore, A. N., Kobori, N., et al. (2007) Molecular activity underlying working memory. Molecular activity underlying working memory. Learn. Mem. (2007), 14, 554–563.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoibGVhcm5tZW0iO3M6NToicmVzaWQiO3M6ODoiMTQvOC81NTQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMS8xNS8yMDIxLjAxLjEyLjIxMjQ5MzE1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

109.109.Nikoletopoulou, V. and Tavernarakis, N. (2012) Calcium homeostasis in aging neurons. Calcium homeostasis in aging neurons. Front. Genet. (2012), 3, 200.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fgene.2012.00200&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23060904&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

110.110.Vangheluwe, P., Sepúlveda, M. R., Missiaen, L., et al. (2009) Intracellular Ca 2+-and Mn 2+-transport ATPases. Chem. Rev., 109, 4733–4759.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/cr900013m&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19678701&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

111.111.He, W. and Hu, Z. (2012) The role of the golgi-resident SPCA Ca 2+/Mn 2+ pump in ionic homeostasis and neural function. The role of the golgi-resident SPCA Ca 2+/Mn 2+ pump in ionic homeostasis and neural function. Neurochem. Res. (2012), 37, 455–468.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11064-011-0644-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22083668&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

112.112.Normandin, L. and Hazell, A. S. (2002) Manganese neurotoxicity: An update of pathophysiologic mechanisms. Metabolic Brain Disease, Vol. 17, pp. 375–387.
    
    

113.113.Harischandra, D. S., Ghaisas, S., Zenitsky, G., et al. (2019) Manganese-induced neurotoxicity: New insights into the triad of protein misfolding, mitochondrial impairment, and neuroinflammation. Manganese-induced neurotoxicity: New insights into the triad of protein misfolding, mitochondrial impairment, and neuroinflammation. Front. Neurosci. (2019), 13, 654.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnins.2019.00654&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31293375&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

114.114.Wojda, U., Salinska, E. and Kuznicki, J. (2008) Calcium ions in neuronal degeneration. IUBMB Life, 60, 575–590.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/iub.91&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18478527&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258893900003&link_type=ISI) 

115.115.Xiang, M., Mohamalawari, D. and Rao, R. (2005) A novel isoform of the secretory pathway Ca2+,Mn(2+)-ATPase, hSPCA2, has unusual properties and is expressed in the brain. J. Biol. Chem., 280, 11608–14.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODAvMTIvMTE2MDgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMS8xNS8yMDIxLjAxLjEyLjIxMjQ5MzE1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

116.116.Ishikawa, K.-I., Nagase, T., Suyama, M., et al. (1998) Prediction of the Coding Sequences of Unidentified Human Genes. X. The Complete Sequences of 100 New cDNA Clones from Brain Which Can Code for Large Proteins in vitro. DNA Res., 5, 169–176.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/dnares/5.3.169&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9734811&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.12.21249315.atom) 

117.117.Aledo, J. C. (2019) Methionine in proteins: The Cinderella of the proteinogenic amino acids. Methionine in proteins: The Cinderella of the proteinogenic amino acids. Protein Sci. (2019), 28, 1785–1796.
    
    

118.118.Bozzi, A. T., Bane, L. B., Weihofen, W. A., et al. (2016) Conserved methionine dictates substrate preference in Nramp-family divalent metal transporters. Proc. Natl. Acad. Sci. U. S. A., 113, 10310–10315.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTEzLzM3LzEwMzEwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDEvMTUvMjAyMS4wMS4xMi4yMTI0OTMxNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)