Next Article in Journal
SAGA-Dependent Histone H2Bub1 Deubiquitination Is Essential for Cellular Ubiquitin Balance during Embryonic Development
Next Article in Special Issue
New Frontiers in Neurodegeneration and Regeneration Associated with Brain-Derived Neurotrophic Factor and the rs6265 Single Nucleotide Polymorphism
Previous Article in Journal
Investigation of the Molecular Mechanisms Underlying the Antiatherogenic Actions of Kaempferol in Human THP-1 Macrophages
Previous Article in Special Issue
Untangling the Conformational Plasticity of V66M Human proBDNF Polymorphism as a Modifier of Psychiatric Disorder Susceptibility
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

NRN1 Gene as a Potential Marker of Early-Onset Schizophrenia: Evidence from Genetic and Neuroimaging Approaches

by
Carmen Almodóvar-Payá
1,2,†,
Maria Guardiola-Ripoll
1,2,†,
Maria Giralt-López
3,4,
Carme Gallego
5,
Pilar Salgado-Pineda
1,2,
Salvador Miret
2,6,7,
Raymond Salvador
1,2,
María J. Muñoz
8,
Luisa Lázaro
2,9,10,11,
Amalia Guerrero-Pedraza
1,8,
Mara Parellada
2,12,13,14,
María I. Carrión
15,
Manuel J. Cuesta
16,17,
Teresa Maristany
18,
Salvador Sarró
1,2,
Lourdes Fañanás
2,19,20,
Luis F. Callado
2,21,22,
Bárbara Arias
2,19,20,
Edith Pomarol-Clotet
1,2,* and
Mar Fatjó-Vilas
1,2,19,*
1
FIDMAG Germanes Hospitalàries Research Foundation, 08830 Sant Boi de Llobregat, Barcelona, Spain
2
Instituto de Salud Carlos III, Biomedical Research Network in Mental Health (CIBERSAM), 28029 Madrid, Madrid, Spain
3
Departament de Psiquiatria, Hospital Universitari Germans Trias i Pujol (HUGTP), 08916 Badalona, Barcelona, Spain
4
Departament de Psiquiatria i Medicina Legal, Universitat Autònoma de Barcelona (UAB), 08193 Bellaterra, Barcelona, Spain
5
Department of Cell Biology, Molecular Biology Institute of Barcelona (IBMB-CSIC), 08028 Barcelona, Barcelona, Spain
6
Centre de Salut Mental d’Adults de Lleida, Servei de Psiquiatria, Salut Mental i Addiccions, Hospital Universitari Santa Maria de Lleida, 25198 Lleida, Lleida, Spain
7
Institut de Recerca Biomèdica (IRB), 25198 Lleida, Lleida, Spain
8
Complex Assistencial en Salut Mental Benito Menni, 08830 Sant Boi de Llobregat, Barcelona, Spain
9
Department of Child and Adolescent Psychiatry and Psychology, Institute of Neurosciences, Hospital Clinic de Barcelona, 08036 Barcelona, Barcelona, Spain
10
Departament de Medicina, Universitat de Barcelona (UB), 08036 Barcelona, Barcelona, Spain
11
Institut d’Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), 08036 Barcelona, Barcelona, Spain
12
Servicio de Psiquiatría del Niño y del Adolescente, Hospital General Universitario Gregorio Marañón, 28007 Madrid, Madrid, Spain
13
Instituto de Investigación Sanitaria del Hospital Gregorio Marañón (IiSGM), 28007 Madrid, Madrid, Spain
14
Departamento de Psiquiatría, Facultad de Medicina, Universidad Complutense, 28040 Madrid, Madrid, Spain
15
Hospital Sant Rafael, 08035 Barcelona, Barcelona, Spain
16
Servicio de Psiquiatría, Hospital Universitario de Navarra, 31008 Pamplona, Navarra, Spain
17
Instituto de Investigación Sanitaria de Navarra (IdiSNA), 31008 Pamplona, Navarra, Spain
18
Departament de Diagnòstic per la Imatge, Hospital Sant Joan de Déu Fundació de Recerca, 08950 Esplugues de Llobregat, Barcelona, Spain
19
Departament de Biologia Evolutiva, Ecología i Ciències Ambientals, Universitat de Barcelona (UB), 08028 Barcelona, Barcelona, Spain
20
Institut de Biomedicina de la Universitat de Barcelona (IBUB), 08028 Barcelona, Barcelona, Spain
21
Department of Pharmacology, University of the Basque Country, UPV/EHU, 48940 Leioa, Bizkaia, Spain
22
Biocruces Bizkaia Health Research Institute, 48903 Barakaldo, Bizkaia, Spain
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(13), 7456; https://doi.org/10.3390/ijms23137456
Submission received: 29 April 2022 / Revised: 27 June 2022 / Accepted: 28 June 2022 / Published: 5 July 2022
(This article belongs to the Special Issue Neurotrophins: Roles and Function in Human Diseases)

Abstract

:
Included in the neurotrophins family, the Neuritin 1 gene (NRN1) has emerged as an attractive candidate gene for schizophrenia (SZ) since it has been associated with the risk for the disorder and general cognitive performance. In this work, we aimed to further investigate the association of NRN1 with SZ by exploring its role on age at onset and its brain activity correlates. First, we developed two genetic association analyses using a family-based sample (80 early-onset (EO) trios (offspring onset ≤ 18 years) and 71 adult-onset (AO) trios) and an independent case–control sample (120 healthy subjects (HS), 87 EO and 138 AO patients). Second, we explored the effect of NRN1 on brain activity during a working memory task (N-back task; 39 HS, 39 EO and 39 AO; matched by age, sex and estimated IQ). Different haplotypes encompassing the same three Single Nucleotide Polymorphisms(SNPs, rs3763180–rs10484320–rs4960155) were associated with EO in the two samples (GCT, TCC and GTT). Besides, the GTT haplotype was associated with worse N-back task performance in EO and was linked to an inefficient dorsolateral prefrontal cortex activity in subjects with EO compared to HS. Our results show convergent evidence on the NRN1 association with EO both from genetic and neuroimaging approaches, highlighting the role of neurotrophins in the pathophysiology of SZ.

1. Introduction

Substantial evidence highlights the importance of the genetic component in the aetiology of schizophrenia (SZ), with an estimated heritability of around 65–79% [1,2]. Indeed, genome-wide association studies (GWAS) have confirmed SZ’s polygenic architecture, resulting from the aggregated effect of low impact variants and reporting an SNP-based heritability of 24% [3]. Moreover, genomic data converge into identifiable biological pathways involved in neurodevelopment, particularly highlighting the mechanism of synaptic plasticity [3,4,5].
However, in the search for specific genetic factors related to SZ, studies face several challenges that arise from the genetic and phenotypic complexity of the disorder [6]. Then, it has been suggested that combining complementary designs, such as family-based and case–control, would clear the way to dissect the genetic influences of the disorder [7]. In this sense, family-based genetic association designs have the advantage of reducing the problem of stratification and spurious association when compared to case–control studies, while the latter usually allow for larger sample sizes [8].
In addition to the design, the approaches to the phenotypic complexity of SZ have also been considered using narrower phenotypes with particular aetiological significance to reduce the heterogeneity and identify specific genetic factors associated with the disorder. One of these phenotypes is the age at onset, which shows a heritability of around 33% [9]. Notably, early age at onset (EO) has captured much attention because it is considered a marker of a higher genetic liability than adult-onset (AO) [10,11]. The EO term includes cases with onset up to 18 years of age and, despite being arbitrary, roughly corresponds with the upper age cut-off in most published studies of child and adolescent psychosis [12]. In this support, EO subjects show a higher familial aggregation of SZ and other mental disorders [13], poorer premorbid adjustment [14] and neurocognitive performance [15], more severe outcomes [16,17] and more prominent alterations in neurodevelopmental trajectories than AO forms [18]. The few GWAS focused on searching for genetic loci associated with age at onset in SZ have confirmed that some variants overlap with those conferring risk for SZ, while others are pure modifiers [19,20,21,22]. Remarkably, EO patients present higher SZ polygenic risk scores than their siblings, with the scores effectively predicting an earlier age at onset [23]. Interestingly, the variants associated with an earlier age at onset converge into molecular networks related to nervous system development, the regulation of axon extension, modulation of glial proliferation, molecular transport, and cell-to-cell signalling and interactions [20,22].
Among genes with pivotal roles through all stages of the brain’s formation, there is the Neuritin 1 gene (NRN1, 6p25.1) (see review [24]), which is highly expressed in the hippocampus, the cerebral cortex and the cerebellum [25,26] in an activity-dependent manner [27,28]. Although the Nrn1 receptor and its downstream signalling effectors are still being studied, it seems that Nrn1 regulates synaptic excitability through the activation of the insulin receptor (IR) and its downstream signalling pathways [29,30]. Consequently, inadequate Nrn1 sustenance could translate into the abnormal formation of synapses, a reduced capacity to perform adaptive responses and, in turn, a higher risk of developing a mental disorder. In fact, the interest in the role of Nrn1 in SZ has been motivated by several studies, which have evidenced its impact on cognitive function through synaptic plasticity mechanisms. From cell- and animal-based approaches, it has been shown that the viral-mediated overexpression of NRN1 in different models (unpredictable stress-induced rat depression model, mice exposed to low-frequency electromagnetic fields and an Alzheimer’s disease model Tg2576 mouse) prevents the atrophy of dendrites and spines and improves associated behaviours, such as anxiety, depression, deficits in novel object recognition, learning and memory [31,32,33]. Additionally, the expression of NRN1 has been shown to increase in the hippocampus of mice exposed to electroconvulsive therapy and fluoxetine administration [32,34]. These studies highlight the potential therapeutic use of NRN1 in disorders associated with loss of cognitive function, such as SZ, and appeal for a better understanding of its molecular mechanisms. From human-based studies, NRN1 has been already defined as a candidate gene for SZ since specific allelic variants have been associated with an incremented risk of developing the disorder. Moreover, NRN1 has also been described as a modifier of the SZ phenotype due to its association with patients’ general cognitive ability and age at onset [35,36]. This suggests that NRN1 may be involved in critical mechanisms of brain development, particularly in those most susceptible to the earlier onset of the symptoms.
Neuroimaging data can provide evidence on how the genetic actors underlying an earlier age at onset contribute to the neurobiology of the disorder [37]. In this sense, functional neuroimaging studies focused on exploring the brain activity during working memory (WM) tasks (related to the capacity to retain and use mental items during a short period) are of particular interest. Subtle WM deviances have been described in the healthy siblings of subjects with SZ compared to healthy subjects (HS) in studies focused on cognition [38], brain activity [39] and connectivity [40]. This suggests that WM alterations in SZ are genetically influenced. Indeed, disabilities in this cognitive domain are considered to be a core feature of SZ [41] and have been reported to be even more severe in EO patients [15].
Several studies based on functional magnetic resonance neuroimaging (fMRI) and exploring brain networks supporting WM have consistently described frontoparietal differences in individuals with SZ when compared to HS. Most of these studies described the decreased activity of the dorsolateral prefrontal cortex (DLPFC), the ventrolateral prefrontal cortex (VLPFC) and anterior cingulate cortex (ACC) as a key mechanism of WM dysfunction [42]. The few functional neuroimaging studies specifically focused on individuals with EO have reported similar patterns of abnormal activations in these regions of the prefrontal cortex (e.g., VLPFC, DLPFC, and ACC) plus some limbic and temporal regions [43,44,45,46,47,48]. However, those studies are scarce, in part, due to the low rate of EO, which represents only about 8% among individuals with SZ [49], and they have reported inconsistent findings regarding the direction of the results. In this context, the study of the genetic WM correlates in individuals with EO forms is particularly pertinent since it could offer insights into the impact of genetic architecture on brain activity and, ultimately, on the clinical manifestation of SZ.
Considering all the above-cited evidence, we hypothesised that the polymorphic variability of NRN1 would be differentially associated with the risk of developing EO forms of SZ compared to AO. We developed this study by combining different designs (family-based and case–control sample approaches) to provide robustness to our findings. Additionally, we hypothesised that those genetic variants conferring risk for EO would differentially impact WM-related brain activity.

