The effects of HIV on fertility by infection duration: evidence from African population cohorts before ART availability: Fertility by duration of HIV infection.

OBJECTIVES
To estimate the relationship between HIV natural history and fertility by duration of infection in East and Southern Africa before the availability of antiretroviral therapy, and assess potential biases in estimates of age-specific sub-fertility when using retrospective birth histories in cross-sectional studies.


DESIGN
Pooled analysis of prospective population-based HIV cohort studies in Masaka (Uganda) Kisesa (Tanzania), and Manicaland (Zimbabwe).


METHODS
Women aged 15-49 who had ever tested for HIV were included. Analyses were censored at antiretroviral treatment roll out. Fertility rate ratios were calculated to see the relationship of duration of HIV infection on fertility, adjusting for background characteristics. Survivorship and misclassification biases on age-specific subfertility estimates from cross-sectional surveys were estimated by reclassifying person time from the cohort data to simulate cross-sectional surveys and comparing fertility rate ratios to true cohort results.


RESULTS
HIV negative and positive women contributed 15,440 births and 86320 person years; and 1,236 births and 11240 thousand person years respectively to the final dataset. Adjusting for age, study site and calendar year, each additional year since HIV sero conversion was associated with a 0.02 (95%CI 0.01-0.03) relative decrease infertility for HIV-positive women. Survivorship and misclassification biases in simulated retrospective birth histories resulted in modest underestimates of sub-fertility by 2-5% for age groups 20-39y.


CONCLUSION
Longer duration of infection is associated with greater relative fertility reduction for HIV-positive women. This should be considered when creating estimates for HIV prevalence among pregnant women and PMTCT need over the course of the HIV epidemic and ART scale-up.


