Value of DNA methylation in predicting curve progression in patients with adolescent idiopathic scoliosis

Identification of Epigenetic Biomarkers of Adolescent Idiopathic Scoliosis Progression: A Workflow to Assess Local Gene Expression

Adolescent idiopathic scoliosis (AIS) is a three-dimensional structural deformity of the spine that affects 2-3% of adolescents under the age of 16. AIS etiopathogenesis is not completely understood; however, the disease phenotype is correlated to multiple genetic loci and results from genetic-environmental interactions. One of the primary, still unresolved issues is the implementation of reliable diagnostic and prognostic markers. For clinical management improvement, predictors of curve progression are particularly needed. Recently, an epigenetic contribution to AIS development and progression was proposed; nevertheless, validation of data obtained in peripheral tissues and identification of the specific mechanisms and genes under epigenetic control remain limited. In this study, we propose a methodological approach for the identification of epigenetic markers of AIS progression through an original workflow based on the preliminary characterization of local expression of candidate genes in tissues directly involved in the pathology. The feasibility of the proposed methodological protocol has been originally tested here in terms of identification of the putative epigenetic markers of AIS progression, collection of the different tissues, retrieval of an appropriate amount and quality of RNA and DNA, and identification of suitable reference genes.

Prevalence and plasma exosome-derive microRNA diagnostic biomarker screening of adolescent idiopathic scoliosis in Yunnan Province, China

Background

Idiopathic scoliosis significantly affects the physical and mental health of children and adolescents, with varying prevalence rates in different regions. The occurrence of idiopathic scoliosis is associated with genetic regulation and biochemical factors, but the changes in exosome-derived miRNA profiles among idiopathic scoliosis patients remain unclear. This study aimed to determine the prevalence of idiopathic scoliosis in Yunnan Province, China, and identify key exosome-derived miRNAs in idiopathic scoliosis through a cohort study.

Methods

From January 2018 to December 2020, a cross-sectional study on idiopathic scoliosis in children and adolescents was conducted in Yunnan Province. A total of 84,460 students from 13 cities and counties in Yunnan Province participated in a scoliosis screening program, with ages ranging from 7 to 19 years. After confirmation through screening and imaging results, patients with severe idiopathic scoliosis and normal control individuals were selected using propensity matching. Subsequently, plasma exosome-derived miRNA sequencing and RT-qPCR validation were performed separately. Based on the validation results, diagnostic performance analysis and target gene prediction were conducted for differential plasma exosome-derived miRNAs.

Results

The overall prevalence of idiopathic scoliosis in children and adolescents in Yunnan Province was 1.10%, with a prevalence of 0.87% in males and 1.32% in females. The peak prevalence was observed at age 13. Among patients diagnosed with idiopathic scoliosis, approximately 12.8% had severe cases, and there were more cases of double curvature than of single curvature, with thoracolumbar curvature being the most common in the single-curvature group. Sequencing of plasma exosome-derived miRNAs associated with idiopathic scoliosis revealed 56 upregulated and 153 downregulated miRNAs. Further validation analysis confirmed that hsa-miR-27a-5p, hsa-miR-539-5p, and hsa-miR-1246 have potential diagnostic value.

Conclusions

We gained insights into the epidemiological characteristics of idiopathic scoliosis in Yunnan Province and conducted further analysis of plasma exosome-derived miRNA changes in patients with severe idiopathic scoliosis. This study has provided new insights for the prevention and diagnosis of idiopathic scoliosis, paving the way for exploring clinical biomarkers and molecular regulatory mechanisms. However, further validation and elucidation of the detailed biological mechanisms underlying these findings will be required in the future.

Integrated DNA methylation analysis reveals a potential role for PTPRN2 in Marfan syndrome scoliosis

Background

Marfan syndrome (MFS) is a rare genetic disorder caused by mutations in the Fibrillin-1 gene (FBN1) with significant clinical features in the skeletal, cardiopulmonary, and ocular systems. To gain deeper insights into the contribution of epigenetics in the variability of phenotypes observed in MFS, we undertook the first analysis of integrating DNA methylation and gene expression profiles in whole blood from MFS and healthy controls (HCs).

Methods

The Illumina 850K (EPIC) DNA methylation array was used to detect DNA methylation changes on peripheral blood samples of seven patients with MFS and five HCs. Associations between methylation levels and clinical features of MFS were analyzed. Subsequently, we conducted an integrated analysis of the outcomes of the transcriptome data to analyze the correlation between differentially methylated positions (DMPs) and differentially expressed genes (DEGs) and explore the potential role of methylation-regulated DEGs (MeDEGs) in MFS scoliosis. The weighted gene co-expression network analysis was used to find gene modules with the highest correlation coefficient with target MeDEGs to annotate their functions in MFS.

