Association Study of a Comprehensive Panel of Neuropeptide-Related Polymorphisms Suggest Potential Roles in Verbal Learning and Memory

Neuropeptides are mostly expressed in regions of the brain responsible for learning and memory and are centrally involved in cognitive pathways. The majority of neuropeptide research has been performed in animal models; with acknowledged differences between species, more research into the role of neuropeptides in humans is necessary to understand their contribution to higher cognitive function. In this study, we investigated the influence of genetic polymorphisms in neuropeptide genes on verbal learning and memory. Variants in genes encoding neuropeptides and neuropeptide receptors were tested for association with learning and memory measures using the Hopkins Verbal Learning Test—Revised (HVLT-R) in a healthy cohort of individuals (n = 597). The HVLT-R is a widely used task for verbal learning and memory assessment and provides five sub-scores: recall, delay, learning, retention, and discrimination. To determine the effect of candidate variants on learning and memory performance, genetic association analyses were performed for each HVLT-R sub-score with over 1300 genetic variants from 124 neuropeptide and neuropeptide receptor genes, genotyped on Illumina OmniExpress BeadChip arrays. This targeted analysis revealed numerous suggestive associations between HVLT-R test scores and neuropeptide and neuropeptide receptor gene variants; candidates include the SCG5, IGFR1, GALR1, OXTR, CCK, and VIPR1 genes. Further characterization of these genes and their variants will improve our understanding of the genetic contribution to learning and memory and provide insight into the importance of the neuropeptide network in humans.


