Mutational Study in the Exon2 of Osteoprotegrin Gene in Osteoporotic Women in Chahar Mahal va Bakhtiari Province of Iran

authors:

avatar Morteza Dehghan 1 , avatar Razieh Pourahmad Jaktaji ORCID 2 , * , avatar Seyedeh Zahra Mousavi ORCID 2

Department of Orthopedics, Shahrekord University of Medical Sciences, Shahrekord, Iran
Department of Genetics, Faculty of Science, Shahrekord University, Shahrekord, Iran

how to cite: Dehghan M, Pourahmad Jaktaji R, Mousavi S Z. Mutational Study in the Exon2 of Osteoprotegrin Gene in Osteoporotic Women in Chahar Mahal va Bakhtiari Province of Iran. Zahedan J Res Med Sci. 2023;25(1):e116972. https://doi.org/10.5812/zjrms-116972.

Abstract

Background:

Osteoporosis is a multifactorial disease. Mutation in cystein-rich domain2 (CRD2) and cystein-rich domain3 (CRD3) of osteoprotegrin protein can prevent its interaction with receptor activator of nuclear factor kappa beta (RANKL) ligand and lead to osteoporosis.

Objectives:

This study aimed to investigate the possible alteration in exon2 of the osteoprotegrin gene encoding CRD2 and CRD3 domains in osteoporotic women in Chahar Mahal va Bakhtiari Province of Iran.

Methods:

The N-terminal region of human osteoprotegerin (OPG) protein was aligned with mice osteoprotegrin proteins using Clustal Omega multiple alignment tool. The genomic DNA of 72 osteoporotic women was extracted by commercial kit. The region of exon2 of the OPG gene encoding CRD2 and CRD3 was amplified by PCR and sequenced by DNA sequencing.

Results:

The result of multiple protein alignment showed dissimilarities among three species in terms of CRD2 and CRD3 domains. The results of polymerase chain reaction (PCR) amplification and DNA sequencing indicated that CRD2 and CRD3 sequence were intact in osteoporotic women.

Conclusion:

Due to the binding possibility of OPG protein with RANKL ligand, it was concluded that the expression of OPG gene may have been different in osteoporotic women. However, it was recommended that further studies should be conducted in order to confirm this finding.

1. Background

Osteoporosis is a multi-factorial disease associated with the combination of reduced bone mass density and deterioration of the bone structure. Today, osteoporosis is considered as an important health problem for societies and is called the silent disease of the century. This disease is asymptomatic and its complications (fracture) can impose great and non-compensable financial and physical detriments on society and patients (1). Epidemiological studies have shown lower bone density in Asian women than in European and American women. Since women have an 8-fold increased rate of osteoporosis than men, more than half of women aged over 50 years are afflicted with this disease (2).

There are many candidate genes that have been examined for their involvement in the pathogenesis of osteoporosis and bone fractures, including collagen 1a1 (COL1A1), vitamin D receptor (VDR), calcitonin receptor, and osteoprotegerin (OPG) (3, 4). Osteoprotegerin is a soluble member of the tumor necrosis factor receptor superfamily (TNFRSF), which is synthesized by osteoblasts; however, it is also expressed in various organs and tissues, including heart, vessel walls, lung and kidney (5), and is secreted in human milk (6). Osteoprotegerin regulates a variety of cellular processes in immune system (maturation of B cells and development of antibody responses) and bone metabolism (inhibition of osteoclastogenesis) (6, 7). Osteoclastogenesis leads to bone resorption. The balance between bone resorption by osteoclast and bone formation by osteoblast is necessary for bone remodeling process (8). Imbalance of bone remodeling results in different bone diseases, including osteoporosis, osteopenia, osteoarthritis, intervertebral disk degeneration, and idiopathic hyperphosphatasia (Juvenile Paget’s disease) (5, 9-12).

The human OPG is a single copy gene located on chromosome 8, consisting of 5 exons and 4 introns. Expression of OPG gene is upregulated by different regulators, including IL-1β, tumor necrosis factors, TGF-β, vitamin D3, and bone morphogenetic protein 2 (BMP2), and is downregulated by various regulators, such as glucocorticoids and prostaglandin E2 (PG-E2) and parathyroid hormone (PTH) (9, 12, 13). Alterations in promoter and OPG gene sequence could change OPG expression and function (11, 14), and may be associated with osteoporosis and bone mineral density (4, 15). Identifying the molecular mechanisms that regulate OPG and RANKL expression may help to find new targets for bone therapy purposes.

