- Research article
- Open Access
- Open Peer Review
Lack of association between the CHL1 gene and adolescent idiopathic scoliosis susceptibility in Han Chinese: a case-control study
BMC Musculoskeletal Disorders volume 15, Article number: 38 (2014)
A previous genome-wide association study (GWAS) suggested a strong association between the single nucleotide polymorphism (SNP) rs10510181 in the proximity of the gene encoding a cell adhesion molecule with homology to L1CAM (CHL1) and adolescent idiopathic scoliosis (AIS) in Caucasians. To clarify the role of CHL1 in the etiopathogenesis of AIS, we performed a case-control replication study in a Han Chinese population.
Five hundred female AIS patients between 10 and 18 years of age, as well as 500 age- and sex-matched controls were included. This study was conducted as a 2-stage case-control analysis: initial screening for the association between AIS and SNPs in and around the CHL1 gene (186 cases and 169 controls) followed by a confirmation test (314 cases and 331 controls). rs10510181 and 4 SNPs (rs2055314, rs331894, rs2272522, and rs2272524) in the CHL1 gene were selected for genotyping.
Putative associations were shown between AIS and rs10510181, rs2055314, and rs2272522 in stage I. However, the associations were not confirmed in stage II. For rs10510181, the genotype frequencies were GG 28.8%, GA 46.2%, and AA 25.0% in AIS patients and GG 29.8%, GA 48.8%, and AA 21.4% in controls. No significant difference was found in genotype distribution between cases and controls (P = 0.39). Similarly, the genotype and allele distribution were comparable between case and control for rs2055314 and rs2272522.
There was no statistical association between polymorphisms of the CHL1 gene and idiopathic scoliosis in a Chinese population.
Adolescent idiopathic scoliosis (AIS) is a complex 3-dimensional deformity of the spine that occurs most commonly in girls during the peripubertal period . Despite decades of research, its specific cause has not been determined . To date, many factors have been proposed to be involved in its etiology, such as genetic factor, nervous system, skeletal growth, hormone, and metabolic dysfunction . It has been strongly suggested that genetic factor plays an important role in the development and progression of AIS [4, 5]. Many AIS susceptibility loci have been identified through genome-wide linkage analyses and genome-wide association study [6–8]. Several promising candidate genes for AIS susceptibility have been reported, including the genes for estrogen receptor , melatonin 1B receptor , matrilin-1  and ladybird homeobox 1 . However, to our knowledge, only few genes can be replicated for their association in different ethnic populations [12–14].
Recently, Sharma et al.  performed a genome-wide association study (GWAS) based on transmission-disequilibrium tests in 419 AIS families, which identified associated single nucleotide polymorphisms (SNPs) in the proximity of the cell adhesion molecule with homology to the L1CAM (CHL1) gene. They genotyped additional SNPs in chromosome 3p26.3 and tested the association in two follow-up case-control cohorts. They obtained the strongest association in rs10510181 with all three cohorts combined (P = 2.58 × 10-8). CHL1 encodes an axon guidance protein related to Robo3 , mutations of which could lead to horizontal gaze palsy with progressive scoliosis (HGPPS), a rare disease marked by severe scoliosis . To confirm the association between rs10510181 and AIS, and to explore whether the CHL1 gene is associated with the occurrence of AIS, we conducted the present genetic association study in a Han Chinese population.
This study initially screened 186 AIS patients and 169 controls for the association between AIS and SNPs in and around the CHL1 gene. Subsequently, the sample size was enlarged to 500 AIS patients and 500 controls for further validation. All patients were ethnically Han Chinese females aged between 11 and 18 years who were seen at our spine center between 2006 and 2010. The diagnosis of AIS was confirmed by a standing X-ray film of the whole spine. The curvature magnitude was measured by Cobb’s method. Only patients with a Cobb angle larger than 20 degrees were included in this study. Patients with congenital scoliosis, neuromuscular scoliosis, scoliosis with skeletal dysplasia, scoliosis with known endocrine and connective tissue abnormalities, or prior treatment for scoliosis were excluded. The control age- and sex-matched adolescent Han Chinese girls were recruited from local secondary schools through health examinations. All controls were examined with either a forward bending test or a radiograph if necessary to rule out any spine deformities. Any potential evidence of bone diseases, metabolic diseases, growth disturbances and other diseases known to affect normal bone metabolism were excluded. Informed consents for DNA analysis were obtained from all subjects or their parents. The protocols and the procedure were approved by the Committee on Medical Ethics of Nanjing Drum Tower Hospital.
The SNP rs10510181, the one that showed strongest association with AIS in the study by Sharma et al. , was selected for the analysis. Furthermore, four SNPs rs2055314, rs331894, rs2272522, and rs2272524 were selected from dbSNP (http://www.ncbi.nlm.nih.gov). These SNPs span the entire region around the CHL1 gene with average intervals of approximately 30.5 kb .
Total genomic DNA was extracted from peripheral blood samples using a DNA extraction kit (Promega, Madison, WI), following to the manufacturer’s instructions. Identification of the polymorphisms was performed using the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. The primers were listed in Table 1. PCR was carried out in a total volume of 20 μL, consisting of 10 μL of 2 × PCR mix, 0.2 μL of each oligonucleotide primer and 8.6 μL of sterile deionized H2O. A typical 35 amplification cycle consisted of 30 seconds at 96°C, 45 seconds at the annealing temperature, and 30 seconds at 72°C (Table 1). The final elongation step was 7 minutes at 72°C. For restriction enzyme digestion, PCR products were subjected to 3 to 5 U of the required enzyme in the presence of the accompanying buffer and incubated overnight. The digested PCR products were visualized using electrophoresis in 2% or 4% agarose gels. A total of 10% of the samples were tested twice to validate the genotyping results with a reproducibility of 100%.
The statistical analyses were performed using the SPSS software (version 13.0, Chicago, USA). The differences in genotype distributions and allelic frequencies between AIS patients and controls were examined using the χ2 test. The odds ratios (OR) and 95% confidence interval (95% CI) ranges were calculated by using logistic regression. Significance was considered at a P-value <0.05. Sample size calculation was done using Quanto (Version 1.2.4; USC, Los Angeles, CA). 473 cases and controls are needed to detect the same OR (1.49) as previously reported for rs10510181  with 80% power assuming MAF (minor allele frequency) of 10% for the log-additive model.
The mean age was 14.7 years for AIS patients and 14.6 years for the controls. No significant different was found for body height between the two groups. However, the body weight was significantly lower in AlS patients than in normal controls. The mean maximum Cobb angle was 33.8° in AIS patients (Table 2).
A total of 1000 subjects (500 cases and 500 controls) were successfully genotyped and subjected to the statistical analysis. The distributions of the alleles and genotypes for the five SNPs in stages I and II were presented in Table 3 and Table 4. No significant deviation of genotype frequencies from the Hardy-Weinberg equilibrium was noted in the AIS and control groups.
In stage I, the genotype frequencies for rs2055314 and rs2272522 were significantly different between case and control groups (P <0.05), and it was marginally different for rs10510181 (P = 0.08). Similarly, the allelic frequencies for rs2055314, rs2272522, and rs10510181 showed significant differences between case and control groups (P <0.05) (Table 3). Otherwise, the genotype and allelic frequencies for the other two SNPs (rs331894 and rs2272524) were comparable between case and control groups (P >0.05) (Table 3).
In stage II, the sample size was enlarged for rs2055314, rs2272522, and rs10510181. The pooled data were shown in Table 4. For rs10510181, the genotype frequencies were GG 28.8%, GA 46.2%, and AA 25.0% in AIS patients and GG 29.8%, GA 48.8%, and AA 21.4% in controls. No significant difference was found in genotype distribution between cases and controls (P = 0.39). In addition, there was no obvious difference in allelic distribution between the two groups (G: 51.9% in cases vs. 54.2% in controls, A: 48.1% in cases vs. 45.8% in controls, P = 0.30). Similarly, the genotype and allelic distributions for rs2055314 and rs2272522 were comparable between case and control groups (Table 4).
It is generally recognized that abnormal growth is associated with the development and progression of the scoliotic curves. After the onset of puberty, significantly longer corrected height, arm span, and various body segments were found. Furthermore, BMI of AIS was significantly lower than that of controls . In the present study, we also found that the body weight was significantly lower in AlS patients than in normal controls.
GWAS is a powerful method for the detection of genetic contributions to polygenic diseases, and it has been used to identify genetic predisposition in AIS, which has been recently regarded as one of the most common complex genetic disorders of the musculoskeletal system [6, 15]. However, this method may produce spurious associations, as the findings could be confounded by systematic biases in comparison groups, such as population stratification and effect heterogeneity [20, 21]. Therefore, replicating the associations in different ethnic groups, as well as studies with large sample sizes, is important to confirm the results of GWAS . Recently, Sharma et al.  reported the association between SNP rs10510181, in the proximity of the CHL1 gene, and AIS by performing GWAS in 419 Caucasian families with AIS. The present study was conducted to validate the association of SNP rs10510181 with AIS. The results showed that rs10510181 was not associated with AIS in a Han Chinese population. Furthermore, 4 SNPs covering the entire region of the CHL1 gene were not associated with the occurrence of AIS either, indicating that the CHL1 gene might not be associated with AIS.
The previously reported association of rs10510181 with AIS predisposition in Caucasians was not replicated in our study. Several factors could lead to this lack of replication, which is not rare in genetic association studies of common diseases. First, the allele and genotype frequencies of rs10510181 could vary with ethnicity. Because of differences in variant frequency and underlying linkage disequilibrium (LD) structure between ancestral populations, the relationship between rs10510181 and a pathogenic variant may be influenced. Second, the inclusion criteria were different between the two studies. The subjects included in Sharma's study  had a curvature larger than 10 degrees, as well as axial rotation toward the side of the deviation. In the present study, patients with a Cobb angle greater than 20 degrees were included.
The CHL1 gene encodes an axon guidance protein that plays an important role in the guidance of thalamocortical axons , as well as the proliferation and differentiation of neural progenitor cells . Mutations in the CHL1 gene may alter the axonal guidance and brain anatomy in mice . Moreover, mutations in the Robo3 gene, which encodes an axon guidance protein related to CHL1, could cause HGPPS . Among the hypotheses of AIS etiology, abnormalities in the central nervous system have long been thought to play a key role [26, 27]. Disturbance of the central nervous system may impair somatosensory function and motor adaptation, which lead to the asymmetry of neuromuscular condition causing AIS [27, 28]. Thus, it seems quite plausible that CHL1 could be involved in the etiology of scoliosis through disturbance of the central nervous system. Although the present study did not find any association between the four SNPs spanning the whole region of the CHL1 gene and AIS, we cannot entirely exclude the role of this gene in the development of AIS because the CHL1 gene spans more than 100 kb and contains thousands of SNPs.
The current study did not find an association between the SNPs in and around the CHL1 gene and AIS predisposition. Considering the power calculation and sample size of the present study, we conclude that there is no association between the CHL1 gene and AIS predisposition in a Han Chinese population. However, the role of the CHL1 gene in other populations cannot be excluded, and replication studies in other ethnic groups are needed to understand the overall significance of the CHL1 gene in AIS.
Asher MA, Burton DC: Adolescent idiopathic scoliosis: natural history and long term treatment effects. Scoliosis. 2006, 1 (1): 2-10.1186/1748-7161-1-2.
Lowe TG, Edgar M, Margulies JY, Miller NH, Raso VJ, Reinker KA, Rivard CH: Etiology of idiopathic scoliosis: current trends in research. J Bone Joint Surg Am. 2000, 82-A (8): 1157-1168.
Wang WJ, Yeung HY, Chu WC, Tang NL, Lee KM, Qiu Y, Burwell RG, Cheng JC: Top theories for the etiopathogenesis of adolescent idiopathic scoliosis. J Pediatr Orthop. 2011, 31 (1 Suppl): S14-27.
Kesling KL, Reinker KA: Scoliosis in twins. A meta-analysis of the literature and report of six cases. Spine (Phila Pa 1976). 1997, 22 (17): 2009-2014. 10.1097/00007632-199709010-00014. discussion 2015
Riseborough EJ, Wynne-Davies R: A genetic survey of idiopathic scoliosis in Boston, Massachusetts. J Bone Joint Surg Am. 1973, 55 (5): 974-982.
Takahashi Y, Kou I, Takahashi A, Johnson TA, Kono K, Kawakami N, Uno K, Ito M, Minami S, Yanagida H, et al: A genome-wide association study identifies common variants near LBX1 associated with adolescent idiopathic scoliosis. Nat Genet. 2011, 43 (12): 1237-1240. 10.1038/ng.974.
Miller NH, Justice CM, Marosy B, Doheny KF, Pugh E, Zhang J, Dietz HC, Wilson AF: Identification of candidate regions for familial idiopathic scoliosis. Spine (Phila Pa 1976). 2005, 30 (10): 1181-1187. 10.1097/01.brs.0000162282.46160.0a.
Ocaka L, Zhao C, Reed JA, Ebenezer ND, Brice G, Morley T, Mehta M, O'Dowd J, Weber JL, Hardcastle AJ, et al: Assignment of two loci for autosomal dominant adolescent idiopathic scoliosis to chromosomes 9q31.2-q34.2 and 17q25.3-qtel. J Med Genet. 2008, 45 (2): 87-92.
Wu J, Qiu Y, Zhang L, Sun Q, Qiu X, He Y: Association of estrogen receptor gene polymorphisms with susceptibility to adolescent idiopathic scoliosis. Spine (Phila Pa 1976). 2006, 31 (10): 1131-1136. 10.1097/01.brs.0000216603.91330.6f.
Qiu XS, Tang NL, Yeung HY, Lee KM, Hung VW, Ng BK, Ma SL, Kwok RH, Qin L, Qiu Y, et al: Melatonin receptor 1B (MTNR1B) gene polymorphism is associated with the occurrence of adolescent idiopathic scoliosis. Spine (Phila Pa 1976). 2007, 32 (16): 1748-1753. 10.1097/BRS.0b013e3180b9f0ff.
Chen Z, Tang NL, Cao X, Qiao D, Yi L, Cheng JC, Qiu Y: Promoter polymorphism of matrilin-1 gene predisposes to adolescent idiopathic scoliosis in a Chinese population. Eur J Hum Genet. 2009, 17 (4): 525-532. 10.1038/ejhg.2008.203.
Takahashi Y, Matsumoto M, Karasugi T, Watanabe K, Chiba K, Kawakami N, Tsuji T, Uno K, Suzuki T, Ito M, et al: Replication study of the association between adolescent idiopathic scoliosis and two estrogen receptor genes. J Orthop Res. 2011, 29 (6): 834-837. 10.1002/jor.21322.
Takahashi Y, Matsumoto M, Karasugi T, Watanabe K, Chiba K, Kawakami N, Tsuji T, Uno K, Suzuki T, Ito M, et al: Lack of association between adolescent idiopathic scoliosis and previously reported single nucleotide polymorphisms in MATN1, MTNR1B, TPH1, and IGF1 in a Japanese population. J Orthop Res. 2011, 29 (7): 1055-1058. 10.1002/jor.21347.
Jiang H, Qiu X, Dai J, Yan H, Zhu Z, Qian B, Qiu Y: Association of rs11190870 near LBX1 with adolescent idiopathic scoliosis susceptibility in a Han Chinese population. Eur Spine J. 2013, 22 (2): 282-286. 10.1007/s00586-012-2532-4.
Sharma S, Gao X, Londono D, Devroy SE, Mauldin KN, Frankel JT, Brandon JM, Zhang D, Li QZ, Dobbs MB, et al: Genome-wide association studies of adolescent idiopathic scoliosis suggest candidate susceptibility genes. Hum Mol Genet. 2011, 20 (7): 1456-1466. 10.1093/hmg/ddq571.
Chen S, Mantei N, Dong L, Schachner M: Prevention of neuronal cell death by neural adhesion molecules L1 and CHL1. J Neurobiol. 1999, 38 (3): 428-439. 10.1002/(SICI)1097-4695(19990215)38:3<428::AID-NEU10>3.0.CO;2-6.
Jen JC, Chan WM, Bosley TM, Wan J, Carr JR, Rub U, Shattuck D, Salamon G, Kudo LC, Ou J, et al: Mutations in a human ROBO gene disrupt hindbrain axon pathway crossing and morphogenesis. Science. 2004, 304 (5676): 1509-1513. 10.1126/science.1096437.
Chen QY, Chen Q, Feng GY, Lindpaintner K, Chen Y, Sun X, Chen Z, Gao Z, Tang J, He L: Case-control association study of the close homologue of L1 (CHL1) gene and schizophrenia in the Chinese population. Schizophr Res. 2005, 73 (2–3): 269-274.
Siu King Cheung C, Tak Keung Lee W, Kit Tse Y, Ping Tang S, Man Lee K, Guo X, Qin L, Chun Yiu Cheng J: Abnormal peri-pubertal anthropometric measurements and growth pattern in adolescent idiopathic scoliosis: a study of 598 patients. Spine (Phila Pa 1976). 2003, 28 (18): 2152-2157. 10.1097/01.BRS.0000084265.15201.D5.
Kingsmore SF, Lindquist IE, Mudge J, Beavis WD: Genome-wide association studies: progress in identifying genetic biomarkers in common, complex diseases. Biomark Insights. 2007, 2: 283-292.
McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JP, Hirschhorn JN: Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008, 9 (5): 356-369. 10.1038/nrg2344.
Siontis KC, Patsopoulos NA, Ioannidis JP: Replication of past candidate loci for common diseases and phenotypes in 100 genome-wide association studies. Eur J Hum Genet. 2010, 18 (7): 832-837. 10.1038/ejhg.2010.26.
Wright AG, Demyanenko GP, Powell A, Schachner M, Enriquez-Barreto L, Tran TS, Polleux F, Maness PF: Close homolog of L1 and neuropilin 1 mediate guidance of thalamocortical axons at the ventral telencephalon. J Neurosci. 2007, 27 (50): 13667-13679. 10.1523/JNEUROSCI.2888-07.2007.
Huang X, Zhu LL, Zhao T, Wu LY, Wu KW, Schachner M, Xiao ZC, Fan M: CHL1 negatively regulates the proliferation and neuronal differentiation of neural progenitor cells through activation of the ERK1/2 MAPK pathway. Mol Cell Neurosci. 2011, 46 (1): 296-307. 10.1016/j.mcn.2010.09.013.
Heyden A, Angenstein F, Sallaz M, Seidenbecher C, Montag D: Abnormal axonal guidance and brain anatomy in mouse mutants for the cell recognition molecules close homolog of L1 and NgCAM-related cell adhesion molecule. Neuroscience. 2008, 155 (1): 221-233. 10.1016/j.neuroscience.2008.04.080.
Potoupnis M, Kapetanos G, Kimiskidis VK, Symeonides PP: Is the central nervous system a causative factor in idiopathic scoliosis?. Stud Health Technol Inform. 2002, 91: 10-11.
Herman R, Mixon J, Fisher A, Maulucci R, Stuyck J: Idiopathic scoliosis and the central nervous system: a motor control problem. The Harrington lecture. Scoliosis Research Society. Spine (Phila Pa 1976) 1985. 1983, 10 (1): 1-14.
Veldhuizen AG, Wever DJ, Webb PJ: The aetiology of idiopathic scoliosis: biomechanical and neuromuscular factors. Eur Spine J. 2000, 9 (3): 178-184. 10.1007/s005860000142.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2474/15/38/prepub
We gratefully recognize the financial support of the Nanjing Medical Science and technique Development Foundation (JQX12005) and the National Natural Science Foundation of China (No. 81171767).
The authors declare that they have no competing interests.
XSQ and FL performed the genetic studies. ZZZ and BPQ recruited patients and were involved in the clinical part of the investigation. BW performed the statistical analysis. YQ conceived the study, participated in its design and coordination and helped to draft the manuscript. All of the authors read and approved the final manuscript.