Article Text

Download PDFPDF

Original research
Prevalence and outcomes of twin pregnancies in Botswana: a national birth outcomes surveillance study
  1. Arielle Isaacson1,2,
  2. Modiegi Diseko2,
  3. Gloria Mayondi2,
  4. Judith Mabuta2,
  5. Sonya Davey2,3,
  6. Mompati Mmalane2,
  7. Joseph Makhema2,
  8. Denise L Jacobson4,
  9. Rebecca Luckett2,5,
  10. Roger L Shapiro2,6,
  11. Rebecca Zash2,5
  1. 1Harvard Medical School, Boston, Massachusetts, USA
  2. 2Botswana-Harvard AIDS Institute Partnership, Gaborone, Botswana
  3. 3Brigham and Women’s Hospital, Boston, Massachusetts, USA
  4. 4Center for Biostatistics and AIDS Research, Harvard University T H Chan School of Public Health, Boston, Massachusetts, USA
  5. 5Beth Israel Deaconess Medical Center, Boston, Massachusetts, USA
  6. 6Harvard University T H Chan School of Public Health, Boston, Massachusetts, USA
  1. Correspondence to Arielle Isaacson; arielleisaacson{at}


Objectives This study aims to evaluate the prevalence and outcome of twin pregnancies in Botswana.

Setting The Tsepamo Study conducted birth outcomes surveillance at 8 government-run hospitals (~45% of all births in Botswana) from August 2014 to June 2018 and expanded to 18 hospitals (~70% of all births in Botswana) from July 2018 to March 2019.

Participants Data were collected for all live-born and stillborn in-hospital deliveries with a gestational age (GA) greater than 24 weeks. This analysis included 117 593 singleton and 3718 twin infants (1859 sets (1.6%)) born to 119 477 women between August 2014 and March 2019 and excluded 73 higher order multiples (23 sets of triplets and 1 set of quadruplets).

Outcomes measured Our primary outcomes were preterm delivery (<37 weeks GA), very preterm delivery (<32 weeks GA) and stillbirth (APGAR (Appearance, Pulse, Grimace, Activity, Respiration) score of 0, 0, 0).

Results Women with twin pregnancies had a similar median number of antenatal care visits (9 vs 10), but were more likely to deliver in a tertiary centre (54.8% vs 45.1%, p<0.001) and more likely to have a cesarean-section (54.6% vs 22.0%, p<0.001) than women with singletons. Compared with singletons, twin pregnancies had a higher risk of preterm delivery (<37 weeks GA) (47.6% vs 16.7%, adjusted risk ratio (aRR) 2.8, 95% CI 2.7 to 2.9) and very preterm delivery (<32 weeks) (11.8% vs 4.0%, aRR 3.0 95% CI 2.6 to 3.4). Among all twin pregnancies, 128 (6.9%) had at least one stillborn infant compared with 2845 (2.4%) stillbirths among singletons (aRR 2.8, 95% CI 2.3 to 3.3).

Conclusion Adverse birth outcomes are common among twins in Botswana, and are often severe. Interventions that allow for earlier identification of twin gestation and improved antenatal management of twin pregnancies may improve infant and child survival.

  • community child health
  • public health
  • obstetrics

Data availability statement

Data are available on reasonable request. Per Botswana IRB regulations, deidentified data from the Tsepamo study can be requested from the principal investigator of the study.

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See:

Statistics from

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • Our study uses a large, nationally representative sample with little missing data.

  • Due to the limited availability of early prenatal ultrasound, gestational age may be less accurate and we are unable to determine the cause of preterm birth or whether preterm birth is spontaneous.

  • We are unable to evaluate the impact of preterm and low birth weight on longer-term outcomes of twins, or estimate perinatal mortality.


Twin pregnancies are universally considered to be ‘high risk’. Maternal complications, including hypertensive disorders, anaemia, postpartum haemorrhage and maternal mortality, are more common among twin pregnancies than singleton pregnancies.1 Twins are also more likely than singletons to be born preterm and to have restricted growth in utero, thereby increasing their risk for intrauterine demise and neonatal mortality.2 The majority of existing research on multigravid pregnancies and birth outcomes among twins is from high-resource settings.3 In these settings, the perinatal mortality rate of twins is 3–7 times higher than that of singletons4 5 and up to 60% of twins are born preterm (before 37 weeks).6 7

