Genetic polymorphisms of GPR126 are functionally associated with PUMC classifications of adolescent idiopathic scoliosis in a Northern Han population

Association between genetic polymorphisms and risk of adolescent idiopathic scoliosis in case-control studies: a systematic review

BACKGROUND

Adolescent idiopathic scoliosis (AIS) is a structural lateral spinal curvature of ≥10° with rotation. Approximately 2%-3% of children across populations are affected with AIS, and this condition is responsible for ~$3 billion in costs within the USA. Although AIS is believed to have a strong genetic contribution, clinical translation of identified genetic variants has stalled.

METHODS

The databases MEDLINE (via PubMed), Embase, Google Scholar and Ovid MEDLINE were searched and limited to articles in English. Title and abstract, full-text and data extraction screening was conducted through Covidence, followed by data transfer to a custom REDCap database. Studies containing variant-level data using genome-wide methodology as well as validation studies of genome-wide methods were considered. Quality assessment was conducted using Q-Genie.

RESULTS

33 studies were included, including 9 genome-wide association studies, 4 whole exome sequencing and 20 validation studies. Combined, these studies included data from >35,000 cases and >67,000 controls, not including validation cohorts. Additionally, results from six meta-analyses containing novel cohorts were also reported. All included study cohorts were from populations of primarily East Asian or Caucasian descent. Quality assessment found that overall study quality was high and control group selection was moderate. The highest number of reported associations were in single nucleotide polymorphisms (SNPs) in or near LBX1, LBX1-AS1, GPR126/ADGRG6 or BNC2.

CONCLUSION

AIS risk may be influenced by specific SNPs, particularly those in/near LBX1 and GPR126. Translatability of study findings is unknown due to an underrepresentation of most ethnic groups as well as few identified genome-wide studies. Further studies may benefit from increased cohort diversity and thorough evaluation of control cohort groups.

Structure, ligands, and roles of GPR126/ADGRG6 in the development and diseases

Adhesion G protein-coupled receptors (aGPCRs) are the second largest diverse group within the GPCR superfamily, which play critical roles in many physiological and pathological processes through cell-cell and cell-extracellular matrix interactions. The adhesion GPCR Adgrg6, also known as GPR126, is one of the better-characterized aGPCRs. GPR126 was previously found to have critical developmental roles in Schwann cell maturation and its mediated myelination in the peripheral nervous system in both zebrafish and mammals. Current studies have extended our understanding of GPR126-mediated roles during development and in human diseases. In this review, we highlighted these recent advances in GPR126 in expression profile, molecular structure, ligand-receptor interactions, and associated physiological and pathological functions in development and diseases.

Analysis of GPR126 polymorphisms and their relationship with scoliosis in Marfan syndrome and Marfan-like syndrome in Mexican patients

Marfan syndrome (MFS) is an inherited connective tissue disorder. As the spinal growth depends on delicate balance of forces, conditions that affect musculoskeletal matrix often lead to spinal deformities. A large cross-sectional study revealed a 63% prevalence of scoliosis among patients with MFS. Multi-ethnic genome-wide association studies and analyses of human genetic mutations showed that variations and mutations of G protein-coupled receptor 126 (GPR126)locus are associated with multiple skeletal defects, including shorter stature and adolescent idiopathic scoliosis. The study included 54 patients with MFS and 196 control patients. The DNA was extracted from peripheral blood using the saline expulsion method and single nucleotide polymorphism (SNP) determination was carried out using TaqMan probes. Allelic discrimination was performed by RT-qPCR. Significant differences in genotype frequencies were found for SNP rs6570507 in relation to MFS and sex (recessive model, OR 2.46, 95% CI 1.03 -5.87; P = 0.03) and rs7755109 (overdominant model, OR 0.39, 95% CI 0.16-0.91; P = 0.03). The most significant association was found in SNP rs7755109, where the frequency of genotype AG was significantly different between MFS patients with scoliosis and those without (OR 5.68, 95% CI 1.09-29.48; P=0.04). This study, for the first time, examined the genetic association of SNP GPR126 with the risk of scoliosis in patients with connective tissue diseases. The study revealed that SNP rs7755109 is associated with the presence of scoliosis in Mexican patients with MFS.

