International Journal of Pediatrics
Volume 2020 (2020), Article ID 3894026, 14 pages
https://doi.org/10.1155/2020/3894026

Research Article

Incidence Density Rate of Neonatal Mortality and Predictors in Sub-Saharan Africa: A Systematic Review and Meta-Analysis

Zebenay Workneh Bitew,¹ Ayinalem Alemu,² Ermias Getaneh Ayele,¹ Desalegn Abebaw Jember,¹ Michael Tamene Haile,¹ and Teshager Worku³

¹St. Paul’s Hospital Millennium Medical College, Nursing Education Directorate, Addis Ababa, Ethiopia
²Ethiopian Public Health Institute, Addis Ababa, Ethiopia
³College of Health and Medical Sciences, School of Nursing and Midwifery, Haramaya University, Harar, Ethiopia

Correspondence should be addressed to Zebenay Workneh Bitew; zedo2015@gmail.com

Received 14 July 2020; Accepted 26 September 2020; Published 17 October 2020

Academic Editor: Asirvatham Alwin Robert

Copyright © 2020 Zebenay Workneh Bitew et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. Neonatal mortality in Sub-Saharan countries is remarkably high. Though there are inconsistent studies about the incidence density rate of neonatal mortalities (IDR) and predictors in Sub-Saharan Africa, they are inconclusive to policymakers and program planners. In this study, the IDR of neonatal mortalities and predictors was determined. Methods. Electronic databases (Web of Science, PubMed, EMBASE (Elsevier), Scopus, CINAHL (EBSCOhost), World Cat, Google Scholar, and Google) were explored. 20 out of 818 studies were included in this study. The IDRs and predictors of neonatal mortality were computed from studies conducted in survival analysis. Fixed and random effect models were used to compute pooled estimates. Subgroup and sensitivity analyses were performed. Results. Neonates were followed for a total of 1,095,611 neonate-days; 67142 neonate-days for neonates treated in neonatal intensive care units and 1,028,469 neonate-days for community-based studies. The IDRs of neonatal mortalities in neonatal intensive care units and in the community were 24.53 and 1.21 per 1000 person-days, respectively. The IDRs of early and late neonatal mortalities neonatal intensive care units were 22.51 and 5.09 per 1000 neonate-days, respectively. Likewise, the IDRs of early and late neonatal mortalities in the community were 0.85 and 0.31, respectively. Not initiating breastfeeding within one hour, multiple births, rural residence, maternal illness, low Apgar score, being preterm, sepsis, asphyxia, and respiratory distress syndrome were independent predictors of time to neonatal mortality in neonatal intensive care units and male gender, perceived small size, multiple births, and ANC were predictors of neonatal mortality in the community. Conclusion. The incidence density rate of neonatal mortality in Sub-Saharan Africa is significantly high. Multiple factors (neonatal and maternal) were found to be independent predictors. Strategies must be designed to address these predictors, and prospective studies could reveal other possible factors of neonatal mortalities.

1. Background

Neonatal mortality rate (NMR) is defined as the death of a newborn within 28 days of birth, and it is expressed per 1000 live births [1]. It is also classified as early (the first 7 days) and late (7 to 28 days) neonatal mortality rates [2, 3]. It is the most vulnerable period in which around one million newborns die in the early neonatal period and 2.8 million die in the late neonatal period [4]. Neonatal mortality rate in the least developed countries (26 per 1000 live births) is significantly lower than NMR in the globe (19 per 1000 live births) [5]. In Sub-Saharan Africa (SSA), NMR was found to be 29 per 1000 live births that accounted for 36% of under-five mortalities [6].

Neonatal mortality (NM) is a major public health problem that endangers the survival of children with a remarkable variation between developed (4 to 46%) and developing countries (0.2 to 64.4%) [7]. It contributes to 44% of under-five mortalities throughout the world and more than 99% of NM took place in low- and middle-income countries, including SSA [2, 3, 8], with a slow progress in decreasing of NM being seen in African and Asian countries [9]. The highest incidence of neonatal mortality was also reported in Sub-Saharan countries, and the burden remained high even after the implementation of the millennium development goals. For instance, NMR in Ethiopia remained stagnant for decades with the current NMR of 30 per 1000 live births [10].

Even though the burden of NMR is considerably high in SSA, there are very few inconsistent studies about the incidence dentistry rate (IDR) of NM. The predictors of time to NM are also very limited. The IDR of NM in SSA ranged from 1.00 per 1000 neonate-days [11] to 39.04 per 1000 neonate-days [12] in Ethiopia and from 1.93 [13] to 30 [14] per 1000 neonate-days in the other Sub-Saharan countries. The burden of IDR of NM is significantly higher in the studies conducted among neonates who were treated in the neonatal intensive care units (NICUs) [12–22] of SSA as compared with community-based studies [11, 23, 24]. The predictors of time to NM were reported by very limited studies that can be classified as service-related (type of delivery, place of delivery, antenatal and postnatal cares, immunization services, and physical distances of health institutions), maternal-related (age, multiple pregnancy, residence, initiation of breast feeding within one hour of birth, delivery complications, and maternal illness), and newborn-related factors (vital signs, low birth weight, gender, and comorbidities like sepsis, hypothermia, asphyxia, respiratory distress syndrome, and jaundice) [11–13, 15–23, 25–30].

All these studies revealed that the IDR of NM remained unacceptably high in SSA. The findings are also inconsistent and inconclusive to policymakers and stakeholders due to variability among the findings. Thus, the pooled IDR of NM and its predictors in SSA was determined in this systematic review and meta-analysis. As to our knowledge, this study is the first of its type to identify the predictors of time to NM using factors that were reported in hazard ratios (HRs). In the original studies, HRs were computed using survival analysis model. The pooled IDRs and predictors were computed for both studies conducted in the NICUs and community-based studies. Hence, the findings of this study will have a vital implication for policymakers, program planners, and to the scientific community. The findings will also have a special importance of depicting gaps hindering the progress of countries towards achieving the sustainable development goals (SDGs), which has a plan of ending preventable causes deaths of newborns by all countries with an aim of reducing neonatal mortality to as low as 12 per 1000 live births by 2030 [31].

2. Methods

2.1. Eligibility Criteria and Information Sources

In this systematic review and meta-analysis, studies conducted in SSA with an objective of assessing incidence density rate (IDR) of neonatal mortality and predictors were included. The studies were assessed using study area, study setups, title, abstract, and full texts prior to inclusion in this study. This study is prepared based on preferred reporting items for systematic reviews and meta-analysis (PRISMA) statement [32]. In the current study, published articles, surveys, and unpublished article that were conducted in English language were explored and included accordingly. Cohort studies that were conducted in clinical and community-based setups reporting the incidence density rates of neonatal mortality and predictors using a hazard function were included in this meta-analysis. Studies that were conducted from inception to June 2020 were searched. However, studies that did not report either IDR of the predictors of neonatal mortality based on the principles of survival analysis or hazard functions were not included in the current meta-analysis. The EndNote X8 reference manager was used to manage retrieved articles.