Osteoprotegerin protects bone from excessive resorption through the interaction and inhibition of its ligand, RANKL. Receptor activator of nuclear factor kappa beta is mainly expressed by osteoblastic lineage cells and activated T cells (16, 17). Parathyroid hormone induces RANKL expression (13). Receptor activator of nuclear factor kappa beta is a critical cytokine for the differentiation, activation, and survival of osteoclasts, and a member of the TNF superfamily. To exert its effect on osteoclast, RANKL must interact with its receptor, RANK, located on the membrane of osteoclast. Similar to OPG, RANK is a member of TNFRSF superfamily. The relative concentrations of RANKL and OPG in bone are important determinants of bone mass and strength, and the blockage of RANK/RANKL signaling has become a therapeutic strategy to cure osteoporosis and other disorders related to increased bone resorption (17).

Osteoprotegerin exists either as a 60 kDa monomer or a 120 kDa disulfide-linked homodimer. Each domain of OPG is composed of 7 structural domains. Domains 1 to 4 of the N-terminal site (cystein-rich domain1 (CRD1), cystein-rich domain2 (CRD2), cystein-rich domain3 (CRD3) and cystein-rich domain4 (CRD4)) are cystein-rich domains and have structural similarity to the extracellular CRD domains of TNFRSP proteins, such as RANK. Osteoprotegerin domains CRD2 and CRD3 participate in the interaction with RANKL and prevention of RANKL binding to RANK. It has been also reported that Glu95 (in CRD2 domain) has interactions with RANKL residues Arg223, Tyr241, and Lys257 (18). Exon 2 of OPG contains 3 CRD domains (i.e., CRD1, CRD2, and CRD3). The location of CRD2 and CRD3 domains in exon2 of OPG is displayed in Figure 1.

The sequence of exon2 of human osteoprotegerin. The arrowhead shows the sequence of forward primer. The red and black rectangles show amino acid sequence of cystein-rich domain2 and cystein-rich domain3 domains, respectively. Star Cys, Glu and Phe are Cys107, Glu95 and Phe117.
The sequence of exon2 of human osteoprotegerin. The arrowhead shows the sequence of forward primer. The red and black rectangles show amino acid sequence of cystein-rich domain2 and cystein-rich domain3 domains, respectively. Star Cys, Glu and Phe are Cys107, Glu95 and Phe117.

2. Objective

Due to the importance of CRD2 and CRD3 domains of OPG in the interaction with RANKL, this study aimed to examine exon2 containing CRD2 and part of CRD3 domains in osteoporotic women aged 45 years and older who lived in Chaharmahal va Bakhtiari Province of Iran and referred to Bone Mineral Density Assay Centers in Shahrekord.

3. Methods

3.1. Amino Acid Alignment

To see whether amino acid sequence of CRD2 and CRD3 domains were conserved among mammalian species, the partial amino acid sequence of OPG protein containing exon2 from 3 mammalian species was aligned via Clustal Omega multiple alignment tool (link: www.ebi.ac.uk/Tools/msa/clustalo/).

3.2. Selection of Osteoporotic Women

Out of all women referring to bone densitometry centers of Shahrekord (Chahar Mahal va Bakhtiari Province, Iran), a total of 144 cases were included in this study in years 2014 - 16. First, the participants were asked to complete informed consent forms before entering the study. Patients with a history of corticosteroid use, ovariectomy or early ovarian failure, thyroid disease, impaired calcium absorption, as well as digestive and kidney diseases were excluded from the study. Our study was prospectively reviewed and approved by the ethics committee of Shahrekord Medical Science University (IR.SKUMS.REC.1394.180), and was conducted in accordance with the guidelines and principles of the Helsinki Declaration (19).

Bone density of lumbar spine and femoral neck were determined by dual energy X-ray absorptionmtery (DEXA) using Hologic QDR device (Germany) and were described based on T score method. Then 72 osteoporotic women were selected for further study based on T score.

3.3. DNA Extraction