Differential Regulation of POC5 by ERα in Human Normal and Scoliotic Cells

Adolescent idiopathic scoliosis (AIS) is a complex three-dimensional spinal deformity. The incidence of AIS in females is 8.4 times higher than in males. Several hypotheses on the role of estrogen have been postulated for the progression of AIS. Recently, Centriolar protein gene POC5 (POC5) was identified as a causative gene of AIS. POC5 is a centriolar protein that is important for cell cycle progression and centriole elongation. However, the hormonal regulation of POC5 remains to be determined. Here, we identify POC5 as an estrogen-responsive gene under the regulation of estrogen receptor ERα in normal osteoblasts (NOBs) and other ERα-positive cells. Using promoter activity, gene, and protein expression assays, we found that the POC5 gene was upregulated by the treatment of osteoblasts with estradiol (E2) through direct genomic signaling. We observed different effects of E2 in NOBs and mutant POC5^A429V AIS osteoblasts. Using promoter assays, we identified an estrogen response element (ERE) in the proximal promoter of POC5, which conferred estrogen responsiveness through ERα. The recruitment of ERα to the ERE of the POC5 promoter was also potentiated by estrogen. Collectively, these findings suggest that estrogen is an etiological factor in scoliosis through the deregulation of POC5.

Genetic animal modeling for idiopathic scoliosis research: history and considerations

PURPOSE

Idiopathic scoliosis (IS) is defined as a structural lateral spinal curvature ≥ 10° in otherwise healthy children and is the most common pediatric spinal deformity. IS is known to have a strong genetic component; however, the underlying etiology is still largely unknown. Animal models have been used historically to both understand and develop treatments for human disease, including within the context of IS. This intended audience for this review is clinicians in the fields of musculoskeletal surgery and research.

METHODS

In this review article, we synthesize current literature of genetic animal models of IS and introduce considerations for researchers.

RESULTS

Due to complex genetic and unique biomechanical factors (i.e., bipedalism) hypothesized to contribute to IS in humans, scoliosis is a difficult condition to replicate in model organisms.

CONCLUSION

We advocate careful selection of animal models based on the scientific question and introduce gaps and limitations in the current literature. We advocate future research efforts to include animal models with multiple characterized genetic or environmental perturbations to reflect current understanding of the human condition.

The N Terminus of Adhesion G Protein–Coupled Receptor GPR126/ADGRG6 as Allosteric Force Integrator

The adhesion G protein–coupled receptor (aGPCR) GPR126/ADGRG6 plays an important role in several physiological functions, such as myelination or peripheral nerve repair. This renders the receptor an attractive pharmacological target. GPR126 is a mechano-sensor that translates the binding of extracellular matrix (ECM) molecules to its N terminus into a metabotropic intracellular signal. To date, the structural requirements and the character of the forces needed for this ECM-mediated receptor activation are largely unknown. In this study, we provide this information by combining classic second-messenger detection with single-cell atomic force microscopy. We established a monoclonal antibody targeting the N terminus to stimulate GPR126 and compared it to the activation through its known ECM ligands, collagen IV and laminin 211. As each ligand uses a distinct mode of action, the N terminus can be regarded as an allosteric module that can fine-tune receptor activation in a context-specific manner.

Genetic variants associated with the occurrence and progression of adolescent idiopathic scoliosis: a systematic review protocol

BACKGROUND

Adolescent idiopathic scoliosis (AIS) is a structural lateral spinal curvature of ≥ 10° with rotation. Approximately 2-3% of children in most populations are affected with AIS, and this condition is responsible for approximately $1.1 billion in surgical costs to the US healthcare system. Although a genetic factor for AIS has been demonstrated for decades, with multiple potentially contributory loci identified across populations, treatment options have remained limited to bracing and surgery.

