The rise of India’s global health diplomacy amid COVID-19 pandemic

The COVID-19 pandemic has highlighted the importance of global health diplomacy (GHD), with India emerging as a key player. India’s commitment to GHD is demonstrated by its active participation in regional and multilateral projects, pharmaceutical expertise, and large-scale manufacturing capabilities, which include the production and distribution of COVID-19 vaccines and essential medicines. India has supported nations in need through bilateral and multilateral platforms, providing vaccines to countries experiencing shortages and offering technical assistance and capacity-building programs to improve healthcare infrastructure and response capabilities. India’s unique approach to GHD, rooted in humanitarian diplomacy, emphasized collaboration and empathy and stressed the well-being of humanity by embracing the philosophy of "Vasudhaiva Kutumbakam," which translates to "the world is one family." Against this background, this paper’s main focus is to analyze the rise of India’s GHD amidst the COVID-19 pandemic and its leadership in addressing various global challenges. India has demonstrated its commitment to global solidarity by offering medical supplies, equipment, and expertise to more than 100 countries. India’s rising global leadership can be attributed to its proactive approach, humanitarian diplomacy, and significant contributions to global health initiatives

Measuring the availability of human resources for health and its relationship to universal health coverage for 204 countries and territories from 1990 to 2019: a systematic analysis for the Global Burden of Disease Study 2019

Background: Human resources for health (HRH) include a range of occupations that aim to promote or improve human health. The UN Sustainable Development Goals (SDGs) and the WHO Health Workforce 2030 strategy have drawn attention to the importance of HRH for achieving policy priorities such as universal health coverage (UHC). Although previous research has found substantial global disparities in HRH, the absence of comparable cross-national estimates of existing workforces has hindered efforts to quantify workforce requirements to meet health system goals. We aimed to use comparable and standardised data sources to estimate HRH densities globally, and to examine the relationship between a subset of HRH cadres and UHC effective coverage performance. Methods: Through the International Labour Organization and Global Health Data Exchange databases, we identified 1404 country-years of data from labour force surveys and 69 country-years of census data, with detailed microdata on health-related employment. From the WHO National Health Workforce Accounts, we identified 2950 country-years of data. We mapped data from all occupational coding systems to the International Standard Classification of Occupations 1988 (ISCO-88), allowing for standardised estimation of densities for 16 categories of health workers across the full time series. Using data from 1990 to 2019 for 196 of 204 countries and territories, covering seven Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) super-regions and 21 regions, we applied spatiotemporal Gaussian process regression (ST-GPR) to model HRH densities from 1990 to 2019 for all countries and territories. We used stochastic frontier meta-regression to model the relationship between the UHC effective coverage index and densities for the four categories of health workers enumerated in SDG indicator 3.c.1 pertaining to HRH: physicians, nurses and midwives, dentistry personnel, and pharmaceutical personnel. We identified minimum workforce density thresholds required to meet a specified target of 80 out of 100 on the UHC effective coverage index, and quantified national shortages with respect to those minimum thresholds. Findings: We estimated that, in 2019, the world had 104·0 million (95% uncertainty interval 83·5–128·0) health workers, including 12·8 million (9·7–16·6) physicians, 29·8 million (23·3–37·7) nurses and midwives, 4·6 million (3·6–6·0) dentistry personnel, and 5·2 million (4·0–6·7) pharmaceutical personnel. We calculated a global physician density of 16·7 (12·6–21·6) per 10 000 population, and a nurse and midwife density of 38·6 (30·1–48·8) per 10 000 population. We found the GBD super-regions of sub-Saharan Africa, south Asia, and north Africa and the Middle East had the lowest HRH densities. To reach 80 out of 100 on the UHC effective coverage index, we estimated that, per 10 000 population, at least 20·7 physicians, 70·6 nurses and midwives, 8·2 dentistry personnel, and 9·4 pharmaceutical personnel would be needed. In total, the 2019 national health workforces fell short of these minimum thresholds by 6·4 million physicians, 30·6 million nurses and midwives, 3·3 million dentistry personnel, and 2·9 million pharmaceutical personnel. Interpretation: Considerable expansion of the world's health workforce is needed to achieve high levels of UHC effective coverage. The largest shortages are in low-income settings, highlighting the need for increased financing and coordination to train, employ, and retain human resources in the health sector. Actual HRH shortages might be larger than estimated because minimum thresholds for each cadre of health workers are benchmarked on health systems that most efficiently translate human resources into UHC attainment

