FANCA Polymorphism Is Associated with the Rate of Proliferation in Uterine Leiomyoma in Korea

Uterine leiomyomas are the most common benign gynecologic tumors. This study was aimed to identify single nucleotide polymorphism of Fanconi anemia complementation group A (FANCA), associated with the rate of proliferation in uterine leiomyomas. In vitro study of patient-derived primary-cultured leiomyoma cells and direct sequencing of fresh frozen leiomyoma from each subject was conducted. Leiomyomas obtained from 44 patients who had underwent surgery were both primary-cultured and freshly frozen. Primary-cultured leiomyoma cells were divided into, according to the rate of proliferation, fast and slow groups. Single nucleotide polymorphism (SNP) of FANCA were determined from fresh frozen tissues of each patient using direct sequencing. Direct sequencing revealed a yet unidentified role of FANCA, a caretaker in the DNA damage-response pathway, as a possible biomarker molecule for the prediction of uterine leiomyoma proliferation. We identified that rs2239359 polymorphism, which causes a missense mutation in FANCA, is associated with the rate of proliferation in uterine leiomyomas. The frequency of C allele in the fast group was 35.29% while that in slow group was 11.11% (odds ratio (OR) 4.036 (1.176–13.855), p = 0.0266). We also found that the TC + CC genotype was more frequently observed in the fast group compared with that in the slow group (OR 6.44 (1.90–31.96), p = 0.0227). Taken together, the results in the current study suggested that a FANCA missense mutation may play an important regulatory role in the proliferation of uterine leiomyoma and thus may serve as a prognostic marker.


Introduction
Uterine leiomyomas, generally known as fibroids, are the most common benign gynecologic tumors. Approximately 20 to 40% of women of reproductive age are reported to have these benign tumors [1]. Although benign in their nature, leiomyomas are socioeconomically and clinically important, since 5 to 10% are estimated to be associated with infertility and 1 to 3% to cause infertility with other possible causes of infertility excluded [2]. With the increasing average age of women at the first childbirth, infertility has become a paramount issue both socioeconomically and clinically [3]. Due to the absence of uterus-sparing medical therapy, the treatment of choice of uterine fibroids is surgical removal via either hysterectomy or myomectomy. Patients with uterine leiomyomas, once diagnosed, have to go through regular check-ups for the progression of the growth of leiomyomas. Sometimes, leiomyoma would unexpectedly grow exponentially in the period of just a month, resulting in unwanted loss of fertility in women of reproductive age. Other times, it remains the same for years. Currently, no therapeutic tool that can predict the growth patterns of leiomyoma is available. Thus, the unpredictability of the growth patterns of leiomyoma is the utmost important unmet medical need regarding the treatment of uterine leiomyoma. Considering the high incidence of leiomyoma, absence of effective medical treatment, and its estimated contribution to infertility, the prognosis of leiomyoma, including growth rate and size variability, is of great importance for women of reproductive age.
Fanconi anemia complementation group A (FANCA), one of Fanconi anemia complementation group (FANC) family genes, codes for a protein that is part of the Fanconi anemia (FA) complex. FA complex mediates the FA pathway which maintains chromosome stabilities via interstrand DNA cross-link repair [4]. FA core complex proteins recruit DNA repair molecules, such as BRCA1, to the area of DNA damage so the error can be repaired. The FA pathway acts, using components of other known DNA repair processes, when DNA damages inhibit replication. Disorders or mutations in this gene cause an autosomal recessive disorder, the main features of which are aplastic anemia in childhood, multiple congenital abnormalities, susceptibility to leukemia and other cancers [4,5]. Evidence indicates that polymorphisms or mutations in the FANCA gene are related to cancer survival and susceptibility [6][7][8][9]. Among numerous polymorphisms of FANCA, rs2239359, a missense single nucleotide polymorphism (SNP) changing serine to glycine, is known to be associated with premature ovarian failure in Korean women [10].
In this study, we discovered the association of rs2239359 polymorphism with the rate of proliferation in leiomyoma in Korea.