METHODS

The databases MEDLINE (via PubMed), Embase, Google Scholar, and Ovid MEDLINE will be searched and limited to articles in English. We will conduct title and abstract, full-text, and data extraction screening through Covidence, followed by data transfer to a custom REDCap database. Quality assessment will be confirmed by multiple reviewers. Studies containing variant-level data (i.e., GWAS, exome sequencing) for AIS subjects and controls will be considered. Outcomes of interest will include presence/absence of AIS, scoliosis curve severity, scoliosis curve progression, and presence/absence of nucleotide-level variants. Analyses will include odds ratios and relative risk assessments, and subgroup analysis (i.e., males vs. females, age groups) may be applied. Quality assessment tools will include GRADE and Q-Genie for genetic studies.

DISCUSSION

In this systematic review, we seek to evaluate the quality of genetic evidence for AIS to better inform research efforts, to ultimately improve the quality of patient care and diagnosis.

SYSTEMATIC REVIEW REGISTRATION

PROSPERO registration #CRD42021243253.

Aberrant interaction between mutated ADAMTSL2 and LTBP4 is associated with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a complex spinal structure deformity with a prevalence of 1%-3%. Genetic and hereditary factors have been associated with the etiology of AIS. However, previous studies mainly focused on common single nucleotide polymorphisms which confer modest disease risk. Recently, rare variants in FBN1 and other extracellular matrix genes have been implicated in AIS, suggesting a potential overlapping disease etiology between AIS and hereditary connective tissue disorders (HCTD). In this study, we systematically analyzed rare variants in a set of HCTD-related genes in 302 AIS patients using exome sequencing. We firstly focused on pathogenic variants based on a monogenic inheritance and identified nine disease-associated variants in FBN1, COL11A1, COL11A2 and TGFBR2. We then explored the potential interactions between variants in different genes based on the case-control statistics. We identified three ADAMTSL2-LTBP4 variant pairs in three AIS patients and none in controls. Furthermore, we revealed that the variant pairs identified in these genes could affect the interaction between ADAMTSL2 and LTBP4 and upregulate TGF-β signaling pathway in human fibroblasts. Our findings supported that the aberrant interaction between mutated ADAMTSL2 and LTBP4 was associated with AIS.

Severity of Idiopathic Scoliosis Is Associated with Differential Methylation: An Epigenome-Wide Association Study of Monozygotic Twins with Idiopathic Scoliosis

Epigenetic mechanisms may contribute to idiopathic scoliosis (IS). We identified 8 monozygotic twin pairs with IS, 6 discordant (Cobb angle difference > 10°) and 2 concordant (Cobb angle difference ≤ 2°). Genome-wide methylation in blood was measured with the Infinium HumanMethylation EPIC Beadchip. We tested for differences in methylation and methylation variability between discordant twins and tested the association between methylation and curve severity in all twins. Differentially methylated region (DMR) analyses identified gene promoter regions. Methylation at cg12959265 (chr. 7 DPY19L1) was less variable in cases (false discovery rate (FDR) = 0.0791). We identified four probes (false discovery rate, FDR < 0.10); cg02477677 (chr. 17, RARA gene), cg12922161 (chr. 2 LOC150622 gene), cg08826461 (chr. 2), and cg16382077 (chr. 7) associated with curve severity. We identified 57 DMRs where hyper- or hypo-methylation was consistent across the region and 28 DMRs with a consistent association with curve severity. Among DMRs, 21 were correlated with bone methylation. Prioritization of regions based on methylation concordance in bone identified promoter regions for WNT10A (WNT signaling), NPY (regulator of bone and energy homeostasis), and others predicted to be relevant for bone formation/remodeling. These regions may aid in understanding the complex interplay between genetics, environment, and IS.

Methylation of Estrogen Receptor 1 Gene in the Paraspinal Muscles of Girls with Idiopathic Scoliosis and Its Association with Disease Severity