Introduction
Cognition and behaviour are the foremost important outputs of human brain functioning and are underpinned by memory functions mediated by a set of neural mechanisms, which allow information to be encoded, consolidated, retained, and retrieved.Memory provides us with a sense of self that links our past, present, and future and is important for individual, social, educational, and vocational functioning as exemplified by the problems experienced by people with memory impairment or amnesic syndromes due to brain injuries [1].
Different types of memory are mediated by different parts of the brain, although they may interact and partly overlap [2].Implicit memory refers to the unconscious facilitation of recall or recognition by previously encountered information and relies on the basal ganglia and cerebellum, while explicit memory refers to the ability to recall or recognise information consciously and is associated with the hippocampus, the neocortex, and the amygdala.Most developed memory tests measure explicit memory, which is subdivided into semantic (knowledge of facts) and episodic (recall or recognition of past events and experiences) [3].Furthermore, there are three subtypes of episodic memory: short-term memory (STM), long-term memory (LTM), and working memory (WM).STM involves processes that last for only a few seconds or minutes and involves hippocampal regions, while LTM entails processes that involve gating of memories formed in the hippocampus through a thalamocortical circuit for longer-term cortical storage [4].WM refers to the ability to hold and actively process information immediately after its presentation and is closely linked to the prefrontal cortex (PFC).Verbal and visuospatial information is also processed in separate domains; verbal tasks have been shown to preferentially involve activity in a network of regions including the PFC and along the inferior prefrontal gyrus, while activity in the occipital-temporal cortex correlates with processing of visual information [5].
Cognitive and behavioural functions are closely entwined and are regulated with overlapping neuronal networks and dependent on chemical neurotransmission.Neurotransmitters can be classified into three categories: fast signalling (including the excitatory amino acids and catecholamines), slow modulatory, and neuropeptides [6].Neuropeptides are 3-to-100-amino-acid-long hormone-like signalling molecules produced by neurons and some other tissues and are important for communication between neurons and other neural cell types as well as non-neuronal target cells [7].It is important to note that the number of known neuropeptides far exceeds that of classical neurotransmitters.Due to the slow and more diffuse release of neuropeptides, they are usually considered to be second messengers and have neurotrophic and neuromodulatory effects.Previous studies have shown the involvement of neuropeptides in numerous neuroendocrine and behavioural processes including metabolism, stress control, learning, and memory [8][9][10].
Neuropeptides localize in regions involved in cognition, learning, and memory in the brain and have been characterized as regulators of these pathways [11,12].To date, research into the function of neuropeptides in learning and memory has mostly been conducted in animal models, with several neuropeptide and receptor families, including the somatostatins, tachykinin, galanin, and vasoactive intestinal peptide (VIP), linked with cognition [12].For example, VIP and its receptors have been implicated in learning and memory processes in mice through investigating gene knockouts [13] and pharmacological interventions which demonstrate a role in synaptic transmission [14].Although there is conservation of the neuropeptide network across mammals [15], studying genetic variation in neuropeptide and neuropeptide receptor genes may allow further dissection and characterization of their roles in human cognitive and behavioural functions.Twin studies have demonstrated that specific cognitive and memory traits are genetically influenced, with heritability estimates ranging from 30 to 60% [16][17][18].Genetic polymorphisms, including single-nucleotide polymorphisms (SNPs) are used in association studies to identify genetic loci that influence human traits and disorders, including complex polygenic cognitive and memory phenotypes [18].
Neuropeptide and neuropeptide receptor gene polymorphisms have been investigated in various human behavioural disorders (e.g., sleep, feeding, and obesity) [19] as well as social and neuropsychological traits [20].There has been much interest in the oxytocin/vasopressin neuropeptide family and oxytocin receptor gene (OXTR); e.g., the common variant rs53576 in OXTR has been investigated in personality traits and behaviour [21][22][23][24].In 2015, Li et al. conducted a meta-analysis of OXTR rs53576 and reported a positive association between the SNP and general sociality but not close relationships (parent-child and/or romantic relationship) or depression [25].SNPs in OXTR have also been studied in relation to autism and intelligence with SNPs and haplotypes in the OXTR gene identified to confer risk for autism spectrum disorders [26].In animal studies, oxytocin and vasopressin have been shown to be associated with learning and memory [27], and in humans, shorter microsatellite repeat polymorphisms in the arginine vasopressin receptor 1A (AVPR1A) gene promoter have been associated with poorer verbal memory performance via modulation of hippocampal structures and connectivity with the thalamus [28].Neuropeptide Y (NPY) Y 2 receptor engagement was shown to be required for learning and memory processing in a mouse knockout model [29].In humans, the NPY -485T promoter allele was shown to be associated with schizophrenia susceptibility due to a decreased level of neuropeptide in the brain [30]; it is worth noting that memory defects are common in schizophrenia.
Although several studies have investigated specific genetic polymorphisms in neuropeptide and neuropeptide receptor genes in relation to social behaviour, learning, and memory, comprehensive studies of neuropeptide-related genetic variants have not been conducted.Here, our aim was to explore associations between an extensive list of neuropeptide and neuropeptide receptor gene variants and verbal learning and memory assessed using the Hopkins Verbal Learning Test-Revised (HVLT-R) in a healthy Australian cohort.The HVLT-R is a widely accepted neurocognitive assessment of verbal learning and memory with measures for five sub-components: recall, delay, learning, retention, and recognition discrimination [31].In this study, we used neuropeptide databases, http://www.neuropeptides.nl/(accessed on 10 January 2017) [7] and NeuroPedia [32], to compile a list of more than 100 human genes encoding classical and candidate neuropeptide genes and their receptors.From these, we identified >1300 SNPs present on Illumina ® Human OmniExpress-24 BeadChip arrays (Illumina, Inc., San Diego, CA, USA) and tested for genetic associations with the measures of verbal learning and memory assessed by the HVLT-R sub-components in a healthy population cohort.

Participants
Participants in this study were individuals recruited through advertisements from the Brisbane and Gold Coast areas of South-East Queensland in Australia.Participation was excluded for individuals with a history of psychiatric disorder or head injury to maintain a representative sample of cognitive and memory ability without additional complications, although it should be noted that other neurological disorders that may impact some memory functions, e.g., autism spectrum disorder, were not specifically screened for.Written informed consent was provided by all participants prior to any study activities.The study was approved by the Griffith University (MSC/01/09/HREC) and Queensland University of Technology Human Research and Ethics (1300000486) Committees.

Phenotyping
All participants were assessed individually and under the same conditions, in a quiet and well-lit room, by the same examiner.

