Skip to main content

Quantitative genome-wide association analyses of receptive language in the Danish High Risk and Resilience Study

Abstract

Background

One of the most basic human traits is language. Linguistic ability, and disability, have been shown to have a strong genetic component in family and twin studies, but molecular genetic studies of language phenotypes are scarce, relative to studies of other cognitive traits and neurodevelopmental phenotypes. Moreover, most genetic studies examining such phenotypes do not incorporate parent-of-origin effects, which could account for some of the heritability of the investigated trait. We performed a genome-wide association study of receptive language, examining both child genetic effects and parent-of-origin effects.

Results

Using a family-based cohort with 400 children with receptive language scores, we found a genome-wide significant paternal parent-of-origin effect with a SNP, rs11787922, on chromosome 9q21.31, whereby the T allele reduced the mean receptive language score by ~ 23, constituting a reduction of more than 1.5 times the population SD (P = 1.04 × 10−8). We further confirmed that this association was not driven by broader neurodevelopmental diagnoses in the child or a family history of psychiatric diagnoses by incorporating covariates for the above and repeating the analysis.

Conclusions

Our study reports a genome-wide significant association for receptive language skills; to our knowledge, this is the first documented genome-wide significant association for this phenotype. Furthermore, our study illustrates the importance of considering parent-of-origin effects in association studies, particularly in the case of cognitive or neurodevelopmental traits, in which parental genetic data are not always incorporated.

Background

Language is one of the most fundamental human traits, and, as such, has been studied from many different angles, in disciplines such as philosophy, linguistics, psychology, anthropology, neuroscience and, more recently, genetics.

In the previous century, many studies showed that linguistic traits were heritable, i.e. that the genetic variation across individuals can explain some of the phenotypic (exhibited trait, in this case, linguistic ability) variation in the studied population. The investigated traits in those studies were either manifestations of atypical language development or of specific language skills in the general population [1], although the genetic factors influencing the two do not have to be the same. These insights were obtained through various study designs, such as twin studies and pedigree studies, which did not examine molecular genetic data.

With advances in molecular genetic techniques and the establishment of samples with available language-related phenotypes, studies were performed whose goal was to identify genes (or genetic variants) associated with linguistic ability. It was a series of studies of one extended pedigree which ultimately reported the first gene implicated in a speech and language disorder, FOXP2 [2,3,4]. Despite the fact that the disorder involving FOXP2 was monogenic (caused by a disruption of a single gene), the precise mechanism through which this gene caused the speech and language disorder is still unknown, and many subsequent studies examined the gene’s molecular function, its molecular evolution and its involvement in vocal communication in other species [5]. The disorder present in the aforementioned extended pedigree was rare and monogenic; non-synonymous variants in FOXP2 were not the cause of a common form of language impairment, namely specific language impairment (SLI) [6]. SLI is diagnosed when language development is deficient in the absence of known neurological or (other) neurodevelopmental abnormities [7]. This includes, for example, a discrepancy between verbal and non-verbal intelligence. However, in recent years, the diagnostic criteria have been changed, and more focus is given to the clinical picture of the child. As a consequence, a new label was adopted: developmental language disorder [8]. A discrepancy between verbal and non-verbal intelligence is not required for the latter diagnosis.

The first genome-wide screen for susceptibility loci for SLI was a non-parametric linkage analysis performed by the SLI Consortium [9], and it was later extended by the addition of more families to the cohort [10]. These studies identified a region on chromosome 16 linked to nonword repetition skills (reflecting phonological short-term memory) and a region on chromosome 19 linked to expressive language skills, in children with SLI. While a receptive language score was one of the quantitative phenotypes examined in the SLI Consortium linkage screen, no chromosomal regions were found to be significantly linked to it. However, in a later multivariate linkage analysis, receptive language traits did contribute to linkage signals e.g. on chromosomes 19 and 10 [11]. Another linkage study identified a region associated with reading impairment in children with SLI on chromosome 13 [12]; however, unlike the SLI Consortium studies, this study used parametric linkage analysis in large pedigrees. The same locus was investigated further in a following study [13]. Using an isolated population with a high rate of language impairment, linkage analyses under several models highlighted a region on chromosome 7 [14]. Regions on chromosomes 10 and 13 were identified in a more recent linkage study [15].

With the advent of dense marker maps and single-nucleotide polymorphism (SNP) arrays towards the end of the 2000s, and in relation to theoretical considerations in complex disease studies, disease-gene mapping strategies in the field of human genetics shifted from linkage studies towards genome-wide association studies (GWAS) [16]. In 2013, two GWAS were published which reported suggestive SNP associations with language-related traits in the general population [17, 18]. In 2014, the first GWAS of SLI reported a significant association exhibiting paternal parent-of-origin effects in a locus on chromosome 14 [19]. The same locus was later associated with reading-related traits, also displaying parent-of-origin effects in two dyslexia cohorts, albeit with discordant trends between the phenotypes [20]. A different approach, also employed in a GWAS design, was to test for association with principal components generated from both language and reading traits [21]. A receptive language phenotype was also investigated in a previous GWAS, although parent-of-origin effects were not modeled [22]. Expressive vocabulary in toddlers was also included in a GWAS design, which found a genome-wide significant association [23]. Complementary to genome-wide studies, other studies examined specific regions and/or variants following prior evidence of relevance to language traits. This included targeted association analyses of linkage regions [24], an exome sequencing study of the aforementioned isolated population [25] and an immunogenetic study of SLI [26], the latter of which also examined parent-of-origin effects.

A parent-of-origin effect denotes a change in the way an allele may influence a trait dependent on which parent it was inherited from. There are good reasons to examine parent-of-origin effects in general, and in studies of language-related traits and disorders in particular: firstly, they may account for some of the so-called “missing heritability” (the difference between disease heritability estimates and the variance explained by significant GWAS associations across published studies) [27]. A recent study illustrated that these effects can be missed in traditional designs, and that, in some cases, the same variant may have opposing effects on a given trait when inherited paternally or maternally [28]. A later study extended this observation to several quantitative traits [29]. Of note, this effect has also been observed in an association between a specific HLA-B allele and receptive language in SLI [26]. Secondly, parent-of-origin effects have been reported for several neurodevelopmental and psychiatric conditions [30]. Importantly, the latter include SLI [19, 26], as previously mentioned, and Williams syndrome [31], both of which involve abnormal development singling out language; while SLI involves impaired language development with other domains relatively unaffected, Williams syndrome is typically presented as involving a cognitive impairment which does not include linguistic ability (people with Williams syndrome are, in fact, hypersocial and hyperverbal), or, in other words, as a “mirror image” of SLI [32], although this view is, perhaps, too simplistic or inaccurate [33,34,35]. A recent family-based GWAS of autism spectrum disorder (ASD) highlighted 18 parent-of-origin effects, but the model employed in that study additionally included child genetic effects and maternal effects, all modeled simultaneously [36].

At this point it would be useful to discuss the expressive language-receptive language divide, as several of the aforementioned studies examined or found significant results for only one type of language trait. Expressive language refers to language production, and receptive language refers to language comprehension. An illustration of what can be measured by scores of expressive or receptive language ability can be found in the Clinical Evaluation of Language Fundamentals [37], which was used by the SLI Consortium in their studies. Expressive language subtests include, for example, formulating sentences (the child needs to formulate a sentence based on visual stimuli) and recalling sentences (the child needs to imitate a sentence produced first by the examiner), while receptive language subtests include, for example, semantic relationships (the child needs to answer questions correctly based on a spoken sentence) and oral directions (the child needs to point to pictured objects in response to directions). While it is easy to outline the different domains of expressive and receptive language, the distinction is not necessarily straightforward when considering real observations. For example, expressive and receptive language scores were found to be highly correlated in the SLI Consortium cohort [9]. Another issue in this regard is with the classification of children into groups based on the type of problems they display. For example, children “moved” from a cluster of expressive language disorder to expressive-receptive language disorder on follow-up only one year later [38]. Nonetheless, many standardized language tests maintain this distinction and provide separate scores for these domains. Leonard discusses this distinction in the context of SLI in more detail [39].