Idiopathic scoliosis (IS) is a multifactorial disease with epigenetic modifications. Tissue dependent and differentially methylated regions (T-DMRs) may regulate tissue-specific expression of the estrogen receptor 1 gene (ESR1). This study aimed to analyze methylation levels within T-DMR1 and T-DMR2 and its concatenation with ESR1 expression of IS patients. The study involved 87 tissue samples (deep paravertebral muscles, both on the convex and the concave side of the curve, and from back superficial muscles) from 29 girls who underwent an operation due to IS. Patient subgroups were analyzed according to Cobb angle ≤70° vs. >70°. Methylation was significantly higher in the superficial muscles than in deep paravertebral muscles in half of the T-DMR1 CpGs and all T-DMR2 CpGs. The methylation level correlated with ESR1 expression level on the concave, but not convex, side of the curvature in a majority of the T-DMR2 CpGs. The T-DMR2 methylation level in the deep paravertebral muscles on the curvature's concave side was significantly lower in patients with a Cobb angle ≤70° in four CpGs. DNA methylation of the T-DMRs is specific to muscle tissue location and may be related to ESR1 expression regulation. Additionally, the difference in T-DMR2 methylation may be associated with IS severity.

The orthopedic characterization of cfap298tm304 mutants validate zebrafish to faithfully model human AIS

Cerebrospinal fluid (CSF) circulation relies on the beating of motile cilia projecting in the lumen of the brain and spinal cord cavities Mutations in genes involved in cilia motility disturb cerebrospinal fluid circulation and result in scoliosis-like deformities of the spine in juvenile zebrafish. However, these defects in spine alignment have not been validated with clinical criteria used to diagnose adolescent idiopathic scoliosis (AIS). The aim of this study was to describe, using orthopaedic criteria the spinal deformities of a zebrafish mutant model of AIS targeting a gene involved in cilia polarity and motility, cfap298tm304. The zebrafish mutant line cfap298tm304, exhibiting alteration of CSF flow due to defective cilia motility, was raised to the juvenile stage. The analysis of mutant animals was based on micro-computed tomography (micro-CT), which was conducted in a QUANTUM FX CALIPER, with a 59 µm-30 mm protocol. 63% of the cfap298tm304 zebrafish analyzed presented a three-dimensional deformity of the spine, that was evolutive during the juvenile phase, more frequent in females, with a right convexity, a rotational component and involving at least one dislocation. We confirm here that cfap298tm304 scoliotic individuals display a typical AIS phenotype, with orthopedic criteria mirroring patient's diagnosis.

From genetics to epigenetics to unravel the etiology of adolescent idiopathic scoliosis

Scoliosis is defined as the three-dimensional (3D) structural deformity of the spine with a radiological lateral Cobb angle (a measure of spinal curvature) of ≥10° that can be caused by congenital, developmental or degenerative problems. However, those cases whose etiology is still unknown, and affect healthy children and adolescents during growth, are the commonest form of spinal deformity, known as adolescent idiopathic scoliosis (AIS). In AIS management, early diagnosis and the accurate prediction of curve progression are most important because they can decrease negative long-term effects of AIS treatment, such as unnecessary bracing, frequent exposure to radiation, as well as saving the high costs of AIS treatment. Despite efforts made to identify a method or technique capable of predicting AIS progression, this challenge still remains unresolved. Genetics and epigenetics, and the application of machine learning and artificial intelligence technologies, open up new avenues to not only clarify AIS etiology, but to also identify potential biomarkers that can substantially improve the clinical management of these patients. This review presents the most relevant biomarkers to help explain the etiopathogenesis of AIS and provide new potential biomarkers to be validated in large clinical trials so they can be finally implemented into clinical settings.

On being the right shape: Roles for motile cilia and cerebrospinal fluid flow in body and spine morphology