Learning and Memory
Verbal learning and memory abilities of the participants were determined using the HVLT-R test [31].The test is based on three trials of free recall of 12 items followed by a 20 to 30 min delayed free recall and subsequent yes and no recognition from a semantically categorized set of words (n = 24).The HVLT-R evaluates learning and memory using 5 sub-components: recall, delay, learning, retention, and recognition discrimination, which are calculated as follows.
The intelligence quotient (IQ) was measured using subsets of the Wechsler Abbreviated Scale of Intelligence (WASI) IQ test.The WASI-IQ test was created in 1955 by David Wechsler and is a well-established IQ test for measuring adult intelligence [33].The vocabulary and matrix reasoning subsets of WASI were completed in this study to estimate the IQ of participants and was used as one of the covariates in the analysis.

Sample Collection and DNA Extraction
Saliva samples were collected from the participants immediately after completion of the memory tests using Oragene ® DNA Self Collection kits (DNA Genotek Inc., Ottawa, ON, Canada).DNA was extracted from whole saliva samples by using the kit and protocol of the same manufacturer.

Genotyping and Gene Selection
Saliva DNA samples were genotyped on Illumina ® Human OmniExpress-24 BeadChip arrays (Illumina, Inc., San Diego, CA, USA) and genotype data for 597 individuals were available after quality control (QC).Genes encoding neuropeptides were identified using two neuropeptide databases.A total of 93 genes were obtained using NeuroPedia [32] and nine additional genes incorporated from the www.neuropeptides.nldatabase (accessed on 10 January 2017) [7].Furthermore, 53 neuropeptide receptor genes from the Illumina OmniExpress 24 v1.1 annotation file were included in the gene set (Supplementary Table S1).From this gene set, we compiled a list of SNPs from the exonic and intronic regions of the loci of interest and extracted the genotype data for the memory cohort.We failed to locate SNPs from 19 of the neuropeptide genes that had been identified from the databases in our genotype data and SNPs from 12 genes did not pass genotyping quality thresholds.Therefore, the final SNP list comprised 1306 SNPs located in 124 neuropeptide and receptor genes which are shown in Table 1.

Statistical Analysis
The R Program for Statistical Computing (v3.3.3)[34] was used to complete the descriptive statistics of the learning and memory phenotypes.Verbal learning and memory scores and IQ were investigated for their correlation using Pearson's r test.Population structure was studied using The R Program and KING (v1.9), a population structure inference tool [35], and principal component analysis (PCA) was undertaken to consider population structure.Self-reported ancestry did not cluster together, indicating an unrecognised mixed background for various individuals, and supported the use of principal components (PCs) instead of self-reported ethnicity in analyses.
QC of the genotype data and association analyses were conducted using PLINK (v1.09) [36].QC thresholds for the analysis were determined as follows: minor allele frequency (MAF) and Hardy Weinberg equilibrium (HWE) set to higher than 0.01 and 0.001, respectively.QC was completed for the selected set of markers and SNPs located in the selected genes that passed thresholds were included in the study.Generalized linear model analysis was carried out for each memory phenotype individually to identify correlations.Age, sex, IQ, PC1, and PC2 were considered as covariates in the analysis.Association analyses were performed for the SNPs within genes of interest from the Illumina OmniExpress BeadChip assays to estimate genotypic effects on learning and memory status.To eliminate type I error due to multiple testing, we calculated the suggestive and significance p-value thresholds using the Genetic type 1 Error Calculator (GEC) (v0.2) [37].A suggestive p-value threshold was set to 0.00113, while the significance p-value was 0.0000564.