2. Results

2.1. Genetic Association Analyses

2.1.1. Family-Based

The genotypes/alleles counts and frequencies of EO/AO offspring and parents are listed in Supplementary Table S1. As shown in Table 1, within EO families, the GCT haplotype including SNP6, SNP7 and SNP8 (HAP678) was significantly under-transmitted from parents to affected offspring (pperm = 0.03). Our analyses did not reveal any association between the genetic variability at NRN1 with the risk for AO SSD, neither in the allelic, genotypic or haplotype approach.

2.1.2. Case–Control

The distribution of the genotypes/alleles in HS, EO and AO subjects is reported in Supplementary Table S2. As exposed in Table 1, we observed a significant association of SNP6 G allele (pperm = 0.02), SNP7 T allele (pperm = 0.03) and SNP8 T allele (pperm = 0.02) with EO SSD under an additive model. We also identified an association of two haplotypes including SNP6, SNP7 and SNP8 (HAP678) with the risk for EO SSD, which was in line with the SNP-based results. The GTT haplotype was significantly more frequent in subjects with EO SSD than in HS (pperm = 0.02), while the TCC was more frequent in HS (pperm = 0.01). Other 2-SNP and 4-SNP haplotypes containing these same variants were also associated with the risk for EO (Supplementary Table S3). Our analyses did not reveal any effect of genetic variability at NRN1 on the risk for AO SSD, in any of the tested models (allelic, genotypic and haplotypic).

2.2. Neuroimaging Genetic Association Analyses

2.2.1. N-Back Functional Response

The three groups (HS, EO and AO) showed typical WM-related activation and deactivation patterns (Supplementary Figures S1–S3). In addition, both EO and AO exhibited a deactivation failure when compared to HS in overlapping regions involving bilateral structures, such as the frontal gyrus (superior, medial and inferior orbital part), the olfactory area, the rectus and the anterior cingulate and the paracingulate gyri, as well as right structures, such as the superior and middle temporal gyrus, the parahippocampal gyrus, the hippocampus, the amygdala, the fusiform gyrus and the caudate nucleus (Supplementary Figure S4).
We detected a significant diagnosis x HAP678 (GTT) interaction for the EO vs. HS comparison in the 2-back vs. 1-back contrast in one cluster located at the superior and middle frontal gyrus, regions of the DLPFC (316 voxels, peak activation at MNI coordinates [−34,42,42], Zmax = 4.54, p = 0.0025, Figure 1A). To further interpret this result, mean activity scores for the 1-back and 2-back contrasts were plotted. As shown in Figure 1B, HS exhibited a cluster mean activity of around zero for the two contrasts, irrespective of their haplotypic profile. Subjects with EO without the risk haplotype showed a pattern towards increased cluster activity from 1-back to 2-back contrasts, whereas those patients carrying the risk haplotype presented a pattern towards decreased cluster activity. We did not observe any significant interaction on brain activity when we compared HS and AO groups.

2.2.2. N-Back Behavioural Response

First, subjects with EO exhibited a globally poorer performance of the N-back task than HS in both difficulty levels (mean (SD) d’1: EO 3.07 (1.16) and HS 4.14 (0.68), F = 13.00, p = 0.001; mean (SD) d’2: EO 2.06 (0.90) and HS 3.41 (0.88), F = 27.52, p < 0.001). While both groups showed different scores at the two levels of the task, their degree of decrease in performance from the 1-back to 2-back was similar (F = 1.55, p = 0.22).
Second, AO and HS exhibited a similar performance in the low memory load condition, but their performance diverged in the high memory load condition (mean (SD) d’1: AO 3.81 (0.90) HS 4.14 (0.68), F = 0.33, p = 0.57; mean (SD) d’2: AO 2.48 (0.82) and HS 3.41 (0.88), F = 15.27, p < 0.001). Then, as the performance of the two groups was similar for the 1-back, the degree of change from the 1-back to -back was more pronounced in subjects with AO than HS (F = 8.11, p = 0.01).
Third, EO performance at the low memory load condition was modulated by NRN1 haplotypic variability. Subjects carrying the HAP678 GTT showed a poorer performance when compared to those without the risk haplotype (mean (SD) d’1: non-carriers 3.67 (0.22) and carriers 2.86 (0.19), F = 5.66, p = 0.02) (Figure 2). No effect of NRN1 haplotypic variability on task performance was detected in either HS or AO subjects.

3. Discussion