Results

Our study identified 1253 DMPs annotated to 236 genes that were primarily associated with scoliosis, cardiomyopathy, and vital capacity. These conditions are typically associated with reduced lifespan in untreated MFS. We calculated correlations between DMPs and clinical features, such as cobb angle to evaluate scoliosis and FEV1% to assess pulmonary function. Notably, cg20223687 (PTPRN2) exhibited a positive correlation with cobb angle of scoliosis, potentially playing a role in ERKs inactivation.

Conclusions

Taken together, our systems-level approach sheds light on the contribution of epigenetics to MFS and offers a plausible explanation for the complex phenotypes that are linked to reduced lifespan in untreated MFS patients.

Knowledge mapping of idiopathic scoliosis genes and research hotspots (2002-2022): a bibliometric analysis

Background

The etiology of idiopathic scoliosis (IS) remains unclear. Gene-based studies on genetic etiology and molecular mechanisms have improved our understanding of IS and guided treatment and diagnosis. Therefore, it is imperative to explicate and demarcate the preponderant areas of inquiry, key scholars, and their aggregate scholarly output, in addition to the collaborative associations amongst publications or researchers.

Methods

Documents were retrieved from the Web of Science Core Collection (WoSCC) with the following criteria: TS = ("idiopathic scoliosis" AND gene) refined by search operators (genomic OR "hereditary substance" OR "germ plasm" OR Cistrons OR genetics OR genetic OR genes OR Polygenic OR genotype OR genome OR allele OR polygenes OR Polygene) AND DOCUMENT TYPES (ARTICLE OR REVIEW), and the timespan of 2002-01-01 to 2022-11-26. The online bibliometric analysis platform (bibliometric), bibliographic item co-occurrence matrix builder (BICOMB), CiteSpace 6.1. R6 and VOS viewer were used to evaluate articles for publications, nations, institutions, journals, references, knowledge bases, keywords, and research hotspots.

Results

A total of 479 documents were retrieved from WoSCC. Fourty-four countries published relevant articles. The country with the most significant number of articles was China, and the institution with the most significant number of articles was Nanjing University. Citation analysis formed eight meaningful clusters and 16 high-frequency keywords. (2) The citation knowledge map included single nucleotide polymorphisms, whole exome sequencing, axonal dynamin, drug development, mesenchymal stem cells, dietary intake, curve progression, zebrafish development model, extracellular matrix, and rare variants were the current research hotspots and frontiers.

Conclusions

Recent research has focused on IS-related genes, whereas the extracellular matrix and unusual variants are research frontiers and hotspots. Functional analysis of susceptibility genes will prove to be valuable for identifying this disease.

Pediatric Sleep Questionnaire for Sleep Apnea in Newly Diagnosed Adolescent Idiopathic Scoliosis Patients

Close association has been established between obstructive sleep apnea (OSA) and adolescent idiopathic scoliosis (AIS), with PSQ being employed as a screening method for OSA. A cross-sectional study was conducted among patients aged from 10 to 16 years who presented to a scoliosis outpatient clinic. Patient demographics, radiological assessments, and PSQ scores were gathered. A total of 299 patients were included in the study, with 28.7% males and 71.2% females. The average Cobb angle was 6.20°. PSQ scores revealed a prevalence of 33.4% for significant obstructive sleep apnea. Patients diagnosed with AIS exhibited a prevalence of 32.9% with positive PSQ results. Among those undergoing adenoid and/or tonsil surgery, 27% had positive PSQ scores. Factors such as genetics, abnormal biomechanical forces, environmental factors including melatonin, and intermittent hypoxia were explored for their potential contribution to AIS etiology. The aim of the study is to underscore the importance of early detection and intervention in OSA cases and highlights the effectiveness of the PSQ, as a screening tool in identifying sleep disorders. The findings underscore the complex relationship between OSA and AIS, and moreover any spinal curvature is in relation with OSA.

Advances in epigenetic research of adolescent idiopathic scoliosis and congenital scoliosis