Study Subjects
Uterine leiomyoma tissues were obtained from patients who underwent myomectomy or hysterectomy at Dongsan medical center (Daegu, Korea) between 2010-2017. Leiomyoma tissues from a total of 44 patients who did not receive preoperative hormonal therapy, including ulipristal acetate and GnRH agonist, were used in this study. Patients in the postmenopausal period were excluded. Written informed consent was obtained from all patients and the study was approved by the institutional review board of Dongsan Medical Center (Daegu, Korea; approval no. 09-156).

Isolation of Uterine Leiomyoma Cells and Culture
The tissue was washed twice in cold phosphate buffered saline before being minced into 5 mm pieces in a sterile culture dish. The minced tissue pieces were transferred into 50 mL conical tubes containing Hanks's balanced salts (Sigma-Aldrich, St Louis, MO, USA) supplemented with 25 mM HEPES (4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid), 100 units/mL antibiotics, 1.5 mg/mL collagenase IV (Sigma-Aldrich) and 0.2 mg/mL of DNase I (Roche Diagnostics, Mannheim, Germany). The tubes were kept at 37 • C in a water bath with gentle agitation for 3 h. Undigested tissues were filtered, and the cells were centrifuged at 1000 rpm for 5 min. The pellet was rinsed once with Hanks's balanced salts and dispersed in a complete medium composed of DMEM (Dulbecco Modified Eagle Medium)/F12 with 10% fetal bovine serum, 1% penicillin/streptomycin. The cells from the second passage were used for the experiments. To assess doubling time, 3 × 10 5 cells/well were plated in 100 mm dishes.

Doubling Time Measurement
Cells were seeded in 100 mm dishes, with 2 × 2 mm crossing lines drawn at the bottom, at a concentration of 10,000 cells/well with a medium containing 10% fetal bovine serum and 1% antibiotics.
Doubling time was determined with the following formula: Doubling time (DT) = TIn2/In(Xe/Xb), where T is incubation time, Xb is the cell number at the beginning of the incubation, and Xe is the cell number at the end of the incubation. Primary-cultured leiomyoma cells were then divided into two groups, fast (DT < 144 h) and slow (DT > 144 h), according to doubling time.

Preparation of Genomic DNA and Identification of SNP
Genomic DNA was isolated from leiomyoma tissue of each subject using GenExTM Tissue, Sx (Geneall, Seoul, Korea) according to the manufacturer's protocol. The genotypes of rs2239359 SNP from each of subject were determined by SolGent Co., Ltd. (Daejeon, Korea), using a BigDye Terminator v3.1 cycle sequencing kit (Thermo Fisher Scientific, Seoul, Korea) and an ABI3730XL DNA analyzer (Thermo Fisher Scientific, Seoul, Korea). The primers for PCR and sequencing were CTCAAGCAACATTACCTCAG (forward) and CCATCTCCTGACCTCGTGAT (reverse) and the genotype was determined using Vector NTI-Contig Express (Thermo Fisher Scientific, Seoul, Korea). Hardy-Weinberg equilibrium (HWE) was evaluated by chi-square test using Haploview 4.2 (HUGO Gene Nomenclature Committee: http://www.genenames.org/ (accessed on 12 June 2018). Each of rs2239359 genotypes is shown in Figure 1. and Xe is the cell number at the end of the incubation. Primary-cultured leiomyoma cells were then divided into two groups, fast (DT < 144 h) and slow (DT > 144 h), according to doubling time.

Preparation of Genomic DNA and Identification of SNP
Genomic DNA was isolated from leiomyoma tissue of each subject using GenExTM Tissue, Sx (Geneall, Seoul, Korea) according to the manufacturer's protocol. The genotypes of rs2239359 SNP from each of subject were determined by SolGent Co., Ltd. (Daejeon, Korea), using a BigDye Terminator v3.1 cycle sequencing kit (Thermo Fisher Scientific, Seoul, Korea) and an ABI3730XL DNA analyzer (Thermo Fisher Scientific, Seoul, Korea). The primers for PCR and sequencing were CTCAAGCAACATTACCTCAG (forward) and CCATCTCCTGACCTCGTGAT (reverse) and the genotype was determined using Vector NTI-Contig Express (Thermo Fisher Scientific, Seoul, Korea). Hardy-Weinberg equilibrium (HWE) was evaluated by chi-square test using Haploview 4.2 (HUGO Gene Nomenclature Committee: http://www.genenames.org/ (accessed on 12 June 2018). Each of rs2239359 genotypes is shown in Figure 1.

