[B]{.smallcaps}ackground {#Sec1}
========================

In sub-Saharan Africa, maternal disorders are the second most important cause of death among women of reproductive age (15-49 years old) \[[@CR1]\]. Disease burden attributable to maternal complications still remains important despite the significant declines in maternal mortality observed in this region since the mid-2000s \[[@CR2], [@CR3]\]. In fact, it is estimated that for every woman dying from maternal complications, another 20 women will have to withstand serious maternal morbidity \[[@CR4]\]. Of all maternal morbidities, obstetric fistula is one of the most debilitating conditions with the immediate consequence of chronic urinary and/or fecal incontinence. Physical comorbidities, psychological distress, and social stigmatization usually follow \[[@CR5]--[@CR9]\].

The etiology of vaginal fistula (VF), an abnormal hole between the bladder (vesico-vaginal fistula) and/or rectum (recto-vaginal fistula) and the reproductive tract of a woman, is divided into two main categories: obstetric and traumatic. VF of obstetric origin are caused by an intertwined set of biological, socio-economic, and cultural factors that favor obstructed labor and triggered by insufficient or delayed access to quality emergency obstetric care \[[@CR7]\]. VF of traumatic origin mostly results from sexual violence. The vast majority of VF in sub-Saharan Africa are of obstetric origins and prevalence of this condition in this region was recently estimated to be between 1.0 and 1.6 per 1,000 women of reproductive age depending on methodology \[[@CR10], [@CR11]\].

The epidemiological evidence regarding risk factors for VF is dominated by facility-based studies \[[@CR12]\]. The numerous clinical series usually report socio-demographic characteristics of VF patients (age of marriage, marital status, literacy, parity, etc.) as well as circumstances of fistula occurrence (duration of labor, type of birth attendance, mode and place of delivery, etc.) \[[@CR13]--[@CR21]\]. This accumulation of hospital-based studies contributed to highlight the diverse characteristics of fistula sufferers who present to facilities \[[@CR12]\]. A few case-control studies tackle individual determinants with the aim to confirm risk factors \[[@CR22]--[@CR24]\] or develop a fistula prevention index \[[@CR25]\]. Other studies, often qualitative, reflect on cultural or health system factors to reduce the three delays causative of obstetric fistula and maternal mortality \[[@CR21], [@CR26]--[@CR28]\]: delay in decision to seek care, delay in reaching care, and delay in receiving adequate care once in the health facility.

Population-based studies could be less susceptible to selection bias than case series from facility and case-control studies but are rarely carried-out \[[@CR11], [@CR29]\]. In sub-Saharan Africa, the main sources of nationally representative health data are Demographic and Health Surveys (DHS) and Multiple Indicator Cluster Surveys (MICS). Since 2004, these surveys progressively began to include questions about VF symptoms. These data sources were recently used to estimate the prevalence of VF in sub-Saharan Africa \[[@CR10]\], while adjusting for the imperfect accuracy of self-reports, but a thorough and systematic examination of individual risk factors has yet to be completed. We are aware of four population-based studies that examined determinants of VF \[[@CR30]--[@CR33]\]. These studies only included a small number of surveys, assumed that self-reports of VF symptoms were perfectly accurate, and none pooled surveys together, severely limiting their statistical power.

The primary objective of this paper is thus to examine the association between selected individual risk factors and lifetime prevalence of self-reported VF symptoms, such as literacy status, education level, female genital mutilation (FGM), sexual violence, short stature, age at first sexual intercourse, age at first birth, and women's ability to get permission to seek health care. By pooling surveys from different countries, we hope to improve the representativeness and precision of the effect size measures for those risk factors.

[M]{.smallcaps}ethods {#Sec2}
=====================

Data sources {#Sec3}
------------

DHS and MICS surveys conducted in sub-Saharan Africa that included questions about VF symptoms were considered for this analysis. A comprehensive overview of DHS and MICS surveys can be found elsewhere \[[@CR34]\]. Briefly, both DHS and MICS are household-based surveys that use a multistage stratified cluster sampling design to select a nationally representative sample of women of reproductive age (15-49 years old). Socio-demographic characteristics and information on selected health indicators are collected through face-to-face interviews by trained personnel and recorded in standard questionnaires. The majority of surveys administered the VF questions to all women of reproductive age but some restricted it to women that were ever married (Mauritania MICS 2011), ever pregnant (Swaziland MICS 2010 and Guinea-Bissau MICS 2010), or that had a live birth in the previous five years (Rwanda DHS 2005).

Procedures {#Sec4}
----------

The specific questions related to vaginal fistula symptoms varied slightly from survey to survey and a contingency question about knowledge of vaginal fistula was sometime incorporated. A full description of the VF and contingency questions (if any), their probes, and the coding of the outcome can be found elsewhere \[[@CR10]\].

Based on previous studies and the information available from DHS/MICS surveys, we estimated the effect of the following risk factors: illiteracy, education level, whether the respondent has experienced female genital mutilation (FGM), short stature, experience of intimate partner sexual violence, young age at first sexual intercourse, young age at first birth, and women's difficulty to get permission to access health care. Literacy status was ascertained in the surveys by asking the interviewee to read a sentence on a card that was handed out to her. If the woman was able to read only part of the sentence, she was considered not being able to read properly. Women who reported having had some secondary education or higher were de facto assumed to be literate. For genital mutilation, we did not stratify our analysis by FGM type as a validation study of the DHS FGM questions in Sierra Leone demonstrated that they were accurate to determine FGM prevalence but inaccurate for determining cutting extent \[[@CR35]\]. Not all surveys recorded information for these risk factors and the list of countries for which such data was collected is presented in Tables [1](#Tab1){ref-type="table"} and [2](#Tab2){ref-type="table"}. As for women's anthropometric measurements, this information is not collected by MICS surveys and the women's height was recorded from a sub-sample of participants in most DHS surveys. Similarly, questions on domestic violence were often administered to a subsample of women, depending on the survey, and the questions about ever having experienced intimate partner sexual violence were only asked to ever married women (or those in a union). As for age at first sexual intercourse, inconsistent responses were disregarded and considered as missing (e.g., a women reporting never having had sexual intercourse but having given birth). Finally, most DHS surveys asked women if getting permission to seek health care was a problem. Those who responded that it was a big problem were considered as having limited ability to seek the care they need.Table 1Number of vaginal fistula (VF) by survey datasets for the following risk factors: literacy status, education level, female genital mutilation (FGM), and short stature (\<150 cm)Country, Survey, and YearLiteratePost-primary educationFGMShort StatureVF (N)/Literate (N)VF (N)/Illiterate (N)VF (N)/Post-primary (N)VF (N)/No Post-primary (N)VF (N)/FGM (N)VF (N)/No FGM (N)VF (N)/Short (N)VF (N)/Tall (N)Benin DHS 2011-1236/418291/1241731/345096/1314910/1370117/1522912/1115111/15061Burkina Faso DHS 20106/306314/139826/220114/1485417/128843/41760/24310/8202Cameroon DHS 201132/921224/616327/714829/8271Not Measured1/41125/7473Chad MICS 20106/233434/134375/185735/1400024/721816/8569Not MeasuredComoros DHS 201250/297459/231949/276160/2547Not Measured16/59092/4593Congo (Brazz.) DHS 2011-1217/674010/405214/629913/4519Not Measured2/51312/5134DRC DHS 200722/501522/495419/410725/5879Not Measured3/78212/3949Ethiopia DHS 200521/393782/1010310/265093/1140569/1001232/36031/85244/5798Guinea DHS 201210/181252/730610/170353/743263/89350/1944/25627/4457Kenya DHS 2008-0954/580924/260921/279657/564019/253959/58915/56972/7742Malawi DHS 201059/1347173/950014/4289118/18721Not Measured7/113628/6496Mali DHS 20060/180918/127300/154618/1303418/117400/282717/1380817/13808Mali DHS 2012-1312/181059/86149/168962/873569/94802/9440/13828/5134Mauritania DHS4/34593/55372/15395/75574/67023/2377Not MeasuredNiger DHS 20062/118518/79851/96719/82222/20918/89520/17912/4351Niger DHS 20120/167216/94320/137316/97591/20915/109240/2147/4912Nigeria DHS 200855/1434586/1880451/1352791/1979037/8452102/2470920/3853120/28572Rwanda DHS 200580/304184/233212/519152/4867Not Measured6/27366/2422Senegal DHS 2010-113/398315/117051/280217/1288611/56897/99990/1268/5633Sierra Leone DHS 201323/541589/1108522/520690/11337105/147736/17604/76446/7185Swaziland MICS 201048/29059/39536/219521/1123Not MeasuredNot MeasuredTanzania DHS 201032/690318/32057/236044/77769/132242/880713/163637/8408Togo MICS 20104/225519/41084/174419/46310/39323/5970Not MeasuredTogo DHS 2013-201435/357958/587427/307066/64046/60287/88612/31450/4517Uganda DHS 200680/3867121/460627/1823174/66532/61199/840311/24863/2596Uganda DHS 201176/429888/430732/2509132/60977/156156/84238/19451/2493Zambia DHS 2013-201441/955449/677434/738657/9001Not Measured12/167779/14562VF = Vaginal Fistula; FGM = Female genital mutilation; DHS = Demographic and Health Survey; MICS = Multiple Indicators Cluster SurveyThe survey-specific total sample sizes can vary by risk factor depending on the number of missing observations and eligibility criteriaTable 2Number of vaginal fistula (VF) by survey datasets for the following risk factors: experience of intimate partner sexual violence (IPSV), young age at first sex (\<14 years old), young age at first birth (\<14 years old), and permission to seek health careIPSV^a^Young age at 1^st^ intercourse^b^Young age at 1^st^ birth^c^Permission to seek health careCountry, Survey, and YearVF (N)/Sex Violence (N)VF (N)/No Violence (N)VF (N)/Young (N)VF (N)/Old (N)VF (N)/Young (N)VF (N)/Old (N)VF (N)/Big problem (N)VF (N)/Not big problem (N)Benin DHS 2011-12Not Measured17/125083/116815/60484/1191855/579772/10802Burkina Faso DHS 20100/14211/98590/52318/139710/10616/131297/385013/13208Cameroon DHS 20115/5749/34253/139649/111482/52647/104927/225122/5195Chad MICS 20101/84433/112835/216533/117295/90831/11933Not MeasuredComoros DHS 20122/3661/24929/24375/31282/11075/281834/240573/2885Congo (Brazz.) DHS 2011-12Not Measured5/151022/84050/25024/853715/475212/6056DRC DHS 20077/76410/20829/95832/72683/25032/689010/211534/7861Ethiopia DHS 2005Not Measured22/189272/81726/36287/896438/399965/10049Guinea DHS 2012Not Measured18/141541/63835/36852/6576Not MeasuredKenya DHS 2008-0918/62639/427310/59761/58902/18973/5906Not MeasuredMalawi DHS 20107/84219/453122/1904107/167725/536121/1749622/2571110/20421Mali DHS 20060/30713/86130/155417/106090/52517/110390/276418/11787Mali DHS 2012-131/37819/27427/110856/71513/48064/800030/306341/7361Mauritania MICS 2011Not MeasuredNot MeasuredNot MeasuredNot MeasuredNiger DHS 2006Not Measured2/129415/62251/25616/69291/86019/8315Niger DHS 2012Not Measured5/116010/83565/33811/88654/242912/8713Nigeria DHS 20085/68888/1850921/3221102/232555/1144114/2255215/4775127/28412Rwanda DHS 200514/25733/18220/88156/51340/17164/53699/102154/5278Senegal DHS 2010-11Not Measured4/150512/97171/60016/100527/257611/13112Sierra Leone DHS 20133/24839/402711/135887/122469/69794/1157030/268682/13831Swaziland MICS 2010Not Measured5/11852/31980/957/3270Not MeasuredTanzania DHS 20107/69525/49912/49840/76391/9642/72271/22950/9899Togo MICS 2010Not Measured1/29922/54172/12418/4592Not MeasuredTogo DHS 2013-1412/42052/49496/49076/74162/15178/678910/117683/8294Uganda DHS 200617/51625/122532/769158/62146/173189/623124/653177/7815Uganda DHS 201115/42321/127122/822127/59978/292137/607612/475152/8127Zambia DHS 2013-1418/150253/78964/89476/122053/19178/122204/49187/15882VF = Vaginal Fistula; IPSV = Intimate Partner Sexual Violence; DHS = Demographic and Health Survey; MICS = Multiple Indicators Cluster Survey^a^Among married and/or ever married women (or those in a union)^b^Among sexually active women^c^Among primi/multiparous womenThe survey-specific total sample sizes can vary by risk factor depending on the number of missing observations and eligibility criteria

The principal threat to the internal validity of our analyses is confounding of the exposure-outcome relationship. The main potential confounders for which information was collected by the survey questionnaires are age, literacy status, location of residence (rural versus urban), gravidity status, and religion. Socio-economic status and marital status were not considered in this analysis because these variables are likely both a cause and an effect of VF. That is, due to the cross-sectional nature of data collection, we do not have information on the temporal sequence in which changes in socio-economic status or marital status occurred. Three surveys (Chad MICS 2010, Mauritania MICS 2011, and Togo MICS 2010) did not record information on gravidity status and we assumed that all nulliparous women were also nulligravid -- a reasonable assumption giving the high correlation observed between these two variables. Finally, four surveys did not record information on religion and these were coded using a missing variable indicator to retain them in the analyses (Mauritania MICS 2011, Niger DHS 2012, Swaziland MICS 2010, and Tanzania DHS 2010).

Statistical analyses {#Sec5}
--------------------

To circumvent the lack of balance and overlap for some of the covariates, matching was used to make the group with the selected risk factor (i.e., exposed) as similar as possible to the group without (i.e., unexposed). By reducing model dependency through this semi-parametric data preprocessing, we aim to produce more robust inferences that are less sensitive to modeling assumptions \[[@CR36]\]. Three of our risk factors are continuous and were dichotomized. Respondents with a height less than 150 cm, a commonly used threshold \[[@CR12], [@CR15]\], were defined as having a short stature. For age at first birth, visual inspection of the exposure-outcome relationship suggested that this variable could be dichotomized at less than 14 years of age at first delivery. This corresponds roughly to the 4^th^ percentile of the distribution of age at first birth. The same threshold of less than 14 years was used to define young age at first sexual intercourse.

All country datasets were pooled together as the low number of VF cases precludes data analysis at the country level for many surveys (i.e., all cases were either exposed or unexposed in these surveys). For the selected risk factors, a nearest neighbor algorithm was used to match women on sampling weight (for sexual violence, the sampling weight from the domestic violence questionnaire was used), age (continuous), and survey identifier. For this latter variable, exact matching was used for risk factors that consistently had more unexposed than exposed observations across surveys: short stature, intimate partner sexual violence, young age at first sexual intercourse, young age at first birth, and problem obtaining permission to seek care (otherwise, nearest neighbor matching was used). The matching ratio of exposed to unexposed units varied for each risk factor and was chosen as to minimize unbalance and maximize statistical power. Matching was implemented using the '*MatchIt*' package \[[@CR37]\] in R. Unmatched women were excluded from the analyses.

Logistic regression models were used on the matched data to estimate the effect of the selected risk factors on lifetime prevalence of VF. Missing values for the selected risk factors and covariates were always less than 1 %, except for height (2.0 %) and age at first sexual intercourse (6.1 % of inconsistent or missing values). Observations with missing values were excluded from the analyses (with the exception of those for religion which were retained using a missing indicator). To provide for additional control of potential confounders, we adjusted for the following covariates: age (15-19, 20-29, and 30-49 years), literacy status (this covariate was not included for literacy status and education level risk factors), gravidity status (not included for age at first birth), location of residence (urban/rural), religion (Christian, Muslim, others, missing), and the survey's country. Such analyses have been described as doubly-robust because statistically consistent inferences can be made "*if either the matching analysis or the analysis model is correct (but not necessarily both)*" \[[@CR37]\]. Surveys that had a different population denominator were included in the analysis since we matched on survey identifier and country fixed effects were included in the parametric analyses. These logistic regressions did not account for the clustered design of surveys as our preliminary analyses have shown that clustering the standard errors had no impact on our conclusions (also discussed in \[[@CR10]\]).

Importantly, women's self-report of vaginal fistula symptoms do not have perfect sensitivity and specificity, as compared to the gold standard of a pelvic examination. In order to account for non-differential misclassification of the self-reported outcome, we used a latent-class Bayesian statistical model \[[@CR10], [@CR38], [@CR39]\]. The underlying assumption being that all surveys have a common sensitivity and specificity (see \[[@CR10]\] for details). This model takes the following form:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \begin{array}{c}\mathrm{Likelihood}:\\ {}{y}_i\sim Binomial\left({p}_i,{N}_i\right)\\ {}{p}_i={\pi}_i(Se)+\left(1-{\pi}_i\right)\left(1-Sp\right)\\ {}\mathrm{logit}\left({\pi}_i\right)=\alpha +\beta {X}_i\end{array} $$\end{document}$$

Because of our very large sample sizes and the computing-intensive nature of Bayesian calculations, we grouped observations with the same covariate patterns and used a binomial likelihood instead of the standard Bernoulli (i.e., grouped logistic regression). In this model, *y*~*i*~ is the total number of women reporting VF symptoms with covariate pattern *i*; *N*~*i*~ is the total number of women with covariate pattern *i*; *p*~*i*~ is the observed probability of reporting VF symptoms, *π*~*i*~ is the true probability of women having ever had VF symptoms; *Se* and *Sp* are the sensitivity and specificity of the survey instrument, respectively; *α* is the model's intercept; *β* is a vector of coefficients for the covariates included in *X*~*i*~. The model's specification is completed using the following prior distributions:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \begin{array}{c}\mathrm{Prior}\ \mathrm{distributions}\ \mathrm{f}\mathrm{o}\mathrm{r}\ \mathrm{model}\ \mathrm{parameters}:\\ {}\alpha \sim Normal\left(0,20\right)\\ {}\beta \sim Normal\left(0,20\right)\\ {}Se\sim Uniform\left(95.10\%,99.90\%\right)\\ {}Sp\sim Uniform\left(99.85\%,99.95\%\right)\end{array} $$\end{document}$$

Both *α* and *β* are given non-informative priors that follow a normal distribution with a mean of zero and standard deviation of 20. For sensitivity and specificity, we used uniform distributions that match the 95 % credible intervals of the posterior distributions of these quantities, as estimated previously \[[@CR10]\]. Posterior distributions were obtained using Markov Chain Monte Carlo sampling, implemented in R using the 'rstan' package \[[@CR40]\]. Samples are obtained using the no-U-turn sampler, a computationally efficient variant of Hamiltonian Monte Carlo \[[@CR41]\]. Inferences were based on three chains of 30,000 samples after an initial warm-up period of 2,500 samples per chain (total of 90,000 iterations used for inferences). Convergence was examined using traceplots and ensuring that the potential scale reduction factor was equal to one. All analyses were performed using the R statistical software \[[@CR42]\].

[R]{.smallcaps}esults {#Sec6}
=====================

Surveys characteristics {#Sec7}
-----------------------

A total of 31 surveys collected information on VF symptoms in sub-Saharan Africa. Of these, individual data records were not available for two surveys (Equatorial Guinea DHS 2011 and Guinea-Bissau MICS 2010), and two other surveys were excluded because the incontinence questions were considered to be non-specific (Côte d'Ivoire DHS 2011-12 and Malawi DHS 2004). Hence, 27 surveys, conducted between 2005 and 2014, informed our analyses. The main characteristics of the interviewees can be found in Additional file [1](#MOESM1){ref-type="media"}: Table S1.

These 27 surveys pooled self-reports from 334,606 eligible women and 2,048 reported having ever experienced VF symptoms (742 had missing information on the outcome (0.2 %)). The specific sample size used in the regression models varied, depending on the considered risk factors, from 332,889 for literacy to 102,928 for intimate partner sexual violence (before matching). Detailed information on the risk factors and number of women reporting VF symptoms, stratified by surveys, can be found in Table [1](#Tab1){ref-type="table"} and Table [2](#Tab2){ref-type="table"}. Briefly, a little over a third of women were able to read (38.6 %), a quarter had completed post-primary education (26.8 %), 42.2 % had experienced FGM, 8.7 % had a height below 150 cm, 9.0 % of ever married women had experienced intimate partner sexual violence, 11.0 % of sexually active women had their first sexual intercourse before the age of 14, 3.8 % of primi/multiparous women had their first birth before the age of 14, and 18.3 % of women reported that obtaining permission to seek health care was a big problem for them.

Risk factors for vaginal fistula {#Sec8}
--------------------------------

The sample size of the pooled datasets before and after matching are presented for each risk factors in a supplementary appendix (Additional file [1](#MOESM1){ref-type="media"}: Table S2). One-to-one matching was used for the risk factors that were most prevalent: being able to read, having a post-primary education, female genital mutilation, and degree of difficulty in obtaining permission to seek health care. For the other risk factors, the ratio was chosen as to minimize imbalances while retaining sufficient statistical power: one-to-two matching for intimate partner sexual violence, one-to-three for young age at first sexual intercourse, one-to-four for short stature, and one-to-eight for young age at first birth.

Results from the matched logistic regressions are presented in Table [3](#Tab3){ref-type="table"}. Preliminary results from the Bayesian models for young age at 1^st^ birth and problem getting permission to seek healthcare suggest convergence issues with the country fixed effects. Since matched logistic regressions with and without country fixed effects for these two risk factors gave very similar results (data not shown), they were omitted from the Bayesian model.Table 3Matched logistic regression results for the different risk factors for vaginal fistula symptomsRisk FactorsBayesian matched logistic regressions adjusting for outcome misclassification**OR (95 % CrI)**Being able to read**0.87 (0.77-0.99)**Post-primary education0.90 (0.76-1.06)Female genital mutilation1.04 (0.82-1.30)Short stature (\<150 cm)**1.31 (1.02-1.68)**Intimate partner sexual violence^a^**2.13 (1.60-2.86)**Young age at 1^st^ intercourse (\<14 years)^b^**1.41 (1.16-1.71)**Young age at 1^st^ birth (\<14 years)^c^**1.39 (1.04-1.82)**Problem with permission to seek care1.20 (0.99-1.47)OR = Odds ratio; 95 % CI = 95 % Confidence Interval; 95 % CrI = 95 % Credible Intervals.Statistically significant results at the alpha = 0.05 level are bolded.The matched logistic regression models adjust for the following covariates: age, literacy status (except for '*Being able to read*' and '*Post-primary education*'), location of residence (urban/rural), gravidity status (except for '*Young age at 1* ^*st*^ *birth*'), religion, and country (country fixed effects were omitted from the Bayesian regressions for '*Young age at 1* ^*st*^ *birth'* and '*Problem getting permission to seek healthcare'*).^a^Among married and/or ever married women (or those in a union).^b^Among sexually active women.^c^Among primi/multiparous women.

Being able to read decreased the odds of VF by 13 % (95 % Credible Intervals (CrI): 1 % to 23 %). The impact of having completed some post-primary education also reduced the odds of VF by 10 % (95 % CrI: -6 % to 24 %) but the effect did not reach statistical significance. For these two determinants, it is likely that gravidity status lies on the causal pathway between literacy/education and occurrence of VF. If that is the case, the effect size measures reported above should be interpreted as the direct effect of literacy/education on VF (i.e., the effect not mediated through gravidity). By not controlling for gravidity status, we can calculate the total effect of literacy/education. The total effect of being literate is a 20 % reduction in the odds of VF (95 % CrI: 10 % to 30 %). For post-primary education, the total effect is a 21 % reduction in the odds of VF (95 % CrI: 7 % to 34 %).

FGM had little effect on the odds of VF, after adjusting for outcome misclassification. Women with a short stature had odds of VF that were 31 % (95 % CrI: 2 % to 68 %) higher than their taller counterparts. Among the sample of ever married women (or in a union), the odds of having had VF for those that experienced intimate partner sexual violence were 2.13 times higher than those that never had (95 % CrI: 1.60-2.86). This finding was confirmed in the subsample of 13 surveys that asked all women (never married and ever/currently married) if they had ever experienced sexual violence (from anyone) with an odds ratio of 2.22 (95 % CrI: 1.72-2.90). Among sexually active women, the odds of VF for those that had sexual intercourse before the age of 14 were 41 % (95 % CrI: 16 % to 71 %) higher than those that had a sexual debut at an older age. Expanding our sample by including women that have not begun their sexually active life had little impact on this effect size estimate (odds ratio (OR) = 1.38; 95 % CrI: 1.14-1.66). Both of these findings are in line with the one from age at first birth. Indeed, having had a first live birth before the age of 14 was associated with odds of VF that were 39 % higher (95 % CrI: 1.04-1.82) than those that had their first birth at older ages. Finally, having difficulty obtaining permission to seek health care was associated with increased odds of reporting VF symptoms but this effect did not reach statistical significance (OR = 1.20; 95 % CrI: 0.99-1.47).

[D]{.smallcaps}iscussion {#Sec9}
========================

Main findings {#Sec10}
-------------

Pooling data from up to 27 population-based surveys conducted in sub-Saharan Africa, we identified the following risk factors for VF: illiteracy, short stature, sexual violence, young age at first sexual intercourse, and young age at first birth. These results corroborate findings from previous studies on the importance of some individual-level risk factors for obstetric fistulas \[[@CR12], [@CR13], [@CR15], [@CR43]\]. Short stature, early sexual debut and young age at first birth are risk factors that are related, among other things (e.g., cultural practices, women's status, access to family planning), to a woman's anthropometry. Women that were young at first sex/birth, with immature pelvic bones, and women of short stature have increased incidence of cephalo-pelvic disproportion, which is a known risk for obstructed labor \[[@CR44]\]. Illiteracy has been found to significantly increase the odds of VF and its effect went beyond that mediated by gravidity status. In contrast, we found no significant direct effect of post-primary education on VF occurrence. This could be explained by the fact that literacy was objectively measured whereas the quality of primary education in sub-Saharan Africa varies widely, even within the same geographical region \[[@CR45]\]. Alternatively, it is possible that education beyond primary school has no impact on VF incidence, suggesting that fistula sufferers are the most disadvantaged of the disadvantaged. We did not evidence any relationship between FGM and VF. The DHS/MICS questionnaire, however, did not enable us to investigate if the most severe forms of FGM, such as infibulation and gishiri cutting, are risk factors for VF. Experience of intimate partner violence had a large effect on VF occurrence, as reported previously \[[@CR31]\]. Taken together, these results suggest that empowerment and improvement of women's status could play a key role in reducing the burden of VF in sub-Saharan Africa.

Strengths and limitations {#Sec11}
-------------------------

A number of strengths characterize this study. First, we have conducted what is believed to be the largest population-based analysis of risk factors for VF, pooling data from up to 23 countries (27 surveys) in sub-Saharan Africa. Second, we explicitly modeled uncertainty of self-reports of VF symptoms using a Bayesian Latent C lass model whereas the few other studies that have examined risk factors for VF using DHS/MICS surveys did not perform such adjustments. Finally, we used a doubly-robust method for inferences as a safeguard to bias of effect size estimates.

This study has some limitations. First, some of the estimates could be affected by reverse causality. This is mostly true for intimate partner sexual violence since the majority of other risk factors are likely to have preceded the occurrence of VF symptoms. Hence, it is possible that women with VF have a higher probability of being affected by intimate partner sexual violence as VF impacts their status within marriage and community \[[@CR46]\], for example by creating financial stress and/or affecting women's economic productivity. Living with fistula was found to interfere with sexual activity for 85.2 % of patients in a multi-country study \[[@CR18]\]. Some physical and psychological consequences of VF persist after repair \[[@CR47]\] and this could influence risk of sexual violence \[[@CR18]\]. Second, we could not exclude from our sample fistulas that were not of obstetric origins as many surveys did not record the cause of VF symptoms. Since more than 90 % of VF in sub-Saharan Africa are from obstetric origins \[[@CR10], [@CR13]\], inclusion of VF from other causes should have little impact on our estimated effect size measures. Third, risk factors like intimate sexual violence, the degree of difficulty of obtaining permission to seek health care, and literacy were measured at the time of interview and we assumed these to be time-invariant. This assumption could be violated if these risk factors have changed since the women's onset of VF symptoms. Finally, the cross-sectional nature of the surveys, coupled with potentially important within-country migration, prohibited us from examining the effect of a number of other risk factors such as access to health services, quality of health services, and coverage of maternal health interventions that may ultimately represent important barriers to the prevention and elimination of obstetric fistula in sub-Saharan Africa.

Interpretation {#Sec12}
--------------

VF embodies many of the challenges of the post-2015 agenda, and, more specifically, of the unfinished reproductive health agenda. Despite a decade of maternal health improvements \[[@CR2], [@CR3]\], poor access to and quality of health services is the norm in most low and middle income countries with antenatal and perinatal care being the least equitable interventions \[[@CR48]\]. The third sustainable development goal (SDG) aims at reducing the maternal mortality ratio to less than 70 per 100,000 live birth and to ensure universal access to sexual and reproductive health-care services, including family planning \[[@CR49]\]. The fifth goal also calls for achieving gender equality and women empowerment, with the elimination of all form of violence against women and girls and of harmful practices such as early and forced marriage \[[@CR49]\]. Attaining these objectives could have important synergistic impacts to reduce incidence of obstetric fistula \[[@CR50]\], but quality of care should be emphasized as poor vulnerable women are often attended by "*the most disenfranchised members of the health-care system*" \[[@CR4]\]. The importance of family planning and antenatal care should also be stressed. Indeed, universal access to sexual and reproductive health is emphasized in both the third and fifth SDG. Alongside, access to comprehensive emergency obstetric care should be viewed as a form of prevention \[[@CR51], [@CR52]\]. Yet, our study highlighted that fistula prevention could be most effective if accompanied with enhanced efforts on education and women empowerment.

[C]{.smallcaps}onclusions {#Sec13}
=========================

Our study confirms a number of important individual-level risk factors for VF, while adding precision to the effect size estimates, using population-based data from a large number of countries in sub-Saharan Africa. Increasing literacy, delaying age at first sex/birth, and preventing sexual violence could contribute to the elimination of obstetric fistula if concomitant improvements in access to quality sexual and reproductive healthcare are ensured.

Ethics approval and consent to participate {#Sec14}
------------------------------------------

Informed consent was provided by all survey participants (or their guardian) before questionnaire administration. Further, all DHS survey protocols have been approved by the Internal Review Board of ICF International in Calverton (USA) and by the relevant country authorities for both DHS and MICS. Further information on the ethics approval can be found in the individual country reports published by DHS and MICS.

Consent for publication {#Sec15}
-----------------------

Not applicable.

Availability of data and materials {#Sec16}
----------------------------------

Datasets containing individual-level records are in the public domain and can be obtained from The DHS Program (DHS surveys) and UNICEF (MICS surveys).

[A]{.smallcaps}dditional file {#Sec17}
=============================

Additional file 1: Table S1.Characteristics of the study population, stratified by survey. **Table S2.** Number of surveys included in the analyses, sample size in the un-matched datasets, matching ratio, and matched sample size for the selected risk factors. (DOCX 22 kb)

CrI

:   credible intervals

DHS

:   Demographic and Health Survey

FGM

:   female genital mutilation

MICS

:   Multiple Indicators Cluster Survey

OR

:   odds ratio

SDG

:   Sustainable Development Goals

VF

:   vaginal fistula

**Competing interests**

The authors declare that they have no competing interests.

**Authors' contributions**

FKS developed the original research idea with contributions by SS and MMG. FKS, MMG, and NMa performed the background literature review for this paper. MMG assembled and managed the databases. MMG also performed the analyses with inputs from VF, MCC, MP and FKS. MMG wrote the manuscript with contributions by FKS and NMa. VF, NMa, SS, MCC, NMe, MP, and FKS contributed intellectual content to the paper and critically reviewed it. All authors have read and approved the final version of this manuscript.

MMG's research is supported by the Bisby prize and a postdoctoral fellowship from the *Canadian Institutes of Health Research*. FKS received a travel award from *Fonds de la Recherche Scientifique-FNRS* to complete this work. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
