ORIGINAL RESEARCH article
Sec. Neurological Biomarkers
Volume 14 - 2023 | https://doi.org/10.3389/fneur.2023.1145737
Evidence of polygenic regulation of the physiological presence of neurofilament light chain in human serum
- 1Department of Genetic Epidemiology, Institute of Human Genetics, University of Münster, Münster, Germany
- 2Department of Neurology, Medical University of Graz, Graz, Austria
- 3Institute for Medical Informatics, Statistics and Documentation, Medical University of Graz, Graz, Austria
- 4Neurologic Clinic and Polyclinic, MS Center and Research Center for Clinical Neuroimmunology and Neuroscience Basel (RC2NB), University Hospital of Basel, Basel, Switzerland
- 5Clinical Trial Unit, Department of Clinical Research, University Hospital of Basel, Basel, Switzerland
- 6Department of Neurology with Institute of Translational Neurology, University Hospital Münster, Münster, Germany
- 7Department of Biochemistry, Genetic Epidemiology and Statistical Genetics, Maastricht University, Maastricht, Netherlands
- 8Institute of Epidemiology and Social Medicine, University of Münster, Münster, Germany
Introduction: The measurement of neurofilament light chain (NfL) in blood is a promising biomarker of neurological injury and disease. We investigated the genetic factors that underlie serum NfL levels (sNfL) of individuals without neurological conditions.
Methods: We performed a discovery genome-wide association study (GWAS) of sNfL in participants of the German BiDirect Study (N = 1,899). A secondary GWAS for meta-analysis was performed in a small Austrian cohort (N = 287). Results from the meta-analysis were investigated in relation with several clinical variables in BiDirect.
Results: Our discovery GWAS identified 12 genomic loci at the suggestive threshold ((p < 1 × 10−5). After meta-analysis, 7 loci were suggestive of an association with sNfL. Genotype-specific differences in sNfL were observed for the lead variants of meta-analysis loci (rs34523114, rs114956339, rs529938, rs73198093, rs34372929, rs10982883, and rs1842909) in BiDirect participants. We identified potential associations in meta-analysis loci with markers of inflammation and renal function. At least 6 protein-coding genes (ACTG2, TPRKB, DMXL1, COL23A1, NAT1, and RIMS2) were suggested as genetic factors contributing to baseline sNfL levels.
Discussion: Our findings suggest that polygenic regulation of neuronal processes, inflammation, metabolism and clearance modulate the variability of NfL in the circulation. These could aid in the interpretation of sNfL measurements in a personalized manner.
Neurofilament light chain (NfL) is a subunit of neurofilaments (NFs), cytoskeletal components found exclusively in neurons and particularly abundant in axons. NfL is a major component of the backbone of NFs in the central and peripheral nervous systems (1). Axonal damage and neuronal death due to neurological diseases, including those of inflammatory, neurodegenerative, traumatic and cerebrovascular nature, result in NfL release into the cerebrospinal fluid (CSF) and blood. Recent technological advances in immunoassay detection have enabled the accurate measurement of the small amounts of NfL that reach the circulation, facilitating its application as a universal peripheral biomarker of the presence and progression of neurological conditions, and of treatment responses (1–3). Therefore, investigating the factors that influence concentrations of NfL in the periphery becomes crucial for the interpretation of results. To date, it has been demonstrated that NfL serum levels (sNfL) increase with age (4) and potential confounding factors, such as body mass index and cardiovascular risk factors, have been suggested (5, 6).
Studies in population-based cohorts have shown a polygenic nature of numerous health-related serum biomarkers, including alanine transaminase (liver function), fibrinogen (clot formation) and glycated hemoglobin (type 2 diabetes mellitus), among many others. These findings can provide novel biological insights and facilitate disease diagnosis and stratification (7). Nevertheless, to our knowledge, no genetic associations with sNfL have been investigated. We hypothesized that the identification of genetic factors that modulate sNfL in physiological conditions will help interpretation on an individual basis, consequently improving the clinical applications of sNfL as a biomarker. To test our hypothesis, we performed a genome-wide association study (GWAS) and meta-analysis of sNfL in a total of 2,186 individuals of European descent without known neurological conditions, and correlated our findings with clinical data to identify potential sources of sNfL variability.
Subjects and methods
The BiDirect Study was initiated in 2009 as a prospective, observational study integrating three cohorts: (1) community-dwelling adults (control cohort), (2) patients with an acute depressive episode (depression cohort), and (3) patients who recently suffered from acute myocardial infarction (MI cohort). The study, whose principal goal is the exploration of the bidirectional relationship between depression and subclinical arteriosclerosis, recruited participants in the district of Münster, Germany, and carried out extensive phenotyping and follow-up of all cohorts in parallel. The study design and methods have been previously described in detail (8). Here, we included 1,899 BiDirect participants (977 males, 922 females; mean age: 52.1 ± 7.9) from the control (763), depression (851) and MI (285) cohorts.
The Austrian Stroke Prevention Family Study (ASPS-Fam) cohort represents an extension of the prospective, population-based ASPS (Austrian Stroke Prevention Study) on the effects of vascular risk factors in normal aging. ASPS was established in 1991 in the city of Graz, Austria (9). For ASPS-Fam, first-degree relatives of ASPS participants were invited to join the study. The study's composition and inclusion criteria have been described elsewhere (10, 11). Here, we included 287 ASPS-Fam participants (115 males, 172 females; mean age: 64.3 ± 10.6).
The basic descriptive information of the BiDirect and ASPS-Fam cohorts are shown in Table 1. Summary information on study design and composition can be found in the Supplementary material 1. All participants of the BiDirect and ASPS-Fam cohorts provided written informed consent. Methods were carried out in accordance with the ethical standards laid down in the updated version of the 1964 Declaration of Helsinki. The BiDirect Study was approved by the Ethics Committee of the University of Münster and the Westphalian Chamber of Physicians in Münster, North-Rhine-Westphalia, Germany. The ASPS-Fam protocol was approved by the Ethics Committee of the Medical University of Graz, Austria.
Serum measurements of NfL
Quantification of sNfL in BiDirect and ASPS-Fam was conducted at the University Hospital Basel, Switzerland, using the single molecule array (Simoa®) HDX analyzer (Quanterix, Lexington, MA, USA). In BiDirect participants, measurements of sNfL were obtained from non-fasting blood samples collected at the first visit, using the Simoa® NF-light Advantage Kit. In ASPS-Fam participants, sNfL measurement (Supplementary material 1) has been previously described in detail (4). The sNfL values obtained at initial assessment were log2-transformed and used for all analyses herein reported. Therefore, with sNfL in our findings, we actually refer to log2 sNfL.
Because it is known that sNfL concentrations increase during aging (4), we tested for age-adjusted sex- and cohort-dependent sNfL differences in BiDirect using analysis of covariance (ANCOVA). We also tested for sNfL correlations, using the Pearson's method, with markers of inflammation, renal and liver function, lipids, hormones and brain volumes derived from magnetic resonance imaging (MRI) data (106 clinical variables in total). All p < 0.05 values were considered statistically significant. Here, age represented the age at participant recruitment, when baseline phenotyping (s0) took place. Clinical variables coming from up to three subsequent follow-up visits were identified as time points s2, s4, and s6.
For BiDirect genotypes, genomic DNA was isolated from whole blood samples with EDTA using standard DNA extraction kits and procedures at the University of Münster. Genome-wide genotyping was performed with the Infinium PsychArray BeadChip v1 (Illumina) at Life&Brain GmbH (Bonn, Germany). Basic quality control (QC) was employed to remove samples and variants with high rates of missing data. This included removal of individuals with genotyping rate < 2%, cryptic relatedness (PI-HAT ≥1/16), sex mismatch and genetic outliers (distance in first two multidimensional scaling components >5 standard deviations from the mean), as well as the removal of variants with call rate < 2% and minor allele frequency (MAF) < 1%. Genotype imputation was performed with SHAPEIT (pre-phasing) (12) and IMPUTE2 (13) using the 1,000 Genomes Project, phase 3, European population reference panel (from here on, 1KG Reference Panel). Imputed variants were filtered for the INFO metric (≥0.8), MAF≥0.01 and Hardy-Weinberg equilibrium (HWE p ≥ 1 × 10−6). Individuals were further removed from the sample based on missing phenotypic data (age and baseline sNfL measurement). The final BiDirect GWAS dataset consisted of 5,597,244 genetic variants and 1,899 individuals.
For ASPS-Fam genotypes, genome-wide genotyping was performed with the Genome-Wide Human SNP Array 6.0 (Affymetrix). During the initial QC, variants with MAF < 0.05, HWE < 5 × 10−6 and low variant call rate (>2%) were excluded. Individuals with sex mismatch, cryptic relatedness, low sample call rate (>2%), a heterozygosity rate exceeding the mean ± 3 standard deviations and erroneous duplicates were removed. No genetic outliers were present. Genotype imputation was performed using the Michigan Imputation Server (14) and the 1KG Reference Panel.
Of note, genetic variants herein comprise single nucleotide polymorphisms (SNPs), as well as small insertions/deletions (indels) present in the datasets.
Screening for genetic associations with sNfL
We conducted a discovery GWAS in the BiDirect dataset under an additive regression model, adjusting for age, sex, cohort and the first 10 principal components. A secondary GWAS in the smaller ASPS-Fam dataset was performed independently at the Medical University of Graz and was adjusted for age, sex and the first 10 principal components. After harmonization of summary statistics from both studies, we performed a weighted meta-analysis of all overlapping variants with Rsq≥0.8 and MAF≥0.01 using Plink 1.9 (15). Variants with high heterogeneity between studies (I>40 and Q < 0.1) were subsequently neglected.
Definition of genomic loci for sNfL
For the discovery GWAS and the meta-analysis, we carried out downstream analyses on the FUMA GWAS platform (16) and defined genomic loci at the suggestive threshold of significance for genome-wide studies (p < 1 × 10−5), obtained variant annotations and identified the level of support for each signal. Linkage disequilibrium (LD) was defined by r2 ≥ 0.6 and a window of 500 kb, according to the 1KG Reference Panel. Subsequently, LD blocks were formed with variants under the suggestive threshold as lead variants, and containing all variants with p < 0.05 in the dataset that were in LD with the corresponding lead variants. Positional (gene) mapping was performed according to a maximum distance of 1 kb for the categories protein-coding, long non-coding RNA (lncRNA), non-coding RNA (ncRNA) and processed transcripts. Expression quantitative trait loci (eQTLs) were mapped using the BRAINEAC and GTEx v8 Brain databases. Only SNP-gene pairs with false discovery rate (FDR) < 0.05 were annotated.
Functional implications of suggested candidate genes
To inform the biological meaning of our findings, we created a protein-protein interaction (PPI) network using our suggested meta-analysis candidate genes as input. The network was generated with the Gene Set analysis tool of the ReactomeFIViz app for Cytoscape v.3.7.1 (17, 18). Linker proteins and functional interaction (FI) annotations were incorporated into the network (version 2018). In addition, we performed clustering of nodes, as well as enrichment analyses of pathways and gene ontology cellular components (GO_CC) for each network cluster. Gene sets with FDR < 0.05 were considered significantly enriched.
SNP heritability (h)
We calculated the proportion of variance in sNfL concentrations explained by our discovery GWAS in BiDirect using the GREML-LDMS (LD- and MAF-stratified GREML) method implemented in GCTA (19, 20). For all autosomal variants with MAF≥0.01 in the imputed dataset, we calculated the 200 kb segment-based LD scores, stratified variants according to LD scores of individual SNPs, computed one genetic relationship matrix for each quartile of the stratified variants, and performed a restricted maximum likelihood analysis using these four matrices. The variance explained was adjusted for the same covariates as the GWAS. SNP heritability from our meta-analysis summary statistics was calculated using LDSC software (21) with LD scores pre-computed in 1KG Reference Panel data, as suggested by the authors.
Screening for associations with clinical variables
For the lead variant of each loci resulting from our meta-analysis, we performed genotype-specific comparisons in BiDirect participants using an ANCOVA model adjusted for age. Moreover, for all variants within meta-analysis loci, we tested for associations with the same set of clinical variables used in the correlation analyses. These association tests were performed in the same manner as for baseline sNfL. The Benjamini-Hochberg method was used to correct for multiple comparisons (adjP).
Basic characterization of sNfL in BiDirect
Our initial characterization of sNfL in BiDirect found similar distributions of sNfL in the three cohorts (sNfL raw mean ± standard deviation values: control 9.49 ± 6.57, depression 9.24 ± 4.99, MI 11.76 ± 11.62; corresponding log2 values: control 2.15 ± 0.44, depression 2.13 ± 0.43, MI 2.29 ± 0.5; Figure 1A) and a positive association with age (p < 2 × 10−16, beta = 0.03), which was independent of the cohort (Figure 1B). Age-adjusted comparisons showed mean differences in sNfL levels between both patient cohorts (depression p = 8.2 × 10−5, MI p = 1.4 × 10−3) and the reference cohort, while no differences could be attributed to sex (p = 0.56) in this dataset (Figure 1C). Moreover, baseline (s0) sNfL correlated well with all other sNfL measurements (i.e., log- and non-transformed values from follow-up visits), and with markers of inflammation, and of the functions of kidneys, liver and thyroid glands (Supplementary material 2, Supplementary Table 1).
Figure 1. Serum neurofilament light chain (NfL) levels (log2-transformed) in BiDirect. Cohorts showed similar distributions of sNfL concentrations (A). There was a positive correlation of sNfL with age (B). Mean differences associated with the diagnostic group (BiDirect cohort: depression, myocardial infarction-MI-, and population-based control individuals), but not with sex, were observed (C). **p < 0.001, ***p < 0.0001.
Genetic associations with sNfL
We identified no genetic associations with sNfL surpassing the desired genome-wide significant threshold (p < 5 × 10−8). But, our observations reached a significance threshold commonly accepted for suggestive associations (p < 1 × 10−5) in GWASs. Therefore, we wished to further explore these suggestive findings from our GWAS and meta-analysis.
With our discovery GWAS in BiDirect (N = 1,899), we observed suggestive signals in 10 chromosomes (Figure 2A). Because the SNP2GENE tool integrates observations coming from GWAS summary statistics with information on LD structure coming from well-established reference panels to define lead variants and genomic loci, and can also be used to annotate an array of functional features for SNPs within the defined loci, we considered this tool to provide an appropriate means for the interpretation of our results. Twelve suggestive genomic loci for sNfL were defined through this analysis. These loci contained 13 lead variants (i.e., identified from independent variants and independent from each other at r2 ≥ 0.1), 14 independent signals (i.e., independent variants at the suggestive p-value threshold and independent from each other at r2 ≥ 0.6), and implicated a total of 246 genetic variants and of 18 mapped genes, from which 7 (CNTNAP5, NAT1, NATP, MTDH, RIMS2, VWA8, and RBFOX1) are protein-coding (Table 2, Supplementary material 2, Supplementary Table 2). The SNP heritability estimation performed with GCTA showed that this GWAS explained about 30% of the variance in sNfL (h= 0.299). However, the analysis also suggested that a larger sample size would be required to confidently detect the genetic component of sNfL (LRT = 2.4, p = 0.061).
Figure 2. Genetic associations with serum neurofilament light chain were only identified at the suggestive level. Manhattan and quantile-quantile (QQ) plots for the discovery genome-wide association study in BiDirect (N = 1,899) (A), and its meta-analysis with ASPS-Fam (N = 2,186) (B). Yellow lines in the Manhattan plots mark the suggestive threshold for genome-wide significance (p < 1 × 10−5).
Because the ASPS-Fam cohort has a small sample size and differences in its composition, in comparison with BiDirect, were evident, we chose not to seek validation of our findings in ASPS-Fam, but to use this cohort to carry out a meta-analysis with the aim to gain statistical power (N = 2,186). After performing a weighted meta-analysis and filtering out heterogeneous variants (i.e., variants with inconsistent effects), we applied again the SNP2GENE approach to extract a relevant interpretation of our results. Even with the addition of the ASPS-Fam cohort, we did not observe genomic variants reaching genome-wide significance (Figure 2B). Nevertheless, we were able to define 7 suggestive meta-analysis loci spanning 5 chromosomes, 144 variants and 8 mapped genes, including 6 protein-coding genes (ACTG2, TPRKB, DMXL1, COL23A1, NAT1, and RIMS2), that associated with sNfL levels in individuals without neurological conditions (Table 2, Supplementary material 2, 3, Supplementary Table 3, Supplementary Figures 1–7). In comparison with our discovery GWAS, meta-analysis loci represented the identification of 4 robust signals (i.e., meta-analysis loci that overlapped GWAS loci; meta-analysis loci #4–7 in chromosomes 8, 9, and 11), as well as the addition of 3 new signals (i.e., meta-analysis loci not found with the discovery GWAS; meta-analysis loci #1–3 in chromosomes 2 and 5). SNP heritability performed with LDSC in our sNfL meta-analysis was estimated to be about 5% (h = 0.0557). Nevertheless, we observed a low Chi2 statistic (mean Chi2 = 1.01) for this analysis, which may be due to the small sample size.
Investigation of biological context
The PPI network created with the protein-coding genes implicated by our meta-analysis loci was able to link 5/6 (exception of NAT1) genes by the incorporation of 9 linker proteins (Figure 3). Four small clusters were defined within this network, which illustrated the differential, yet interconnected functional properties between clusters. The most prominent pathways enriched in each cluster (Supplementary material 2, Supplementary Table 4) were related to cell signaling and organization of the extracellular matrix (lilac module: ACTG2, COL23A1, FURIN, MMP13, MMP16), senescence, inflammation and cell death (green module: AKT1, TP53, TP53RK, TPRKB), glucose and insulin metabolism (magenta module: MYH9, RAB8A, RIMS2), and immune processes (olive module: DMXL1, RICTOR). These pathways showed consistency with the associations observed between sNfL and clinical variables, including not only inflammation but also those related to thyroid and renal functions, and to blood lipids (e.g., Parathyroid hormone synthesis, secretion and action-FDR = 0.0086 in lilac module-; Thyroid hormone signaling pathway-FDR = 0.0052 in green module-; Plasma lipoprotein assembly, remodeling, and clearance-FDR = 0.03 in lilac module). Additionally, network modules were enriched for distinct cellular compartments (Supplementary material 2, Supplementary Table 5), mainly: extracellular matrix and Golgi (lilac module), cytoplasm and nucleus (green module), presynaptic cytoskeleton and transport vesicles (magenta module), and the RAVE (regulator of ATPase of vacuoles and endosomes) and TORC2 (target of rapamycin complex 2) complexes (olive module).
Figure 3. Protein-protein interaction network of mapped and brain expression quantitative trail loci genes implicated by the identified (meta-analysis) suggestive loci for serum neurofilament light chain. Circles denote input genes. Diamonds denote linker proteins. Colors denote network clusters, whose enrichments for pathways and gene ontology cellular compartments can be found in the Supplementary material 2, Supplementary Tables 4, 5.
Because none of the variants in our GWAS reached the common threshold accepted for genome-wide significance, we also dissected these associations. For all lead variants from our meta-analysis loci (rs34523114, rs114956339, rs529938, rs73198093, rs34372929, rs10982883 and rs1842909), we found significant differences in sNfL levels from BiDirect participants with different genotypes, particularly in those individuals with two copies of the effect/minor allele (AA genotype), as compared to those homozygous for the non-effect/major allele (BB genotype) (Table 3). With the exception of rs114956339 (p = 0.0016), we found no interactions for sNfL measurements between the genotypes of these variants and the diagnostic group (i.e., depression, MI and control).
Finally, we tested the associations of meta-analysis loci with clinical variables. None of these survived correction for multiple comparisons (adjP > 0.05), therefore we focused on the top signals (p < 0.05) from these tests only. At this threshold, we found evidence suggesting associations of meta-analysis loci with several clinical variables (Supplementary material 2, Supplementary Table 1). When prioritizing these by the integration of our results from genetic association and sNfL correlation tests, we identified overlaps for 18 variables from the clinical phenotypes (Table 4). These included markers of inflammation (interferon-α, and interleukins 6 and 1α), renal function (cystatin, creatinine, albumin and urea), liver and muscle function (lactate dehydrogenase and lipase), thyroid function (free thyroxine and free triiodothyronine), and blood lipids (HDL cholesterol and triglycerides). Noticeably, the index of comorbidity (which included stroke, leg thrombosis, peripheral artery disease, hypertension, MI, diabetes, depression, cancer, kidney and lung diseases, chronic arthritis, and Parkinson's disease) and gray matter volume (relative to total brain, coming from magnetic resonance imaging data) were also prioritized. Moreover, the associations with all sNfL measurements from follow-up visits remained suggested (Supplementary material 2, Supplementary Table 1).
Table 4. Prioritized clinical variables in BiDirect showed significant correlations with sNfL and the GWAS meta-analysis loci.
With the increasing interest in the clinical use of sNfL as a peripheral biomarker for the presence, progression and treatment response of neurological conditions in general, there is a need to define which biological factors contribute to physiological variations in sNfL concentrations. Previous studies have reported age, body mass index, blood volume, renal function (as measured by serum creatinine levels), hypertension and pregnancy may act as determinants of sNfL (3–6, 22). To some extent, we corroborated the association of sNfL with aging and renal function, and observed other physiological variables potentially associated with sNfL in the BiDirect study. Nevertheless, because of the small-effect interactions and overlaps at the genetic level that we observed, more studies will be necessary to clarify whether these findings may represent true determinants of serum sNfL levels or an epiphenomenon.
As our primary goal was to determine genetic factors that contribute to modulate sNfL concentrations, we performed a discovery GWAS and meta-analysis study in Europeans. Although we report here the findings from both analyses, we focused on the 7 suggestive loci resulting from our meta-analysis of the BiDirect and ASPS-Fam study populations to gain some biological insights on the implicated genomic regions. Results from our network analysis and overlapping genetic associations with a set of clinical variables show consistency. These highlighted particularly important roles for inflammation, lipids, thyroid hormones and vesicular transport. We also found in the literature, for all protein-coding mapped and/or any-tissue eQTL genes for variants in all of our meta-analysis loci, functions that are relevant for neuronal development and function. As neuronal processes may impact the release of NfL into the CSF and, consequently, its dissemination into peripheral blood, we focused on identifying potential roles of our meta-analysis loci in neuronal functions. However, as suggested by our analyses, it is possible that some variants contribute to regulate sNfL levels through effects on the body's metabolism and renal clearance.
In our study, NAT1, RIMS2 and DEC1 (meta-analysis loci #4–6, respectively) were the more robustly suggested candidate genes. The NAT1 (N-Acetyltransferase 1) protein forms an enzymatic complex with ARD1 (N-Alpha-Acetyltransferase 10, NatA Catalytic Subunit; NAA10 gene) that is required for neuronal differentiation and dendritic arborization (23, 24). The product of RIMS2 (Regulating Synaptic Membrane Exocytosis 2) functions as a Rab effector involved in synaptic membrane exocytosis (25). DEC1 (deleted in esophageal cancer 1, DELEC1), a lncRNA gene, is a candidate tumor suppressor (26), which means that it may regulate the cell cycle and other fundamental cellular processes.
Moreover, meta-analysis locus #1 mapped to ACTG2 (Actin Gamma 2) and implicated TPRKB (TP53RK Binding Protein) as a brain eQTL gene. Although the ACTG2 protein primarily localizes to the cytoskeleton of enteric smooth muscle, this gene has also been found downregulated during the chemical conversion of cultured human cortical astrocytes into neurons by treatment with small molecules (27), suggesting a role for ACTG2 in neuronal development. TPRKB is a subunit of the KEOPS (Kinase, Endopeptidase and Other Proteins of small Size) complex, which is required for the threonyl carbamoyl adenosine (t6A) transfer (t)RNA modification (28). An increasing number of reports link defects in these modifications to various neurodevelopmental disorders, suggesting a role in the development of the nervous system (29, 30). Additionally, when looking at any-tissue eQTL effects, genetic variants in meta-analysis locus #1 were found to regulate the expression of DCTN1 (Dynactin Subunit 1) and DGUOK (Deoxyguanosine Kinase). The product of DCTN1 is essential for the retrograde transport of vesicles and organelles along microtubules mediated by dynein. In neurons, it activates retrograde axonal transport and regulates microtubule stability (31, 32). On the other hand, DGUOK is a mitochondrial protein that may be involved in neuronal differentiation, as suggested by experiments in retinoic acid-induced differentiated neuronal-like cells (33).
Meta-analysis locus #2 mapped to DMXL1 (Dmx Like 1). In ngr1−/− mice, this gene was upregulated in axotomized corticospinal motor neurons 4 weeks after pyramidotomy (34), suggesting a role in axonal repair. Meta-analysis locus #3 mapped to COL23A1 (Collagen Type XXIII Alpha 1 Chain), whose dysregulated expression has been reported in different brain regions of mice with repeated experience of agonistic interactions (35). The work suggested the involvement of extracellular matrix remodeling (and of COL23A1) in the development of experimental psychopathologies. Although meta-analysis locus #7 did not map to protein-coding genes or showed eQTL effects on any in the brain datasets, we found variants in this locus with any-tissue eQTL effects on PTPN5 (Protein Tyrosine Phosphatase Non-Receptor Type 5). This gene regulates synaptic plasticity, and has been implicated in diverse neurological and psychiatric disorders (36–38).
We acknowledge important limitations of our study. First, the relatively small sample size limited the power to detect genetic associations at the genome-wide level and, therefore, to estimate SNP heritability. This was indeed reflected by the statistics from our heritability analyses. Second, serum samples from non-fasting study participants were used to determine sNfL concentrations. However, it is unknown if fasting status influences sNFL levels. Future assessments of sNfL levels in fasting and non-fasting blood should clarify whether this is a relevant factor for sNfL measurement. And, third, the nature of the study design of the sample populations included in the present study derived in an enrichment of individuals with depression, cardiovascular risk factors and cardiovascular disease. While most prior research focused on neurological conditions, recent studies have shown increased levels of sNfL in patients with cardiovascular or metabolic conditions and multimorbidity (39). In fact, we also showed this to be the case in the BiDirect study. To overcome this, we adjusted for these conditions and other confounding factors, including age. We expect that this is sufficient to adequately address condition-induced biases. Finally, we did not perform analyses within each condition cohort due to their limited sample sizes. Overall, we are positive that the future inclusion of appropriate population-based cohorts will help establish these and other genomic regions as genetic drivers of sNfL variations in individuals without neurological conditions. Further bioinformatics and functional studies should help to elucidate the biological relevance of our findings for sNfL measurements. The potential genetic and physiological factors associated with sNfL that were identified by our study warrant future investigations that will pave the way for an optimal application of sNfL as a marker of neuronal conditions.
Data availability statement
The data analyzed in this study is subject to the following licenses/restrictions: The summary statistics datasets generated in this study are available from the authors on reasonable request. The derived data supporting the conclusions presented in this article are included within the article and the corresponding additional files. Requests to access these datasets should be directed to MH-R, firstname.lastname@example.org.
The studies involving human participants were reviewed and approved by the Ethics Committee of the University of Münster and the Westphalian Chamber of Physicians in Münster, North Rhine-Westphalia, Germany and the Ethics Committee of the Medical University of Graz, Austria. The patients/participants provided their written informed consent to participate in this study.
MH-R: project design, data analysis, interpretation, and manuscript preparation. EH, MK, and RS: GWAS in ASPS-Fam. MS and KB: project design and critical revisions. KB: coordination of the BiDirect study. HW, AM, DL, PB, and JK: measurements of NfL. All authors contributed to the article and approved the submitted version.
ASPS-Fam was funded by the Austrian Science Fund (FWF, grant number PI904). The Medical University of Graz and the Steiermärkische Krankenanstalten Gesellschaft support the ASPS-Fam databank. The BiDirect Study was funded by the German Federal Ministry of Education and Research (BMBF, grant numbers FKZ-01ER0816 and -01ER1506).
The authors thank Till Andlauer for his support with genotype imputation. An early version of this manuscript was shared as a pre-print under a CC BY 4.0 License. This is accessible at: https://doi.org/10.21203/rs.3.rs-422221/v1.
Conflict of interest
DL is Chief Medical Officer at GeNeuro.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2023.1145737/full#supplementary-material
1. Yuan A, Rao MV, Veeranna Nixon RA. Neurofilaments and neurofilament proteins in health and disease. Cold Spring Harb Perspect Biol. (2017) 9:a018309. doi: 10.1101/cshperspect.a018309
2. Khalil M, Teunissen CE, Otto M, Piehl F, Sormani MP, Gattringer T, et al. Neurofilaments as biomarkers in neurological disorders. Nat Rev Neurol. (2018) 14:577–89. doi: 10.1038/s41582-018-0058-z
3. Barro C, Chitnis T, Weiner HL. Blood neurofilament light: a critical review of its application to neurologic disease. Ann Clin Transl Neurol. (2020) 7:2508–23. doi: 10.1002/acn3.51234
4. Khalil M, Pirpamer L, Hofer E, Voortman MM, Barro C, Leppert D, et al. Serum neurofilament light levels in normal aging and their association with morphologic brain changes. Nat Commun. (2020) 11:812. doi: 10.1038/s41467-020-14612-6
5. Manouchehrinia A, Piehl F, Hillert J, Kuhle J, Alfredsson L, Olsson T, et al. Confounding effect of blood volume and body mass index on blood neurofilament light chain levels. Ann Clin Transl Neurol. (2020) 7:139–43. doi: 10.1002/acn3.50972
6. Korley FK, Goldstick J, Mastali M, Van Eyk JE, Barsan W, Meurer WJ, et al. Serum NfL (Neurofilament Light Chain) levels and incident stroke in adults with diabetes mellitus. Stroke. (2019) 50:1669–75. doi: 10.1161/STROKEAHA.119.024941
7. Prins BP, Kuchenbaecker KB, Bao Y, Smart M, Zabaneh D, Fatemifar G, et al. Genome-wide analysis of health-related biomarkers in the UK Household Longitudinal Study reveals novel associations. Sci Rep. (2017) 7:11008. doi: 10.1038/s41598-017-10812-1
8. Teismann H, Wersching H, Nagel M, Arolt V, Heindel W, Baune BT, et al. Establishing the bidirectional relationship between depression and subclinical arteriosclerosis—rationale, design, and characteristics of the BiDirect study. BMC Psychiatry. (2014) 14:174. doi: 10.1186/1471-244X-14-174
9. Schmidt R, Lechner H, Fazekas F, Niederkorn K, Reinhart B, Grieshofer P, et al. Assessment of cerebrovascular risk profiles in healthy persons: definition of research goals and the Austrian stroke prevention study (ASPS). Neuroepidemiology. (1994) 13:308–13. doi: 10.1159/000110396
10. Seiler S, Pirpamer L, Hofer E, Duering M, Jouvent E, Fazekas F, et al. Magnetization transfer ratio relates to cognitive impairment in normal elderly. Front Aging Neurosci. (2014) 6:263. doi: 10.3389/fnagi.2014.00263
11. Hilal S, Tan CS, Adams HHH, Habes M, Mok V, Venketasubramanian N, et al. Enlarged perivascular spaces and cognition: a meta-analysis of 5 population-based studies. Neurology. (2018) 91:e832–42. doi: 10.1212/WNL.0000000000006079
12. Delaneau O, Zagury JF, Marchini J. Improved whole-chromosome phasing for disease and population genetic studies. Nat Methods. (2013) 10:5–6. doi: 10.1038/nmeth.2307
13. Howie BN, Donnelly P, Marchini JA. Flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. (2009) 5:e1000529. doi: 10.1371/journal.pgen.1000529
14. Das S, Forer L, Schönherr S, Sidore C, Locke AE, Kwong A, et al. Next-generation genotype imputation service and methods. Nat Genet. (2016) 48:1284–7. doi: 10.1038/ng.3656
15. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ, et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. (2015) 4:7. doi: 10.1186/s13742-015-0047-8
16. Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. (2017) 8:1826. doi: 10.1038/s41467-017-01261-5
17. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303
18. Wu G, Feng X, Stein L. A human functional protein interaction network and its application to cancer data analysis. Genome Biol. (2010) 11:R53. doi: 10.1186/gb-2010-11-5-r53
19. Yang J, Lee SH, Goddard ME, Visscher PMGCTA. A tool for genome-wide complex trait analysis. Am J Hum Genet. (2011) 88:76–82. doi: 10.1016/j.ajhg.2010.11.011
20. Yang J, Bakshi A, Zhu Z, Hemani G, Vinkhuyzen AAE, Lee SH, et al. Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index. Nat Genet. (2015) 47:1114–20. doi: 10.1038/ng.3390
21. Bulik-Sullivan BK Loh PR Finucane HK Ripke S Yang J Schizophrenia Schizophrenia Working Group of the Psychiatric Genomics Consortium . LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. (2015) 47:291–5. doi: 10.1038/ng.3211
22. Akamine S, Marutani N, Kanayama D, Gotoh S, Maruyama R, Yanagida K, et al. Renal function is associated with blood neurofilament light chain level in older adults. Sci Rep. (2020) 10:20350. doi: 10.1038/s41598-020-76990-7
23. Sugiura N, Adams SM, Corriveau RA. An evolutionarily conserved N-terminal acetyltransferase complex associated with neuronal development. J Biol Chem. (2003) 278:40113–20. doi: 10.1074/jbc.M301218200
24. Ohkawa N, Sugisaki S, Tokunaga E, Fujitani K, Hayasaka T, Setou M, et al. N-acetyltransferase ARD1-NAT1 regulates neuronal dendritic development. Genes Cells. (2008) 13:1171–83. doi: 10.1111/j.1365-2443.2008.01235.x
25. Wang Y, Südhof TC. Genomic definition of RIM proteins: evolutionary amplification of a family of synaptic regulatory proteins. Genomics. (2003) 81:126–37. doi: 10.1016/S0888-7543(02)00024-1
26. Nishiwaki T, Daigo Y, Kawasoe T, Nakamura Y. Isolation and mutational analysis of a novel human cDNA, DEC1 (deleted in esophageal cancer 1), derived from the tumor suppressor locus in 9q32. Genes Chromosomes Cancer. (2000) 27:169–76. doi: 10.1002/(SICI)1098-2264(200002)27:2<169::AID-GCC8>3.0.CO;2-M
27. Ma NX, Yin JC, Chen G. Transcriptome analysis of small molecule-mediated astrocyte-to-neuron reprogramming. Front Cell Dev Biol. (2019) 7:82. doi: 10.3389/fcell.2019.00082
28. Srinivasan M, Mehta P, Yu Y, Prugar E, Koonin EV, Karzai AW, et al. The highly conserved KEOPS/EKC complex is essential for a universal tRNA modification, t6A. EMBO J. (2011) 30:873–81. doi: 10.1038/emboj.2010.343
29. Ramos J, Fu D. The emerging impact of tRNA modifications in the brain and nervous system. Biochim Biophys Acta Gene Regul Mech. (2019) 1862:412–28. doi: 10.1016/j.bbagrm.2018.11.007
30. Arrondel C, Missoury S, Snoek R, Patat J, Menara G, Collinet B, et al. Defects in t6A tRNA modification due to GON7 and YRDC mutations lead to Galloway-Mowat syndrome. Nat Commun. (2019) 10:3967. doi: 10.1038/s41467-019-11951-x
31. Lazarus JE, Moughamian AJ, Tokito MK, Holzbaur EL. Dynactin subunit p150(Glued) is a neuron-specific anti-catastrophe factor. PLoS Biol. (2013) 11:e1001611. doi: 10.1371/journal.pbio.1001611
32. Ayloo S, Lazarus JE, Dodda A, Tokito M, Ostap EM, Holzbaur EL, et al. Dynactin functions as both a dynamic tether and brake during dynein-driven motility. Nat Commun. (2014) 5:4807. doi: 10.1038/ncomms5807
33. Moutaoufik MT, Malty R, Amin S, Zhang Q, Phanse S, Gagarinova A, et al. Rewiring of the human mitochondrial interactome during neuronal reprogramming reveals regulators of the respirasome and neurogenesis. iScience. (2019) 19:1114–32. doi: 10.1016/j.isci.2019.08.057
34. Fink KL, López-Giráldez F, Kim IJ, Strittmatter SM, Cafferty WB. Identification of intrinsic axon growth modulators for intact CNS neurons after injury. Cell Rep. (2017) 18:2687–2701. doi: 10.1016/j.celrep.2017.02.058
35. Smagin DA, Galyamina AG, Kovalenko IL, Babenko VN, Kudryavtseva NN. Aberrant expression of collagen gene family in the brain regions of male mice with behavioral psychopathologies induced by chronic agonistic interactions. Biomed Res Int. (2019) 2019:7276389. doi: 10.1155/2019/7276389
36. Olausson P, Venkitaramani DV, Moran TD, Salter MW, Taylor JR, Lombroso PJ, et al. The tyrosine phosphatase STEP constrains amygdala-dependent memory formation and neuroplasticity. Neuroscience. (2012) 225:1–8. doi: 10.1016/j.neuroscience.2012.07.069
37. Karasawa T, Lombroso PJ. Disruption of striatal-enriched protein tyrosine phosphatase (STEP) function in neuropsychiatric disorders. Neurosci Res. (2014) 89:1–9. doi: 10.1016/j.neures.2014.08.018
38. Jang S-S, Royston SE, Xu J, Cavaretta JP, Vest MO, Lee KY, et al. Regulation of STEP61 and tyrosine-phosphorylation of NMDA and AMPA receptors during homeostatic synaptic plasticity. Mol Brain. (2015) 8:55. doi: 10.1186/s13041-015-0148-4
39. Abu-Rumeileh S, Abdelhak A, Foschi M, D'Anna L, Russo M, Steinacker P, et al. The multifaceted role of neurofilament light chain protein in non-primary neurological diseases. Brain. (2023) 146:421–37. doi: 10.1093/brain/awac328
Keywords: GWAS, neurofilament light chain, serum biomarkers, neuropathology, genetics
Citation: Herrera-Rivero M, Hofer E, Maceski A, Leppert D, Benkert P, Kuhle J, Schmidt R, Khalil M, Wiendl H, Stoll M and Berger K (2023) Evidence of polygenic regulation of the physiological presence of neurofilament light chain in human serum. Front. Neurol. 14:1145737. doi: 10.3389/fneur.2023.1145737
Received: 16 January 2023; Accepted: 20 February 2023;
Published: 08 March 2023.
Edited by:Samir Abu Rumeileh, Martin Luther University of Halle-Wittenberg, Germany
Reviewed by:Tommaso Schirinzi, University of Rome Tor Vergata, Italy
Lorenzo Barba, University Hospital in Halle, Germany
Copyright © 2023 Herrera-Rivero, Hofer, Maceski, Leppert, Benkert, Kuhle, Schmidt, Khalil, Wiendl, Stoll and Berger. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Marisol Herrera-Rivero, email@example.com
†Present address: Marisol Herrera-Rivero, Department of Psychiatry, University of Münster, Münster, Germany