- Open Access
Genomic and proteomic profiling I: Leiomyomas in African Americans and Caucasians
Reproductive Biology and Endocrinology volume 5, Article number: 34 (2007)
Clinical observations indicate that leiomyomas occur more frequently in African Americans compared to other ethnic groups with unknown etiology. To identify the molecular basis for the difference we compared leiomyomas form A. Americans with Caucasians using genomic and proteomic strategies.
Microarray, realtime PCR, 2D-PAGE, mass spectrometry, Western blotting and immunohistochemistry.
Using Affymetrix U133A array and analysis based on P ranking (P < 0.01) 1470 genes were identified as differentially expressed in leiomyomas compared to myometrium regardless of ethnicity. Of these, 268 genes were either over-expressed (177 genes) or under-expressed (91 genes) based on P < 0.01 followed by 2-fold cutoff selection in leiomyomas of A. Americans as compared to Caucasians. Among them, the expression E2F1, RUNX3, EGR3, TBPIP, ECM2, ESM1, THBS1, GAS1, ADAM17, CST6, CST7, FBLN5, ICAM2, EDN1 and COL18 was validated using realtime PCR low-density arrays. 2D PAGE coupled with image analysis identified 332 protein spots of which the density/volume of 31 varied by greater than or equal to 1.5 fold in leiomyomas as compared to myometrium. The density/volume of 34 protein-spots varied by greater than or equal to 1.5 fold (26 increased and 8 decreased) in leiomyomas of A. Americans as compared to Caucasians. Tandem mass spectrometric analysis of 15 protein spots identified several proteins whose transcripts were also identified by microarray, including 14-3-3 beta and mimecan, whose expression was confirmed using western blotting and immunohistochemistry.
These findings imply that the level rather than the ethnic-specific expression of a number of genes and proteins may account for the difference between leiomyomas and possibly myometrium, in A. Americans and Caucasians. Further study using larger sample size is required to confirm these findings.
It has been estimated that 70% of women have a life-long risk of developing leiomyomas, with symptomatic tumors accounting for a 1/3 of all the hysterectomies preformed annually in the United States alone. Such estimates are even higher among African Americans with leiomyomas that develop earlier, become larger and more symptomatic as compared to other ethnic groups [1–3]. Since leiomyomas develop during the reproductive years and regress with menopause, ovarian steroids are critical to their growth, however the molecular environment to explain these ethnic differences is unknown.
Various conventional and recent large-scale gene expression-profiling studies have provided valuable information regarding the molecular environment of leiomyomas [4–11]. Most of these studies were based on comparing leiomyomas with myometrium obtained from different ethnic groups without examining the influence of ethnicity on overall gene expression [9, 10, 12]. Using tissue microarray a recent immunohistochemistry study has reported the expression of a number of proteins in leiomyomas with some differences in their expression intensities in leiomyomas from African Americans as compared to other ethnic groups . In addition, an increased risk of developing leiomyomas, including among African Americans, has been examined in a number of genetic studies, although in a majority of the cases the evidence of genomic instability is inconsistent [14–17].
Considering the above information and availability of limited data differentiating leiomyomas microenvironment in African Americans in the present study we combined genomic and proteomic strategies comparing the molecular environment of paired leiomyoma and myometrium from African Americans with Caucasians. Furthermore, from microarray data sets we selected 15 differentially expressed genes from different functional categories rather than levels of their expression and validated their expression using realtime PCR in a microfluidic card format. Additionally, 15 protein spots identified by proteomic were isolated and subjected to tandem mass spectrometry (MS/MS) for identification of their content and verified two of these proteins by Western blotting and immunohistochemistry.
Portions of leiomyoma and matched myometrium were collected from six women, ranging from 29 to 38 years old who were scheduled to undergo hysterectomy for indications related to symptomatic leiomyomas. Three of the patients were Caucasians and three African Americans. These women were not taking any medication including hormonal therapy for the pervious 3 months prior to surgery and based on their last menstrual period and endometrial histology they were from early-mid secretory phase of the menstrual cycle. The myometrium used in this study was collected from regions distal from leiomyomas. All the leiomyomas were subserosal/intramural and 3 to 4 cm in diameter. These tissues were collected at the University of Florida affiliated Shands Hospital with prior approval from the Institutional Review Board without requiring to obtain informed consent. The tissues were snapped frozen and stored in liquid nitrogen until used.
Microarray and gene expression profiling
Small portions of the above tissues were used to isolate total RNA using Trizol (Invitrogen, Carlsbad, CA). RNA quality and yield was analyzed with Agilent 2100 Bioanalyzer (Agilent Technologies, Foster City, CA). RNA amplification was carried out as previously described , and following second-strand cDNA synthesis 5 μg of purified cDNA was reverse transcribed using Enzo BioArray high yield RNA transcript labeling kit (Affymetrix, Santa Clara, CA). After purifying the product, 20 μg of cRNA (0.5 μg/μl) was fragmented. The fragment then mixed with 300 μl of hybridization mixture and 200 μl of the mixture was hybridized to human U133A Affymetrix GeneChip which consists of 22,277 oligonucleotide probe sets representing 18,400 transcripts and variants, including 14,500 known genes. The chips were processed after meeting recommended criteria for use of the expression arrays as previously described  and according to manufactures protocol.
Microarray data analysis
The chips were scanned using Affymetrix Genepix 5000A scanner and the resulting images were assessed using Genepix software with a manual supervision to detect any inaccuracies. The net hybridization values were determined and subjected to global normalization and transformed expression values were subjected to Affymetrix Analysis Suite V 5.0 to remove any probe sets that were flagged as absent on all arrays using default settings. The expression value of the remaining probe sets was then subjected to unsupervised and supervised learning, discrimination analysis, and cross validation as previously described . After variation filtering, the coefficient of variation was calculated for each probe set across all chips, ranked and the expression values of the selected genes were then subjected to statistical analysis in "R" programming as previously described . The gene expression values having a statistical significance of p ≤ 0.01 (ANOVA, Turkey test) were selected and subjected to 2-fold cutoff. Hierarchical clustering and TreeView analysis was carried out and the selected genes were subjected to functional annotation and visualization using Database for Annotation, Visualization, and Integrated Discovery  software. The integrated GoCharts assigns genes to specific ontology functional categories based on selected classifications.
Low-density microfluidic card and realtime PCR
Total RNA isolated from the same paired leiomyomas and myometrium from African Americans and Caucasians used for microarray analysis was also subjected to realtime PCR using microfluidic cards TaqMan Low Density Arrays (LDAs) assessing the level of expression of 23 selected genes, and the house-keeping gene, G3PDH. These genes were selected based on their functional categories rather than the level of their expression from the list of differentially expressed genes identified in the present and our previous [6, 7, 19] microarray studies. The products of these 23 genes have been implicated in several biological processes involved in pathogenesis of leiomyoma and other fibrotic disorders. Pre-designed TaqMan probe and primer sets for these genes were selected from the on-line catalogue and were factory-loaded into customized 384-well LDA plates (Applied Biosystem, Foster City, CA). The LDA format was arranged on-line with each array containing two replicates for each 23-target gene and GAPDH. Since these LDAs were used to determine the expression of these genes in several different tissues for establishing the least variable housekeeping gene as standard prior to these experiments we utilized a factory-loaded LDA (384-well plate) specific for 12 housekeeping genes. The results indicated that GAPDH expression displayed the least variation among the reference genes in these tissues as well as the tissues used in the accompanied manuscript . Based on these results we selected GAPDH as a housekeeping gene during the development of the LDA and for normalization of the expression of other target genes on LDA. The exact locations and the sequences of the oligonucleotides used in all assays can be downloaded from the Applied Biosystems website  by selecting the Assays IDs. All genes span exonintron boundaries and cover the major transcript forms.
The LDAs were analyzed using the 7900HT system with a TaqMan LDA Upgrade (Applied Biosystems). In brief, 5 μl of single-stranded cDNA at final concentration of 2 μg starting RNA was combined with 45 μl water and 50 μl TaqMan Universal PCR Master Mix. The mixture was injected into selected sample ports and the cards were centrifuged twice for 1 min at 1200 rpm and sealed to prevent well-to-well contamination. Thermal cycling conditions were as follows: 50°C for 2 min, 95°C for 10 min, 97°C for 30s and 60°C for 1 min for 40 cycles. In all the experiments total RNA isolated from each sample were assayed in duplicate and one sample was used in all cards to assess inter- and intra-assay variability. Since each LDA card has a unique barcode the Sequence Detection System plate documents stores the information on the plate type, detector, sample/target gene configurations, thermal cycling conditions, data collection, and raw fluorescence data at each cycle. The data obtained from each card was analyzed using the comparative method and following normalization of expression values to GAPDH expression using relative quantity (RQ) and Sequence Detection Software 2.2.1 for automated data analysis according to the manufacturer's guidelines . For calculating the RQ of selected mRNA in leiomyomas compared with myometrium, or tissues from African Americans compared to Caucasians, myometrial RNA used in each card was designated as calibrator. To analyze the RQ amounts of selected genes as an effect of ethnicity, Caucasian myometrium was used as calibrators. The Ct value >35 was selected as cutoff for the absence of gene expression. The expression of 15 of these genes is reported here.
Proteomic: protein isolation and CyDye labeling
Total protein was isolated from 100 to 300 mg of paired frozen leiomyoma and myometrium from African Americans and Caucasians. The tissues were grinded in liquid nitrogen until fine powder, 5 ml of Trizol reagent and 5 to 15 μl of protease inhibitor cocktail (Sigma-Aldrich, St. Louis, MO) was added and centrifuged at 10,000 × g for 20 min at 4°C. Following phase extraction with chloroform the mixture was centrifuged at 5000 RPM at 4°C for 20 min, the inter- and bottom phases were extracted with 100% ethanol and clarified by centrifugation. The protein was precipitated with isopropanol overnight at -20°C and the pellets were washed twice with ice cold 80% ethanol, air dried and dissolved in 300 μl of CyDye labeling buffer (8 M urea, 2 M thiourea, 4% CHAPS, 20 mM Tris pH 8.5) and centrifuged at 40,000 × g for 30 min. The supernatants were transferred to a clean tube and adjusted to pH of 8.5 and protein concentration was determined using EZQ® (Invitrogen). Fifty μg of myometrial and leiomyomas protein mixtures was labeled with 400 pM of Cy3 and Cy5 respectively, for 30 min on ice in dark and the reaction was stopped and quenched with 1 μl of 10 mM Lysin.
For the first dimension, the Cy3- and Cy5-labeled proteins were combined with 450 μg each of unlabeled protein and 100 mM of DTT ampholyte buffer pH 3 to 11, to 0.5% and adjusted the final volume to 400 μl. An 18 cm of pH 3 to 11 immobilized pH gradient (IPG) strip was passively rehydrated with above sample solution overnight in the dark. Isoelectrofocusing (IEF) was carried in IPGphor (Amersham, GE Healthcare Bio-Sciences Corp. Piscataway, NJ) and the proteins co-migrated together was focused at 8000 V at 20°C for 60 kV hr. For the second dimension SDS-PAGE, an 8 to 16% precast Tris Glycine gel (BioRad, CA) was used. Prior to mounting, the IPG strips were reduced and equilibrated in 15 ml of 50 mM Tris-HCl pH 6.8, 6 M Urea, 30% (v/v) glycerol 2% (w/v) SDS 100 mM DTT for 15 min, then was alkylated with 15 ml of 50 mM Tris-HCl pH 6.8 6 M Urea 30% (v/v) glycerol 2% SDS 2.5% idoacetamide for 15 min. 2nd dimension was carried in the dark at 20°C for 4 to 5 hrs at 24 mA/gel. MW markers (New England Biolabs, Beverly, MA) were added to the 2D-PAGE.
Image acquisition and spot quantification
The gels were scanned with Typhoon 9400 Variable Mode Imager (Amersham) and following image acquisition the protein spots were identified, their volumes were determined and analyzed using computer-assisted gel analysis and Progenesis PG220 2D software (Nonlinear Dynamics, Durham, NC). The analysis allows automatic detection and quantification of protein spots as well as resizing, alignment, and matching between different 2D gel images. A database of all protein spots obtained from digital images was created, and integrated intensity of each spot was calculated based on spot area and optical density and expressed as percent volume, namely, spot volume/volumes of all spots resolved in the gel. The comparative program was used to geometrically correct spatial differences and allow comparison between the images. The information from comparison of the image databases was then used to establish an average gel. Each protein spot volume was normalized against the total protein spots volume present in the gel and the fold difference of each protein was calculated for the paired myometrium and leiomyoma sample. The list of protein spots with assigned identification number was prepared and arranged according to the fold difference and the values were statistically analyzed using ANOVA, Turkey test.
Protein identification by mass spectrometric analysis
The gels were Coomassie blue stained and the spots of interest were identified and manually excised with a disposable gel puncher and stored in 25% methanol until tyrptic digestion according to a protocol used by Proteomics Core Facility of the Inter-disciplinary Cores for Biomedical Research (ICBR) at the University of Florida. Peptides derived from tryptic digestion were analyzed by reversed-phase HPLC-tandem mass spectrometry on a hybrid quadrupole time-of-flight instrument (QSTAR XL, Applied Biosystems) equipped with a nano-electrospray source. Solvent delivery at 200 nl/min was provided by an integrated capillary HPLC system (Ultimate, LC Packings, Sunnyvale, CA) in which a 30 min gradient from 3% to 50% acetonitrile in 0.1% acetic acid was employed. Each information dependent acquisition (IDA) cycle consisted of a survey scan from m/z 400 to 1500 and three MS/MS scans obtained by collision-induced dissociation of ions that demonstrated the largest signal intensity at a given chromatographic time point. Survey and MS/MS scans were accumulated for 1 and 2 sec, respectively.
Tandem mass spectrometric data were searched against the IPI human protein database using the Mascot search algorithm. Carbamidomethylation of cysteine was included as a fixed modification whereas oxidation of methionine, deamidation of asparagine and glutamine, and pyro-glu formation from N-terminal glutamine or glutamic acid were included as variable modifications. Probability-based MOWSE scores that exceeded the value corresponding to p < 0.05 were considered for protein identification. After inputting the pI and MW of proteins Swiss-Prot Database was searched by TagIdent software (Swiss Institute of Bioinformatics, Basel, Switzerland) for matching proteins.
Western blot analysis
Briefly, small portions of paired leiomyoma and myometrial tissues were homogenized in homogenizing buffer containing protease inhibitor cocktail . The laysets were centrifuged, supernatants protein content was determined and an equal amount of proteins was subjected to SDS-PAGE, transfer by electroblotting and the blots were exposed to polyclonal antibodies generated against 14-3-3β and mimecan (Santa Cruz Biotechnology, Santa Cruz, CA), respectively and monoclonal antibody to β-actin (Sigma) to normalize for protein loading. The blots were exposed to corresponding HRP-conjugated IgG and visualized using enhanced chemiluminesence reagents (Amersham). The band intensity was determined using Kodak gel analysis software (Eastman Kodak, Rochester, New York).
For immunohistochemical localization of 14-3-3β and mimecan, fixed and paraffin embedded leiomyomas and myometrial tissue sections 3 to 5 μm thick were prepared. Following standard procedures, the sections were incubated with polyclonal antibodies at 5 μg of IgG/ml prepared in phosphate buffered saline, pH 7.4 . The sections were than exposed to biotinylated second antibodies and avidin horseradish peroxidase (Vector Laboratories, Burlingame, CA) and chromogenic reaction was developed using 3, 3'diaminobenzidine and counter sections were counter staining with hemotoxalin as indicated. Tissue sections incubated with normal IgG instead of the primary antibodies, or deletion of the primary antibodies during immunostaining served as controls.
Gene expression profiling of leiomyomas from African Americans and Caucasians
Gene expression values obtained from microarray analysis of paired leiomyoma and myometrium from African Americans and Caucasians were subjected to unsupervised and supervised learning. Based on statistical analysis with p ranking of P ≤ 0.01 we identified 1470 transcripts, or 8% of the genes on the array, as differentially expressed among these tissues regardless of ethnicity. Further analysis based on 2-fold cutoff resulted in identification of 268 genes with 177 genes over-expressed and 91 genes under-expressed in leiomyomas of African Americans as compared to Caucasians (Table 1; for the complete list of genes see Table 4, Additional file 1). Hierarchical clustering separated these genes into several distinctive clusters (Fig. 1) and functional pathway analysis indicated that majority of their products serve as transcriptional, translational and signal transduction mediators, cell cycle regulation, ECM turnover, cell-cell communication and metabolic activities etc.
Gene expression using TaqMan LDA
We used customized 384-well LDA cards to determine the expression of 15 genes in leiomyomas and myometrium, with each card containing two replicates of the target genes and GAPDH as a housekeeping gene (Fig. 2). Since these cards were also used to assess the expression of these genes in tissues reported in the accompanied manuscript  we first established the least variable housekeeping gene to use as a standard among these tissues. Using factory-loaded LDA cards representing 12 housekeeping genes GAPDH expression displayed the least variation among these tissues (data not shown). Based on these results we selected GAPDH as a housekeeping gene and along with 23 target genes for preparation of customized LDA for our study. The selection of the 23 genes was based on the result of the present and our previous microarray studies [6, 7, 19]. Figure 2 show the level of expression of E2F1, RUNX3, EGR3, TBPIP, ECM2, ESM1, THBS1, GAS1, ADAM17, CST6, CST7, FBLN5, ICAM2, EDN1 and COL18, representing transcription factors, cell cycle and apoptosis regulators, matrix remodeling and cell adhesion signaling, respectively in matched leiomyoma and myometrium from African Americans and Caucasians. The level of expression of these genes varied significantly in leiomyoma and myometrium as well as tissues from African Americans as compared to Caucasians (Fig. 2; P < 0.05), some displaying a similar pattern of expression observed with microarray.
Proteomic analysis of leiomyomas and myometrium
A representative of total protein isolated from a paired myometrium and leiomyoma, labeled with Cy3 or Cy5 respectively, and co-separated in 2D PAGE is shown in Figure 3A. Using image analysis 332 protein spots were identified and their density/volume determined with two spots (#107 and #192) later extracted for identification shown in graphic presentation (Fig. 3B). The average density of each protein spot in leiomyomas and myometrium was normalized against their density in one of the myometrium from a Caucasian patient. Of the 332 protein spots identified 28 and 31 spots were unique to leiomyomas and myometrium, respectively with the density/volume of 31 spots varied by at least 1.5 fold (12 were overexpressed and 19 underexpressed) in leiomyomas as compared to myometrium regardless of ethnicity. 109 protein spots identified in leiomyomas the average density/volume of 34 spots differed by at least 1.5 fold, of which 26 spots were overexpressed and 8 were underexpressed in African American as compared to Caucasians (Table 2). Mass spectrometric analysis of 15 selected spots resulted in identification of 137 proteins. A database search against Swiss-Prot resulted in identification of several proteins and the results are summarized in Table 3, with MS and MS/MS scores showing the reliability of identification. Among the content of these protein spots included several keratins, Annexin A1 and V, transgelin (SM22), galgizzarin (S100A11) and EF-hand domin-containing protein 2 (S100A12), vimentin, retinoic acid binding protein II, cofilin, several isoforms of 14-3-3 and mimecan. Transcripts corresponding to these proteins were also identified among differentially expressed genes in leiomyomas and myometrium of both ethnic groups (Table 3).
Western blot analysis and immunohistochemical localization
Western blot analysis indicated that leiomyoma and myometrium from both ethnic groups express 14-3-3β with a considerable variation among these tissues (Figs. 4). We were unable to detect mimecan immunoreactive protein in these tissue extracts by western blotting possibly because of the nature of the antibody for this application.
Immunohistochemically, 14-3-3β (Fig. 5A and 5B) and mimecan (Fig. 5C and 5D) were localized in leiomyoma and myometrial smooth muscle cells, connective tissue fibroblasts and vascular cells, with a considerable variability in their intensity among these cells within the same tissue (Figs. 5). Incubation of the tissue sections with normal rabbit and/or goat serum instead of the primary antibodies resulted in considerable reduction in staining intensity (Fig. 5E and 5F; the insert images show a higher magnification portion of these figures indicated by the arrows).
Using genomic and proteomic strategies the present study provided further insight into the molecular environments of leiomyoma and myometrium in Caucasians and African Americans. At genomic level 1470 genes or about 10% of the transcripts on the array were identified as differentially expressed in leiomyomas as compared to myometrium regardless of ethnicity. Of these genes 268 were identified as either over-expressed (177 genes) or under-expressed (91 genes) in leiomyomas of African Americans as compared to Caucasians. Hierarchical cluster analysis separated these genes into several clusters reflecting their tissue as well as ethnic-dependent and -independent association. However, the profile of some of the genes in these clusters displayed an area of relatedness between myometrium and leiomyomas within each ethnic group. This suggests that some of the differences in leiomyoma's gene expression might be attributed to differences in myometrial gene expression between the ethnic groups, as well as the differences in leiomyomas vs. myometrium. Functional analysis indicated that the majority of these genes serve as transcriptional, translational and signal transduction mediators, cell cycle regulation, ECM turnover, cell-cell communication and transport/metabolic activities.
At proteomic level we identified several proteins displaying both ethnic-dependent and -independent profiles, with 34, or 10% of total protein spots identified, displaying altered intensity in leiomyomas of African Americans as compared to Caucasians. Tandem mass spectrometry analysis of 15 selected protein spots revealed an array of proteins as part of their content (Table 2) of which the transcripts of some of them, or their related proteins, were among the genes identified in myometrium and leiomyomas. Although the genomic and proteomic results of our study in part confirmed the previous reports involving leiomyomas and myometrium [8–10, 12, 13, 19, 23, 24], previous microarray analysis did not consider the influence of ethnicity as part of their comparative assessment. We recognize that low sample size is a limitation of our study, however validating the expression of a large number of transcripts simultaneously in the same samples by quantitative realtime PCR (low density array) and identification and confirmation of two proteins, mimecan (osteoglycin) and 14-3-3β in paired leiomyoma and myometrium compensated for this limitation.
Of the genes functionally relevant to processes that contribute to characteristic and pathogenesis of leiomyomas was the identification of several components of TGF-β system and their elevated expression in leiomyomas of African Americans as compared to Caucasians. Based on this and our previous studies we consider over-expression of TGF-β and TGF-β receptors in leiomyomas of African Americans result in further amplification of their signaling, targeted transcription factors and downstream genes expression, which are already altered in leiomyomas as compared to myometrium [7, 19, 22, 25–28]. In contrast, leiomyomas of African Americans expressed lower levels of LTBP-1, a member of fibrillin super-family that associates and stores TGF-β into the ECM. LTBPs consist of LTBP-1 to LTBP-4  with LTBP-1 and -2 expression reported in leiomyomas [30, 31]. Commencing our study a recent report provided evidence for the expression of LTBP-1 and FBN-1 in leiomyomas of three different sizes with higher levels of expression in the medium-sized (3–5 cm) tumors compared with myometrium in the proliferative phase, while FBN-1 mRNA expression was size-independent . Although the results suggested that ECM turnover in leiomyomas might be size-dependent, the larger tumors undergo significant alteration in their tissue structure, specifically at the center that must be taken into account during such comparative analysis. Our results that LTBP-1 is expressed at lower levels in leiomyomas of African Americans as compared to Caucasians suggest that TGF-β upon secretion and activation becomes readily available for binding rather than being stored into the ECM, a mechanism that regulates local availability of many growth factors and cytokines, including TGF-β.
Plasmin, integrins, thrombospondin-1 (TSP-1), mannose 6-phosphate, and decorin are known to activate TGF-β [32–34] and their differential expression suggests the presence of regulatory mechanism that control TGF-β action in leiomyoma where it regulates the expression of several genes involved in tissue fibrosis [19, 25–28, 35]. Additionally, we identified TCF8/ZB-1, a member of zinc finger transcription factors that include Smad interacting protein 1 (SIP1/ZEB-2), as differentially expressed gene in leiomyomas of African Americans. TCF8 through interaction with Smad and recruitment co-activator (p300) and co-repressor (CtBP), is involved in ER and TGF-β receptors signaling [35–37]. Leiomyomas express these and several other transcription factors that interact with TCF8 including Runx, CITED, EGR1, EGR3, E2F1, Nurr77, c-Myc, Max, and Mad, whose expression is validated in this and our previous study [7, 19]. These transcription factors regulate the expression of group of genes involved in cell cycle regulation, apoptosis, angiogenesis and inflammation mediated by the action of TGF-β receptor and ER activation [38–43]. Because leiomyomas grow more rapidly in African Americans as compared to other ethnic groups, altered expression of these genes regulated by TGF-β and ovarian steroids could influence the outcome of this process.
A balance between cell growth and apoptosis is critical in progression of tissue fibrosis and the expression of several genes functionally related to these categories, including Bcl-XL, Bad, Bax, p27Kip1, p57Kip2, Gas1, Gas7, CST6, CST7, caspases, etc., were identified in leiomyomas and myometrium. Differential expression of some of these genes in leiomyomas of African Americans most likely contributes toward their rapid growth as compared to other ethnic group. Local apoptotic and non-apoptotic cell-death also results in regulation of multiple genes involved in inflammation, angiogenesis, fibrogenesis and tissue turn over. Among these genes involved in these processes are collagens, versican, fibromodulin, syndecan 4, TSP-1, tenascin-C, osteonectin/SPARC and ECM2, endothelial cell specific molecule-1 (ESM-1), ICAM2, EDN1, FZD7, β1 catenin, CST6, CST7 as well as several member of integrins, MMPs, TIMPs and ADAMs [44–49] that are expressed in leiomyomas [7, 26, 27] with altered expression in African Americans. The identification of 14-3-3β and mimecan as differentially expressed genes in leiomyomas of African Americans is also of interest. 14-3-3s are regulatory proteins that bind to a variety of cellular targets, including Raf kinase, cell cycle-dependent phosphatase Cdc25, pro-apoptotic protein Bad and many others proteins and are considered to regulate hypertrophic scar formation through regulation of MMPs expression [50–52]. The ability of 14-3-3 proteins to bind and regulate various oncogenic gene products as well as various tumor suppressor gene products points to a potential role in cancer As such 14-3-3 has been shown to inhibit TSC1/TSC2 complex functions and overexpression of either TSC1 or TSC2 in Hela cells has been reported to increase the expression of various 14-3-3 isoforms  suggesting that deregulated expression of 14-3-3 and TSC can be associated with leiomyoma pathogenesis. Mimecan is a member of small leucine-rich proteoglycans family, which includes lumican, fibromodulin, decorin and biglycan, and has been implicated in collagen fibrillogenesis, cellular growth, and migration [54–56]. Mimecan is secreted and proteolytically cleaved by a serine protease to form the 25–30 kDa form found in association with ECM of connective tissue and considered to play a structural role by maintaining the tensile strength and hydrated nature of the tissue [55, 56]. Mimecan is expressed in mouse pituitary gland and in human anterior pituitary gland and Pit-1 is reported to activate the human mimecan promoter through Pit-1 response element sites . As such mime can and 14-3-3 as well as other genes products in these categories may have a similar biological activity in leiomyoma pathogenesis, specifically in African Americans with more symptomatic tumors.
Because the ovarian steroids are critical to fibroid growth it is essential to assess their relationship with genomic and proteomic profiling presented here since the expression of some of these genes are the target of ovarian steroids regulatory functions. Since the tissues used in our study come from the early-mid secretary phase of the menstrual cycle it is possible to assume that both estrogen and progesterone-dependent and independent mechanism influence the expression of these genes. Additionally, we only identified the content of 15 of the protein spots isolated from the proteomic protein profiling as compared with a large number of genes identified by microarry. This limited our ability to assess any overlapping expression between genomic and proteomic results which determine whether different genes or proteins are regulated at different molecular level since some genes are regulated at transcription level, while others at protein level.
In conclusion, subjecting paired leiomyoma and myometrium from African Americans and Caucasians to genomic and proteomic analysis we identified a considerable similarity between their molecular environments with differences seen at the level rather than ethnic-specific expression of a number of genes. The area of relatedness among some of the gene clusters between myometrial and leiomyomas within each ethnic group suggests that some of the differences in leiomyoma's gene expression might be attributed to differences in myometrial gene expression between the ethnic groups, in addition to differences in leiomyomas vs. myometrium. Because many of the differentially expressed genes identified in these cohorts are know to regulate inflammatory response, angiogenesis, cell cycle, apoptosis and ECM matrix turnover, their products may account for rapid growth and associated symptoms in African Americans as compared to other ethnic groups.
Flake GP, Andersen J, Dixon D: Etiology and pathogenesis of uterine leiomyomas: a review. Environ Health Perspect. 2003, 111: 1037-1054.
Wise LA, Palmer JR, Harlow BL, Spiegelman D, Stewart EA, Adams-Campbell LL, Rosenberg L: Reproductive factors, hormonal contraception, and risk of uterine leiomyomata in African-American women: a prospective study. Am J Epidemiol. 2004, 159: 113-123. 10.1093/aje/kwh016.
Wise LA, Palmer JR, Stewart EA, Rosenberg L: Age-specific incidence rates for self-reported uterine leiomyomata in the Black Women's Health Study. Obstet Gynecol. 2005, 105: 563-568.
N C, R.Stan.Williams: Implication of growth factor and cytokine networks in leiomyomas. Cytokines in human reproduction. Edited by: Hill J. 2000, New York:, Wiley & Sons, 133-162.
N C: Gene expression and hormonal response. Uterine Leiomyomata:Pathogenesis and Management. Edited by: Ivo B. 2005, London, Taylor &Francis Group, 41-67.
Chegini N, Verala J, Luo X, Xu J, Williams RS: Gene expression profile of leiomyoma and myometrium and the effect of gonadotropin releasing hormone analogue therapy. J Soc Gynecol Investig. 2003, 10: 161-171. 10.1016/S1071-5576(03)00004-2.
Luo X, Ding L, Xu J, Williams RS, Chegini N: Leiomyoma and myometrial gene expression profiles and their responses to gonadotropin-releasing hormone analog therapy. Endocrinology. 2005, 146: 1074-1096. 10.1210/en.2004-1384.
Quade BJ, Wang TY, Sornberger K, Dal Cin P, Mutter GL, Morton CC: Molecular pathogenesis of uterine smooth muscle tumors from transcriptional profiling. Genes Chromosomes Cancer. 2004, 40: 97-108. 10.1002/gcc.20018.
Tsibris JC, Segars J, Coppola D, Mane S, Wilbanks GD, O'Brien WF, Spellacy WN: Insights from gene arrays on the development and growth regulation of uterine leiomyomata. Fertil Steril. 2002, 78: 114-121. 10.1016/S0015-0282(02)03191-6.
Wang H, Mahadevappa M, Yamamoto K, Wen Y, Chen B, Warrington JA, Polan ML: Distinctive proliferative phase differences in gene expression in human myometrium and leiomyomata. Fertil Steril. 2003, 80: 266-276. 10.1016/S0015-0282(03)00730-1.
Zaitseva M, Vollenhoven BJ, Rogers PA: In vitro culture significantly alters gene expression profiles and reduces differences between myometrial and fibroid smooth muscle cells. Mol Hum Reprod. 2006, 12: 187-207. 10.1093/molehr/gal018.
Ahn WS, Kim KW, Bae SM, Yoon JH, Lee JM, Namkoong SE, Kim JH, Kim CK, Lee YJ, Kim YW: Targeted cellular process profiling approach for uterine leiomyoma using cDNA microarray, proteomics and gene ontology analysis. Int J Exp Pathol. 2003, 84: 267-279. 10.1111/j.0959-9673.2003.00362.x.
Wei JJ, Chiriboga L, Arslan AA, Melamed J, Yee H, Mittal K: Ethnic differences in expression of the dysregulated proteins in uterine leiomyomata. Hum Reprod. 2006, 21: 57-67. 10.1093/humrep/dei309.
Al Hendy A, Salama SA: Catechol-O-methyltransferase polymorphism is associated with increased uterine leiomyoma risk in different ethnic groups. J Soc Gynecol Investig. 2006, 13: 136-144. 10.1016/j.jsgi.2005.10.007.
Denschlag D, Bentz EK, Hefler L, Pietrowski D, Zeillinger R, Tempfer C, Tong D: Genotype distribution of estrogen receptor-alpha, catechol-O-methyltransferase, and cytochrome P450 17 gene polymorphisms in Caucasian women with uterine leiomyomas. Fertil Steril. 2006, 85: 462-467. 10.1016/j.fertnstert.2005.07.1308.
Ligon AH, Morton CC: Leiomyomata: heritability and cytogenetic studies. Hum Reprod Update. 2001, 7: 8-14. 10.1093/humupd/7.1.8.
Villanova FE, Andrade PM, Otsuka AY, Gomes MT, Leal ES, Castro RA, Girao MJ, Nishimura E, Baracat EC, Silva ID: Estrogen receptor alpha polymorphism and susceptibility to uterine leiomyoma. Steroids. 2006, 71: 960-965. 10.1016/j.steroids.2006.07.005.
Luo X, Ding L, Xu J, Chegini N: Gene expression profiling of leiomyoma and myometrial smooth muscle cells in response to transforming growth factor-beta. Endocrinology. 2005, 146: 1097-1118. 10.1210/en.2004-1377.
Luo Xiaoping, Pan Qun, Chegini Nasser: Genomic and proteomic profiling II: Comparative assessment of gene expression profiles in leiomyomas, keloids, and surgically-induced scars. Reprod Biol & Endocrinol. 2007, BioMed Central, 5 (1): 35-10.1186/1477-7827-5-35.
Chegini N, Luo X, Ding L, Ripley D: The expression of Smads and transforming growth factor beta receptors in leiomyoma and myometrium and the effect of gonadotropin releasing hormone analogue therapy. Mol Cell Endocrinol. 2003, 209: 9-16. 10.1016/j.mce.2003.08.007.
Kanamori T, Takakura K, Mandai M, Kariya M, Fukuhara K, Kusakari T, Momma C, Shime H, Yagi H, Konishi M, Suzuki A, Matsumura N, Nanbu K, Fujita J, Fujii S: PEP-19 overexpression in human uterine leiomyoma. Mol Hum Reprod. 2003, 9: 709-717. 10.1093/molehr/gag088.
Kanamori T, Takakura K, Mandai M, Kariya M, Fukuhara K, Sakaguchi M, Huh NH, Saito K, Sakurai T, Fujita J, Fujii S: Increased expression of calcium-binding protein S100 in human uterine smooth muscle tumours. Mol Hum Reprod. 2004, 10: 735-742. 10.1093/molehr/gah100.
Ding L, Xu J, Luo X, Chegini N: Gonadotropin releasing hormone and transforming growth factor beta activate mitogen-activated protein kinase/extracellularly regulated kinase and differentially regulate fibronectin, type I collagen, and plasminogen activator inhibitor-1 expression in leiomyoma and myometrial smooth muscle cells. J Clin Endocrinol Metab. 2004, 89: 5549-5557. 10.1210/jc.2004-0161.
Levens E, Luo X, Ding L, Williams RS, Chegini N: Fibromodulin is expressed in leiomyoma and myometrium and regulated by gonadotropin-releasing hormone analogue therapy and TGF-beta through Smad and MAPK-mediated signalling. Mol Hum Reprod. 2005, 11: 489-494. 10.1093/molehr/gah187.
Luo X, Ding L, Chegini N: CCNs, fibulin-1C and S100A4 expression in leiomyoma and myometrium: inverse association with TGF-beta and regulation by TGF-beta in leiomyoma and myometrial smooth muscle cells. Mol Hum Reprod. 2006, 12: 245-256. 10.1093/molehr/gal015.
Xu J, Luo X, Chegini N: Differential expression, regulation, and induction of Smads, transforming growth factor-beta signal transduction pathway in leiomyoma, and myometrial smooth muscle cells and alteration by gonadotropin-releasing hormone analog. J Clin Endocrinol Metab. 2003, 88: 1350-1361. 10.1210/jc.2002-021325.
Hyytiainen M, Penttinen C, Keski-Oja J: Latent TGF-beta binding proteins: extracellular matrix association and roles in TGF-beta activation. Crit Rev Clin Lab Sci. 2004, 41: 233-264. 10.1080/10408360490460933.
Wu X, Blanck A, Norstedt G, Sahlin L, Flores-Morales A: Identification of genes with higher expression in human uterine leiomyomas than in the corresponding myometrium. Mol Hum Reprod. 2002, 8: 246-254. 10.1093/molehr/8.3.246.
Zhao Y, Wen Y, Polan ML, Qiao J, Chen BH: Increased expression of latent TGF-beta binding protein-1 and fibrillin-1 in human uterine leiomyomata. Mol Hum Reprod. 2007, 13: 343-349. 10.1093/molehr/gam007.
Crawford SE, Stellmach V, Murphy-Ullrich JE, Ribeiro SM, Lawler J, Hynes RO, Boivin GP, Bouck N: Thrombospondin-1 is a major activator of TGF-beta1 in vivo. Cell. 1998, 93: 1159-1170. 10.1016/S0092-8674(00)81460-9.
Ghahary A, Tredget EE, Shen Q: Insulin-like growth factor-II/mannose 6 phosphate receptors facilitate the matrix effects of latent transforming growth factor-beta1 released from genetically modified keratinocytes in a fibroblast/keratinocyte co-culture system. J Cell Physiol. 1999, 180: 61-70. 10.1002/(SICI)1097-4652(199907)180:1<61::AID-JCP7>3.0.CO;2-R.
Godar S, Horejsi V, Weidle UH, Binder BR, Hansmann C, Stockinger H: M6P/IGFII-receptor complexes urokinase receptor and plasminogen for activation of transforming growth factor-beta1. Eur J Immunol. 1999, 29: 1004-1013. 10.1002/(SICI)1521-4141(199903)29:03<1004::AID-IMMU1004>3.0.CO;2-Q.
Chegini N, Ma C, Tang XM, Williams RS: Effects of GnRH analogues, 'add-back' steroid therapy, antiestrogen and antiprogestins on leiomyoma and myometrial smooth muscle cell growth and transforming growth factor-beta expression. Mol Hum Reprod. 2002, 8: 1071-1078. 10.1093/molehr/8.12.1071.
Dillner NB, Sanders MM: Transcriptional activation by the zinc-finger homeodomain protein delta EF1 in estrogen signaling cascades. DNA Cell Biol. 2004, 23: 25-34. 10.1089/104454904322745907.
Postigo AA, Depp JL, Taylor JJ, Kroll KL: Regulation of Smad signaling through a differential recruitment of coactivators and corepressors by ZEB proteins. EMBO J. 2003, 22: 2453-2462. 10.1093/emboj/cdg226.
Du B, Fu C, Kent KC, Bush H, Schulick AH, Kreiger K, Collins T, McCaffrey TA: Elevated Egr-1 in human atherosclerotic cells transcriptionally represses the transforming growth factor-beta type II receptor. J Biol Chem. 2000, 275: 39039-39047. 10.1074/jbc.M005159200.
Ghosh R, Nadiminty N, Fitzpatrick JE, Alworth WL, Slaga TJ, Kumar AP: Eugenol causes melanoma growth suppression through inhibition of E2F1 transcriptional activity. J Biol Chem. 2005, 280: 5812-5819. 10.1074/jbc.M411429200.
Khachigian LM: Early growth response-1: blocking angiogenesis by shooting the messenger. Cell Cycle. 2004, 3: 10-11.
Lee CG, Cho SJ, Kang MJ, Chapoval SP, Lee PJ, Noble PW, Yehualaeshet T, Lu B, Flavell RA, Milbrandt J, Homer RJ, Elias JA: Early growth response gene 1-mediated apoptosis is essential for transforming growth factor beta1-induced pulmonary fibrosis. J Exp Med. 2004, 200: 377-389. 10.1084/jem.20040104.
Nishimura G, Manabe I, Tsushima K, Fujiu K, Oishi Y, Imai Y, Maemura K, Miyagishi M, Higashi Y, Kondoh H, Nagai R: DeltaEF1 mediates TGF-beta signaling in vascular smooth muscle cell differentiation. Dev Cell. 2006, 11: 93-104. 10.1016/j.devcel.2006.05.011.
Thiel G, Cibelli G: Regulation of life and death by the zinc finger transcription factor Egr-1. J Cell Physiol. 2002, 193: 287-292. 10.1002/jcp.10178.
Albig AR, Neil JR, Schiemann WP: Fibulins 3 and 5 antagonize tumor angiogenesis in vivo. Cancer Res. 2006, 66: 2621-2629. 10.1158/0008-5472.CAN-04-4096.
Bodner-Adler B, Nather A, Bodner K, Czerwenka K, Kimberger O, Leodolter S, Mayerhofer K: Expression of thrombospondin 1 (TSP 1) in patients with uterine smooth muscle tumors: an immunohistochemical study. Gynecol Oncol. 2006, 103: 186-189. 10.1016/j.ygyno.2006.02.023.
Brekken RA, Sullivan MM, Workman G, Bradshaw AD, Carbon J, Siadak A, Murri C, Framson PE, Sage EH: Expression and characterization of murine hevin (SC1), a member of the SPARC family of matricellular proteins. J Histochem Cytochem. 2004, 52: 735-748. 10.1369/jhc.3A6245.2004.
Kuznetsova SA, Issa P, Perruccio EM, Zeng B, Sipes JM, Ward Y, Seyfried NT, Fielder HL, Day AJ, Wight TN, Roberts DD: Versican-thrombospondin-1 binding in vitro and colocalization in microfibrils induced by inflammation on vascular smooth muscle cells. J Cell Sci. 2006, 119: 4499-4509. 10.1242/jcs.03171.
Nishiu J, Tanaka T, Nakamura Y: Identification of a novel gene (ECM2) encoding a putative extracellular matrix protein expressed predominantly in adipose and female-specific tissues and its chromosomal localization to 9q22.3. Genomics. 1998, 52: 378-381. 10.1006/geno.1998.5455.
Siegel PM, Shu W, Massague J: Mad upregulation and Id2 repression accompany transforming growth factor (TGF)-beta-mediated epithelial cell growth suppression. J Biol Chem. 2003, 278: 35444-35450. 10.1074/jbc.M301413200.
Fu H, Subramanian RR, Masters SC: 14-3-3 proteins: structure, function, and regulation. Annu Rev Pharmacol Toxicol. 2000, 40: 617-647. 10.1146/annurev.pharmtox.40.1.617.
Sarrazin S, Adam E, Lyon M, Depontieu F, Motte V, Landolfi C, Lortat-Jacob H, Bechard D, Lassalle P, Delehedde M: Endocan or endothelial cell specific molecule-1 (ESM-1): a potential novel endothelial cell marker and a new target for cancer therapy. Biochim Biophys Acta. 2006, 1765: 25-37.
Zhai J, Lin H, Shamim M, Schlaepfer WW, Canete-Soler R: Identification of a novel interaction of 14-3-3 with p190RhoGEF. J Biol Chem. 2001, 276: 41318-41324. 10.1074/jbc.M107709200.
Tzivion G, Gupta VS, Kaplun L, Balan V: 14-3-3 proteins as potential oncogenes. Semin Cancer Biol. 2006, 16: 203-213. 10.1016/j.semcancer.2006.03.004.
Ghaffari A, Li Y, Karami A, Ghaffari M, Tredget EE, Ghahary A: Fibroblast extracellular matrix gene expression in response to keratinocyte-releasable stratifin. J Cell Biochem. 2006, 98: 383-393. 10.1002/jcb.20782.
Tasheva ES, Conrad GW: Interferon-gamma regulation of the human mimecan promoter. Mol Vis. 2003, 9: 277-287.
Tasheva ES, Klocke B, Conrad GW: Analysis of transcriptional regulation of the small leucine rich proteoglycans. Mol Vis. 2004, 10: 758-772.
Hu SM, Li F, Yu HM, Li RY, Ma QY, Ye TJ, Lu ZY, Chen JL, Song HD: The mimecan gene expressed in human pituitary and regulated by pituitary transcription factor-1 as a marker for diagnosing pituitary tumors. J Clin Endocrinol Metab. 2005, 90: 6657-6664. 10.1210/jc.2005-0322.
We thank Drs. Mick Popp and Stanley Steven at Interdisciplinary Center for Biotechnology Research at the University of Florida for their assistance with microarray and proteomic analysis. The work presented here is supported by a grant HD37432 from the National Institute of Health. The work was presented in part at 53rd Annual Meeting of the Society for Gynecological Investigation, Reno NA, March 2007.
The author(s) declare that they have no competing interests.
QP, XL and NC participated in all aspect of the experiments presented here with the final microarray gene chips and two-d gels were performed at Interdisciplinary Center for Biotechnology Research at the University of Florida. Gene expression analysis performed by XL, proteomic analysis and immunohistochemistry by NC, realtime PCR and Western blot by XL and QP. All the authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table 4 - Differentially expressed genes in leiomyoma of African Americans and Caucasians. The complete list of differentially expressed genes identified in leiomyomas from African Americans and Caucasians as illustrated in Figure 1 and reported in part in Table 1. The genes were selected based on p ≤ 0.001 and 2-fold cutoff change (F. Change) as described in materials and methods. (DOC 51 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Pan, Q., Luo, X. & Chegini, N. Genomic and proteomic profiling I: Leiomyomas in African Americans and Caucasians. Reprod Biol Endocrinol 5, 34 (2007). https://doi.org/10.1186/1477-7827-5-34
- Protein Spot
- Ovarian Steroid
- Smad Interact Protein
- Proteomic Strategy