Genomic DNA was extracted from blood samples using DNA extraction kit (Cinaclone, Iran) according to supplier manual and was stored at -20°C. The quality and quantity of extracted DNA were evaluated using the spectrophotometer (Ultrospec 1100, Amersham Pharmacia Biothech, Australia) and agarose gel electrophoresis (horizontal gel electrophoresis set, paya pajohesh, Iran).

3.4. Polymerase Chain Reaction Amplification

The sequence of exon 2 OPG gene (Gene ID 4982) containing CRD2 and CRD3 domains was amplified using primers designed through Gene Runner v4.0.9.67 Beta and PCR method. The primers used were: Forward, 5’-GTACAGCAAAGTGGAAGA-3’ (as shown by arrowhead in Figure 1) and reverse, 5’-CTCCTAAACTGTCACAACTA-3’ (in intron). As for PCR reaction, 2.5 μL buffer 10X, 0.5 mL MgCl2 50 mM, 0.5 μL Mix 40 mM dNTP, 0.5 μL forward primer 10 pM, 0.5 μL reverse primer 10 pM, 0.2 μL Taq polymerase (0.5 U/μL), and 2.3 μL DNA template (50 ng) were used. The final reaction volume was reached to 25 μL using ddH2O. Thermal condition for PCR reaction was: 95°C for 5 min followed by 35 cycles of 94°C for 30 s, 45°C for 30 s, 72°C for 30 s with final extension at 72°C for 10 min.

3.5. DNA Sequencing

Amplified fragments were further analyzed by adopting DNA sequencing method to show the possible change in exon2 of OPG.

4. Results

The result of amino acid sequence alignment of OPG protein from 3 mammal species using Clustal Omega server showed that human OPG has approximately 81 percent similarity with mouse OPGs, and most of dissimilarities were outside of CRD2 and CRD3 in exon2 (Figure 2). The cystein residues and Glu95 were in the same locations in 3 species.

Partial amino acid sequence alignment of osteoprotegerin from 3 mammalian species. NP_002537.3, NP_032790.3 and NP037002.1 represent amino acid sequence of Homo sapiens, Mus muscullus and Rattus norvegicus, respectively. Black and grey arrowheads show start and end of exon2 and cystein-rich domain2, respectively. White arrowhead shows the start of cystein-rich domain3. Vertical rectangles show Cys107, Glu95, Phe117 Cys118 residues.
Partial amino acid sequence alignment of osteoprotegerin from 3 mammalian species. NP_002537.3, NP_032790.3 and NP037002.1 represent amino acid sequence of Homo sapiens, Mus muscullus and Rattus norvegicus, respectively. Black and grey arrowheads show start and end of exon2 and cystein-rich domain2, respectively. White arrowhead shows the start of cystein-rich domain3. Vertical rectangles show Cys107, Glu95, Phe117 Cys118 residues.

After conducting PCR amplification and gel electrophoresis, the PCR product was observed on agarose gel (Figure 3). The size of product was 303 bp. Polymerase chain reaction products were sequenced and compared with human OPG gene sequence retrieved from NCBI. The results of DNA sequencing of CRD2 and part of CRD3 in exon2 in 72 osteoporotic women did not show any changes in these domains (Figures 4 and 5). As displayed in figures, all cystein residues and Glu95 were intact.

Agarose gel electerophoresis of polymerase chain reaction (PCR) products. The last lane from right shows size marker (1 kb). Rest of lanes present PCR products with size of 303 bp.
Agarose gel electerophoresis of polymerase chain reaction (PCR) products. The last lane from right shows size marker (1 kb). Rest of lanes present PCR products with size of 303 bp.
DNA Sequencing of cystein-rich domain2 (rectangle). Star Glu is Glu95.
DNA Sequencing of cystein-rich domain2 (rectangle). Star Glu is Glu95.
DNA sequencing of part of cystein-rich domain3. Star Cyss and Phe are Cys107, Phe117 and Cys118.
DNA sequencing of part of cystein-rich domain3. Star Cyss and Phe are Cys107, Phe117 and Cys118.

5. Discussion