Scoliosis is a three-dimensional structural deformity of the spine; more than 80% of scoliosis has no specific pathogenesis but is understood to be closely related to genetic, hormonal, and environmental factors. In recent years, the epigenetic alterations observed in scoliosis have been analyzed in numerous studies to determine the pathogenesis and progression of this condition, however, there is currently no comprehensive review of the epigenetic factors to date. We searched PubMed, Embase, and Web of Science databases for relative studies without language and date restrictions in March 2023. Twenty-five studies were included in this review and analyzed from the four main aspects of epigenetic alteration: DNA methylation, non-coding RNAs, histone modifications, and chromatin remodeling. The relationship between DNA methylation, non-coding RNAs, and scoliosis was considerably reported in the literature, and the corresponding related signaling pathways and novel biomarkers observed in scoliosis provide insights into innovative prevention and treatment strategies. However, the role of histone modifications is rarely reported in scoliosis, and few studies have investigated the relationship between scoliosis and chromatin remodeling. Therefore, these related fields need to be further explored to elucidate the overall effects of epigenetics in scoliosis.

The asymmetrical ESR1 signaling in muscle progenitor cells determines the progression of adolescent idiopathic scoliosis

Adolescent Idiopathic Scoliosis (AIS) is a common pediatric skeletal disease highly occurred in females. The pathogenesis of AIS has not been fully elucidated. Here, we reveal that ESR1 (Estrogen Receptor 1) expression declines in muscle stem/progenitor cells at the concave side of AIS patients. Furthermore, ESR1 is required for muscle stem/progenitor cell differentiation and disrupted ESR1 signaling leads to differentiation defects. The imbalance of ESR1 signaling in the para-spinal muscles induces scoliosis in mice, while reactivation of ESR1 signaling at the concave side by an FDA approved drug Raloxifene alleviates the curve progression. This work reveals that the asymmetric inactivation of ESR1 signaling is one of the causes of AIS. Reactivation of ESR1 signaling in para-spinal muscle by Raloxifene at the concave side could be a new strategy to treat AIS.

Identification of a Functional Susceptibility Variant for Adolescent Idiopathic Scoliosis that Upregulates Early Growth Response 1 (EGR1)-Mediated UNCX Expression

Adolescent idiopathic scoliosis (AIS) is a serious health problem affecting 3% of live births all over the world. Many loci associated with AIS have been identified by previous genome wide association studies, but their biological implication remains mostly unclear. In this study, we evaluated the AIS-associated variants in the 7p22.3 locus by combining in silico, in vitro, and in vivo analyses. rs78148157 was located in an enhancer of UNCX, a homeobox gene and its risk allele upregulated the UNCX expression. A transcription factor, early growth response 1 (EGR1), transactivated the rs78148157-located enhancer and showed a higher binding affinity for the risk allele of rs78148157. Furthermore, zebrafish larvae with UNCX messenger RNA (mRNA) injection developed body curvature and defective neurogenesis in a dose-dependent manner. rs78148157 confers the genetic susceptibility to AIS by enhancing the EGR1-regulated UNCX expression. © 2022 American Society for Bone and Mineral Research (ASBMR).

Where is bone science taking us?

Bone science has over the last decades unraveled many important pathways in bone and mineral metabolism and the interplay between genetic factors and the environment. Some of these discoveries have led to the development of pharmacological treatments of osteoporosis and rare bone diseases. Other scientific avenues have uncovered a role for the gut microbiome in regulating bone mass, which have led to investigations on the possible therapeutic role of probiotics in the prevention of osteoporosis. Huge advances have been made in identifying the genes that cause rare bone diseases, which in some cases have led to therapeutic interventions. Advances have also been made in understanding the genetic basis of the more common polygenic bone diseases, including osteoporosis and Paget's disease of bone (PDB). Polygenic profiles are used for establishing genetic risk scores aiming at early diagnosis and intervention, but also in Mendelian randomization (MR) studies to investigate both desired and undesired effects of targets for drug design.

Biological effect of dysregulated LBX1 on adolescent idiopathic scoliosis through modulating muscle carbohydrate metabolism

BACKGROUND CONTEXT

Abnormal energy metabolism such as lower body weight and body mass index (BMI) and less fat mass is widely reported in patients with adolescent idiopathic scoliosis (AIS) and has been implicated in deformity development. However, the underlying mechanism is largely unclear. LBX1 is one of the promising AIS predisposing genes validated by multicenter studies.

PURPOSE

This study aimed to identify differentially expressed proteins (DEPs) relating to energy metabolism in AIS by using proteomic and metabolic analysis and to explore if the expression of these DEPs is associated with clinical parameters and modulated by LBX1.

STUDY DESIGN

This is a cross-sectional study using clinical data and biological samples followed by basic study using a cellular model.

PATIENT SAMPLE

Plasma samples were collected from Chinese girls with nonprogressive and progressive AIS (N=7 and 8, respectively) and age-matched healthy girls (N=50). Paraspinal muscle tissues were collected intraoperatively from concave and convex side of the apex of the major spinal curve in AIS (N=24) and either side from nonscoliosis patients (N=14).