Demographics and Phenotype Analysis
A cohort of healthy participants recruited from the Brisbane and Gold Coast areas of Australia undertook a battery of memory tests, including the HVLT-R to assess verbal learning and memory, and donated a saliva sample to extract DNA for genetic analysis.Statistical analysis was performed for the 597 participants for whom genotyping data were obtained and the demographics of this memory cohort are presented in Table 2. Two-thirds of the cohort were female (71%) with an age range of 16 to 65 years old (M = 20, SD = 8.57).The majority of the cohort identified as Caucasian (75%).Due to the large number of individuals with other ethnicities (n = 151, mostly Asian), population structure was studied using PCA, and the first two PCs were added into analyses as a covariate, along with age and gender.IQ was normally distributed in our cohort and ranged from 78 to 138 (µ = 108.1).In this study, verbal learning and memory performance was measured using the HVLT-R using five sub-components: recall, delay, learning, retention, and discrimination.When we further examined these sub-components, we found correlations between recall HVLT-R and delay HVLT-R measures (0.64) as well as between delay HVLT-R and retention HVLT-R measures (0.72).Analysis of HVLT-R subtest scores (recall, delay, learning, retention, discrimination) by age using Pearson's correlation showed a moderate positive correlation between score and age (r = 0.2, p < 0.001), with no other significant correlations observed.Analysis by sex using an independent T-test showed that, in this cohort, males had higher average IQ scores compared to females (p < 0.001), while comparisons for HVLT-R scores were all non-significant.Finally, analysis by ethnic group (Caucasian versus non-Caucasian) using an independent-samples T-test showed higher average recall HVLT-R , delay HVLT-R , and IQ scores, compared to non-Caucasians (p = 0.002, p = 0.005, and p < 0.001, respectively), while other comparisons were non-significant.Significant results may reflect the relatively small sample sizes for some of the groups and could suggest some impact of English as a second language on some of the verbal memory scores.The results justify adding age, sex, IQ, and PCs 1 and 2 (accounting for population structure) as covariates while subsequent analyses were performed on the entire cohort to maintain statistical power.