In most high-resource settings, twin pregnancies are managed by routine prenatal monitoring and neonatal intensive care services, including early identification of twin pregnancies, serial antenatal ultrasound, referral to specialised centres and prenatal counselling.8 In lower-resource settings, data are limited, but twin pregnancies may be at particularly high risk for adverse birth outcomes because of the lack of routine intensive prenatal monitoring and neonatal intensive care services.1 9 Additionally, maternal outcomes with multigravid pregnancies may also be worse because of limitations in management of maternal conditions more common in twin pregnancies, such as postpartum haemorrhage and pre-eclampsia.3

In Southern Africa, where assisted reproductive technology is not commonly available, the incidence of naturally occurring twins is estimated to be high (12–18 per 1000 births) compared with other low-income and middle-income countries (LMIC) in East Asia and Latin America where the incidence of twins is as low as 6–9 per 1000 births.10 This relatively large number of twins may contribute substantially to perinatal mortality in the region.11 Increased prevalence of preterm delivery and low birth weight (LBW) among twins leads to increased risk of under-5 mortality due to malnutrition, respiratory disorders, vulnerability to infection and developmental delays.12 13 It is estimated that one in five twins born in sub-Saharan Africa dies before the age of 5.11 Prior studies of twins in sub-Saharan Africa predate implementation of the United Nations Sustainable Development Goals in 2015 and do not evaluate specific birth outcomes.9–12 This study focuses on Botswana, a country with a rapidly growing economy, strong investment in healthcare, a national programme to prevent maternal mortality,14 15 and a large, nationally representative birth outcomes surveillance study (Tsepamo). We aim to provide the first published data on the prevalence and outcome of twin pregnancies in Botswana.


The Tsepamo Study

The Tsepamo Study is a birth outcomes surveillance study in Botswana, where >95% of women deliver in healthcare settings (not at home).16 Details of methodology have been previously published.17 18 In summary, deidentified information was abstracted from obstetric cards (antenatal care records used throughout the pregnancy) at the time of discharge from the postnatal ward from women who deliver live-born or stillborn infants at select government maternity hospitals in Botswana. From August 2014 to June 2018, data were collected from 8 sites across the country (~45% of births in Botswana) and from July 2018 to March 2019 data were collected from up to 18 sites (~72% of births in Botswana).

At each site, data were collected for all in-hospital deliveries with a gestational age (GA) greater than 24 weeks. Information included maternal demographics, antenatal care visits, ultrasound reports with date of ultrasound, HIV status, method of delivery and infant delivery characteristics (including the number of infants delivered, GA at delivery, birth weight and vital status at birth and discharge).


The primary adverse outcomes assessed were stillbirth Appearance, Pulse, Grimace, Activity, Respiration (APGAR) scores 0, 0, 0), preterm birth (<37 weeks gestation), very preterm birth (<32 weeks gestation), LBW (<2500 g) and very LBW (VLBW) (<1500 g). While APGAR scores and birth weights were recorded for each individual twin, neonatal death status was only collected on the first twin per the original protocol of our study. Therefore, neonatal deaths in twins could not be analysed and were not included as a primary outcome assessed. We chose not to provide data on neonatal deaths in twin 1 without data on twin 2 because it would only provide an incomplete comparison to singletons, which could over or underestimate the total Neonatal deaths (NND) among twins. The GA was documented by midwives at the time of delivery using the estimated delivery date, which is calculated at the first antenatal care visit based on the last reported menstrual period and confirmed by ultrasound when available. If the last menstrual period was unknown or suspected to be incorrect, and if no ultrasound data were available, midwives occasionally used fundal height measurements to estimate GA.

Statistical analysis