The above studies provided insight into the genetic basis of language, but genome-wide studies of language-related traits and language disorders remain scarce, relative to studies of other cognitive phenotypes and neurodevelopmental disorders. Studies in which parent-of-origin effects are modeled are even scarcer, as parental genotype data are not always available. In this study, we attempted to fill in some of these gaps by examining receptive language (which has not been the focus in follow-up studies to the first genome-wide screen of SLI) in a relatively large family cohort, and by examining parent-of-origin effects as well as child genetic effects in a GWAS design.

Results

In total, 400 children had Test for Reception of Grammar 2 (TROG-2) scores in our dataset. The distribution of the TROG-2 scores is shown in Fig. 1. The mean score across all 400 children was 101.2 with a standard deviation (SD) of 15.6, close to the population values/norms of 100 and 15, respectively. The distribution of these scores was not completely Gaussian. The means (SDs) within the ASD, attention deficit/hyperactivity disorder (ADHD) and index family subgroups were lower at 94.9 (± 16.8), 92.7 (± 17.3) and 99.8 (± 16.5, based on 243 children, as one child did not have a TROG-2 score), respectively. It should be mentioned that, in some cases, children who received an ASD diagnosis might have received it in part due to having low language ability, as mentioned in the Methods section; however, not all children with ASD had TROG-2 scores below the population mean for their age, and some indeed had scores above the population mean.

Fig. 1
figure1

The distribution of the standardized TROG-2 scores across all children in the study

In the discovery stage, we found a genome-wide significant association in the paternal parent-of-origin test; the top hits in the other two tests showed only suggestive association: rs12297354 in the general test (P = 7 × 10−6), and rs9397338 in the maternal parent-of-origin test (P = 8 × 10−6). Manhattan and QQ plots for the discovery analyses are shown in Additional file 1: Figure S1. Associations that are at least suggestive (P ≤ 10−5) from all three tests can be found in Additional file 2: Table S1.

The top hit in the paternal parent-of-origin test was with an intergenic SNP on chromosome 9:82,788,836, rs11787922 (Fig. 2). In our analysis, when examining alleles inherited from the father, the T allele (the minor allele), was associated with a reduction of 23.05 in the TROG-2 score (P = 1.04 × 10−8), which is more than 1.5 SD away from the population mean. The frequency of the T allele in our entire sample was 7.54%. It was 5.88%, 10.64% and 6.98% in the ASD, ADHD and index family children subgroups, respectively. The T allele was not genome-wide significantly associated in the general test (its effect was − 9.46, P = 9 × 10−6) or in the maternal parent-of-origin test (its effect was − 0.287, P = 0.927), but it showed suggestive association in the former.

Fig. 2
figure2

A regional association plot showing the top SNP in the paternal parent-of-origin analysis and surrounding SNPs and genes

To assess the evidence for a paternal association at the above locus further, we performed five post hoc tests. The first test assessed whether there was a significant difference between the maternal and paternal alleles; it is possible that, when looking only at paternally inherited alleles, the association is significant, and that, when looking at maternally inherited alleles, it is not (or vice versa). However, that alone does not mean that the difference itself between the associations with paternal and maternal alleles is significant. QTDT has a special test for that, and our first post hoc test confirmed that the difference was indeed significant in this case (P = 5 × 10−6). In our second post hoc test, we repeated the association test having added covariates for ASD, ADHD and for whether the family was an index family or a control family, as we had observed lower scores (on average) in children from all three groups compared to all children in the study. While the covariates explained some of the reduction in the TROG-2 score, the genetic effect was still large (effect of T allele on TROG-2 score = − 21.3, P = 3.25 × 10−8). When using a covariate for sex, the effect of T allele on TROG-2 score was − 22.9, (P = 1.13 × 10−8).

In the fourth post hoc test, where scores were transformed to approach normality and scaled, the association signal between rs11787922 and receptive language remained significant (P = 8.33 × 10−7), with the T allele being associated with a reduced score, but the effect was somewhat diminished at − 1.3 SD of the transformed score. Variance component models are sensitive to departure from normality; we used them in our model to account for relatedness in 11 families with a sibship of 2. If we remove one sibling from each of these families (keeping children with the T allele, if only one sibling had it, and otherwise randomly choosing one) and do not include variance components in the model, the results are similar to the original, with an effect of the T allele = − 22.9 (P = 2 × 10−8).

Lastly, we used an independent method of assessing parent-of-origin effects in unrelated individuals. Using this method, we observe an association signal for rs11787922 (P = 0.0016). This method relies on unrelated individuals and was designed with larger samples in mind, and we also lost power by excluding individuals. Nonetheless, even with this method, we observed some signal at this locus. Table 1 includes a summary of all post hoc tests.

Table 1 Results of the post hoc tests for rs11787922 (the original test is included as a point of reference)

Discussion

Our GWAS identified a genome-wide significant association for receptive language with a paternal parent-of-origin effect of a SNP on chromosome 9, rs11787922, supported by both a family-based association method and an independent method of assessing parent-of-origin effects which is based on the variance of the trait in unrelated individuals. This association was not driven by the presence of two other neurodevelopmental disorders in which language may be impaired, namely ASD and ADHD, nor was it driven by the parents’ having a diagnosis of schizophrenia or bipolar disorder, despite the fact that the receptive language scores across individuals with ASD, ADHD or whose parents had the above psychiatric diagnoses were lower on average compared to the entire cohort. This suggests that the association of this SNP with receptive language is not part of a larger behavioral phenotype encompassing ASD, ADHD, familial risk of schizophrenia or bipolar disorder. In this context, it would be useful to note that a study of polygenic risk which used the same cohort used in this study [40] found that a polygenic risk score (PRS) trained on a previous GWAS of SLI was not predictive of the risk of ASD or ADHD in this cohort, but it was predictive of the risk of SLI. The latter resulted in a Nagelkerke’s R2 (adjusted for prevalence and the proportion of cases in the sample) of > 5%, compared to values of close to 0% for ASD and ADHD (see footnote in the Methods section for more information about this). A different way of estimating additive genetic overlaps between disorders is estimating genetic correlations between them. This can be achieved in several ways e.g. by using family data or through examining genetic and phenotypic variation across unrelated individuals, the latter being more similar to the PRS analyses. Earlier studies, which examined genetic correlations across ASD and ADHD or specific traits associated with the two disorders, reported moderate to high significant genetic correlations [41, 42], but SNP-based methods did not obtain similar estimates: they were either lower [43] or non-significant [44].

To our knowledge, this is the first genome-wide significant association reported for receptive language; a previous GWAS of receptive language (which did not model parent-of-origin effects) found only suggestive associations [22].

Our top SNP has not been highlighted in past genetic studies of language or other cognitive traits. Therefore, we surveyed the literature and examined the region surrounding the SNP for possible involvement in relevant traits. Interestingly, the protein-coding gene that is closest to the SNP is TLE4 (Fig. 2), which is involved in neural development [45, 46]. The gene is involved in cell determination and differentiation in neurogenesis, together with other genes of the same family, working in different combinations of expression patterns during neurogenesis. TLE4, specifically, is expressed in neural progenitor cells in the early embryonic stage.

Through PhenoScanner [47], we found a proxy SNP for rs11787922, namely rs10512097 (r2 = 0.85 in Europeans, not included in our GWAS), which was a suggestive expression quantitative trait locus (eQTL) for IMP5 (also called SPPL2C) (P = 7.18 × 10−6) [48]. Deletions of this gene have been reported for some forms of intellectual disability and developmental delay [49, 50], although de novo loss-of-function variants in a nearby gene were enough to cause a similar phenotype [51]. Recently, SPPL2C has been highlighted in a G × E study which examined the influence of common variants on language through low-frequency hearing ability [52].