OUTCOME MEASURES

Radiological Cobb angle and basic anthropometric data of recruited subjects were measured. The DEPs and metabolites were compared in plasma using proteomics and metabolomics technique. The relative expression of selected genes was measured in muscles.

METHODS

Plasma samples from AIS were collected at first clinical visit and were further divided into nonprogressive or progressive groups according to Cobb angle changes in 6-year follow-up. Age-matched healthy girls were recruited as control. High-performance liquid chromatography-mass spectrometry based proteomic analysis was carried out in three groups to identify DEPs and their annotated metabolic pathways. An independent cohort was used for validation by gas chromatography-mass spectrometry based metabolomic analysis. Paraspinal muscles were subjected to quantitative polymerase chain reaction (qPCR) followed by correlation analysis. Human skeletal muscle myoblast (HSMM) was used as the cellular model.

RESULTS

The likelihood of aberrant galactose metabolism and glycolysis was found to be associated with AIS curve progression as evidenced by the thirteen DEPs and seven related metabolites according to proteomic and metabolomic analysis. Some of the DEPs showed significantly altered expression in AIS concave and convex sides paraspinal muscles compared with those in nonscoliosis control. Four DEPs were found significantly and negatively correlated with LBX1 in AIS convex side paraspinal muscles. Overexpressing LBX1 in HSMM cells led to increased expression of three DEPs and decreased expression of three DEPs, respectively.

CONCLUSIONS

This is the first integrated proteomic and metabolomic analysis on AIS. Our findings show dysregulated galactose metabolism and glycolysis pathways in progressive group of AIS, suggesting the presence of abnormal energy metabolism at early stage of this disease, and their association with higher risk of progressing into more severe curvature. Evidence from ex vivo study with human muscle biopsies and in vitro study with human myoblast cells propose the possible effect of LBX1 on these two pathways in skeletal muscles. The present study provides new evidence of LBX1 function in AIS via modulating effect on the expression of energy metabolism related genes. This study might provide new insights into etiopathogenesis and development of novel treatment strategy targeting on abnormal body weight and BMI in patients with AIS. Additionally, the plasma proteomic and metabolomic studies suggested new candidates as biomarkers for establishing predictive model for AIS onset/progression.

Association of LBX1 Gene Methylation Level with Disease Severity in Patients with Idiopathic Scoliosis: Study on Deep Paravertebral Muscles

Idiopathic scoliosis (IS) is a multifactorial disease with a genetic background. The association of Ladybird Homeobox 1 (LBX1) polymorphisms with IS has been proven in multiple studies. However, the epigenetic mechanisms have not been evaluated. This study aimed to evaluate the LBX1 methylation level in deep paravertebral muscles in order to analyze its association with IS occurrence and/or IS severity. Fifty-seven IS patients and twenty non-IS patients were examined for the paravertebral muscles’ methylation level of the LBX1 promoter region. There was no significant difference in methylation level within paravertebral muscles between patients vs. controls, except for one CpG site. The comparison of the paravertebral muscles’ LBX1 promoter region methylation level between patients with a major curve angle of ≤70° vs. >70° revealed significantly higher methylation levels in 17 of 23 analyzed CpG sequences at the convex side of the curvature in patients with a major curve angle of >70° for the reverse strand promoter region. The association between LBX1 promoter methylation and IS severity was demonstrated. In patients with severe IS, the deep paravertebral muscles show an asymmetric LBX1 promoter region methylation level, higher at the convex scoliosis side, which reveals the role of locally acting factors in IS progression.

Abnormal TNS3 gene methylation in patients with congenital scoliosis

Background

Congenital scoliosis (CS) is a congenital deformity of the spine resulting from abnormal and asymmetrical development of vertebral bodies during pregnancy. However, the etiology and mechanism of CS remain unclear. Epigenetics is the study of heritable variations in gene expression outside of changes in nucleotide sequence. Among these, DNA methylation was described first and is the most characteristic and most stable epigenetic mechanism. Therefore, in this study, we aim to explore the association between genome methylation and CS which are not been studied before.

Methods

Two pairs of monozygotic twins were included, with each pair involving one individual with and one without CS. Agilent SureSelect XT Human Methyl-Sequencing was used for genome methylation sequencing. MethylTarget was used to detect methylation levels in target regions. Immunohistochemistry was performed to visualize expression of associated genes in candidate regions.

Results