How developing and growing organisms attain their proper shape is a central problem of developmental biology. In this review, we investigate this question with respect to how the body axis and spine form in their characteristic linear head-to-tail fashion in vertebrates. Recent work in the zebrafish has implicated motile cilia and cerebrospinal fluid flow in axial morphogenesis and spinal straightness. We begin by introducing motile cilia, the fluid flows they generate and their roles in zebrafish development and growth. We then describe how cilia control body and spine shape through sensory cells in the spinal canal, a thread-like extracellular structure called the Reissner fiber, and expression of neuropeptide signals. Last, we discuss zebrafish mutants in which spinal straightness breaks down and three-dimensional curves form. These curves resemble the common but little-understood human disease Idiopathic Scoliosis. Zebrafish research is therefore poised to make progress in our understanding of this condition and, more generally, how body and spine shape is acquired and maintained through development and growth.

Estrogen Receptors (ESRs) Mutations in Adolescent Idiopathic Scoliosis: A Cross-Sectional Study

Background

Adolescent idiopathic scoliosis (AIS) is the most common spinal deformity, but its etiology is unclear. Multiple genetic mutations have been reported to be associated with AIS.

Material/Methods

We enrolled a cohort of 113 surgically treated AIS patients with available parental subjects from the Peking Union Medical College Hospital. We performed whole-exome sequencing in 10 trio families and whole-genome sequencing in 103 singleton patients. Luciferase assay was used to detect the functional alterations of candidate ESR1 and ESR2 variants.

Results

Using a de novo strategy, a missense variant in ESR1 (c.868A>G) was selected as a candidate gene for AIS. The main Cobb angle of this patient was 41° (T6–T10). Another potential pathogenic variant in ESR2 (c.236T>C) was identified. The main curve of the patient was 45° at T10–L3. The transactivation capacities of the mutated ESR1 and ESR2 protein were both significantly decreased (p=0.026 and 0.014, respectively).

Conclusions

Potential pathogenic variants in ESR1 and ESR2 were identified in 113 AIS patients, suggesting that genetic mutations in ESR1/2 were associated with the risk of AIS.

Whole-Genome Methylation Analysis of Phenotype Discordant Monozygotic Twins Reveals Novel Epigenetic Perturbation Contributing to the Pathogenesis of Adolescent Idiopathic Scoliosis

Background: Adolescent idiopathic scoliosis (AIS) is a complex disease affecting a large number of teenagers, especially in female. This study reveals novel epigenetic perturbation to the pathogenesis of AIS.

Methods: A female monozygotic (MZ) twin pair discordant for AIS were examined for whole-exome sequencing and epigenome difference. Sets of differentially methylated regions (DMRs) were validated using MethylTarget™ method in 20 AIS female patients and 20 healthy female controls.

Results: Few exome difference but several potential DMRs were found between the MZ twins. We identified 313 hypermethylated DMRs and 397 hypomethylated DMRs, respectively. Most of them were enriched in the MAPK and PI3K-Akt signaling pathway, which may contribute to the discordance of AIS. Several DMRs related to scoliosis genes were tested, and the NDN: TSS-DMR (chr15:23932133-23932304, hg19) was confirmed in additional samples. The methylation level of this DMR was significantly higher in the AIS group than in the control group (p = 0.04).

Conclusions: We described the epigenome difference in an AIS female discordant MZ twin pair using Whole Genome Bisulfite Sequencing (WGBS). The NDN: TSS-DMR had higher methylation level in female AIS, which can help elucidate the potential etiology of AIS.

Paraspinal muscle ladybird homeobox 1 (LBX1) in adolescent idiopathic scoliosis: a cross-sectional study

BACKGROUND CONTEXT

Adolescent idiopathic scoliosis (AIS) is the leading cause of spinal deformity in adolescents globally. Recent evidence from genome-wide association studies has implicated variants in or near the ladybird homeobox 1 (LBX1) gene, encoding the ladybird homeobox 1 transcription factor, in AIS development. This gene plays a critical role in guiding embryonic neurogenesis and myogenesis and is vital in muscle mass determination. Despite the confirmation of the role for LBX1 gene variants in the development of AIS, the biological basis of LBX1 contribution to AIS remains mostly unknown.

