Genetic Determinants of Pelvic Organ Prolapse among African American and Hispanic Women in the Women’s Health Initiative

Current evidence suggests a multifactorial etiology to pelvic organ prolapse (POP), including genetic predisposition. We conducted a genome-wide association study of POP in African American (AA) and Hispanic (HP) women from the Women’s Health Initiative Hormone Therapy study. Cases were defined as any POP (grades 1–3) or moderate/severe POP (grades 2–3), while controls had grade 0 POP. We performed race-specific multiple logistic regression analyses between SNPs imputed to 1000 genomes in relation to POP (grade 0 vs 1–3; grade 0 vs 2–3) adjusting for age at diagnosis, body mass index, parity, and genetic ancestry. There were 1274 controls and 1427 cases of any POP and 317 cases of moderate/severe POP. Although none of the analyses reached genome-wide significance (p<5x10-8), we noted variants in several loci that met p<10−6. In race-specific analysis of grade 0 vs 2–3, intronic SNPs in the CPE gene (rs28573326, OR:2.14; 95% CI 1.62–2.83; p = 1.0x10-7) were associated with POP in AAs, and SNPs in the gene AL132709.5 (rs1950626, OR:2.96; 95% CI 1.96–4.48, p = 2.6x10-7) were associated with POP in HPs. Inverse variance fixed-effect meta-analysis of the race-specific results showed suggestive signals for SNPs in the DPP6 gene (rs11243354, OR:1.36; p = 4.2x10-7) in the grade 0 vs 1–3 analyses and for SNPs around PGBD5 (rs740494, OR:2.17; p = 8.6x10-7) and SHC3 (rs2209875, OR:0.60; p = 9.3x10-7) in the grade 0 vs 2–3 analyses. While we did not identify genome-wide significant findings, we document several SNPs reaching suggestive statistical significance. Further interrogation of POP in larger minority samples is warranted.


Introduction
Pelvic organ prolapse (POP) is a highly prevalent condition associated with significant morbidity that affects up to 40% of postmenopausal women [1], with 20% of women opting for surgical management by age 80 [2]. Although several risk factors for POP have been identified, including age [3], race [4;5], parity [1;6] and obesity [1;7;8], the underlying pathophysiology of POP remains unknown.
Epidemiologic studies have reported family history as a significant risk factor [9][10][11] and familial forms of POP have been reported [12;13]. In addition, several candidate gene studies [14][15][16][17][18] and one genome wide association study (GWAS) have identified promising single nucleotide polymorphisms (SNPs) associated with POP [19]. A recently published genomewide linkage analysis provided evidence for two additional loci in relation to symptomatic POP [20]. However, a majority of the existing studies have focused on women of European ancestry and there are few validated loci known for POP [21].
Given the limitations in the current literature, the Women's Health Initiative SNP Health Associated Resources (SHARe) dataset provides a unique opportunity to explore the genetic susceptibility to POP in minority women, as this dataset includes extensive phenotypic and genotypic information on African American (AA) and Hispanic (HP) postmenopausal women. Investigating the genetic epidemiology of POP in these populations is particularly important, as epidemiologic data suggest that the risk for POP varies by race. Prior studies have reported that AA women have the lowest risk [1;5], while white and Hispanic women have the highest risks [1;5]. Thus, the objective of our study was to identify loci associated with POP in AA and HP women using the Women's Health Initiative-SHARe dataset.