Our top SNP, rs11787922, lies in a region encompassed by several pathogenic copy-number variants (CNVs) (found with the ClinGen track on UCSC Genome Browser) associated with various forms of developmental delay, including a CNV identified in an individual with delayed speech and language development (ClinVar accession number: VCV000148901.1). This CNV, however, was quite large (chromosome 9:71,130,848-86,285,142) and contains many genes, and so it could, at most, provide suggestive evidence for the region around the SNP’s being associated with a language phenotype. A survey of the literature found that the chromosomal band encompassing the SNP on chromosome 9 (i.e. 9q21.31), as well as adjacent bands, were implicated in several studies of phenotypes such as developmental delay, language delay and intellectual disability. At times these were part of a broad syndrome, namely Gorlin syndrome, which typically involves 9q22 and the PTCH1 gene [53], but there were cases of individuals with CNVs in that region, not encompassing PTCH1, that were still associated with language problems, including associations with receptive language [54]. A study of nine individuals of whom eight showed severe speech and language delay identified deletions on 9q21 (mean size = 7.14 Mbp), of which six involved 9q21.31 [55]. Among the genes deleted across all deletion cases in that study were RORB, TRPM6, NMRK1 and OSTF1. Of note, RORB was reported to be involved in verbal intelligence [56].

While there seems to be some evidence for the involvement of the locus on chromosome 9 in language-related phenotypes, one further point to consider is that our top association displayed a paternal parent-origin effect. Information regarding imprinted regions in the human genome is not readily available in online databases (imprinting, or silencing of one parental allele, being one of the major mechanisms underlying parent-of-origin effects). A “Catalogue of Imprinted Genes and Parent-of-origin Effects” exists [57], but it is quite small and has not been updated since 2016. Another database, Geneimprint [58], also contains a small number of entries, and they are limited to specific genes only. However, it is sometimes possible to find information in cytogenetic studies which could provide some insight regarding imprinting of genomic regions. For example, for 9q22, some deletions (not all involving PTCH1) associated with language delay or intellectual disability are of paternal origin, suggesting that the paternal allele is expressed in that region [59, 60]. This does not mean that the maternal allele is not expressed. Further insight can be gained by considering uniparental disomy (UPD) cases, i.e. scenarios in which both chromosomes of the same kind (or chromosomal segments) in the child come from one parent. If a region is paternally or maternally imprinted, then a UPD involving it may result in an abnormal phenotype, either due to having too much gene expression from two non-imprinted chromosomes, or insufficient gene expression, if both chromosomes are imprinted. One study [61] reported the case of a child with neurodevelopmental problems who had a paternal UPD from 9q21.33 onwards (a triplication of 9q21.11–q21.33, a region which included our top SNP, was also observed). This UPD involved inheriting two copies of a segment of one of the homologous chromosome 9’s from the father (isodisomy), which means that the phenotype may potentially be caused by either homozygosity or imprinting. Thus, there is some indirect evidence suggesting that the implicated region on chromosome 9 may be imprinted. It should also be noted that loci may display parent-of-origin effects if they interact with imprinted regions or genes, even if they are not imprinted themselves [62], which implies that even if the top SNP in our analyses is not in an imprinted region, it might still display parent-of-origin effects through interaction with imprinted genes or loci elsewhere in the genome.

Some of the suggestive associations (Additional file 2: Table S1) are also worth mentioning. For example, rs4632856 (P = 4 × 10−6 in the paternal parent-of-origin test) was shown to be an eQTL for PHIP [63] (P = 2.36 × 10−8), a gene which has been implicated in intellectual disability and developmental delay [64, 65]; rs3731769 (P = 0.00001 in the maternal parent-of-origin test) was suggestive (P = 4.95 × 10−6) in a study of autism [66]. Interestingly, it exhibited a paternal parent-of-origin effect. This opposite parental effect is similar to the one observed across the SLI and dyslexia studies, as mentioned earlier.

Our study had several limitations. Firstly, the genotyping array used in this study included both common and rare variants, the latter having been identified through studies of psychiatric disorders or included from an exome chip. Unfortunately, due to these variants being rare, most of them were filtered out for having low minor allele frequencies during QC; nearly half of the genotyped variants were removed and thus not included in downstream analyses. This also means that the genome-wide coverage is lower than in studies which used genotyping arrays with more common variants. Also related to this is the fact that we used software not designed to work with imputed dosage data or extremely large marker datasets. Secondly, cohorts with data on receptive language and genotyped parents (both mothers and fathers, to be able to test for both parent-of-origin effects and the difference between them) are scarce. Specifically, even though samples meeting those criteria do exist, they are often ascertained for having language disorders or related neurodevelopmental conditions, meaning that the receptive language scores are already likely to be quite below the population mean, which could potentially mask an effect observed in non-ascertained samples, and/or they may also comprise fewer families, as is the case for some of the previously mentioned studies. Our association thus requires replication in a suitable, independent sample. Lastly, the effect size of the top SNP in our study is quite large (the T allele was associated with a reduction of 23.05 in the TROG-2 score). In the general test of association (without parent-of-origin effects), the effect of the T allele was a reduction of 9.46 in the TROG-2 score, and it was not genome-wide significant. While we could not find relevant examples of studies employing parent-of-origin tests, the result in the general test is of the same order of magnitude (i.e. 1 < X < 10 for standardized scores with a mean of 100 and a SD of 15 in some studies, or 0.5 SD < X < 1 SD in other studies) as observed in associations between CNTNAP2 variants and language traits in children with SLI [67], intelligence in a general population twin cohort [68], and handedness in individuals with dyslexia [69]. The fact that this effect appears smaller and less significant in the general test, assuming a true parent-of-origin effect is present, is in line with the results of the aforementioned studies which investigated parent-of-origin effects across several quantitative traits. However, we cannot rule out an overestimation of this effect or the winner’s curse [70, 71].

Conclusions

We report a genome-wide significant association with receptive language in a sample not selected for language impairment. The association displays a parent-of-origin effect, whereby the minor allele reduces the language test score by more than 1.5 times the population SD, when inherited from the father. Our results contribute to the scientific literature and shed light on a relatively understudied research topic–the genetic underpinnings of linguistic ability. Moreover, our study illustrates the importance of considering parent-of-origin effects in genetic studies, which are often not modeled, in part due to constraints in data availability. In the same vein, we hope that our study will encourage further investigations into this and similar phenotypes.

Methods

Participants

The cohort used in this study comes from the Danish High Risk and Resilience Study–the VIA 7 study [72]. The VIA 7 study recruited children around age 7 and their biological parents. Families were selected on account of having one or two parents with a diagnosis of either schizophrenia spectrum psychosis or bipolar disorder (hereafter index families) or as control families, in which neither parent had schizophrenia or bipolar disorder (hereafter control families). NB: in other VIA 7 publications “index” may refer to the biological parents themselves who were selected for having a psychiatric diagnosis as well as to matched control parents (the non-index parent being the other biological parent). However, here, “index” refers to the families of parents who have a psychiatric diagnosis, and “control” refers to families of parents who do not have a psychiatric diagnosis. The families in the study were recruited from across Denmark. Just under half of the families lived in densely populated areas, with no significant differences in urbanization between index families (neither schizophrenia nor bipolar disorder) and control families. The parents in the schizophrenia group were younger and had lower levels of education compared to the other two groups. Both biological parents in the index families had lower levels of functioning and were less often in employment or enrolled in a study program [73]. In the genetic analyses, we removed individuals of non-European ancestry and their relatives, as explained below. The children were administered numerous tests and interviews, focusing on cognitive, behavioral, social and psychomotor measures and assessed for psychopathology.

A subset of the families in VIA 7 provided DNA samples, which were subsequently genotyped on the Illumina PsychChip v1.1. This subset was used in this study. Following genetic quality control (QC), individuals from 429 families remained, consisting of child-parents trios, child-parent duos, or only children, and, in a minority of families, a sibling as well (also, occasionally children might have been removed due to QC or not genotyped, in which case the parents might have been kept in the pedigree if they themselves passed QC, e.g. for more accurate allele frequency calculations). Due to the nature of the analyses, the precise number of children used in the association test with a given marker may differ across markers, as explained later in this section, but we note the precise number of informative children in the test for the top association observed and for the suggestive associations in Additional file 2: Table S1.