Osteoporosis is a common and serious disease worldwide. Environmental and genetic factors control the BMD value and, thereby, the osteoporotic fractures. Previous studies have shown that alterations in human OPG gene could be associated with osteoporosis or its phenotypes such as osteoporotic fractures and decrease in BMD (4, 15). Chong et al. suggested the involvement of OPG mutations with the idiopathic hyperphosphatasia (IH) juvenile Paget’s disease phenotype (11). This congenital bone disease is a rare autosomal recessive disorder that causes deformity and physical disability. They identified 6 recessive IH juvenile-associated mutations in OPG gene, three of which were located in exon 2, including Cys65Arg, Cys87Tyr, and Phe117Leu. Alteration of these cysteins to other amino acids destroys disulfide bonds in the second CRD2 of OPG which is the second most effective domain in OPG interaction with RANKL. Furthermore, alteration of Phe117 in CRD3 to other amino acid may affect the formation of disulfide bond between Cys118 and Cys107 (11) and decreases OPG binding to RANKL (20). In this study, therefore, it was decided to examine part of exon2 of OPG containing CRD2 and part of CRD3 in order for identifying possible mutation in 72 osteoporotic women in Chahar Mahal va Bakhtiari Province of Iran. It was found that majority of amino acids -cysteins, in particular- were conserved in CRD2 and CRD3 in comparison with mice OPG, and that these regions were intact in osteoporotic women. This may have implied that there was no defect in the interaction between OPG and RANKL in these osteoporotic women. Moreover, it was indicated that women who were homozygous for Phe117 had higher spine and neck hip BMD in comparison with those who were heterozygous (Phe/Leu) or homozygous (Leu/Leu) for position 117 (21). The same situation was found for women who were homozygous for Cys87 in comparison with those who were heterozygous (Cys/Tyr) or homozygous (Tyr/Tyr) for position 87 (22).

Since osteoporosis is a multifactorial disease, there might be another reason for having low BMD and frequent bone fracture in present cases. It has been argued that the ratio of RANKL to OPG is important for the regulation of osteoclast differentiation and function (8). It has been also indicated that the level of OPG in patients with a long bone fracture is higher than that in healthy controls. In contrast, the level of RANKL decreases during regular fracture repair (23). Alterations in promoter of OPG gene could change OPG expression (15) and may be associated with osteoporosis. Alteration in four locations in OPG promoter, -960, -946, -900, and -864, either alone or in combination, produces different haplotypes and affects the ability of transcription factors (e.g., specificity protein1 (Sp1) and NF-1) to attach to the region containing these sites, leading to various levels of transcription (14).

Moreover, the degree of methylation of CpG sites in the CpG island of OPG/RANKL promoters affects these genes expression, so that higher methylation degree causes lower gene expression (9).

In addition to OPG, an association was detected between some single nucleotide polymorphisms (SNPs) in COL1A1, VDR, and calcitonin receptor genes with bone density in Iranian osteoporotic women living in Chahar Mahal va Bakhtiari Province, Iran (24-26). Taq1 and Bsm1 polymorphisms of VDR as well as Sp1 polymorphism of COL1A1 were also found to have a significant association with bone density (P < 0.05) (18, 19). However, it was recommended that further studies should be carried out to examine OPG expression in these women and extra cases and show the role of OPG in order for finding more information about the pathogenesis of osteoporosis in osteoporotic women.

5.1. Conclusions

Since the majority of amino acids and specially cysteins were intact in the studied osteoporotic women, it was concluded that there was no defect in interaction between OPG and RANKL. Due to the importance of OPG protein in bone remodeling, however, it was recommended OPG gene expression in osteoporotic women should be further investigated.

5.2. Limitations and Suggestions

Due to time and financial limitations, performing extra PCR amplification and DNA sequencing was not possible. Therefore, it was suggested that a larger sample of osteoporotic women should be selected and investigated in order for identifying possible mutation in CRD2 and CRD3 domains.