In this study, we combined genetic association and neuroimaging approaches to deepen into the role of NRN1 in the age at onset of SZ. Regarding the genetic association approach, our study adds to the only two previous studies on the association of the NRN1 gene with SZ and other disorders within the spectrum [35,36], but it is the first to be developed through family-based and case–control designs in two independent samples. Our results, derived from the two samples, suggest that the variability at NRN1 may explain a modest proportion of the risk of EO. Concerning the neuroimaging approach, our study represents the first to explore WM neural correlates of NRN1. Our findings indicate that NRN1 variants conferring risk for SZ also have an effect on the performance of the N-back task, specifically within EO subjects. Additionally, we report brain activity differences between EO subjects and HS located at the DLPFC conditional to the same genetic variants.
Our genetic association analyses identified different SNPs and haplotypes at NRN1 associated with EO in SSD. We detected a significant under-transmission of the HAP678 GCT from parents to affected offspring, specifically in EO families. In parallel, through a case–control approach, we identified the effect of SNP6, 7 and 8 on the risk for EO SSD and the association of a risk haplotype encompassing these same polymorphisms HAP678 GTT. On the contrary, we did not detect any association with AO SSD. Therefore, our data converge into the view of polymorphisms at NRN1 (SNP6, SNP7, SNP8) as a relevant genetic variability source in modifying the neurodevelopment processes related to the earlier emergence of these disorders. It is of note that these SNPs have been previously associated with the risk for SZ in the two-preceding works [35,36]. Moreover, one of these studies also identified a role of NRN1 in age at onset of SSD [36]. However, our results should be interpreted in the context of the polygenic architecture of these disorders, as the effect of the SNPs and haplotypes is small (see the corresponding ORs). Still, this evidence suggests that those genes that influence brain development, such as NRN1, may modify illness traits, such as age at onset, and ultimately affect the risk for these disorders.
Due to the few studies focused on examining the association of the NRN1 gene with SZ, data from whole-genome approaches must be taken into consideration for the further interpretation of our results. First, different genetic linkage studies mapping SZ to a genomic location pointed towards the association of chromosome region 6p24-25 and highlighted NRN1 as a positional candidate gene [50,51,52,53]. However, as far as we know, NRN1 has not appeared as a significant locus in the latest genome-wide association studies [3]. These negative results could be explained due to the modifier properties of NRN1, which means that, as our results suggest, NRN1 modulates SSD phenotype through its impact on age at onset. Some linkage and whole-genome studies that specifically aimed to identify modifier loci related to the age at onset in SZ have highlighted the chromosome region 6p24 and some NRN1 neighbouring intergenic variants with putative regulatory roles on its expression [21,54]. Additionally, whole-genome approaches have also linked NRN1 and SZ through epigenetic mechanisms. In this respect, Pidsley et al., 2014 [55] identified, through a methylomic approach in human post-mortem prefrontal cortex samples, a wide genetic region that is hypomethylated in patients compared to controls, spanning the body of NRN1. This result suggests that NRN1 could be differentially expressed in the prefrontal cortex of subjects with SZ, a brain region repetitively described to be altered in this disorder [56].
To explore the neurobiological translation of the observed genetic variants conferring a higher risk for EO SZ, we developed a neuroimaging genetic study in a matched case–control sub-set. Our functional data suggest that the risk haplotype (HAP678 GTT) that is associated with the earlier emergence of the disorder is also associated with DLPFC activity changes within this group of patients. Concretely, through the analysis of differences between the two levels of the N-back task (2-back vs. 1-back contrast), we observed that EO subjects not carrying the risk haplotype changed DLPFC activity towards activation in response to the task’s increasing difficulty. At the same time, those carrying the risk haplotype were prone to decreased activity. The previous few studies exploring whole-brain activity differences between subjects with EO and HS have reported inconsistent findings regarding the implicated regions. Some reported reduced activation of the left VLPFC and extrastriate visual cortex [43], while others described VLPFC hyperactivation [47]. Other works reported the reduced engagement of the DLPFC, the ACC, frontal operculum and inferior and posterior parietal and caudate [44,45,46], contrary to other investigations that suggested increased activations in the ACC, medial temporal lobe structures, the insula and bilateral lateral temporal lobes [48]. In this respect, our results shed light on those controversial findings, as they provide evidence that genetic factors, in this case, the NRN1 gene, could be underlying these differences in DLPFC activity.
To interpret our functional results in the DLPFC, it is important to integrate brain activity findings with N-back behavioural data. On the one hand, the sustained activation of the prefrontal circuits is considered a key mechanism for executing high-memory-load tasks [57]. On the other hand, it is also known that the degree of change in DLPFC activity is related to the cognitive effort needed to perform the task. In other words, if the computational cost is unlikely to result in the accurate performance of the task, the prefrontal resources get disengaged [58]. In this sense, EO subjects not carrying the risk haplotype displayed a better performance at low memory load and exhibited a higher degree of DLPFC modulation towards activation in response to a higher memory load level. This result suggests that EO subjects without the risk haplotype may use greater prefrontal resources in response to task difficulty increase than those with the risk haplotype, who seem to reach activation and performance peaks at a lower processing load.
On the whole, our functional and behavioural results align with the preceding evidence linking NRN1 and cognitive performance in SZ [35,36] and executive function in HS [59]. Furthermore, they suggest that these prefrontal networks of sustained activation during WM, in which NRN1 appears to have a relevant role, might be especially sensitive to the earlier onset of psychosis. The specific effects in EO forms of SZ seem reasonable since several investigations have described adolescence as a crucial period for the development of the prefrontal cortex [60] and the reorganisation of the WM network [61,62]. This is also supported by several studies showing the greater recruitment of WM regions in adults than in children [63]. In this view, the earlier onset of the disorder might strongly impact the neural trajectories associated with WM development, potentially leading to WM-characteristic impairments in EO compared to AO [64]. Interestingly, regarding the specific role of DLPFC activity and NRN1 as potential markers of EO, a recent study using an innovative transcriptomic approach defined two molecularly distinct subgroups of subjects with SZ [65]. The first presented a DLPFC transcriptome very similar to that of HS, while the second exhibited a strikingly different DLPFC transcriptome, with the NRN1 gene included among the differentially expressed genes. These data suggest that fundamental biologic differences exist between subjects diagnosed with SZ. Thus, our results on the modulation effect of NRN1 haplotypic variability on brain function and performance contribute to bridging the gap between the role of NRN1 in synaptic plasticity processes and the pathophysiological mechanisms underlying SZ.
Towards a further understanding of such mechanisms, considering the putative effects of the analysed polymorphic sites on gene expression regulatory mechanisms represents a valuable resource to provide additional meaning and importance to our association data. Among the three variants encompassed by the HAP678 risk haplotype (LD = 0.99 in both samples), data from the RegulomeDB and Haploreg highlight the functional effects of SNP6 (rs3763180). There is evidence that this variant could modify the histone enhancer and promoter marks in the brain, contributing to the chromatin state at this locus. Moreover, this variant is predicted to alter motifs that overlap the recognition sequence of different transcription factors, such as the alpha isoform of the CCAAT-enhancer binding proteins (C/EBP), PBX homeobox 3 (PBX3) and the Neuron-Restrictive Silencing Factor (NRSF). Interestingly, the change of a T to a G in that position is linked to increased NRSF affinity, implicated in the programming of stress-sensitive neurons by neonatal experience through epigenetic mechanisms, promoting resilience to stress-related emotional disorders [66]. The other two variants included in the HAP678 presented a lower functionality score; still, data show their putative modulatory effects on the affinity of some transcription factors. For instance, the SNP7 (rs10484320) is suggested to modify the binding of the TATA-binding protein (TBP) and PU.1. The TBP has been associated with the risk for SZ, age at onset and prefrontal function [67]. Additionally, higher levels of the transcription factor PU.1, required for the development of the immune system, have been detected in post-mortem brain samples from individuals diagnosed with SZ compared to HS [68]. Additionally, accordingly to Brainiac data, when the effect of these three SNPs stratifies the expression of NRN1 transcripts, genotype-based differences emerge in the hippocampus, and a trend effect is detected in the cortex. These lines of evidence suggest putative molecular mechanisms by which the SNPs included in the HAP678 may affect the complex phenotype of SZ. Nevertheless, further functional data on these SNPs are needed to fully characterise their impact on the underlying mechanisms that connect NRN1 and age at onset of psychosis.
Finally, our study should be interpreted in the context of some limitations. First, regarding our genetic association approach, the samples could be considered to be relatively small. However, according to the statistical power of our analyses and after multiple testing correction procedures, we concurrently identified, in two independent samples, the impact of NRN1 genetic variants on the risk for the earlier onset of SZ. Second, all the variants included in the present study are polymorphic; however, it is known that a certain proportion of the variance in genetic liability of SZ is also accounted for by rare variants [69,70]. Therefore, different approaches analysing the combined role of common and low-frequency variants along NRN1 gene on SZ would be of potential interest. Third, while the present study has not directly analysed the functional consequences of the NRN1 variants associated with EO, our results and the available functional data suggest the need for cell-based studies integrating genetic variability information. Fourth, in the case of neuroimaging approaches, although we compared EO subjects and HS, patients were scanned in their adulthood, years after the onset of the illness. Therefore, illness duration and related clinical variables could have affected the results. Based on this, we checked the possible impact of illness duration or medication on the mean activity and the d’ scores through regressions. While we cannot completely rule out the effect of these variables the lack of significance suggests that our results are not modulated by them. Additionally, it should be underlined that activation differences at prefrontal regions have been observed in unaffected first relatives of SZ patients [39], individuals at clinical high risk for psychosis [71] and individuals with treatment-naïve first episode psychosis [72], suggesting that this pattern may represent an intrinsic feature of SZ rather than a medication effect. Lastly, the absence of representation of diverse ethnic groups and the low proportion of females within our EO group hampers the extrapolation of our results and demands the need for new studies in larger samples with equal representation of those populations.

4. Materials and Methods

4.1. Sample

This study included 798 individuals (Table 2). Two independent samples were used to develop separate genetic association analyses: (i) Sample 1 comprised 151 trios (with an offspring diagnosed with schizophrenia spectrum disorders (SSD) plus 302 healthy parents), (ii) Sample 2 consisted of 225 independent patients diagnosed with SSD and 120 HS. Also, from Sample 2, a sub-set of cases with SZ (39 EO and 39 AO) and HS (39) (matched by sex, age and estimated IQ) was selected to develop a neuroimaging genetic analysis (Sample 3). Participants were drawn from admissions to both Child and Adolescent and Adult Psychiatric Units. All HS were recruited from non-medical staff working in the hospital, their relatives and acquaintances, plus independent sources in the community.
All patients were evaluated by experienced psychiatrists and met the DSM-IV-TR criteria for SSD, including schizophrenia, schizophreniform disorder, schizoaffective disorder and psychosis disorder not otherwise specified (Table 2). Patients up to 17 years old were diagnosed following Kiddie Schedule for Affective Disorders and Schizophrenia (KSDAS, [73]), while the Comprehensive Assessment of Symptoms and History (CASH, [74]) or the Structured Clinical Interview for DSM Disorders (SCID, [75]) was used for adult patients. Age at onset of the first episode was determined using these clinical schedules and/or the Symptom Onset in Schizophrenia inventory (SOS, [76]). Following previous studies [12], subjects with SSD were classified as either EO when the first episode occurred before or at 18 years, or as adult-onset AO when presented at age 19 or older.
The general exclusion criteria included an age above 65 years, major medical illnesses that could affect brain functions, substance-induced psychotic disorder, neurological conditions and having had at least one parent not from European ancestry. Moreover, all the relatives and HS underwent a clinical interview on personal and/or familial psychiatric history using Family Interview for Genetic Studies (FIGS) [77] and those who reported a personal history of mental illness or treatment with psychotropic medication were excluded.
For the subjects included in the neuroimaging study, the exclusion criteria also included an estimated IQ under 70, left manual dominance, and a history of head trauma with loss of consciousness. The evaluation of patients comprised the Positive and Negative Symptoms Scale (PANSS), while the estimated IQ of both patients and controls, was assessed using the Word Accentuation Test [78], which requires the pronunciation of 30 low-frequency Spanish words whose accents were removed.
All participants provided written consent after being informed about the study procedures and implications. In the case of patients below the age of 18, written consent was also obtained from their parents. The study was performed following the guidelines of the institutions involved and was approved by the local ethics committee of the centre. All procedures were carried out according to the Declaration of Helsinki.

4.2. Genotyping