2.2. Search Strategy and Selection of Studies

Comprehensive searching was conducted to identify the studies. Electronic databases, grey literature sources, and reference lists of articles were explored by three investigators (ZWB, AA, and TW), independently. Web of Science, PubMed, EMBASE (Elsevier), Scopus, CINAHL (EBSCOhost), World Cat, Google Scholar, and Google were explored to find out articles. Searching was conducted using key terms: (a) population (neonates, newborns, and infants), (b) exposure (associated factors, risk factors, determinants, and predictors), (c) outcome (incidence density rate of neonatal mortality, death, mortality, time to death, survival status, and treatment outcomes), (d) study design (cohort studies, epidemiology, observational, and national demographic and health surveys), (e) study setting (hospitals, neonatal intensive care unit, NICU, community-based surveys, and health centers), and (f) location (Sub-Saharan Africa, SSA, and Sub-Saharan countries). The Boolean search operators such as “OR” and “AND” were used during the searching process. Key terms were verified for appropriateness prior to actual searching. Example of search strategy in PubMed: ((((((((((((((((Incidence density rate) AND (“Infant, Newborn”[Mesh])) AND (“Mortality”[Mesh])) OR (“Death”[Mesh])) OR (Time to death)) OR (survival status)) OR (“Treatment Outcome”[Mesh]) AND (predictors))) OR (Determinants)) OR (Associated factors)) OR (“Risk Factors”[Mesh])) AND (“Intensive Care Units, Neonatal”[Mesh])) AND (“Africa South of the Sahara”[Mesh]))).

2.3. Data Extraction Process

Structured and pretested data extraction checklist was used to collect information by three authors (ZWB, DAJ, and EGA). Name of the author(s), publication year, study country, sample sizes, study population, treatment outcome, type of neonatal mortality, follow-up in neonate-days, the incidence density rate of death, and factor in HRs were extracted for this study. During the extraction process, disagreements between the three were solved by the other author (MTH). The data for subgroup analysis based on the type neonatal mortalities and the HRs of the studies were extracted from the included studies. The appropriateness of each datum was verified before the analyses (Additional file 1).

2.4. Quality Assessment of Studies

The qualities of the included studies were independently appraised by two authors (ZWB and MTH) using Joanna Briggs Institute (JBI) Critical Appraisal Checklist for Cohort Studies [33]. The tool has Yes, No, Unknown, and Not Applicable options; 1 is given for Yes and 0 for other options. The minimum score was 0, and 11 was the maximum one. The scores were summed up and changed to percentages. Studies with >50% quality scores were included in this meta-analysis (Additional file 2). The mean scores of the two reviewers were used for the final decision of inclusion of the studies in this systematic review and meta-analysis. During critical appraisal, the third author (TW) participated actively in solving disagreements between two authors.

2.5. Summary Measures

The primary outcome of this study was the incidence density rate of neonatal mortality in SSA in both neonatal intensive care units and in the community setups. The IDR was also computed for early neonatal mortality (birth to 6 days) and late neonatal mortality (7 to 28 days). Time to neonatal mortality was calculated from admission to NICU to death or discharge and from birth to death (for community-based studies). Neonates were followed in days, and the IDRs were calculated by dividing the number of neonates died to total follow-up periods in person-days (neonate-days) and multiplying it by 1000. The IDRs were presented per 1000 per person-days (neonate-days). The second outcomes were predictors of time to neonatal mortality that were computed from studies reporting predictors in the form of HRs. The pooled HRs determining time to death were computed using the crude hazard ratios with 95% confidence intervals (CI) that were reported in the original studies. This was done to control confounders. STATA (version 15) software was used to compute the pooled estimates. The “metan” command was used to compute the pooled hazards ratios of each predictor variable. The pooled hazard ratios (PHRs) were presented with 95% CI. Independent predictors were declared when the 95% CI did not include one. The effect sizes were IDR of neonatal mortality and hazard ratios predicting time to neonatal mortality.

2.6. Statistical Methods and Analysis

In the present meta-analysis, STATA Version 15 (STATA Corporation, College Station Texas) software was used for computing the pooled estimates of both IDR of neonatal mortality and predictors of NM. The pooled estimates were computed using both random and fixed effect models. Results were presented by using either the fixed (when there is no significant heterogeneity among studies) and random effect models (when there is heterogeneity among studies). Subgroup analyses were performed using the type of neonatal mortality (early NM or late NM), study setting (NICU or community based), and country the original studies were conducted. The predictors of neonatal mortality were identified for both clinical-based and community-based studies. The pooled estimates of predictors were also presented separately. The pooled estimates were presented with 95% CI. Meta-analyses were presented using forest plot, summery tables, and texts.

2.7. Publication Bias and Heterogeneity

Publication bias was assessed by looking at the asymmetry of the funnel plot and/or statistical significance of Egger’s regression test. Publication bias was declared when Egger’s regression test was significant () [34]. Heterogeneities among studies were explored using forest plots, test, and Cochrane statistics [35]. The values of 25%, 50%, and 75% were interpreted as low, medium, and high heterogeneity, respectively [36]. In this study, the presence of heterogeneity was declared and justified when with value < 0.05. The sources of possible significant heterogeneities were explored through subgroup analyses and sensitivity analysis (using Duval and Tweedie’s Trim and Fill analysis in the random effect model).

3. Results

3.1. Selection of Eligible Studies

In the initial search, 818 studies were found. These studies were retrieved from electronic databases and other sources. Of these studies, 425 were duplicate files and 342 studies were removed after screening using titles and abstracts. The full texts of 51 studies were reviewed. Finally, 20 studies [11–30] were included in the final analysis of this systematic review and meta-analysis (Figure 1).

Figure 1: PRISMA flow chart showing the literature search process.

3.2. Characteristics of the Studies

All studies included in this study were cohort studies. Regarding the study setting, 12 studies [12–22, 28] were conducted in NICUs and three [11, 23, 24] were community-based prospective cohort studies. The others [25–27, 29, 30] were reported from retrospective analysis of demographic and health surveys of the respective countries. From the included studies to this systematic review and meta-analysis, the majorities (15 articles) [11, 12, 15–24, 28–30] were conducted in Ethiopia, three were performed in Nigeria [25–27], and others were done in Burkina Faso [13] and Uganda [14]. In this study, neonates were followed in days and the follow-up times were reported in neonate-days. Neonates were followed for a total of 1,095,611 neonate days, of which neonates admitted in NICUs of SSA were followed for 67142 neonate-days and neonates in community-based studies were followed for a total of 1,028,469 neonate-days. The IDR of NM in NICU ranged from 1.93 per 1000 neonate-days [13] to 39.04 per 1000 neonate days [12]. Likewise, it ranged from 1.00 [11] to 1.51 [23] per 1000 neonate-days in the community-based studies (Table 1).

Table 1: Incidence density rates (IDRs) of neonatal mortality in Sub-Saharan Africa.

Author year
Country
Sample size
Study population
Treatment outcome
Type of NM
Neonate days
IDR/1000 neonate days (95% CI)
Quality score

Preterm
Term
Dead
Censored
ENM
LNM

Dessu et al., 2019 [15]
Ethiopia
332
37
30
67
265
—
—
2121
31.59 (24.14, 39.04)
High

Gudayu et al., 2020 [16]
Ethiopia
504
—
—
87
417
73
14
2508
34.69 (27.54, 41.84)
High

Kayode et al., 2016 [37]
Ethiopia
1152
—
—
68
1084
50
18
27357.5
2.49 (1.90, 3.08)
High

Orsido et al., 2019 [18]
Ethiopia
964
84
75
159
804
132
27
5889
27.00 (22.86, 31.14)
High

Asmare, 2018 [12]
Ethiopia
571
170
—
170
401
103
67
4354
39.04 (33.28, 44.80)
Medium

Yismaw et al., 2019 [19]
Ethiopia
516
149
—
149
—
127
22
4527
32.91 (27.72, 38.10)
High

Alebel et al., 2020 [20]
Ethiopia
513
42
67
109
404
91
18
4223
25.81 (21.03, 30.59)
High

Mengistu et al., 2020 [21]
Ethiopia
612
69
45
114
498
103
11
4177.8
27.29 (22.35, 32.23)
High

Kahsay et al., 2020 [22]
Ethiopia
253
—
—
32
221
32
—
1425.32
22.45 (14.77, 30.13
Medium

Coulibaly et al., 2016 [13]
Burkina Faso
341
—
—
18
323
—
—
9326.4
1.93 (1.05, 2.81)
Low

Musooko et al., 2014 [14]
Uganda
341
—
—
37
304
37
—
1233.3
30.00 (20.47, 39.53)
Medium

Wosenu et al., 2017 [28]
Ethiopia
490
—
—
171
319
171
—
—
—
Medium

Limaso et al., 2020 [23]
Ethiopia
584
—
584
24
—
15
11
15894
1.51 (0.90, 2.12)
High

Gizaw et al., 2014 [24]
Ethiopia
—
—
—
1055
—
768
287
803370
1.31 (1.23, 1.39)
Medium

Desta et al., 2016 [11]
Ethiopia
7367
—
—
209
7158
144
20
209205
1.00 (0.86, 1.14)
High

Wakgari and Wencheko, 2013 [30]
Ethiopia
8651
—
—
517
—
367
150
—
—
Medium

Ezeh et al., 2014 [26]
Nigeria
27147
—
—
996
26151
—
—
—
—
Medium

Dahiru, 2017 [27]
Nigeria
244836
—
—
8176
236660
8176
—
—
—
Low

Dahiru, 2015 [25]
Nigeria
119024
—
—
3772
115,252
3772
—
—
—
Medium

Mekonen et al., 2013 [29]
Ethiopia
32 042
—
—
—
—
—
—
—
—
Low

In the current study, the predictors of time to NM were presented in the form of hazard ratios. Out of 20 articles included in this systematic review and meta-analysis, 11 articles [12, 13, 15–22, 28] reported the predictors of time to NM among neonates admitted to NICU. Similarly, the predictors of time to NM were reported by seven [11, 23, 25–27, 29, 30] community-based studies, studies conducted through the analysis of secondary data obtained from demographic and health surveys. Regarding the quality scores of articles, nine [11, 15–21, 23], eight [12, 14, 22, 24–26, 28, 30], and three [13, 27, 29] were classified under high, medium, and low qualities, respectively.

3.3. The Pooled Incidence Rate of Neonatal Mortality in Sub-Saharan Africa

Fourteen articles were included to estimate the pooled IDR of NM in SSA. Of these, 11 articles [12–22] were used to compute the pooled IDR of NM among neonates admitted to NICU. Three articles [11, 23, 24] were included to estimate the pooled IDR of NM in the community. The pooled IDR of NM in SSA was found to be 9.85 per 1000 neonate-days (95% CI: 8.65, 11.05; , ). The IDR of NM among neonates admitted to NICUs of SSA was 24.53 per 1000 person-days (95% CI: 18.68, 30.38; , ), whereas the pooled IDR of NM in the community was 1.21 per 1000 person-days (95% CI: 0.94, 1.48; , ) in the random effect model (Figure 2).

Figure 2: The incidence density rate of neonatal mortality in SSA.

3.4. Subgroup Analysis of Incidence Rate of Neonatal Mortality in Sub-Saharan Africa

Due to the presence of higher heterogeneity both in the clinical-based and community-based studies, subgroup analysis was done on the type on NM (ENM and LNM). Twelve studies were eligible to estimate the pooled IDR of ENM in SSA, of which nine [12, 14, 16–22] were used to compute IDR in NICU and three [11, 23, 24] were used for community-based studies. The overall pooled IDR of ENM in SSA was 6.96 per 1000 neonate-days (95% CI: 5.96, 7.96; , ), of which IDR of ENM in NICU was 22.51 per 1000 neonate-days (95% CI: 12.11, 32.91; , ) and 0.85 per 1000 neonate-days (95% CI; 0.62, 1.08; , ) in the community. Likewise, the pooled IDR of LNM (0.92 per1000 neonate-days; 95% CI: 0.62, 1.22; , ) was computed using 10 eligible articles. Of these, the IDR of LNM in NICU was 5.09 per 1000 neonate-days (95% CI: 2.52, 7.65; , ) that was computed from seven eligible studies [12, 16–21]. The IDR of LNM in the community [11, 23, 24] was found to be 0.31 per 1000 neonate-days (95% CI: 0.08, 0.55; , ) (Figure 3).

Figure 3: The IDR of early neonatal mortality (ENM) and late neonatal mortality (LNM) in SSA.

The IDR of NM among neonates admitted to NICUs in SSA was further analyzed based on countries, where the original studies were performed. This was done due to significant heterogeneities among the included studies and to verify the possible source of heterogeneities. Hence, the IDR of NM in Ethiopia, Burkina Faso, and Uganda were 26.96 (95% CI: 14.39, 39.53; , ), 1.93 (95% CI: 1.05, 2.81), and 30.00 (95% CI: 20.47, 39.53) per 1000 neonate-days, respectively (Figure 4). However, heterogeneity among studies remained higher and the possibility of its association with publication bias was explored using funnel plot and Egger’s regression test. The funnel plot revealed that there was publication bias, and the asymmetry of funnel plot was confirmed by a significant Egger’s test () (Figure 5). Finally, sensitivity analysis was done and it was found that no single study affected the pooled IDR. However two studies [13, 17] could have a slight contribution to heterogeneity among the studies used to compute the pooled IDR of NM among neonates treated with in NICUs of countries in Sub-Saharan Africa (Figure 6).

Figure 4: The IDR of neonatal mortality among neonates treated in the NICUs of SSA.

Figure 5: Funnel plot showing publication bias among studies used to compute the IDR of NM among neonates treated in the NICUs of SSA.

Figure 6: Sensitivity analysis for studies included to estimate the pooled incidence density rate of neonatal mortality in NICUs of SSA.

3.5. Predictors of Neonatal Mortality in Sub-Saharan Africa

The predictors of time to NM were computed for both studies conducted in NICUs and studies conducted through analysis of secondary data (community-based studies). For a predictor to be included in this meta-analysis, it has to be statistically significant at least in two studies. Eleven studies [12, 13, 15–22, 28] were included in the analysis of the predictors of time to NM among neonates treated in NICUs of SSA. Fifteen predictors were included in this meta-analysis: low Apgar score, not initiating breastfeeding within one hour of birth, multiple births, CS delivery, home delivery, low birth weight, being preterm, being female, no ANC follow-up, rural residence, sepsis, asphyxia, respiratory distress syndrome (RDS), hypothermia, and maternal illness. Of these CS delivery, home delivery, low birth weight, being female, no ANC follow-up, and hypothermia were found to be statistically insignificant. Low Apgar score, not breastfeeding within one hour of birth, multiple births, being preterm, rural residence, sepsis, asphyxia, RDS, and maternal illness were independent predictors of time to NM. Neonates with low first and 5th minute Agar score were 3.47 (PHR: 3.47, 95% CI: 2.53, 4.41) times more likely to die as compared to neonates who were born with normal Apgar score. The hazard of death of neonates who did not initiate breastfeeding within one hour of birth was 3.77 (PHR: 3.77, 95% CI: 1.44, 6.10) times than the counterparts. Neonates who were part of multiple births were four times (PHR: 3.92, 95% CI: 2.64, 5.37) more likely to die compared to neonates born from single pregnancies. Likewise, preterm neonates were 3.78 (PHR: 3.78, 95% CI: 2.28, 5.27) times more likely to die than term neonates. Neonates who were from a rural residence were nearly two times (PHR: 1.91, 95% CI: 1.30, 2.52) more likely to die than neonates who were from urban residence. The likelihood of death of neonates whose mother had a history of illness during pregnancy was higher than the counterparts (PHR: 1.91, 95% CI: 1.41, 2.41). Moreover, comorbidities (sepsis, asphyxia, and RDS) increased the hazards of death of neonates who were treated in the NICUs of SSA. The hazard of death of neonates with sepsis was 1.94 (PHR: 1.94, 95% CI: 1.53, 2.34) times than neonates free from sepsis. Similarly, the likelihood of death among neonates with perinatal asphyxia were two (PHR: 2.03, 95% CI: 1.32, 2.75) times as compared to neonates who had no perinatal asphyxia. Besides, neonates who had RDS were 2.49 (PHR: 2.49, 95% CI: 1.98, 3.00) times more likely to die than neonates who did not have RDS (Table 2).

Table 2: Predictors of neonatal mortality in Sub-Saharan Africa.

Predictors of neonatal mortality among neonates treated in the neonatal intensive care units (NICUs) of SSA.

Predictors
Included studies
HR (95% CI)
Pooled HR (95% CI)
Heterogeneity

Low Apgar score
Dessu et al., 2019
4.8 (3.61, 9.87)
3.47 (2.53, 4.41)
,

Asmare & Mekonnen, 2018
3.83 (2.70, 5.44)

Alebel et al., 2020
2.70 (1.50, 4.80)

Kahsay et al., 2020
3.13 (1.43, 6.86)

No BF in 1 hour
Dessu et al., 2019
1.79 (1.09, 2.96)
3.77 (1.44, 6.10)
,

Mengesha et al., 2016
21.63 (12.82, 36.49)

Orsido et al., 2019
7.42 (4.73, 11.64)

Alebel et al., 2020
2.00 (1.30, 3.20

Multiple births
Dessu et al., 2019
2.40 (1.19, 4.93
3.92 (2.64, 5.37)
,

Orsido et al., 2019
4.43 (3.24, 6.07)

Mengistu et al., 2020
5.68 (2.95, 6.73)

Kahsay et al., 2020
2.73 (1.18, 6.32)

CS delivery
Dessu et al., 2019
1.90 (1.08, 3.34)
0.99 (0.04, 1.95)
,

Orsido et al., 2019
0.26 (0.15, 0.46)

Mengesha et al., 2016
1.17 (0.65, 2.10)

Home delivery
Gudayu et al., 2020
2.26 (1.04, 4.90)
2.19 (0.97, 3.410)
,

Yismaw et al., 2019
2.14 (1.08, 4.23)

Low birth weight
Mengesha et al., 2016
11.49 (4.49, 29.39)
4.85 (0.68, 10.37)
,

Kahsay et al., 2020
7.30 (3.41, 15.63)

Yismaw et al., 2019
1.56 (1.10, 2.21)

Preterm
Gudayu et al., 2020
2.55 (1.59, 4.08)
3.78 (2.28, 5.27)
,

Orsido et al., 2019
3.41 (2.46, 4.73)

Asmare & Mekonnen, 2018
6.31 (3.89, 10.24)

Kahsay et al., 2020
5.83 (2.52, 13.49)

Coulibaly et al., 2016
17.60 (6.90, 44.8)

Female gender
Orsido et al., 2019
0.59 (0.42, 0.84)
0.87 (0.32, 1.41)
,

Mengistu et al., 2020
1.55 (0.96, 2.01)

Kahsay et al., 2020
0.56 (0.26, 1.19)

No ANC follow-up
Orsido et al., 2019
12.96 (7.86, 21.39)
7.24 (-2.86, 17.33)
,

Alebel et al., 2020
2.60 (1.50, 4.60)

Rural residence
Mengistu et al., 2020
1.98 (1.27, 2.66)
1.91 (1.30, 2.52)
,

Kahsay et al., 2020
1.68 (0.84, 3.36)

Sepsis
Gudayu et al., 2020
1.95 (1.06, 3.59)
1.94 (1.53, 2.34)
,

Asmare & Mekonnen, 2018
2.21 (1.57, 3.12)

Yismaw et al., 2019
1.68 (1.04, 2.72)

Mengistu et al., 2020
2.50 (1.81, 5.27)

Wosenu et al., 2017
1.80 (1.20, 2.60)

Asphyxia
Gudayu et al., 2020
1.23 (0.79, 1.92)
2.03 (1.32, 2.75)
,

Orsido et al., 2019
1.73 (1.25, 2.39)

Yismaw et al., 2019
2.18 (1.58, 3.03)

Wosenu et al., 2017
1.60 (1.05, 2.56)

Mengistu et al., 2020
4.84 (2.64, 5.90)

Respiratory distress syndrome
Orsido et al., 2019
4.60 (3.10, 6.82)
2.49 (1.98, 3.00)
,

Asmare & Mekonnen, 2018
2.27 (1.56, 3.28)

Alebel et al., 2020
2.50 (1.70, 3.70)

Wosenu et al., 2017
2.22 (1.50, 3.30)

Hypothermia
Orsido et al., 2019
4.75 (3.33, 6.78)
3.12 (0.12, 6.13)
,

Yismaw et al., 2019
1.68 (1.00, 2.82)

Maternal illness
Asmare & Mekonnen, 2018
2.38 (1.63, 3.46)
1.91 (1.41, 2.41)
,

Yismaw et al., 2019
1.71 (1.21, 2.40)

Predictors of neonatal mortality in SSA using community-based studies

Male gender
Limaso et al., 2020
2.61 (1.08, 6.29)
1.31 (1.22, 1.41)
,

Wakgari & Wencheko, 2013
1.26 (1.06, 1.50)

Ezeh et al., 2014
1.29 (1.11, 1.51)

Dahiru, 2017
1.25 (1.10, 1.43)

Mekonen et al., 2013
1.42 (1.26, 1.59)

Perceived small size
Desta et al., 2016
5.89 (2.44, 14.22)
1.85 (1.05, 2.65)
,

Wakgari & Wencheko, 2013
1.38 (1.05, 1.82)

Ezeh et al., 2014
2.17 (1.82, 2.58)

Multiple births
Limaso et al., 2020
3.48 (1.04, 11.68)
3.71 (1.72, 5.69)
,

Desta et al., 2016
6.76 (4.64, 9.81)

Wakgari & Wencheko, 2013
3.73 (2.81, 4.94)

Dahiru, 2017
1.82 (1.43, 2.31)

Rural residence
Ezeh et al., 2014
1.36 (1.11, 1.66)
1.30 (1.27, 1.32)
,

Dahiru, 2015
1.30 (1.28, 1.33)

Mekonen et al., 2013
1.22 (1.03, 1.46)

Home delivery
Limaso et al., 2020
2.55 (1.01, 6.42)
1.01 (0.84, 1.17)
,

Ezeh et al., 2014
1.00 (0.84, 1.17)

CS delivery
Ezeh et al., 2014
2.33 (1.54, 3.51)
1.58 (0.91, 2.20)
,

Dahiru, 2017
2.09 (1.59, 2.74)

Dahiru, 2015
0.93 (0.81, 1.07)

Mekonen et al., 2013
1.26 (0.86, 1.85)

Pregnancy complication
Wakgari & Wencheko, 2013
Dahiru, 2015
1.73 (1.27, 2.24)
1.35 (0.71, 1.99)
,

1.07 (1.01, 1.13)

ANC follow-up
Wakgari & Wencheko, 2013
0.72 (0.59, 0.89)
0.82 (0.66, 0.97)
,

Dahiru, 2015
0.88 (0.85, 0.92)

No PNC
Limaso et al., 2020
2.799 (1.16, 6.75)
—
—

BF: breastfeeding; CS: cesarean section; ANC: antenatal care; PNC: postnatal care.

In the current study, seven community-based studies [11, 23, 25–27, 29, 30] were included to compute the PHRs predicting time to death of neonates in the community. From the eligible predictors, three variables (home delivery, CS delivery, and pregnancy complications) were not significantly associated with time to NM, whereas male gender, perceived small size, multiple births, and ANC follow-up were significantly associated with time to NM. Perceived small size is a substitute of low birth weight in NICU-based studies due to the fact that birth weights and weeks of delivery are usually unclear in the community-based studies. Male neonates were 1.31 (PHR: 1.31, 95% CI: 1.22, 1.41) times more likely to die than females. Neonates who were perceived as small in size by their mother were 1.85 (PHR: 1.85, 95% CI: 1.05, 2.65) times highly likely to die as compared to neonates who were perceived to be normal in size at birth. Similarly, multiple births increased the likelihood of death by 3.71 (PHR: 3.71, 95% CI 1.72, 5.69) times than single births. Finally, this study revealed that ANC follow-up could decrease the hazard NM by 12% (PHR: 0.82, 95% CI: 0.66, 0.97). Only one study [23] reported that no postnatal care increased the likelihood of neonatal mortality in the community (HR: 2.80, 95% CI: 1.16, 6.75) (Table 2).

4. Discussion

This study is the first of its type in investigating the pooled incidence density rate of neonatal mortality and its predictors among neonate treated in the neonatal intensive care units and community setups of SSA. Though there are limited studies in SSA about IDR of NM and its predictors, 20 studies were found to fulfill the eligibility criteria of this study. Thus, 14 studies were used to compute the pooled IDR of NM and 18 studies were used to identify the predictors of NM both in the clinical setups and in the community.

The pooled IDRs of NM were computed using fixed and random effect models. However, the final estimates were reported using random effect model, due to the presence of remarkable heterogeneity among the included studies (, ). The pooled IDR of NM in SSA is 9.85 per 1000 neonate-days, of which the IDR of neonates who were treated in NICUs is 24.53 per 1000 person-days and 1.21 per 1000 person-days in the community. Due to high heterogeneity among studies, these pooled estimates were presented with a random effect model. In this study, subgroup analyses were performed to compute the pooled IDR of NM. Based on the type of NM, IDR of ENM is 6.96 per 1000 neonate-days, of which IDR of ENM in NICU is 22.51 per 1000 neonate-days and 0.85 per 1000 neonate-days in the community. Likewise, the IDR of LNM is 0.92 per1000 neonate-days. Of these, IDR in NICU is 5.09 per 1000 neonate-days and the IDR of LNM is 0.31 per 1000 neonate-days. The IDR of NM in NICUs of SSA was further analyzed based on study country. The IDR of NM in Uganda, Ethiopia, and Burkina Faso were 30.00, 26.96, and 1.93 per 1000 neonate-days, respectively. Finally, trim and fill analysis (sensitivity analysis) was performed; however, no significant change was observed in the pooled estimates and heterogeneities between the studies included in this meta-analysis. The possible source of heterogeneity could be best explained by the presence of a wide variation in NM among Sub-Saharan countries that might be related to differences in the quality of healthcare governance, prevalence of HIV, and socioeconomic deprivation [37]. From these findings, it can be deduced that the IDR of NM in most of studies in NICUs is significantly higher than IDR in the community. The pooled IDRs are also remarkably higher NICUs. This could be attributed to the high burden of comorbidities in neonates admitted to NICUs, and most of neonates admitted to NICUs are preterm newborns which could uplift the likelihood of NM [38].

Regarding the predictors of NM, the predictors were classified into two. The first were those predicting time to NM of neonates who were treated in the NICUs of SSA. The second were predictors that determine the time to NM of neonates from birth to death using community-based studies. The independent predictors of NM among neonates who were treated in NICUs were further classified as maternal-related and newborn-related factors. The maternal-related factors were not initiating breastfeeding within one hour of birth, multiple births, rural residence, and maternal illness. The likelihood of NM during delayed initiation of breastfeeding is significantly higher as compared to the early initiation (with in 1 hour of birth) and this is supported by previous meta-analyses [39–41]. Similarly, the hazard of death among neonates born from multiple pregnancies is considerably higher than single pregnancies. The possible elucidation for this is that multiple births are usually associated with preterm births that can amplify the risk of NM [42, 43]. Rural residence is also an independent predictor of NM and the possible rational might be lack of access to health care services during antenatal and postnatal periods which could increase NM. Besides, maternal illness during pregnancy is found to have a significant association with NM and this is supported by a previous study finding [44].

The newborn-related factors were low Apgar score in the first and fifth minute of birth, being preterm, sepsis, asphyxia, and RDS. Low Apgar score in the first and 5th minute of birth is found to increase the likelihood of NM, and it is corroborated by a previous finding that depicted low Apgar scores have a higher predictive value of NM [45]. The hazard of death of preterm babies is remarkably higher than neonates born on the due date. This could be due to the fact that preterm newborns are at greater risk of death, which could be resulted from physical and physiologic immaturities [46, 47]. In addition, comorbidities (sepsis, asphyxia, and RDS) increased the hazards on NM. It is clear that neonatal sepsis is one of the public health problems affecting the survival of children, primarily in countries included in this study [48]. Neonatal sepsis is frequently reported by prior studies as a predictor of NM [2, 49, 50]. Likewise, perinatal asphyxia is found to have a main predictor of NM. This could be due the fact that asphyxia is caused by interruption of blood flow to the placenta that could further decrease blood flow to the vital organs like the kidney, brain, heart, and lung [51, 52]. Another possible explanation for this is that the presence of high prevalence of perinatal asphyxia in East and Central Africa could be associated with increased incidence of NM in SSA [53]. Respiratory distress syndrome is also found to increase the hazard of NM of neonates treated in the NICUs of SSA. The mechanism of NM secondary to RDS is nearly similar with asphyxia, and RDS is the main predictor of NM which is supported by other findings [54, 55].

In the current study, the predictors of NM among neonates born in SSA were male gender, perceived small size, multiple births, and ANC follow-up. These predictors were computed from demographic and health surveys of the respective countries. The likelihood of NM in males is significantly higher than females. This could be associated with physiologic or hormonal differences, and the other possible rationales need to be explained by future researches. Likewise, perceived small size as reported by mothers of neonates is significantly associated with NM. This could be due to the fact that it is highly likely for perceived small sized babies to be either low birth weight or preterm and the likelihood of NM is highly likely in either of the two cases [56]. Multiple births are also independent predictors of NM in the community-based studies similar to studies conducted in the NICUs. In the present study, access to antenatal care is another independent determinate of prohibiting death neonates in SSA. This is supported a prior meta-analytic finding which revealed utilization of at least one antenatal visit could decrease NM by up to 39% [57]. Besides, only one study reported that no postnatal follow-up increased the hazard of NM [23]. There were very limited studies about the effect of postnatal care on the survival of newborns in SSA and it is a potential area of investigation.

In general, the findings of this study pinpointed that NM in SSA is caused by multiple predators, which are informative to program planners and policymakers of perinatal cares. The policymakers need to focus on promotion breastfeeding in the first one hour of birth. The NICUs should be organized to the standard. This can mitigate the three leading causes (asphyxia, preterm-related complications, and sepsis) NM in SSA. Besides, designing policies focusing on universal access of antenatal and postnatal cares may minimize the maternal causes of NM, such as maternal illnesses in SSA.

4.1. Strengths and Limitations

The main strength of this study is that the reputable databases were explored to find all possible articles. This study is the first of its type in determining the IDR of NM in SSA and unfolding the possible predictors of NM using hazard ratios. The predictors were separately estimated from studies in NICUs and community-based studies. This finding will have a paramount importance for program planners and policymakers in SSA, where the burden of NM is considerably high. Despite these strengths, the presence of scant studies in countries other than Ethiopia could obscure some other predictors of time to NM in SSA. Lack of reports on IDR of NM in community-based studies could also be associated with under estimation of IDR of NM in SSA. Prospective cohort studies in countries other than Ethiopia and community-based prospective cohort studies could unfold all possible causes of NM in the future.

5. Conclusion

The incidence density rate of neonatal mortality among neonates treated in neonatal intensive care units of SSA is significantly high. It is also considerable in the community. The hazard of neonatal mortality in NICUs is increased by comorbidities (sepsis, asphyxia, and RDS), low Apgar scores, being preterm, not initiating breastfeeding within one hour of birth, multiple births, rural residence, and maternal illness. Likewise, male gender, perceived small size, multiple births, and ANC follow-up were independent predictors of time NM in the communities of SSA. Designing strategies to address all these predictors could decrease the high burden of NM in SSA. We also recommend prospective longitudinal researches to be done in the respective Sub-Saharan countries to identify all possible bottlenecks of neonatal survival in SSA.

Abbreviations

ANC: | Antenatal care
CI: | Confidence interval
ENM: | Early neonatal mortality
JBI: | Juana Brigg’s Institute
HH: | Hazard ratio
LNM: | Late neonatal mortality
IDR: | Incidence density rate
NICU: | Neonatal intensive care unit
NM: | Neonatal mortality
NMR: | Neonatal mortality rate
PHR: | Pooled hazard ratio
SDGs: | Sustainable development goals
SSA: | Sub-Saharan Africa
RDS: | Respiratory distress syndrome.

Data Availability

The data that support the review findings of this study are included in the manuscript, and the supporting files are submitted as additional files (Additional file 1 & Additional file 2) and PRISMA checklists and search strings.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

All authors contributed equally in designing the study; data extraction, analysis, and interpretation; and final approval of the manuscript.

Acknowledgments

We are thankful to the authors of the original studies included in this systematic review and meta-analysis and to those who contributed a lot in this work.

Supplementary Materials

Supplementary 1. Extraction checklist of the included studies.

Supplementary 2. The critical appraisal of cohort studies.

References

V. Y. Yu, “Global, regional and national perinatal and neonatal mortality,” Journal of Perinatal Medicine, vol. 31, no. 5, pp. 376–379, 2003.

M. J. Sankar, C. K. Natarajan, R. R. Das, R. Agarwal, A. Chandrasekaran, and V. K. Paul, “When do newborns die? A systematic review of timing of overall and cause-specific neonatal deaths in developing countries,” Journal of Perinatology, vol. 36, no. S1, pp. S1–S11, 2016.

H. T. Tran, L. W. Doyle, K. J. Lee, and S. M. Graham, “A systematic review of the burden of neonatal mortality and morbidity in the ASEAN region,” WHO South-East Asia Journal of Public Health, vol. 1, no. 3, pp. 239–248, 2012.

IGME U, Levels & trends in child mortality: report 2017, Estimates developed by the UN Inter-Agency Group for child mortality estimation, United Nations Children’s Fund, New York, 2017.

UNICEF, <SOWC-Statistical-Tables-2017.pdf>, 2018, http://www.unicef.org/sowc.

UN IGME, Levels & trends in child mortality: estimates developed by the UN Interagency Group for child mortality estimation. Report 2015, World Health Organization, Geneva, Switherland, 2015.

S. Chow, R. Chow, M. Popovic et al., “A selected review of the mortality rates of neonatal intensive care units,” Frontiers in Public Health, vol. 3, 2015.

M. Z. Oestergaard, M. Inoue, S. Yoshida et al., “Neonatal mortality levels for 193 countries in 2009 with trends since 1990: a systematic analysis of progress, projections, and priorities,” PLoS Medicine, vol. 8, no. 8, article e1001080, 2011.

The partnership for maternal and child health: newborn death and illness, 2015, http://www.who.int/pmnch/media/press_materials/fs/fs_newborndealth_illness/en/.

Institute EPH and ICF, Ethiopia mini demographic and health survey 2019: key indicators, EPHI and ICF Rockville, Maryland, USA, 2019.

B. N. Desta, N. Assefa, T. D. Damte, and L. O. Hordofa, “Neonatal mortality and its risk factors in eastern Ethiopia: a prospective cohort study in Kersa health and demographic surveillance system (Kersa HDSS),” Mortality, vol. 13, no. 18, 2016.

Y. Asmare, Survival status and predictor of mortality among premature neonate admitted to neonatal intensive care unit from 2013-2017 in Tikur Anbesa Specialized Hospital, Addis Ababa University, Addis Ababa, Ethiopia, 2018.

A. Coulibaly, A. Baguiya, T. Millogo, I. B. Meda, F. Koueta, and S. Kouanda, “Predictors of mortality of low birth weight newborns during the neonatal period: a cohort study in two health districts of Burkina Faso,” International Journal of Gynecology & Obstetrics, vol. 135, no. S1, pp. S89–S92, 2016.

M. Musooko, O. Kakaire, A. Nakimuli et al., “Incidence and risk factors for early neonatal mortality in newborns with severe perinatal morbidity in Uganda,” International Journal of Gynecology & Obstetrics, vol. 127, no. 2, pp. 201–205, 2014.

S. Dessu, M. Kote, F. Gebremeskel, and T. Girum, “Predictors of neonatal mortality among neonates who admitted in neonatal intensive care unit at Arba Minch General Hospital,” Ethiopian Journal of Health Development, vol. 33, no. 1, 2019.

T. W. Gudayu, E. G. Zeleke, and A. M. Lakew, “Time to death and its predictors among neonates admitted in the intensive care unit of the University of Gondar Comprehensive Specialized Hospital, Northwest Ethiopia,” Research and Reports in Neonatology, vol. Volume 10, no. 1, pp. 1–10, 2020.

H. G. Mengesha, A. D. Wuneh, W. T. Lerebo, and T. H. Tekle, “Survival of neonates and predictors of their mortality in Tigray region, northern Ethiopia: prospective cohort study,” BMC Pregnancy and Childbirth, vol. 16, no. 1, p. 202, 2016.

T. T. Orsido, N. A. Asseffa, and T. M. Berheto, “Predictors of neonatal mortality in neonatal intensive care unit at referral hospital in Southern Ethiopia: a retrospective cohort study,” BMC Pregnancy and Childbirth, vol. 19, no. 1, p. 83, 2019.

A. E. Yismaw, A. A. Gelagay, and M. M. Sisay, “Survival and predictors among preterm neonates admitted at University of Gondar comprehensive specialized hospital neonatal intensive care unit, Northwest Ethiopia,” Italian Journal of Pediatrics, vol. 45, no. 1, p. 4, 2019.

A. Alebel, F. Wagnew, P. Petrucka et al., “Neonatal mortality in the neonatal intensive care unit of Debre Markos referral hospital, Northwest Ethiopia: a prospective cohort study,” BMC Pediatrics, vol. 20, no. 1, p. 72, 2020.

B. A. Mengistu, A. E. Yismaw, Z. N. Azene, and M. S. Mihret, “Incidence and predictors of neonatal mortality among neonates admitted in Amhara regional state referral hospitals, Ethiopia: prospective follow up study,” BMC Pediatrics, vol. 20, no. 1, p. 142, 2020.

A. Kahsay, H. Abebe, L. Gebretsadik, and T. Tekle, Survival and predictors of early neonatal death in neonatal intensive care unit of Mekelle general and Ayder comprehensive specialized hospitals, 2018: Prospective Cohort Study, Northern Ethiopia, 2019.

A. A. Limaso, M. H. Dangisso, and D. T. Hibstu, “Neonatal survival and determinants of mortality in Aroresa district, southern Ethiopia: a prospective cohort study,” BMC Pediatrics, vol. 20, no. 1, p. 33, 2020.

M. Gizaw, M. Molla, and W. Mekonnen, “Trends and risk factors for neonatal mortality in Butajira District, South Central Ethiopia, (1987-2008): a prospective cohort study,” BMC Pregnancy and Childbirth, vol. 14, no. 1, 2014.

T. Dahiru, “Determinants of early neonatal mortality in Nigeria: results from 2013 Nigeria DHS,” Journal of Pediatrics & Neonatal Care, vol. 2, no. 5, 2015.

O. K. Ezeh, K. E. Agho, M. J. Dibley, J. Hall, and A. N. Page, “Determinants of neonatal mortality in Nigeria: evidence from the 2008 demographic and health survey,” BMC Public Health, vol. 14, no. 1, 2014.

T. Dahiru, “First-day and early neonatal mortality in Nigeria: a pooled cross-sectional analysis of Nigeria DHS data,” Journal of Advances in Medicine and Medical Research, vol. 19, no. 9, pp. 1–12, 2017.

M. Wesenu, S. Kulkarni, and T. Tilahun, “Modeling determinants of time-to-death in premature infants admitted to neonatal intensive care unit in Jimma University Specialized Hospital,” Annals of Data Science, vol. 4, no. 3, pp. 361–381, 2017.

Y. Mekonnen, B. Tensou, D. S. Telake, T. Degefie, and A. Bekele, “Neonatal mortality in Ethiopia: trends and determinants,” BMC Public Health, vol. 13, no. 1, 2013.

N. Wakgari and E. Wencheko, “Risk factors of neonatal mortality in Ethiopia,” Ethiopian Journal of Health Development, vol. 27, no. 3, pp. 192–199, 2013.

WHO, Sustainable development goal 3, Gütersloh: WHO Reginal Office, 2015, http://www.who.int/sdg/targets/en/.

D. Moher, A. Liberati, J. Tetzlaff, D. G. Altman, and PRISMA Group, “Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement,” Annals of Internal Medicine, vol. 151, no. 4, pp. 264–9, W64, 2009.

Z. Munn, S. Moola, K. Lisy, D. Riitano, and C. Tufanaru, “Methodological guidance for systematic reviews of observational epidemiological studies reporting prevalence and cumulative incidence data,” International Journal of Evidence-Based Healthcare, vol. 13, no. 3, pp. 147–153, 2015.

J. A. C. Sterne and M. Egger, “Funnel plots for detecting bias in meta-analysis: guidelines on choice of axis,” Journal of Clinical Epidemiology, vol. 54, no. 10, pp. 1046–1055, 2001.

G. Rücker, G. Schwarzer, J. R. Carpenter, and M. Schumacher, “Undue reliance on I 2 in assessing heterogeneity may mislead,” BMC Medical Research Methodology, vol. 8, no. 1, 2008.

J. P. T. Higgins and S. G. Thompson, “Quantifying heterogeneity in a meta-analysis,” Statistics in Medicine, vol. 21, no. 11, pp. 1539–1558, 2002.

G. A. Kayode, D. E. Grobbee, M. Amoakoh-Coleman, E. Ansah, O. A. Uthman, and K. Klipstein-Grobusch, “Variation in neonatal mortality and its relation to country characteristics in sub-saharan Africa: an ecological study,” BMJ Global Health, vol. 2, no. 4, p. e000209, 2017.

L. D. Liang, N. Kotadia, L. English et al., “Predictors of mortality in neonates and infants hospitalized with sepsis or serious infections in developing countries: a systematic review,” Frontiers in Pediatrics, vol. 6, 2018.

E. R. Smith, L. Hurt, R. Chowdhury et al., “Delayed breastfeeding initiation and infant survival: a systematic review and meta-analysis,” PLoS One, vol. 12, no. 7, article e0180722, 2017.

A. K. Debes, A. Kohli, N. Walker, K. Edmond, and L. C. Mullany, “Time to initiation of breastfeeding and neonatal mortality and morbidity: a systematic review,” BMC Public Health, vol. 13, Supplement 3, p. S19, 2013.

A. Alebel, G. Dejenu, G. Mullu, N. Abebe, T. Gualu, and S. Eshetie, “Timely initiation of breastfeeding and its association with birth place in Ethiopia: a systematic review and meta-analysis,” International Breastfeeding Journal, vol. 12, no. 1, 2017.

H. S. Ko, J. H. Wie, S. K. Choi, I. Y. Park, Y. G. Park, and J. C. Shin, “Multiple birth rates of Korea and fetal/neonatal/infant mortality in multiple gestation,” PLoS One, vol. 13, no. 8, article e0202318, 2018.

S. Bellizzi, H. Sobel, A. P. Betran, and M. Temmerman, “Early neonatal mortality in twin pregnancy: findings from 60 low- and middle-income countries,” Journal of Global Health, vol. 8, no. 1, article 010404, 2018.

C. V. Ananth and O. Basso, “Impact of pregnancy-induced hypertension on stillbirth and neonatal mortality,” Epidemiology, vol. 21, no. 1, pp. 118–123, 2010.

S. Cnattingius, M. Norman, F. Granath, G. Petersson, O. Stephansson, and T. Frisell, “Apgar score components at 5 minutes: risks and prediction of neonatal mortality,” Paediatric and Perinatal Epidemiology, vol. 31, no. 4, pp. 328–337, 2017.

A. Sania, E. R. Smith, K. Manji et al., “Neonatal and infant mortality risk associated with preterm and small for gestational age births in Tanzania: individual level pooled analysis using the intergrowth standard,” The Journal of Pediatrics, vol. 192, pp. 66–72.e4, 2018.

J. G. Ray, A. L. Park, and D. B. Fell, “Mortality in infants affected by preterm birth and severe small-for-gestational age birth weight,” Pediatrics, vol. 140, no. 6, p. e20171881, 2017.

M. A. Assemie, M. Alene, L. Yismaw et al., “Prevalence of neonatal sepsis in Ethiopia: a systematic review and meta-analysis,” International Journal of Pediatrics, vol. 2020, Article ID 6468492, 9 pages, 2020.

J. C. Reyes, R. O. Ramírez, L. L. Ramos, L. M. Ruiz, E. A. Vázquez, and V. R. Patiño, “Neonatal mortality and associated factors in newborn infants admitted to a neonatal care unit,” Archivos Argentinos de Pediatria, vol. 116, no. 1, pp. 42–48, 2018.

Z. Akbarian-Rad, S. M. Riahi, A. Abdollahi et al., “Neonatal sepsis in Iran: a systematic review and meta-analysis on national prevalence and causative pathogens,” PLoS One, vol. 15, no. 1, article e0227570, 2020.

M. A. Rainaldi and J. M. Perlman, “Pathophysiology of birth asphyxia,” Clinics in Perinatology, vol. 43, no. 3, pp. 409–422, 2016.

J. J. Volpe, Hypoxic-ischemic injury in the term infant, Volpe’s Neurology of the Newborn, pp. 500–509, 2018.

Y. Workineh, A. Semachew, E. Ayalew, W. Animaw, M. Tirfie, and M. Birhanu, “Prevalence of perinatal asphyxia in East and Central Africa: systematic review and meta-analysis,” Heliyon, vol. 6, no. 4, article e03793, 2020.

J. J. M. Wong, M. Jit, R. Sultana et al., “Mortality in pediatric acute respiratory distress syndrome: a systematic review and meta-analysis,” Journal of Intensive Care Medicine, vol. 34, no. 7, pp. 563–571, 2017.

L. R. Schouten, F. Veltkamp, A. P. Bos et al., “Incidence and mortality of acute respiratory distress syndrome in children: a systematic review and meta-analysis,” Critical Care Medicine, vol. 44, no. 4, pp. 819–829, 2016.

T. Marchant, B. Willey, J. Katz et al., “Neonatal mortality risk associated with preterm birth in East Africa, adjusted by weight for gestational age: individual participant level meta-analysis,” PLoS Medicine, vol. 9, no. 8, article e1001292, 2012.

T. Tekelab, C. Chojenta, R. Smith, and D. Loxton, “The impact of antenatal care on neonatal mortality in sub-Saharan Africa: a systematic review and meta-analysis,” PLoS One, vol. 14, no. 9, article e0222566, 2019.