A total of 75 differentially methylated regions were identified, including 24 with an increased methylation level and 51 with a decreased methylation level in the CS group. Nine of the differentially methylated regions were selected (TNS3, SEMAC3, GPR124, MEST, DLK1, SNTG1, PPIB, DEF8, and GRHL2). The results showed that the methylation level of the promoter region of TNS3 was 0.72 ± 0.08 in the CS group and 0.43 ± 0.06 in the control group (p = 0.00070 < 0.01). There was no significant difference in the degree of methylation of SEMAC3, GPR124, MEST, DLK1, SNTG1, PPIB, DEF8, or GRHL2 between the two groups. Immunohistochemistry showed significantly decreased TNS3 expression in the cartilage of the articular process in CS (CS: 0.011 ± 0.002; control: 0.018 ± 0.006, P = 0.003 < 0.01).

Conclusion

Compared with the control group, high-level methylation of the TNS3 promoter region and low TNS3 expression in the cartilage layer of the articular process characterize CS. Thus, DNA methylation and TNS3 may play important roles in the pathogenesis of CS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12891-022-05730-x.

Adolescent idiopathic scoliosis : a review of aetiological theories of a multifactorial disease

Adolescent idiopathic scoliosis (AIS), defined by an age at presentation of 11 to 18 years, has a prevalence of 0.47% and accounts for approximately 90% of all cases of idiopathic scoliosis. Despite decades of research, the exact aetiology of AIS remains unknown. It is becoming evident that it is the result of a complex interplay of genetic, internal, and environmental factors. It has been hypothesized that genetic variants act as the initial trigger that allow epigenetic factors to propagate AIS, which could also explain the wide phenotypic variation in the presentation of the disorder. A better understanding of the underlying aetiological mechanisms could help to establish the diagnosis earlier and allow a more accurate prediction of deformity progression. This, in turn, would prompt imaging and therapeutic intervention at the appropriate time, thereby achieving the best clinical outcome for this group of patients. Cite this article: Bone Joint J 2022;104-B(8):915-921.

Upregulation of microRNA-96-5p is associated with adolescent idiopathic scoliosis and low bone mass phenotype

Bone densitometry revealed low bone mass in patients with adolescent idiopathic scoliosis (AIS) and its prognostic potential to predict curve progression. Recent studies showed differential circulating miRNAs in AIS but their diagnostic potential and links to low bone mass have not been well-documented. The present study aimed to compare miRNA profiles in bone tissues collected from AIS and non-scoliotic subjects, and to explore if the selected miRNA candidates could be useful diagnostic biomarkers for AIS. Microarray analysis identified miR-96-5p being the most upregulated among the candidates. miR-96-5p level was measured in plasma samples from 100 AIS and 52 healthy girls. Our results showed significantly higher plasma levels of miR-96-5p in AIS girls with an area under the curve (AUC) of 0.671 for diagnostic accuracy. A model that was composed of plasma miR-96-5p and patient-specific parameters (age, body weight and years since menarche) gave rise to an improved AUC of 0.752. Ingenuity Pathway Analysis (IPA) indicated functional links between bone metabolic pathways and miR-96-5p. In conclusion, differentially expressed miRNAs in AIS bone and plasma samples represented a new source of disease biomarkers and players in AIS etiopathogenesis, which required further validation study involving AIS patients of both genders with long-term follow-up.

Epigenetic and Genetic Factors Related to Curve Progression in Adolescent Idiopathic Scoliosis: A Systematic Scoping Review of the Current Literature

Adolescent idiopathic scoliosis (AIS) is a progressive deformity of the spine. Scoliotic curves progress until skeletal maturity leading, in rare cases, to a severe deformity. While the Cobb angle is a straightforward tool in initial curve magnitude measurement, assessing the risk of curve progression at the time of diagnosis may be more challenging. Epigenetic and genetic markers are potential prognostic tools to predict curve progression. The aim of this study is to review the available literature regarding the epigenetic and genetic factors associated with the risk of AIS curve progression. This review was carried out in accordance with Preferential Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) guidelines. The search was carried out in January 2022. Only peer-reviewed articles were considered for inclusion. Forty studies were included; fifteen genes were reported as having SNPs with significant association with progressive AIS, but none showed sufficient power to sustain clinical applications. In contrast, nine studies reporting epigenetic modifications showed promising results in terms of reliable markers. Prognostic testing for AIS has the potential to significantly modify disease management. Most recent evidence suggests epigenetics as a more promising field for the identification of factors associated with AIS progression, offering a rationale for further investigation in this field.

Polygenic Risk Score of Adolescent Idiopathic Scoliosis for Potential Clinical Use