Genomic DNA was obtained for all individuals either from buccal mucosa through cotton swabs or from peripheral blood cells by punction and extracted using an ATP Genomic DNA Mini Kit Tissue (Teknokroma Analítica, S.A., Sant Cugat del Vallès, Barcelona, Spain) or using a Realpure SSS Kit for DNA Extraction (Durviz, S.L.U, Valencia, Spain), respectively.
All Single Nucleotide Polymorphisms (SNPs) were determined via a fluorescence-based allelic discrimination procedure (Applied Biosystems Taqman 5′-exonuclease assays) using standard conditions.
The information about the SNPs is given in Table 3. The SNPs were selected based on two previous studies [35,36]. All SNPs had a minor allele frequency above 5% and were non-coding. As previous evidence suggests that non-coding variants exert important regulatory effects [79] and that such effects are particularly important in SZ [80], the functional consequences of the analysed SNPs were evaluated using different resources. First, HaploReg was used to obtain information about the impact of non-coding variants on chromatin state, protein binding, sequence conservation across mammals, regulatory motifs and expression (https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php, [81], accessed on 1 June 2022). It showed that several SNPs (from rs12333117 to rs3763180) are classified as genetic promoters or enhancers in the brain tissue, based on histone marker data from the Epigenetic Roadmap. Moreover, all the variants, except for one (rs582186), are predicted to change the affinity of multiple regulatory motifs based on data from the ENCODE project. Second, the Regulome DataBase (http://www.regulomedb.org/, [82], accessed on 1 June 2022), a model integrating functional genomic features, was used to obtain a functional probability score for each SNP. This score ranges from 0 to 1, with 1 being most likely to be a regulatory variant [83]. As shown in Table 3, several of the selected SNPs had a score above 0.61 (from rs12333117 to rs3763180). Third, the Brain eQTL Almanac (Braineac), which is a web-based (http://www.braineac.org, accessed on 1 June 2022) resource to access the UK Brain Expression Consortium (UKBEC) dataset, showed that the assessed NRN1 transcripts are mainly expressed in the cortex, hippocampus and cerebellum. In addition, this tool was also used to evaluate the effect of the SNPs associated with EO forms of SZ on brain expression patterns (see Discussion).
All markers were in Hardy–Weinberg equilibrium in the two samples: the family-based (parents and offspring) and the case–control (subjects with SSD and HS). The total genotypic call rate was 95.48 and 98.85%, respectively.

4.3. fMRI Task Description and Acquisition Parameters

4.3.1. N-Back Task

Functional images were acquired while participants performed a sequential-letter version of the N-back task, which engages storage and executive processes related to attention and WM. Briefly, in this task, letters were presented sequentially in a random way, and the participants were required to press a button when the letter shown on the screen matched the one presented one step prior in the sequence (condition 1-back) or the one from two steps before in the sequence (condition 2-back). The two levels of memory load were presented in a block design manner. Each block consisted of 24 letters shown every 2 s (1 s on, 1 s off), and all blocks contained 5 letter repetitions located randomly within the blocks. Four 1-back and four 2-back blocks were presented in an interleaved way, and between them, a baseline stimulus (an asterisk flashing with the same frequency as the letters) was presented for 16 s. Letters were displayed in green for 1-back blocks and in red for 2-back blocks to identify which condition had to be performed. All participants went through a training session before entering the scanner.

4.3.2. N-Back Performance Data

To measure the behavioural performance of the task, we used the signal detection theory index sensitivity, d’ scores [84]. Higher values of d’ indicate a better ability to discriminate between targets and distractors, while negative values indicate that subjects are not performing the task. Then, all the individuals included in the analyses had positive d’ values (d’1 for 1-back and d’2 for 2-back).

4.3.3. fMRI Acquisition Parameters

The fMRI data acquisition was performed with a GE Sigma 1.5-T scanner (General Electric Medical Systems, Milwaukee, WI, USA) at Hospital Sant Joan de Déu (Barcelona, Spain). Functional images included 266 volumes for each individual and a gradient echo-planar imaging sequence depicting the blood oxygen level-dependent (BOLD) signal. Each volume contained 16 axial planes acquired with the following parameters: repetition time = 2000 ms, echo time = 20 ms, flip angle = 70°, section thickness = 7 mm, section skip = 0.7 mm, in-plane resolution = 3 × 3 mm. The first 10 volumes were discarded to avoid T1-saturation effects.

4.4. Statistical Analyses

4.4.1. Design

First, family-based genetic association analyses were conducted within EO families and AO families separately. Second, case–control genetic association analyses were tested by comparing HS to EO and AO patients independently. Third, neuroimaging analyses were developed to determine whether brain activity differences exist between HS and subjects with EO and AO SZ depending on the NRN1 genetic variability.

4.4.2. Genetic Association Analyses

Haploview v4.1 [85] was employed to estimate the linkage disequilibrium (LD) between NRN1 SNPs. Three haplotype blocks were identified (Block 1: SNP1-SNP3, Block 2: SNP4-SNP5 and Block 3: SNP6-SNP11) in both the family-based and the case–control sample (Supplementary Figure S5). Hardy–Weinberg and genetic association analyses between NRN1 SNPs/haplotypes and SZ risk/age at onset were conducted using PLINK-v1.07 software [86].
For the family-based analyses, SNP and haplotype associations were tested using the Transmission Disequilibrium Test (TDT). This test evaluates whether the transmission frequency of alleles/haplotypes from heterozygous parents to their affected children deviates from the expected Mendelian frequency by comparing the transmitted and not transmitted alleles/haplotypes.
The case–control analyses were conducted using logistic regressions under different models of inheritance (allelic, genotypic/additive, recessive and dominant), all adjusted by sex.
In both genetic association approaches, a cut-off threshold for rare haplotypes of 1% and a sliding window approach were applied to the haplotype analyses. The odds ratios (OR) were estimated either from the absolute number of alleles/haplotypes transmitted and not transmitted from parents to affected offspring or from the absolute number of alleles/haplotypes estimated in patients and controls. Multiple testing corrections (10,000 permutations procedure) were applied to all analyses, and all the reported p-values are those obtained with this correction (pperm). As haplotype TDT implemented in PLINK does not include the permutation procedure, to confirm our associations, all the possible haplotypes were reconstructed based on the most likely expectation maximisation (EM) phase and once reconstructed, haplotypic associations were tested using a simple TDT with 10,000 permutations.
The statistical power was calculated, in the case of the family-based sample, using the ‘trio’ R package version 3.1.2 available at http://www.bioconductor.org, accessed on 1 April 2022. As the sample consisted of 151 trios and the MAF of selected SNPs ranged from 0.20 to 0.50 (Table 1), by assuming an allelic TDT model and a statistical power of 0.80, the smaller detectable relative risk is 1.75. In the case of the case–control sample, the statistical power was calculated using the ‘genpwr’ R package version 1.0.2. As the sample comprised 345 subjects and the MAF of selected SNPs ranged from 0.16 to 0.49, assuming a logistic model and a power of 0.80, the smaller detectable odds ratio is 1.57.

4.4.3. Neuroimaging Association Study

Based on our genetic association results highlighting the genetic region spanning SNP 6, 7 and 8, we performed the neuroimaging analyses with the haplotype significantly associated with the risk for EO (the HAP678 GTT) in the matched sub-set (39 HS, 39 EO and 39 AO subjects). For these analyses, individuals’ possible haplotype phases were estimated using PLINK and only those with a probability ≥ 95% were included. Because of the haplotypic frequencies in our sample, the analyses were conducted considering the haplotype as a dichotomous variable, and each subject was classified as a non-carrier (0 copies of the risk haplotype) or carrier (1 or 2 copies of the risk haplotype).
Functional MRI pre-processing and analyses were performed with the FEAT tool, from FSL software (FMRIB Software Library, University of Oxford, Oxford, UK; [84]). Pre-processing included motion correction (using the MCFLIRT algorithm with 6 degrees of freedom) and co-registration and normalisation to a common stereotactic space (Montreal Neurological Institute [MNI] template with 2 × 2 × 2 mm resolution) using linear transformations with 12 degrees of freedom. Before group analyses, normalised images were spatially filtered with a Gaussian filter (FWHM = 5mm). To minimise unwanted movement-related effects, individuals with an estimated maximum absolute movement > 3.0 mm or an average absolute movement > 0.3 mm were excluded from analyses.
At the single-subject level analysis, General Linear Models (GLMs) were fitted to generate individual activation maps for each condition of interest compared to baseline and for the comparison between conditions (1-back vs. baseline, 2-back vs. baseline and 2-back vs. 1-back). Temporal derivatives for each condition of interest, as well as movement parameters (six in total, three rotations and three translations), were also included as additional regressors. Fixation periods were not modelled and thus acted as an implicit baseline (i.e., to compare a condition of interest of any given task with its baseline periods, the average BOLD signal from all the baseline periods across the whole task is subtracted from that of the blocks corresponding to the condition of interest). Images were high-pass filtered with a 130 s cut-off. All statistical tests were performed at the cluster level with a corrected p-value of 0.05 and an initial height threshold of 3.1 (equivalent to an uncorrected p-value of 0.001), using the Standard Field Theory correction implemented in FSL [87].
At the group-level analyses, we studied brain activations and deactivations associated with the execution of the N-back task within each group for all the contrasts, as well as the differences between groups using ANOVA models (two comparisons whole-brain corrected: HS vs. EO/HS vs. AO adjusted for age, sex, and premorbid-IQ).
Since the 1-back requires the maintenance of the target in the memory (keeping track of the target when the consecutive letter is represented) and the 2-back demands both maintenance and target switching (updating the target identity with the appearance of each new letter), we decided to focus our interaction analyses on the 2-back vs. 1-back contrast as it highlights regional responses specific to a higher WM capacity [88,89]. Accordingly, our group differences between HS and patients (both EO and AO), were more pronounced in the 2-back vs. 1-back contrast (Supplementary Figure S5).
The interaction effect between the diagnosis and the risk haplotype HAP678 GTT was investigated using a regression model, which tests whether the slope between groups and haplotype differs (two models, whole-brain corrected: HS and EO/HS and AO, adjusted for age, sex, and premorbid-IQ). Four contrasts were explored (EO > HS and the reverse contrast; AO > HS and the reverse contrast). Therefore, to control for all the comparisons, the significance threshold used was set to p < 0.05/4 = 0.0125. To interpret the direction of the interaction results, we estimated individual mean activity scores from the areas where a significant interaction was detected using the FSLSTATS tool in FSL and afterwards, these values were plotted using SPSS.
Analyses of the behavioural data (d’1 and d’2) were carried out using SPSS and ANOVA models adjusted for age, sex and estimated IQ. First, N-back task performance was compared between HS and patients (two comparisons: HS vs. EO/HS vs. AO). Second, we explored the effect of the risk haplotype HAP678 GTT on N-back task performance in each group (HS, EO and AO). Third, we tested the interaction between diagnosis and the risk haplotype HAP678 GTT (two models: HS vs. EO/HS vs. AO). Additionally, the change between d’1 and d’2 was explored using repeated measures ANOVA models adjusted for age, sex and estimated IQ. These repeated measures models were conducted, first, between groups (two comparisons: HS vs. EO/HS vs. AO) and, second, between groups depending on the haplotype (two models: HS vs. EO/HS vs. AO).

5. Conclusions

Our results contribute to the understanding of the molecular mechanisms underlying SZ early age at onset. Specifically, our work suggests that studying the role of neurotrophins (such as NRN1), in specific phenotypes with particularly etiological underpinnings (such as early-onset) and their effect on intermediate phenotypes (such as functional neuroimaging data), helps to elucidate the impact of common genetic variability on biological networks underlying mental disorders.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23137456/s1.

Author Contributions

Conceptualisation, C.A.-P., C.G. and M.F.-V.; Data curation, C.A.-P., M.G.-R. and P.S.-P.; Formal analysis, C.A.-P., M.G.-R. and M.F.-V.; Funding acquisition, C.G., L.F., L.F.C., E.P.-C. and M.F.-V.; Investigation, C.A.-P., M.G.-R., M.G.-L., B.A. and M.F.-V.; Methodology, C.A.-P., R.S., A.G.-P., E.P.-C. and M.F.-V.; Project administration, M.F.-V.; Resources, M.G.-L., P.S.-P., S.M., M.J.M., L.L., M.P., M.I.C., M.J.C., T.M., S.S., E.P.-C. and M.F.-V.; Supervision, M.F.-V.; Visualisation, C.A.-P. and M.G.-R.; Writing—original draft, C.A.-P., M.G.-R. and M.F.-V.; Writing—review & editing, C.A.-P., M.G.-R., M.G.-L., C.G., P.S.-P., S.M., R.S., M.J.M., L.L., A.G.-P., M.P., M.I.C., M.J.C., T.M., S.S., L.F., L.F.C., B.A., E.P.-C. and M.F.-V. All authors have read and agreed to the published version of the manuscript.

Funding

This study received funding provided by: (i) Fundación Alicia Koplowitz; (ii) Acadèmia de les Ciències Mèdiques i de la Salut de Catalunya i de Balears (predoctoral contract to C.A.-P.); (iii) the Instituto de Salud Carlos III through a PFIS predoctoral contract to M.G.-R. (FI19/0352) and a Miguel Servet contract to M.F.-V. (CP20/00072), co-funded by European Regional Development Fund (ERDF)/European Social Fund “Investing in your future”; (iv) the Comissionat per a Universitats i Recerca del DIUE of the Generalitat de Catalunya (Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR), 2017SGR1271 and 2017SGR1577).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and ethical approval was obtained from the research ethical committee “CEIC Hermanas Hospitalarias del Sagrado Corazón de Jesús” (protocol code PR-2016-07 and date of approval 31 March 2016).