Mapping development and health effects of cooking with solid fuels in low-income and middle-income countries, 2000-18 : a geospatial modelling study

Background More than 3 billion people do not have access to clean energy and primarily use solid fuels to cook. Use of solid fuels generates household air pollution, which was associated with more than 2 million deaths in 2019. Although local patterns in cooking vary systematically, subnational trends in use of solid fuels have yet to be comprehensively analysed. We estimated the prevalence of solid-fuel use with high spatial resolution to explore subnational inequalities, assess local progress, and assess the effects on health in low-income and middle-income countries (LMICs) without universal access to clean fuels.Methods We did a geospatial modelling study to map the prevalence of solid-fuel use for cooking at a 5 km x 5 km resolution in 98 LMICs based on 2.1 million household observations of the primary cooking fuel used from 663 population-based household surveys over the years 2000 to 2018. We use observed temporal patterns to forecast household air pollution in 2030 and to assess the probability of attaining the Sustainable Development Goal (SDG) target indicator for clean cooking. We aligned our estimates of household air pollution to geospatial estimates of ambient air pollution to establish the risk transition occurring in LMICs. Finally, we quantified the effect of residual primary solid-fuel use for cooking on child health by doing a counterfactual risk assessment to estimate the proportion of deaths from lower respiratory tract infections in children younger than 5 years that could be associated with household air pollution.Findings Although primary reliance on solid-fuel use for cooking has declined globally, it remains widespread. 593 million people live in districts where the prevalence of solid-fuel use for cooking exceeds 95%. 66% of people in LMICs live in districts that are not on track to meet the SDG target for universal access to clean energy by 2030. Household air pollution continues to be a major contributor to particulate exposure in LMICs, and rising ambient air pollution is undermining potential gains from reductions in the prevalence of solid-fuel use for cooking in many countries. We estimated that, in 2018, 205000 (95% uncertainty interval 147000-257000) children younger than 5 years died from lower respiratory tract infections that could be attributed to household air pollution.Interpretation Efforts to accelerate the adoption of clean cooking fuels need to be substantially increased and recalibrated to account for subnational inequalities, because there are substantial opportunities to improve air quality and avert child mortality associated with household air pollution. Copyright (C) 2022 The Author(s). Published by Elsevier Ltd.Peer reviewe

Global, regional, and national sex-specific burden and control of the HIV epidemic, 1990-2019, for 204 countries and territories: the Global Burden of Diseases Study 2019

Background: The sustainable development goals (SDGs) aim to end HIV/AIDS as a public health threat by 2030. Understanding the current state of the HIV epidemic and its change over time is essential to this effort. This study assesses the current sex-specific HIV burden in 204 countries and territories and measures progress in the control of the epidemic. Methods: To estimate age-specific and sex-specific trends in 48 of 204 countries, we extended the Estimation and Projection Package Age-Sex Model to also implement the spectrum paediatric model. We used this model in cases where age and sex specific HIV-seroprevalence surveys and antenatal care-clinic sentinel surveillance data were available. For the remaining 156 of 204 locations, we developed a cohort-incidence bias adjustment to derive incidence as a function of cause-of-death data from vital registration systems. The incidence was input to a custom Spectrum model. To assess progress, we measured the percentage change in incident cases and deaths between 2010 and 2019 (threshold >75% decline), the ratio of incident cases to number of people living with HIV (incidence-to-prevalence ratio threshold <0·03), and the ratio of incident cases to deaths (incidence-to-mortality ratio threshold <1·0). Findings: In 2019, there were 36·8 million (95% uncertainty interval [UI] 35·1–38·9) people living with HIV worldwide. There were 0·84 males (95% UI 0·78–0·91) per female living with HIV in 2019, 0·99 male infections (0·91–1·10) for every female infection, and 1·02 male deaths (0·95–1·10) per female death. Global progress in incident cases and deaths between 2010 and 2019 was driven by sub-Saharan Africa (with a 28·52% decrease in incident cases, 95% UI 19·58–35·43, and a 39·66% decrease in deaths, 36·49–42·36). Elsewhere, the incidence remained stable or increased, whereas deaths generally decreased. In 2019, the global incidence-to-prevalence ratio was 0·05 (95% UI 0·05–0·06) and the global incidence-to-mortality ratio was 1·94 (1·76–2·12). No regions met suggested thresholds for progress. Interpretation: Sub-Saharan Africa had both the highest HIV burden and the greatest progress between 1990 and 2019. The number of incident cases and deaths in males and females approached parity in 2019, although there remained more females with HIV than males with HIV. Globally, the HIV epidemic is far from the UNAIDS benchmarks on progress metrics. Funding: The Bill & Melinda Gates Foundation, the National Institute of Mental Health of the US National Institutes of Health (NIH), and the National Institute on Aging of the NIH