Statistical Analysis
Data were analyzed with SPSS software, version 19.0 (IBM Analytics, Seoul, Korea). Differences in continuous variables were determined using parametric Student's t-test. Categorical variables were compared with a chi-square test or Fisher's exact test. The association of rs2239359 with the proliferation of leiomyoma was analyzed using binary logistic regression and adjusted for age, body mass index (BMI), smoking, drinking, recurrence, hypertension, and the presence of diabetes mellitus. Odds ratios (ORs) with 95% confidence intervals (95% CI) were calculated. To investigate whether the rs2239359 polymorphism is associated with the clinical parameters of the study subjects, we performed a statistical analysis using a general linear model adjusted for age, BMI, smoking, and drinking status. Statistical significance was set at p < 0.05.

Anthropometrics, Comorbidities, and Blood Chemistry
The general anthropometric characteristics in the slow and fast groups are shown in Table 1. Body weight, height, and age were not different between the two groups. No differences in the number and size of leiomyomas were observed between the two groups. Also, the presence of comorbidities (hypertension, diabetes mellitus, thyroid disease, and hyperlipidemia) was not different. We did not find any association between factors in blood chemistry and the rate of proliferation of leiomyoma (Table 2). Of note, hematocrit showed a tendency to decrease in the fast group (p = 0.085).

Statistical Analysis
Data were analyzed with SPSS software, version 19.0 (IBM Analytics, Seoul, Korea). Differences in continuous variables were determined using parametric Student's t-test. Categorical variables were compared with a chi-square test or Fisher's exact test. The association of rs2239359 with the proliferation of leiomyoma was analyzed using binary logistic regression and adjusted for age, body mass index (BMI), smoking, drinking, recurrence, hypertension, and the presence of diabetes mellitus. Odds ratios (ORs) with 95% confidence intervals (95% CI) were calculated. To investigate whether the rs2239359 polymorphism is associated with the clinical parameters of the study subjects, we performed a statistical analysis using a general linear model adjusted for age, BMI, smoking, and drinking status. Statistical significance was set at p < 0.05.

Anthropometrics, Comorbidities, and Blood Chemistry
The general anthropometric characteristics in the slow and fast groups are shown in Table 1. Body weight, height, and age were not different between the two groups. No differences in the number and size of leiomyomas were observed between the two groups. Also, the presence of comorbidities (hypertension, diabetes mellitus, thyroid disease, and hyperlipidemia) was not different. We did not find any association between factors in blood chemistry and the rate of proliferation of leiomyoma (Table 2). Of note, hematocrit showed a tendency to decrease in the fast group (p = 0.085).

Association of rs2239359 SNP in FANCA with the Proliferation Rate of Leiomyoma
To find the possible association between the missense SNP rs2239359 in FANCA and the proliferation rate of leiomyoma, we performed direct sequencing ( Figure 1). The rs2239359 SNP in FANCA satisfied the Hardy-Weinberg equilibrium (HWE) (p = 0.658), and minor allele frequency (MAF) was 0.205, which was slightly lower than in a previous report [10] and NCBI information [11]. Table 3 shows the allele and genotype distributions of rs2239359 in the fast and slow groups. In the fast group, C allele frequency was 35.29%, which was significantly higher than that in the slow group (11.11%, OR = 4.036 (1.176-13.855), p = 0.0266). The frequency of subjects with C allele in the fast group was 64.7%, which was also higher than 22.2% of the slow group (OR = 6.44 (1.90-31.96), p = 0.023).