PURPOSE

To investigate the potential role of LBX1 in driving spinal curving, curve laterality, and progression through muscle-based mechanisms in AIS patients by analyzing its gene and protein expression.

STUDY DESIGN

This is a cross-sectional study using clinical data and biological samples from the Immunometabolic CONnections to Scoliosis study (ICONS study).

PATIENT SAMPLE

Twenty-five patients with AIS provided informed consent. Paraspinal muscle biopsies from the maximal points of concavity and convexity for gene expression and protein analysis were obtained at the start of corrective spinal surgery.

OUTCOME MEASURES

The outcome measures included the detection of paraspinal muscle LBX1 mRNA abundance and LBX1 protein expression and the correlation of the latter with age, sex, and curve severity.

METHODS

The measurement of mRNA abundance was done using quantitative real-time polymerase chain reaction (qRT-PCR). Additionally, protein lysates from the biopsied muscle samples were probed with a monoclonal LBX1 antibody to compare the muscle protein levels on either side of the scoliotic curve by western blot. This study received funding from the Division of Orthopedics, Department of Surgery, McMaster University, Hamilton, Ontario, Canada ($39,900 CAN for 2 years). The authors have no conflicts of interest to disclose.

RESULTS

LBX1 mRNA abundance (concave 2.98±0.87, convex 3.40±1.10, p value 0.73) and protein expression (concave 1.20±0.13, convex 1.21±0.10, p value 0.43) were detected on both sides of the scoliotic curve at equivalent levels. The expression of LBX1 protein did not correlate with age (concave: correlation coefficient 0.32, p value 0.12; convex: correlation coefficient 0.08, p value 0.69), sex (concave: correlation coefficient -0.03, p value 0.08; convex: correlation coefficient 0.07, p value 0.72), or the severity of spinal curving measured using the Cobb angle (concave: correlation coefficient -0.16, p value 0.45; convex: correlation coefficient -0.08, p value 0.69).

CONCLUSIONS

LBX1 is expressed in erector spinae muscles, and its levels are equal in muscles on both sides of the scoliotic curve in AIS. The expression of LBX1 on the convex and concave sides of the scoliotic curve did not correlate with age, sex, or the severity of spinal curving. The molecular mechanisms by which LBX1contributes to the development and propagation of AIS need to be explored further in muscle and other tissues.

Adhesion G protein-coupled receptors: opportunities for drug discovery

Adhesion G protein-coupled receptors (aGPCRs) - one of the five main families in the GPCR superfamily - have several atypical characteristics, including large, multi-domain N termini and a highly conserved region that can be autoproteolytically cleaved. Although GPCRs overall have well-established pharmacological tractability, currently no therapies that target any of the 33 members of the aGPCR family are either approved or in clinical trials. However, human genetics and preclinical research have strengthened the links between aGPCRs and disease in recent years. This, together with a greater understanding of their functional complexity, has led to growing interest in aGPCRs as drug targets. A framework for prioritizing aGPCR targets and supporting approaches to develop aGPCR modulators could therefore be valuable in harnessing the untapped therapeutic potential of this family. With this in mind, here we discuss the unique opportunities and challenges for drug discovery in modulating aGPCR functions, including target identification, target validation, assay development and safety considerations, using ADGRG1 as an illustrative example.

The Progress of CRISPR/Cas9-Mediated Gene Editing in Generating Mouse/Zebrafish Models of Human Skeletal Diseases