Measuring routine childhood vaccination coverage in 204 countries and territories, 1980-2019 : a systematic analysis for the Global Burden of Disease Study 2020, Release 1

Background Measuring routine childhood vaccination is crucial to inform global vaccine policies and programme implementation, and to track progress towards targets set by the Global Vaccine Action Plan (GVAP) and Immunization Agenda 2030. Robust estimates of routine vaccine coverage are needed to identify past successes and persistent vulnerabilities. Drawing from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2020, Release 1, we did a systematic analysis of global, regional, and national vaccine coverage trends using a statistical framework, by vaccine and over time. Methods For this analysis we collated 55 326 country-specific, cohort-specific, year-specific, vaccine-specific, and dosespecific observations of routine childhood vaccination coverage between 1980 and 2019. Using spatiotemporal Gaussian process regression, we produced location-specific and year-specific estimates of 11 routine childhood vaccine coverage indicators for 204 countries and territories from 1980 to 2019, adjusting for biases in countryreported data and reflecting reported stockouts and supply disruptions. We analysed global and regional trends in coverage and numbers of zero-dose children (defined as those who never received a diphtheria-tetanus-pertussis [DTP] vaccine dose), progress towards GVAP targets, and the relationship between vaccine coverage and sociodemographic development. Findings By 2019, global coverage of third-dose DTP (DTP3; 81.6% [95% uncertainty interval 80.4-82 .7]) more than doubled from levels estimated in 1980 (39.9% [37.5-42.1]), as did global coverage of the first-dose measles-containing vaccine (MCV1; from 38.5% [35.4-41.3] in 1980 to 83.6% [82.3-84.8] in 2019). Third- dose polio vaccine (Pol3) coverage also increased, from 42.6% (41.4-44.1) in 1980 to 79.8% (78.4-81.1) in 2019, and global coverage of newer vaccines increased rapidly between 2000 and 2019. The global number of zero-dose children fell by nearly 75% between 1980 and 2019, from 56.8 million (52.6-60. 9) to 14.5 million (13.4-15.9). However, over the past decade, global vaccine coverage broadly plateaued; 94 countries and territories recorded decreasing DTP3 coverage since 2010. Only 11 countries and territories were estimated to have reached the national GVAP target of at least 90% coverage for all assessed vaccines in 2019. Interpretation After achieving large gains in childhood vaccine coverage worldwide, in much of the world this progress was stalled or reversed from 2010 to 2019. These findings underscore the importance of revisiting routine immunisation strategies and programmatic approaches, recentring service delivery around equity and underserved populations. Strengthening vaccine data and monitoring systems is crucial to these pursuits, now and through to 2030, to ensure that all children have access to, and can benefit from, lifesaving vaccines. Copyright (C) 2021 The Author(s). Published by Elsevier Ltd.Peer reviewe

Global injury morbidity and mortality from 1990 to 2017 : results from the Global Burden of Disease Study 2017

Correction:Background Past research in population health trends has shown that injuries form a substantial burden of population health loss. Regular updates to injury burden assessments are critical. We report Global Burden of Disease (GBD) 2017 Study estimates on morbidity and mortality for all injuries. Methods We reviewed results for injuries from the GBD 2017 study. GBD 2017 measured injury-specific mortality and years of life lost (YLLs) using the Cause of Death Ensemble model. To measure non-fatal injuries, GBD 2017 modelled injury-specific incidence and converted this to prevalence and years lived with disability (YLDs). YLLs and YLDs were summed to calculate disability-adjusted life years (DALYs). Findings In 1990, there were 4 260 493 (4 085 700 to 4 396 138) injury deaths, which increased to 4 484 722 (4 332 010 to 4 585 554) deaths in 2017, while age-standardised mortality decreased from 1079 (1073 to 1086) to 738 (730 to 745) per 100 000. In 1990, there were 354 064 302 (95% uncertainty interval: 338 174 876 to 371 610 802) new cases of injury globally, which increased to 520 710 288 (493 430 247 to 547 988 635) new cases in 2017. During this time, age-standardised incidence decreased non-significantly from 6824 (6534 to 7147) to 6763 (6412 to 7118) per 100 000. Between 1990 and 2017, age-standardised DALYs decreased from 4947 (4655 to 5233) per 100 000 to 3267 (3058 to 3505). Interpretation Injuries are an important cause of health loss globally, though mortality has declined between 1990 and 2017. Future research in injury burden should focus on prevention in high-burden populations, improving data collection and ensuring access to medical care.Peer reviewe

