Association of MTHFR rs1801133 and homocysteine with Legg–Calvé–Perthes disease in Mexican patients
Orphanet Journal of Rare Diseases volume 17, Article number: 123 (2022)
Legg–Calvé–Perthes disease (LCPD) is the avascular osteonecrosis of the proximal femoral epiphysis. It is a rare disease of unclear etiology in children, although alterations in coagulation or the collagen gene have been described and could be associated with its etiology. Our objective was to evaluate the following alterations: COL1A1 (rs1107946, rs2412298), COL2A1 (rs121912891 and rs387106558), MTHFR rs1801133, CBS rs115742905, and PT rs1799963 and their relationship with LCPD.
DNA was obtained and genotyped by real-time PCR with TaqMan probes. Prothrombin (FII) and homocysteine (Hcy) were determined by a coagulometric method. The variables were described as mean and standard deviation or percentages, and genotypic and allelic distributions were analyzed using the Student's t-test. The Hardy–Weinberg equilibrium and OR were also used.
We studied 23 patients with LCPD and 46 controls. We did not find any association of the MTHFR, CBS, PT, COL1A1, and COL2A1 genetic variants with LCPD. However, when adjusting the data with the Hcy values for the MTHFR C677T polymorphism, the C/C genotypes showed an association with the recessive model (p = 0.038), with susceptibility to LCPD.
No association was found with the CBS, PT, COL1A1, and COL2A1 genes. Nevertheless, our results suggest a significant link between moderately elevated Hcy levels and the MTHFR C677T polymorphism in a cohort of Mexican children with LCPD.
Due to its low prevalence, Legg–Calvé–Perthes disease (LCPD) is rare, particularly in Latin America. However, the avascular osteonecrosis of the femoral head does occur in the pediatric population [1, 2]. The interruption of the blood supply causes fragmentation and deformation of the femoral head, which leads to mechanical hip alterations . Although there are several proposed causes, its etiology is still unknown. Thrombophilia is one of the most widely accepted since the tendency to develop states of hypercoagulability could lead to the development of necrosis [4,5,6]. On the other hand, described collagen gene alterations could also be related to the etiology of the disease .
The methylenetetrahydrofolate reductase enzyme (MTHFR), together with the Cystathionine β-synthase (CBS), regulates the serum homocysteine (Hcy) level. MTHFR C677T and CBS T833C polymorphisms, Hcy levels, and other risk factors like polymorphisms of prothrombin G20210A (PT G20210A) have been linked to the development of prothrombotic states . Various studies have suggested that these conditions could be related to the etiology of LCPD [9, 10].
LCPD consists of sequential stages of involvement of the proximal femoral epiphyses, including subchondral fracture, fragmentation, reossification, and healing with residual deformity. Even though genetic factors have been linked to the etiology of LCPD, the causal gene has not been identified. Nonetheless, COL2A1 gene mutations have been implicated in the etiology of LCPD because collagen is the essential matrix protein of all connective tissues. Some of these alterations could cause structural defects of type II collagen related to the major clinical aspects [11, 12].
On the other hand, the COL1A1 gene encodes the alpha-1 protein chain of type I collagen, the primary protein of bones. Some research focused on the COL1A1 polymorphisms associates them significantly with low bone mineral density, osteoporosis, increased fracture risk, and osteonecrosis [13, 14]. Also, the harmful effect of Hcy on bones may be mainly mediated by the binding of the molecule to collagen, which interferes with the formation of type I collagen cross-links . Furthermore, several studies report a significant relationship between hypercoagulability, elevated levels of Hcy, and osteonecrosis [16, 17]. Because of the above, our objective is to evaluate the following alterations involved with coagulation genes, which could be associated with LCPD etiology: MTHFR C677T rs1801133, CBS T833C rs115742905, and PT G20210A rs1799963; the following collagen-related polymorphisms: COL1A1 G1997T rs1107946, COL1A1 C1663T rs2412298, COL2A1 G2306A rs121912891, and COL2A1 G3665A rs387106558; and their relationship with Mexican patients suffering from LCPD.
We carried out a case–control study. The cases were first-time or recurrent patients of all ages, of both sexes, with a clinical and radiological diagnosis of LCPD, without other bone diseases or diseases related to coagulation abnormalities, and without any pharmacological treatment. Additionally, we matched healthy controls with the patients in a 2:1 proportion by age and sex, with radiological studies showing that the controls did not have any alterations in the femur or hip, had no history of hematological or thrombotic pathologies, and were not under pharmacological treatment.
All patients were recruited from the Orthopedic Service at Instituto Nacional de Rehabilitación Luis Guillermo Ibarra Ibarra (INR-LGII) in Mexico City. In addition, we obtained weight (kg), height (m), and body mass index (BMI).
We obtained DNA from peripheral blood leukocytes. Genotyping was performed by real-time PCR with TaqMan probes in the conditions recommended by the manufacturer (Applied Biosystems, Foster City, CA). Briefly, reactions were performed in 25 ml volume, containing TaqMan PCR master mix (1X), TaqMan probe (100 nm), primers (900 nm of each), and DNA (25 ng). Samples were run in a StepOne Real-Time PCR System (Applied Biosystems). Cycling conditions: one cycle at 95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s, and a final cycle at 60 °C for 1 min.
Determination of coagulation factor II (prothrombin) and Hcy
Plasma was obtained from peripheral blood in tubes with 3.8% sodium citrate. The hemosiLTM 0020012200 chromogenic kit and the HemosiLTM 0020007800 kit for Hcy were used in ACL-7000 elite pro equipment to determine the prothrombin or coagulation factor II (FII). The results were expressed as a percentage of activity, reference range: 98 to 136% for FII; and in µmol/L, reference range: 10 to 14.4 for µmol/L Hcy.
The variables are described as mean ± standard deviation or percentages (%) where appropriate. Differences in genotypic and allelic distribution between case and control groups were analyzed using Fisher's exact test or Student's t-test. The Hardy–Weinberg equilibrium (HWE) was calculated for the genotype frequencies in cases and controls using a Chi-square goodness of fit test in SNPStats. The odds ratio (OR) and the confidence interval (95% CI) were calculated to assess the relative risk of having the disease according to the relative frequency of different genotypes between cases and controls with PAST 4.03. Furthermore, the inheritance patterns for SNPs and associations were estimated using SNPStats (http://www.snpstats.net/start.htm).
The cases were patients with LCPD. Controls were selected under the official Mexican standard NOM-253-SSA1-2012 guidelines, all participants received oral and written information, and parents or guardians signed a letter of consent. This study was approved by the INR-LGII research and ethics committees.
This study included 69 subjects: 23 patients (21 men and 2 women) of 3 to 18 years of age and 46 controls (42 men and 4 women). The age of all participants ranged from 3 to 18 years, with a mean (X) of 16.8 and a standard deviation (SD) of ± 11.6 (patients X = 16.8 ± 11.3, and controls X = 16.3 ± 11.0). The average height in patients was 1.48 ± 0.23 m compared to 1.49 ± 0.23 m in controls. The average weight was 46.21 ± 15.02 kg in patients, compared to 53.83 ± 24.06 kg in controls, without significant difference between patients and controls. Similarly, BMI presented no significant differences between patients and controls.
Prothrombin activity did not significantly differ between the patients' mean and the controls. However, regarding the levels of Hcy, there were substantial differences in the mean of the patients compared to the controls (Fig. 1).
The genotype frequencies of the MTHFR rs1801133 polymorphism were in Hardy–Weinberg equilibrium (HWE) because the mutant allele is present in controls and patients (p = 0.76 and p = 0.40, respectively). The CBS rs115742905 and PT rs1799963 variants were not found in HWE because the mutant allele is absent in either case. The wild genotype resulted in 100% in both populations (Table 1). The approach to the disease is complicated because its etiology is not yet clear.
The genotype frequency of the COL1A1 rs1107946 variant was not found in HWE with respect to the controls. No significant differences were observed for the genotype and allelic distribution for the polymorphism between the patients with LCPD and the controls. The COL1A1 variant rs24122948 was found in HWE, but no differences were observed. In the case of COL2A1 rs24122948 and rs106558, we did not find the presence of the mutant alleles (Table 1).
Table 1 presents the link of MTHFR to age, which showed a significant association with the recessive inheritance model and was detected among carriers of the MTHFR CC compared to carriers of MTHFR TT + CT, regardless of age. However, apparently, a worse outcome is more closely related to an older age.
The prevalence of LCPD is highly variable in different populations. Since it is a rare disease, few reports in Mexico and Latin America describe it [18,19,20]. Children with LCPD begin showing pain in one or both legs and a slight limp, sometimes attributed to falls. On occasion, the pain goes away temporarily, but the lameness continues. Generally, the disease encompasses a broad spectrum of symptoms, from mild with no long-term sequelae to severe with a permanent degenerative change of the hip joint. Obtaining a timely and accurate diagnosis for a rare disease like LCPD can often be challenging due to the lack of experienced healthcare providers. In addition, research activities are less common, and, in developing countries, the problems are compounded by other resource limitations [21, 22]. Therefore, the correct diagnosis is only given when patients finally go to specialized pediatric orthopedic institutions, such as the INR-LGII.
The approach to the disease is complicated because its etiology is not yet clear. However, the main priority is to make an early diagnosis and offer a conservative or surgical treatment that keeps the head of the femur and the hip joint in the best conditions. For this reason, it is vital that primary care physicians know and understand this disease.
Genes involved in bone remodeling and architecture have been studied, and elevated levels of Hcy have been linked to osteonecrosis . Because of this, we decided to explore variants in the COL1A1 gene that have been studied in relation to osteoporosis and fractures [23, 24]. However, none of these variants was informative in the studied sample of Mexican patients with LCPD. It has also been suggested that COL2A1, involved in osteonecrosis of the femoral head, is related to skeletal dysplasia due to failed cartilage development and growth . However, none of these alterations had significant differences between the patients and the controls studied. Additionally, COL2A1 is suggested to be one of the representative causal genes, as Li et al. found a mutation in COL2A1 G1888A in a Chinese family affected by LCPD and osteonecrosis of the femoral head [12, 26], and Miyamoto et al. identified a mutation in COL2A1 G3508A in a Japanese family with LCPD. However, these findings differ from our results since we found no relationship with LCPD .
The PT G20210A mutation increases the risk of thrombosis . López et al. studied 90 children with LCPD, and no patient had a family or personal history of early thrombotic events. However, compared to controls, four children with LCPD (4.4%) were heterozygous for the G20210A polymorphism but without significant association . Vosmaer et al. found an increase in the incidence of LCPD in the presence of the prothrombin mutation when the FVL polymorphism was also present . Our study did not find significant differences that relate this polymorphism, PT rs1799963, with LCPD.
Common alterations in some genes related to Hcy metabolism, such as MTHFR C677T and CBS T833C, have been shown to cause increased plasma homocysteine levels, thus bringing about a predisposition to thrombosis . Our study did not find these variants, CBS T833C rs5742905, neither in controls nor in patients. Azarpira et al. showed that the MTHFR C677T polymorphism was not associated with the risk of LCPD in Iranian children . When we analyzed the results for the MTHFR C677T rs1801133 polymorphism, they did not present a significant association with the risk of susceptibility to LCPD. However, the mutant allele frequency was higher than in the Iranian population, even though several epidemiological case–control studies have found that MTHFR polymorphisms could not play a significant role in the susceptibility to the development of osteonecrosis on the femoral head [33, 34]. In this study, we also found that the levels of Hcy between cases and controls present a significant difference (p < 0.05). In addition, when adjusting the value of the MTHFR C677T polymorphism with age, we found a significant association in the recessive model (OR 5.62 (1.02–30.93), p = 0.038). Age is one of the main factors that can influence a poor prognosis for patients.
In general, hyperhomocysteinemia is a risk factor for various diseases: vascular, neurological, diabetes, psoriasis, cancers, osteoporosis, among others . In addition, Hcy is thought to alter bone remodeling . Furthermore, it has been described that short-term moderate hyperhomocysteinemia affects bone and cartilage characteristics. Trabecular bone microarchitecture is especially sensitive to hyperhomocysteinemia. It shows clearly negative bone balance, mainly due to a decrease in trabecular numbers and markedly reduced trabecular connections, which indicates multiple alterations in collagen due to homocysteine accumulation in bones, indicative of broken collagenous cross-links .
One of the weaknesses of this work was that it required considerable time to capture this small population of patients, which poses difficulties in research. However, the results strengthen and support the theories that describe thrombosis and homocysteine as risk factors in developing LCPD. In addition, it contributes to the scientific knowledge and information about the disease in the Mexican population.
This is the first study to evaluate genetic aspects of coagulation mediators and type I and II collagen in Mexican patients with LCPD. We did not find a relationship of collagen alterations (COL1A1 rs1107946, COL1A1 rs2412298, COL2A1 rs121912891, and COL2A1 rs387106558) or coagulation alterations (CBS T833C rs115742905 and PT G20210A rs1799963) with LCPD. However, we found an association of the MTHFR rs1801133 polymorphism in a recessive model adjusted by age and elevated levels of Hcy in a population of Mexican children with LCPD.
Availability of data and materials
All relevant data used in this study are included in the manuscript. The corresponding author can be contacted if any further information is needed.
Collagen type I alpha 1 chain
Collagen type II alpha 1 chain
Instituto Nacional de Rehabilitación “Luis Guillermo Ibarra Ibarra”
Body mass index
Coagulation factor II
Single nucleotide polymorphisms
Norma Oficial Mexicana
Jaramillo D. What is the optimal imaging of osteonecrosis, Perthes, and bone infarcts? Pediatr Radiol. 2009;39:216.
Srzentić S, Nikčević G, Spasovski D, Baščarević Z, Živković Z, Terzic-Šupić Z, Matanović D, Djordjević V, Pavlović S, Spasovski V. Predictive genetic markers of coagulation, inflammation, and apoptosis in Perthes disease—Serbian experience. Eur J Pediatr. 2015;174:1085–92.
Chen L, Hong G, Zhang Q, Chen Z, Fang B, He W. Waldenström classification: better evaluation for the intervention time of Scottish rite brace treatment in Legg–Calvé–Perthes disease. Int J Clin Exp Med. 2017;10:724–31.
Kealey WD, Mayne EE, McDonald W, Murray P, Cosgrove AP. The role of coagulation abnormalities in the development of Perthes’ disease. J Bone Jt Surg Br. 2000;82:744–6.
Arruda VR, Belangero WD, Ozelo MC, Oliveira GB, Pagnano RG, Volpon JB, Annichino-Bizzacchi JM. Inherited risk factors for thrombophilia among children with Legg–Calvé–Perthes disease. J Pediatr Orthop. 1999;19:84–7.
Perry DC, Hall AJ. The epidemiology and etiology of Perthes disease. Orthop Clin North Am. 2011;42:279–83.
Poort SR, Rosendaal FR, Reitsma PH, Bertina RM. A common genetic variation in the 3′-untranslated region of the prothrombin gene is associated with elevated plasma prothrombin levels and an increase in venous thrombosis. Blood. 1996;88:3698–703.
Zavala-Hernández C, Hernández-Zamora E, Martínez-Murillo C, Majluf-Cruz A, Vela-Ojeda J, García-Chávez J, Reyes-Maldonado E. Risk factors for thrombosis development in Mexican patients. Ann Vasc Surg. 2015;29:1625–32.
Loder RT, Skopelja EN. The epidemiology and demographics of Legg–Calvé–Perthes disease. ISRN Orthop. 2011;2011:504393.
Chaudhry S, Phillips D, Feldman D. Legg–Calvé–Perthes disease: an overview with recent literature. Bull Hosp Jt Dis. 2014;72(1):18–27.
Kannu P, Bateman J, Savarirayan R. Clinical phenotypes associated with type II collagen mutations. J Paediatr Child Health. 2012;48:E38–43.
Li N, Yu J, Cao X, Wu QY, Li WW, Li TF, Zhang C, Cui YX, Li XJ, Yin ZM, Xia XY. A novel p.Gly630Ser mutation of COL2A1 in a Chinese family with presentations of Legg–Calve–Perthes disease or avascular necrosis of the femoral head. PLoS ONE. 2014;9:e100505.
Urreizti R, Garcia-Giralt N, Riancho JA, González-Macías J, Civit S, Güerri R, Yoskovitz G, Sarrion P, Mellivobsky L, Díez-Pérez A, Nogués X, Balcells S, Grinberg D. COL1A1 haplotypes and hip fracture. J Bone Miner Res. 2012;27:950–3.
Erdem M, Tüfekçi Ö, Kızıldağ S, Yılmaz Ş, Kızmazoğlu D, Eroğlu Filibeli B, Ören H. Investigation of the relationship between Fok1 and Col1A1 gene polymorphisms and development of treatment-related bone complications in children with acute lymphoblastic leukemia. Turk J Haematol. 2019;36:12–8.
Jackson SH. The reaction of homocysteine with aldehyde: an explanation of the collagen defects in homocystinuria. Clin Chim Acta. 1973;45:215–7.
Azarpira MR, Ghilian MM, Sobhan MR, Mehdinezhad-Yazdi M, Aghili K, Miresmaeili SM, Neamatzadeh H. Association of MTHFR and TNF-α genes polymorphisms with susceptibility to Legg–Calve–Perthes disease in Iranian children: a case-control study. J Orthop. 2018;15:984–7.
Glueck CJ, Freiberg RA, Fontaine RN, Tracy T, Wang P. Hypofibrinolysis, thrombophilia, osteonecrosis. Clin Orthop Relat Res. 2001;386:19–33.
Nava Zarco NA, Mijares Morán J, Hernández MA. Delay in skeletal maturity in Mexican patients with Legg–Calvé–Perthes. Acta Ortop Mex. 2008;22:157–61.
Castañeda P, Haynes R, Mijares J, Quevedo H, Cassis N. Varus-producing osteotomy for patients with lateral pillar type B and C Legg–Calvé–Perthes disease followed to skeletal maturity. J Child Orthop. 2008;2:373–9.
Sarassa CA, Herrera AM, Carvajal J, Gomez LF, Lopez CA, Rojas AF. Early clinical and radiological outcomes after double osteotomy in patients with late presentation Legg–Calvé–Perthes disease. J Child Orthop. 2008;2:425–9.
Divi SN, Bielski RJ. Legg–Calvé–Perthes disease. Pediatr Ann. 2016;45:e144–9.
Forman J, Taruscio D, Llera VA, Barrera LA, Coté TR, Edfjäll C, Gavhed D, Haffner ME, Nishimura Y, Posada M, Tambuyzer E, Groft SC, Henter JI, International Conference for Rare Diseases and Orphan Drugs (ICORD). The need for worldwide policy and action plans for rare diseases. Acta Paediatr. 2012;101:805–7.
Moradifard S, Hoseinbeyki M, Emam MM, Parchiniparchin F, Ebrahimi-Rad M. Association of the Sp1 binding site and-1997 promoter variations in COL1A1 with osteoporosis risk: the application of meta-analysis and bioinformatics approaches offers a new perspective for future research. Mutat Res. 2020;786:108339.
Jin H, Evangelou E, Ioannidis JP, Ralston SH. Polymorphisms in the 5′ flank of COL1A1 gene and osteoporosis: a meta-analysis of published studies. Osteoporos Int. 2011;22:911–21.
Sakamoto Y, Yamamoto T, Miyake N, Matsumoto N, Iida A, Nakashima Y, Research Committee on Idiopathic Osteonecrosis of the Femoral Head of the Ministry of Health, Labour and Welfare of Japan, Iwamoto Y, Ikegawa S. Screening of the COL2A1 mutation in idiopathic osteonecrosis of the femoral head. J Orthop Res. 2017;35:768–74.
Al-Omran AK, Sadat-Ali M. Legg–Calvé–Perthes disease in two generations of male family members: a case report. J Orthop Surg (Hong Kong). 2013;21:258–61.
Miyamoto Y, Matsuda T, Kitoh H, Haga N, Ohashi H, Nishimura G, Ikegawa S. A recurrent mutation in type II collagen gene causes Legg–Calvé–Perthes disease in a Japanese family. Hum Genet. 2007;121:625–9.
Carter AM, Sachchithananthan M, Stasinopoulos S, Maurer F, Medcalf RL. Prothrombin G20210A is a bifunctional gene polymorphism. Thromb Haemost. 2002;87:846–53.
López-Franco M, González-Morán G, De Lucas JC, de Velasco LPJrJF, Vivancos JC, Epeldegui-Torre T. Legg–Perthes disease and heritable thrombophilia. J Pediatr Orthop. 2005;25:456–9.
Vosmaer A, Pereira RR, Koenderman JS, Rosendaal FR, Cannegieter SC. Coagulation abnormalities in Legg–Calvé–Perthes disease. J Bone Jt Surg Am. 2010;92:121–8.
Akar N, Akar E, Misirlioğlu M, Avcu F, Yalçin A, Cin S. Search for genetic factors favoring thrombosis in Turkish population. Thromb Res. 1998;92:79–82.
Azarpira MR, Ghilian MM, Sobhan MR, Mahdinezhad-Yazdi M, Aghili K, Ahrar H, Neamatzadeh H. Association of eNOS 27-bp VNTR, 894G>T and 786T>C polymorphisms with susceptibility to Legg–Calve–Perthes disease in Iranian children. J Orthop. 2019;16:137–40.
Chai W, Zhang Z, Ni M, Geng P, Lian Z, Zhang G, Shi LL, Chen J. Genetic association between methylenetetrahydrofolate reductase gene polymorphism and risk of osteonecrosis of the femoral head. Biomed Res Int. 2015;2015:196495.
Kim TH, Hong JM, Kim HJ, Park EK, Kim SY. Lack of association of MTHFR gene polymorphisms with the risk of osteonecrosis of the femoral head in a Korean population. Mol Cells. 2010;29:343–8.
Liew SC, Gupta ED. Methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism: epidemiology, metabolism, and the associated diseases. Eur J Med Genet. 2015;58:1–10.
Vacek TP, Kalani A, Voor MJ, Tyagi SC, Tyagi N. The role of homocysteine in bone remodeling. Clin Chem Lab Med. 2013;51:579–90.
Milovanovic P, Hrncic D, Radotic K, Stankovic M, Mutavdzic D, Djonic D, Rasic-Markovic A, Djuric D, Stanojlovic O, Djuric M. Moderate hyperhomocysteinemia induced by short-term dietary methionine overload alters bone microarchitecture and collagen features during growth. Life Sci. 2017;191:9–16.
Thanks to proofreader Alejandra Tapia Alcazar for her assistance in the preparation of this manuscript. In loving memory of M.D. Antonio Redón Tavera, head of the Pediatric Hip Service, INR-LGII.
This study had no specific funding.
Ethics approval and consent to participate
All participants received oral and written information about the study and signed a letter of consent. The INR-LGII Research and Ethics Committees reviewed and approved the study protocol.
Consent for publication
All authors consent to publish.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Buendía-Pazarán, J.G., Hernández-Zamora, E., Rodríguez-Olivas, A.O. et al. Association of MTHFR rs1801133 and homocysteine with Legg–Calvé–Perthes disease in Mexican patients. Orphanet J Rare Dis 17, 123 (2022). https://doi.org/10.1186/s13023-022-02264-2