rs2239359 in FANCA and Thyroid Function
Although we found no significant difference between thyroid functional markers (T3, TSH, and fT4) and the proliferation rate of leiomyoma (Table 2), we found an association of rs2239359 in FANCA with the serum fT4 level (Figure 2). The average serum fT4 level in subjects with C allele was 1.69 ng/dL, which was significantly higher than 1.26 ng/mL of subjects with the TT type (p = 0.040). TSH also showed an increased tendency in subjects with C allele (p = 0.061). T3 level showed no difference between the two groups. We observed no other association of rs2239359 in FANCA with any factors in blood parameters (Table 4).

Discussion
In this study, we found that FANCA missense SNP rs2239359 was associated with increased proliferation rate of leiomyoma in Korean women. We also found the association of rs2239359 with serum fT4 level.
FANCA is a major member of the Fanconi anemia (FA) core complex genes, mutations of which are responsible for most cases of FA, a heritable autosomal recessive disease characterized by high susceptibility to chromosomal breakage, leading to progressive bone marrow failure, congenital abnormalities, and cancer predisposition [12]. Clinical manifestations of FA are highly variable, implicating genetic heterogeneity. FANCA mutant mice exhibited growth retardation, microphthalmia, and craniofacial malformations [13]. One of the various clinical manifestations of FA is impaired fertility. Women with FA present hypoplastic ovaries with menstrual irregularities, secondary amenorrhea, anovulation, and premature senescence of reproduction [14,15]. Mutant FANCA mice displayed reduced ability to reproduce and gonadal defects [13,16,17].
Of the numerous clinical manifestations that are related to FANCA mutations, genetic association of FANCA mutations with uterine leiomyoma, especially in relation to the proliferation rate, has not been reported. Given the extent of genetic heterogeneity of the FANCA mutation and its involvement in diverse cellular functions, finding the association of FANCA missense SNP rs2239359 with the proliferation rate of leiomyoma in the current study was not surprising. To support a highly possible role of FANCA in the growth of leiomyoma, a study reported protein-protein interactions between FANCA and BRCA1 in leiomyomas [18]. FA core complex recruits key molecules of DNA repair system through the FA/BRCA pathway [9]. BRCA1 also participates in the homologous recombination process and double-strand break repair. Therefore, the loss of activity of FANCA could result in DNA damage and contribute to the differential growth rate of leiomyoma.
Further, a very recent study revealed that genes involved in genome stability are associated with the development of uterine leiomyoma [19]. They showed the association of ATM and TP53 with uterine neoplasia, possibly through genetic instability. Interestingly, like FANCA, ATM and TP53 are regulatory molecules involved in the DNA damage response [20]. They have also been found to be recurrently mutated in leiomyosarcoma [21].
FANCA is believed to induce the synthesis of GnRH [22], which in turn stimulates gonadotropin and sex hormones, including estrogen. Estrogen is a firmly established inducer of leiomyoma growth [23]. It also affects thyroid functions [24,25]. Thyroid dysfunction in patients with leiomyoma after administration of a GnRH agonist has been reported [26,27]. Given the result in this study and the above-mentioned information, there might be a possible, not yet identified relationship between FANCA and thyroid hormones in relation to the growth of leiomyoma. Further studies establishing the mechanistic details to determine the effect of FANCA on thyroid function in relation to the growth of leiomyoma will be of great importance.
This study has limitations such as sample size and absence of other SNP information positioned in FANCA. We also have not undertaken the diepoxybutane (DEB) and mitomycin-C (MMC) stress test to rule out FA. However, we excluded all the anthropomorphic associations among the patients to rule out FA. Nevertheless, we found the relation between rs2239359 of FANCA and proliferation of leiomyoma in Korean women. We also found the association of rs2239359 with the serum fT4 level. The significance of this study lies in the finding of FANCA SNP rs2239359 as the possible determining factor for the differential growth rate of leiomyoma. Since FANCA mutations are highly common in gynecological malignant tumors such as breast cancer, further studies including DEB and MMC stress test to determine the factors that render leiomyoma benign and biomarkers that predict the growth patterns of leiomyoma will be of great benefit.