International Journal of
Reproductive Biomedicine
Sun, May 19, 2024
[Archive]
Volume 20, Issue 2 (February 2022)
IJRM 2022, 20(2): 111-122 |
Back to browse issues page
Download citation:
BibTeX | RIS | EndNote | Medlars | ProCite | Reference Manager | RefWorks
Send citation to:
BibTeX | RIS | EndNote | Medlars | ProCite | Reference Manager | RefWorks
Send citation to:
Nemati Vahedi S, Kheirkhah B, Malekirad A A, Hosseini S M. Association of selected polymorphisms in GPX4, COMT, pre-miR-125a, pre-miR-10a, and pre-miR-323b genes in Iranian women with idiopathic recurrent pregnancy loss: A case-control study. IJRM 2022; 20 (2) :111-122
URL: http://ijrm.ir/article-1-1968-en.html
URL: http://ijrm.ir/article-1-1968-en.html
1- Department of Cellular and Molecular Biology, Faculty of Basic Sciences, Islamic Azad University Sirjan Branch, Sirjan, Iran.
2- Department of Microbiology, Faculty of Veterinary Medicine, Islamic Azad University Baft Branch, Baft, Iran. , BabakKheirkhah@yahoo.com
3- Department of Biology, Faculty of Basic Sciences, Payame Noor University, Tehran, Iran.
4- Human Genetic Research Center, Baqiyatallah University of Medical Sciences, Tehran, Iran. Department of Molecular Genetics, Baqiyatallah Hospital, Iran University of Medical Sciences, Tehran, Iran.
2- Department of Microbiology, Faculty of Veterinary Medicine, Islamic Azad University Baft Branch, Baft, Iran. , BabakKheirkhah@yahoo.com
3- Department of Biology, Faculty of Basic Sciences, Payame Noor University, Tehran, Iran.
4- Human Genetic Research Center, Baqiyatallah University of Medical Sciences, Tehran, Iran. Department of Molecular Genetics, Baqiyatallah Hospital, Iran University of Medical Sciences, Tehran, Iran.
Full-Text [PDF 1935 kb]
(678 Downloads)
| Abstract (HTML) (1094 Views)
1. Introduction
Recurrent pregnancy loss (RPL) is defined as “two or more abortions before 24 weeks of gestation” (1, 2). According to previous studies (3), multiple factors are involved in this condition, i.e. chromosomal abnormalities, genetic disorders, uterine structural abnormalities, immune diseases, thrombophilia, thyroid abnormalities, infectious diseases, high parental age, and lifestyle. However, the reason for more than 50% of miscarriages is still unknown (4). Among the factors affecting RPL, genes responsible for fetal growth and development have considerable roles.
Single nucleotide polymorphisms (SNPs) in several genes for proteins or fetal development are associated with RPL. Phospholipid hydroperoxide glutathione peroxidase (GPX4) is involved in glutathione synthesis (5) and may be associated with RPL in women; however, the exact mechanism is not clear. Catechol-O-methyl transferase protein (COMT) is now increasingly reported to have a role in pregnancy. In particular, recent reports have noted its expression in fetal membranes. This enzyme is also active both in the placenta and in the decidua (6). Reports have demonstrated that hypertension in pregnancy can reduce the placental COMT activity (7). Recent studies have shown that women with severe preeclampsia have reduced placental COMT protein expression (8). Besides, micro RNAs (miRs) are among the newly discovered genetic factors associated with RPL. miRs are small non-protein coding RNAs causing gene silencing at the post-transcription level. They target about one-third of human gene expressions and also play a critical role in various gestational events such as implantation, ectopic pregnancy, pre-term labor, pre-eclampsia, low birth weight, RPL, and anti-inflammatory responses. Recent investigations on women with RPL have shown that the overexpression of miR-133a can lead to a decline in human leukocyte antigen G (HLA-G) translation by joining to unranslated regions (UTR-3'). This protein has an essential role in the mother's immune system tolerating the fetus, and the down-regulation of this protein leads to RPL (9). miR-125 is a group of miRs which consists of two subgroups: miR-125a and miR-125b. The expression of these miRs occurs in different organs like the stomach, liver, lungs, rectum, breast, prostate, and ovaries. Research on these organs suggests that the malignant process may be controlled by preventing gene expression. Embryo development may also be affected by low secretion of the mother’s miR-125a polymorphism. Previous study have also shown that women with RPL have downregulation of 41 miRs and overexpression of four miRs (10). Although the critical role of miR-125a in RPL is unclear, recent evidence has proposed that its down-regulation may induce RPL in multiple ways (11).
By considering the wide range of miR activities in messaging processes, especially in the female genital system, studies on these regulatory RNAs seem to be necessary. The importance of miscarriage necessitates more comprehensive studies on the role of miRs in RPL. Therefore, the current study aimed to use the SNaPshot method to investigate the association between polymorphisms of pre-miR-125a (rs41275794), pre-miR-10a (rs3809783), pre-miR-323b (rs56103835), GPX4 (rs4680) and GPX4 (rs713041), and RPL in Iranian women with idiopathic recurrent miscarriage.
2. Materials and Methods
2.1. Sample collection
In the present case-control study, 116 women with idiopathic RPL and 89 healthy women as controls who had previously had at least two successful pregnancies (a total of 205 individuals) were included and referred to a Baqiyatallah hospital (Tehran, Iran) for sampling during 2017-2019.
A diagnosis of idiopathic RPL was confirmed by the expert gynecologist after evaluation of hormonal levels, immunological and hematological deficiencies, karyotyping, and uterine anomalies. Uterine anomalies were evaluated through hysterosalpingography, sonhysterography and hysteroscopy techniques. Healthy controls visited the hospital for checkups and no abnormalities were found in their physical or laboratory examinations. They did not have any history of RPL or previous related medical diseases.
Inclusion criteria were: (i) women with idiopathic RPL; and (ii) who had undergone complete clinical and pathologic examinations. Cases who met the following criteria were excluded from the study: (i) history of autoimmune or hormonal diseases; (ii) chromosomal aneuploidies or genetic abnormalities; and/or (iii) immunological or hematological disorders.
2.2. DNA extraction
In the first examination, after an overnight fast, 5 ml peripheral blood samples were collected and placed in the EDTA-containing tubes. Then, the rapid genomic DNA extraction (RGDE) method was used to extract genomic DNA (12).
Subsequently, a Nanodrop ND-1000 spectrophotometer (Thermo Sci., Newington, NH, USA) was used to determine the quantity and quality of the extracted DNA. A fraction of each DNA sample was placed on 1.5% denaturing agarose gel stained with Red Safe and was then subjected to electrophoresis to assess the DNA quality. Following that, the isolated DNA was stored at -20°C until further analysis.
2.3. Uniplex polymerase chain reaction (PCR) amplification
After DNA extraction, the functional SNPs pre-miR-125a (rs41275794), pre-miR-10a (rs3809783), pre-miR-323b (rs56103835), GPX4 (rs4680) and GPX4 (rs713041) were selected in each gene for genotyping based on published data and Primer-BLAST in the National Center for Biotechnology Information dbSNP database.
Oligo analyzer software was used to consider or predict the formation of dimer or secondary structures of primer. Primer sequences and their amplicon size are shown in table I. Gene runner software was used to design appropriate single-base extension (SBE) primers, as mentioned in table II. All individuals were analyzed for the presence of these polymorphisms (13).
2.4. Multiplex PCR amplification
To conduct the multiplex PCR, the following chemicals were prepared: 1 μM of each primer, 0.17 mM of dNTPs, 1.5 mM of MgCl2, one unit of Taq DNA polymerase, and 20 ng of template DNA in a 20 μl volume flask containing 10 X PCR buffer. The three conditions applied to the thermal cycle were: 1) 10 min at 95°C; 2) 30 sec at 95°C for 35 cycles; 3) 30 sec at 60°C, one min at 72°C and a final extension of 10 min at 72°C. Then, 3 μl of PCR product was run in 12% polyacrylamide gel to check for the quality and yield of the multiplex PCR. Next, the excess dNTPs and primers were removed by treating the 3 μl of remaining PCR products with one unit (1 μl) of exonuclease I. Then, three units of shrimp alkaline phosphatase and two units of exonuclease I (ExoSAP/USB, Cleveland, OH, USA) were used to purify the remaining PCR products, which was followed by incubation for 15 min at 37oC and 15 min at 80oC to remove excess dNTPs and primers, and to deactivate the Taq enzyme (14).
2.5. SNaPshot method
The SNaPshot assay applied five SBE primers which simultaneously annealed adjacent to SNP variant in our reaction. As shown in table II, poly (dC) tails of different lengths were attached to the SBE primers.
The SNaPshot solutions for multiplex PCR contained PCR product (1 μl), master mix SNaPshot (1 μl), each SBE primer (1 μl) and DDW (2 μl) in a final volume of 5 μl. The three thermal protocols applied to the SNaPshot were: two min at 96°C; 10 sec at 96°C; five sec at 50°C and 30 sec at 60°C. dNTPs were removed from the second PCR product by treating the remaining PCR products with three units (1 μl) of shrimp alkaline phosphatase and then the following two thermal protocols: 45 min at 37°C and 15 min at 75°C. An ABI PRISM 3130 Genetic analyzer was used to conduct capillary electrophoresis after which Gene Mapper idx, Ver 4.0 (Applied Biosystems, Foster City, CA, USA) was used to analyze the results (14).
2.6. Ethical considerations
This case-control study was approved by the institutional review board and ethical consideration committee of Baqiyatallah University of Medical Sciences (Code: IR.BMSU.REC.1398.054). Written informed consent was obtained from all individuals.
2.7. Statistical analysis
Statistical analysis was applied on SPSS v. 21 statistics software. Data were showed as percentages of the mean or allele frequency. Pearson’s chi-square test was used to calculate inter-group significance. As well as, the homozygous and heterozygous genotypes of each group were unified as carriers and the odds ratios (OR) and 95% confidence intervals (CI) were obtained. Statistical significance was set at p < 0.05.
3. Results
The present study was performed on 116 cases with RPL and 89 healthy women without any history of RPL. There was no significant difference between their ages (38.93 ± 10.45 years. vs. 41.31 ± 10.74 years, respectively; p = 0.40). After collection of blood samples, the DNA was extracted with the rapid genomic DNA extraction method, as shown in figure 1.
The Uniplex and Multiplex PCR techniques were applied to detect G>A polymorphism in pre-miR-125a (rs41275794), A>T polymorphism in pre-miR-10a (rs3809783), T>A polymorphism in pre-miR-323b (rs56103835), A>G polymorphism in GPX4 (rs4680) and C>T polymorphism in GPX4 (rs713041) genes, as shown in figure 2 and figure 3, respectively.
After multiplex PCR, the SNaPshot kit and specific SBE primers were used to detect rs713041, rs4680, rs41275794, rs3809783, and rs56103835 SNPs with the genetic analyzer machine, as shown in figure 4.
The comparison of SNP genotype distribution for all of the studied genes in the control and case groups is summarized in table III. The frequencies of AG, GG and AA genotypes in the polymorphism rs41275794 of pre-miR-125a (p < 0.001) and of TT, AT and AA genotypes in the polymorphism rs3809783 of pre-miR-10a (p = 0.04) were calculated for the women in the case and control groups. The findings showed significant differences and an association of these polymorphisms with the symptoms of recurrent miscarriage (p < 0.05). There was no significant difference (p = 0.09) in the frequency of TT, AT and AA genotypes in the GPX4 gene (rs713041 polymorphism) between the case and control groups. The frequency of TT, AT and AA genotypes in the miR-323b gene (rs56103835 polymorphism) showed no significant difference (p = 0.14) between the case and control groups. Also, no significant difference (p = 0.21) was observed in the frequency of GG, AG or AA genotypes in the COMT gene (rs4680 polymorphism) between the two groups (Table III). In the allele frequency study for alleles A and G of polymorphism rs41275794, a significant difference was found between the case and control groups (p < 0.001). Similarly, there was a significant difference (p = 0.02) in the frequency of alleles A and T of polymorphism rs3809783 between the two groups of cases and controls (Table IV). The assessment of risk at a 95% confidence level for allele A relative to G in the rs41275794 polymorphism showed that by changing the G allele to A, the risk of the RPL phenotype was increased by 1.912 times. Similarly, for the rs3809783 polymorphism, the changing of allele A to T increased the risk of the RPL phenotype by 1.704 times, showing a greater association of risk. To evaluate the genotypic inheritance of significant polymorphisms, we studied and calculated the risk of RPL in three models: dominant, recessive, and obvious models. The results showed that the dominant hereditary model had a significant correlation based on the genotypic frequency of the variant (Table V).
4. Discussion
RPL is a clinical challenge that affects approximately 1-5% of women of childbearing age (11). The underlying mechanism of RPL in 37-79% of cases is unclear. Recent evidence has revealed that molecular and genetic disorders, including various polymorphisms, may be the main cause of RPL in these cases (15).
In this research, we compared for the first time five polymorphisms of pre-miR-125a (rs41275794), pre-miR-10a (rs3809783), pre-miR-323b (rs56103835), GPX4 (rs4680) and GPX4 (rs713041) between healthy women and women with RPL using the SNaPshot genotyping method. Our findings showed a significant difference in the miR-125a gene polymorphism rs41275794 between the control and case groups. GG was found to be the most common genotype in the control subjects (60.7%), while GA (50.9%) was the most common genotype in the case group. These data indicated that the polymorphism in miR-125a may be a risk factor for RPL among the Iranian female population. In our study, allele A was significantly associated with an increased risk of RPL indicating that this variant allele in pre-miR125a can change precursor and mature miR-125a expression.
Our data revealed that the dominant genotype was significantly associated with an increased risk of RPL. In other study, expression of mature miR-125a in the CC/AA genotype was ~2.1-fold higher than that in the TT/AA genotype (16). They found that the A>T haplotype in rs41275794 affected miR-125a expression and caused RPL in women. In another study, Li and colleagues compared the expression pattern of six miRs, including miR-125a, miR-125b, miR-34a, miR-155, miR-24 and miR-141 between healthy women and those with RPL. They found that miR-125a, miR-125b, miR-34a, miR-155 and miR-141 were significantly overexpressed in women with RPL while miR-24 was significantly down-regulated in these cases, compared to the healthy group. Additionally, they showed that the A>T polymorphism in pri-miR-10a (rs3809783) was significantly associated with an increased risk of RPL in women.
More studies are being conducted nowadays to find relationships between SNPs and the miRNA system's critical components, biological functions and disease (4, 16, 17). The current study indicated that Iranian women with rs3809783 A>T in pri-miR-10a were significantly associated with an increased risk of RPL. This study’s findings also suggested that the expression level of the miR-10a target gene, Bim, might have been reduced because of miR-10a rs3809783. Additionally, this study showed a significant correlation between the miR-10a rs3809783 polymorphism and the risk of RPL. AT was the most common genotype in the case group while the control group most commonly exhibited the AA genotype. The frequency of the T allele in the case group (30.2%) was significantly greater than that of the control group (20.2%). On the other hand, the control group had a significantly higher frequency (79.8%) of the A allele compared to the case group (69.8%).
The current study also found a significant association between the dominant genotype and an increased risk of RPL. Our data showed that polymorphisms in miR-125a and miR-10a can be associated with the risk of RPL in women; however, further studies are needed to confirm this. Several studies have considered the role of different miRs in RPL (17, 18). In a case-control study among Korean women, RPL cases had significantly different frequencies of the miR-499TC + CC and miR-196a2CC genotypes than the control group. It was shown that these two polymorphisms and their combinations were significantly associated with the incidence of idiopathic RPL in the study population (19). In addition, in another study conducted in north Indian females show the miR-196a2CC, miR-499TC + CC, and miR-196a2CC/miR-499TC + CC were associated with an increased risk of RPL (20).
Furthermore, cumulative gene risk score analysis were calculated mentioning oxidative stress-related genes suggested that, common variations in individual genes like GPX4 and to a greater extent in the multi-locus model (ABCB1, GPX4, GPX4, OGG1 genes) are associated with RPL (21). In the current study, we did not find a significant difference in miR-323b and GPX4 polymorphisms between the healthy and case groups. Previous studies have reported that SNPs in the GPX4 gene significantly affect enzyme activity (6, 22, 23). The evidence shows that pre-eclampsia pathogenesis is associated with SNPs in this gene (6). A recent cohort study on a Korean population of 164 women with preeclampsia and 182 normotensive cases showed that preeclampsia was associated with the Val108/158Met polymorphism (6). Nevertheless, the variation in GPX4 enzyme activity is likely not due to the Val108/158Met polymorphism alone. Additionally, various vascular health impairments can be due to reduced levels of estrogen metabolites, such as 2-methoxyestradiol, because of a decrease in GPX4 activity (6, 24). Hypoxia-inducible factor-1a (HIF-1a) plays an essential role in angiogenesis and studies have found that this factor can be suppressed by anti-angiogenic 2-methoxyestradiol (6). Additionally, this transcription factor can induce genes such as soluble fms-like tyrosine kinase-1 (sFlt-1) which facilitates the adaptation and survival of cells during low-oxygen levels. Studies have also shown the stimulation of trophoblast invasion due to the combination of 2-methoxyestradiol and low levels of oxygen (25-27). Placental pathology can also be affected throughout the different stages of pregnancy if the homeostasis and regulation of GPX4 and 2-methoxyestradiol are disrupted. Also, it has been hypothesized that hypoxia-driven trophoblast invasion and vascular remodeling can be disrupted with a premature increase in 2-methoxyestradiol, thereby contributing to preeclampsia pathogenesis (6, 25). Decreased GPX4 activity in late pregnancy can cause reduced levels of 2-methoxyestradiol and thus reduce HIF-1a inhibition, resulting in possible vascular pathology and inflammation (6). Therefore, future studies on the expression and activity of the GPX4 enzyme throughout pregnancy can pinpoint the effect of GPX4/2-methoxyestradiol at different stages of miscarriage.
Since there is no similar report on the frequency of selected polymorphisms in the Iranian female population, one of the aims of the current study was to investigate the possible role of these polymorphisms in the occurrence of spontaneous miscarriage. A strength of this study was the careful selection of eligible individuals within the case group, i.e. women who had had at least two recurrent miscarriages, were completely healthy in terms of anatomical and cytogenetic studies and did not have any significant clinical problems in their medical records. The findings indicated that genetic polymorphisms can influence RPL, but the lack of significant differences in the GPX4, COMT and miR-323b polymorphisms in the Iranian females also reinforced the hypothesis that a single miR or polymorphism is unlikely to be solely responsible for RPL. One of the limitation of our study was the diversity of ethnicity in the Iranian population of the individuals who participate in investigation and lack of indigenous genetically database which harbor allele frequencies and other genomic data for estimation calculation samples collection and etc.
It is recommended that the expression of these polymorphisms in the blood as well as in the uterine tissue of women with RPL be measured to create a better understanding of the molecular mechanisms leading to miscarriage and to compare different conditions. Moreover, this research can be utilized to provide a clearer picture of the mechanisms involved in causing a miscarriage in women and to analyze a complete profile of polymorphisms involved in the various pathways leading to miscarriage.
Bioinformatics studies on the secondary and tertiary protein structures of selected genes using related databases can provide more comprehensive information on recurrent miscarriages.
5. Conclusion
In conclusion, our data illustrated that a polymorphism in pre-miR-125a (rs41275794) and in pre-miR-10a (rs3809783) can be associated with RPL in women. SNaPshot was shown to be a suitable method for investigation of the role of polymorphisms in the incidence of health conditions.
Acknowledgments
We thank Dr. Sayed Mostafa Hosseini from Baqiyatallah University of Medical Sciences for their technical support and Mahshid Shahrzad Kavkani for their assistance. The author(s) received no financial support for the research, authorship, and/or publication of this article.
Conflict of Interest
The authors declare that there is no conflict of interest.
Full-Text: (301 Views)
1. Introduction
Recurrent pregnancy loss (RPL) is defined as “two or more abortions before 24 weeks of gestation” (1, 2). According to previous studies (3), multiple factors are involved in this condition, i.e. chromosomal abnormalities, genetic disorders, uterine structural abnormalities, immune diseases, thrombophilia, thyroid abnormalities, infectious diseases, high parental age, and lifestyle. However, the reason for more than 50% of miscarriages is still unknown (4). Among the factors affecting RPL, genes responsible for fetal growth and development have considerable roles.
Single nucleotide polymorphisms (SNPs) in several genes for proteins or fetal development are associated with RPL. Phospholipid hydroperoxide glutathione peroxidase (GPX4) is involved in glutathione synthesis (5) and may be associated with RPL in women; however, the exact mechanism is not clear. Catechol-O-methyl transferase protein (COMT) is now increasingly reported to have a role in pregnancy. In particular, recent reports have noted its expression in fetal membranes. This enzyme is also active both in the placenta and in the decidua (6). Reports have demonstrated that hypertension in pregnancy can reduce the placental COMT activity (7). Recent studies have shown that women with severe preeclampsia have reduced placental COMT protein expression (8). Besides, micro RNAs (miRs) are among the newly discovered genetic factors associated with RPL. miRs are small non-protein coding RNAs causing gene silencing at the post-transcription level. They target about one-third of human gene expressions and also play a critical role in various gestational events such as implantation, ectopic pregnancy, pre-term labor, pre-eclampsia, low birth weight, RPL, and anti-inflammatory responses. Recent investigations on women with RPL have shown that the overexpression of miR-133a can lead to a decline in human leukocyte antigen G (HLA-G) translation by joining to unranslated regions (UTR-3'). This protein has an essential role in the mother's immune system tolerating the fetus, and the down-regulation of this protein leads to RPL (9). miR-125 is a group of miRs which consists of two subgroups: miR-125a and miR-125b. The expression of these miRs occurs in different organs like the stomach, liver, lungs, rectum, breast, prostate, and ovaries. Research on these organs suggests that the malignant process may be controlled by preventing gene expression. Embryo development may also be affected by low secretion of the mother’s miR-125a polymorphism. Previous study have also shown that women with RPL have downregulation of 41 miRs and overexpression of four miRs (10). Although the critical role of miR-125a in RPL is unclear, recent evidence has proposed that its down-regulation may induce RPL in multiple ways (11).
By considering the wide range of miR activities in messaging processes, especially in the female genital system, studies on these regulatory RNAs seem to be necessary. The importance of miscarriage necessitates more comprehensive studies on the role of miRs in RPL. Therefore, the current study aimed to use the SNaPshot method to investigate the association between polymorphisms of pre-miR-125a (rs41275794), pre-miR-10a (rs3809783), pre-miR-323b (rs56103835), GPX4 (rs4680) and GPX4 (rs713041), and RPL in Iranian women with idiopathic recurrent miscarriage.
2. Materials and Methods
2.1. Sample collection
In the present case-control study, 116 women with idiopathic RPL and 89 healthy women as controls who had previously had at least two successful pregnancies (a total of 205 individuals) were included and referred to a Baqiyatallah hospital (Tehran, Iran) for sampling during 2017-2019.
A diagnosis of idiopathic RPL was confirmed by the expert gynecologist after evaluation of hormonal levels, immunological and hematological deficiencies, karyotyping, and uterine anomalies. Uterine anomalies were evaluated through hysterosalpingography, sonhysterography and hysteroscopy techniques. Healthy controls visited the hospital for checkups and no abnormalities were found in their physical or laboratory examinations. They did not have any history of RPL or previous related medical diseases.
Inclusion criteria were: (i) women with idiopathic RPL; and (ii) who had undergone complete clinical and pathologic examinations. Cases who met the following criteria were excluded from the study: (i) history of autoimmune or hormonal diseases; (ii) chromosomal aneuploidies or genetic abnormalities; and/or (iii) immunological or hematological disorders.
2.2. DNA extraction
In the first examination, after an overnight fast, 5 ml peripheral blood samples were collected and placed in the EDTA-containing tubes. Then, the rapid genomic DNA extraction (RGDE) method was used to extract genomic DNA (12).
Subsequently, a Nanodrop ND-1000 spectrophotometer (Thermo Sci., Newington, NH, USA) was used to determine the quantity and quality of the extracted DNA. A fraction of each DNA sample was placed on 1.5% denaturing agarose gel stained with Red Safe and was then subjected to electrophoresis to assess the DNA quality. Following that, the isolated DNA was stored at -20°C until further analysis.
2.3. Uniplex polymerase chain reaction (PCR) amplification
After DNA extraction, the functional SNPs pre-miR-125a (rs41275794), pre-miR-10a (rs3809783), pre-miR-323b (rs56103835), GPX4 (rs4680) and GPX4 (rs713041) were selected in each gene for genotyping based on published data and Primer-BLAST in the National Center for Biotechnology Information dbSNP database.
Oligo analyzer software was used to consider or predict the formation of dimer or secondary structures of primer. Primer sequences and their amplicon size are shown in table I. Gene runner software was used to design appropriate single-base extension (SBE) primers, as mentioned in table II. All individuals were analyzed for the presence of these polymorphisms (13).
2.4. Multiplex PCR amplification
To conduct the multiplex PCR, the following chemicals were prepared: 1 μM of each primer, 0.17 mM of dNTPs, 1.5 mM of MgCl2, one unit of Taq DNA polymerase, and 20 ng of template DNA in a 20 μl volume flask containing 10 X PCR buffer. The three conditions applied to the thermal cycle were: 1) 10 min at 95°C; 2) 30 sec at 95°C for 35 cycles; 3) 30 sec at 60°C, one min at 72°C and a final extension of 10 min at 72°C. Then, 3 μl of PCR product was run in 12% polyacrylamide gel to check for the quality and yield of the multiplex PCR. Next, the excess dNTPs and primers were removed by treating the 3 μl of remaining PCR products with one unit (1 μl) of exonuclease I. Then, three units of shrimp alkaline phosphatase and two units of exonuclease I (ExoSAP/USB, Cleveland, OH, USA) were used to purify the remaining PCR products, which was followed by incubation for 15 min at 37oC and 15 min at 80oC to remove excess dNTPs and primers, and to deactivate the Taq enzyme (14).
2.5. SNaPshot method
The SNaPshot assay applied five SBE primers which simultaneously annealed adjacent to SNP variant in our reaction. As shown in table II, poly (dC) tails of different lengths were attached to the SBE primers.
The SNaPshot solutions for multiplex PCR contained PCR product (1 μl), master mix SNaPshot (1 μl), each SBE primer (1 μl) and DDW (2 μl) in a final volume of 5 μl. The three thermal protocols applied to the SNaPshot were: two min at 96°C; 10 sec at 96°C; five sec at 50°C and 30 sec at 60°C. dNTPs were removed from the second PCR product by treating the remaining PCR products with three units (1 μl) of shrimp alkaline phosphatase and then the following two thermal protocols: 45 min at 37°C and 15 min at 75°C. An ABI PRISM 3130 Genetic analyzer was used to conduct capillary electrophoresis after which Gene Mapper idx, Ver 4.0 (Applied Biosystems, Foster City, CA, USA) was used to analyze the results (14).
2.6. Ethical considerations
This case-control study was approved by the institutional review board and ethical consideration committee of Baqiyatallah University of Medical Sciences (Code: IR.BMSU.REC.1398.054). Written informed consent was obtained from all individuals.
2.7. Statistical analysis
Statistical analysis was applied on SPSS v. 21 statistics software. Data were showed as percentages of the mean or allele frequency. Pearson’s chi-square test was used to calculate inter-group significance. As well as, the homozygous and heterozygous genotypes of each group were unified as carriers and the odds ratios (OR) and 95% confidence intervals (CI) were obtained. Statistical significance was set at p < 0.05.
3. Results
The present study was performed on 116 cases with RPL and 89 healthy women without any history of RPL. There was no significant difference between their ages (38.93 ± 10.45 years. vs. 41.31 ± 10.74 years, respectively; p = 0.40). After collection of blood samples, the DNA was extracted with the rapid genomic DNA extraction method, as shown in figure 1.
The Uniplex and Multiplex PCR techniques were applied to detect G>A polymorphism in pre-miR-125a (rs41275794), A>T polymorphism in pre-miR-10a (rs3809783), T>A polymorphism in pre-miR-323b (rs56103835), A>G polymorphism in GPX4 (rs4680) and C>T polymorphism in GPX4 (rs713041) genes, as shown in figure 2 and figure 3, respectively.
After multiplex PCR, the SNaPshot kit and specific SBE primers were used to detect rs713041, rs4680, rs41275794, rs3809783, and rs56103835 SNPs with the genetic analyzer machine, as shown in figure 4.
The comparison of SNP genotype distribution for all of the studied genes in the control and case groups is summarized in table III. The frequencies of AG, GG and AA genotypes in the polymorphism rs41275794 of pre-miR-125a (p < 0.001) and of TT, AT and AA genotypes in the polymorphism rs3809783 of pre-miR-10a (p = 0.04) were calculated for the women in the case and control groups. The findings showed significant differences and an association of these polymorphisms with the symptoms of recurrent miscarriage (p < 0.05). There was no significant difference (p = 0.09) in the frequency of TT, AT and AA genotypes in the GPX4 gene (rs713041 polymorphism) between the case and control groups. The frequency of TT, AT and AA genotypes in the miR-323b gene (rs56103835 polymorphism) showed no significant difference (p = 0.14) between the case and control groups. Also, no significant difference (p = 0.21) was observed in the frequency of GG, AG or AA genotypes in the COMT gene (rs4680 polymorphism) between the two groups (Table III). In the allele frequency study for alleles A and G of polymorphism rs41275794, a significant difference was found between the case and control groups (p < 0.001). Similarly, there was a significant difference (p = 0.02) in the frequency of alleles A and T of polymorphism rs3809783 between the two groups of cases and controls (Table IV). The assessment of risk at a 95% confidence level for allele A relative to G in the rs41275794 polymorphism showed that by changing the G allele to A, the risk of the RPL phenotype was increased by 1.912 times. Similarly, for the rs3809783 polymorphism, the changing of allele A to T increased the risk of the RPL phenotype by 1.704 times, showing a greater association of risk. To evaluate the genotypic inheritance of significant polymorphisms, we studied and calculated the risk of RPL in three models: dominant, recessive, and obvious models. The results showed that the dominant hereditary model had a significant correlation based on the genotypic frequency of the variant (Table V).
4. Discussion
RPL is a clinical challenge that affects approximately 1-5% of women of childbearing age (11). The underlying mechanism of RPL in 37-79% of cases is unclear. Recent evidence has revealed that molecular and genetic disorders, including various polymorphisms, may be the main cause of RPL in these cases (15).
In this research, we compared for the first time five polymorphisms of pre-miR-125a (rs41275794), pre-miR-10a (rs3809783), pre-miR-323b (rs56103835), GPX4 (rs4680) and GPX4 (rs713041) between healthy women and women with RPL using the SNaPshot genotyping method. Our findings showed a significant difference in the miR-125a gene polymorphism rs41275794 between the control and case groups. GG was found to be the most common genotype in the control subjects (60.7%), while GA (50.9%) was the most common genotype in the case group. These data indicated that the polymorphism in miR-125a may be a risk factor for RPL among the Iranian female population. In our study, allele A was significantly associated with an increased risk of RPL indicating that this variant allele in pre-miR125a can change precursor and mature miR-125a expression.
Our data revealed that the dominant genotype was significantly associated with an increased risk of RPL. In other study, expression of mature miR-125a in the CC/AA genotype was ~2.1-fold higher than that in the TT/AA genotype (16). They found that the A>T haplotype in rs41275794 affected miR-125a expression and caused RPL in women. In another study, Li and colleagues compared the expression pattern of six miRs, including miR-125a, miR-125b, miR-34a, miR-155, miR-24 and miR-141 between healthy women and those with RPL. They found that miR-125a, miR-125b, miR-34a, miR-155 and miR-141 were significantly overexpressed in women with RPL while miR-24 was significantly down-regulated in these cases, compared to the healthy group. Additionally, they showed that the A>T polymorphism in pri-miR-10a (rs3809783) was significantly associated with an increased risk of RPL in women.
More studies are being conducted nowadays to find relationships between SNPs and the miRNA system's critical components, biological functions and disease (4, 16, 17). The current study indicated that Iranian women with rs3809783 A>T in pri-miR-10a were significantly associated with an increased risk of RPL. This study’s findings also suggested that the expression level of the miR-10a target gene, Bim, might have been reduced because of miR-10a rs3809783. Additionally, this study showed a significant correlation between the miR-10a rs3809783 polymorphism and the risk of RPL. AT was the most common genotype in the case group while the control group most commonly exhibited the AA genotype. The frequency of the T allele in the case group (30.2%) was significantly greater than that of the control group (20.2%). On the other hand, the control group had a significantly higher frequency (79.8%) of the A allele compared to the case group (69.8%).
The current study also found a significant association between the dominant genotype and an increased risk of RPL. Our data showed that polymorphisms in miR-125a and miR-10a can be associated with the risk of RPL in women; however, further studies are needed to confirm this. Several studies have considered the role of different miRs in RPL (17, 18). In a case-control study among Korean women, RPL cases had significantly different frequencies of the miR-499TC + CC and miR-196a2CC genotypes than the control group. It was shown that these two polymorphisms and their combinations were significantly associated with the incidence of idiopathic RPL in the study population (19). In addition, in another study conducted in north Indian females show the miR-196a2CC, miR-499TC + CC, and miR-196a2CC/miR-499TC + CC were associated with an increased risk of RPL (20).
Furthermore, cumulative gene risk score analysis were calculated mentioning oxidative stress-related genes suggested that, common variations in individual genes like GPX4 and to a greater extent in the multi-locus model (ABCB1, GPX4, GPX4, OGG1 genes) are associated with RPL (21). In the current study, we did not find a significant difference in miR-323b and GPX4 polymorphisms between the healthy and case groups. Previous studies have reported that SNPs in the GPX4 gene significantly affect enzyme activity (6, 22, 23). The evidence shows that pre-eclampsia pathogenesis is associated with SNPs in this gene (6). A recent cohort study on a Korean population of 164 women with preeclampsia and 182 normotensive cases showed that preeclampsia was associated with the Val108/158Met polymorphism (6). Nevertheless, the variation in GPX4 enzyme activity is likely not due to the Val108/158Met polymorphism alone. Additionally, various vascular health impairments can be due to reduced levels of estrogen metabolites, such as 2-methoxyestradiol, because of a decrease in GPX4 activity (6, 24). Hypoxia-inducible factor-1a (HIF-1a) plays an essential role in angiogenesis and studies have found that this factor can be suppressed by anti-angiogenic 2-methoxyestradiol (6). Additionally, this transcription factor can induce genes such as soluble fms-like tyrosine kinase-1 (sFlt-1) which facilitates the adaptation and survival of cells during low-oxygen levels. Studies have also shown the stimulation of trophoblast invasion due to the combination of 2-methoxyestradiol and low levels of oxygen (25-27). Placental pathology can also be affected throughout the different stages of pregnancy if the homeostasis and regulation of GPX4 and 2-methoxyestradiol are disrupted. Also, it has been hypothesized that hypoxia-driven trophoblast invasion and vascular remodeling can be disrupted with a premature increase in 2-methoxyestradiol, thereby contributing to preeclampsia pathogenesis (6, 25). Decreased GPX4 activity in late pregnancy can cause reduced levels of 2-methoxyestradiol and thus reduce HIF-1a inhibition, resulting in possible vascular pathology and inflammation (6). Therefore, future studies on the expression and activity of the GPX4 enzyme throughout pregnancy can pinpoint the effect of GPX4/2-methoxyestradiol at different stages of miscarriage.
Since there is no similar report on the frequency of selected polymorphisms in the Iranian female population, one of the aims of the current study was to investigate the possible role of these polymorphisms in the occurrence of spontaneous miscarriage. A strength of this study was the careful selection of eligible individuals within the case group, i.e. women who had had at least two recurrent miscarriages, were completely healthy in terms of anatomical and cytogenetic studies and did not have any significant clinical problems in their medical records. The findings indicated that genetic polymorphisms can influence RPL, but the lack of significant differences in the GPX4, COMT and miR-323b polymorphisms in the Iranian females also reinforced the hypothesis that a single miR or polymorphism is unlikely to be solely responsible for RPL. One of the limitation of our study was the diversity of ethnicity in the Iranian population of the individuals who participate in investigation and lack of indigenous genetically database which harbor allele frequencies and other genomic data for estimation calculation samples collection and etc.
It is recommended that the expression of these polymorphisms in the blood as well as in the uterine tissue of women with RPL be measured to create a better understanding of the molecular mechanisms leading to miscarriage and to compare different conditions. Moreover, this research can be utilized to provide a clearer picture of the mechanisms involved in causing a miscarriage in women and to analyze a complete profile of polymorphisms involved in the various pathways leading to miscarriage.
Bioinformatics studies on the secondary and tertiary protein structures of selected genes using related databases can provide more comprehensive information on recurrent miscarriages.
5. Conclusion
In conclusion, our data illustrated that a polymorphism in pre-miR-125a (rs41275794) and in pre-miR-10a (rs3809783) can be associated with RPL in women. SNaPshot was shown to be a suitable method for investigation of the role of polymorphisms in the incidence of health conditions.
Acknowledgments
We thank Dr. Sayed Mostafa Hosseini from Baqiyatallah University of Medical Sciences for their technical support and Mahshid Shahrzad Kavkani for their assistance. The author(s) received no financial support for the research, authorship, and/or publication of this article.
Conflict of Interest
The authors declare that there is no conflict of interest.
Type of Study: Original Article |
Subject:
Cellular and Molecular Biology of Reproduction
References
1. Youssef A, Lashley L, Dieben S, Verburg H, van der Hoorn ML. Defining recurrent pregnancy loss: associated factors and prognosis in couples with two versus three or more pregnancy losses. Reprod BioMed Online 2020; 41: 679-685. [DOI:10.1016/j.rbmo.2020.05.016] [PMID]
2. Azani A, Hosseinzadeh A, Azadkhah R, Poursadegh Zonouzi AA, Poursadegh Zonouzi A, Aftabi Y, et al. Association of endothelial nitric oxide synthase gene variants (-786 T>C, intron 4 b/a VNTR and 894 G>T) with idiopathic recurrent pregnancy loss: A case-control study with haplotype and in silico analysis. Eur J Obstet Gynecol Reprod Biol 2017; 215: 93-100. [DOI:10.1016/j.ejogrb.2017.05.024] [PMID]
3. Turienzo A, Lledó B, Ortiz JA, Morales R, Sanz J, Llácer J, et al. Prevalence of candidate single nucleotide polymorphisms on p53, IL-11, IL-10, VEGF and APOE in patients with repeated implantation failure (RIF) and pregnancy loss (RPL). Hum Fertil 2020; 23: 117-122. [DOI:10.1080/14647273.2018.1524935] [PMID]
4. Li Y, Wang XQ, Zhang L, Lv XD, Su X, Tian S, et al. A SNP in pri-miR-10a is associated with recurrent spontaneous abortion in a Han-Chinese population. Oncotarget 2016; 7: 8208-8222. [DOI:10.18632/oncotarget.7002] [PMID] [PMCID]
5. Flohé L. The impact of thiol peroxidases on redox regulation. Free Radic Res 2016; 50: 126-142. [DOI:10.3109/10715762.2015.1046858]
6. Roten LT, Fenstad MH, Forsmo S, Johnson MP, Moses EK, Austgulen R, et al. A low COMT activity haplotype is associated with recurrent preeclampsia in a Norwegian population cohort (HUNT2). Mol Hum Reprod 2011; 17: 439-446. [DOI:10.1093/molehr/gar014] [PMID] [PMCID]
7. Barnea ER, MacLusky NJ, DeCherney AH, Naftolin F. Catechol-o-methyl transferase activity in the human term placenta. Am J Perinatol 1988; 5: 121-127. [DOI:10.1055/s-2007-999669] [PMID]
8. Kanasaki K, Palmsten K, Sugimoto H, Ahmad S, Hamano Y, Xie L, et al. Deficiency in catechol-O-methyltransferase and 2-methoxyoestradiol is associated with pre-eclampsia. Nature 2008; 453: 1117-1121. [DOI:10.1038/nature06951] [PMID]
9. Riesco-Eizaguirre G, Wert-Lamas L, Perales-Patón J, Sastre-Perona A, Fernández LP, Santisteban P. The miR-146b-3p/PAX8/NIS regulatory circuit modulates the differentiation phenotype and function of thyroid cells during carcinogenesis. Cancer Res 2015; 75: 4119-4130. [DOI:10.1158/0008-5472.CAN-14-3547] [PMID]
10. Barchitta M, Maugeri A, Quattrocchi A, Agrifoglio O, Agodi A. The role of miRNAs as biomarkers for pregnancy outcomes: A comprehensive review. Int J Genomics 2017; 2017: 8067972. [DOI:10.1155/2017/8067972] [PMID] [PMCID]
11. Hu Y, Huo ZhH, Liu ChM, Liu ShG, Zhang N, Yin KL, et al. Functional study of one nucleotide mutation in pri-miR-125a coding region which related to recurrent pregnancy loss. PloS One 2014; 9: e114781. [DOI:10.1371/journal.pone.0114781] [PMID] [PMCID]
12. Saremi MA, Saremi M, Tavallaei M. Rapid genomic DNA extraction (RGDE). Forensic Science International: Genetics Supplement Series 2008; 1: 63-65. [DOI:10.1016/j.fsigss.2007.12.001]
13. Nemati Vahedi S, Kheirkhah B, Malekirad AA. Association between polymorphism of has-miR-125 (rs12976445) and susceptibility to idiopathic recurrent pregnancy loss in Iranian women. Iran J Obstet Gynecol Infertil 2020; 23: 37-48.
14. Daniel R, Santos C, Phillips C, Fondevila M, van Oorschot RAH, Carracedo A, et al. A SNaPshot of next generation sequencing for forensic SNP analysis. Forensic Sci Int Genet 2015; 14: 50-60. [DOI:10.1016/j.fsigen.2014.08.013] [PMID]
15. Moghbeli M. Genetics of recurrent pregnancy loss among Iranian population. Mol Genet Genomic Med 2019; 7: e891. [DOI:10.1002/mgg3.891] [PMID] [PMCID]
16. Hu Y, Liu ChM, Qi L, He TZh, Shi-Guo L, Hao CJ, et al. Two common SNPs in pri-miR-125a alter the mature miRNA expression and associate with recurrent pregnancy loss in a Han-Chinese population. RNA Biol 2011; 8: 861-872. [DOI:10.4161/rna.8.5.16034] [PMID]
17. Jung YW, Jeon YJ, Rah HCh, Kim JH, Shin JE, Choi DH, et al. Genetic variants in microRNA machinery genes are associate with idiopathic recurrent pregnancy loss risk. PloS One 2014; 9: e95803. [DOI:10.1371/journal.pone.0095803] [PMID] [PMCID]
18. Rah HCh, Chung KW, Ko KH, Kim ES, Kim JO, Sakong JH, et al. miR-27a and miR-449b polymorphisms associated with a risk of idiopathic recurrent pregnancy loss. PLoS One 2017; 12: e0177160. [DOI:10.1371/journal.pone.0177160] [PMID] [PMCID]
19. Ahn DH, Rah H, Choi YK, Jeon YJ, Min KT, Kwack K, et al. Association of the miR‐146aC> G, miR‐149T> C, miR‐196a2T> C, and miR‐499A> G polymorphisms with gastric cancer risk and survival in the korean population. Mol Carcino 2013; 52 (Suppl.): 39-51.
20. Parveen F, Agrawal S. Recurrent miscarriage and micro-RNA among north Indian women. Reprod Sci 2015; 22: 410-415. [DOI:10.1177/1933719114529376] [PMID] [PMCID]
21. Khadzhieva MB, Lutcenko NN, Volodin IV, Morozova KV, Salnikova LE. Association of oxidative stress-related genes with idiopathic recurrent miscarriage. Free Radic Res 2014; 48: 534-541. [DOI:10.3109/10715762.2014.891735] [PMID]
22. Meplan C, Dragsted LO, Ravn-Haren G, Tjønneland A, Vogel U, Hesketh J. Association between polymorphisms in glutathione peroxidase and selenoprotein P genes, glutathione peroxidase activity, HRT use and breast cancer risk. PLoS One 2013; 8: e73316. [DOI:10.1371/journal.pone.0073316] [PMID] [PMCID]
23. Méplan C, Crosley LK, Nicol F, Horgan GW, Mathers JC, Arthur JR, et al. Functional effects of a common single-nucleotide polymorphism (GPX4c718t) in the glutathione peroxidase 4 gene: Interaction with sex. Am J Clin Nutr 2008; 87: 1019-1027. [DOI:10.1093/ajcn/87.4.1019] [PMID]
24. Dubey RK, Imthurn B, Jackson EK. 2-Methoxyestradiol: A potential treatment for multiple proliferative disorders. Endocrinology 2007; 148: 4125-4127. [DOI:10.1210/en.2007-0514] [PMID]
25. Lee SB, Wong AP, Kanasaki K, Xu Y, Shenoy VK, McElrath TF, et al. Preeclampsia: 2-methoxyestradiol induces cytotrophoblast invasion and vascular development specifically under hypoxic conditions. Am J Pathol 2010; 176: 710-720. [DOI:10.2353/ajpath.2010.090513] [PMID] [PMCID]
26. Shen Zh, Wu Y, Chen X, Chang X, Zhou Q, Zhou J, et al. Decreased maternal serum 2-methoxyestradiol levels are associated with the development of preeclampsia. Cell Physiol Biochem 2014; 34: 2189-2199. [DOI:10.1159/000369662] [PMID]
27. Perez-Sepulveda A, España-Perrot PP, Norwitz ER, Illanes SE. Metabolic pathways involved in 2-methoxyestradiol synthesis and their role in preeclampsia. Reprod Sci 2013; 20: 1020-1029. [DOI:10.1177/1933719113477483] [PMID] [PMCID]
Send email to the article author
| Rights and permissions | |
 |
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License. |