Study Population
Data used in this study were obtained from AA and HP women who were enrolled in the Women's Health Initiative hormone therapy (HT) trial (registered in ClinicalTrials.gov; registration number: NCT00000611) and for whom genotyping data was available through Women's Health Initiative-SHARe. The Women's Health Initiative is a large prospective study which recruited 161,861 post-menopausal women between 50-79 years of age from 40 clinical centers throughout the US from 1993-1998. Eligible participants were enrolled to the observational study or one or more of the three clinical trials: HT, dietary modification and/or calcium and Vitamin D supplementation trial [22]. Briefly, post-menopausal women who were unlikely to move and had predicted survival for three or more years, who were not using hormone therapy or were willing to stop, and who were currently not participating in any other clinical trial were eligible to participate.
Information regarding demographics, clinical, behavioral characteristics, medical history, and lifestyle/behavioral factors, among other risk factors, was obtained by standardized selfadministered questionnaires at baseline. Information regarding POP was ascertained through standardized pelvic exams performed on women participating in the Women's Health Initiative-HT trial. Baseline pelvic exams were performed using standardized procedures by a gynecologist, experienced nurse, or physician assistant. Pelvic exams included an evaluation of the presence of uterine prolapse, cystocele, and rectocele using Women's Health Initiative Prolapse Classification System (grades: 0, no prolapse; 1, prolapse in vagina; 2, prolapse to introitus; and 3, prolapse outside vagina). Examination for POP was performed in a supine lithotomy position while participant was asked to perform the Valsalva maneuver. Annual pelvic exams were also performed during follow-up physical exams with participants having between one to 10 visits. Further details regarding Women's Health Initiative protocols and ascertainment for Women's Health Initiative have been previously described [23;24].
For this study, cases were defined as women with a cystocele, rectocoele, or uterine prolapse of any severity (grades 1-3) at baseline or at one of the 10 follow-up visits. Since multiple assessments of POP were available for women in this study, the first occurrence of a cystocele, rectocele or uterine prolapse of any severity (grades 1-3) was coded as any POP, and the first occurrence of a cystocele, rectocele or uterine prolapse of grade II or higher was coded as moderate-severe POP. For women with multiple assessments of POP through time, the highest grade of POP across visits was used as their POP status for the purposes of our analyses. Controls were women who did not have any of the three forms of POP: rectocele, cystocele and uterine prolapse (grade 0) at baseline or at follow-up visits when POP was assessed. Our primary analyses compared controls to women with any POP (grade 0 vs. grade 1-3) or moderate/severe POP (grade 0 vs. grade 2-3). Variables examined in this study included age at diagnosis and variables collected at baseline: age at baseline interview, parity, body mass index (BMI), waist circumference, cigarette smoking history, prior hysterectomy, and menopausal hormone use status. De-identified data were assessed through authorized and secure methods after Institutional Review Board approval at Vanderbilt University and approval from the Women's Health Initiative, and the database of Genotypes and Phenotypes (dbGaP).

Genotyping
The Affymetrix Human SNP Array 6.0 (Affymetrix1, Inc Santa Clara, CA) was used for genome wide SNP genotyping. Genomic DNA was quantitated via an ND-8000 spectrophotometer and DNA quality was evaluated via gel electrophoresis. The genomic DNA samples were processed according to standard Affymetrix procedures for processing of the assay. The data were processed for genotype calling using the Affymetrix1 Genotypic Console software using the Birdseed calling algorithm version 2.0 (Affymetrix1, Inc., Santa Clara, CA) [25].

Quality Control (QC)
Data on 909,622 SNPs and 12,007 individuals were available prior to implementation of quality control. 177 individuals were removed after excluding individuals with low genotyping quality (<98%). 188 individuals with first degree or higher relatedness were identified using identityby-descent sharing from a random selection of 100,000 autosomal SNPs, and removed from analysis. Finally, we removed 140 subjects with inconsistent reported versus genetically determined sex. This resulted in 8,180 AA and 3,322 HP women remaining after sample QC.
All SNPs were tested for deviation from Hardy-Weinberg equilibrium (HWE) using PLINK software, stratified by race [26]. We excluded SNPs with HWE P10 −6 , low genotyping quality (<98%), a minor allele frequency (MAF) < 0.01, non-autosomal SNPs, and SNPs that did not map to a chromosomal position. Furthermore, we plotted the allele frequencies for the AA and HP datasets against the allele frequencies corresponding to the African reference population and American Indian reference population from the 1000 genomes build 37, respectively. We created these plots to identify and exclude SNPs which had discordant allele frequencies (absolute value difference > 0.1) compared with their respective reference populations, and to exclude palindromic SNPs with MAF > 0.4 for which the reference strand could not be determined easily. This resulted in a total of 777,060 SNPs in AAs and 730,985 SNPs in HPs available for imputation. Quality control procedures are detailed in Fig 1. Imputation was conducted using IMPUTE2, using samples and genotypes that passed QC and all 1000 genomes reference panels (build 37, 2013) to increase imputation accuracy [27;28]. Imputation was limited to 1769 AAs and 932 HPs, as data on POP was limited to these individuals. Prior to imputation, genotype data were phased using SHAPEIT software [29].