The prevalence of each birth outcome was calculated for singletons and twins and higher order multiples were excluded from analyses. Prevalence of preterm birth and very preterm birth were calculated per pregnancy as both twins were born at the same GA. Prevalence of stillbirth, LBW and VLBW among twins were calculated by pregnancy (defined as either twin with the outcome), and by infant (defined as the outcome in the total number of individual infants). When comparing birth outcomes between singletons vs twins, we used pregnancy as the unit of analysis (eg, how many pregnancies ended in at least one stillbirth). Log binomial regression models were fit to determine the relative risk (RR), adjusted risk ratio (aRR) and 95% CIs of adverse birth outcomes among twin compared with singleton pregnancies. Multivariable models were adjusted for maternal age, gravida, educational attainment and maternal HIV status, which were chosen a priori based on prior analysis of risk factors for adverse birth outcomes in Tsepamo.17–19 Statistical analyses were performed using SAS software (Version 9.2).

Patient and public involvement

Patients and the public were not involved in the design or conduct of this study.


Study population

Between August 2014 and March 2019, there were 121 385 infants born to 119 477 women who delivered in the Tsepamo Study, including 117 593 (98.4%) singletons, 3718 twin infants (1859 sets of twins) (1.6%), 69 triplet infants (23 sets of triplets) (0.00%), 4 quadruplets (1 set) (0.00%) and 1 with missing data for delivery number. GA was known in 98.8% of singletons and 98.7% of twins.

Maternal characteristics and obstetric care

Maternal demographics and obstetric care parameters are shown in table 1. Women with twin pregnancies were older, less likely to be primigravid and more likely to have had >4 prior pregnancies compared with women with singleton pregnancies. The median number of antenatal care visits was similar among women with singletons (10) and women with twins (9). Prenatal ultrasound was performed in 78.5% (980/1248) of women with twin pregnancies at a median GA of 25 (IQR 19–31) weeks (20.3% prior to 20 weeks GA). In comparison, 66.9% (52 012/77 786) of women with singleton pregnancies received an ultrasound scan, at a median GA of 27 (IQR 20–33) weeks (14.6% prior to 20 weeks GA). Twins were more likely to be delivered in a tertiary hospital (54.8% vs 45.1%) and also more likely to deliver via caesarean section (C-section) (54.6% vs 22.0%).

Table 1

Maternal characteristics and obstetric care

Birth outcomes

The median GA was 39 weeks (IQR 37–40), and the median birth weight was 3080 g (IQR 2750–3400) among all singletons with 114 749 (97.6%) infants live-born. Among twin pregnancies, the median GA was 37 weeks (IQR 34–38) and the median birth weight was 2330 g (IQR 1920–2675) for the first twin and 2290 g (IQR 1830–2610) for the second twin with 3552 (95.5%) infants live-born. The mean head circumference was 32.4 cm for both the first and second twin, and slightly lower than the mean head circumference for singletons (34.2 cm). Mean length at birth was similar for the first and second twin (45.7 cm and 45.4 cm, respectively) but lower than the mean length of singletons (50.1 cm).

Twin pregnancies were more likely to result in preterm birth than singleton pregnancies (table 2). Compared with singletons, twin pregnancies had a higher risk of preterm birth (47.6% vs 16.7%, aRR 2.8, 95% CI 2.7 to 3.2.9) and very preterm birth (11.8% vs 4.0%, aRR 3.0, 95% CI 2.6 to 3.4). C-sections were more common in preterm twins than preterm singletons (47.7% vs 23.0%), however, among very preterm births (<32 weeks GA), C-section rate was similar (24.9% vs 22.7%).

Table 2

Preterm birth among twin and singleton pregnancies

Among all 3718 individual twin infants, 2397 (64.5%) were LBW, 508 (13.7%) were VLBW and 166 (4.5%) were stillborn. Both twins were LBW in 52% of pregnancies, VLBW in 10.6% of pregnancies and stillborn in 2.0%. The first born twin and the second born twin had similar rates of low and VLBW, but twin 2 was more likely to be stillborn than twin 1 (5.2% vs 3.8%) (table 3). The risk of at least one infant being LBW (aRR 4.5, 95% CI 4.3 to 4.6),VLBW (aRR 5.2, 95% CI 4.7 to 5.8) or stillbirth (aRR 2.8, 95% CI 2.3 to 3.3) was higher among twin pregnancies than singleton pregnancies (table 3).