Introduction
The effects of HIV infection on fertility have been extensively studied in generalized HIV epidemic settings in sub-Saharan Africa [1][2][3][4][5][6]. This was of interest for two reasons: first, to forecast the demographic impacts of hyperendemic HIV [7,8] and, second, because HIV prevalence among pregnant women was widely used for estimating general population HIV prevalence levels and trends [9][10][11]. More recently, the need to plan and evaluate prevention of mother-to-child transmission (PMTCT) programmes has further increased the importance of accurate predictions of fertility of HIVpositive women and changes therein.
Existing literature, largely based on analysis of crosssectional data, has demonstrated that the relationship between HIV infection and fertility depends strongly on age. Among young women (age 15-19 years) antenatal care prevalence is higher than general population prevalence because both pregnancy and HIV risk occur among the subset of women who are sexually active, but among older age groups the fertility rate ratio (FRR) among HIV-positive women becomes increasingly lower relative to HIV-negative women [1,12,13].
Currently, the Spectrum model (Avenir Health, Connecticut, USA) uses estimates of the FRR for HIVpositive to HIV-negative women by age-group estimated by Chen and Walker [1] to generate estimates of HIV prevalence among pregnant women and need for PMTCT. However, rather than a direct effect of age, the lower prevalence among older pregnant women may primarily be associated with reduced fertility during later stages of HIV infection [5,[14][15][16][17]. This distinction is potentially important because of its interaction with the stages of the HIVepidemic -during the early exponential growth period of the epidemic, many more women are recently infected, and so HIV-related subfertility will be lower than later in the epidemic, even among older women. Moreover, antiretroviral treatment (ART) is disproportionately provided to those infected the longest and experiencing the most serious clinical symptomsthose who are expected to experience the greatest fertility reductions. If the effects of HIV on fertility are strongly related to the duration of infection, then these two effects may contribute to biased predictions about need for PMTCT services as ART programmes scale up.
Finally, the hypothesized relationship between duration of HIV infection and fertility may influence our ability to estimate the relationship between HIV and fertility. Widely used estimates of age-specific FRRs by HIV status rely on cross-sectional Demographic and Health Survey data to compare fertility over the previous 3 years among HIV-positive and HIV-negative women [1]. This poses two potential biases (Fig. 1). First, it excludes women who do not survive the 3-year period preceding the survey. If duration of infection influences fertility, then this survivorship bias would exclude women with the lowest fertility, resulting in an underestimate of subfertility based on cross-sectional surveys. Second, retrospective analyses assume that the HIV status at the time of the survey is unchanged over the previous 3 years. For women who seroconverted during the 3 years prior to the survey, this misclassifies some HIV-negative person-time as HIV-positive, again potentially overestimating the true fertility of HIV-positive women.
In this analysis, we estimate the relationship between the imputed duration of HIV infection and fertility using data from three prospective general population open cohorts in Uganda, Tanzania and Zimbabwe -all members of the ALPHA network (London School of Hygiene & Tropical Medicine, London, UK) [18]. The objective of this analysis is to estimate the relationship of HIV natural history and fertility in the absence of treatment, and as such we censor the cohort data at the time when ART became available in the population (population-wide fertility trends in these cohorts since ART scale up have been described elsewhere [19]). We use the prospective demographic and HIV surveillance data to empirically quantify the expected magnitude of survivorship and misclassification biases on age-specific subfertility from cross-sectional surveys.

Methods
Sites and setting Data come from three community-based demographic and HIV open-cohort studies. Kisesa (managed by the National Institute for Medical Research Mwanza) located in northwestern Tanzania, was established in 1994 and has a population of around 34 000. It is predominately rural with a small trading centre on the main road. The average HIV prevalence between 1994 and 2010 was 6% [20]. The Manicaland study (managed by the Biomedical Research and Training Institute and Imperial College London) in Zimbabwe was established in 1998. A prospective household census (population size approximately 37 000) and general population cohort survey (10 000-12 000) were initiated in 12 geographically distinct study sites spread across three districts, with follow-up rounds conducted every 2 or 3 years. The Manicaland study sites comprise two small towns, four S70 AIDS 2017, Vol 31 (Suppl 1) Copyright © 2017 Wolters Kluwer Health, Inc. All rights reserved.
agricultural estates, two roadside settlements and four subsistence farming areas. Overall adult HIV prevalence was around 25% in the late 1990s and has declined steadily to around 15% in 2012-2013 [21]. Masaka (managed by MRC/UVRI Uganda Research Unit on AIDS) is situated in rural southwest Uganda and was established in 1989. Its initial population was around 10 000 which then increased to 18 000 when 10 villages were added to the census area in 2000. Average HIV prevalence between 1989 and 2011 was 8% [22].

Fertility data
In Kisesa, there are two sources of data that are used to estimate fertility. At each demographic surveillance round conducted one to two times per year, a proxy respondent is asked whether each woman in the household gave birth since the previous round and the birth outcome. Also, all new members of the household, including newborns, are linked to their mother if she lives in the household. These two pieces of information are reconciled to give the date of delivery of each birth observed in the demographic surveillance site.
In Masaka, there are four sources of data for estimating fertility. At each annual census, women of child-bearing age are asked whether they were pregnant in the previous 12 months and the birth outcome. The names and identification number of the child are recorded on the mother's record. Second, each new member of the household is enumerated during the annual census and the reasons for joining obtained. If the reason is new born, the mother's identification number is recorded on the child's census record. Third, village leaders are asked to report all births in their village on a monthly basis to the study clerks. This information is entered, and any child reported by these recorders but not on census is added to the census file. Fourth, every 3 years, all children aged less than 18 years are asked about their parents to establish/ confirm who they are and their vital status.
In the Manicaland study, survey rounds are conducted every 2-3 years. At each survey round, eligible women are enumerated in a household census and invited to participate in an open-cohort study. Participants report all births since the previous survey round through a structured questionnaire. For women who die between survey rounds, any births occurring since the previous survey round are recorded in a verbal autopsy interview with the next of kin.

HIV data
In Kisesa, the HIV surveys were carried out separately to the demographic surveillance rounds every 2-3 years, and data were linked afterwards using unique personal identifiers. In Masaka, HIV testing was done immediately after demographic surveillance rounds that were used to list those eligible for HIV testing. HIV testing took place in the home for all sites apart from Kisesa where temporary village clinics are used, to which people are transported from their homes. Prior to the availability of antiretroviral therapy, testing protocols used informed consent without disclosure, so that participants did not learn the results of the HIV research tests. In Manicaland, following household census enumeration, research assistants interview eligible individual participants to collect dried blood spot samples, which are transported to and analysed in an offsite laboratory.

Statistical analysis
Imputation of date of seroconversion Calculating the fertility rate by duration of HIV infection requires data about when a woman seroconverted, which is not exactly observed. We generated 100 imputations for the date of seroconversion for each HIV-positive woman. For women who are observed HIV-negative in one survey round and HIV-positive in a subsequent round (seroconverters), we imputed dates of seroconversion from a uniform distribution between the dates of the last negative and first HIV-positive test.
For women who were already HIV positive the first time they were tested in the cohort (prevalent cases), we imputed 100 seroconversion dates from a distribution determined by the convolution of the age-specific HIV incidence rates and the probability of surviving from seroconversion until the woman's latest age at interview.

Fertility rate ratio by duration of infection
Person-time and live births of women of reproductive age (15-49 years) who had ever tested for HIV in the studies were eligible for inclusion in the analysis. HIV-negative person-time for women with no subsequent positive test was assumed to last for up to 5 years past their last negative test, the exact cut-off point was determined by the HIV incidence rates in the sites, defined as the time at which the cumulated probability of becoming infected following the last negative test reached 5%. Data for each cohort were censored at the start of ART combined using Rubin's rules [23] to give confidence intervals (CIs) that reflect the uncertainty about the exact date of seroconversion. Older age at infection pre-ART is associated with a shorter survival time [24] independent of current age [25]. We investigated whether this could also have an effect on subfertility classified by duration of infection (model not shown).
The effects of survivorship and misclassification bias in retrospective survey analysis We quantified the potential magnitude of survivorship and misclassification biases when estimating age-specific subfertility from cross-sectional surveys by using the population cohort data to simulate the 3-year retrospective fertility history analysis and compared the resulting agespecific FRRs with the true FRRs observed in the cohorts.

Ethics statement
Each of the three sites contributing data to the pooled analysis received ethical clearance from the appropriate local ethics review bodies and from the corresponding Institutional Review Boards at relevant collaborating partner universities.

Estimates of HIV subfertility by duration of infection
The dataset compiled for women aged 15-49 years contained 15 451 births and 86 280 person-years to HIV-negative women, 993 births and 9580 person-years to HIV-positive women and 315 births and 2510 personyears with HIV status unknown. Prior to imputation, the latter group comprised the time before a first positive test and person-time in the seroconversion interval (Table 1). Kisesa contributed the most births (54%) and personyears (42%) ( Table 1) (Table 2).
Compared with HIV-negative women, the relative fertility of HIV-positive 20-24-year-olds was 0.72 (95% CI 0.66-0.79), and relative fertility further reduced with age (Table 3, Model 1). The 15-19-year-old HIV-positive women have higher fertility compared with those who are uninfected due to the fact that many women in this age group are not sexually active and therefore are not exposed to HIV.
Including duration of infection in the model showed that each additional year since seroconversion was associated with a 0.979 (95% CI 0.965-0.995) times reduction in fertility for HIV-positive women, adjusted for age, the effect of age at seroconversion, study site and calendar year (Model 2, Table 3). Accounting for duration attenuated the relative fertility of positive women compared with negative women to 0.78 (95% CI 0.70-0.88) and similarly for other age groups (Model 2, Table 3).
Restricting the model to HIV-positive women (not shown) shows that with increasing year of age at seroconversion, there is an increase in the effect of duration on subfertility (FFR 0.997, 95% CI 0.994-0.999).

Estimates of survivorship bias in retrospective surveys
Age-specific subfertility was larger in the ALPHA sites compared with that found by Chen and Walker [1] apart from the 15-19-year age group (Fig. 2a). The reduction in fertility was 3-12% greater in the age groups 20-34 years and somewhat larger at the oldest age groups, for example 41% lower in the 40-44-year age group. However, CIs encompassed Chen and Walker estimates apart from the 40-44-year-old age group. Figure 2a compares the observed subfertility by age in the cohorts (red dots) with the subfertility estimates when analysed using the assumptions of a retrospective cross-sectional survey (blue triangles). Estimates with simulated misclassification and survivorship bias attenuated the subfertility by age by between 2 and 5% in the age groups between 20 and 39 years and 22% in the 40-44-year age group.
There was some evidence for variation of age-specific subfertility by study site with subfertility in Manicaland lower than in Masaka and Kisesa (Fig. 2b).

Discussion
These data show that longer duration of HIV infection is associated with increased subfertility. Estimating agespecific HIV subfertility using retrospective cross-sectional surveys underestimates subfertility, particularly for older ages, due to survivorship bias being more important at longer duration of infection, which corresponds to greater fertility-reducing effects of HIV infection.
Many studies have documented the effect of HIV on fertility and on age-specific subfertility [1,4,12,13] at the  Copyright © 2017 Wolters Kluwer Health, Inc. All rights reserved.
population level during the pre-ART period. A number of studies in sub-Saharan Africa have looked at disease progression in relation to fertility, a case-control study in Uganda found that high viral load was associated with reduced rates of pregnancy and a reduction in live births [5], despite being sexually active and not using contraception. Also, a clinical cohort found that fertility is reduced from the earliest stage of HIV infection with a large reduction in fertility following the progression to AIDS [16] -this finding was adjusted for sexual activity but not for contraceptive use. A clinical cohort study in Tanzania also found reduced fertility related to clinical stage of HIV [17] adjusting for social and demographic characteristics. A multisite HIV care and treatment programme analysis showed a strong association between disease progression and a reduction in the incidence of pregnancy [15].
Increased subfertility by duration of infection at the population level could have both biological and behavioural factors. Biologically, as well as increases in viral load or decreases in CD4 þ cell count as explanatory factors, the semen quality of HIV-positive partners could be reduced over the time of their infection [26][27][28] or their increased illness could impact on their sexual activity. In terms of behaviour, HIV-positive women are more likely to be widowed [6,29,30] due to having had an HIV-positive partner. Although voluntary testing and counselling was rare in these sites prior to ART introduction, suspicion of HIV status or illness in a partner with HIV may reduce the desire for more pregnancies [31], which may be more obvious at longer durations of infection and it may also increase divorce or separation [6,30].
Increased age at seroconversion accelerated the effects of infection duration on subfertility. Older age at infection leads to shorter survival postinfection [24,25], so a shorter duration to low CD4 þ cell count and higher viral load have been shown to reduce fertility. Also, at older ages of seroconversion, it is more likely that the partner (who is more likely to be older) has been infected for a longer duration; therefore, there is a higher chance of widowhood early on in the women's HIV infection lowering her changes of pregnancy. Finally, older women are likely to have higher parity and therefore may have lower desires for more children than a younger woman who has none or few children.
Compared with the demographic and health survey (DHS) analysis by Chen and Walker [1], ALPHA cohorts showed greater fertility reductions among HIV-positive women by 5-year age group, particularly in the older age groups. Around half of this discrepancy was explained by biases inherent in estimating subfertility from cross-sectional data due to not including the person-years and births of those who died prior to interview and classifying all person-years according to the HIV status at time of interview.

S74
AIDS 2017, Vol 31 (Suppl 1)  [32], which may contribute to these differences [6]. Deliveries and the deaths of children dying in early infancy (particularly in the neonatal period) could be underreported in the ALPHA studies due to recall bias or lack of knowledge on the part of a proxy respondent, which would affect HIV-positive women disproportionally due to the high infant mortality of children infected through vertical transmission [33]. This could artificially increase subfertility estimates in the cohort studies. The DHS will be prone to more recall bias than the cohort studies; however, if analysis is limited to the first few years prior to the interview and the respondent is the women rather than a proxy, it is possible that this will lead to less bias in reporting of births to infants who have died in DHS compared with ALPHA studies. We find that subfertility increases with duration of HIV infection in the absence of ART. This has two important implications that should be considered in future HIVepidemic estimates and the estimates of need for PMTCT. First, over the course of the epidemic, the distribution of duration of infection changes. During the exponential growth phase, a higher proportion of women will be recently infected, and as incidence declines average duration of infection will become longer. This means that the population-level effects of HIV on fertility, and hence the relationship between HIV prevalence measured among pregnant women and general population prevalence, will change.
Second, initiation of antiretroviral treatment has been disproportionately among women in later stages of infection who might be expected to have the lowest fertility rates. Thus, following ART scale up, not only might women on ART have increased fertility [34], but also the fertility of untreated HIV-positive women may be higher because those who would have the lowest fertility are selectively removed into the treatment group. Implementation of Option Bþ over the past several years, in which all pregnant women are initiated on lifelong ART, will further change these dynamics. In light of the demonstrated association between duration of infection and fertility reduction, we recommend that model-based approaches account for not only age but also stage of infection and ART status when estimating HIV prevalence among pregnant women and PMTCT need.
Our results also imply that there are differences in fertility by setting. This underscores that, where possible, locally available data such as prevalence from routine HIV testing of pregnant women should be used in place of default model values to inform appropriate model assumptions about subfertility when generating estimates of PMTCT need.
Finally, it is worth noting that survivorship bias will be less important in the era of ART, as HIV mortality is lower.
The assumption that women who are HIV positive at the time of interview have been infected for at least 3 years will also become more realistic as longer durations of infection become more common in the era of ART. These factors should also be considered when interpreting changes over time in the relationship between HIV and fertility from cross-sectional surveys.