Estimating global injuries morbidity and mortality : methods and data used in the Global Burden of Disease 2017 study

Background While there is a long history of measuring death and disability from injuries, modern research methods must account for the wide spectrum of disability that can occur in an injury, and must provide estimates with sufficient demographic, geographical and temporal detail to be useful for policy makers. The Global Burden of Disease (GBD) 2017 study used methods to provide highly detailed estimates of global injury burden that meet these criteria. Methods In this study, we report and discuss the methods used in GBD 2017 for injury morbidity and mortality burden estimation. In summary, these methods included estimating cause-specific mortality for every cause of injury, and then estimating incidence for every cause of injury. Non-fatal disability for each cause is then calculated based on the probabilities of suffering from different types of bodily injury experienced. Results GBD 2017 produced morbidity and mortality estimates for 38 causes of injury. Estimates were produced in terms of incidence, prevalence, years lived with disability, cause-specific mortality, years of life lost and disability-adjusted life-years for a 28-year period for 22 age groups, 195 countries and both sexes. Conclusions GBD 2017 demonstrated a complex and sophisticated series of analytical steps using the largest known database of morbidity and mortality data on injuries. GBD 2017 results should be used to help inform injury prevention policy making and resource allocation. We also identify important avenues for improving injury burden estimation in the future.Peer reviewe

Predicting the environmental suitability for onchocerciasis in Africa as an aid to elimination planning

