Incidence and predictors of mortality among neonates admitted with perinatal asphyxia at west Oromia tertiary hospitals, Ethiopia, 2022

Background Perinatal asphyxia is one of the preventable and treatable causes of neonatal mortality. However, it is the fifth-largest cause of under-five mortality. Even with management advancements, it remains one of the key public health issues in underdeveloped countries, including Ethiopia. Comorbidities are also understated; therefore, adequate information regarding the incidence of death and its predictors is required. Methods A four-year retrospective follow-up study was conducted from October 3 to November 2, 2022. From a total sample size, of 655, 616 data were collected by nurse through follow-up reviews charts using Kobo Toolbox software. The data was exported to STATA Version 14 for analysis. The Cox proportional hazard assumption was checked, and the model for the data was selected using Akaike Information Criteria. Finally, an adjusted hazard ratio with 95% CI was computed, and variables with a P-value < 0.05 in the multivariable analysis were taken as significant predictors of death. Result The overall incidence of mortality was 38.86/1000 (95% CI: 33.85–44.60). The median time of follow-up was 15 days (95% CI: 14–20). The proportion of deaths was 202 (32.79%, 95% CI: 29.18–36.61) among neonates with perinatal asphyxia. While the distance from health facility > 10 km is (AHR: 2.25; 95% CI: 1.60–3.17), direct oxygen (AHR: 1.83; 95% CI: 1.35–2.48), APGAR score (Appearance, Pulse, Grimace, Activity, and Respiration) < 3 at the fifth minute (AHR: 2.63; 95% CI: 1.03–6.73), prolonged rupture of membrane (AHR: 1.41; 95% CI: 1.02–1.94), and stage III hypoxic ischemic encephalopathy (AHR: 2.02; 95% CI: 1.18–3.47) were predictors of mortality among neonates with perinatal asphyxia. Conclusion According to this study’s findings, high neonatal mortality due to perinatal asphyxia requires proper intervention regarding membrane rupture, APGAR score (Appearance, Pulse, Grimace, Activity, and Respiration), oxygen use, stage III hypoxic-ischemic encephalopathy, and residence distance.


Background
Perinatal asphyxia (PNA) is the failure of a neonate to establish and sustain normal respiration after birth [1].Perinatal asphyxia causes progressive hypoxemia, hypercarbia, and metabolic acidosis, as well as multi-organ failure [1,2].Perinatal asphyxia is an injury to the fetus or neonate caused by inadequate oxygenation (hypoxia) or inadequate blood flow (ischemia) to numerous organs and systems, including the central nervous system, cardiovascular system, respiratory system, and hematologic system [3,4].
Globally, 2.4 million infant deaths occurred in the first month of life in 2020 [5].Neonatal mortality declined less quickly between 1990 and 2020 than under-five mortality, and the highest neonatal mortality rate in 2020 was on the African continent, with 27 deaths per 1000 live births, and perinatal asphyxia was the major cause of neonatal death on this continent [6,7].After pneumonia, diarrhea, neonatal infections, and preterm birth problems, perinatal asphyxia ranks fifth in terms of the causes of under-five infant mortality, and it accounts for about 23% of neonatal deaths [8,9].Nearly 60% of all neonatal deaths occur in the first three days of life, and the majority of preterm, malformation, and asphyxia-related deaths occur in the first week [10].Perinatal asphyxia rates are higher in developing countries, ranging from 4.6 to 1000 in Cape Town to 7-26 per 1000 in Nigeria, and mortality rates may be as high as 40%, with the precise burden of severe neurological disability being higher in underdeveloped countries [5,11,12].
Low birth weight, preterm deliveries, malpresentation, mode of delivery, meconium-stained amniotic fluid, chorioamnionitis, and prolonged rupture of membranes were shown to be risk factors for poor outcomes for neonates who were born asphyxiated [13,14].Numerous clinical, pathological, biochemical, and metabolic alterations result from perinatal hypoxia, which increases the risk of death [15].The survival status of neonates with birth asphyxia was low, particularly in the primary care facilities in Ethiopia [16].According to research, stage III birth asphyxia is the leading cause of perinatal and neonatal death [17].
Despite improvements in the management of perinatal care, ANC follow-up strategies, and the accessibility of NICU (neonatal intensive care unit) care, perinatal asphyxia continues to account for the majority of neonatal intensive care unit (NICU) admissions and is one of the leading causes of neonatal deaths in underdeveloped countries [18][19][20].A study conducted in Vietnam identified that prematurity, asphyxia, and designated congenital malformations are underestimated because most of the research did not consider the comorbidities [21,22].A retrospective cohort study establishes perinatal asphyxia as a substantial factor in illness and mortality among neonates, and the incidence of mortality from PNA was not adequately reported, but it is high according to some studies, particularly in developing countries [13,23,24].In Ethiopian healthcare institutions, midwives' and nurses' retention of their knowledge of neonatal resuscitation practices could contribute to an increase in the number of neonates dying from asphyxia [25].
The neonatal mortality rate is still high and linked to some predictors, and the leading cause of death was perinatal asphyxia [26].Neonatal mortality is unexpectedly growing in Ethiopia, while it has been stagnating in many low-income countries at 28-33 per 1000 live births [27].Perinatal asphyxia was the main reason for admission (55.3%), which increases the economic impact and the increased rate of neonate death, according to a facilitybased cross-sectional study, but lacks the incidences of death and related factors among neonates admitted with perinatal asphyxia [28,29].However, perinatal asphyxia is acknowledged as the greatest contributor to neonatal mortality, but previous research in this field has focused on the prevalence of PNA and missed some important variables.The goal of this study is to determine the incidence of PNA-related mortality and identify the factors that predict mortality in neonates with perinatal asphyxia.

