Xiao-Lin Yang1,2†
Shao-Yan Zhang2,3†Hong Zhang1,2Xin-Tong Wei1,2Gui-Juan Feng1,2Yu-Fang Pei2,3*‡
Lei Zhang1,2*‡- 1Center for Genetic Epidemiology and Genomics, School of Public Health, Medical College, Soochow University, Jiangsu, China
- 2Jiangsu Key Laboratory of Preventive and Translational Medicine for Geriatric Diseases, School of Public Health, Medical College, Soochow University, Jiangsu, China
- 3Department of Epidemiology and Health Statistics, School of Public Health, Medical College, Soochow University, Jiangsu, China
As an important trait at birth, infant head circumference (HC) is associated with a variety of intelligence- and mental-related conditions. Despite being dominated by genetics, the mechanism underlying the variation of HC is poorly understood. Aiming to uncover the genetic basis of HC, we performed a genome-wide joint association analysis by integrating the genome-wide association summary statistics of HC with that of its two related traits, birth length and birth weight, using a recently developed integrative method, multitrait analysis of genome-wide association (MTAG), and performed in silico replication in an independent sample of intracranial volume (N = 26,577). We then conducted a series of bioinformatic investigations on the identified loci. Combining the evidence from both the MTAG analysis and the in silico replication, we identified three novel loci at the genome-wide significance level (α = 5.0 × 10−8): 3q23 [lead single nucleotide polymorphism (SNP) rs9846396, pMTAG = 3.35 × 10−8, preplication = 0.01], 7p15.3 (rs12534093, pMTAG = 2.00 × 10−8, preplication = 0.004), and 9q33.3 (rs7048271 pMTAG = 9.23 × 10−10, preplication = 1.14 × 10−4). Each of the three lead SNPs was associated with at least one of eight brain-related traits including intelligence and educational attainment. Credible risk variants, defined as those SNPs located within 500 kb of the lead SNP and with p values within two orders of magnitude of the lead SNP, were enriched in DNase I hypersensitive site region in brain. Nine candidate genes were prioritized at the three novel loci using multiple sources of information. Gene set enrichment analysis identified one associated pathway GO:0048009, which participates in the development of nervous system. Our findings provide useful insights into the genetic basis of HC and the relationship between brain growth and mental health.
Introduction
The size of infant head circumference (HC) is thought to affect brain development in later life (Ivanovic et al., 2004). It is associated with a variety of intelligence- and mental-related conditions, such as intelligence quotient (IQ) score (Ivanovic et al., 2004), autism (Aylward et al., 2002), and Alzheimer’s disease (Schofield et al., 1997). Its relation to nonmental-related disorders, such as cardiovascular disease, is also reported (Barker et al., 1993). Therefore, understanding the mechanism of HC development is important in relation to its impact on phenotypic outcomes in later life in order to better prevent and intervene the related disorders.
Despite being influenced by maternal nutrition during gestation, infant HC is dominantly determined by inheritance (Smit et al., 2010). A previous twin study estimated its heritability to be as high as 0.7–0.9 (Smit et al., 2010). Early candidate gene studies have identified multiple rare genomic mutations with large effects (Bond et al., 2005; Kumar et al., 2009). For common variants, Taal et al. (2012) and Haworth et al. (2019)performed two genome-wide association study (GWAS) meta-analyses in European subjects and identified five loci (4q28.1, 6p21.21, 12q15, 12q24, and 17p13.1) in total. Nonetheless, variants identified by these studies explain only a small proportion of the total HC heritability, and the vast majority of missing heritability is yet to be discovered.
As a measure at birth, HC is related to birth length (BL) and birth weight (BW). They have robust links with each other in epidemiology and are all associated with common factors such as fetal sex, maternal factors, and gestational age (Yokoyama et al., 2005; Rydbeck et al., 2014). Many studies attempt to define some conceptual growth parameters by integrating the three traits to comprehensively reflect children’s health status (Sankilampi et al., 2013; Haksari et al., 2016). Their high phenotypic correlation (Scheffler et al., 2017) implies that these traits may share some common genetic architecture.
Recently, a new statistical method multitrait analysis of genome-wide association (MTAG) was proposed, which integrates the information contained in correlated traits to boost genetic association signals for each individual trait (Turley et al., 2018). In the present study, we will perform a joint genetic association analysis of HC, BL, and BW with MTAG. The main purpose is to identify novel genomic loci for HC by integrating the GWAS summary results from BL and BW. We will perform independent in silico replication in large-scale GWAS summary results of adult intracranial volume. We will then conduct genomic feature enrichment analysis and candidate gene prioritization. We will also perform gene set enrichment analysis to highlight relevant pathways.
Materials and Methods
Study Design
In this study, we performed a joint GWAS analysis by integrating the information carried in the GWAS summary results of HC, BL, and BW traits, with the main purpose of identifying responsible genomic loci for HC.
Gwas Summary Results
All GWAS meta-analyses for the three individual traits were carried out by the Early Growth Genetics (EGG) consortium, and basic characteristics of individual studies are detailed in respective publications (Table 1) (Taal et al., 2012; van der Valk et al., 2015; Horikoshi et al., 2016). In brief, these three studies conducted meta-analyses in 28,459, 10,678, and 143,677 European participants for BL, HC, and BW, respectively. We accessed GWAS summary results from the consortium’s website (http://egg-consortium.org) and retrieved the part of results that included participants with European ancestry only for each trait.
Table 1 Basic characteristics of the three studied traits.
Single Nucleotide Polymorphism Inclusion Criteria
A series of sequential quality control (QC) criteria were applied to retain eligible single nucleotide polymorphisms (SNPs): 1) Duplicated SNPs were dropped within each trait (0 for BL, 2 for HC, 30,565 for BW); 2) SNP genomic coordinate for BL and HC was transformed from GRCH 36 to GRCh 37 with LiftOver (http://genome.ucsc.edu/cgi-bin/hgLiftOver), and deprecated SNPs were dropped (213 for BL and 256 for HC); 3) for BL summary results with no allele frequency information, allele frequency was estimated based on the HapMap Phase 2 + 3 European individual genotypes, and SNPs with missing allele frequency were excluded (Ivanovic et al., 2004, 748); 4) rare SNPs of minor allele frequency (MAF) below 1% were excluded within each trait (2,876 for BL and 6,665,395 for BW); 5) SNPs on which the sample size was smaller than 75% of 90th percentile of SNP-level sample size distribution were excluded within each trait (3,809 for BL, 36,610 for HC, and 10,607 for BW) (Turley et al., 2018); and 6) SNPs having nonconcordant alleles (e.g., A/G vs. A/C, 61) or with ambiguous strands (e.g., A/T and C/G, 168,789) across the three traits were excluded. After the above QC steps, there were 1,977,400 qualified SNPs in common to the three traits. All subsequent analyses were based on these SNPs.
LDSC Analysis
We applied the linkage disequilibrium score (LDSC)-based regression analysis to the summary results of each trait to evaluate its genetic architecture. LDSC takes GWAS summary statistics as input and partitions overall inflated association statistics into one part attributable to polygenic architecture and another part due to cryptic relatedness and population stratification (Bulik-Sullivan et al., 2015b). LDSC can also estimate genetic heritability attributable to the GWAS SNPs. Reference LD scores for the European population were downloaded from the software website (https://github.com/bulik/ldsc).
We also applied bivariate LDSC regression analysis implemented in LD-Hub (http://ldsc.broadinstitute.org/ldhub/) to the GWAS summary results of HC and BL traits, and of HC and BW traits, to study the genetic correlations between HC and the other two traits (Bulik-Sullivan et al., 2015a). In brief, this method is an extension of single-trait LD score regression and requires only GWAS summary statistics from individual studies, being robust to sample overlap.
Joint Analysis by MTAG
We applied the MTAG method to the three traits for a joint association analysis (Turley et al., 2018). MTAG estimates per SNP effect size for each trait by incorporating information contained in the other correlated traits and therefore has potential to improve statistical power of association test. MTAG takes summary statistics from multiple studies as input. We used the 1,000 genomes project European individuals as reference panel for LD estimation for all the three traits. The effect size re-estimated by MTAG is a generalized estimate of inverse-variance-weighted meta-analysis by integrating association summary statistics from different traits, where p value is derived from the re-estimated effect size (Turley et al., 2018). A distinct locus was defined as a genomic region encompassing the lead SNP and its flanking 500-kb region in either direction.
Variance Explained by SNPs
All the three traits are standardized (i.e., mean 0 and variance 1) in their original studies. Therefore, we estimated individual variant effect size as the variance it explains with the formula 2f(1 − f)β2, where f is allele frequency and β is regression coefficient.
In Silico Replication in Intracranial Volume Summary Results
Each lead SNP from the newly identified loci was subject to independent in silico replication in intracranial volume. This is because HC is both phenotypically and genetically highly correlated with intracranial volume (Adams et al., 2016; Breakey et al., 2018). Their genetic correlation coefficient is reported to be as high as 0.75, indicating a shared genetic background (Adams et al., 2016). Some previous studies integrate these two traits together for a direct meta-analysis (Adams et al., 2016; Haworth et al., 2019).
The summary results of intracranial volume come from a large-scale GWAS meta-analysis conducted by the CHARGE and the ENIGMA consortia in 26,577 subjects of European ancestry and were downloaded through the ENIGMA’s website (http://enigma.ini.usc.edu/research/download-enigma-gwas-results/) with permission.
The effect direction was first compared between the discovery and the replication stages. Given consistent effect direction, one-sided p value from replication stage was then reported. To account for multiple testing problem, the Bonferroni correction was adopted to declare statistical significance in the replication stage.
In Silico Replication of BL and BW
Although the focus of the present study was HC, it did produce re-estimated MTAG results for BL and BW. We performed in silico replication of these two traits to strengthen the validity of our analytical approach. For BL, we replicated in the large-scale GWAS summary results released by the GIANT consortium for adult height (Yengo et al., 2018) (http://portals.broadinstitute.org/collaboration/giant/index.php/GIANT_consortium_data_files), where adult height is believed to closely link to the BL (van der Valk et al., 2015). For BW, we replicated in BW GWAS summary results conducted in ∼330,000 UK Biobank participants. The analysis was performed by Dr. Neale’s lab, and the results were kindly released to be publicly available at the lab website (http://www.nealelab.is/uk-biobank). We kept in mind that about 67,000 participants overlapped between this analysis and the BW results used in the present study. Again, the replication significance was evaluated in terms of both effect direction and Bonferroni corrected one-sided p value.
Association With Brain/Mental-Related Traits
We explored the association of the SNPs identified in the present study for HC with multiple mental health and intelligence-related traits. A total of eight traits were collected, which could be classified into two categories. The first category contains two normal traits intelligence (Savage et al., 2018) and educational attainment (Lee et al., 2018). The second category contains the remaining six disease-related traits, including subjective wellbeing (SWB) (Okbay et al., 2016), broad depression (Howard et al., 2018), autism spectrum disorder (ASD) (Grove et al., 2017), loneliness (http://www.nealelab.is/uk-biobank), cognitive impairment (CI) (http://www.nealelab.is/uk-biobank) and Parkinson disease (PD) (http://www.nealelab.is/uk-biobank).
Functional Enrichment Analysis
We assessed the enrichment of the HC SNPs into different genomic features with the GWAS Analysis of Regulatory or Functional Information Enrichment with LD correction (GARFIELD) (Iotchkova et al., 2016). GARFIELD employs an adaptive permutation procedure that creates and utilizes a set of null variant that accounts for systematic differences in minor allele frequency, gene distance, and number of proxies in the test variant set (Iotchkova et al., 2016).
Genomic feature information was accessed from the software website (https://www.ebi.ac.uk/birney-srv/GARFIELD/). Genomic features include DNase I hypersensitive sites and genic features from three tissues (fetal brain/brain, fetal hippocampus, and cerebellum) to capture the cell-type-specific feature enrichment. The assessment was conducted at four p value thresholds (1.0 × 10−5, 1.0 × 10−6, 1.0 × 10−7, and 5.0 × 10−8). To account for correlation between genomic features, the effective number of genomic features was first estimated by GARFIELD, which was 35.32. Statistical significance was then declared at a Bonferroni-corrected threshold, i.e., 1.42 × 10−3 (0.05/35.32).
Candidate Gene Prioritization
Within each locus, we first defined a set of credible risk variants (CRVs) as those located within 500 kb of the lead SNP, and with p values within two orders of magnitude of the lead SNP (Michailidou et al., 2017). This is approximately equivalent to defining CRVs whose posterior probability of causality is within two orders of magnitude of the value of the lead SNP (Udler et al., 2010).
We incorporated five sources of information to seek evidence of a gene’s causality: i) being the gene nearest the lead CRV, ii) containing a missense coding CRV, (iii) being the target gene for any cis-eQTL CRV, iv) being prioritized by DEPICT, or v) being prioritized by SMR analysis.
Cis-eQTL information was obtained from two datasets: GTEx (v7) (https://gtexportal.org/home/) (GTEx, 2017) and the study of Westra et al. (Westra et al., 2013). Cis-eQTL information for GTEx is available for over 50 tissues, and for Westra et al., it is only available for whole blood. We explored all tissues to avoid missing a tissue-specific eQTL. We used the web tool HaploReg v4.1 (Ward et al., 2012) to annotate missense coding variants and cis-eQTLs (https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php). Significance threshold was set at 0.05.
DEPICT is an integrative tool that employs predicted gene functions to systematically prioritize the most likely causal genes at associated loci (Pers et al., 2015). Significant genes were declared at a false discovery rate (FDR) = 0.05.
SMR is another prioritization method that integrates summary-level data from GWAS with cis-eQTL data to identify genes whose expression levels are associated with trait because of pleiotropy (Zhu et al., 2016). Three cis-eQTL datasets were analyzed. The first was the study of Westra et al. (Westra et al., 2013), and the second and third ones were two brain-related tissues from Qi et al. brain eQTL results and the PsychENCODE dataset (Gandal et al., 2018; Qi et al., 2018). All cis-eQTL datasets were downloaded from the SMR website. The LD reference panel required by SMR was formed by the 503 European subject haplotypes of the 1000G project. Specifically, the length of the reference panel was 2 mega-base (MB) and was centered at the target lead SNP. Significance threshold was set at 0.05.
Gene Set Enrichment Analysis
We performed gene set enrichment analysis with MAGMA (v1.6) (de Leeuw et al., 2015). MAGMA simultaneously analyzes multiple genetic markers incorporating LD between markers to determine their joint effect utilizing either summary data or individual-level data (de Leeuw et al., 2015). The total SNPs were mapped to 16,645 genes according to de Leeuw et al. (2015)and these genes were mapped to 9,467 predefined gene sets retrieved from GO consortium (http://www.geneontology.org/). A gene-based association test was examined to determine individual gene association signals. A gene set association signal was then calculated by integrating all gene signals within the set. We used competitive model to test whether genes in a gene set are more strongly associated with the phenotype than other genes and multitesting problem was corrected by permutation test.
Results
Joint Association Analysis
After a series of inclusion criteria, we finally included 1,977,400 SNPs for analysis. Using HC summary association statistics at these SNPs, the LDSC estimated a SNP attributable heritability 0.23 (s.e. = 0.05) for the trait. Meanwhile, the SNP attributable heritability for BL and BW was 0.169 and 0.098, respectively. There are two loci 12q15 and 12q24 being significant at the genome-wide significance level in the original HC summary results. After excluding the two lead SNPs (rs1042725 and rs7980687) and their neighboring SNPs within 500 kb to either direction, the heritability is 0.22 (s.e. = 0.04), showing that most of the SNP-attributable heritability is hidden in the current GWAS setting.
When applied a LDSC-based genetic correlation analysis, the genetic correlation coefficient between HC and BL is estimated to be as high as 0.55 (s.e. = 0.12) and that between HC and BW is 0.37 (s.e. = 0.08), in line with previous research (Horikoshi et al., 2016) (Table 2). This implies that the three traits have shared genetic background.
Table 2 Genetic correlations of head circumference (HC) and other two related traits.
We then applied the MTAG method to the three studies for a joint analysis. Briefly, MTAG jointly analyzes GWAS summary statistics from different traits, possibly from overlapping samples, and reports association results for each trait. In the MTAG analysis of the HC trait, the LDSC intercept and mean chi-square are 0.99 and 1.04, suggesting that the vast majority of inflation in the mean chi-squared statistic is from polygenic architecture rather than from population stratification.
The two originally significant loci remain significant at the genome-wide significance level (GWS, 5 × 10−8), with the exact same two lead SNPs (rs1042725 and rs7980687) and stronger signals. To search for additional loci, we evaluated MTAG results of all SNPs excluding those within 500 kb of a locus that was reported previously. These identified 76 SNPs in 9 distinct loci that are significant at the GWS level. Of the nine lead SNPs, four are not significant in the original HC study at the nominal level (α = 0.05) but extremely significant for BW and BL, indicating that the significant signals are likely to be driven solely by BL and/or BW instead by HC. As a QC criterion, we discarded these four loci.
We then sought for in silico replication of the identified SNPs in an independent large-scale GWAS meta-analysis of intracranial volume (N = 26,577). Both previously identified SNPs are significant (one-sided p = 3.22 × 10−9 and 0.03) in the replication results, demonstrating their robustness. Of the five newly identified loci, all are consistent in effect direction between HC and intracranial volume. Among them, one-sided p value at three SNPs further achieves the Bonferroni-corrected significance level (α = 0.01), while rs3795128 has a weaker but nonsignificant signal (p = 0.02). Taken together, SNPs from three loci were convincingly associated with HC in both discovery and replication stages. The three lead SNPs explain 0.6% of HC phenotypic variance. Results are listed in Table 3 and Supplementary Table 1, and the Manhattan plot is displayed in Figure 1.
Table 3 The main results of the identified loci for HC.
Figure 1 Manhattan plot of head circumference (HC). Known loci were retrieved from the EBI GWAS catalog website. Novel single nucleotide polymorphisms (SNPs) below the significance threshold 5.0 × 10−8 are marked in green.
For BL and BW, new loci beyond the originally significant ones are identified, too. For BL, a total of 539 SNPs from 27 distinct loci are significant at the GWS level. Of them, only one was reported by the original study, while the remaining 26 loci are novel (Supplementary Table 2). For BW, a total of 731 SNPs from 40 distinct loci are significant at the GWS level. Of them, 39 loci are significant in the original BW results, and the remaining one locus 10q22.1 is novel (Supplementary Table 2). We also performed in silico replications in adult height and BW, respectively. As shown in Supplementary Table 2, up to 21 of the 26 newly identified loci for BL are significant in adult height, and the only one newly identified locus for BW is significant in the UK Biobank BW. These largely replicated results further strength the validity of the present integrative analysis.
Association With Brain/Mental Related Traits
We next evaluated the relation of the five identified HC SNPs (two known + three new) to several mental health and intelligence-related traits and diseases. We collected GWAS summary results for eight traits. p values of the five lead SNPs in these traits are listed in Table 4.
Table 4 The association of head circumference (HC) single nucleotide polymorphisms (SNPs) with other mental traits.
At the Bonferroni corrected significance level (0.05/40 = 1.25 × 10−3), any of the five SNPs is associated with at least one trait. The two known SNPs rs1042725 and rs7980687 are both associated with educational attainment (rs1042725 p = 7.63 × 10−6 and rs7980687 p = 3.01 × 10−25). In addition, rs1042725 is associated with subjective well-being (p = 6.63 × 10−4), and rs7980687 is associated with intelligence (p = 1.44 × 10−6). The three newly identified SNPs are associated with educational attainment (rs9846396 p = 1.44 × 10−7), intelligence (rs12534093 p = 3.25 × 10−4) or loneliness (rs7048271 p = 1.82 × 10−4). There is also a hint for the association of rs12534093 with loneliness, although the signal (p = 9.74 × 10−3) does not reach the significance threshold.
Notably, the three loci (3q23, 9q33.3, and 12q24.31) associated with educational attainment have been identified to be significant at the GWS level (Savage et al., 2018), with different lead SNPs. At 3q23 and 12q24.31, the lead SNPs in the present study are in high linkage disequilibrium (0.81–0.98) with the GWS lead SNPs, respectively, indicating that both signals arise from the same association signal at respective locus.
SNP Enrichment in Genomic Features
We investigated the enrichment of the identified SNPs into a variety of genomic features in three brain-related tissues including brain/fetal brain, fetal hippocampus, and cerebellum. The results are shown in Supplementary Table 3. At the most stringent threshold (T = 5.0 × 10−8), significant enrichment is observed in up to 13 features. Notably, SNPs are enriched in the DNase I hypersensitive site peak and a broader DNase I hypersensitive site category hotspot in only one (brain/fetal brain, OR: 10.71–30.44, p: 1.66 × 10−6–2.79 × 10−4) of the three tissues being investigated, but not in fetal hippocampus and cerebellum. At footprint feature, which is defined as a paucity of cleavage because protein-bound DNA is typically protected from DNase I cleavage, SNPs are observed to be enriched in brain (OR = 12.11, p = 7.62 × 10−4) and brain hippocampus (OR = 12.26, p = 3.39 × 10−4). SNPs are also enriched in exon (OR = 10.94, p = 1.21 × 10−3) and intergenic regions (OR = 10.70, p = 1.33 × 10−3) with lower ORs. At T = 1.0 × 10−7, significant enrichment is found in the same 13 features with similar ORs. At T = 1.0 × 10−6 and 1.0 × 10−5, however, enrichment is observed in a fewer number of features (9 and 0, respectively).
Candidate Gene Prioritization
We combined five sources of information to prioritize responsible candidate genes. The results are listed in Table 5. Sample size for the several tissues being analyzed varies from 175 to 491, all greater than the median across the total 53 tissues. A total of nine genes are prioritized. Among them, two genes (ZBTB38 and STRBP) are prioritized by more than one source. Both genes are closest to respectively lead SNPs and are the cis-eQTL target gene of the latter. What is more, previous research demonstrates that STRBP at 9q33.3 is differentially expressed between patients with Down’s syndrome and controls (Salemi et al., 2012).
Table 5 Prioritized candidate genes at the identified novel loci.
Gene Set Enrichment Analysis
Based on genome-wide SNP associations, we performed gene set enrichment analysis. Two GO terms, GO:0048009 and GO:0003680, are significant after permutation-based multiple-testing correction (Table 6). The term GO:0003680 consists of two genes KMT2A and HMGA2 whose p values are 0.15 and 2.55 × 10−11 in gene-based association results. Therefore, the signal of this gene set is presumably dominated by HMGA2. GO:0048009 is the series of molecular signals generated as a consequence of the insulin-like growth factor (IGF) receptor binding to one of its physiological ligands (http://amigo.geneontology.org/amigo/). The role of this gene set includes regulating the development and growth of central nervous system (O’Kusky et al., 2012). It may also contribute to the development of cognitive disorders (Xie et al., 2014; Pomytkin et al., 2018).
Table 6 Main results of gene set enrichment for head circumference (HC).
Discussion
In this study, we have performed a joint association analysis of HC and its two related traits BL and BW. We have identified three novel loci that are successfully replicated by intracranial volume in an independent large-scale study, increasing the number of associated loci from five to eight. We have also evaluated the enrichment of the associated SNPs into multiple genomic features, prioritized candidate genes, and identified responsible gene sets.
Our MTAG analysis identified four HC loci at the GWS level that are not significant in the original HC study even at the nominal level (α = 0.05). As a QC criterion, we discarded these four loci because they might represent false positive signals. As stated by the authors of MTAG, this false positive finding could happen when SNP is nonsignificant for one trait but significant for others in finite samples. Indeed, none of these SNPs could be replicated by intracranial volume (Supplementary Table 2), strengthening the possibility of false positive signals.
We could not replicate rs3795128 (20q12) in intracranial volume at Bonferroni-corrected threshold, although it was nominally significant (one-sided p = 0.02) in this dataset. We prioritized four genes (PLCG1, LPIN3, ZHX3, and TOP1) for this SNP using same gene prioritization strategy. Among them, PLCG1 is the most promising candidate gene at 20q12 detected by multiple sources because it is the closest to lead CRV (rs3795128) and contains a missense CRV (rs753381) that is in high linkage equilibrium (r2 = 0.7, 23 kb) with lead CRV. In a previous case–control study, missense mutations in PLCG1 were reported to participate in the pathogenesis of cutaneous T-cell lymphomas by enhancing downstream signaling towards nuclear factor of activated T-cell activation. We therefore hypothesized that rs3795128 was linked to T-cell lymphomas through its strong LD with the missense mutation in PLCG1 (Vaque et al., 2014). Evident impact of PLCG1 is repeatedly reported in learning and memory through various pathways in hippocampal tissue in mouse model (Gruart et al., 2007; Yang et al., 2017). These clues linking HC with neurological diseases imply that early brain growth may have a connection with cognitive ability. Therefore, PLCG1 may have a comprehensive regulatory effect on fetal growth and affect brain development in later life. Another candidate gene STRBP at 9q33.3 is reported to have a slight association with late-onset AD (Herold et al., 2016).
Given the connection between HC and cognitive ability and intelligence in epidemiology (Mosing et al., 2018), we explored this relation by investigating the association of the SNPs identified in the present study with brain health and intelligence-related traits/diseases. We observed an excess of significant SNPs in several traits including intelligence and educational attainment, in consistent with previous studies (Ivanovic et al., 2004).
The associated gene set GO:0048009 is an IGF receptor signaling pathway. In conjunction with its glycoprotein receptor (IGF-1R), IGF-1 is a phylogenetically ancient neurotrophic hormone and plays crucial roles in central nervous system development and maturation (Dyer et al., 2016). In astrocyte from mice, two previous studies demonstrated that cognition is influenced by IGF-1 by exhibiting expression of IGF-1R (Logan et al., 2018; Chen et al., 2019), as supported by both in animal and in clinical experiments (Baker et al., 2012; Andersson et al., 2018).
Certain limitations exist in the present study. First, the identified novel SNPs lack independent replication in the same trait. This is partly due to the relatively limited GWAS for HC in the community. As a compromise, we performed in silico replication in a genetically highly correlated trait intracranial volume. We believe that the replication results imply important clues towards true associations. Second, the unbalanced sample sizes used in the MTAG analysis may elevate false positive rate, as described by the MTAG authors (Scheffler et al., 2017). We tackled this issue by monitoring both the original and the integrative association signals for any unusual difference (e.g., non- to extreme significance). We also conducted an in-depth evaluation towards the relevance of the identified loci from a variety of sources, such as bioinformatical and functional investigation.
In summary, by integrating the GWAS summary statistics of HC, BL, and BW, we have performed a joint association analysis and identified three novel loci for HC. Not only that, 21 novel loci for BL and 1 for BW were also identified and replicated. We have also prioritized candidate genes and performed gene set enrichment analyses. Our results provide useful insights into the genetic architecture underlying HC and the development of its related mental conditions.
Data Availability Statement
Publicly available datasets were analyzed in this study. These data can be found here: http://www.gefos.org and http://egg-consortium.org.
Author Contributions
In this paper, X-LY is in charge of conceptualization/design, methodology, investigation, formal analysis, and resources. S-YZ contributed to methodology, formal analysis and funding for this study. HZ, X-TW and G-JF, are also co-authors in charge of conceptualization/design and data curation. The two corresponding authors, Y-FP and LZ are responsible for methodology, formal analysis, supervision/oversight, funding, and a final approval of the vision to be published.
Funding
This study was partially supported by the funding from National Natural Science Foundation of China (31571291 to LZ, 31771417 and 31501026 to Y-FP), the Natural Science Foundation of Jiangsu Province of China (BK20150323 to Y-FP), a project funded by the Priority Academic Program Development of Jiangsu higher education institutions (PAPD), and the Undergraduate Innovation Program of Jiangsu Province (KY2019881B to X-LY).
Conflict of Interest
The 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.
Acknowledgments
We are grateful to the Early Growth Genetics consortium for releasing genome-wide summary statistics for the three traits being analyzed. We are also grateful to the PGC, GIANT, CHARGE, and ENIGMA consortia and Dr. Neale for releasing summary results for replication. We appreciated efforts of Lu liu and Zhi-Zhen Cui to this work. The numerical calculations in this paper have been done on the supercomputing system of the National Supercomputing Center in Changsha.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00947/full#supplementary-material
References
Adams, H.H., Hibar, D.P., Chouraki, V., Stein, J.L., Nyquist, P.A., Rentería, M. E., et al. (2016). Novel genetic loci underlying human intracranial volume identified through genome-wide association. Nat. Neurosci. 19 (12), 1569–1582. doi: 10.1038/nn.4398
Andersson, K. M. E., Wasén, C., Juzokaite, L., Leifsdottir, L., Erlandsson, M. C., Silfverswärd, S. T., et al. (2018). Inflammation in the hippocampus affects IGF1 receptor signaling and contributes to neurological sequelae in rheumatoid arthritis. Proc. Natl. Acad. Sci. U. S. A. 115 (51), E12063–E12072. doi: 10.1073/pnas.1810553115
Aylward, E. H., Minshew, N. J., Field, K., Sparks, B. F., Singh, N. (2002). Effects of age on brain volume and head circumference in autism. Neurology 59 (2), 175–183. doi: 10.1212/WNL.59.2.175
Baker, L. D., Barsness, S. M., Borson, S., Merriam, G. R., Friedman, S. D., Craft, S., et al. (2012). Effects of growth hormone–releasing hormone on cognitive function in adults with mild cognitive impairment and healthy older adults: results of a controlled trial. Arch. Neurol. 69 (11), 1420–1429. doi: 10.1001/archneurol.2012.1970
Barker, D. J., Osmond, C., Simmonds, S. J., Wield, G. A. (1993). The relation of small head circumference and thinness at birth to death from cardiovascular disease in adult life. BMJ 306 (6875), 422–426. doi: 10.1136/bmj.306.6875.422
Bond, J., Roberts, E., Springell, K., Lizarraga, S., Scott, S., Higgins, J., et al. (2005). A centrosomal mechanism involving CDK5RAP2 and CENPJ controls brain size. Nat. Genet. 37 (4), 353–355. doi: 10.1038/ng1539
Breakey, R. W., Knoops, P. G., Borghi, A., Rodriguez-Florez, N., O’Hara, J., James, G., et al. (2018). Intracranial volume and head circumference in children with unoperated syndromic craniosynostosis. Plast. Reconstr. Surg. 142 (5), 708e–717e. doi: 10.1097/PRS.0000000000004843
Bulik-Sullivan, B., Finucane, H. K., Anttila, V., Gusev, A., Day, F. R., Loh, P-R., et al. (2015a). An atlas of genetic correlations across human diseases and traits. Nat. Genet. 47 (11), 1236–1241. doi: 10.1038/ng.3406
Bulik-Sullivan, B. K., Loh, P-R., Finucane, H., Ripke, S., Yang, J., Schizophrenia Working Group of the Psychiatric Genomics Consortium, et al. (2015b). LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. 47 (3), 291–295. doi: 10.1038/ng.3211
Chen, W., He, B., Tong, W., Zeng, J., Ping Zheng (2019). Astrocytic insulin-like growth factor-1 protects neurons against excitotoxicity. Front. Cell Neurosci. 13, 298. doi: 10.3389/fncel.2019.00298
de Leeuw, C. A., Mooij, J. M., Heskes, T., Posthuma, D. (2015). MAGMA: generalized gene-set analysis of GWAS data. PLoS Comput. Biol. 11 (4), e1004219. doi: 10.1371/journal.pcbi.1004219
Dyer, A. H., Vahdatpour, C., Sanfeliu, A., Tropea, D. (2016). The role of insulin-like growth factor 1 (IGF-1) in brain development, maturation and neuroplasticity. Neuroscience 325, 89–99. doi: 10.1016/j.neuroscience.2016.03.056
Gandal, M. J., Zhang, P., Hadjimichael, E., Walker, R. L., Chen, C., Liu, S., et al. (2018). Transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, and bipolar disorder. Science 362 (6420), eaat8127. doi: 10.1126/science.aat8127
Grove, J., Ripke, S., Als, T. D., Mattheisen, M., Walters, R., Won, H., et al. (2017). Common risk variants identified in autism spectrum disorder. bioRxiv p, 224774. doi: 10.1101/224774
Gruart, A., Sciarretta, C., Valenzuela-Harrington, M., Delgado-García, J. M., Minichiello, L. (2007). Mutation at the TrkB PLC{gamma}-docking site affects hippocampal LTP and associative learning in conscious mice. Learn Mem. 14 (1), 54–62. doi: 10.1101/lm.428307
GTEx Consortium (2017). Genetic effects on gene expression across human tissues. Nature 550 (7675), 204–213. doi: 10.1038/nature24277
Haksari, E. L., Lafeber, H. N., Hakimi, M., Pawirohartono, E. P., Nyström, L. (2016). Reference curves of birth weight, length, and head circumference for gestational ages in Yogyakarta, Indonesia. BMC Pediatr. 16 (1), 188. doi: 10.1186/s12887-016-0728-1
Haworth, S., Shapland, C. Y., Hayward, C., Prins, B. P., Felix, J. F., Medina-Gomez, C., et al. (2019). Low-frequency variation in TP53 has large effects on head circumference and intracranial volume. Nat. Commun. 10 (1), 357. doi: 10.1038/s41467-018-07863-x
Herold, C., Hooli, B. V., Mullin, K., Liu, T., Roehr, J. T., Mattheisen, M., et al. (2016). Family-based association analyses of imputed genotypes reveal genome-wide significant association of Alzheimer’s disease with OSBPL6, PTPRG, and PDCL3. Mol. Psychiatry 21 (11), 1608–1612. doi: 10.1038/mp.2015.218
Horikoshi, M., Beaumont, R. N., Day, F. R., Warrington, N. M., Kooijman, M. N., Fernandez-Tajes, J., et al. (2016). Genome-wide associations for birth weight and correlations with adult disease. Nature 538 (7624), 248–252. doi: 10.1038/nature19806
Howard, D. M., Adams, M. J., Shirali, M., Clarke, T-K., Marioni, R. E., Davies, G., et al. (2018). Genome-wide association study of depression phenotypes in UK Biobank identifies variants in excitatory synaptic pathways. Nat. Commun. 9 (1), 1470. doi: 10.1038/s41467-018-03819-3
Iotchkova, V., Ritchie, G. R. S., Geihs, M. G., Morganella, S., Min, J. L., Walter, K., et al. (2016). GARFIELD classifies disease-relevant genomic features through integration of functional annotations with association signals. Nat. Genet. 51 (2), 343–353.doi: 10.1038/s41588-018-0322-6
Ivanovic, D. M., Leiva, B. P., Pérez, H. T., Olivares, M. G., Díaz, N. S., Urruti, M. S. C., et al. (2004). Head size and intelligence, learning, nutritional status and brain development. Head, IQ, learning, nutrition and brain. Neuropsychologia 42 (8), 1118–1131. doi: 10.1016/j.neuropsychologia.2003.11.022
Kumar, A., Girimaji, S. C., Duvvari, M. R., Blanton, S. H. (2009). Mutations in STIL, encoding a pericentriolar and centrosomal protein, cause primary microcephaly. Am. J. Hum. Genet. 84 (2), 286–290. doi: 10.1016/j.ajhg.2009.01.017
Lee, J. J., Wedow, R., Okbay, A., Kong, E., Maghzian, O., Zacher, M., et al. (2018). Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat. Genet. 50 (8), 1112–1121. doi: 10.1038/s41588-018-0147-3
Logan, S., Pharaoh, G. A., Marlin, M. C., Masser, D. R., Matsuzaki, S., Wronowski, B., et al. (2018). Insulin-like growth factor receptor signaling regulates working memory, mitochondrial metabolism, and amyloid-beta uptake in astrocytes. Mol. Metab. 9, 141–155. doi: 10.1016/j.molmet.2018.01.013
Michailidou, K., Lindström, S., Dennis, J., Beesly, J., Hui, S., Kar, S., et al. (2017). Association analysis identifies 65 new breast cancer risk loci. Nature 551 (7678), 92–94. doi: 10.1038/nature24284
Mosing, M. A., Lundholm, C., Cnattingius, S., Gatz, M., Pedersen, N. L.(2018). Associations between birth characteristics and age-related cognitive impairment and dementia: a registry-based cohort study. PLoS Med. 15 (7), e1002609. doi: 10.1371/journal.pmed.1002609
O’Kusky, J., Ye, P. (2012). Neurodevelopmental effects of insulin-like growth factor signaling. Front. Neuroendocrinol. 33 (3), 230–251. doi: 10.1016/j.yfrne.2012.06.002
Okbay, A., Baselmans, B. M. L., De Neve, J. M., Turley, P., Nivard, M. G., Fontana, M. A., et al. (2016). Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses. Nat. Genet. 48 (6), 624–633. doi: 10.1038/ng.3552
Pers, T. H., Karjalainen, J. M., Chan, Y., Westra, H-J., Wood, A. R., Yang, J., et al. (2015). Biological interpretation of genome-wide association studies using predicted gene functions. Nat. Commun. 6, 5890. doi: 10.1038/ncomms6890
Pomytkin, I., Costa-Nunes, J. P., Kasatkin, V., Veniaminova, E., Lyundup, A. D. A., Lesch, K-P., et al. (2018). Insulin receptor in the brain: Mechanisms of activation and the role in the CNS pathology and treatment. CNS Neurosci. Ther. 24 (9), 763–774. doi: 10.1111/cns.12866
Qi, T., Wu, Y., Zeng, J., Zhang, F., Xue, A., Jiang, L., et al. (2018). Identifying gene targets for brain-related traits using transcriptomic and methylomic data from blood. Nat. Commun. 9 (1), 2282. doi: 10.1038/s41467-018-04558-1
Rydbeck, F., Rahman, A., Grandér, M., Ekström, E-C., Vahter, M., Kippler, M. (2014). Maternal urinary iodine concentration up to 1.0 mg/L is positively associated with birth weight, length, and head circumference of male offspring. J. Nutr 144 (9), 1438–1444. doi: 10.3945/jn.114.193029
Salemi, M., La Vignera, S., Castiglione, R., Condorelli, R. A., Cimino, L., Bosco, P., et al. (2012). Expression of STRBP mRNA in patients with cryptorchidism and Down’s syndrome. J. Endocrinol. Invest. 35 (1), 5–7. doi: 10.1007/BF03345414
Sankilampi, U., Hannila, M-L., Saari, A., Gissler, M., Dunkel, L. (2013). New population-based references for birth weight, length, and head circumference in singletons and twins from 23 to 43 gestation weeks. Ann. Med. 45 (5–6), 446–454. doi: 10.3109/07853890.2013.803739
Savage, J. E., Jansen, P. R., Stringer, S., Watanabe, K., Bryois, J., de Leeuw, C. A., et al. (2018). Genome-wide association meta-analysis in 269,867 individuals identifies new genetic and functional links to intelligence. Nat. Genet. 50 (7), 912–919. doi: 10.1038/s41588-018-0152-6
Scheffler, C., Greil, H., Hermanussen, M. (2017). The association between weight, height, and head circumference reconsidered. Pediatr. Res. 81 (5), 825–830. doi: 10.1038/pr.2017.3
Schofield, P. W., Logroscino, G., Andrews, H. F., Stern, A. Y. (1997). An association between head circumference and Alzheimer’s disease in a population-based study of aging and dementia. Neurology 49 (1), 30–37. doi: 10.1212/WNL.49.1.30
Smit, D. J. A., Luciano, M., Bartels, M., van Beijsterveldt, C. E. M., Wright, M. J., Hansell, N. K., et al. (2010). Heritability of head size in Dutch and Australian twin families at ages 0–50 years. Twin Res. Hum. Genet. 13 (4), 370–380. doi: 10.1375/twin.13.4.370
Taal, H. R., St Pourcain, B., Thiering, E., Das, S., Mook-Kanamori, D. O., Warrington, N. M., et al. (2012). Common variants at 12q15 and 12q24 are associated with infant head circumference. Nat. Genet. 44 (5), 532–538. doi: 10.1038/ng.2238
Turley, P., Walters, R. K., Maghzian, O., Okbay, A., Lee, J. J., Fontana, M. A., et al. (2018). Multi-trait analysis of genome-wide association summary statistics using MTAG. Nat. Genet. 50 (2), 229–237. doi: 10.1038/s41588-017-0009-4
Udler, M. S., Tyrer, J., Easton, D. F. (2010). Evaluating the power to discriminate between highly correlated SNPs in genetic association studies. Genet. Epidemiol. 34 (5), 463–468. doi: 10.1002/gepi.20504
van der Valk, R. J. P., Kreiner-Møller, E., Kooijman, M. N., Guxens, M., Stergiakouli, E., Sääf, A., et al. (2015). A novel common variant in DCST2 is associated with length in early life and height in adulthood. Hum. Mol. Genet. 24 (4), 1155–1168. doi: 10.1093/hmg/ddu510
Vaqué, J. P., Gómez-López, G., Monsálvez, V., Varela, I., Martínez, N., Pérez, C., et al. (2014). PLCG1 mutations in cutaneous T-cell lymphomas. Blood 123 (13), 2034–2043. doi: 10.1182/blood-2013-05-504308
Ward, L. D., Kellis, M. (2012). HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 40 (Database issue), D930–D934. doi: 10.1093/nar/gkr917
Westra, H. J., Peters, M. J., Esko, T., Yaghootkar, H., Schurmann, C., Kettunen, J., et al. (2013). Systematic identification of trans eQTLs as putative drivers of known disease associations. Nat. Genet. 45 (10), 1238–1243. doi: 10.1038/ng.2756
Xie, J., Wei, Q., Deng, H., Li, G., Ma, L., Zeng, H. (2014). Negative regulation of Grb10 Interacting GYF Protein 2 on insulin-like growth factor-1 receptor signaling pathway caused diabetic mice cognitive impairment. PLoS One 9 (9), e108559. doi: 10.1371/journal.pone.0108559
Yang, Y. R., Jung, J. H., Kim, S. J., Hamada, K., Suzuki, A., Kim, H. J. (2017). Forebrain-specific ablation of phospholipase Cgamma1 causes manic-like behavior. Mol. Psychiatry 22 (10), 1473–1482. doi: 10.1038/mp.2016.261
Yengo, L., Sidorenko, J., Kemper, K. E., Zheng, Z., Wood, A. R., Weedon, M. N., et al. (2018). Meta-analysis of genome-wide association studies for height and body mass index in approximately 700000 individuals of European ancestry. Hum. Mol. Genet. 27 (20), 3641–3649. doi: 10.1093/hmg/ddy271
Yokoyama, Y., Sugimoto, M., Ooki, S. (2005). Analysis of factors affecting birthweight, birth length and head circumference: study of Japanese triplets. Twin Res. Hum. Genet. 8 (6), 657–663. doi: 10.1375/twin.8.6.657
Keywords: infant head circumference, genome-wide association study, 3q23, 7p15.3, 9q33.3, birth weight, birth length
Citation: Yang X-L, Zhang S-Y, Zhang H, Wei X-T, Feng G-J, Pei Y-F and Zhang L (2019) Three Novel Loci for Infant Head Circumference Identified by a Joint Association Analysis. Front. Genet. 10:947. doi: 10.3389/fgene.2019.00947
Received: 02 May 2019; Accepted: 05 September 2019;
Published: 11 October 2019.
Edited by:
Marika Kaakinen, University of Surrey, United KingdomReviewed by:
Alexander Teumer, University of Greifswald, GermanyWouter Van Rheenen, University Medical Center Utrecht, Netherlands
Copyright © 2019 Yang, Zhang, Zhang, Wei, Feng, Pei and Zhang. 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: Yu-Fang Pei, eXBlaUBzdWRhLmVkdS5jbg==; Lei Zhang, bHpoYW5nNkBzdWRhLmVkdS5jbg==
†These authors have contributed equally to this work
‡These authors have jointly supervised this work