Table 3

Birth outcomes among twins by birth order


We performed the first published analysis of the prevalence and outcomes of twin births in Botswana using nationally representative data from 2014 to 2019. Because of the absence of medically assisted reproduction, Botswana represents a particularly valuable opportunity to study naturally occurring rates and outcomes of twinning. We found that twin births occurred in 16 per 1000 pregnancies, a prevalence much higher than reported twinning rates in LMIC outside of sub-Saharan Africa. We also found that adverse birth outcomes were high among twin pregnancies, including a 6.9% stillbirth prevalence, almost 3-fold higher than stillbirths in singleton pregnancies (2.4%).

The rate of twinning in our study, 16/1000 pregnancies, is consistent with previously reported high rates of naturally occurring twins in Southern Africa (12–15/1000 pregnancies in Namibia, South Africa and Lesotho and 15–18/1000 pregnancies in Zimbabwe).10 In contrast, very low twinning rates (6–9/1000 pregnancies) have been documented in LMIC across East Asia and Central and South America.10 Geographical differences in naturally occurring twin birth rates are primarily due to genetic, racial and ethnic differences in predisposition to dizygotic (DZ) twin birth.10 While the rate of twinning in Botswana is high compared with other LMIC, it is lower than that in many high-income countries (HICs), where the introduction of medically assisted reproduction has led to dramatic increases in twinning rates over the past several decades. In the USA, for example, the rate of twin births rose from 18.9 per 1000 births in 1980–33.9 per 1000 births in 2014.20 21

The prevalence of stillbirths among twin pregnancies in Botswana (6.9%) is well above the WHO’s target of under 12 stillbirths per 1000 births.22 23 However, our prevalence is similar to that reported across sub-Saharan Africa, and lower than a reported 10.2% stillbirth prevalence among twins in one Nigerian study.24 In contrast, the prevalence of stillbirth in twins in Botswana is higher than in HICs such as the USA (0.35%)25 and South Korea (0.48%).26 The prevalence of stillbirth among singletons in our study is 2.4%, which is nearly 10 times greater than that among singletons in many HICs.27 28 While twin pregnancies are clearly particularly high-risk and contribute substantially to perinatal mortality, the high prevalence of stillbirths among all pregnancies in Botswana highlights a strong need for interventions to decrease stillbirths in the entire population.

The prevalence of preterm birth among twin pregnancies (46.%) in Botswana is substantially higher than that among singleton pregnancies (16.8%) and slightly higher than reported prevalence of preterm birth among twins in other sub-Saharan countries.24 However, prevalence of preterm birth with twin pregnancy in Botswana is lower than the USA, where 60% of twin pregnancies have preterm births.6 28 While the high prevalence of preterm births among twins in the USA may be partially attributed to the increased risk for preterm birth in multiple gestation, it may also be explained by clearer antenatal care guidelines and closer prenatal monitoring among twin pregnancies in the USA compared with Botswana. The high prevalence of preterm birth among twins in the USA is in line with the American College of Obstetrics and Gynecology’s guidelines, which recommend timing delivery based on chorionicity (38 weeks for dichorionic-diamniotic, 34–36 6/7 weeks for monochorionic (MC)-diamniotic, 32–34 for MC-monoamniotic) and delivery by 37–38 week’s gestation in uncomplicated twin pregnancies to reduce the risk of stillbirth.29–32

The relatively low rate of preterm birth and high rate of stillbirth among twin pregnancies in Botswana raises the question of whether twin stillbirths in Botswana could be decreased by increasing delivery of high risk (MC-monoamniotic, MC-diamniotic) twin pregnancies between 34 and 37 weeks.31–33 To implement this type of change, improvements in proportion of ultrasound in the first trimester of twin pregnancy to identify chorionicity and to improve accuracy of GA dating would be needed. Additionally, preterm delivery may carry increased risk in low-resource settings like Botswana where intensive neonatal care services are limited. Though delivery by 37 weeks may not be appropriate in all settings, closer surveillance of twin pregnancies after 34 weeks GA, at which time the risk of stillbirth and maternal complications begins to rise, may improve twin outcomes in Botswana.26 Close monitoring for fetal growth concordance,34 screening for aneuploidy, which is more common among twin pregnancies and associated with higher risk of stillbirth,35 36 and clear delivery planning during antenatal care may also help to reduce adverse perinatal outcomes among twins. In addition to increased antenatal care, improvements in maternal nutrition may help improve outcomes in twin pregnancies in Botswana. Despite the relatively lower prevalence of preterm births in Botswana, we found very high rates of LBW and VLBW among twins (77.2% and 16.8%, respectively, compared with 55.6% and 9.1% in the USA),7 which is likely multifactorial, but suggests that women with twin pregnancies in Botswana may not receive adequate nutrition to support multiple gestations.37 38