References

  • 1.

    Munch S, Shapiro S. The silent thief: osteoporosis and women's health care across the life span. Health Soc Work. 2006;31(1):44-53. [PubMed ID: 16550847]. https://doi.org/10.1093/hsw/31.1.44.

  • 2.

    Khoury MJ. Genetic and epidemiologic approaches to the search for gene-environment interaction: the case of osteoporosis. Am J Epidemiol. 1998;147(1):1-2. [PubMed ID: 9440392]. https://doi.org/10.1093/oxfordjournals.aje.a009360.

  • 3.

    Marozik P, Mosse I, Ameliyanovich M, Rudenka E, Alekna V, Tamulaitienė M. Molecular and genetic mechanisms of predisposition to osteoporosis. Gerontologija. 2011;12(4):250-8.

  • 4.

    Langdahl BL, Carstens M, Stenkjaer L, Eriksen EF. Polymorphisms in the osteoprotegerin gene are associated with osteoporotic fractures. J Bone Miner Res. 2002;17(7):1245-55. [PubMed ID: 12096838]. https://doi.org/10.1359/jbmr.2002.17.7.1245.

  • 5.

    Xue JB, Zhan XL, Wang WJ, Yan YG, Liu C. OPG rs2073617 polymorphism is associated with upregulated OPG protein expression and an increased risk of intervertebral disc degeneration. Exp Ther Med. 2016;12(2):702-10. [PubMed ID: 27446264]. [PubMed Central ID: PMC4950610]. https://doi.org/10.3892/etm.2016.3342.

  • 6.

    Vidal K, van den Broek P, Lorget F, Donnet-Hughes A. Osteoprotegerin in human milk: a potential role in the regulation of bone metabolism and immune development. Pediatr Res. 2004;55(6):1001-8. [PubMed ID: 15155868]. https://doi.org/10.1203/01.pdr.0000127014.22068.15.

  • 7.

    Sano T, Akeda K, Yamada J, Takegami N, Sudo T, Sudo A. Expression of the RANK/RANKL/OPG system in the human intervertebral disc: implication for the pathogenesis of intervertebral disc degeneration. BMC Musculoskelet Disord. 2019;20(1):225. [PubMed ID: 31101043]. [PubMed Central ID: PMC6525349]. https://doi.org/10.1186/s12891-019-2609-x.

  • 8.

    Lv WT, Du DH, Gao RJ, Yu CW, Jia Y, Jia ZF, et al. Regulation of Hedgehog signaling Offers A Novel Perspective for Bone Homeostasis Disorder Treatment. Int J Mol Sci. 2019;20(16). [PubMed ID: 31426273]. [PubMed Central ID: PMC6719140]. https://doi.org/10.3390/ijms20163981.

  • 9.

    Wang P, Cao Y, Zhan D, Wang D, Wang B, Liu Y, et al. Influence of DNA methylation on the expression of OPG/RANKL in primary osteoporosis. Int J Med Sci. 2018;15(13):1480-5. [PubMed ID: 30443169]. [PubMed Central ID: PMC6216050]. https://doi.org/10.7150/ijms.27333.

  • 10.

    Takegami N, Akeda K, Yamada J, Sano T, Murata K, Huang J, et al. RANK/RANKL/OPG system in the intervertebral disc. Arthritis Res Ther. 2017;19(1):121. [PubMed ID: 28576140]. [PubMed Central ID: PMC5457592]. https://doi.org/10.1186/s13075-017-1332-y.

  • 11.

    Chong B, Hegde M, Fawkner M, Simonet S, Cassinelli H, Coker M, et al. Idiopathic hyperphosphatasia and TNFRSF11B mutations: relationships between phenotype and genotype. J Bone Miner Res. 2003;18(12):2095-104. [PubMed ID: 14672344]. https://doi.org/10.1359/jbmr.2003.18.12.2095.

  • 12.

    Wan M, Shi X, Feng X, Cao X. Transcriptional mechanisms of bone morphogenetic protein-induced osteoprotegrin gene expression. J Biol Chem. 2001;276(13):10119-25. [PubMed ID: 11139569]. https://doi.org/10.1074/jbc.M006918200.

  • 13.

    Fu Q, Jilka RL, Manolagas SC, O'Brien CA. Parathyroid hormone stimulates receptor activator of NFkappa B ligand and inhibits osteoprotegerin expression via protein kinase A activation of cAMP-response element-binding protein. J Biol Chem. 2002;277(50):48868-75. [PubMed ID: 12364326]. https://doi.org/10.1074/jbc.M208494200.

  • 14.

    Hagedorn C, Telgmann R, Dordelmann C, Schmitz B, Hasenkamp S, Cambien F, et al. Identification and functional analyses of molecular haplotypes of the human osteoprotegerin gene promoter. Arterioscler Thromb Vasc Biol. 2009;29(10):1638-43. [PubMed ID: 19661483]. https://doi.org/10.1161/ATVBAHA.109.193532.

  • 15.

    Arko B, Prezelj J, Komel R, Kocijancic A, Hudler P, Marc J. Sequence variations in the osteoprotegerin gene promoter in patients with postmenopausal osteoporosis. J Clin Endocrinol Metab. 2002;87(9):4080-4. [PubMed ID: 12213850]. https://doi.org/10.1210/jc.2002-020124.

  • 16.

    Soufi M, Schoppet M, Sattler AM, Herzum M, Maisch B, Hofbauer LC, et al. Osteoprotegerin gene polymorphisms in men with coronary artery disease. J Clin Endocrinol Metab. 2004;89(8):3764-8. [PubMed ID: 15292302]. https://doi.org/10.1210/jc.2003-032054.

  • 17.

    Delgado-Calle J, Sanudo C, Fernandez AF, Garcia-Renedo R, Fraga MF, Riancho JA. Role of DNA methylation in the regulation of the RANKL-OPG system in human bone. Epigenetics. 2012;7(1):83-91. [PubMed ID: 22207352]. [PubMed Central ID: PMC3337833]. https://doi.org/10.4161/epi.7.1.18753.

  • 18.

    Luan X, Lu Q, Jiang Y, Zhang S, Wang Q, Yuan H, et al. Crystal structure of human RANKL complexed with its decoy receptor osteoprotegerin. J Immunol. 2012;189(1):245-52. [PubMed ID: 22664871]. https://doi.org/10.4049/jimmunol.1103387.

  • 19.

    World Medical Association. World Medical Association Declaration of Helsinki: ethical principles for medical research involving human subjects. JAMA. 2013;310(20):2191-4. [PubMed ID: 24141714]. https://doi.org/10.1001/jama.2013.281053.

  • 20.

    Nelson CA, Warren JT, Wang MW, Teitelbaum SL, Fremont DH. RANKL employs distinct binding modes to engage RANK and the osteoprotegerin decoy receptor. Structure. 2012;20(11):1971-82. [PubMed ID: 23039992]. [PubMed Central ID: PMC3607351]. https://doi.org/10.1016/j.str.2012.08.030.

  • 21.

    Liu S, Yi Z, Ling M, Shi J. Association between g.19163A>G and g.23298T>C genetic variants of the osteoprotegerin gene and bone mineral density in Chinese women. Hormones (Athens). 2013;12(4):578-83. [PubMed ID: 24457406]. https://doi.org/10.14310/horm.2002.1446.

  • 22.

    Zhang YD, Zhang Z, Zhou NF, Jia WT, Cheng XG, Wei XJ. Association of the g.19074G>A genetic variant in the osteoprotegerin gene with bone mineral density in Chinese postmenopausal women. Genet Mol Res. 2014;13(3):6646-52. [PubMed ID: 25177945]. https://doi.org/10.4238/2014.August.28.9.

  • 23.

    Kottstorfer J, Thomas A, Gregori M, Kecht M, Kaiser G, Eipeldauer S, et al. Are OPG and RANKL involved in human fracture healing? J Orthop Res. 2014;32(12):1557-61. [PubMed ID: 25212894]. https://doi.org/10.1002/jor.22723.

  • 24.

    Dehghan M, Pourahmad-Jaktaji R. Sp1 binding site polymorphism of a collagen gene (rs 1800012) in women aged 45 and over and its association with bone density. Turk J Med Sci. 2015;45(3):644-50. [PubMed ID: 26281333]. https://doi.org/10.3906/sag-1405-80.

  • 25.

    Dehghan M, Pourahmad-Jaktaji R. The Effect of Some Polymorphisms in Vitamin D Receptor Gene in Menopausal Women with Osteoporosis. J Clin Diagn Res. 2016;10(6):RC06-10. [PubMed ID: 27504361]. [PubMed Central ID: PMC4963721]. https://doi.org/10.7860/JCDR/2016/17147.8006.

  • 26.

    Dehghan M, Pourahmad-Jaktaji R, Farzaneh Z. Calcitonin Receptor AluI (rs1801197) and Taq1 Calcitonin Genes Polymorphism in 45-and Over 45-year-old Women and their Association with Bone Density. Acta Inform Med. 2016;24(4):239-43. [PubMed ID: 27708484]. [PubMed Central ID: PMC5037997]. https://doi.org/10.5455/aim.2016.24.239-243.