Adolescent idiopathic scoliosis (AIS) is a common disease causing three-dimensional spinal deformity in as many as 3% of adolescents. Development of a method that can accurately predict the onset and progression of AIS is an immediate need for clinical practice. Because the heritability of AIS is estimated as high as 87.5% in twin studies, prediction of its onset and progression based on genetic data is a promising option. We show the usefulness of polygenic risk score (PRS) for the prediction of onset and progression of AIS. We used AIS genomewide association study (GWAS) data comprising 79,211 subjects in three cohorts and constructed a PRS based on association statistics in a discovery set including 31,999 female subjects. After calibration using a validation data set, we applied the PRS to a test data set. By integrating functional annotations showing heritability enrichment in the selection of variants, the PRS demonstrated an association with AIS susceptibility (p = 3.5 × 10^-40 with area under the receiver-operating characteristic [AUROC] = 0.674, sensitivity = 0.644, and specificity = 0.622). The decile with the highest PRS showed an odds ratio of as high as 3.36 (p = 1.4 × 10^-10 ) to develop AIS compared with the fifth in decile. The addition of a predictive model with only a single clinical parameter (body mass index) improved predictive ability for development of AIS (AUROC = 0.722, net reclassification improvement [NRI] 0.505 ± 0.054, p = 1.6 × 10^-8 ), potentiating clinical use of the prediction model. Furthermore, we found the Cobb angle (CA), the severity measurement of AIS, to be a polygenic trait that showed a significant genetic correlation with AIS susceptibility (rg = 0.6, p = 3.0 × 10^-4 ). The AIS PRS demonstrated a significant association with CA. These results indicate a shared polygenic architecture between onset and progression of AIS and the potential usefulness of PRS in clinical settings as a predictor to promote early intervention of AIS and avoid invasive surgery. © 2021 American Society for Bone and Mineral Research (ASBMR).

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.

Whole-genome methylation analysis reveals novel epigenetic perturbations of congenital scoliosis

Congenital scoliosis (CS) is a congenital disease caused by malformations of vertebrae. Recent studies demonstrated that DNA modification could contribute to the pathogenesis of disease. This study aims to identify epigenetic perturbations that may contribute to the pathogenesis of CS. Four CS patients with hemivertebra were enrolled and underwent spine correction operations. DNA was extracted from the hemivertebrae and spinal process collected from the specimen during the hemivertebra resection. Genome-wide DNA methylation profiling was examined at base-pair resolution using whole-genome bisulfite sequencing (WGBS). We identified 343 genes with hyper-differentially methylated regions (DMRs) and 222 genes with hypo-DMRs, respectively. These genes were enriched in the mitogen-activated protein kinase (MAPK) signaling pathway, calcium signaling pathway, and axon guidance in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and were enriched in positive regulation of cell morphogenesis involved in differentiation, regulation of cell morphogenesis involved in differentiation, and regulation of neuron projection development in Biological Process of Gene Ontology (GO-BP) terms. Hyper-DMR-related genes, including IGHG1, IGHM, IGHG3, RNF213, and GSE1, and hypo DMR-related genes, including SORCS2, COL5A1, GRID1, RGS3, and ROBO2, may contribute to the pathogenesis of hemivertebra. The aberrant DNA methylation may be associated with the formation of hemivertebra and congenital scoliosis.

The Susceptibility and Potential Functions of the LBX1 Gene in Adolescent Idiopathic Scoliosis

Genome-wide association studies have identified many susceptibility genes for adolescent idiopathic scoliosis (AIS). However, most of the results are hard to be replicated in multi-ethnic populations. LBX1 is the most promising candidate gene in the etiology of AIS. We aimed to appraise the literature for the association of LBX1 gene polymorphisms with susceptibility and curve progression in AIS. We also reviewed the function of the LBX1 gene in muscle progenitor cell migration and neuronal determination processes. Three susceptibility loci (rs11190870, rs625039, and rs11598564) near the LBX1 gene, as well as another susceptibility locus (rs678741), related to LBX1 regulation, have been successfully verified to have robust associations with AIS in multi-ethnic populations. The LBX1 gene plays an essential role in regulating the migration and proliferation of muscle precursor cells, and it is known to play a role in neuronal determination processes, especially for the fate of somatosensory relay neurons. The LBX1 gene is the most promising candidate gene in AIS susceptibility due to its position and possible functions in muscle progenitor cell migration and neuronal determination processes. The causality between susceptibility loci related to the LBX1 gene and the pathogenesis of AIS deserves to be explored with further integrated genome-wide and epigenome-wide association studies.

Methylation of estrogen receptor 2 (ESR2) in deep paravertebral muscles and its association with idiopathic scoliosis