Language data and other phenotypic measures

The test used in this study was the Test for Reception of Grammar (TROG-2) [74]. This test measures receptive language skills by presenting children with blocks of four sets of pictures. In each set of pictures, only one picture fully corresponds to a sentence uttered by the examiner. The child must choose the correct picture in every set to have “passed” the block correctly. Scores comprising the number of correct blocks were age-standardized according to the norms from the Danish manual. Children also underwent screening with the Danish version of the Schedule for Affective Disorders and Schizophrenia for School-Age Children (K-SADS) [75]. Children who received a probable or a definite diagnosis for ASD or ADHD based on the K-SADS were defined as cases for their respective diagnosis/diagnoses (children who did not have a probable or a definite diagnosis were defined as controls for each diagnosis separately) [76]. All diagnostic assessments of children who were suspected of having ASD or ADHD were reviewed by a specialist in child and adolescent psychiatry, and data from other measures available in the VIA 7 study might have been taken into consideration when making the diagnosis, including data on language ability. We used these indications to construct covariates in a post hoc test as described later. Out of the entire sample of children used in this study, 17 were diagnosed with ASD, 47 were diagnosed with ADHD and 244 came from an index family, as defined above.

Genetic data and quality control

The QC steps for the genetic data are described in detail in a recent study of PRS for language impairment [40]. QC was performed with PLINK v1.9b5.2 [77]. The main differences between the dataset used in the PRS study and the one used in the present study are (i) that the present study did not exclude parents or siblings from the analyses, as the analyses were family-based, unlike the PRS analyses, (ii) that markers in the major histocompatibility region were not removed, and (iii) that duplicate markers (based on position) were removed (in the PRS analyses these would have been removed during the clumping procedure, if still present after earlier filtering steps and in the unlikely event that the same duplicate set of markers was included in both the base and target datasets). Due to the requirements of the software used in the present study, “dummy” individuals (i.e. individuals for whom genetic data were not available) were added to the pedigree so that all children had two parents in each nuclear family. Since by default PLINK checks for Mendelian errors only in trios, this resulted in a small number of additional errors found in some child-parent duos. No child-parent duos not otherwise removed during other QC steps showed an excess of errors large enough (> 1%) to have been removed for this alone, and only a small number of markers that were not removed previously were flagged for excess errors. While we confirmed that this issue did not have any major impact on the results of the PRS analyses in the previous study,Footnote 1 we chose to remove the offending genotypes and markers from the present study, as it examined one marker at a time. Apart from the individuals removed for excess Mendelian errors (N = 10), individuals and markers were also removed following the other main steps, which are repeated here in brief form for easy reference: initial QC on raw genetic data: individuals with low call rates or discordant sex information were removed in the first step (N = 18), as were SNPs with a Gentrain score < 0.3. QC with PLINK: SNPs with > 5% missing data were removed (all remaining individuals had < 5% missing data). Individuals with extreme heterozygosity rates (with a threshold of ± 3 SD from the sample mean) were removed (N = 21). Genetic ancestry was estimated in a principal component analysis (PCA). The threshold for the exclusion of samples was 2 SD above or below the VIA 7 mean for either PC1 or PC2, using the VIA 7 samples and the CEU, CHB, JPT and YRI HAPMAP samples to create the PC space [78]. Individuals of divergent ancestry were removed along with their relatives (N = 36, the remaining individuals clustered near the European (CEU) reference individuals), as were individuals who exhibited cryptic relatedness (the Pi-hat threshold for the exclusion of individuals not known to be related was 0.185) or who were less related to family members than expected from pedigree information (N = 13). A Hardy–Weinberg Equilibrium P value threshold of 1 × 10−6 was employed for QC-passing SNPs, as well as a minor allele frequency threshold of 1%. In total, 299,604 autosomal SNPs passed all QC steps. Positions in the text and tables are in hg19.

Discovery analyses

The software package QTDT v2.6.1 [79] was used for conducting the statistical genetic analyses (QTDT stands for quantitative transmission-disequilibrium test, although the tests employed here were not TDTs, as explained below). MERLIN v1.1.2 [80] was used for estimating identity by descent scores for each marker to be used by QTDT. Three types of analyses were performed: a total test of association using all family data (qtdt -at), a paternal parent-of-origin analysis (qtdt -at -op) and a maternal parent-of-origin analysis (qtdt -at -om). The “total association” model was used for all the above tests, as it is more powerful in the absence of population stratification (individuals of divergent ancestry had been removed from our dataset in the PCA step during QC). In this model, a combined between/within family component “X”, denoting the between/within effect on the means, is tested (X is the effect size reported for the QTDT analyses). In the null model this component is fixed to 0, and in the full model it is estimated from the data. The likelihoods of these two models are then compared resulting in a χ2 statistic (in most cases). We included variance components in all models (-wega), incorporating an environmental component, a polygenic component and an additive major locus component. This allowed for the use of families with multiple children. In the two parent-of-origin analyses, the same model was used, except that only paternally inherited alleles and maternally inherited alleles were used in the paternal parent-of-origin analysis and maternal parent-of-origin analysis, respectively. The phenotype used in the discovery analyses was the standardized TROG-2 score. LocusZoom [81] was used to plot the region surrounding the top association.

The power of a QTDT analysis

The power of a QTDT analysis depends on factors such as: the marker allele frequencies, the effect size, the linkage disequilibrium between the marker and the quantitative trait locus, the number of genotypes in the analysis (determined by the number of children/siblings), whether parental genotypes are included or not, and the significance level. Power analyses for QTDT are carried out through simulation and were performed in several studies focusing on the method itself or on how it compares with other association methods. In the original QTDT paper, assuming a maximum D’, h2 of 0.1, risk allele frequency of 0.5, significance level of 0.001 and including parental genotypes, a sample of 480 children (families with a sibship of 1) resulted in a power estimate of 97.4% [79]. Another study reported a power of 74% with N = 200, h2 of 0.1, and a risk allele frequency of 0.3 [82]. However, it is important to keep in mind two main ways in which our analyses differ from these studied examples: (i) studies which looked into the power of a typical QTDT analysis employed the orthogonal model and not the total association model, the latter of which was used in our study, and (ii) they did not examine parent-of-origin models. These two factors are expected to suggest that our analyses should have greater power, all other things being equal, because: (i) the power to detect a parent-of-origin effect in trios tends be higher than the power to detect a child genetic effect, given equal sample sizes, if the model is correct [83], and (ii) in the absence of population stratification, the QTDT total association model has greater power than the orthogonal model [84].

Post hoc tests

For the top association from the discovery analyses, we performed five additional tests: four tests were performed with QTDT as well: the first test checked whether the difference between the paternal and maternal alleles was significant (qtdt -at -ot).

The second test repeated the original analysis in which the top SNP was highlighted, incorporating covariates for a diagnosis of ASD, a diagnosis of ADHD, and whether the family was an index family or a control family. These were included as “dummy variables” having a value of 0 (control for ASD, ADHD/control family) or 1 (case for ASD, ADHD/index family). We also repeated the test with a covariate for sex (third test). This was done separately due to the sex bias observed in ASD and ADHD (to avoid correlation between the covariates).