Recent evidence suggests that, in some foci, elimination of onchocerciasis from Africa may be feasible with mass drug administration (MDA) of ivermectin. To achieve continental elimination of transmission, mapping surveys will need to be conducted across all implementation units (IUs) for which endemicity status is currently unknown. Using boosted regression tree models with optimised hyperparameter selection, we estimated environmental suitability for onchocerciasis at the 5 × 5-km resolution across Africa. In order to classify IUs that include locations that are environmentally suitable, we used receiver operating characteristic (ROC) analysis to identify an optimal threshold for suitability concordant with locations where onchocerciasis has been previously detected. This threshold value was then used to classify IUs (more suitable or less suitable) based on the location within the IU with the largest mean prediction. Mean estimates of environmental suitability suggest large areas across West and Central Africa, as well as focal areas of East Africa, are suitable for onchocerciasis transmission, consistent with the presence of current control and elimination of transmission efforts. The ROC analysis identified a mean environmental suitability index of 071 as a threshold to classify based on the location with the largest mean prediction within the IU. Of the IUs considered for mapping surveys, 502% exceed this threshold for suitability in at least one 5 × 5-km location. The formidable scale of data collection required to map onchocerciasis endemicity across the African continent presents an opportunity to use spatial data to identify areas likely to be suitable for onchocerciasis transmission. National onchocerciasis elimination programmes may wish to consider prioritising these IUs for mapping surveys as human resources, laboratory capacity, and programmatic schedules may constrain survey implementation, and possibly delaying MDA initiation in areas that would ultimately qualify.SUPPORTING INFORMATION : FIGURE S1. Data coverage by year. Here we visualise the volume of data used in the analysis by country and year. Larger circles indicate more data inputs. ‘NA’ indicates records for which no year was reported (eg, ‘pre-2000’). https://doi.org/10.1371/journal.pntd.0008824.s001FIGURE S2. Illustration of covariate values for year 2000. Maps were produced using ArcGIS Desktop 10.6. https://doi.org/10.1371/journal.pntd.0008824.s002FIGURE S3. Environmental suitability of onchocerciasis including locations that have received MDA for which no pre-intervention data are available. This plot shows suitability predictions from green (low = 0%) to pink (high = 100%), representing those areas where environmental conditions are most similar to prior pathogen detections. Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s003FIGURE S4. Environmental suitability prediction uncertainty including locations that have received MDA for which no pre-intervention data are available. This plot shows uncertainty associated with environmental suitability predictions colored from blue to red (least to most uncertain). Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s004FIGURE S5. Environmental suitability of onchocerciasis excluding morbidity data. This plot shows suitability predictions from green (low = 0%) to pink (high = 100%), representing those areas where environmental conditions are most similar to prior pathogen detections. Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s005FIGURE S6. Environmental suitability prediction uncertainty excluding morbidity data. This plot shows uncertainty associated with environmental suitability predictions colored from blue to red (least to most uncertain). Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. https://doi.org/10.1371/journal.pntd.0008824.s006FIGURE S7. Covariate Effect Curves for all onchocerciasis occurrences (measures of infection prevalence and disability). On the right set of axes we show the frequency density of the occurrences taking covariate values over 20 bins of the horizontal axis. The left set of axes shows the effect of each on the model, where the mean effect is plotted on the black line and its uncertainty is represented by the upper and lower confidence interval bounds plotted in dark grey. The figures show the fit per covariate relative to the data that correspond to specific values of the covariate. https://doi.org/10.1371/journal.pntd.0008824.s007FIGURE S8. Covariate Effect Curves for all onchocerciasis occurrences (measures of infection prevalence and disability). On the right set of axes we show the frequency density of the occurrences taking covariate values over 20 bins of the horizontal axis. The left set of axes shows the effect of each on the model, where the mean effect is plotted on the black line and its uncertainty is represented by the upper and lower confidence interval bounds plotted in dark grey. https://doi.org/10.1371/journal.pntd.0008824.s008FIGURE S9. ROC analysis for threshold. Results of the area under the receiver operating characteristic (ROC) curve analysis are presented below, with false positive rate (FPR) on the x-axis and true positive rate (TPR) on the y-axis. The red dot on the curve represents the location on the curve that corresponds to a threshold that most closely agreed with the input data. For each of the 100 BRT models, we estimated the optimal threshold that maximised agreement between occurrence inputs (considered true positives) and the mean model predictions as 0·71. https://doi.org/10.1371/journal.pntd.0008824.s009TABLE S1. Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) checklist. https://doi.org/10.1371/journal.pntd.0008824.s010TABLE S2. Total number of occurrence data classified as point and polygon inputs by diagnostic. We present the total number of occurrence points extracted from the input data sources by diagnostic type. ‘Other diagnostics’ include: DEC Patch test; Knott’s Method (Mazotti Test); 2 types of LAMP; blood smears; and urine tests. https://doi.org/10.1371/journal.pntd.0008824.s011TABLE S3. Total number of occurrence data classified as point and polygon inputs by location. https://doi.org/10.1371/journal.pntd.0008824.s012TABLE S4. Covariate information. https://doi.org/10.1371/journal.pntd.0008824.s013TEXT S1. Details outlining construction of occurrence dataset. https://doi.org/10.1371/journal.pntd.0008824.s014TEXT S2. Covariate rationale. https://doi.org/10.1371/journal.pntd.0008824.s015TEXT S3. Boosted regression tree methodology additional details. https://doi.org/10.1371/journal.pntd.0008824.s016APPENDIX S1. Country-level maps and data results. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s017This work was primarily supported by a grant from the Bill & Melinda Gates Foundation OPP1132415 (SIH). Financial support from the Neglected Tropical Disease Modelling Consortium (https://www.ntdmodelling.org/), which is funded by the Bill & Melinda Gates Foundation (grants No. OPP1184344 and OPP1186851), and joint centre funding (grant No. MR/R015600/1) by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement which is also part of the EDCTP2 programme supported by the European Union (MGB).The Neglected Tropical Disease Modelling Consortium which is funded by the Bill & Melinda Gates Foundation, the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement which is also part of the EDCTP2 programme supported by the European Union (MGB).http://www.plosNTDS.orgam2022Medical Microbiolog

Predicting the environmental suitability for onchocerciasis in Africa as an aid to elimination planning