Study design and period
Institution-based retrospective follow-up study design was employed from October 3 to November 2, 2022.

Study setting
The study was conducted at five tertiary hospitals in the west Oromia region.Those hospitals were Mettu Karl Referral Hospital (MKRH), Jimma University Medical Centre (JUMC), Wollega University Referral Hospital (WURH), Nekemte Specialized Hospital (NSH), and Ambo University Referral Hospital (AURH).Each hospital is far from Addis Ababa, the Ethiopian capital, at 600 km, 420 km, 325 km, and 125 km, respectively.All of them provide tertiary NICU service when there are 10 nurses (4 neonatal nurses), 1 pediatrician at MKRH, 11 nurses (2 neonatal nurses), and 2 pediatricians at NSRH, 22 nurses (8 neonatal nurses), 8 pediatricians, and 35 residents at JUMC.Eleven nurses (6 neonatal nurses), 5 pediatricians, and 13 residents at WURH, and 14 nurses (0 neonatal nurses), 7 pediatricians, and 12 residents at AURH are actively working in each hospital.Average annual ANC services were 6100 in WURH, 7500, 7900, 5700, and 6300 in NSRH, JUMC, AURH, and MKRH, respectively.The client transfer referral system with a neonatal referral form is used in these hospitals through three tiers of health facilities.Annually, in all hospitals, about 5,200 neonates are admitted, and in the previous four years, in all hospitals, 3,934 perinatally asphyxiated neonates were admitted.

Source population
The source populations were all neonates diagnosed with perinatal asphyxia and admitted to the NICU wards of the tertiary hospitals of the west Oromia region.

Study population
All neonates diagnosed with perinatal asphyxia and admitted to the NICU of the Western Oromia tertiary hospitals from January 1, 2018, to December 31, 2021, were included in the study population.

Inclusion and exclusion criteria Inclusion criteria
All neonates in the study area had PNA diagnoses and were admitted in the NICU between January 1, 2018, and December 31, 2021.

Exclusion criteria
Records of neonates whose admission date, discharge date, and outcome were not recorded on the chart and home delivery were excluded from the study.

Sample size determination
For the first objective, the sample size was calculated by using the single population proportion formula with consideration of the following statistical assumptions based on previous study [30]: a 7.85 proportion of deaths (p = 0.0785 and q = 0.9215), a 95% confidence level, and a 5% margin of error.
The sample size was also calculated by using the second objective to check the adequacy of the sample size used for a survival sample size calculation power approach using STATA 14 software with Cox proportional assumptions.The four predictive variables included were cord prolapse, pregnancy-induced hypertension, maternal iron deficiency anemia, and convulsions from the previous study conducted in southern Ethiopia [30].By using the Schoenfeld formula (Schoenfeld DA, sample-size formula for the proportional-hazards regression model, biometrics 1983;39 (2): 499-503.)[31].
Where E = number of events.N = number of sample sizes.P1 = proportion of the event among exposed.Prev.= Probability of an event.P2 = proportion of event among non-exposed.
From those predictors, the sample size obtained from maternal iron deficiency anemia was considered the final sample size of the study because it gives a maximum sample size.Using the probability of events (7.85), and a crude hazard ratio of (2.28), proportional withdrawals (10%), 95% level of confidence & and power of (80%); the final sample size becomes 655 (Table 1).

Sampling technique and procedure
A sample of 655 charts of neonates with PNA was selected from 3,934 neonates' charts admitted in the same case, and the subsequent unique chart numbers from the registration file were extracted using a computer-generated random number approach.Estimated previous four-year PNA neonates admitted to these hospitals were: Mettu Karl Referral Hospital (MKRH) (762), Jimma University Medical Centre (JUMC) (898), Wollega University Referral Hospital (WURH) (558), Nekemte Specialized Hospital (NSH) (970), and Ambo University Referral Hospital (AURH) (746) between January 1, 2018, and December 31, 2021.
Since the final sample size was (655), proportional allocation was done for each hospital.The sampling frame was prepared by collecting the identification numbers of PNA patients from the registration book.After identifying the patients who fulfilled the inclusion criteria, study subjects were selected by a simple random sampling technique using computer-generated random numbers (Fig. 1).Defaulted PNA cases that are signed (by parents on behalf of their child) against treatment to leave treatment before the cure are completed.
Stages of HIE were determined based on saranat's classifications of the clinical diagnosis made by a health care provider (stage I, stage II, and stage III) [2].

Data collection tool and procedure
A data extraction format was adapted from the peerreviewed articles [16,32,33]

Data quality control assurance
To ensure the quality of the data, the data extraction tool was checked for the existence of variables in the registration format on the patient's chart via a preliminary chart follow-up review of 5% (33 charts) of the sample, and the tool validity was checked by experts (a pediatrician and a neonatal nurse).Data collectors were nurses who had experience working in the NICU ward.In addition, data collectors were trained for one day at each hospital before data collection started.The data retrieval process was closely monitored by the principal investigator and/or two supervisors through a well-prepared and restricted data collection tool by Kobo Toolbox software.

Data processing and analysis
Data inconsistency, coding errors, completeness, clarity, missing values, and data cleaning were checked before being exported.The data were exported to the STATA 14 version for further analysis.Descriptive statistics of numeric variables were carried out, and categorical variables were presented using frequency, tables, and percentages.The incidence rate of mortality was also calculated by dividing the total number of new occurrences of death among PNA by the total number of person-days of follow-up.The Kaplan-Meier curve was used to estimate survival time.The log-rank test was used to compare survival experiences between independent groups.The proportional hazard assumption was checked both graphically and using a Schoenfeld residual test, which assesses the relationship between the scaled Schoenfeld residuals and time.Multicollinearity was checked between independent variables by a variance inflation factor.The frailty model was taken into account to handle variations among hospitals (Table 2).
The log-likelihood and Akaike Information Criteria (AIC) were applied to select the best-fitted model, and a model with a minimum AIC value was considered the best-fitted model.The goodness of fit of the model was assessed by Cox-Snell residuals and the Nelson-Aalen cumulative hazard function graph.A variable with a P < 0.2 in the bivariable analyses was considered a potential candidate variable for the multivariable analysis.A variable with a p-value ≤ 0.05 in the multivariable model was declared a significant predictor of the outcome of interest.Finally, the crude and adjusted hazard ratio (HR) with a 95% confidence interval (CI) was computed.

Sociodemographic characteristics
Almost half of the mothers (326, or 52.925%) came from rural areas and came from an area with a distance of > 10 km (149, or 24.19%), and (464, or 75.32%) were aged between 21 and 34 years, and the median age of mothers was 24 years (IQR: 7).About (358, or 58.12%) of the neonates were male, and (270, or 43.83%) of them were admitted after 60 min of age.Most of the mothers (476, or 77.27%) were given birth in hospital (Table 3).

Maternal and obstetric-related factors
Half of the mothers (314, or 50.97%) attended the ANC three times or more during their current pregnancy.More than half of the mothers (365, or 59.25%) were multi-gravida, and about two-thirds (474, or 76.95%) of them had no history of pregnancy-induced hypertension.About (357, or 57.95%) of mothers had a history of non-prolonged rupture of membranes before delivery; more than half of them (366, or 59.42%) had no history of obstructed labor, and the majority of mothers (436, or 70.78%) were in labor for less than 18 h.Approximately half of the mothers, 303 (49.19%), had spontaneous vaginal deliveries during the night (328, or 53.25%).During the current pregnancy, the majority of the mothers (477, or 77.44%), (569, or 92.37%), (494, or 80.19%), and (604, or 98.05%) had no history of anemia, DM, abortion, or HIV infection, respectively (Table 4).

Medical care-related factors
Two-thirds (387 or 62.82%) of the neonates' orders were revised every 12 h.About half (328 or 53.25%) and (329 or 53.41%) of the asphyxiated neonates were treated with CPAP and initiated for enteral feeding after 48 h, respectively.The majority of the participants (471 or 76.46%) were given antibiotics, and more than one-third (197 or 31.98%) of the neonates were given anti-epileptic drugs (Table 6).

Overall neonatal outcome
From the total 28 days followed, 616 neonates (202, or 32.79%) developed the outcome variable (death) after being diagnosed with perinatal asphyxia and admitted to the NICU (Fig. 2).

Accessing proportional hazard assumption
Accessing the proportional hazard assumption was tested by using Schoenfeld residual proportional hazard assumption, and all the covariates fulfilled the assumption, and the overall global model result is described (Table 7).

Model comparison
Both semi-parametric and parametric proportional hazard models were fitted to determine the survival time to death and identify its determinants among neonates with perinatal asphyxia after the proportional hazard assumption was verified.The model comparison was graphically  evaluated by Cox-Snell residual.Furthermore, model comparison was statistically examined using information criteria (AIC) and log-likelihood to determine which model was the best fit for the given data set.As a result, the Weibull regression model with an AIC of 895.47 was the model that suited the data (Table 8; Fig. 3).
With a shape parameter (P) greater than zero, which indicates the hazard exponentially increases over time.The Weibull distribution was used to match the baseline hazard (the effect of time when all categorical variables are in a reference category) of experiencing an event or death (P = 1.71 95% CI: 1.53-1.90).The shared frailty was checked according to estimates, demonstrating that there is no difference between the hospitals in the distribution of unmeasured variables.

Incidence of mortality
Of the total 616 neonates that were admitted with perinatal asphyxia (PNA), 202 (32.79%) (95% CI: 29.18-36.61) of them experienced the event of interest (death).The total neonate days observed over the full follow-up period were 5198 person-days, comprised of 883 from WURH, 1400 from NSH, 912 from MKSH, 1044 from JUMC, and 959 from AURH, respectively, with minimum and maximum follow-up times of 1 and 27 days.The 95% confidence interval for the median follow-up time was 14-20 days, or 15 days.Throughout the whole period of follow-up, the overall incidence of mortality was

Overall survival function
The overall Kaplan-Meier failure function demonstrated that neonates with PNA had a higher risk of dying during the follow-up period.On the first day after admission, there was a 1.78% chance of dying.The cumulative likelihood of dying at the end of 5, 10, and 15 days was 0.59 at five days, 0.57 at 10 days, and 0.50 at 15 days, respectively, and the least probability of death was observed at the end of 15 days of follow-up time (Fig. 4).

Predictors of mortality among neonates with perinatal asphyxia
With Kaplan-Meier and the log-rank test, the survival distribution for several categories of variables was tested for equality.In general, it was evident from the pattern that one group's survivorship function was positioned above another group and that the lower curve had a higher chance of death than the upper curve (Figs. 5 and  6).The upper groups, on the other hand, as depicted by the Kaplan-Meier survival curve, had a high survival rate.Additionally, the log-rank test determined whether the observed difference was a statistical difference that could be seen on the KM graph.
In the current study, neonates whose families lived more than 10 km from a health facility had shorter Fig. 3 Summary of model comparison using Cox-Snell residual graph for neonates with perinatal asphyxia, West Oromia Tertiary Hospitals Ethiopia 2022 A residual is the difference between an observed data point and a predicted or fitted value and as the following graph indicates the Weibull regression fits than other models with hazards greatly followed the predicted hazard at 45' survival times than their counterparts.Moreover, neonates whose mothers had a history of prolonged membrane rupture during their current pregnancy had a lower survival probability than those whose mothers had not.
This study also revealed that neonates who had low APGAR scores (< 3) at the fifth minute of birth were more likely to die when compared with those with scores greater than three.Furthermore, asphyxiated neonates who were treated with CPAP had a better probability of survival compared with neonates that were treated with direct oxygen.Moreover, asphyxiated neonates who had been diagnosed with stage III hypoxic-ischemic encephalopathy had a lower survival probability compared with those diagnosed with stage I and stage II HIE.These differences were statistically significant with a p-value of < 0.05.

Bivariable and multivariable Weibull regression analysis
In the bivariable Weibull proportional hazard model, onset of labor, stage of HIE, obstructed labor, history of cry at birth, frequency of revision of order, type of oxygen used, antibiotics, antiepileptics, thrombocytopenia, resuscitation, preeclampsia, birth type, duration of labor, congenital anomaly, history of prolonged ROM, neonatal age at admission, APGAR score at the fifth minute, and distance from the health facility were associated with the incidence of mortality among neonates admitted with PNA.However, in the multivariable analysis, distance from HF, APGAR score at the fifth minute, stage of HIE, type of oxygen used, and prolonged rupture of the membrane were significant predictors of mortality among neonates admitted with PNA (Table 9).
Keeping other variables constant, the hazard of mortality was about 2 fold higher in asphyxiated neonates with stage III HIE (AHR: 2.02; 95% CI: 1.18-3.47)than in stage I and stage II HIE.Neonates whose mothers had a history of prolonged membrane rupture had a 1.41-fold higher hazard of death than those whose mothers had not (AHR: 1.41; 95% CI: 1.02-1.94).Neonates with an APGAR score of less than 3 at the fifth minute were about 3 times more likely to die (AHR: 2.63; 95% CI: 1.03-6.73)than neonates with an APGAR score of greater than 3 at the fifth minute by keeping out other variables constant.Moreover, neonates who traveled more than 10 km from HF had a 2-fold hazard of death than those who traveled less than 10 km (AHR: 2.25; 95% CI: 1.60-3.17)while other variables remained constant.Furthermore, when all other variables were held constant, neonates on direct oxygen had about a 2-fold higher hazard of death (AHR: 1.83; 95% CI: 1.35-2.48)than asphyxiated neonates on CPAP.

Discussion
The primary objective of this study was to estimate the incidence rate of death and assess its predictors among asphyxiated neonates admitted to the NICU.From the total of 616 asphyxiated babies enrolled in this study, 202 (32.79%) (95% CI: 29.18-36.61)neonates experienced the event of interest (died).Throughout the full follow-up period, the overall incidence of death was 38.86 per 1000 neonates (95% CI: 33.85-44.60).
This finding is lower than that of the study conducted in northwest Ethiopia, which was (53.49 per 1000) [33].This marked difference might be due to the study area, as the current study is multi-intuitional, while the northwest Ethiopia study was conducted monocentrically as a single hospital incidence represents a limited area of the population.Another possible justification could be due to the difference in overall included data, as this study followed data for four years, while the northwest Ethiopia study followed data three years.Meanwhile, different from the period of data incorporation, neonatal care is improving in the current period, as this study used more recent data than the northwest Ethiopia study [33].
According to this study's findings the incidence of mortality is higher than the study conducted in Brazil 1.71 deaths per 1,000 [34].This great difference might be due to the difference in study methods, as this study used a prospective cohort while our study used a retrospective follow-up because the prospective follow-up study follows the actual patients up to the occurrence of the event.Another reason could be due to differences in geographical area, health facilities relative to the total population, the skill of neonatal resuscitation among midwives and nurses, the structure of the NICU, and the advanced technology used in neonatal care [34].
The proportion of deaths in this finding is greater than in the studies conducted in southern Ethiopia (7.85%) [30], Addis Ababa (24.09%) [32], Nigeria (18%) [23], and Tanzania (23%) [7].The difference between the current study findings and the study of southern Ethiopia might be due to the difference in study design, as this study used a retrospective follow-up and the study of southern Ethiopia was a prospective follow-up.In addition, it might be due to differences in health care provider level; in Addis Ababa, even if it is not adequate, there are neonatologists, but not in our study area, which may make a difference in neonatal care [30].
The variability of results obtained in this study and a study conducted in Addis Ababa might be due to the discrepancy in sample size, the discrepancy in health facilities and health professionals in these study areas, and political instability in the current study area.In other words, the study conducted in Addis Ababa did not consider the comorbidities like sepsis, respiratory distress, and medical care that have been carried out for neonates, on which the current study was focused [32].
To justify the difference, the study conducted in Nigeria used clinical laboratory investigations as inclusion criteria, but in this study, clinical diagnosis was the only criterion for the diagnosis of PNA, and there might be differences in the study population and treatment modalities in two different countries [23].Meanwhile, the reason why the study conducted in Tanzania was lower might be because the study conducted in Tanzania was monocenteric and followed two years of PNA data, while this study followed four years of multi-center data.In addition to the aforementioned factors, the variability could be explained by the fact that the Tanzanian study excluded neonates who died within the first 30 min of admission, whereas the current study included all types of deaths [7].
The proportion of deaths in the current study is lower than the studies conducted in northern Ethiopia (37.5%) [35], and eastern Ethiopia (55.3%) [28].The difference could be attributed to the study period, as the northern Ethiopia study only looked at two years of data, whereas the current study looked at four years.Another reason might be that the northern Ethiopia study included all NICU-admitted neonates, whereas this study evaluated neonates with PNA.Meanwhile, there could be differences in sample size and study design [35].
The discrepancy between this study and the study conducted in eastern Ethiopia might be due to the monocenteric study and the difference in sample size.Meanwhile, in the Eastern Ethiopian study, all neonates admitted to the NICU were included, whereas in the current study, only perinatally asphyxiated neonates were evaluated [28].
In this study, severe asphyxia or less than 3 APGAR scores at the fifth minute is an independent predictor and increases the hazards of death 3 times among neonates with PNA.This is supported by the study conducted in Addis Ababa [36], Nigeria [37], Tanzania [7], and Cameroon [38].The reason for this finding might be perinatal asphyxia, which is a lack of blood flow or gas exchange to or from the fetus in the period immediately before, during, or after the birth process.The APGAR score is a numeric assessment of the neonatal condition at birth and can be affected by numerous factors, including gestational age, medication taken by the mother, resuscitation efforts, and cardiac, respiratory, or neurological abnormalities in the baby.If the neonate is severely asphyxiated at the golden times of the first and fifth minutes, the tissues and vital organs (muscle, liver, heart, and ultimately the brain) will develop an oxygen deficit.Anaerobic glycolysis and lactic acidosis will result and also lead to death [1,39].Improving the resuscitation skills of NICU and labor and delivery health care providers, prompt referral to the hospital, and strict follow-up during labor and delivery are needed to alleviate the occurrence of low APGAR scores in the perinatal period.
In the current study, using direct oxygen rather than CPAP is another predictor of mortality and increases the hazard of death almost 2 times.This finding is supported by different studies conducted in Ugandan [40], and China [41].Justification for this finding could be that CPAP is a simple and effective respiratory support modality used to support neonates with respiratory failure by providing constant pressure maintained throughout the respiratory cycle but with no additional inspiratory pressure support.CPAP keeps alveoli open, improves oxygenation by reducing the amount of blood shunted through atelectasis while the infant breathes spontaneously, and reduces the hazards of mortality from asphyxia [42].As this predictive variable indicates, healthcare providers should be aware of the use of CPAP to minimize the incidence of death among PNA neonates.
A neonate who develops stage III HIE is an independent factor that predicts mortality and increases the hazards of death by two.This finding is supported by the study conducted in northwest Ethiopia [33], India [18,[43][44][45], Nigeria [23,46], and Cameroon [38].It could be due to the fact that, hypoxic ischemic brain injury is an ongoing process composed of several different phases.When the primary critical energy failure occurs, an uncontrolled release of excitatory neurotransmitters begins, starting the ischemic cascade that damages neuronal cells (both at the cytoplasmic and mitochondrial levels), disrupts the brain-blood barrier, and activates an important inflammatory response.These result in failure of oxidative metabolism, cytotoxic edema, and the accumulation of excitotoxins, which may lead the neonates to death [47].Because these neonates had multi-organ failure, we need to pay special attention to them while providing cares to reduce mortality.
In this study, PNA neonates of mothers who had a history of PROM in their current pregnancy are another factor that predicts mortality among PNA and increases the hazards of death by about 1.41 times.This finding is supported by the studies conducted in India [48], and Cameroon [38].The reasons could be due to the complications of prolonged rupture of the membrane (amnionitis and endometritis), which both put the fetus at a high risk of developing an overwhelming infection in the bloodstream (sepsis), which can exacerbate hypoxia and lead to death [2].Minimizing the time of labor after rupture of the membrane could ameliorate this condition, which is a predictor of death among PNA.
According to the findings of this study, distance from a health facility is a significant predictor of death from PNA and doubles the hazards of death.This finding might be due to a delay in seeking health care, which can prolong the initiation of care for asphyxiated neonates, increase the severity of the problem, and lead to death.It is also supported by the study conducted in Uganda, as delay in seeking health care is directly associated with the distance from a health facility [49].Providing improved community-based newborn care and providing health education for pregnant mothers on how to plan where to give birth and when to seek care can change the effect of distance from a health facility on PNA neonates' ability to survive.

Limitations of the study
The primary weakness of the study could be incomplete medical records because it was based on secondary data.In this regard, medically recorded data did not provide information on certain important factors, such as mothers' educational attainment, nutritional status, and monthly income.Furthermore, because the study was conducted in a hospital, neonatal mortality may have been underestimated, as home deliveries that ended in neonatal deaths were not properly reported.

Conclusion
The study found that neonatal mortality among perinatal asphyxia patients remains high.Besides, the overall incidence rate of mortality was 38.86 per 1000 neonates.The type of oxygen used, residence distance, history of prolonged ROM, and stages of HIE were significant predictors of mortality.Therefore, pregnant mothers should get information on when to seek care following their pregnancy.Providing appropriate follow-up for PNA neonates with stage III HIE, minimizing the time of labor after rupture of membrane, and evaluating the treatment options regarding the type of oxygen for PNA neonates are also mandatory to reduce mortality from perinatal asphyxia.Moreover, future researchers should conduct an interventional study to address further significant predictors of mortality in asphyxiated neonates.

Fig. 1
Fig. 1 Schematic presentation of sampling procedure for incidence and predictors of mortality among neonates with perinatal asphyxia admitted at West Oromia tertiary hospitals from January 1, 2018 to December 31, 2021

Fig. 5 Fig. 4
Fig. 5 Kaplan-Meier inter-group comparison of survival function among neonates with PNA at West Oromia Tertiary Hospitals Ethiopia, 2022

Fig. 6
Fig. 6 Kaplan-Meier inter-group comparison of survival function among neonates with PNA at West Oromia Tertiary Hospitals Ethiopia, 2022

Table 1
Summary of sample size calculation to access incidence and predictors of mortality among neonates with PNA admitted at West Oromia Tertiary hospitals Ethiopia, 2022 Censored neonates with PNA who did not develop the outcome of interest (death) until the end of the followup period or lost to follow-up, recovered from illness, discharged against medical advice, or transferred out to other health institutions without knowing the outcome.

Table 2
The frailty test procedure to handle variations between hospitals, among neonates with PNA, at West Oromia Tertiary Hospitals Ethiopia, 2022

Table 3
Socio-demographic characteristics of the mother and neonate with perinatal asphyxia at West Oromia Tertiary Hospitals, Ethiopia, 2022

Table 4
Maternal and obstetric related characteristics of newborn with perinatal asphyxia at, West Oromia Tertiary Hospitals Ethiopia, 2022

Table 5
Clinical characteristics of the new-born with perinatal asphyxia at, West Oromia Tertiary Hospitals Ethiopia, 2022

Table 7
Global test for proportional hazard assumption using schoenfeld residuals for neonates with perinatal asphyxia, West Oromia Tertiary Hospitals Ethiopia 2022 All authors have read and approved the tables

Table 9
Bivariate and multivariate Weibull regression analysis for predictors of death among neonates with PNA, at West Oromia Tertiary Hospitals Ethiopia, 2022