Informed Consent Statement

Informed consent statement was obtained from all the participants involved in the study. Written informed consent has been obtained from all the participants to publish this paper.

Conflicts of Interest

The authors declare to have no competing interest to declare.

References

  1. Hilker, R.; Helenius, D.; Fagerlund, B.; Skytthe, A.; Christensen, K.; Werge, T.M.; Nordentoft, M.; Glenthøj, B. Heritability of schizophrenia and schizophrenia spectrum based on the nationwide Danish Twin register. Biol. Psychiatry 2018, 83, 492–498. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Sullivan, P.F.; Kendler, K.S.; Neale, M.C. Schizophrenia as a complex trait: Evidence from a meta-analysis of twin studies. Arch. Gen. Psychiatry 2003, 60, 1187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Trubetskoy, V.; Pardiñas, A.F.; Qi, T.; Panagiotaropoulou, G.; Awasthi, S.; Bigdeli, T.B.; Bryois, J.; Chen, C.-Y.; Dennison, C.A.; Hall, L.S.; et al. Mapping genomic Loci implicates genes and synaptic biology in schizophrenia. Nature 2022, 604, 502–508. [Google Scholar] [CrossRef] [PubMed]
  4. Pardiñas, A.F.; Holmans, P.; Pocklington, A.J.; Escott-Price, V.; Ripke, S.; Carrera, N.; Legge, S.E.; Bishop, S.; Cameron, D.; Hamshere, M.L.; et al. Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection. Nat. Genet. 2018, 50, 381–389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ripke, S.; Neale, B.M.; Corvin, A.; Walters, J.T.R. Biological insights from 108 schizophrenia-associated genetic Loci. Nature 2014, 511, 421–427. [Google Scholar] [CrossRef] [Green Version]
  6. Liang, S.G.; Greenwood, T.A. The impact of clinical heterogeneity in schizophrenia on genomic analyses. Schizophr. Res. 2015, 161, 490–495. [Google Scholar] [CrossRef] [Green Version]
  7. Glahn, D.C.; Nimgaonkar, V.L.; Raventós, H.; Contreras, J.; McIntosh, A.M.; Thomson, P.A.; Jablensky, A.; McCarthy, N.S.; Charlesworth, J.C.; Blackburn, N.B.; et al. Rediscovering the value of families for psychiatric genetics research. Mol. Psychiatry 2019, 24, 523–535. [Google Scholar] [CrossRef]
  8. Lange, E.M.; Sun, J.; Lange, L.A.; Zheng, S.L.; Duggan, D.; Carpten, J.D.; Gronberg, H.; Isaacs, W.B.; Xu, J.; Chang, B.-L. Family-based samples can play an important role in genetic association studies. Cancer Epidemiol. Prev. Biomark. 2008, 17, 2208–2214. [Google Scholar] [CrossRef] [Green Version]
  9. Hare, E.; Glahn, D.C.; Dassori, A.; Raventos, H.; Nicolini, H.; Ontiveros, A.; Medina, R.; Mendoza, R.; Jerez, A.; Muñoz, R.; et al. Heritability of age of onset of psychosis in schizophrenia. Am. J. Med. Genet. B Neuropsychiatr. Genet. 2010, 153, 298–302. [Google Scholar] [CrossRef]
  10. Esterberg, M.L.; Trotman, H.D.; Holtzman, C.; Compton, M.T.; Walker, E.F. The impact of a family history of psychosis on age-at-onset and positive and negative symptoms of schizophrenia: A meta-analysis. Schizophr. Res. 2010, 120, 121–130. [Google Scholar] [CrossRef]
  11. Hilker, R.; Helenius, D.; Fagerlund, B.; Skytthe, A.; Christensen, K.; Werge, T.M.; Nordentoft, M.; Glenthøj, B. Is an early age at illness onset in schizophrenia associated with increased genetic susceptibility? Analysis of data from the nationwide Danish Twin register. EBioMedicine 2017, 18, 320–326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Hollis, C. Adult outcomes of child- and adolescent-onset schizophrenia: Diagnostic stability and predictive validity. Am. J. Psychiatry 2000, 157, 1652–1659. [Google Scholar] [CrossRef] [PubMed]
  13. Forsyth, J.K.; Asarnow, R.F. Genetics of childhood-onset schizophrenia: 2019 update. Child Adolesc. Psychiatr. Clin. N. Am. 2020, 29, 157–170. [Google Scholar] [CrossRef] [PubMed]
  14. Larsen, T.K.; Friis, S.; Haahr, U.; Johannessen, J.O.; Melle, I.; Opjordsmoen, S.; Rund, B.R.; Simonsen, E.; Vaglum, P.; McGlashan, T.H. Premorbid adjustment in first-episode non-affective psychosis: Distinct patterns of pre-onset course. Br. J. Psychiatry 2004, 185, 108–115. [Google Scholar] [CrossRef] [PubMed]
  15. Rajji, T.K.; Ismail, Z.; Mulsant, B.H. Age at onset and cognition in schizophrenia: Meta-analysis. Br. J. Psychiatry 2009, 195, 286–293. [Google Scholar] [CrossRef] [Green Version]
  16. Clemmensen, L.; Vernal, D.L.; Steinhausen, H.-C. A systematic review of the long-term outcome of early onset schizophrenia. BMC Psychiatry 2012, 12, 150. [Google Scholar] [CrossRef] [Green Version]
  17. Rapoport, J.L.; Addington, A.M.; Frangou, S.; Psych, M.R.C. The neurodevelopmental model of schizophrenia: Update. Mol. Psychiatry 2005, 10, 434–449. [Google Scholar] [CrossRef]
  18. Rapoport, J.L.; Gogtay, N. Childhood onset schizophrenia: Support for a progressive neurodevelopmental disorder. Int. J. Dev. Neurosci. Off. J. Int. Soc. Dev. Neurosci. 2011, 29, 251–258. [Google Scholar] [CrossRef] [Green Version]
  19. Bergen, S.E.; O’Dushlaine, C.T.; Lee, P.H.; Fanous, A.H.; Ruderfer, D.M.; Ripke, S.; Sullivan, P.F.; Smoller, J.W.; Purcell, S.M.; Corvin, A. Genetic modifiers and subtypes in schizophrenia: Investigations of age at onset, severity, sex and family history. Schizophr. Res. 2014, 154, 48–53. [Google Scholar] [CrossRef] [Green Version]
  20. Guo, S.; Liu, J.; Li, W.; Yang, Y.; Lv, L.; Xiao, X.; Li, M.; Guan, F.; Luo, X.-J. Genome wide association study identifies four loci for early onset schizophrenia. Transl. Psychiatry 2021, 11, 248. [Google Scholar] [CrossRef]
  21. Wang, K.-S.; Liu, X.; Zhang, Q.; Aragam, N.; Pan, Y. Genome-wide association analysis of age at onset in schizophrenia in a European-American sample. Am. J. Med. Genet. B Neuropsychiatr. Genet. 2011, 156, 671–680. [Google Scholar] [CrossRef] [PubMed]
  22. Woolston, A.L.; Hsiao, P.-C.; Kuo, P.-H.; Wang, S.-H.; Lien, Y.-J.; Liu, C.-M.; Hwu, H.-G.; Lu, T.-P.; Chuang, E.Y.; Chang, L.-C.; et al. Genetic Loci associated with an earlier age at onset in multiplex schizophrenia. Sci. Rep. 2017, 7, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Ahn, K.; An, S.S.; Shugart, Y.Y.; Rapoport, J.L. Common polygenic variation and risk for childhood-onset schizophrenia. Mol. Psychiatry 2016, 21, 94–96. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhou, S.; Zhou, J. Neuritin, a Neurotrophic factor in nervous system physiology. Curr. Med. Chem. 2014, 21, 1212–1219. [Google Scholar] [CrossRef]
  25. Aid, T.; Kazantseva, A.; Piirsoo, M.; Palm, K.; Timmusk, T. Mouse and rat BDNF gene structure and expression revisited. J. Neurosci. Res. 2007, 85, 525–535. [Google Scholar] [CrossRef]
  26. Nedivi, E.; Fieldust, S.; Theill, L.E.; Hevron, D. A Set of genes expressed in response to light in the adult cerebral cortex and regulated during development. Proc. Natl. Acad. Sci. USA 1996, 93, 2048–2053. [Google Scholar] [CrossRef] [Green Version]
  27. Nedivi, E.; Hevroni, D.; Naot, D.; Israeli, D.; Citri, Y. Numerous candidate plasticity-related genes revealed by differential cDNA cloning. Nature 1993, 363, 718–722. [Google Scholar] [CrossRef]
  28. Tongiorgi, E. Activity-Dependent expression of brain-derived neurotrophic factor in dendrites: Facts and open questions. Neurosci. Res. 2008, 61, 335–346. [Google Scholar] [CrossRef]
  29. Lu, J.-M.; Liu, D.-D.; Li, Z.-Y.; Ling, C.; Mei, Y.-A. Neuritin enhances synaptic transmission in medial prefrontal cortex in mice by increasing CaV3.3 surface expression. Cereb. Cortex 2017, 27, 3842–3855. [Google Scholar] [CrossRef] [Green Version]
  30. Yao, J.-J.; Gao, X.-F.; Chow, C.-W.; Zhan, X.-Q.; Hu, C.-L.; Mei, Y.-A. Neuritin activates insulin receptor pathway to up-regulate Kv4.2-mediated transient outward K + current in rat cerebellar granule neurons. J. Biol. Chem. 2012, 287, 41534–41545. [Google Scholar] [CrossRef] [Green Version]
  31. An, K.; Jung, J.H.; Jeong, A.Y.; Kim, H.G.; Jung, S.Y.; Lee, K.; Kim, H.J.; Kim, S.-J.; Jeong, T.-Y.; Son, Y.; et al. Neuritin can normalize neural deficits of Alzheimer’s disease. Cell Death Dis. 2014, 5, e1523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Son, H.; Banasr, M.; Choi, M.; Chae, S.Y.; Licznerski, P.; Lee, B.; Voleti, B.; Li, N.; Lepack, A.; Fournier, N.M.; et al. Neuritin produces antidepressant actions and blocks the neuronal and behavioral deficits caused by chronic stress. Proc. Natl. Acad. Sci. USA 2012, 109, 11378–11383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Choi, Y.; Lee, K.; Ryu, J.; Kim, H.G.; Jeong, A.Y.; Woo, R.-S.; Lee, J.-H.; Hyun, J.W.; Hahn, S.; Kim, J.-H.; et al. Neuritin attenuates cognitive function impairments in Tg2576 mouse model of Alzheimer’s Disease. PLoS ONE 2014, 9, e104121. [Google Scholar] [CrossRef]
  34. Park, K.W.; Kim, I.H.; Sun, W.; Kim, H. Sustained expression of neuritin MRNA after repeated electroconvulsive stimulations in the rat hippocampal formation. Exp. Neurobiol. 2009, 18, 62. [Google Scholar] [CrossRef]
  35. Chandler, D.; Dragović, M.; Cooper, M.; Badcock, J.C.; Mullin, B.H.; Faulkner, D.; Wilson, S.G.; Hallmayer, J.; Howell, S.; Rock, D.; et al. Impact of neuritin 1 (NRN1) polymorphisms on fluid intelligence in schizophrenia. Am. J. Med. Genet. B Neuropsychiatr. Genet. 2010, 153, 428–437. [Google Scholar] [CrossRef]
  36. Fatjó-Vilas, M.; Prats, C.; Pomarol-Clotet, E.; Lázaro, L.; Moreno, C.; González-Ortega, I.; Lera-Miguel, S.; Miret, S.; Muñoz, M.J.; Ibáñez, I.; et al. Involvement of NRN1 gene in schizophrenia-spectrum and bipolar disorders and its impact on age at onset and cognitive functioning. World J. Biol. Psychiatry 2016, 17, 129–139. [Google Scholar] [CrossRef] [Green Version]
  37. Birnbaum, R.; Weinberger, D.R. Functional neuroimaging and schizophrenia: A view towards effective connectivity modeling and polygenic risk. Dialogues Clin. Neurosci. 2013, 15, 279–289. [Google Scholar] [CrossRef]
  38. Bora, E.; Lin, A.; Wood, S.J.; Yung, A.R.; McGorry, P.D.; Pantelis, C. Cognitive deficits in youth with familial and clinical high risk to psychosis: A systematic review and meta-analysis. Acta Psychiatr. Scand. 2014, 130, 1–15. [Google Scholar] [CrossRef]
  39. Zhang, R.; Picchioni, M.; Allen, P.; Toulopoulou, T. Working memory in unaffected relatives of patients with schizophrenia: A meta-analysis of functional magnetic resonance imaging studies. Schizophr. Bull. 2016, 42, 1068–1077. [Google Scholar] [CrossRef] [Green Version]
  40. Repovš, G.; Barch, D.M. Working memory related brain network connectivity in individuals with schizophrenia and their siblings. Front. Hum. Neurosci. 2012, 6, 137. [Google Scholar] [CrossRef] [Green Version]
  41. Silver, H.; Feldman, P.; Bilker, W.; Gur, R.C. Working memory deficit as a core neuropsychological dysfunction in schizophrenia. Am. J. Psychiatry 2003, 160, 1809–1816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Van Erp, T.G.M.; Walton, E.; Hibar, D.P.; Schmaal, L.; Jiang, W.; Glahn, D.C.; Pearlson, G.D.; Yao, N.; Fukunaga, M.; Hashimoto, R.; et al. Cortical brain abnormalities in 4474 individuals with schizophrenia and 5098 control subjects via the enhancing neuro imaging genetics through meta analysis (ENIGMA) consortium. Biol. Psychiatry 2018, 84, 644–654. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Bittner, R.A.; Linden, D.E.J.; Roebroeck, A.; Härtling, F.; Rotarska-Jagiela, A.; Maurer, K.; Goebel, R.; Singer, W.; Haenschel, C. The when and where of working memory dysfunction in early-onset schizophrenia—A functional magnetic resonance imaging study. Cereb. Cortex 2015, 25, 2494–2506. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kyriakopoulos, M.; Dima, D.; Roiser, J.P.; Corrigall, R.; Barker, G.J.; Frangou, S. Abnormal functional activation and connectivity in the working memory network in early-onset schizophrenia. J. Am. Acad. Child Adolesc. Psychiatry 2012, 51, 911–920.e2. [Google Scholar] [CrossRef] [PubMed]
  45. Loeb, F.F.; Zhou, X.; Craddock, K.E.S.; Shora, L.; Broadnax, D.D.; Gochman, P.; Clasen, L.S.; Lalonde, F.M.; Berman, R.A.; Berman, K.F.; et al. Reduced functional brain activation and connectivity during a working memory task in childhood-onset schizophrenia. J. Am. Acad. Child Adolesc. Psychiatry 2018, 57, 166–174. [Google Scholar] [CrossRef] [PubMed]
  46. Pauly, K.; Seiferth, N.Y.; Kellermann, T.; Backes, V.; Vloet, T.D.; Shah, N.J.; Schneider, F.; Habel, U.; Kircher, T.T. Cerebral dysfunctions of emotion—Cognition interactions in adolescent-onset schizophrenia. J. Am. Acad. Child Adolesc. Psychiatry 2008, 47, 1299–1310. [Google Scholar] [CrossRef] [PubMed]
  47. Thormodsen, R.; Jensen, J.; Holmèn, A.; Juuhl-Langseth, M.; Emblem, K.E.; Andreassen, O.A.; Rund, B.R. Prefrontal hyperactivation during a working memory task in early-onset schizophrenia spectrum disorders: An FMRI study. Psychiatry Res. Neuroimaging 2011, 194, 257–262. [Google Scholar] [CrossRef]
  48. White, T.; Hongwanishkul, D.; Schmidt, M. Increased anterior cingulate and temporal lobe activity during visuospatial working memory in children and adolescents with schizophrenia. Schizophr. Res. 2011, 125, 118–128. [Google Scholar] [CrossRef] [Green Version]
  49. Solmi, M.; Radua, J.; Olivola, M.; Croce, E.; Soardo, L.; Salazar de Pablo, G.; Il Shin, J.; Kirkbride, J.B.; Jones, P.; Kim, J.H.; et al. Age at onset of mental disorders worldwide: Large-scale meta-analysis of 192 epidemiological studies. Mol. Psychiatry 2021, 27, 1–15. [Google Scholar] [CrossRef]
  50. Lindholm, E.; Ekholm, B.; Shaw, S.; Jalonen, P.; Johansson, G.; Pettersson, U.; Sherrington, R.; Adolfsson, R.; Jazin, E. A schizophrenia-susceptibility locus at 6q25, in one of the world’s largest reported pedigrees. Am. J. Hum. Genet. 2001, 69, 96–105. [Google Scholar] [CrossRef] [Green Version]
  51. Maziade, M.; Roy, M.-A.; Rouillard, É.; Bissonnette, L.; Fournier, J.-P.; Roy, A.; Garneau, Y.; Montgrain, N.; Potvin, A.; Cliche, D.; et al. A search for specific and common susceptibility loci for schizophrenia and bipolar disorder: A linkage study in 13 target chromosomes. Mol. Psychiatry 2001, 6, 684–693. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Moises, H.W.; Yang, L.; Kristbjarnarson, H.; Wiese, C.; Byerley, W.; Macciardi, F.; Arolt, V.; Blackwood, D.; Liu, X.; Sjögren, B.; et al. An international two–stage genome–wide search for schizophrenia susceptibility genes. Nat. Genet. 1995, 11, 321–324. [Google Scholar] [CrossRef] [PubMed]
  53. Straub, R.E.; MacLean, C.J.; Ma, Y.; Webb, B.T.; Myakishev, M.V.; Harris-Kerr, C.; Wormley, B.; Sadek, H.; Kadambi, B.; O’Neill, F.A.; et al. Genome-wide scans of three independent sets of 90 Irish Multiplex schizophrenia families and follow-up of selected regions in all families provides evidence for multiple susceptibility genes. Mol. Psychiatry 2002, 7, 542–559. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Fanous, A.H.; Neale, M.C.; Webb, B.T.; Straub, R.E.; Amdur, R.L.; O’Neill, F.A.; Walsh, D.; Riley, B.P.; Kendler, K.S. A genome-wide scan for modifier Loci in schizophrenia. Am. J. Med. Genet. B Neuropsychiatr. Genet. 2007, 144, 589–595. [Google Scholar] [CrossRef]
  55. Pidsley, R.; Viana, J.; Hannon, E.; Spiers, H.; Troakes, C.; Al-Saraj, S.; Mechawar, N.; Turecki, G.; Schalkwyk, L.C.; Bray, N.J.; et al. Methylomic profiling of human brain tissue supports a neurodevelopmental origin for schizophrenia. Genome Biol. 2014, 15, 483. [Google Scholar] [CrossRef]
  56. Bähner, F.; Meyer-Lindenberg, A. Hippocampal–prefrontal connectivity as a translational phenotype for schizophrenia. Eur. Neuropsychopharmacol. 2017, 27, 93–106. [Google Scholar] [CrossRef]
  57. Riley, M.R.; Constantinidis, C. Role of prefrontal persistent activity in working memory. Front. Syst. Neurosci. 2016, 9, 181. [Google Scholar] [CrossRef] [Green Version]
  58. Thomas, M.L.; Duffy, J.R.; Swerdlow, N.; Light, G.A.; Brown, G.G. Detecting the inverted-U in FMRI studies of schizophrenia: A comparison of three analysis methods. J. Int. Neuropsychol. Soc. 2022, 28, 258–269. [Google Scholar] [CrossRef]
  59. Prats, C.; Arias, B.; Ortet, G.; Ibáñez, M.I.; Moya, J.; Pomarol-Clotet, E.; Fañanás, L.; Fatjó-Vilas, M. Role of neurotrophins in depressive symptoms and executive function: Association analysis of NRN1 gene and its interaction with BDNF gene in a non-clinical sample. J. Affect. Disord. 2017, 211, 92–98. [Google Scholar] [CrossRef] [Green Version]
  60. Lenroot, R.K.; Giedd, J.N. Brain development in children and adolescents: Insights from anatomical magnetic resonance imaging. Neurosci. Biobehav. Rev. 2006, 30, 718–729. [Google Scholar] [CrossRef]
  61. Fair, D.A.; Dosenbach, N.U.F.; Church, J.A.; Cohen, A.L.; Brahmbhatt, S.; Miezin, F.M.; Barch, D.M.; Raichle, M.E.; Petersen, S.E.; Schlaggar, B.L. Development of distinct control networks through segregation and integration. Proc. Natl. Acad. Sci. USA 2007, 104, 13507–13512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Finn, A.S.; Sheridan, M.A.; Kam, C.L.H.; Hinshaw, S.; D’Esposito, M. Longitudinal evidence for functional specialization of the neural circuit supporting working memory in the human brain. J. Neurosci. 2010, 30, 11062–11067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Andre, J.; Picchioni, M.; Zhang, R.; Toulopoulou, T. Working memory circuit as a function of increasing age in healthy adolescence: A systematic review and meta-analyses. NeuroImage Clin. 2015, 12, 940–948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Hartshorne, J.K.; Germine, L.T. When does cognitive functioning peak? The asynchronous rise and fall of different cognitive abilities across the lifespan. Psychol. Sci. 2015, 26, 433–443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Bowen, E.F.W.; Burgess, J.L.; Granger, R.; Kleinman, J.E.; Rhodes, C.H. DLPFC Transcriptome defines two molecular subtypes of schizophrenia. Transl. Psychiatry 2019, 9, 1–10. [Google Scholar] [CrossRef] [Green Version]
  66. Singh-Taylor, A.; Molet, J.; Jiang, S.; Korosi, A.; Bolton, J.L.; Noam, Y.; Simeone, K.; Cope, J.; Chen, Y.; Mortazavi, A.; et al. NRSF-Dependent epigenetic mechanisms contribute to programming of stress-sensitive neurons by neonatal experience, promoting resilience. Mol. Psychiatry 2018, 23, 648–657. [Google Scholar] [CrossRef] [Green Version]
  67. Ohi, K.; Hashimoto, R.; Yasuda, Y.; Kiribayashi, M.; Iike, N.; Yoshida, T.; Azechi, M.; Ikezawa, K.; Takahashi, H.; Morihara, T.; et al. TATA box-binding protein gene is associated with risk for schizophrenia, age at onset and prefrontal function. Genes Brain Behav. 2009, 8, 473–480. [Google Scholar] [CrossRef]
  68. Zhang, J.; Chang, L.; Pu, Y.; Hashimoto, K. Abnormal expression of colony stimulating factor 1 receptor (CSF1R) and transcription factor PU.1 (SPI1) in the spleen from patients with major psychiatric disorders: A role of brain–spleen axis. J. Affect. Disord. 2020, 272, 110–115. [Google Scholar] [CrossRef]
  69. Lee, S.H.; DeCandia, T.R.; Ripke, S.; Yang, J.; Sullivan, P.F.; Goddard, M.E.; Keller, M.C.; Visscher, P.M.; Wray, N.R. Estimating the proportion of variation in susceptibility to schizophrenia captured by common SNPs. Nat. Genet. 2012, 44, 247–250. [Google Scholar] [CrossRef] [Green Version]
  70. Singh, T.; Poterba, T.; Curtis, D.; Akil, H.; Al Eissa, M.; Barchas, J.D.; Bass, N.; Bigdeli, T.B.; Breen, G.; Bromet, E.J.; et al. Rare coding variants in ten genes confer substantial risk for schizophrenia. Nature 2022, 604, 509–516. [Google Scholar] [CrossRef]
  71. Fusar-Poli, P. Voxel-wise meta-analysis of FMRI studies in patients at clinical high risk for psychosis. J. Psychiatry Neurosci. 2012, 37, 106–112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Radua, J.; Borgwardt, S.; Crescini, A.; Mataix-Cols, D.; Meyer-Lindenberg, A.; McGuire, P.K.; Fusar-Poli, P. Multimodal meta-analysis of structural and functional brain changes in first episode psychosis and the effects of antipsychotic medication. Neurosci. Biobehav. Rev. 2012, 36, 2325–2333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Geller, B.; Zimerman, B.; Williams, M.; Bolhofner, K.; Craney, J.L.; Delbello, M.P.; Soutullo, C. Reliability of the Washington University in St. Louis kiddie schedule for affective disorders and schizophrenia (WASH-U-KSADS) mania and rapid cycling sections. J. Am. Acad. Child Adolesc. Psychiatry 2001, 40, 450–455. [Google Scholar] [CrossRef]
  74. Andreasen, N.C. The comprehensive assessment of symptoms and history (CASH): An instrument for assessing diagnosis and psychopathology. Arch. Gen. Psychiatry 1992, 49, 615. [Google Scholar] [CrossRef]
  75. First, M.; Spitzer, R.; Gibbon, M.; Williams, J. Structured clinical interview for DSM-IV-TR Axis I disorders, research version, non-patient edition. In (SCID-I/P); New York State Psychiatric Institute: New York, NY, USA, 2002. [Google Scholar]
  76. Perkins, D.O.; Leserman, J.; Jarskog, L.F.; Graham, K.; Kazmer, J.; Lieberman, J.A. Characterizing and dating the onset of symptoms in psychotic illness: The symptom onset in schizophrenia (SOS) inventory. Schizophr. Res. 2000, 44, 1–10. [Google Scholar] [CrossRef]
  77. Gershon, E.S.; DeLisi, L.E.; Hamovit, J.; Nurnberger, J.I., Jr.; Maxwell, M.E.; Schreiber, J.; Dauphinais, D.; Dingman, C.W., II; Guroff, J.J. A controlled family study of chronic psychoses: Schizophrenia and schizoaffective disorder. Arch. Gen. Psychiatry 1988, 45, 328–336. [Google Scholar] [CrossRef] [PubMed]
  78. Gomar, J.J.; Ortiz-Gil, J.; McKenna, P.J.; Salvador, R.; Sans-Sansa, B.; Sarró, S.; Guerrero, A.; Pomarol-Clotet, E. Validation of the word accentuation test (TAP) as a means of estimating premorbid IQ in spanish speakers. Schizophr. Res. 2011, 128, 175–176. [Google Scholar] [CrossRef]
  79. Zhang, F.; Lupski, J.R. Non-coding genetic variants in human disease. Hum. Mol. Genet. 2015, 24, R102–R110. [Google Scholar] [CrossRef] [Green Version]
  80. Fullard, J.F.; Giambartolomei, C.; Hauberg, M.E.; Xu, K.; Voloudakis, G.; Shao, Z.; Bare, C.; Dudley, J.T.; Mattheisen, M.; Robakis, N.K.; et al. Open chromatin profiling of human postmortem brain infers functional roles for non-coding schizophrenia Loci. Hum. Mol. Genet. 2017, 26, 1942–1951. [Google Scholar] [CrossRef] [Green Version]
  81. Ward, L.D.; Kellis, M. HaploReg: A resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2012, 40, D930–D934. [Google Scholar] [CrossRef] [Green Version]
  82. Boyle, E.A.; Li, Y.I.; Pritchard, J.K. An expanded view of complex traits: From polygenic to omnigenic. Cell 2017, 169, 1177–1186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Dong, S.; Boyle, A.P. Predicting functional variants in enhancer and promoter elements using regulomeDB. Hum. Mutat. 2019, 40, 1292–1298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Smith, S.M.; Jenkinson, M.; Woolrich, M.W.; Beckmann, C.F.; Behrens, T.E.J.; Johansen-Berg, H.; Bannister, P.R.; De Luca, M.; Drobnjak, I.; Flitney, D.E.; et al. Advances in functional and structural mr image analysis and implementation as FSL. NeuroImage 2004, 23, S208–S219. [Google Scholar] [CrossRef] [Green Version]
  85. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Worsley, K.J. Statistical Analysis of Activation Images. In Functional Magnetic Resonance Imaging; Oxford University Press: Oxford, UK, 2001; ISBN 978-0-19-263071-1. [Google Scholar]
  88. Egli, T.; Coynel, D.; Spalek, K.; Fastenrath, M.; Freytag, V.; Heck, A.; Loos, E.; Auschra, B.; Papassotiropoulos, A.; de Quervain, D.J.-F.; et al. Identification of two distinct working memory-related brain networks in healthy young adults. eNeuro 2018, 5, ENEURO.0222-17.2018. [Google Scholar] [CrossRef]
  89. Wilhelm, O.; Hildebrandt, A.; Oberauer, K. What is working memory capacity, and how can we measure it? Front. Psychol. 2013, 4, 433. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (A) Axial view of the brain showing the significant cluster derived from the diagnosis x NRN1 HAP678 GTT analysis in the 2-back vs. 1-back contrast. A sagittal view with the marks of the cross slices is also included. The right side of the image represents the right side of the brain. The MNI coordinates are given for the shown slices. Units of the bar correspond to the β values of the regression, standardised to Z scores. (B) Bar plots with the cluster mean activity (estimated marginal means and ±2 standard errors (se)) for healthy subjects (HS; left, non-carriers: n = 27, carriers: n = 10) and subjects with early-onset schizophrenia (EO; right, non-carriers: n = 19, carriers: n = 20).