Genotype Association
We compiled a comprehensive list of 80 human neuropeptide and 44 neuropeptide receptor genes from databases (Table 1) and identified SNPs in their exonic and intronic regions that were present on Illumina ® Human OmniExpress-24 BeadChip arrays.We extracted genotypes for these SNPs from the microarray data of the 597 individuals of the memory cohort.Due to the selected quality control thresholds (minor allele frequency (MAF) > 0.01; Hardy Weinberg equilibrium (HWE) > 0.001), eight neuropeptide genes and four neuropeptide receptors were excluded from the analysis.Association analyses with the test scores from each of the HVLT-R sub-components were performed for the 1306 remaining SNPs using a generalized linear model (GLM).As our memory cohort comprised mostly females (71%), young adults (78% were 16 to 25 years old), and Caucasians (75%) (Table 2), we included age, gender, and population structure (PC1 and PC2) as well as IQ to the analysis as covariates in order to minimize their effect on the study and to focus on the association of the loci of interest and learning and memory.To consider and eliminate type I error due to multiple testing, the p-value thresholds were calculated and adjusted to be able to declare the significant SNPs.To calculate our adjusted p-value thresholds, SNPs in linkage disequilibrium (LD) were considered and after estimating the effective number of independent markers (n = 887), we set the significance threshold to 0.0000564 and a suggestive p-value as 0.00113.
While 298 markers were below an α level of 0.05 in one phenotype or more, none of the markers in the loci of interest were significantly associated with verbal learning and memory after the multiple testing corrections (Supplementary Table S2).However, nine SNPs were found to be suggestively associated with HVLT-R phenotypes, with rs1997317 in the Secretogranin V gene (SCG5) approaching significance with the delay HVLT-R component (Table 3).Box plots showing the allelic distributions of the suggestive SNPs presented in Table 3 are illustrated in Figure 1.We then used the GTEx (Genotype-Tissue Expression) expression quantitative trait loci (eQTL) Browser (http://www.gtexportal.org/home/)to further investigate eQTLs of the SNPs potentially linked to learning and memory traits in Table 3. rs11634086 is a significant eQTL for SCG5 in testis (p = 1.47 × 10 −6 ), but it is also an eQTL in various human  We then used the GTEx (Genotype-Tissue Expression) expression quantitative trait loci (eQTL) Browser (http://www.gtexportal.org/home/)to further investigate eQTLs of the SNPs potentially linked to learning and memory traits in Table 3. rs11634086 is a significant eQTL for SCG5 in testis (p = 1.47 × 10 −6 ), but it is also an eQTL in various human brain regions including the cortex (p = 6.5 × 10 −5 ).In addition to the pancreas and pituitary, SCG5 shows high expression levels in brain-related tissues (Figure 2), which may support the supposition that the SNPs we have identified may impact or are in high LD with SNPs that affect memory performance.rs2684788 is an eQTL for IGF1R in brain tissues, including the cerebellum (p = 5.1 × 10 −15 ) and hippocampus (2.6 × 10 −4 ), and shows expression in these regions as well as high expression in arteries and thyroid tissue.GALR1 has higher expression only in the pitui- rs2684788 is an eQTL for IGF1R in brain tissues, including the cerebellum (p = 5.1 × 10 −15 ) and hippocampus (2.6 × 10 −4 ), and shows expression in these regions as well as high expression in arteries and thyroid tissue.GALR1 has higher expression only in the pituitary according to GTEx.The SNP rs8192472 is reported as an eQTL in GTEx for CCK in the traverse colon (p = 4.3 × 10 −5 ), with CCK being highly expressed in some tissue regions of the brain (Figure 3).rs437876 is an eQTL for VIPR1 in whole blood (1.4 × 10 −7 ), with some signal in brain regions such as the cortex (p = 2.9 × 10 −3 ) and amygdala (p = 0.03).The other suggestively associated SNPs in Table 3 were not present in GTEx.Both VIPR1 and OXTR are expressed in some brain regions, including areas important for cognition, but VIPR1 is most highly expressed in the lungs and OXTR in breast tissue according to GTEx.

Discussion
In the present study, we performed targeted analysis of SNPs in neuropeptide-and neuropeptide receptor-encoding genes, to test for associations between genotypes and verbal learning and memory performances as assessed by the HVLT-R, in a healthy Australian cohort.We found numerous SNPs in various neuropeptide and neuropeptide receptor genes with a p-value < 0.05 (Supplementary Table S2).However, after multiple-testing correction and adjustment to the significance threshold (0.0000564), no SNPs were identified as significant.Nevertheless, suggestive associations (below the suggestive p-value of 0.00113) remained for nine SNPs in six neuropeptide and neuropeptide receptor genes (SCG5, IGF1R, GALR1, OXTR, CCK, VIPR1) with specific HVLT-R components of recall, delay, learning, retention, and discrimination (Table 3).
Three suggestively associated SNPs were identified in the SCG5 gene that encodes Secretogranin V, belonging to the chromogranin and secretogranin family of proteins; rs11634086 and rs10519737 (both located in intron 2) were associated with the recall HVLT-R score and rs1997317 (located in the final intron) was associated with the delay HVLT-R score.Minor alleles of these SNPs were positively associated with higher memory performance in the recall HVLT-R and delay HVLT-R sub-components.The recall and delay sub-components

Discussion
In the present study, we performed targeted analysis of SNPs in neuropeptide-and neuropeptide receptor-encoding genes, to test for associations between genotypes and verbal learning and memory performances as assessed by the HVLT-R, in a healthy Australian cohort.We found numerous SNPs in various neuropeptide and neuropeptide receptor genes with a p-value < 0.05 (Supplementary Table S2).However, after multiple-testing correction and adjustment to the significance threshold (0.0000564), no SNPs were identified as significant.Nevertheless, suggestive associations (below the suggestive p-value of 0.00113) remained for nine SNPs in six neuropeptide and neuropeptide receptor genes (SCG5, IGF1R, GALR1, OXTR, CCK, VIPR1) with specific HVLT-R components of recall, delay, learning, retention, and discrimination (Table 3).
Three suggestively associated SNPs were identified in the SCG5 gene that encodes Secretogranin V, belonging to the chromogranin and secretogranin family of proteins; rs11634086 and rs10519737 (both located in intron 2) were associated with the recall HVLT-R score and rs1997317 (located in the final intron) was associated with the delay HVLT-R score.Minor alleles of these SNPs were positively associated with higher memory performance in the recall HVLT-R and delay HVLT-R sub-components.The recall and delay sub-components are the only memory phenotypes assessed by the test, with recall HVLT-R assessing immediate verbal memory, and longer-term verbal memory assessed by the delay HVLT-R , using the same word list but measured after a 20 to 30 min interval.This study presents the first suggestion that polymorphisms in SCG5 may be associated with human memory and is supported by in silico analysis of SCG5 expression showing high levels in the brain and its sub-regions important in learning and memory (Figure 2).While these SNPs have not been previously associated with any phenotypes, Cao-Lei et al. previously reported a behavioural association for SCG5, showing that DNA methylation changes in SCG5 are highly correlated with maternal objective stress [38].
The insulin-like growth factor 1 receptor (IGFR1) is a member of the insulin gene family.Two SNPs at this gene locus, rs867431 and rs2684788, showed associations with the learning HVLT-R and retention HVLT-R scores, respectively.These SNPs have not been previously linked to human behaviour or learning and memory.rs2684788 is a 3 UTR SNP and was previously reported to be associated with susceptibility to idiopathic short stature [39].The other SNP that showed correlation with the retention HVLT-R score was the galanin receptor 1 (GALR1) polymorphism rs9959924.GALR1 belongs to the galanin neuropeptide family; galanin has been studied in Alzheimer's disease, epilepsy, and memory consolidation studies [40][41][42][43] and is known for its potential role in learning and memory [11].The GALR1 SNP rs9959924 we identified to be potentially associated with verbal retention was previously studied in relation to smoking and nicotine dependence and found to be significantly associated with smoking quantity [44].While not directly related to memory, it suggests the SNP may have a functional effect.
SNPs rs237888 (OXTR), rs8192472 (CCK), and rs437876 (VIPR1) showed suggestive associations with the descrimination HVLT-R score, an index obtained by subtracting the number of objects recognized falsely from a list from the number of true positives.OXTR is an oxytocin receptor gene in which a number of polymorphisms have been widely studied in autism, cognition, and behaviour [26,[45][46][47].While we did not find high LD between rs237888 and these other OXTR SNPs reported in the literature, however, OXTR SNP rs237888 has previously been found to be associated with recognition and decision making [48], consistent with our findings.The cholecystokinin gene (CCK) has been studied in relation to anxiety and cognition [49] and more recently has been associated with body mass index, similar to other neuropeptides that have been found to regulate eating behaviour and energy homeostasis [50].Furthermore, in murine models, CCK has been found to modulate hippocampal long-term potentiation (LTP) and spatial learning and memory by facilitating neuroplasticity of hippocampal CA3-CA1 synapses [51].
VIPR1, vasoactive intestinal polypeptide receptor 1 gene, encodes the receptor for VIP, which is known to play a modulatory role in learning and memory [52].VIP is a neuropeptide released from dentate gyrus interneurons and enhances the proliferative and proneurogenic effect of microglia via IL-4 release by the VPAC1 receptor [53].Furthermore, Cunha-Reis et al. studied long-term depression, focusing on VIPR1 for possible pharmacological targets for treatment of cognitive dysfunction [54].This supports our finding that VIPR1 polymorphisms may play a role in learning and memory.Although VIPR1 shows only moderate expression in tissues related to the brain or cognition in GTEx, VIP and its receptor gene VIPR1 are expressed in the postnatal dentate gyrus where they direct neural stem/precursor cells towards a neuronal fate [52].
There are a number of limitations in our study.Although we identified suggestive associations for neuropeptide-related genes in learning and memory phenotypes, none were significant.The relatively small sample size may have contributed to SNPs failing to pass the significance threshold.Furthermore, we did not impute our genotype data due to the considerable proportion of non-Caucasian individuals in the cohort; imputation may have increased power to detect possible associated variants and would have provided a higher resolution of association signals.Another limitation is that the SNPs were not validated in an independent cohort.Since the HVLT-R is a widely used assessment of learning and memory, it would be of interest to investigate these SNPs and genes in other cohorts phenotyped using either this specific test or other verbal learning and memory measures.We also only included SNPs that were in the coding and intronic regions of the neuropeptide and neuropeptide receptor genes in the study.SNPs in the upstream or downstream regulatory regions of these genes may also affect their expression levels and may, therefore, also be relevant for future studies.Finally, this was a targeted study and genome-wide approaches are necessary for more comprehensive association testing, although such approaches require larger sample sizes, particularly if cohorts or phenotyping is heterogeneous.Notably, the well-known Alzheimer's disease APOE/APOC1/TOMM40 region has been associated with verbal memory phenotypes in a number of studies [55][56][57], and a large meta-analysis additionally identified significant signals at CDH18, NT5DC2, STAB1, ITIH1, ITIH4, and PBRM1 and implicated synaptic and neurodevelopmental genes [57].

Conclusions
Our study is the first to examine verbal learning and memory in a healthy population assessed using the HVLT-R with respect to a comprehensive list of neuropeptide and neuropeptide receptor gene variants.In summary, we found a number of suggestive associations between SNPs in genes related to neuropeptide networks with verbal learning and memory.Our findings establish new potential targets for future studies and enhance our knowledge on the importance of neuropeptides in learning and memory.Interestingly, similar to work in animal studies, treatment of patients with memory deficits (e.g., with schizophrenia) with exogenous neuropeptides such as arginine vasopressin has been shown to have potential therapeutic benefits for cognitive and memory functions [58].Further work investigating the role of neuropeptide and neuropeptide receptor genes in learning and memory is warranted.

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes15010030/s1,Table S1: SNP list obtained from neuropeptide and neuropeptide receptor genes.;Table S2: Suggestively significant association results for neuropeptide and neuropeptide gene receptor SNPs with HVLT-R subcomponent performance phenotypes.
Author Contributions: N.A. analysed and interpreted the data, generated figures, and co-wrote the paper; H.G.S. Sutherland prepared the samples for genotyping, provided supervision, and co-wrote the paper; R.A.L. was involved in study design and funding acquisition, performed analyses, and supervised the study; L.M.H. was involved in supervision and editing of the manuscript; D.H.K.S. was involved in funding acquisition, designed the memory performance evaluations, oversaw collection of the phenotype data, and provided supervision for the study; L.R.G. was involved in design, funding acquisition, and overall supervision of the study.All authors contributed to the editing of the final manuscript.All authors have read and agreed to the published version of the manuscript.Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

17 Figure 1 .
Figure 1.Hopkins Verbal Learning Test-Revised (HVLT-R) phenotypes and allelic distributions of suggestively associated SNPs from neuropeptide and neuropeptide receptor genes.Genotypes homozygous for the minor allele are bolded.Circles represent values of outlier samples.

Figure 1 .
Figure 1.Hopkins Verbal Learning Test-Revised (HVLT-R) phenotypes and allelic distributions of suggestively associated SNPs from neuropeptide and neuropeptide receptor genes.Genotypes homozygous for the minor allele are bolded.Circles represent values of outlier samples.

Figure 3 .
Figure 3. GTEx CCK gene expression illustrated across brain regions.(RPKM: reads per kilo base transcript per million).Circles represent values of outlier samples.

Funding:
This work was supported by an Australian Research Council grant-The Genetic Basis of Human Memory (Project ID: DP130101921).This work was also supported by infrastructure purchased with Australian Government EIF Super Science Funds as part of the Therapeutic Innovation Australia-Queensland Node project.Institutional Review Board Statement: The study was conducted in accordance with the Declaration of Helsinki, and approved by the Griffith University Human Research and Ethics Committee (MSC/01/09/HREC [2009]) and Queensland University of Technology Human Research and Ethics Committee (1300000486 [2013] and 2000000491 [22/09/2020]).

Table 1 .
Overview of neuropeptide-encoding genes and neuropeptide receptor genes that were analysed in this study.

Table 2 .
Demographics of the memory cohort.

Table 3 .
Suggestive genetic associations for nine SNPs in neuropeptide-related genes with HVLT-R subcomponent performance phenotypes.