Recent evidence suggests that, in some foci, elimination of onchocerciasis from Africa may be feasible with mass drug administration (MDA) of ivermectin. To achieve continental elimination of transmission, mapping surveys will need to be conducted across all implementation units (IUs) for which endemicity status is currently unknown. Using boosted regression tree models with optimised hyperparameter selection, we estimated environmental suitability for onchocerciasis at the 5 × 5-km resolution across Africa. In order to classify IUs that include locations that are environmentally suitable, we used receiver operating characteristic (ROC) analysis to identify an optimal threshold for suitability concordant with locations where onchocerciasis has been previously detected. This threshold value was then used to classify IUs (more suitable or less suitable) based on the location within the IU with the largest mean prediction. Mean estimates of environmental suitability suggest large areas across West and Central Africa, as well as focal areas of East Africa, are suitable for onchocerciasis transmission, consistent with the presence of current control and elimination of transmission efforts. The ROC analysis identified a mean environmental suitability index of 0.71 as a threshold to classify based on the location with the largest mean prediction within the IU. Of the IUs considered for mapping surveys, 50.2% exceed this threshold for suitability in at least one 5×5-km location. The formidable scale of data collection required to map onchocerciasis endemicity across the African continent presents an opportunity to use spatial data to identify areas likely to be suitable for onchocerciasis transmission. National onchocerciasis elimination programmes may wish to consider prioritising these IUs for mapping surveys as human resources, laboratory capacity, and programmatic schedules may constrain survey implementation, and possibly delaying MDA initiation in areas that would ultimately qualify

Global, regional, and national sex differences in the global burden of tuberculosis by HIV status, 1990–2019: results from the Global Burden of Disease Study 2019

Background Tuberculosis is a major contributor to the global burden of disease, causing more than a million deaths annually. Given an emphasis on equity in access to diagnosis and treatment of tuberculosis in global health targets, evaluations of differences in tuberculosis burden by sex are crucial. We aimed to assess the levels and trends of the global burden of tuberculosis, with an emphasis on investigating differences in sex by HIV status for 204 countries and territories from 1990 to 2019. Methods We used a Bayesian hierarchical Cause of Death Ensemble model (CODEm) platform to analyse 21 505 site-years of vital registration data, 705 site-years of verbal autopsy data, 825 site-years of sample-based vital registration data, and 680 site-years of mortality surveillance data to estimate mortality due to tuberculosis among HIV-negative individuals. We used a population attributable fraction approach to estimate mortality related to HIV and tuberculosis coinfection. A compartmental meta-regression tool (DisMod-MR 2.1) was then used to synthesise all available data sources, including prevalence surveys, annual case notifications, population-based tuberculin surveys, and tuberculosis cause-specific mortality, to produce estimates of incidence, prevalence, and mortality that were internally consistent. We further estimated the fraction of tuberculosis mortality that is attributable to independent effects of risk factors, including smoking, alcohol use, and diabetes, for HIV-negative individuals. For individuals with HIV and tuberculosis coinfection, we assessed mortality attributable to HIV risk factors including unsafe sex, intimate partner violence (only estimated among females), and injection drug use. We present 95% uncertainty intervals for all estimates. Findings Globally, in 2019, among HIV-negative individuals, there were 1.18 million (95% uncertainty interval 1.08-1.29) deaths due to tuberculosis and 8.50 million (7.45-9.73) incident cases of tuberculosis. Among HIV-positive individuals, there were 217 000 (153 000-279 000) deaths due to tuberculosis and 1.15 million (1.01-1.32) incident cases in 2019. More deaths and incident cases occurred in males than in females among HIV-negative individuals globally in 2019, with 342 000 (234 000-425 000) more deaths and 1.01 million (0.82-1.23) more incident cases in males than in females. Among HIV-positive individuals, 6250 (1820-11 400) more deaths and 81 100 (63 300-100 000) more incident cases occurred among females than among males in 2019. Age-standardised mortality rates among HIV-negative males were more than two times greater in 105 countries and age-standardised incidence rates were more than 1.5 times greater in 74 countries than among HIV-negative females in 2019. The fraction of global tuberculosis deaths among HIV-negative individuals attributable to alcohol use, smoking, and diabetes was 4.27 (3.69-5.02), 6.17 (5.48-7.02), and 1.17 (1.07-1.28) times higher, respectively, among males than among females in 2019. Among individuals with HIV and tuberculosis coinfection, the fraction of mortality attributable to injection drug use was 2.23 (2.03-2.44) times greater among males than females, whereas the fraction due to unsafe sex was 1.06 (1.05-1.08) times greater among females than males. Interpretation As countries refine national tuberculosis programmes and strategies to end the tuberculosis epidemic, the excess burden experienced by males is important. Interventions are needed to actively communicate, especially to men, the importance of early diagnosis and treatment. These interventions should occur in parallel with efforts to minimise excess HIV burden among women in the highest HIV burden countries that are contributing to excess HIV and tuberculosis coinfection burden for females. Placing a focus on tuberculosis burden among HIV-negative males and HIV and tuberculosis coinfection among females might help to diminish the overall burden of tuberculosis. This strategy will be crucial in reaching both equity and burden targets outlined by global health milestone