Idiopathic scoliosis (IS) is one of the most common spinal disorders in adolescents. Despite many studies, the etiopathogenesis of IS is still poorly understood. In recent years, the role of epigenetic factors in the etiopathogenesis of IS has been increasingly investigated. It has also been postulated that the development and progression of the disease is related to gender and puberty, and could be associated with estrogen action. Estrogen hormones act via estrogen receptor 1 (ESR1) and estrogen receptor 2 (ESR2). It has been suggested that ESR2 expression is dependent on methylation within its gene promoter. So far, no studies have evaluated local, tissue-specific DNA methylation in patients with IS. Thus, our study aimed to analyze the methylation and expression level of ESR2 in the paraspinal muscles of the convex and concave side of the IS curvature. The methylation level within ESR2 promoter 0N, but not exon 0N, was significantly higher on the concave side of the curvature compared to the convex side. There was no significant correlation between ESR2 expression and methylation level in the promoter 0N on the convexity of thoracic scoliosis, whereas, on the concave side of the curvature, we observed a moderate negative correlation. There was no difference in the methylation levels of the ESR2 promoter and exon 0N between groups of patients with Cobb angle ≤ 70° and > 70° on the concave and convex side of the curvature. We also found no statistically significant correlation between the Cobb angle value and the mean methylation level in either the ESR2 promoter or exon 0N on the convex or concave side of the curvature. Our findings demonstrate that DNA methylation at the ESR2 promoter in deep paravertebral muscle tissue is associated with the occurrence but not with the severity of idiopathic scoliosis.

Research progress on the etiology and pathogenesis of adolescent idiopathic scoliosis

Etiology of adolescent idiopathic scoliosis (AIS), a complicated three-dimensional spinal deformity with early-onset, receives continuous attention but remains unclear. To gain an insight into AIS pathogenesis, this review searched PubMed database up to June 2019, using key words or medical subject headings terms including "adolescent idiopathic scoliosis," "scoliosis," "pathogenesis," "etiology," "genetics," "mesenchymal stem cells," and their combinations, summarized existing literatures and categorized the theories or hypothesis into nine aspects. These aspects include bone marrow mesenchymal stem cell studies, genetic studies, tissue analysis, spine biomechanics measurements, neurologic analysis, hormone studies, biochemical analysis, environmental factor analysis, and lifestyle explorations. These categories could be a guidance for further etiology or treatment researches to gain inspiration.

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.

DNA methylation biomarkers of future health outcomes in children

Biomarkers which predict future health outcomes are key to the goals of precision health. Such biomarkers do not have to be involved in the causal pathway of a disease, and their performance is best assessed using statistical tests of clinical performance and evaluation of net health impact. DNA methylation is the most commonly studied epigenetic process and represents a potential biomarker of future health outcomes. We review 25 studies in non-oncological paediatric conditions where DNA methylation biomarkers of future health outcomes are assessed. Whilst a number of positive findings have been described, the body of evidence is severely limited by issues with outcome measures, tissue-specific samples, accounting for sample cell type heterogeneity, lack of appropriate statistical testing, small effect sizes, limited validation, and no assessment of net health impact. Future studies should concentrate on careful study design to overcome these issues, and integration of DNA methylation data with other 'omic', clinical, and environmental data to generate the most clinically useful biomarkers of paediatric disease.

High methylation of lysine acetyltransferase 6B is associated with the Cobb angle in patients with congenital scoliosis

Background

The etiology of congenital scoliosis (CS) is complex and uncertain. Abnormal DNA methylation affects the growth and development of spinal development. In this study, we investigated the role of DNA methylation in CS.

Methods

The target region DNA methylation level in the peripheral blood of patients with CS was analyzed. Through in-depth analysis, genes closely related to the growth and development of the vertebra were identified. EdU staining was performed to verify the role of differentially expressed genes in chondrocyte proliferation.

Results

The hypermethylated KAT6B gene was observed in patients with CS, and was positively correlated with the Cobb angle. KAT6B was primarily expressed on chondrocytes. The promoter of KAT6B in CS patients was hypermethylated, and its expression was significantly reduced. Further mechanistic studies revealed that EZH2 mediated trimethylation of lysine 27 on histone H3 of the KAT6B promoter. Overexpression of KAT6B in CS-derived primary chondrocytes can significantly promote chondrocyte proliferation, which may be related to activation of the RUNX2/Wnt/β-catenin signaling pathway.

Conclusion

Epigenetic modification of KAT6B may be a cause of CS. If similar epigenetic modification abnormalities can be detected through maternal liquid biopsy screening, they may provide useful biomarkers for early screening and diagnosis of CS.

Exome sequencing analysis identifies frequent oligogenic involvement and FLNB variants in adolescent idiopathic scoliosis

Background