Strengths of this study include the large, nationally representative sample with little missing data. Our study also has several limitations. Due to the limited availability of early prenatal ultrasound, GA may be less accurate. Also, our study did not collect information necessary to determine the aetiology of preterm delivery and stillbirth among twins, such as induction status, spontaneity of preterm delivery, whether c-sections were planned or emergent, indication for c-section (including presentation of infants), or fetal heartbeat on admission. Understanding the causes of adverse outcomes is clearly necessary before effective interventions can be designed and implemented. However, we hope our findings will be the catalyst for further research to elucidate these causes and ultimately lead to fewer adverse birth outcomes among twins. We are also unable to distinguish the proportion of twins that are MC, which is a known risk factor for stillbirth and neonatal death.9 33 However, the prevalence of MC twins is thought to be relatively constant worldwide (prevalence of DZ twins varies), so this is unlikely to explain the high prevalence of stillbirth that we found among twin pregnancies in Botswana.39 40 Due to our study design, we are unable to evaluate the impact of preterm and LBW on longer-term outcomes of twins, or estimate perinatal mortality. Finally, we do not have data from deliveries occurring outside the hospital, though this is rare (<5%) in Botswana.15 16


In summary, we found that the prevalence of twinning is high in Botswana and adverse birth outcomes are common among twin pregnancies and are often severe. The high rate of stillborn and LBW twins despite a relatively low preterm birthrate compared with twins in HICs points to the need for interventions that include early ultrasound and identification of twin gestation and type of twin gestation, improved maternal nutrition, close antepartum surveillance and advanced delivery planning in order to reduce morbidity and mortality among twin pregnancies in Botswana.

Data availability statement

Data are available on reasonable request. Per Botswana IRB regulations, deidentified data from the Tsepamo study can be requested from the principal investigator of the study.

Ethics statements

Patient consent for publication


We thank our research assistants Cynthia Dube, Daphne Segobye, Gosego Legase, Keemenao France, Mmapula Ofhentse, Naledi Kamanga, Onkabetse Mokgosi, Rosemary Moremi, Shally Morgan, Tsaone Gaonakala, Tshepang Motlotlegi, Edith Moseki, Patricia Mophutegi, Keba Rabasiako, Nametsang Tshosa, Maipelo Kegakilwe, Masego Kgafela, Tshegofato Motladile, Tsholofelo Tsokunyane, Kealeboga Mmokele, Obakeng Makalane, Thuto Rabana, Seele Mafokate, Annah Bojang, Tlhabologo Baitsemi, Priscilla Mashona and Bathoba Mabiletsa; the maternity staff and administrators at the 18 participating hospitals; the members of the Botswana Ministry of Health and Wellness and the Department of Maternal and Child Health.



  • Contributors AI and RZ designed and performed the analysis and drafted the initial manuscript. RL, JoM, MM, RZ and DLJ contributed to the design, data collection and study oversight of the parent study. GM, MD, JuM, SD and AI oversaw data collection and data cleaning. All authors (AI, MD, GM, JuM, SD, MM, JoM, DLJ, RL, RLS and RZ) contributed substantive feedback, revised and approved the final manuscript. RZ is responsible for the overall content as the guarantor.

  • Funding Supported by grants from the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), National Institutes of Health (R01 HD080471 and R01 HD095766, to RLS; and K23 HD088230, to RZ).

  • Competing interests None declared.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Provenance and peer review Not commissioned; externally peer reviewed.