The fourth test repeated the original analysis but used Z-scores of transformed TROG-2 standardized scores (Z-scores were used, because the transformed scores could have extreme values). The method for the transformation, the Yeo-Johnson transformation [85], and the parameter lambda = 2.762526, were selected with the bestNormalize package v1.4.3 [86] and employed with the VGAM package v1.1-2 [87] in R v3.6.3 [88]. This test was performed because QTDT assumes normality of the phenotype scores, and the distribution of the scores in our sample was not completely Gaussian. We therefore transformed the scores to approach normality. It is worth noting that the original QTDT method paper as well as several publications mentioning the assumption of normality and/or the effect of violating it discussed in this context the orthogonal model presented by Abecasis et al. in the QTDT paper [79, 89, 90]. As previously mentioned, the model we used was the total association model. In fact, the QTDT program provides a way to account for non-normality using a permutation procedure, but this is not implemented for the total association model in QTDT. Some differences between the two models should be noted: the total association model is not a TDT and makes use of more children (the difference in sample size between the orthogonal model and the total association model may be quite large, in favor of the latter). In particular, in addition to this model using information from more children, the total association model for parent-of-origin effects specifically is less constricted in terms of parental genotypes: the orthogonal parent-of-origin model uses children whose both parents are genotyped and where one parent is homozygous, or whose mother and father have different genotypes. In addition, when paternal parent-of-origin effects are tested, the father must be heterozygous and, when maternal effects are tested, the mother must be heterozygous. The total association parent-of-origin model uses, in addition to the above group of children, all children with at least one homozygous parent, even if the other parent has a missing genotype [91]. With regards to the top SNP in our analyses, the orthogonal model could use only 23 children (fewer than the default minimal number of informative children), whereas the total association model used 353 children, making it much more powered. This also illustrates the fact that individuals, and thus phenotype distributions, could vary greatly for each tested marker, especially when testing for parent-of-origin effects. Nonetheless, since the effect of a violation of normality in the total association model is not known, we decided to repeat the test following the transformation of the scores across all children. It should be mentioned that transforming trait values does not necessarily improve the way in which the model behaves and may limit the interpretation of the results [92,93,94], but we thought it best to include both estimates in the paper for the sake of completeness.

Additionally, we employed a newer approach for the detection of parent-of-origin effects for quantitative traits in unrelated individuals (fifth test). To this end we used the “POE method” as implemented in QUICKTEST v0.99b [95]. In this analysis, 391 unrelated children were used, the same sample as used in the PRS study [40] (TROG-2 scores were available for 389 children of those 391 unrelated children). The underlying principle of this approach is the following: assuming that a locus with alleles A and B has an additive effect on a phenotype (the two alleles have different effects on the quantitative trait), the three genotype groups AA, AB and BB will differ in their mean phenotype scores. If a parent-of-origin effect is present, then the heterozygous group (individuals with genotype AB) will consist of two subgroups, namely individuals with paternal A and maternal B and individuals with maternal A and paternal B, which should have different means, thereby increasing the variance of the heterozygous group as a whole. If there is no parent-of-origin effect at this locus, then the two heterozygous subgroups should not be distinguishable in this respect. This method provides a formal test for this by comparing the variance of the heterozygous group to those of the homozygous groups. It should be noted that his method does not indicate what type of parental effect (i.e. paternal or maternal) is present, because parental genotypes are not available to it; rather, it provides statistical evidence for the presence of such an effect.

Availability of data and materials

Access to the dataset used in the current study are available from the corresponding authors upon reasonable request.

Notes

  1. 1.

    For an example of this effect: the adjusted Nagelkerke’s R2 for the model for the SLI phenotype was 5.21% (with P < 0.05), for ASD it was 0.015%, and for ADHD it was 0.009%. The R2 for height was 0.019%; following the removal of some SNPs and Mendelian errors in duos, the values were 5.09% (with P < 0.05), 0.017%, 0.009% and 0.01%, respectively (using the dataset termed Correction 1 in the paper).

Abbreviations

ADHD:

Attention deficit/hyperactivity disorder

ASD:

Autism spectrum disorder

CNV:

Copy-number variant

eQTL:

Expression quantitative trait locus

GWAS:

Genome-wide association study

K-SADS:

Schedule for Affective Disorders and Schizophrenia for School-Age Children

PCA:

Principal component analysis

PRS:

Polygenic risk score

QC:

Quality control

QTDT:

Quantitative transmission-disequilibrium test

SD:

Standard deviation

SLI:

Specific language impairment

SNP:

Single-nucleotide polymorphism

TROG-2:

Test for Reception of Grammar 2

UPD:

Uniparental disomy