Adolescent idiopathic scoliosis (AIS) is a genetically heterogeneous disease characterised by three-dimensional deformity of the spine in the absence of a congenital spinal anomaly or neurological musculoskeletal disorder. The clinical variability and incomplete penetrance of some genes linked with AIS indicate that this disease constitutes an oligogenic trait.

Objective

We aimed to explore the oligogenic nature of this disease and identify novel AIS genes.

Methods

We analysed rare damaging variants within AIS-associated genes by using exome sequencing in 40 AIS trios and 183 sporadic patients.

Results

Multiple variants within AIS-associated genes were identified in eight AIS trios, and five individuals harboured rare damaging variants in the FLNB gene. The patients showed more frequent oligogenicity than the controls. In the gene-based burden test, the top signal resided in FLNB. In functional studies, we found that the AIS-associated FLNB variants altered the protein’s conformation and subcellular localisation and its interaction with other proteins (TTC26 and OFD1) involved in AIS. The most compelling evidence of an oligogenic basis was that the number of rare damaging variants was recognised as an independent prognostic factor for curve progression in Cox regression analysis.

Conclusion

Our data indicate that AIS is an oligogenic disease and identify FLNB as a susceptibility gene for AIS.

CHD7 gene polymorphisms in female patients with idiopathic scoliosis

BACKGROUND

The CHD7 (chromosome domain helicase DNA binding protein 7) gene has been associated with familial idiopathic scoliosis (IS) in families of European descent. The CHD7 single-nucleotide polymorphisms have never been studied in Polish Caucasian IS patients.

METHODS

The aim of this study was to investigate the relationship of CHD7 gene polymorphisms with susceptibility to or progression of IS in Polish Caucasian females. The study group comprised 211 females who underwent clinical, radiological and genetic examination. The study group was analyzed in three subgroups according to: (1) Cobb angle (Cobb angle ≤30° vs. Cobb angle ≥35°), (2) age of diagnosis (adolescent IS vs. early-onset IS) and (3) rate of progression (non-progressive vs. slowly progressive vs. rapidly progressive IS). The control group comprised 83 females with no scoliosis and with a negative family history who underwent clinical and genetic examination. In total six CHD7 gene polymorphisms were examined. Three polymorphisms (rs1017861, rs13248429, and rs4738813) were examined by RFLP (restriction fragment length polymorphism) analysis, and three were quantified by Sanger sequencing (rs78874766, rs4738824, and rs74797613).

RESULTS

In rs13248429, rs78874766, and rs74797613 polymorphisms only the wild allele was present. The rs1017861 polymorphism demonstrated an association with IS susceptibility (p < 0.01). Two polymorphisms, rs1017861 and rs4738813, were associated with curve severity and progression rate (p < 0.05). None of the evaluated polymorphisms in CHD7 gene showed any association with the age of IS onset.

CONCLUSIONS

The polymorphism rs1017861 in CHD7 gene showed an association with IS susceptibility. Two polymorphisms (rs1017861 and rs4738813) were associated with curve severity and progression rate. None of the evaluated polymorphisms in CHD7 gene showed any association with the age of IS onset. Further evaluation of CHD7 gene should be considered as IS modifying factor.

A validated composite model to predict risk of curve progression in adolescent idiopathic scoliosis

BACKGROUND

In adolescent idiopathic scoliosis (AIS), the continuous search for effective prognostication of significant curve progression at the initial clinical consultation to inform decision for timely treatment and to avoid unnecessary overtreatment remains a big challenge as evidence of the multifactorial etiopathogenic nature is increasingly reported. This study aimed to formulate a composite model composed of clinical parameters and circulating markers in the prediction of curve progression.

METHOD

This is a two-phase study consisting of an exploration cohort (120 AIS, mean Cobb angle of 25°± 8.5 at their first clinical visit) and a validation cohort (51 AIS, mean Cobb angle of 23° ± 5.0° at the first visit). Patients with AIS were followed-up for a minimum of six years to formulate a composite model for prediction. At the first visit, clinical parameters were collected from routine clinical practice, and circulating markers were assayed from blood.

FINDING

We constructed the composite predictive model for curve progression to severe Cobb angle > 40° with a high HR of 27.9 (95% CI of 6.55 to 119.16). The area under curve of the composite model is higher than that of individual parameters used in current clinical practice. The model was validated by an independent cohort and achieved a sensitivity of 72.7% and a specificity of 90%.

INTERPRETATION

This is the first study proposing and validating a prognostic composite model consisting of clinical and circulating parameters which could quantitatively evaluate the probability of curve progression to a severe curvature in AIS at the initial consultation. Further validation in clinic will facilitate application of composite model in assisting objective clinical decision.

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.