Genetic factors play a substantial role in the etiology of skeletal diseases, which involve 1) defects in skeletal development, including intramembranous ossification and endochondral ossification; 2) defects in skeletal metabolism, including late bone growth and bone remodeling; 3) defects in early developmental processes related to skeletal diseases, such as neural crest cell (NCC) and cilia functions; 4) disturbance of the cellular signaling pathways which potentially affect bone growth. Efficient and high-throughput genetic methods have enabled the exploration and verification of disease-causing genes and variants. Animal models including mouse and zebrafish have been extensively used in functional mechanism studies of causal genes and variants. The conventional approaches of generating mutant animal models include spontaneous mutagenesis, random integration, and targeted integration via mouse embryonic stem cells. These approaches are costly and time-consuming. Recent development and application of gene-editing tools, especially the CRISPR/Cas9 system, has significantly accelerated the process of gene-editing in diverse organisms. Here we review both mice and zebrafish models of human skeletal diseases generated by CRISPR/Cas9 system, and their contributions to deciphering the underpins of disease mechanisms.

Replication Study for the Association of GWAS-associated Loci With Adolescent Idiopathic Scoliosis Susceptibility and Curve Progression in a Chinese Population

STUDY DESIGN

A genetic association (replication) study.

OBJECTIVE

The aim of this study was to replicate and further evaluate the association among seven genome-wide association studies (GWAS)-identified single nucleotide polymorphisms (SNPs) in Chinese girls with adolescent idiopathic scoliosis (AIS) with disease onset, curve types, and progression.

SUMMARY OF BACKGROUND DATA

AIS is the most common pediatric spinal deformity with a strong genetic predisposition. Recent GWAS identified 10 new disease predisposition loci for AIS.

METHODS

Three hundred nineteen female AIS patients with Cobb angle ≥ 10 and 201 healthy controls were studied for the association with disease onset. Seven GWAS-identified SNPs (rs11190870 in LBX1, rs12946942 in SOX9/KCNJ2, rs13398147 in PAX3/EPH4, rs241215 in AJAP1, rs3904778 in BNC2, rs6570507 in GPR126, and rs678741 in LBX1-AS1) were analyzed. In subgroup analysis, AIS patients were subdivided by curve types and disease progression to examine for genotype association.

RESULTS

We replicated the association with disease onset in four common SNPs rs11190870, rs3904778, rs6570507, and rs678741. In addition, rs1190870 and rs678741 remained significantly associated in the right thoracic curves only subgroup. However, no significant difference was observed with both clinical curve progression or Cobb angle.

CONCLUSION

This study replicated the associations of four GWAS-associated SNPs with occurrence of AIS in our Chinese population. However, none of these SNPs was associated with curve severity and progression. The results suggest that curve progression may be determined by environmental (nongenetic) factor, but further study with a larger sample size is required to address this issue.

LEVEL OF EVIDENCE

4.

A Genetic Variant in GPR126 Causing a Decreased Inclusion of Exon 6 Is Associated with Cartilage Development in Adolescent Idiopathic Scoliosis Population

Adolescent idiopathic scoliosis (AIS) is the most common spinal deformity disease in adolescents but its etiology and pathogenesis are still unclear. The current study aims to identify the relationship between single nucleotide polymorphisms (SNPs) of G protein-coupled receptor 126 (GPR126) gene and AIS predisposition. GPR126 contains 26 exons and alternative splicing of exon 6 and exon 25 produces 4 protein-coding transcripts. We genotyped SNPs of GPR126 gene around exon 6 and exon 25 in 131 Chinese AIS patients and 132 healthy controls and provided evidence that SNP rs41289839 G>A is strongly associated with AIS susceptibility. Linkage disequilibrium analysis suggests that rs41289839 and other AIS-related SNPs were in strong LD. Next, we demonstrated that rs41289839 G>A inhibits the inclusion of exon 6 during alternative splicing, resulting in a decreased expression level of exon 6-included transcript (GPR126-exon6ⁱⁿ) relative to the exon 6 excluded transcript (GPR126-exon6^ex) by minigene assay. Chondrogenic differentiation experiment showed that GPR126-exon6ⁱⁿ has a high expression level relative to GPR126-exon6^ex during chondrogenic differentiation of hMSCs. Our findings indicate that newly discovered SNP is related to cartilage development and may provide valuable insights into the etiology and pathogenesis of adolescent idiopathic scoliosis.