Analytic Sample
Of the 8180 AA and 3322 HP women with genetic data, 1,769 AA and 932 HP women had 1 pelvic exams and thus contributed information on POP. Of the 1,769 AA and 932 HP women, 514 AA and 437 HP women had POP (considering any severity of rectocele, cystocele and uterine prolapse) at baseline. During follow-up, 292 AA women and 184 HP women additionally developed POP. Our final analytical sample included 793 grades 1-3 cases, 154 grades 2-3 cases, and 948 controls in AA women, and 606 grades 1-3 cases, 163 grades 2-3 cases, and 305 controls in HP women who had no missing data on important covariates such as age at diagnosis, body mass index, hysterectomy status and parity.

Population Stratification
In order to assess population stratification among AAs and HP samples multi-dimensional scaling (MDS) was employed using PLINK software to estimate continuous axes of ancestry [26]. The top four MDS components were extracted for AAs and HP groups individually and used as covariates in SNP association analyses to adjust for potential confounding due to ancestry. A plot of the top two MDS components from the sample data (AAs and HP), and the International HapMap Project Phase 3 populations (ASW: African ancestry in Southwest USA;

Power
We performed power calculations a priori using the Quanto (Version 1.2.4) software. Assuming a meta-analysis sample size of approximately 1400 cases, a genome-wide p-value threshold of 5 x 10 −8 , an approximate 1:1 case to control ratio, a baseline prevalence of 40% and a logadditive model for SNP effect, we estimated 80% power to detect minimum odds ratios of 1. 48

Statistical Analyses
Participant and demographic characteristics were analyzed with two-sample t-tests assuming unequal variance to compare between racial groups when variables were continuous, and Chisquared tests when variables were categorical using Stata, version 11 (StataCorp, College Station, TX, USA).
The association between each genetic marker and POP status (binary) was assessed using logistic regression stratified by race adjusting for established POP risk factors identified in the literature including: age at POP ascertainment (continuous), BMI (continuous), parity (continuous) and MDS-derived axes of ancestry (continuous) using SNPTEST software [30]. Interpretation of analyses from logistic regression models were limited to genotyped and imputed SNPs with a post-imputation information score of 0.4, Hardy-Weinberg equilibrium P value >1x10 -8 , and minor allele frequency >0.05. Quantile-quantile plots for SNP association analyses for AA and HP subjects considering grade -0 vs. 1-3 and grade 0 vs. 2-3 POP are shown in weighted fixed-effect meta-analysis between AA and HP results using the PLINK software. We resorted to using a fixed-effect meta-analysis technique instead of more sophisticated transethnic meta-analysis techniques such as new random-effect models [31] or MANTRA [32] since only two independent datasets were available.
We performed several sensitivity analyses. We compared the results of six genome-wide significant SNPs reported by a previous GWAS study on POP with our results [19]. In the WHI-HT a large proportion of women had hysterectomy at baseline, with no indication for the reason for hysterectomy and the possibility that some of these women may have had POP previously could not be ruled out. Consequently, in the WHI-HT hysterectomy status was inversely associated with POP. However, women without a uterus may still develop other forms of POP such as rectocele and cystocele. In order to minimize the possibility that the associations observed with SNPs were not directed by any potential bias introduced by this misclassification, top hits (p < 10 −5 ) from primary models were additionally adjusted for hysterectomy status. We also performed logistic regression analyses by strata of hysterectomy status for the top-SNPs from primary models. All p-values presented are two-sided.

Results
In this case-control study, the mean ages at diagnosis for AA and HP women were 62 and 60 years respectively ( Table 1). Characteristics of AA and HP participants from this sub-study are similar to those in the parent Women's Health Initiative-HT study. Compared with HP women, AA women were on average more likely to be older, have higher BMI (29.8 vs. 31.6), more likely to be current smokers, and more likely to have had prior hysterectomy at baseline (40.7% vs. 57.7%). Advancing age, increasing BMI and parity were all associated with increased odds of POP across all grades (S1 Table).Prior hysterectomy was inversely associated with POP (OR: 0.24; 95% CI: 0.20, 0.29). HP women were approximately twice as likely as AA women to have POP (S1 Table). The racial disparity in POP existed irrespective of hysterectomy status.
We first evaluated the relationships between SNPs and any severity of POP in AA and HP populations separately using multiple logistic regression models adjusted for age at ascertainment, BMI, parity and four MDS components of genetic ancestry (Table 2). While no SNPs reached genome-wide significance, we noted those that reached a P value threshold of <10 −6 . In analyses specific to AA women considering any POP, SNP rs7035589 (OR: 1.44; 95% CI: 1.25-1.65; P = 2.92x10 -7 ), located near the gene ATP-binding cassette, subfamily-A (ABC1), member 1 (ABCA1) had the smallest P value. Similarly, in analyses specific to HP women considering any POP, SNP rs144039930 (OR: 0.15; 95% CI: 0.07-0.31; P = 3.80x10 -7 ) located in the intronic region of gene ankyrin repeat and sterile alpha motif domain containing 4B (ANKS4B) had the smallest P value. Regional association plots corresponding to these loci are presented in Fig 3. In addition to evaluating women with any POP, we also analyzed AA and HP women with moderate-severe POP (grades 2-3), adjusting for the same variables described above (Table 3). Again, in race-ethic stratified analyses, no loci reached genome-wide significance, but we identified SNPs in seven loci in the AA population and SNPs in three loci in the HP population which were associated with grades 2-3 POP at a P value threshold of <1.00x10 -6 . In analyses specific to AA women, the most significant SNP associated with moderate/severe POP was SNP rs28573326 (OR: 2.14; 95% CI 1.62-2.83; P = 1.04x10 -7 ) which is located in the intronic region of the carboxypeptidase E (CPE) gene. Similarly, in analyses specific to HP women we found several SNPs in and around the gene AL132709.5 to be associated with grades 2-3 POP, with rs1950626 (OR: 2.96; 95% CI: 1.96-4.48; P = 2.64x10 -7 ) as the most significant SNP in the region. Regional association plots for top loci from the grade 0 vs. 2-3 analyses are shown in Fig 4. for AA subgroup and Fig 5. for HP subgroup. Race/ethnicity specific correlation plots for top associated SNPs (P < 10 −5 ) from analyses of either any POP or moderate/severe POP not only showed that the effect size estimates were in the same direction, but that the magnitude of the effect size was larger in the moderate/severe POP analyses than in the any POP analyses (S6 Fig). We then aggregated the signals from the AA and HP analyses by performing inverse variance weighted fixed effects meta-analysis for both definitions of POP (Table 4): any POP   To evaluate the robustness of SNPs with the lowest P values in our study (P<10 −6 ), we performed several sensitivity analyses. To evaluate if there was more than one independent SNP contributing to the signals in a given region, we regressed SNPs in the region while conditioning on the SNP with the smallest P value for that respective region. We did not identify any additional SNPs that had signals <10 −6 after conditional analyses. Furthermore, since women who had previously had hysterectomy at baseline assessment were less likely to have/develop POP, we reevaluated the top SNPs (P<10 −5 ) from our original models while additionally adjusting for hysterectomy status (S7 Fig). The effect estimates we observed from these models  Finally, we provide effect estimates for six previously reported genome-wide significant SNPs [19]. We found suggestive evidence for rs1810636 (OR: 1.29; 95% CI: 0.99-1.69; P = 0.06) in our AA analyses when analyzing grade 0 vs. grades 2-3 POP. The magnitude and direction of OR for other SNPs were not consistent across ethnicities and close to the null (Table 5).

Discussion
We observed several potential loci which may be associated with POP among AA and HP women in the Women's Health Initiative. While strong associations were observed, these  associations were not genome-wide significant. However, we note that several of the associations were consistent with regards to the direction and magnitude across populations. Additionally, the magnitudes of associations were generally larger for the grade 0 vs 2-3 analyses than for grade 0 vs. 1-3 analyses. This is consistent with the hypothesis that women with a strong genetic risk for POP are more likely to have advanced POP. The strongest observed associations (smallest p-values) for grade 0 vs 2&3 were in both AAs (CPE gene) and HPs (nearby AL132709.5, MIR370 and RTL1). Meta-analysis of grade 0 vs. 1-3 analysis identified a SNP in the gene DPP6 had the strongest association with POP. CPE is a secretory protein involved in the biosynthesis of peptide hormones and neurotransmitters, such as insulin, involved in energy balance, nutrient partition, and satiety [33]. Animal model studies have associated mutations in the CPE gene with obesity, diabetes, infertility, and hyperinsulinemia in mice [34;35]. Although, there is no known direct relationship between CPE and POP, high BMI is an important risk factor for POP. Furthermore, metabolic syndrome, which involves obesity, cardiac risk factors, insulin resistance, systemic inflammation, atherogenic dyslipidemia and endothelial dysfunction, has associated with multiple urologic conditions including POP [36]. The strongest associated variant in HPs was not in a gene, but was located near AL132709.5, MIR370, and RTL1. There is no established relationship between these genes and POP or POP-related phenotypes. Another locus surrounding SNP rs7187140, identified in the AA grade 0 vs. 2-3 analysis, is also noteworthy in relation to POP. Although the SNP is intergenic, it is within 60 kb of the smooth muscle myosin heavy chain 11 (MYH11) gene. In a small gene expression study in tissue from anterior vaginal walls, Bortolini and colleagues showed MYH11 gene expression was down regulated in pre-menopausal women with POP compared with premenopausal women without POP [37].
Meta-analyses identified strong associations on the gene DPP6. The direction of the association for the most statistically significant SNP (rs11273354) was consistent across AAs and HPs. The SNP was also common in both populations (minor allele frequency >0.05). DPP6 is a single transmembrane protein an important protein for voltage-gated ion channels in determining cellular excitability. To date, there have not been any reports of an association of DPP6 with POP or any POP-related phenotypes. It is of interest that base-pair variation in the SNP rs11273354 alters binding motif of transcription factors containing the forkhead box D1 (FOXD1) gene [38]. Evaluation of expression quantitative trait loci (eQTL) using the GTEx portal [39] showed several SNPs in the DPP6 gene were associated with subcutaneous adiposity, potentially suggesting indirect ties between obesity and POP which is not fully understood. DPP6 has been associated with pancreatic cancer, tardive dyskinesia and amyotrophic lateral sclerosis in prior studies [40][41][42][43][44]. The piggyBac transposable element derived 5 (PGBD5) gene, identified in the grade 0 vs. 2-3 meta-analysis may also be a potential gene of interest in relation to POP. Although there have been no direct ties between PGBD5 and POP in the literature, it is of interest that several eQTL have been shown to alter expression of this gene in the tibial nerve [39]. The tibial nerve originates from the same lumbar and sacral nerves (L5 through S3) which also innervate the pelvic floor [45]. Studies suggest that neuro-modulation through percutaneous or posterior tibial nerve stimulation, may be beneficial in alleviating symptoms associated with pelvic floor disorders including fecal and urinary incontinence [46][47][48].
We took several measures to ensure the validity of the observed associations. In the WHI-HT a large proportion of women had undergone hysterectomy prior to study enrollment, with no information on the indication for hysterectomy. Considering that POP is the 3 rd leading indication for hysterectomy in women, it is then possible that some women might have had POP previously but was not detected at baseline. At the same time, women without a uterus can still develop other forms of prolapse such as rectocele and cystocele, as was the case during follow-up. To address this concern we performed analyses adjusting for hysterectomy and by strata of hysterectomy for our top hits and showed that these estimates were highly correlated. Therefore, bias introduced by hysterectomy status is not a likely explanation for associations observed in this study. Additionally, we compared the magnitudes of associations for our top hits in the any POP and moderate-severe POP analyses to show that not only were these associations in the same direction, but that the magnitude of effect estimates were larger for analyses considering moderate-severe POP than any POP. Effect estimates in the any POP analyses were likely driven by moderate-severe POP cases. One additional advantage of using the WHI-HT dataset is that information on important risk factors for POP was collected uniformly across all patients. In our analyses we adjusted for age, BMI and parity as these have been shown to be important risk factors for POP. Ideally, adjustment for mode of delivery would have provided for a more robust analysis, however, the WHI-HT did not collect this information.
We were not able to replicate the associations observed in a previously published GWAS by Allen-Brady et al., in a EA population with 115 cases and 2,976 controls [19]. Although cases were well-characterized, the study utilized general-population controls for whom the phenotype status was not verified. The study identified associations in six chromosomal regions 4q21, 8q24, 9q22, 15q11, 20p13, and 21q22. Two of the strongest associated SNPs were located within a gene: rs2236479 (COL18A1) and rs1036819 (ZFAT). The odds ratios we observed for these loci were close to the null and corresponding P values were >0.05. Potential reasons we were unable to replicate these findings are that genetic variants for POP may differ across racial/ethnic populations or power, as the minor allele frequency for these previously associated variants varies significantly across geographic populations.
To our knowledge, this is the second GWAS that evaluates POP risk and the first GWAS in AA and HP women. Majority of other genetic investigations evaluating POP risk have been limited to relatively smaller candidate gene studies that primarily focused on women of European or Asian descent [21]. The lack of GWAS evaluating POP, especially in minority populations and the lack of replication of candidate gene signals across existing studies are likely due to several factors and barriers, as highlighted by Wu et al. [49]. First, POP is a difficult phenotype to characterize as it involves a specialized pelvic exam and is rarely evaluated in research studies as routine practice. Secondly, a well conducted genetic epidemiologic study also equally rests on proper characterization of controls to reduce outcome misclassification. Prevalence of POP increases with age, especially after menopause and is highly prevalent in this population. The timing of control classification is thus of crucial importance, with ideal preference for including older post-menopausal women with no POP over younger pre-menopausal women to potentially reduce outcome misclassification.
Using the WHI-HT data allowed us to overcome these two barriers to a certain degree to conduct the study at hand. The WHI-HT used centrally validated standardized procedures performed by trained medical staff to collect information on POP at baseline and at up to 10-annual visits per individual (all of whom were post-menopausal). This not only allowed us to properly identify POP cases who developed POP during follow-up visits, but also allowed us to classify women as controls considering multiple visits, when possible. Even though is a marked improvement over previous studies that have investigated POP, there is likely some degree of misclassification in our study as well. While information on POP was available for multiple visits, the loss-to-follow-up rate in the WHI-HT was extremely high. Consequently, for controls identified at baseline that were lost to follow-up, we could not determine POP status through time. For majority of the cases who had two or more assessments of POP (612 AA POP cases and 498 HP POP cases), POP could be verified at two or more visits for 70% of these AA POP cases and 80% of these HP POP cases. For some women POP grade fluctuated between visits; this speaks to a degree of subjectivity in POP measurement even when using standardized procedures. For women with information on POP for multiple visits, we noted the visit at the first ascertainment of POP when considering any POP and the first ascertainment of grade 2 or higher POP when considering moderate/severe POP. It should also be noted that POP was measured in the supine lithotomy position in the WHI-HT, and therefore, the maximal extant of prolapse was likely underestimated compared with a semi-upright or upright position of assessment, which is ideal for determining the maximum extent of prolapse.
Thirdly, resources for large scale genetic investigations for phenotypes such as POP are extremely limited, making replication of findings a difficult task to achieve. This especially rings true for investigation in minority populations such as AAs and HPs, who have historically been under-represented in medical research in general. We were unable to replicate our findings in an independent cohort of AA and HP women, and stress the need for an independent replication of our findings. To address this, we performed trans-ethnic replication across AAs and HPs. Although we have a strong sample size of 2,701 women, the number of women with POP was limited (any POP: n = 1,382). However, a small proportion of women (n = 317) had grades 2-3 POP and even fewer women had POP which protruded beyond the introitus. This resulted in our analysis being underpowered to detect smaller associations. Given the highly stringent p-value threshold due to multiple testing, the possibility that our findings are due to chance alone cannot be ruled out. Additionally, since we had only two datasets for meta-analysis, we resorted to using a fixed effects meta-analysis framework, instead of newer trans-ethnic meta-analysis techniques such as new-random-effects model [31] or MANTRA [32]. Both these methods have been shown to boost power especially in scenarios with allele frequency differences and/or extreme allelic heterogeneity [50]. Simulations of a considerably larger number of studies (N = 30) have demonstrated considerably improved power for these techniques, however, these techniques are also underpowered as the number of meta-analysis studies decrease (N = 10) in the presence of effect-size heterogeneity [50]. The added benefit from utilizing these techniques when the number of studies is small (N = 2) is not clear. By design our meta-analysis is most sensitive to identifying associations for variants with similar allele frequencies and effect estimates across populations. Considering lower prevalence of POP in AA women compared with HP women it is likely that there may be ancestry specific loci which we were not able to detect through this meta-analytic technique. Future meta-analysis of a greater number of studies with diverse populations, coupled with meta-analysis techniques designed to leverage heterogeneity associated with between population differences may provide insight into these loci.
To our knowledge, the WHI-HT data is the only resource available which has validated detailed information on POP and simultaneous availability of genetic data. The surge of interest in utilizing Electronic medical record (EMR) databases for research investigation provides future hope for replicating these results. However, constructing a well characterized set of cases and controls requires building algorithms that can identify cases and controls with a high degree of sensitivity and specificity, respectively. To this end, our group is currently in the process of developing algorithms using EMR data from the BioVU DNA repository at Vanderbilt University, which we hope to apply in larger EMR database networks with available GWAS data for replication in the future.
Little is known about POP pathophysiology or genetic risk factors, particularly among AA and HP women. Our study is the first POP GWAS in minority women. Our findings suggest that common germ-line variations may contribute to increased risk for POP among AAs and HPs; however, further research with larger minority sample sizes is necessary to validate our study findings. any POP models to grade 0 vs. moderate/severe POP models. This figure compares beta-estimates of top hits (p <10-5) originating from either any POP or moderate/severe POP models where all models were adjusted for age, BMI, parity and ancestry components. X-axis and Yaxis represents natural log transformed odds ratios from any POP and moderate/severe POP models, respectively. left plot: AAs; right plot: HPs. (TIF) S7 Fig. Summary of log(OR) of single SNP association analyses adjusted for hysterectomy compared to no adjustment for hysterectomy. This figure compares beta-estimates of top hits (p < 10 −5 ) from models adjusted for age, BMI, parity and ancestry components to models adjusted for hysterectomy in addition to the aforementioned variables. Original refers to betas from models adjusted for age, BMI, parity and ancestry components; Hysterectomy adjusted refers to betas from models adjusted for hysterectomy status in addition to factors in the original models. (TIF) S8 Fig. Summary of log(OR) estimates from single SNP association analyses by strata of hysterectomy status. This figure compares beta-estimates of top hits (p < 10 −5 ) from models adjusted for age, BMI, parity and ancestry components in women with and without a hysterectomy at baseline. X-axis and Y-axis represents natural log transformed odds ratios from analyses in women without hysterectomy and women with hysterectomy, respectively. Top 2 plots are for AAs (Any POP and moderate/severe POP, left to right). Bottom two plots are for HPs (Any POP and moderate/severe POP, left to right). (TIF) S1 Table. Associations between key risk factors for any pelvic organ prolapse and grades 2-3 pelvic organ prolapse in the WHI Hormone Therapy trial. The following table presents the associations between key risk factors in relation to any POP (grades 1-3) and severe/moderate POP (grades 2-3). (DOCX) S2 Table. Results from meta-analysis of Grade 0 vs. Grade 1-3 POP across African American (AA) and Hispanic (HP) Women from the Women's Health Initiative. The following table presents results on single nucleotide polymorphisms (SNPs) associated with Grade 1-3 prolapse, which had p < 0.00001 in the fixed effects meta-analysis. Results from randomeffects models and Heterogeneity Scores (I 2 ) are also presented. (DOCX) S3 Table. Results from meta-analysis of Grade 0 vs. Grade 2-3 POP across African American (AA) and Hispanic (HP) women from the Women's Health Initiative. The following table presents results on single nucleotide polymorphisms (SNPs) associated with Grade 2-3 prolapse, which had p < 0.00001 in the fixed effects meta-analysis. Results from randomeffects models and Heterogeneity Scores (I 2 ) are also presented. (DOCX)