Figure 1. (A) Axial view of the brain showing the significant cluster derived from the diagnosis x NRN1 HAP678 GTT analysis in the 2-back vs. 1-back contrast. A sagittal view with the marks of the cross slices is also included. The right side of the image represents the right side of the brain. The MNI coordinates are given for the shown slices. Units of the bar correspond to the β values of the regression, standardised to Z scores. (B) Bar plots with the cluster mean activity (estimated marginal means and ±2 standard errors (se)) for healthy subjects (HS; left, non-carriers: n = 27, carriers: n = 10) and subjects with early-onset schizophrenia (EO; right, non-carriers: n = 19, carriers: n = 20).
Ijms 23 07456 g001
Figure 2. Bar plots with mean performance (estimated marginal means and ±2 standard errors(se)) for healthy subjects (HS; left, non-carriers: n = 27, carriers: n = 10) and subjects with early-onset schizophrenia (EO; right, non-carriers: n = 19, carriers: n = 20) by NRN1 HAP678 GTT.
Figure 2. Bar plots with mean performance (estimated marginal means and ±2 standard errors(se)) for healthy subjects (HS; left, non-carriers: n = 27, carriers: n = 10) and subjects with early-onset schizophrenia (EO; right, non-carriers: n = 19, carriers: n = 20) by NRN1 HAP678 GTT.
Ijms 23 07456 g002
Table 1. Significant genetic association results within early-onset (EO) family-based and EO case–control samples. At the top row, for the family-based approach, there are the transmitted and not-transmitted haplotype (HAP) counts from heterozygous parents to affected offspring for family-based analyses. Below, for the case–control approach, the frequency (%) in healthy subjects (HS) and schizophrenia spectrum disorders (SSD) are given with the risk genotype placed last. The odds ratio (OR) associated with the genotype and the confidence interval (CI 95%) are also reported. The empirical p-values obtained after 10,000 permutation procedures (pperm) for the Transmission Disequilibrium Test (TDT) or the logistic regression (additive model) are shown.
Table 1. Significant genetic association results within early-onset (EO) family-based and EO case–control samples. At the top row, for the family-based approach, there are the transmitted and not-transmitted haplotype (HAP) counts from heterozygous parents to affected offspring for family-based analyses. Below, for the case–control approach, the frequency (%) in healthy subjects (HS) and schizophrenia spectrum disorders (SSD) are given with the risk genotype placed last. The odds ratio (OR) associated with the genotype and the confidence interval (CI 95%) are also reported. The empirical p-values obtained after 10,000 permutation procedures (pperm) for the Transmission Disequilibrium Test (TDT) or the logistic regression (additive model) are shown.
SNPsHaplotypeTransmitted EO SSDNot Transmitted EO SSDOR (CI 95%)TDT; pperm
HAP678GCT13270.48 (0.25–0.93)4.90; 0.03
SNPsGenotypes
Haplotypes
Frequency EO SSDFrequency HSOR (CI 95%)Wald; pperm
SNP6TT/TG/GG11 (0.13)/40 (0.48)/33 (0.39)31 (0.26)/58 (0.49)/30 (0.25)1.68 (1.01–2.57)2.39; 0.02 a,b
SNP7CC/CT/TT43 (0.50)/35 (0.41)/8 (0.09)78 (0.66)/35 (0.29)/6 (0.05)1.69 (1.06–2.71)2.19; 0.03 b
SNP8CC/CT/TT14 (0.17)/41 (0.49)/29 (0.35)33 (0.29)/58 (0.51)/23 (0.20)1.66 (1.08–2.56)2.31; 0.02 a,c
HAP678TCC0.370.500.59 (0.39–0.89)6.44; 0.01
HAP678GTT0.300.201.70 (1.08–2.67)5.28; 0.02
a The genotypic model was also significant (pperm < 0.05). b The dominant model was also significant (pperm < 0.05). c The recessive model was also significant (pperm < 0.05).
Table 2. Sociodemographic and clinical information for the family-based and case-control samples included in the genetic association and neuroimaging analyses. Data for patients are given separately for subjects with early-onset (EO) and adult-onset (AO) schizophrenia spectrum disorders (SSD). Number (percentage in brackets) are shown for qualitative variables. Mean scores (standard deviation in brackets) are provided for quantitative variables. Illness duration refers to years. The psychopathology was assessed using the Positive and Negative Symptoms Scale (PANSS). Treatment was defined by chlorpromazine equivalence (CPZE). Those not significant values are not reported (n.s.).
Table 2. Sociodemographic and clinical information for the family-based and case-control samples included in the genetic association and neuroimaging analyses. Data for patients are given separately for subjects with early-onset (EO) and adult-onset (AO) schizophrenia spectrum disorders (SSD). Number (percentage in brackets) are shown for qualitative variables. Mean scores (standard deviation in brackets) are provided for quantitative variables. Illness duration refers to years. The psychopathology was assessed using the Positive and Negative Symptoms Scale (PANSS). Treatment was defined by chlorpromazine equivalence (CPZE). Those not significant values are not reported (n.s.).
Sample 1:
Family-Based a
(n = 453)
AO Offspring
(n = 71)
EO Offspring
(n = 80)
AO Parents
(n = 142)
EO Parents
(n = 160)
Male 58 (81.70)54 (67.50)n.s.71 (50.00)80 (50.00)n.s.
Age at interview27.45 (5.03)18.04 (4.94)t = −11.51, p < 0.00150.04 (7.97)58.34 (8.44)t = −7.27, p < 0.001
Age at onset23.24 (4.24)15.41 (2.12) dt = −13.35, p < 0.001
Sample 2:
Case-control b
(n = 345)
AO Subjects
(n = 138)
EO Subjects
(n = 87)
Healthy Subjects
(n = 120)
Male93 (67.40)67 (77.00)n.s.60 (50.00)χ2 = 17.13, p < 0.001
Age at interview41.97 (10.03)39.79 (10.87)n.s.38.24 (11.21)F = 3.45; p = 0.033 e
Age at onset25.12 (5.86)16.38 (2.00)t = 16.11, p < 0.001
Sample 3: Neuroimaging c
(n = 117)
AO Subjects
(n = 39)
EO Subjects
(n = 39)
Healthy Subjects
(n = 39)
Male37 (94.87)37 (94.87)n.s.37 (94.87)n.s.
Age at interview39.49 (1.90)39.30 (1.87)n.s.38.43 (1.78)n.s.
Age at onset24.56 (0.80)16.85 (0.26)t = 9.17; p < 0.001 --
Illness duration14.92 (11.01)22.46 (11.31)t = −2.98; p = 0.004--
PANSS total68.72 (20.46)80.05 (21.11)t = −2.60; p = 0.011--
CPZE367.01 (188.83)633.66 (304.39)t = −4.28; p < 0.001--
a The different SSD diagnoses were equally distributed between the AO and EO offspring (χ2 = 1.271, p = 0.736): Schizophrenia (AO n = 41 (57.75%); EO n = 47 (58.75%)), Schizophreniform (AO n = 13 (18.31%); EO n = 10 (12.50%)), Schizoaffective (AO n = 6 (8.45%); EO 7 (8.75%)) and Psychosis not otherwise specified (AO n = 11 (15.4%); EO n = 16 (20.00%)). b The different diagnoses were equally distributed between the AO and EO cases (χ2 = 1.558, p = 0.212): Schizophrenia (AO n = 130 (84.42%); EO n = 84 (90.32%)) and Schizoaffective (AO n = 24 (15.58%); EO n = 9 (9.80%)). c Sub-set of individuals coming from Sample 3. Patients included in the neuroimaging sample were all diagnosed with schizophrenia. d Information was available for 60% of patients but all were drawn from child and adolescent units which allowed their classification as EO. e The post-hoc analyses showed significant differences only between healthy subjects and subjects with AO SSD (p = 0.010).
Table 3. Information on the Single Nucleotide Polymorphisms (SNPs) at Neuritin 1 gene included in this study (NRN1, chromosome 6p25.1, from 5,997,999 to 6,007,605 bp, UCSC Genome Browser on Human Assembly GRCh38/hg38, http://genome.ucsc.edu/cgi-bin/hgTracks, accessed on 1 April 2022). The table includes dbSNP number, the chromosome and gene position, the alleles of each SNP, the minor allele frequency (MAF; described for all and EUR populations in the 1000 Genomes Project and the MAF observed in each sample included in the present study) and the functional score according to the Regulome Database.
Table 3. Information on the Single Nucleotide Polymorphisms (SNPs) at Neuritin 1 gene included in this study (NRN1, chromosome 6p25.1, from 5,997,999 to 6,007,605 bp, UCSC Genome Browser on Human Assembly GRCh38/hg38, http://genome.ucsc.edu/cgi-bin/hgTracks, accessed on 1 April 2022). The table includes dbSNP number, the chromosome and gene position, the alleles of each SNP, the minor allele frequency (MAF; described for all and EUR populations in the 1000 Genomes Project and the MAF observed in each sample included in the present study) and the functional score according to the Regulome Database.
SNPsChromosome
Position
Gene
Position
Alleles
(Minor/Major)
MAF
(All/Eur)
Family-Based
MAF
Case-Control
MAF
RegulomeDB Score a
SNP1rs22088705,992,257intergenicG/A0.33/0.340.330.330.61
SNP2rs123331175,994,759intergenicT/C0.35/0.400.430.380.61
SNP3rs5821866,001,148downstreamA/G0.45/0.620.610.400.61
SNP4rs6456496,004,726intronicC/G0.45/0.640.640.380.61
SNP5rs5822626,007,758intronicC/G0.30/0.480.280.270.70
SNP6rs37631806,009,615upstreamT/G0.40/0.460.450.430.63
SNP7rs104843206,010,204upstreamT/C0.15/0.220.260.240.16
SNP8rs49601556,010,306upstreamT/C0.43/0.490.500.490.13
SNP9rs93790026,012,158intergenicG/T0.29/0.420.240.260.13
SNP10rs94058906,012,488intergenicC/T0.31/0.380.280.330.18
SNP11rs14751576,016,936intergenicG/A0.16/0.170.160.160.18
a This score ranges from 0 to 1, with 1 being most likely to be a regulatory variant.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Almodóvar-Payá, C.; Guardiola-Ripoll, M.; Giralt-López, M.; Gallego, C.; Salgado-Pineda, P.; Miret, S.; Salvador, R.; Muñoz, M.J.; Lázaro, L.; Guerrero-Pedraza, A.; et al. NRN1 Gene as a Potential Marker of Early-Onset Schizophrenia: Evidence from Genetic and Neuroimaging Approaches. Int. J. Mol. Sci. 2022, 23, 7456. https://doi.org/10.3390/ijms23137456

AMA Style

Almodóvar-Payá C, Guardiola-Ripoll M, Giralt-López M, Gallego C, Salgado-Pineda P, Miret S, Salvador R, Muñoz MJ, Lázaro L, Guerrero-Pedraza A, et al. NRN1 Gene as a Potential Marker of Early-Onset Schizophrenia: Evidence from Genetic and Neuroimaging Approaches. International Journal of Molecular Sciences. 2022; 23(13):7456. https://doi.org/10.3390/ijms23137456

Chicago/Turabian Style

Almodóvar-Payá, Carmen, Maria Guardiola-Ripoll, Maria Giralt-López, Carme Gallego, Pilar Salgado-Pineda, Salvador Miret, Raymond Salvador, María J. Muñoz, Luisa Lázaro, Amalia Guerrero-Pedraza, and et al. 2022. "NRN1 Gene as a Potential Marker of Early-Onset Schizophrenia: Evidence from Genetic and Neuroimaging Approaches" International Journal of Molecular Sciences 23, no. 13: 7456. https://doi.org/10.3390/ijms23137456

APA Style

Almodóvar-Payá, C., Guardiola-Ripoll, M., Giralt-López, M., Gallego, C., Salgado-Pineda, P., Miret, S., Salvador, R., Muñoz, M. J., Lázaro, L., Guerrero-Pedraza, A., Parellada, M., Carrión, M. I., Cuesta, M. J., Maristany, T., Sarró, S., Fañanás, L., Callado, L. F., Arias, B., Pomarol-Clotet, E., & Fatjó-Vilas, M. (2022). NRN1 Gene as a Potential Marker of Early-Onset Schizophrenia: Evidence from Genetic and Neuroimaging Approaches. International Journal of Molecular Sciences, 23(13), 7456. https://doi.org/10.3390/ijms23137456

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop