Predicting COVID-19 outcomes from clinical and laboratory parameters in an intensive care facility during the second wave of the pandemic in South Africa

Background: The second wave of coronavirus disease 2019 (COVID-19) in South Africa was caused by the Beta variant of severe acute respiratory syndrome coronavirurus-2. This study aimed to explore clinical and biochemical parameters that could predict outcome in patients with COVID-19. Methods: A prospective study was conducted between 5 November 2020 and 30 April 2021 among patients with confirmed COVID-19 admitted to the intensive care unit (ICU) of a tertiary hospital. The Cox proportional hazards model in Stata 16 was used to assess risk factors associated with survival or death. Factors with P<0.05 were considered significant. Results: Patients who died were found to have significantly lower median pH (P<0.001), higher median arterial partial pressure of carbon dioxide (P<0.001), higher D-dimer levels (P=0.001), higher troponin T levels (P=0.001), higher N-terminal-prohormone B-type natriuretic peptide levels (P=0.007) and higher C-reactive protein levels (P=0.010) compared with patients who survived. Increased standard bicarbonate (HCO3std) was associated with lower risk of death (hazard ratio 0.96, 95% confidence interval 0.93–0.99). Conclusions: The mortality of patients with COVID-19 admitted to the ICU was associated with elevated D-dimer and a low HCO3std level. Large studies are warranted to increase the identification of patients at risk of poor prognosis, and to improve the clinical approach

Latent class analysis: an innovative approach for identification of clinical and laboratory markers of disease severity among COVID-19 patients admitted to the intensive care unit

Objective: The aim of this study was to identify clinical and laboratory phenotype distribution patterns and their usefulness as prognostic markers in COVID-19 patients admitted to the intensive care unit (ICU) at Tygerberg Hospital, Cape Town. Methods and results: A latent class analysis (LCA) model was applied in a prospective, observational cohort study. Data from 343 COVID-19 patients were analysed. Two distinct phenotypes (1 and 2) were identified, comprising 68.46% and 31.54% of patients, respectively. The phenotype 2 patients were characterized by increased coagulopathy markers (D-dimer, median value 1.73 ng/L vs 0.94 ng/L; p < 0.001), end-organ dysfunction (creatinine, median value 79 µmol/L vs 69.5 µmol/L; p < 0.003), under-perfusion markers (lactate, median value 1.60 mmol/L vs 1.20 mmol/L; p < 0.001), abnormal cardiac function markers (median N‐terminal pro‐brain natriuretic peptide (NT-proBNP) 314 pg/ml vs 63.5 pg/ml; p < 0.001 and median high‐sensitivity cardiac troponin (Hs-TropT) 39 ng/L vs 12 ng/L; p < 0.001), and acute inflammatory syndrome (median neutrophil-to-lymphocyte ratio 15.08 vs 8.68; p < 0.001 and median monocyte value 0.68 × 109/L vs 0.45 × 109/L; p < 0.001). Conclusion: The identification of COVID-19 phenotypes and sub-phenotypes in ICU patients could help as a prognostic marker in the day-to-day management of COVID-19 patients admitted to the ICU

Comparison of patients with severe COVID-19 admitted to an intensive care unit in South Africa during the first and second wave of the COVID-19 pandemic

BACKGROUND: The second wave of coronavirus disease 2019 (COVID‑19), dominated by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Beta variant, has been reported to be associated with increased severity in South Africa (SA). OBJECTIVES: To describe and compare clinical characteristics, management and outcomes of COVID‑19 patients admitted to an intensive care unit (ICU) in SA during the first and second waves. METHODS: In a prospective, single-centre, descriptive study, we compared all patients with severe COVID‑19 admitted to ICU during the first and second waves. The primary outcomes assessed were ICU mortality and ICU length of stay (LOS). RESULTS: In 490 patients with comparable ages and comorbidities, no difference in mortality was demonstrated during the second compared with the first wave (65.9% v. 62.5%, p=0.57). ICU LOS was longer in the second wave (10 v. 6 days, p<0.001). More female admissions (67.1% v. 44.6%, p<0.001) and a greater proportion of patients were managed with invasive mechanical ventilation than with non-invasive respiratory support (39.0% v. 14%, p<0.001) in the second wave. CONCLUSIONS: While clinical characteristics were comparable between the two waves, a higher proportion of patients was invasively ventilated and ICU stay was longer in the second. ICU mortality was unchanged

Prognostic value of biochemical parameters among severe COVID-19 patients admitted to an intensive care unit of a tertiary hospital in South Africa

Background: Data on biochemical markers and their association with mortality rates in patients with severe coronavirus disease 2019 (COVID-19) admitted to intensive care units (ICUs) in sub-Saharan Africa are scarce. An evaluation of baseline routine biochemical parameters was performed in COVID-19 patients admitted to the ICU, in order to identify prognostic biomarkers. Methods: Demographic, clinical, and laboratory data were collected prospectively from patients with PCR-confirmed COVID-19 admitted to the adult ICU of a tertiary hospital in Cape Town, South Africa, between October 2020 and February 2021. Robust Poisson regression methods and the receiver operating characteristic (ROC) curve were used to explore the association of biochemical parameters with severity and mortality. Results: A total of 82 patients (median age 53.8 years, interquartile range 46.4–59.7 years) were enrolled, of whom 55 (67%) were female and 27 (33%) were male. The median duration of ICU stay was 10 days (interquartile range 5–14 days); 54/82 patients died (66% case fatality rate). Baseline lactate dehydrogenase (LDH) (adjusted relative risk 1.002, 95% confidence interval 1.0004–1.004; P = 0.016) and N-terminal pro B-type natriuretic peptide (NT-proBNP) (adjusted relative risk 1.0004, 95% confidence interval 1.0001–1.0007; P = 0.014) were both found to be independent risk factors of a poor prognosis, with optimal cut-off values of 449.5 U/l (sensitivity 100%, specificity 43%) and 551 pg/ml (sensitivity 49%, specificity 86%), respectively. Conclusions: LDH and NT-proBNP appear to be promising predictors of a poor prognosis in COVID-19 patients in the ICU. Studies with a larger sample size are required to confirm the validity of this combination of biomarkers

Haematological predictors of poor outcome among COVID-19 patients admitted to an intensive care unit of a tertiary hospital in South Africa

BACKGROUND: Studies from Asia, Europe and the USA indicate that widely available haematological parameters could be used to determine the clinical severity of Coronavirus disease 2019 (COVID-19) and predict management outcome. There is limited data from Africa on their usefulness in patients admitted to Intensive Care Units (ICUs). We performed an evaluation of baseline haematological parameters as prognostic biomarkers in ICU COVID-19 patients. METHODS: Demographic, clinical and laboratory data were collected prospectively on patients with confirmed COVID-19, admitted to the adult ICU in a tertiary hospital in Cape Town, South Africa, between March 2020 and February 2021. Robust Poisson regression methods and receiver operating characteristic (ROC) curves were used to explore the association of haematological parameters with COVID-19 severity and mortality. RESULTS: A total of 490 patients (median age 54.1 years) were included, of whom 237 (48%) were female. The median duration of ICU stay was 6 days and 309/490 (63%) patients died. Raised neutrophil count and neutrophil/lymphocyte ratio (NLR) were associated with worse outcome. Independent risk factors associated with mortality were age (ARR 1.01, 95%CI 1.0–1.02; p = 0.002); female sex (ARR 1.23, 95%CI 1.05–1.42; p = 0.008) and D-dimer levels (ARR 1.01, 95%CI 1.002–1.03; p = 0.016). CONCLUSIONS: Our study showed that raised neutrophil count, NLR and D-dimer at the time of ICU admission were associated with higher mortality. Contrary to what has previously been reported, our study revealed females admitted to the ICU had a higher risk of mortality

Cardiopulmonary disease as sequelae of long-term COVID-19: Current perspectives and challenges

COVID-19 infection primarily targets the lungs, which in severe cases progresses to cytokine storm, acute respiratory distress syndrome, multiorgan dysfunction, and shock. Survivors are now presenting evidence of cardiopulmonary sequelae such as persistent right ventricular dysfunction, chronic thrombosis, lung fibrosis, and pulmonary hypertension. This review will summarize the current knowledge on long-term cardiopulmonary sequelae of COVID-19 and provide a framework for approaching the diagnosis and management of these entities. We will also identify research priorities to address areas of uncertainty and improve the quality of care provided to these patients

Global, regional, and national burden of disorders affecting the nervous system, 1990–2021: a systematic analysis for the Global Burden of Disease Study 2021

BackgroundDisorders affecting the nervous system are diverse and include neurodevelopmental disorders, late-life neurodegeneration, and newly emergent conditions, such as cognitive impairment following COVID-19. Previous publications from the Global Burden of Disease, Injuries, and Risk Factor Study estimated the burden of 15 neurological conditions in 2015 and 2016, but these analyses did not include neurodevelopmental disorders, as defined by the International Classification of Diseases (ICD)-11, or a subset of cases of congenital, neonatal, and infectious conditions that cause neurological damage. Here, we estimate nervous system health loss caused by 37 unique conditions and their associated risk factors globally, regionally, and nationally from 1990 to 2021.MethodsWe estimated mortality, prevalence, years lived with disability (YLDs), years of life lost (YLLs), and disability-adjusted life-years (DALYs), with corresponding 95% uncertainty intervals (UIs), by age and sex in 204 countries and territories, from 1990 to 2021. We included morbidity and deaths due to neurological conditions, for which health loss is directly due to damage to the CNS or peripheral nervous system. We also isolated neurological health loss from conditions for which nervous system morbidity is a consequence, but not the primary feature, including a subset of congenital conditions (ie, chromosomal anomalies and congenital birth defects), neonatal conditions (ie, jaundice, preterm birth, and sepsis), infectious diseases (ie, COVID-19, cystic echinococcosis, malaria, syphilis, and Zika virus disease), and diabetic neuropathy. By conducting a sequela-level analysis of the health outcomes for these conditions, only cases where nervous system damage occurred were included, and YLDs were recalculated to isolate the non-fatal burden directly attributable to nervous system health loss. A comorbidity correction was used to calculate total prevalence of all conditions that affect the nervous system combined.FindingsGlobally, the 37 conditions affecting the nervous system were collectively ranked as the leading group cause of DALYs in 2021 (443 million, 95% UI 378–521), affecting 3·40 billion (3·20–3·62) individuals (43·1%, 40·5–45·9 of the global population); global DALY counts attributed to these conditions increased by 18·2% (8·7–26·7) between 1990 and 2021. Age-standardised rates of deaths per 100 000 people attributed to these conditions decreased from 1990 to 2021 by 33·6% (27·6–38·8), and age-standardised rates of DALYs attributed to these conditions decreased by 27·0% (21·5–32·4). Age-standardised prevalence was almost stable, with a change of 1·5% (0·7–2·4). The ten conditions with the highest age-standardised DALYs in 2021 were stroke, neonatal encephalopathy, migraine, Alzheimer's disease and other dementias, diabetic neuropathy, meningitis, epilepsy, neurological complications due to preterm birth, autism spectrum disorder, and nervous system cancer.InterpretationAs the leading cause of overall disease burden in the world, with increasing global DALY counts, effective prevention, treatment, and rehabilitation strategies for disorders affecting the nervous system are needed

Global, regional, and national burden of diabetes from 1990 to 2021, with projections of prevalence to 2050: a systematic analysis for the Global Burden of Disease Study 2021

Background: Diabetes is one of the leading causes of death and disability worldwide, and affects people regardless of country, age group, or sex. Using the most recent evidentiary and analytical framework from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD), we produced location-specific, age-specific, and sex-specific estimates of diabetes prevalence and burden from 1990 to 2021, the proportion of type 1 and type 2 diabetes in 2021, the proportion of the type 2 diabetes burden attributable to selected risk factors, and projections of diabetes prevalence through 2050. Methods: Estimates of diabetes prevalence and burden were computed in 204 countries and territories, across 25 age groups, for males and females separately and combined; these estimates comprised lost years of healthy life, measured in disability-adjusted life-years (DALYs; defined as the sum of years of life lost [YLLs] and years lived with disability [YLDs]). We used the Cause of Death Ensemble model (CODEm) approach to estimate deaths due to diabetes, incorporating 25 666 location-years of data from vital registration and verbal autopsy reports in separate total (including both type 1 and type 2 diabetes) and type-specific models. Other forms of diabetes, including gestational and monogenic diabetes, were not explicitly modelled. Total and type 1 diabetes prevalence was estimated by use of a Bayesian meta-regression modelling tool, DisMod-MR 2.1, to analyse 1527 location-years of data from the scientific literature, survey microdata, and insurance claims; type 2 diabetes estimates were computed by subtracting type 1 diabetes from total estimates. Mortality and prevalence estimates, along with standard life expectancy and disability weights, were used to calculate YLLs, YLDs, and DALYs. When appropriate, we extrapolated estimates to a hypothetical population with a standardised age structure to allow comparison in populations with different age structures. We used the comparative risk assessment framework to estimate the risk-attributable type 2 diabetes burden for 16 risk factors falling under risk categories including environmental and occupational factors, tobacco use, high alcohol use, high body-mass index (BMI), dietary factors, and low physical activity. Using a regression framework, we forecast type 1 and type 2 diabetes prevalence through 2050 with Socio-demographic Index (SDI) and high BMI as predictors, respectively. Findings: In 2021, there were 529 million (95% uncertainty interval [UI] 500-564) people living with diabetes worldwide, and the global age-standardised total diabetes prevalence was 6·1% (5·8-6·5). At the super-region level, the highest age-standardised rates were observed in north Africa and the Middle East (9·3% [8·7-9·9]) and, at the regional level, in Oceania (12·3% [11·5-13·0]). Nationally, Qatar had the world's highest age-specific prevalence of diabetes, at 76·1% (73·1-79·5) in individuals aged 75-79 years. Total diabetes prevalence-especially among older adults-primarily reflects type 2 diabetes, which in 2021 accounted for 96·0% (95·1-96·8) of diabetes cases and 95·4% (94·9-95·9) of diabetes DALYs worldwide. In 2021, 52·2% (25·5-71·8) of global type 2 diabetes DALYs were attributable to high BMI. The contribution of high BMI to type 2 diabetes DALYs rose by 24·3% (18·5-30·4) worldwide between 1990 and 2021. By 2050, more than 1·31 billion (1·22-1·39) people are projected to have diabetes, with expected age-standardised total diabetes prevalence rates greater than 10% in two super-regions: 16·8% (16·1-17·6) in north Africa and the Middle East and 11·3% (10·8-11·9) in Latin America and Caribbean. By 2050, 89 (43·6%) of 204 countries and territories will have an age-standardised rate greater than 10%. Interpretation: Diabetes remains a substantial public health issue. Type 2 diabetes, which makes up the bulk of diabetes cases, is largely preventable and, in some cases, potentially reversible if identified and managed early in the disease course. However, all evidence indicates that diabetes prevalence is increasing worldwide, primarily due to a rise in obesity caused by multiple factors. Preventing and controlling type 2 diabetes remains an ongoing challenge. It is essential to better understand disparities in risk factor profiles and diabetes burden across populations, to inform strategies to successfully control diabetes risk factors within the context of multiple and complex drivers. Funding: Bill & Melinda Gates Foundation

Burden of disease scenarios for 204 countries and territories, 2022–2050: a forecasting analysis for the Global Burden of Disease Study 2021

BackgroundFuture trends in disease burden and drivers of health are of great interest to policy makers and the public at large. This information can be used for policy and long-term health investment, planning, and prioritisation. We have expanded and improved upon previous forecasts produced as part of the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) and provide a reference forecast (the most likely future), and alternative scenarios assessing disease burden trajectories if selected sets of risk factors were eliminated from current levels by 2050.MethodsUsing forecasts of major drivers of health such as the Socio-demographic Index (SDI; a composite measure of lag-distributed income per capita, mean years of education, and total fertility under 25 years of age) and the full set of risk factor exposures captured by GBD, we provide cause-specific forecasts of mortality, years of life lost (YLLs), years lived with disability (YLDs), and disability-adjusted life-years (DALYs) by age and sex from 2022 to 2050 for 204 countries and territories, 21 GBD regions, seven super-regions, and the world. All analyses were done at the cause-specific level so that only risk factors deemed causal by the GBD comparative risk assessment influenced future trajectories of mortality for each disease. Cause-specific mortality was modelled using mixed-effects models with SDI and time as the main covariates, and the combined impact of causal risk factors as an offset in the model. At the all-cause mortality level, we captured unexplained variation by modelling residuals with an autoregressive integrated moving average model with drift attenuation. These all-cause forecasts constrained the cause-specific forecasts at successively deeper levels of the GBD cause hierarchy using cascading mortality models, thus ensuring a robust estimate of cause-specific mortality. For non-fatal measures (eg, low back pain), incidence and prevalence were forecasted from mixed-effects models with SDI as the main covariate, and YLDs were computed from the resulting prevalence forecasts and average disability weights from GBD. Alternative future scenarios were constructed by replacing appropriate reference trajectories for risk factors with hypothetical trajectories of gradual elimination of risk factor exposure from current levels to 2050. The scenarios were constructed from various sets of risk factors: environmental risks (Safer Environment scenario), risks associated with communicable, maternal, neonatal, and nutritional diseases (CMNNs; Improved Childhood Nutrition and Vaccination scenario), risks associated with major non-communicable diseases (NCDs; Improved Behavioural and Metabolic Risks scenario), and the combined effects of these three scenarios. Using the Shared Socioeconomic Pathways climate scenarios SSP2-4.5 as reference and SSP1-1.9 as an optimistic alternative in the Safer Environment scenario, we accounted for climate change impact on health by using the most recent Intergovernmental Panel on Climate Change temperature forecasts and published trajectories of ambient air pollution for the same two scenarios. Life expectancy and healthy life expectancy were computed using standard methods. The forecasting framework includes computing the age-sex-specific future population for each location and separately for each scenario. 95% uncertainty intervals (UIs) for each individual future estimate were derived from the 2·5th and 97·5th percentiles of distributions generated from propagating 500 draws through the multistage computational pipeline.FindingsIn the reference scenario forecast, global and super-regional life expectancy increased from 2022 to 2050, but improvement was at a slower pace than in the three decades preceding the COVID-19 pandemic (beginning in 2020). Gains in future life expectancy were forecasted to be greatest in super-regions with comparatively low life expectancies (such as sub-Saharan Africa) compared with super-regions with higher life expectancies (such as the high-income super-region), leading to a trend towards convergence in life expectancy across locations between now and 2050. At the super-region level, forecasted healthy life expectancy patterns were similar to those of life expectancies. Forecasts for the reference scenario found that health will improve in the coming decades, with all-cause age-standardised DALY rates decreasing in every GBD super-region. The total DALY burden measured in counts, however, will increase in every super-region, largely a function of population ageing and growth. We also forecasted that both DALY counts and age-standardised DALY rates will continue to shift from CMNNs to NCDs, with the most pronounced shifts occurring in sub-Saharan Africa (60·1% [95% UI 56·8–63·1] of DALYs were from CMNNs in 2022 compared with 35·8% [31·0–45·0] in 2050) and south Asia (31·7% [29·2–34·1] to 15·5% [13·7–17·5]). This shift is reflected in the leading global causes of DALYs, with the top four causes in 2050 being ischaemic heart disease, stroke, diabetes, and chronic obstructive pulmonary disease, compared with 2022, with ischaemic heart disease, neonatal disorders, stroke, and lower respiratory infections at the top. The global proportion of DALYs due to YLDs likewise increased from 33·8% (27·4–40·3) to 41·1% (33·9–48·1) from 2022 to 2050, demonstrating an important shift in overall disease burden towards morbidity and away from premature death. The largest shift of this kind was forecasted for sub-Saharan Africa, from 20·1% (15·6–25·3) of DALYs due to YLDs in 2022 to 35·6% (26·5–43·0) in 2050. In the assessment of alternative future scenarios, the combined effects of the scenarios (Safer Environment, Improved Childhood Nutrition and Vaccination, and Improved Behavioural and Metabolic Risks scenarios) demonstrated an important decrease in the global burden of DALYs in 2050 of 15·4% (13·5–17·5) compared with the reference scenario, with decreases across super-regions ranging from 10·4% (9·7–11·3) in the high-income super-region to 23·9% (20·7–27·3) in north Africa and the Middle East. The Safer Environment scenario had its largest decrease in sub-Saharan Africa (5·2% [3·5–6·8]), the Improved Behavioural and Metabolic Risks scenario in north Africa and the Middle East (23·2% [20·2–26·5]), and the Improved Nutrition and Vaccination scenario in sub-Saharan Africa (2·0% [–0·6 to 3·6]).InterpretationGlobally, life expectancy and age-standardised disease burden were forecasted to improve between 2022 and 2050, with the majority of the burden continuing to shift from CMNNs to NCDs. That said, continued progress on reducing the CMNN disease burden will be dependent on maintaining investment in and policy emphasis on CMNN disease prevention and treatment. Mostly due to growth and ageing of populations, the number of deaths and DALYs due to all causes combined will generally increase. By constructing alternative future scenarios wherein certain risk exposures are eliminated by 2050, we have shown that opportunities exist to substantially improve health outcomes in the future through concerted efforts to prevent exposure to well established risk factors and to expand access to key health interventions.FundingBill & Melinda Gates Foundation.</p

Global fertility in 204 countries and territories, 1950–2021, with forecasts to 2100: a comprehensive demographic analysis for the Global Burden of Disease Study 2021

BackgroundAccurate assessments of current and future fertility—including overall trends and changing population age structures across countries and regions—are essential to help plan for the profound social, economic, environmental, and geopolitical challenges that these changes will bring. Estimates and projections of fertility are necessary to inform policies involving resource and health-care needs, labour supply, education, gender equality, and family planning and support. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2021 produced up-to-date and comprehensive demographic assessments of key fertility indicators at global, regional, and national levels from 1950 to 2021 and forecast fertility metrics to 2100 based on a reference scenario and key policy-dependent alternative scenarios.MethodsTo estimate fertility indicators from 1950 to 2021, mixed-effects regression models and spatiotemporal Gaussian process regression were used to synthesise data from 8709 country-years of vital and sample registrations, 1455 surveys and censuses, and 150 other sources, and to generate age-specific fertility rates (ASFRs) for 5-year age groups from age 10 years to 54 years. ASFRs were summed across age groups to produce estimates of total fertility rate (TFR). Livebirths were calculated by multiplying ASFR and age-specific female population, then summing across ages 10–54 years. To forecast future fertility up to 2100, our Institute for Health Metrics and Evaluation (IHME) forecasting model was based on projections of completed cohort fertility at age 50 years (CCF50; the average number of children born over time to females from a specified birth cohort), which yields more stable and accurate measures of fertility than directly modelling TFR. CCF50 was modelled using an ensemble approach in which three sub-models (with two, three, and four covariates variously consisting of female educational attainment, contraceptive met need, population density in habitable areas, and under-5 mortality) were given equal weights, and analyses were conducted utilising the MR-BRT (meta-regression—Bayesian, regularised, trimmed) tool. To capture time-series trends in CCF50 not explained by these covariates, we used a first-order autoregressive model on the residual term. CCF50 as a proportion of each 5-year ASFR was predicted using a linear mixed-effects model with fixed-effects covariates (female educational attainment and contraceptive met need) and random intercepts for geographical regions. Projected TFRs were then computed for each calendar year as the sum of single-year ASFRs across age groups. The reference forecast is our estimate of the most likely fertility future given the model, past fertility, forecasts of covariates, and historical relationships between covariates and fertility. We additionally produced forecasts for multiple alternative scenarios in each location: the UN Sustainable Development Goal (SDG) for education is achieved by 2030; the contraceptive met need SDG is achieved by 2030; pro-natal policies are enacted to create supportive environments for those who give birth; and the previous three scenarios combined. Uncertainty from past data inputs and model estimation was propagated throughout analyses by taking 1000 draws for past and present fertility estimates and 500 draws for future forecasts from the estimated distribution for each metric, with 95% uncertainty intervals (UIs) given as the 2·5 and 97·5 percentiles of the draws. To evaluate the forecasting performance of our model and others, we computed skill values—a metric assessing gain in forecasting accuracy—by comparing predicted versus observed ASFRs from the past 15 years (2007–21). A positive skill metric indicates that the model being evaluated performs better than the baseline model (here, a simplified model holding 2007 values constant in the future), and a negative metric indicates that the evaluated model performs worse than baseline.FindingsDuring the period from 1950 to 2021, global TFR more than halved, from 4·84 (95% UI 4·63–5·06) to 2·23 (2·09–2·38). Global annual livebirths peaked in 2016 at 142 million (95% UI 137–147), declining to 129 million (121–138) in 2021. Fertility rates declined in all countries and territories since 1950, with TFR remaining above 2·1—canonically considered replacement-level fertility—in 94 (46·1%) countries and territories in 2021. This included 44 of 46 countries in sub-Saharan Africa, which was the super-region with the largest share of livebirths in 2021 (29·2% [28·7–29·6]). 47 countries and territories in which lowest estimated fertility between 1950 and 2021 was below replacement experienced one or more subsequent years with higher fertility; only three of these locations rebounded above replacement levels. Future fertility rates were projected to continue to decline worldwide, reaching a global TFR of 1·83 (1·59–2·08) in 2050 and 1·59 (1·25–1·96) in 2100 under the reference scenario. The number of countries and territories with fertility rates remaining above replacement was forecast to be 49 (24·0%) in 2050 and only six (2·9%) in 2100, with three of these six countries included in the 2021 World Bank-defined low-income group, all located in the GBD super-region of sub-Saharan Africa. The proportion of livebirths occurring in sub-Saharan Africa was forecast to increase to more than half of the world's livebirths in 2100, to 41·3% (39·6–43·1) in 2050 and 54·3% (47·1–59·5) in 2100. The share of livebirths was projected to decline between 2021 and 2100 in most of the six other super-regions—decreasing, for example, in south Asia from 24·8% (23·7–25·8) in 2021 to 16·7% (14·3–19·1) in 2050 and 7·1% (4·4–10·1) in 2100—but was forecast to increase modestly in the north Africa and Middle East and high-income super-regions. Forecast estimates for the alternative combined scenario suggest that meeting SDG targets for education and contraceptive met need, as well as implementing pro-natal policies, would result in global TFRs of 1·65 (1·40–1·92) in 2050 and 1·62 (1·35–1·95) in 2100. The forecasting skill metric values for the IHME model were positive across all age groups, indicating that the model is better than the constant prediction.InterpretationFertility is declining globally, with rates in more than half of all countries and territories in 2021 below replacement level. Trends since 2000 show considerable heterogeneity in the steepness of declines, and only a small number of countries experienced even a slight fertility rebound after their lowest observed rate, with none reaching replacement level. Additionally, the distribution of livebirths across the globe is shifting, with a greater proportion occurring in the lowest-income countries. Future fertility rates will continue to decline worldwide and will remain low even under successful implementation of pro-natal policies. These changes will have far-reaching economic and societal consequences due to ageing populations and declining workforces in higher-income countries, combined with an increasing share of livebirths among the already poorest regions of the world.</p