References

  1. 1.

    Stromswold K. The heritability of language: a review and metaanalysis of twin, adoption, and linkage studies. Language. 2001;77(4):647–723.

    Google Scholar 

  2. 2.

    Hurst JA, Baraitser M, Auger E, Graham F, Norell S. An extended family with a dominantly inherited speech disorder. Dev Med Child Neurol. 1990;32(4):352–5 Epub 1990/04/01.

    CAS  PubMed  Google Scholar 

  3. 3.

    Fisher SE, Vargha-Khadem F, Watkins KE, Monaco AP, Pembrey ME. Localisation of a gene implicated in a severe speech and language disorder. Nat Genet. 1998;18(2):168–70 Epub 1998/02/14.

    CAS  PubMed  Google Scholar 

  4. 4.

    Lai CS, Fisher SE, Hurst JA, Vargha-Khadem F, Monaco AP. A forkhead-domain gene is mutated in a severe speech and language disorder. Nature. 2001;413(6855):519–23 Epub 2001/10/05.

    CAS  PubMed  Google Scholar 

  5. 5.

    Nudel R, Newbury DF. Foxp2. Wiley interdisciplinary reviews Cognitive science. 2013;4(5):547-60. Epub 2014/04/26.

  6. 6.

    Newbury DF, Bonora E, Lamb JA, Fisher SE, Lai CS, Baird G, et al. FOXP2 is not a major susceptibility gene for autism or specific language impairment. Am J Hum Genet. 2002;70(5):1318–27 Epub 2002/03/15.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Bishop DVM. What causes specific language impairment in children? Curr Dir Psychol Sci. 2006;15(5):217–21 Epub 2008/11/15.

    PubMed  PubMed Central  Google Scholar 

  8. 8.

    Bishop DVM, Snowling MJ, Thompson PA, Greenhalgh T. Phase 2 of CATALISE: a multinational and multidisciplinary Delphi consensus study of problems with language development: terminology. J Child Psychol Psychiatry. 2017;58(10):1068–80 Epub 2017/04/04.

    PubMed  PubMed Central  Google Scholar 

  9. 9.

    The SLI Consortium. A genomewide scan identifies two novel loci involved in specific language impairment. Am J Hum Genet. 2002;70(2):384–98 Epub 2002/01/16.

    PubMed Central  Google Scholar 

  10. 10.

    The SLI Consortium. Highly significant linkage to the SLI1 locus in an expanded sample of individuals affected by specific language impairment. Am J Hum Genet. 2004;74(6):1225–38 Epub 2004/05/11.

    Google Scholar 

  11. 11.

    Monaco AP. The SLI consortium. Multivariate linkage analysis of specific language impairment (SLI). Ann Human Genet. 2007;71(5):660–73.

    CAS  Google Scholar 

  12. 12.

    Bartlett CW, Flax JF, Logue MW, Vieland VJ, Bassett AS, Tallal P, et al. A major susceptibility locus for specific language impairment is located on 13q21. Am J Hum Genet. 2002;71(1):45–55 Epub 2002/06/06.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Bartlett CW, Flax JF, Logue MW, Smith BJ, Vieland VJ, Tallal P, et al. 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. 2004;57(1):10–20.

    PubMed  PubMed Central  Google Scholar 

  14. 14.

    Villanueva P, Newbury DF, Jara L, De Barbieri Z, Mirza G, Palomino HM, et al. Genome-wide analysis of genetic susceptibility to language impairment in an isolated Chilean population. Eur J Human Gene. 2011;19(6):687–95 Epub 2011/01/21.

    CAS  Google Scholar 

  15. 15.

    Evans PD, Mueller KL, Gamazon ER, Cox NJ, Tomblin JB. A genome-wide sib-pair scan for quantitative language traits reveals linkage to chromosomes 10 and 13. Genes Brain Behav. 2015;14(5):387–97 Epub 2015/05/23.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JP, et al. Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008;9(5):356–69 Epub 2008/04/10.

    CAS  PubMed  Google Scholar 

  17. 17.

    Luciano M, Evans DM, Hansell NK, Medland SE, Montgomery GW, Martin NG, et al. A genome-wide association study for reading and language abilities in two population cohorts. Genes Brain Behav. 2013;12(6):645–52 Epub 2013/06/07.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Eicher JD, Powers NR, Miller LL, Akshoomoff N, Amaral DG, Bloss CS, et al. Genome-wide association study of shared components of reading disability and language impairment. Genes Brain Behav. 2013;12(8):792–801 Epub 2013/09/13.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Nudel R, Simpson NH, Baird G, O’Hare A, Conti-Ramsden G, Bolton PF, et al. Genome-wide association analyses of child genotype effects and parent-of-origin effects in specific language impairment. Genes Brain Behav. 2014;13(4):418–29 Epub 2014/02/28.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Pettigrew KA, Frinton E, Nudel R, Chan MTM, Thompson P, Hayiou-Thomas ME, et al. Further evidence for a parent-of-origin effect at the NOP9 locus on language-related phenotypes. J Neurodev Disorders. 2016;8:24 Epub 2016/06/17.

    Google Scholar 

  21. 21.

    Gialluisi A, Newbury DF, Wilcutt EG, Olson RK, DeFries JC, Brandler WM, et al. Genome-wide screening for DNA variants associated with reading and language traits. Genes Brain Behav. 2014;13(7):686–701 Epub 2014/07/30.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Harlaar N, Meaburn EL, Hayiou-Thomas ME, Davis OS, Docherty S, Hanscombe KB, et al. Genome-wide association study of receptive language ability of 12-year-olds. JSLHR. 2014;57(1):96–105 Epub 2014/04/02.

    PubMed  Google Scholar 

  23. 23.

    St Pourcain B, Cents RAM, Whitehouse AJO, Haworth CMA, Davis OSP, O’Reilly PF, et al. Common variation near ROBO2 is associated with expressive vocabulary in infancy. Nat Commun. 2014;5(1):4831.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Newbury DF, Winchester L, Addis L, Paracchini S, Buckingham LL, Clark A, et al. CMIP and ATP2C2 modulate phonological short-term memory in language impairment. Am J Hum Genet. 2009;85(2):264–72 Epub 2009/08/04.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Villanueva P, Nudel R, Hoischen A, Fernandez MA, Simpson NH, Gilissen C, et al. Exome sequencing in an admixed isolated population indicates NFXL1 variants confer a risk for specific language impairment. PLoS Genet. 2015;11(3):e1004925 Epub 2015/03/18.

    PubMed  PubMed Central  Google Scholar 

  26. 26.

    Nudel R, Simpson NH, Baird G, O’Hare A, Conti-Ramsden G, Bolton PF, et al. Associations of HLA alleles with specific language impairment. J Neurodev Disorders. 2014;6(1):1.

    Google Scholar 

  27. 27.

    Kong A, Steinthorsdottir V, Masson G, Thorleifsson G, Sulem P, Besenbacher S, et al. Parental origin of sequence variants associated with complex diseases. Nature. 2009;462(7275):868–74 Epub 2009/12/18.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Benonisdottir S, Oddsson A, Helgason A, Kristjansson RP, Sveinbjornsson G, Oskarsdottir A, et al. Epigenetic and genetic components of height regulation. Nature Commun. 2016;7:13490 Epub 2016/11/17.

    CAS  Google Scholar 

  29. 29.

    Mozaffari SV, DeCara JM, Shah SJ, Sidore C, Fiorillo E, Cucca F, et al. Parent-of-origin effects on quantitative phenotypes in a large Hutterite pedigree. Commun Biol. 2019;2:28 Epub 2019/01/25.

    PubMed  PubMed Central  Google Scholar 

  30. 30.

    Davies W, Isles AR, Wilkinson LS. Imprinted genes and mental dysfunction. Ann Med. 2001;33(6):428–36 Epub 2001/10/05.

    CAS  PubMed  Google Scholar 

  31. 31.

    Collette JC, Chen XN, Mills DL, Galaburda AM, Reiss AL, Bellugi U, et al. William’s syndrome: gene expression is related to parental origin and regional coordinate control. J Hum Genet. 2009;54(4):193–8 Epub 2009/03/14.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Tomblin JB. Children with specific language impairment. In: Bavin EL, editor. The Cambridge Handbook of Child Language. Cambridge: Cambridge University Press; 2009. p. 417–32.

    Google Scholar 

  33. 33.

    Stojanovik V, Perkins M, Howard S. Williams syndrome and specific language impairment do not support claims for developmental double dissociations and innate modularity. J Neurolinguistics. 2004;17(6):403–24.

    Google Scholar 

  34. 34.

    Laws G, Bishop D. Pragmatic language impairment and social deficits in Williams syndrome: a comparison with Down’s syndrome and specific language impairment. Int J Lang Commun Disorders Royal College Speech Lang Ther. 2004;39(1):45–64.

    Google Scholar 

  35. 35.

    Rice ML, Warren SF, Betz SK. Language symptoms of developmental language disorders: an overview of autism, Down syndrome, fragile X, specific language impairment, and Williams syndrome. Appl Psycholinguistics. 2005;26(1):7–27.

    Google Scholar 

  36. 36.

    Connolly S, Anney R, Gallagher L, Heron EA. A genome-wide investigation into parent-of-origin effects in autism spectrum disorder identifies previously associated genes including SHANK3. EJHG. 2017;25(2):234–9 Epub 2016/11/24.

    CAS  PubMed  Google Scholar 

  37. 37.

    Semel EM, Wiig EH, Secord W. Clinical evaluation of language fundamentals-revised. San Antonio: Phychol Corporat; 1992.

    Google Scholar 

  38. 38.

    Conti-Ramsden G, Botting N. Classification of children with specific language impairment: longitudinal considerations. JSLHR. 1999;42(5):1195–204 Epub 1999/10/09.

    CAS  PubMed  Google Scholar 

  39. 39.

    Leonard LB. Children with specific language impairment. 2nd ed. Cambridge: MIT press; 2014.

    Google Scholar 

  40. 40.

    Nudel R, Christiani CAJ, Ohland J, Uddin MJ, Hemager N, Ellersgaard DV, et al. Language deficits in specific language impairment, attention deficit/hyperactivity disorder, and autism spectrum disorder: an analysis of polygenic risk. Autism Res. 2020;13(3):369–81 Epub 2019/10/03.

    PubMed  Google Scholar 

  41. 41.

    Reiersen AM, Constantino JN, Grimmer M, Martin NG, Todd RD. Evidence for shared genetic influences on self-reported ADHD and autistic symptoms in young adult Australian twins. Twin Res Human Gene. 2008;11(6):579–85 Epub 2008/11/20.

    Google Scholar 

  42. 42.

    Pinto R, Rijsdijk F, Ronald A, Asherson P, Kuntsi J. The genetic overlap of attention-deficit/hyperactivity disorder and autistic-like traits: an investigation of individual symptom scales and cognitive markers. J Abnorm Child Psychol. 2016;44(2):335–45 Epub 2015/05/30.

    PubMed  Google Scholar 

  43. 43.

    Schork AJ, Won H, Appadurai V, Nudel R, Gandal M, Delaneau O, et al. A genome-wide association study of shared risk across psychiatric disorders implicates gene regulation during fetal neurodevelopment. Nat Neurosci. 2019;22(3):353–61 Epub 2019/01/30.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Lee SH, Ripke S, Neale BM, Faraone SV, Purcell SM, Perlis RH, et al. Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs. Nat Genet. 2013;45(9):984–94 Epub 2013/08/13.

    CAS  PubMed  Google Scholar 

  45. 45.

    Yao J, Liu Y, Husain J, Lo R, Palaparti A, Henderson J, et al. Combinatorial expression patterns of individual TLE proteins during cell determination and differentiation suggest non-redundant functions for mammalian homologs of Drosophila Groucho. Dev Growth Differ. 1998;40(2):133–46 Epub 1998/05/08.

    CAS  PubMed  Google Scholar 

  46. 46.

    Shen Q, Wang Y, Dimos JT, Fasano CA, Phoenix TN, Lemischka IR, et al. The timing of cortical neurogenesis is encoded within lineages of individual progenitor cells. Nat Neurosci. 2006;9(6):743–51 Epub 2006/05/09.

    CAS  PubMed  Google Scholar 

  47. 47.

    Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, et al. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. 2019;35(22):4851–3 Epub 2019/06/25.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Zeller T, Wild P, Szymczak S, Rotival M, Schillert A, Castagne R, et al. Genetics and beyond–the transcriptome of human monocytes and disease susceptibility. PLoS ONE. 2010;5(5):e10693 Epub 2010/05/27..

    PubMed  PubMed Central  Google Scholar 

  49. 49.

    Nascimento GR, Pinto IP, de Melo AV, da Cruz DM, Ribeiro CL, da Silva CC, et al. Molecular characterization of koolen de vries syndrome in two girls with idiopathic intellectual disability from Central Brazil. Mol Syndromol. 2017;8(3):155–60 Epub 2017/06/08.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Mc Cormack A, Taylor J, Te Weehi L, Love DR, George AM. A case of 17q21.31 microduplication and 7q31.33 microdeletion, associated with developmental delay, microcephaly, and mild dysmorphic features. Case Rep Gene. 2014;2014:658570 Epub 2014/03/22.

    Google Scholar 

  51. 51.

    Zollino M, Orteschi D, Murdolo M, Lattante S, Battaglia D, Stefanini C, et al. Mutations in KANSL1 cause the 17q21.31 microdeletion syndrome phenotype. Nat Gene. 2012;44(6):636–8 Epub 2012/05/01.

    CAS  Google Scholar 

  52. 52.

    Perrino PA, Talbot L, Kirkland R, Hill A, Rendall AR, Mountford HS, et al. Multi-level evidence of an allelic hierarchy of USH2A variants in hearing, auditory processing and speech/language outcomes. Commun Biol. 2020;3(1):180 Epub 2020/04/22.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Hahn H, Wicking C, Zaphiropoulous PG, Gailani MR, Shanley S, Chidambaram A, et al. Mutations of the human homolog of Drosophila patched in the nevoid basal cell carcinoma syndrome. Cell. 1996;85(6):841–51 Epub 1996/06/14.

    CAS  PubMed  Google Scholar 

  54. 54.

    Wang JC, Boyar FZ. Chromosomal microarray analysis as the first-tier test for the identification of pathogenic copy number variants in chromosome 9 pericentric regions and its challenge. Mol Cytogene. 2016;9:64 Epub 2016/08/16.

    Google Scholar 

  55. 55.

    Boudry-Labis E, Demeer B, Le Caignec C, Isidor B, Mathieu-Dramard M, Plessis G, et al. A novel microdeletion syndrome at 9q21.13 characterised by mental retardation, speech delay, epilepsy and characteristic facial features. Eur J Med Gene. 2013;56(3):163–70 Epub 2013/01/03.

    Google Scholar 

  56. 56.

    Ersland KM, Christoforou A, Stansberg C, Espeseth T, Mattheisen M, Mattingsdal M, et al. Gene-based analysis of regionally enriched cortical genes in GWAS data sets of cognitive traits and psychiatric disorders. PLoS ONE. 2012;7(2):e31687 Epub 2012/03/03.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Catalogue of Imprinted Genes and Parent-of-origin Effects. igc.otago.ac.nz. Accessed 6 Jan 2020.

  58. 58.

    Geneimprint. geneimprint.com. Accessed 6 Jan 2020.

  59. 59.

    Shimojima K, Adachi M, Tanaka M, Tanaka Y, Kurosawa K, Yamamoto T. Clinical features of microdeletion 9q22.3 (pat). Clin Genet. 2009;75(4):384–93 Epub 2009/03/27.

    CAS  PubMed  Google Scholar 

  60. 60.

    Siggberg L, Peippo M, Sipponen M, Miikkulainen T, Shimojima K, Yamamoto T, et al. 9q22 Deletion–first familial case. Orphanet J Rare Dis. 2011;6:45 Epub 2011/06/23.

    PubMed  PubMed Central  Google Scholar 

  61. 61.

    Sahoo T, Wang JC, Elnaggar MM, Sanchez-Lara P, Ross LP, Mahon LW, et al. Concurrent triplication and uniparental isodisomy: evidence for microhomology-mediated break-induced replication model for genomic rearrangements. EJHG. 2015;23(1):61–6 Epub 2014/04/10.

    CAS  PubMed  Google Scholar 

  62. 62.

    Mott R, Yuan W, Kaisaki P, Gan X, Cleak J, Edwards A, et al. The architecture of parent-of-origin effects in mice. Cell. 2014;156(1–2):332–42 Epub 2014/01/21.

    CAS  PubMed  PubMed Central  Google Scholar 

  63. 63.

    Võsa U, Claringbould A, Westra H-J, Bonder MJ, Deelen P, Zeng B, et al. Unraveling the polygenic architecture of complex traits using blood eQTL metaanalysis. bioRxiv. 2018:447367.

  64. 64.

    Webster E, Cho MT, Alexander N, Desai S, Naidu S, Bekheirnia MR, et al. De novo PHIP-predicted deleterious variants are associated with developmental delay, intellectual disability, obesity, and dysmorphic features. Cold Spring Harbor Mol Case Stud. 2016;2(6):a001172 Epub 2016/12/03.

    Google Scholar 

  65. 65.

    Jansen S, Hoischen A, Coe BP, Carvill GL, Van Esch H, Bosch DGM, et al. A genotype-first approach identifies an intellectual disability-overweight syndrome caused by PHIP haploinsufficiency. EJHG. 2018;26(1):54–63 Epub 2017/12/07.

    CAS  PubMed  Google Scholar 

  66. 66.

    Anney R, Klei L, Pinto D, Regan R, Conroy J, Magalhaes TR, et al. A genome-wide scan for common alleles affecting risk for autism. Hum Mol Genet. 2010;19(20):4072–82 Epub 2010/07/29.

    CAS  PubMed  PubMed Central  Google Scholar 

  67. 67.

    Vernes SC, Newbury DF, Abrahams BS, Winchester L, Nicod J, Groszer M, et al. A functional genetic link between distinct developmental language disorders. New Engl J Med. 2008;359(22):2337–45 Epub 2008/11/07.

    CAS  PubMed  Google Scholar 

  68. 68.

    Gosso MF, de Geus EJ, van Belzen MJ, Polderman TJ, Heutink P, Boomsma DI, et al. The SNAP-25 gene is associated with cognitive ability: evidence from a family-based study in two independent Dutch cohorts. Mol Psychiatry. 2006;11(9):878–86 Epub 2006/06/28.

    CAS  PubMed  Google Scholar 

  69. 69.

    Scerri TS, Brandler WM, Paracchini S, Morris AP, Ring SM, Richardson AJ, et al. PCSK6 is associated with handedness in individuals with dyslexia. Hum Mol Genet. 2011;20(3):608–14 Epub 2010/11/06.

    CAS  PubMed  Google Scholar 

  70. 70.

    Ioannidis JP. Why most discovered true associations are inflated. Epidemiology. 2008;19(5):640–8 Epub 2008/07/18.

    PubMed  Google Scholar 

  71. 71.

    Xiao R, Boehnke M. Quantifying and correcting for the winner’s curse in genetic association studies. Genet Epidemiol. 2009;33(5):453–62 Epub 2009/01/14.

    PubMed  PubMed Central  Google Scholar 

  72. 72.

    Thorup AA, Jepsen JR, Ellersgaard DV, Burton BK, Christiani CJ, Hemager N, et al. The danish high risk and resilience study–VIA 7–a cohort study of 520 7-year-old children born of parents diagnosed with either schizophrenia, bipolar disorder or neither of these two mental disorders. BMC Psychiatry. 2015;15:233 Epub 2015/10/04.

    PubMed  PubMed Central  Google Scholar 

  73. 73.

    Ellersgaard DV. Psychopathology, psychotic-like experiences, quality of life, and self-perception in seven-year-old children with familial high risk of schizophrenia or bipolar disorder 2018.

  74. 74.

    Bishop DVM. Test for reception of grammar: TROG-2. London: Pearson Assessment; 2003.

    Google Scholar 

  75. 75.

    Kaufman J, Birmaher B, Brent D, Rao U, Flynn C, Moreci P, et al. Schedule for affective disorders and schizophrenia for school-age children-present and lifetime version (K-SADS-PL): initial reliability and validity data. J Am Acad Child Adolesc Psychiatry. 1997;36(7):980–8 Epub 1997/07/01.

    CAS  PubMed  Google Scholar 

  76. 76.

    Ellersgaard D, Jessica Plessen K, Richardt Jepsen J, Soeborg Spang K, Hemager N, Klee Burton B, et al. Psychopathology in 7-year-old children with familial high risk of developing schizophrenia spectrum psychosis or bipolar disorder—The Danish High Risk and Resilience Study—VIA 7, a population-based cohort study. World Psychiatry. 2018;17(2):210–9 Epub 2018/06/02.

    PubMed  PubMed Central  Google Scholar 

  77. 77.

    Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.

    CAS  PubMed  PubMed Central  Google Scholar 

  78. 78.

    Anderson CA, Pettersson FH, Clarke GM, Cardon LR, Morris AP, Zondervan KT. Data quality control in genetic case-control association studies. Nat Protoc. 2010;5(9):1564–73 Epub 2010/11/19.

    CAS  PubMed  PubMed Central  Google Scholar 

  79. 79.

    Abecasis GR, Cardon LR, Cookson WO. A general test of association for quantitative traits in nuclear families. Am J Hum Genet. 2000;66(1):279–92 Epub 2000/01/13.

    CAS  PubMed  Google Scholar 

  80. 80.

    Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Merlin–rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002;30(1):97–101 Epub 2001/12/04.

    CAS  PubMed  Google Scholar 

  81. 81.

    Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26(18):2336–7 Epub 2010/07/17.

    CAS  PubMed  PubMed Central  Google Scholar 

  82. 82.

    Lange C, DeMeo DL, Laird NM. Power and design considerations for a general class of family-based association tests: quantitative traits. Am J Hum Genet. 2002;71(6):1330–41 Epub 2002/11/28.

    CAS  PubMed  PubMed Central  Google Scholar 

  83. 83.

    Weinberg CR, Wilcox AJ, Lie RT. A log-linear approach to case-parent-triad data: assessing effects of disease genes that act either directly or through maternal effects and that may be subject to parental imprinting. Am J Hum Genet. 1998;62(4):969–78 Epub 1998/06/13.

    CAS  PubMed  PubMed Central  Google Scholar 

  84. 84.

    van der Sluis S, Posthuma D. Single-locus association models. In: Neale BM, Ferreira MAR, Medland SE, Posthuma D, editors. Statistical Genetics: Gene Mapping through Linkage and Association. Abingdon: Taylor & Francis; 2008.

    Google Scholar 

  85. 85.

    Yeo I-K, Johnson RA. A new family of power transformations to improve normality or symmetry. Biometrika. 2000;87(4):954–9.

    Google Scholar 

  86. 86.

    Peterson RA, Cavanaugh JE. Ordered quantile normalization: a semiparametric transformation built for the cross-validation era. J Appl Stat. 2019:1–16.

  87. 87.

    Yee TW. The VGAM package for categorical data analysis. J Stat Softw. 2010;2010(10(2010)):1–34.

    Google Scholar 

  88. 88.

    R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2014.

  89. 89.

    Saint Pierre A, Vitezica Z, Martinez M. A comparative study of three methods for detecting association of quantitative traits in samples of related subjects. BMC Proceed. 2009;3(Suppl 7):122 Epub 2009/12/19.

    Google Scholar 

  90. 90.

    Iles MM. Linkage and Association: The Transmission/Disequilibrium Test for QTLs. Quantitative Trait Loci: Methods and Protocols. Totowa: Humana Press; 2002.

    Google Scholar 

  91. 91.

    QTDT—Online Reference. https://csg.sph.umich.edu/abecasis/QTDT/docs/index.htmlAccessed 30 Mar 2020.

  92. 92.

    Diao G, Lin DY. Improving the power of association tests for quantitative traits in family studies. Genet Epidemiol. 2006;30(4):301–13 Epub 2006/04/12.

    CAS  PubMed  Google Scholar 

  93. 93.

    Beasley TM, Erickson S, Allison DB. Rank-based inverse normal transformations are increasingly used, but are they merited? Behav Genet. 2009;39(5):580–95 Epub 2009/06/16.

    PubMed  PubMed Central  Google Scholar 

  94. 94.

    Lo S, Andrews S. To transform or not to transform: using generalized linear mixed models to analyse reaction time data. Front Psychol. 2015;6:1171 Epub 2015/08/25.

    PubMed  PubMed Central  Google Scholar 

  95. 95.

    Hoggart CJ, Venturini G, Mangino M, Gomez F, Ascari G, Zhao JH, et al. Novel approach identifies SNPs in SLC2A10 and KCNK9 with evidence for parent-of-origin effect on body mass index. PLoS Genet. 2014;10(7):e1004508 Epub 2014/08/01.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to express their gratitude to the dedicated families participating in the study. Figure 1 was exported with Daniel’s XL Toolbox; Manhattan and QQ plots were generated with an R script by Stephen Turner and Daniel Capurso; P-values for these plots were extracted from the QTDT output files using a Perl script by Thomas Scerri. We thank Mette Falkenberg Krantz for her assistance with the description of the VIA 7 cohort.

Funding

The VIA 7 project is supported by the Mental Health Services of the Capital Region of Denmark (Region Hovedstadens Psykiatri), the Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH), Aarhus University, the Tryg Foundation and the Beatrice Surovell Haskell Fund for Child Mental Health Research of Copenhagen. RN is supported by a postdoctoral grant from Region Hovedstadens Psykiatri. The funding bodies had no influence on the study design, analysis, or interpretation or on the writing of the manuscript.

Author information

Affiliations

Authors

Contributions

RN conceived the study, performed the QC of the genetic data, performed the genetic and statistical analyses, analyzed the results, wrote the paper; MJU performed the transformation of language test scores to norm-based scores; JO performed data management for VIA 7 and assisted with the QC of the pedigree information; CAJC, NH, DE, KSS, BKB, ANG, DLG contributed to the VIA 7 data collection and/or pilot study; J-BG oversaw sample preparation and genotyping and performed initial QC on the raw genetic data; AAET, JRMJ, OM, MN contributed to the conception of the VIA 7 study and its design, coordination and funding applications; TW designed and oversaw the genetic part of the VIA 7 study. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Thomas Werge or Merete Nordentoft.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Danish Data Protection Agency and follows all laws concerning the processing of personal data. Permission to draw data from registers was granted by the Danish Ministry of Health. The study protocol was sent to the Danish Committee on Health Research Ethics, who decided that ethical approval was not needed due to the observational nature of the study. The genetic part of the study obtained ethical approval from the outset of the study and The Danish High Risk and Resilience Study –VIA 7 was later incorporated into the protocol (Arv og Miljø–genetics and environment) as an appendix, which has then been approved by the ethics committee (ARV OG MILJØ: betydning for psykisk sygdom hos børn og unge (H-B-2009-026)). Written informed consent was obtained from all adult participants and from the legal guardians of participating children.

Consent for publication

Not applicable.

Competing interests

The authors have no competing interests to declare, but TW states that he has acted as a lecturer and scientific counselor to H. Lundbeck A/S. DE has been employed by H. Lundbeck A/S since March 3rd 2020.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1: Figure S1.

Manhattan and QQ plots for the discovery analyses.

Additional file 2: Table S1.

Suggestive associations (P ≤ 0.00001) from the discovery analyses.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Nudel, R., Christiani, C.A.J., Ohland, J. et al. Quantitative genome-wide association analyses of receptive language in the Danish High Risk and Resilience Study. BMC Neurosci 21, 30 (2020). https://doi.org/10.1186/s12868-020-00581-5

Download citation

Keywords

  • Linguistics
  • Genetics
  • GWAS
  • Parent-of-origin
  • Receptive language
  • Danish High Risk and Resilience Study