This paper is in the following e-collection/theme issue:

Longitudinal and Cohort Studies in Public Health (266) JMIR Theme Issue: COVID-19 Special Issue (2394) Health Services Research (413) Outbreak and Pandemic Preparedness and Management (1679) Clinical Informatics (918) Decision Support for Health Professionals (1145)

Published on in Vol 10 (2024)

Preprints (earlier versions) of this paper are available at https://preprints.jmir.org/preprint/44062, first published .
The Use of ICD-9-CM Coding to Identify COVID-19 Diagnoses and Determine Risk Factors for 30-Day Death Rate in Hospitalized Patients in Italy: Retrospective Study

The Use of ICD-9-CM Coding to Identify COVID-19 Diagnoses and Determine Risk Factors for 30-Day Death Rate in Hospitalized Patients in Italy: Retrospective Study

The Use of ICD-9-CM Coding to Identify COVID-19 Diagnoses and Determine Risk Factors for 30-Day Death Rate in Hospitalized Patients in Italy: Retrospective Study

Authors of this article:

Barbara Giordani1 Author Orcid Image ;   Alessandra Burgio2 Author Orcid Image ;   Francesco Grippo2 Author Orcid Image ;   Alessandra Barone1 Author Orcid Image ;   Erica Eugeni1 Author Orcid Image ;   Giovanni Baglio1 Author Orcid Image
Article Authors Cited by Tweetations (3) Metrics

    Original Paper

    1Research, National Outcomes Evaluation Programme (PNE) and International Relations Unit, Italian National Agency for Regional Healthcare Services, Rome, Italy

    2Italian National Institute of Statistics, Rome, Italy

    Corresponding Author:

    Barbara Giordani, MSc

    Research, National Outcomes Evaluation Programme (PNE) and International Relations Unit

    Italian National Agency for Regional Healthcare Services

    Via Piemonte 60

    Rome, 00187

    Italy

    Phone: 39 06 42749713

    Email: giordani@agenas.it


    Background: In Italy, it has been difficult to accurately quantify hospital admissions of patients with a COVID-19 diagnosis using the Hospital Information System (HIS), mainly due to the heterogeneity of codes used in the hospital discharge records during different waves of the COVID-19 pandemic.

    Objective: The objective of this study was to define a specific combination of codes to identify the COVID-19 hospitalizations within the HIS and to investigate the risk factors associated with mortality due to COVID-19 among patients admitted to Italian hospitals in 2020.

    Methods: A retrospective study was conducted using the hospital discharge records, provided by more than 1300 public and private Italian hospitals. Inpatient hospitalizations were detected by implementing an algorithm based on specific International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) code combinations. Hospitalizations were analyzed by different clinical presentations associated with COVID-19 diagnoses. In addition, 2 multivariable Cox regression models were performed among patients hospitalized “due to COVID-19” from January 1 to December 31, 2020, to investigate potential risk factors associated with 30-day death and the temporal changes over the course of the pandemic; in particular, the 30-day death rates during the first and the second waves were analyzed across 3 main geographical areas (North, Center, and South and Islands) and by discharge wards (ordinary and intensive care).

    Results: We identified a total of 325,810 hospitalizations with COVID-19–related diagnosis codes. Among these, 73.4% (n=239,114) were classified as “due to COVID-19,” 14.5% (n=47,416) as “SARS-CoV-2 positive, but not due to COVID-19,” and 12.1% (n=39,280) as “suspected COVID-19” hospitalizations. The cohort of patients hospitalized “due to COVID-19” included 205,048 patients, with a median age of 72 years and a higher prevalence of male patients (n=124,181, 60.6%). The overall 30-day death rate among hospitalized patients due to COVID-19 was 9.9 per 1000 person-days. Mortality was lower for women (hazard ratio [HR]=0.83; P<.001) and for patients coming from high migration pressure countries, especially Northern Africans (HR=0.65; P<.001) and Central and Eastern Europeans (HR=0.66; P<.001), compared to patients coming from Italy and high-income countries. In the southern regions and the Islands, mortality was higher compared to the northern regions (HR=1.17; P<.001), especially during the second wave of COVID-19 among patients with a transfer to intensive care units (HR=2.52; P<.001).

    Conclusions: To our knowledge, the algorithm is the first attempt to define, at a national level, selection criteria for identifying COVID-19 hospitalizations within the HIS. The implemented algorithm will be used to monitor the pandemic over time, and the patients selected in 2020 will be followed up in the next years to assess the long-term effects of COVID-19.

    JMIR Public Health Surveill 2024;10:e44062

    doi:10.2196/44062

    Keywords

    COVID-19; ICD-9-CM coding; hospitalizations; SARS-CoV-2; coronavirus; risk factor; Italy; death rate; monitoring; hospital records; coding; algorithm 


    The COVID-19 pandemic affected global health care systems, leading to an increase in hospitalizations and intensive care admissions and a high risk of severe clinical outcomes, particularly in older adults [1].

    SARS-CoV-2 infection led to clinical manifestations ranging from asymptomatic or paucisymptomatic forms (medical conditions in which patients experience only a small number of symptoms or mild symptoms) to a wide spectrum of symptoms, such as pneumonia, upper respiratory tract infection, or respiratory failures [2,3].

    Some studies showed that the most critically ill patients with COVID-19 in intensive care units (ICUs) required invasive mechanical ventilation and reported a high mortality rate [4]. A meta-analysis of observational studies also found that older age (≥65 years), male gender, and comorbidities (eg, hypertension, cardiovascular diseases, diabetes, chronic obstructive pulmonary diseases, and malignancies) were associated with a greater risk of death from COVID-19 [5]. Although environmental, clinical, and social characteristics have a role as risk factors for SARS-CoV-2 infection and the severity of the disease [6-8], evidence is arising about the importance of human genetics in the onset and course of COVID-19 [9,10].

    In Italy, despite early data provided by the COVID-19 Integrated Surveillance System [11], it has been difficult to accurately quantify hospital admissions of patients with a COVID-19 diagnosis using the Hospital Information System (HIS), mainly due to the heterogeneity of codes used in the hospital discharge records (HDRs) during different waves of the COVID-19 pandemic.

    Indeed, with the beginning of the COVID-19 pandemic, the World Health Organization provided, in response to member state requests, codes and instructions for COVID-19 coding in the International Classification of Diseases, 10th Revision (ICD-10) and International Classification of Diseases, 11th Revision (ICD-11) [12], but not in the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM), which is the international classification of diseases currently used in Italy for coding in the HDRs.

    Therefore, on March 20, 2020, the Italian Ministry of Health (MoH) published the first guidelines to ensure the needed homogeneity in the criteria and methods of HDR coding throughout the country [13]. Such guidelines provide instructions to identify COVID-19 and related clinical conditions by adapting unspecific codes already existing in the ICD-9-CM. Only with the Ministerial Decree of October 28, 2020 [14], specific ICD-9-CM codes were introduced for COVID-19. Annex 1 of the Ministerial Decree and the guidelines issued on February 19, 2021, provided new instructions for both the principal and secondary diagnoses, aiming to overcome differences in coding among different Italian regions. However, the harmonization of COVID-19 information on HDRs has been gradual in different Italian regions, and sometimes, it was done retrospectively by updating and changing codes in HDRs.

    The algorithm presented in this paper, based on the COVID-19 codes included in different MoH guidelines, the ICD-9-CM SARS codes, and the codes that reported a dramatic increase in 2020, has made it possible to quantify patients with COVID-19–related diagnoses admitted to Italian hospitals. Moreover, the algorithm made it possible to classify hospitalizations as “due to COVID-19,” “SARS-CoV-2 positive, but not due to COVID-19,” and “suspected COVID-19,” and investigate the risk factors associated with mortality “due to COVID-19” among patients admitted to Italian hospitals in 2020.

    Moreover, the algorithm is important because other researchers could use it for potential investigations on specific issues related to the pandemic.


    Source

    A nationwide retrospective study was conducted using the HDRs on health care services in 2020, provided by more than 1300 public and private Italian hospitals. Hospital discharge data are routinely collected by the Italian MoH and contain patient demographic information (eg, gender and age), admission and discharge dates, up to 6 discharge diagnoses (ICD-9-CM), up to 11 medical procedures or surgical interventions, and status at discharge (eg, alive, deceased, or transferred to another hospital). In addition, the National Tax Registry was used to determine vital status or death following hospitalization. HDRs were linked with National Tax Registry records using deterministic record linkage.

    Ethical Considerations

    This study is conducted in accordance with the Ministerial Decree, “Evaluation in terms of quality, safety and appropriateness of activities provided for accreditation and contractual agreements with health facilities,” released on December 19, 2022, Article 4 [15]. This decree grants a mandate, inter alia, to the Italian National Agency for Regional Healthcare Services to monitor hospital activity and publish reports based on anonymized Italian hospitalizations. We did not include any identifiable data of patients’ personal information, including name, identity information, address, and telephone number; therefore, ethics approval was not required.

    Algorithm to Identify COVID-19 Hospitalizations

    Inpatient hospitalizations were detected by implementing an algorithm based on specific ICD-9-CM code combinations. This was carried out through the integration of 2 distinct approaches. One approach, deductive (ie, from the general rule to practical use), was based on the ICD-9-CM codes reported in the official documents of the MoH for COVID-19 case identifications. In addition, we used codes referring to SARS (though not specifically due to SARS-CoV-2), assuming that these codes, mainly at the beginning of the pandemic and before the MoH guidelines publication, were used to indicate COVID-19 hospitalizations. These codes are presented in Table 1.

    Table 1. Codes from the deductive approach.
    Codes sourcesCodes
    First guidelines by MoHa (March 20, 2020)
    • 078.89 (other specified diseases due to viruses)
    • V01.79 (contact with or exposure to other viral diseases)
    • V07.0 (need for isolation)
    • V71.83 (observation and evaluation for suspected exposure to other biological agents)
    MoH Decree (October 28, 2020)
    • 043 (COVID-19 disease)
    • 480.4 (COVID-19 pneumonia)
    • 518.9 (COVID-19 acute respiratory distress syndrome)
    • 519.7 (COVID-19 other respiratory infections)
    • V01.85 (exposure to SARS-CoV-2)
    • V01.82 (exposure to SARS-associated coronavirus)
    • V07.00 (need for isolation after contact with SARS-COV-2)
    • V12.04 (personal history of COVID-19)
    • V71.84 (observation and evaluation for suspected exposure to SARS-COV-2)
    Use of ICD-9-CMb codes of SARS
    • 079.82 (SARS-associated coronavirus)
    • 480.3 (pneumonia due to SARS-associated coronavirus)
    • V01.82 (exposure to SARS-associated coronavirus)

    aMoH: Ministry of Health.

    bICD-9-CM: International Classification of Diseases, Ninth Revision, Clinical Modification.

    The other approach, inductive (ie, from an empirical observation to a general rule) was used to confirm and integrate the algorithm. In detail, comparing the number of hospitalizations between 2019 and 2020, a dramatic increase was found in the principal diagnosis with the ICD-9-CM code 484.8 for “pneumonia in other infectious diseases classified elsewhere.” This code was presumed to be associated with COVID-19 and was integrated into the algorithm.

    The hospitalizations were then grouped into clinical profiles according to a hierarchical approach based on 3 main criteria, as follows:

    • The main reason for hospitalization: admissions categorized as “due to COVID-19” (codes identifying COVID-19 in the principal diagnosis) or “SARS-CoV-2 positive, but not due to COVID-19” (codes identifying COVID-19 only in secondary diagnosis).
    • Degree of certainty: diagnosis confirmed or categorized as “suspected COVID-19”; the latter included hospitalizations with nonspecific diagnoses or V codes (exposure to SARS-CoV-2, observation, or isolation), for which it is difficult to ascertain if they refer to COVID-19 cases.
    • Clinical presentation: pneumonia or acute respiratory distress syndrome, other acute respiratory infections, full-blown COVID-19 without respiratory symptoms, and asymptomatic or pauciasymptomatic COVID-19.

    All diagnosis codes and their combined use (algorithm) in the principal and secondary diagnoses used to select and classify COVID-19 hospitalizations are shown in Table S1 in Multimedia Appendix 1.

    Statistical Analysis

    Hospitalizations were analyzed with respect to the distribution of different clinical presentations associated with diagnoses of COVID-19.

    A further analysis was performed only among patients hospitalized “due to COVID-19” from January 1 to December 31, 2020, to investigate potential risk factors associated with 30-day mortality. The cohort included all patients with at least 1 “due to COVID-19” hospitalization, excluding repeated admissions due to transfers to other hospitals or readmissions.

    Time-to-event techniques were used to analyze survival. The follow-up of patients began from the hospital admission date until either the ascertainment of death (outcome) or up to 30 days after admission (censoring).

    The 30-day death rate was calculated as the total number of all-cause deaths up to 30 days after admission divided by the total follow-up days of hospitalized patients in the cohort and was reported per 1000 patient-days. To investigate the association of one or more risk factors with the 30-day death rate among patients hospitalized “due to COVID-19,” a multivariable Cox proportional-hazards regression model (Cox, 1972) was used, after having checked the proportional hazard assumption.

    Sociodemographic data (eg, age, gender, area of residence, or nationality) were considered in the analysis. Patients were classified according to the following geographic areas of nationality: Italy or other high-income countries; migrants from high migration pressure countries (HMPC), which include Central and Eastern European countries (including those belonging to the European Union) and Malta; Northern African countries; other African countries; Central and South American countries, Asian countries (excluding Japan, Israel, and South Korea); and Oceania (excluding Australia and New Zealand) [16].

    The transfer to ICU has been used as a proxy for the clinical severity of the disease.

    The Charlson Comorbidity Index (CCI) [17,18] was used to account for current comorbidities or comorbidities over the previous 4 years in patients that could affect the death rate. To calculate the CCI, 17 medical conditions were considered [19,20]. The CCI was categorized into the following scores: CCI=0, CCI=1 or 2, and CCI≥3.

    To detect differences in the 30-day death rate during the year 2020, three pandemic periods were considered, as follows:

    • First wave, including the Italian lockdown (January 1 to May 5, 2020)
    • Postlockdown period (May 6 to October 7, 2020)
    • Second wave (October 8 to December 31, 2020)

    Two further multivariable Cox regression models were performed to assess temporal changes over the course of the pandemic; in particular, the 30-day death rates during the first and second waves were analyzed across 3 main geographical areas (ie, North, Center, and South and Islands) and by discharge wards (ie, ordinary and intensive care).

    The selection of patients and hospital admissions was performed using SAS Studio (version 3.81; SAS Institute; Enterprise edition). The Cox proportional-hazards regression was performed using STATA (version 14.0; StataCorp). Tests were 2-sided, and statistical significance was set at P<.05.


    COVID-19 Hospitalizations Analysis

    We identified a total of 325,810 hospitalizations with COVID-19-related diagnosis codes within the public and private accredited Italian hospitals throughout the year 2020 (Table 2).

    A total of 239,114 (73.4%) out of the initial HDR selection were classified as “due to COVID-19” hospitalizations (in most cases, with a diagnosis of pneumonia and acute respiratory distress syndrome), and 47,416 (14.5%) were classified as “SARS-CoV-2 positive, but not due to COVID-19” hospitalizations, mainly in association with diseases of the circulatory system, injuries, chronic diseases of respiratory system, and complication of pregnancy, childbirth, and the puerperium.

    The remaining 12.1% (n=39,280) were “suspected COVID-19 hospitalizations.” A high proportion of such hospitalizations (n=29,327, 9%) were also found in association with other pathological conditions reported in principal diagnosis.

    There was a significant increase in the length of ICU stay out of the overall length of hospital stay (calculated in days), compared with the prepandemic period (1,031,037/54,170,360, 1.9% in 2019 vs 1,204,216/46,113,945, 2.6% in 2020; P<.001). In addition, this proportion was 2.6% (33,509/1,295,932) in 2020 for the “SARS-CoV-2 positive, but not due to COVID-19” and “suspected COVID-19” hospitalizations, while it increased to 8.4% (281,783/3,355,386) among the “due to COVID-19” hospitalizations (data are not shown in the Tables).

    Table 2. The number of admissions categorized by reason for hospitalization and clinical presentation of COVID-19 in Italy, 2020.
    Reason for hospitalization and clinical presentationValues (N=325,810), n (%)
    Due to COVID-19, n=239,114 (73.4%)

    		Pneumonia or acute respiratory distress syndrome222,655 (93.1)

    		Other acute respiratory infections2275 (1)

    		Full-blown COVID-19 (without respiratory symptoms)5718 (2.4)

    		Pauciasymptomatic COVID-198466 (3.5)

    		Total239,114 (100)
    SARS-CoV-2–positive but not due to COVID-19 (with secondary diagnoses of COVID-19),n=47,416 (14.5%)

    		Diseases of the circulatory system10,318 (21.8)

    		Injuries6038 (12.7)

    		Chronic diseases of the respiratory system4424 (9.3)

    		Complication of pregnancy, childbirth, and the puerperium4089 (8.6)

    		Diseases of the digestive system3208 (6.8)

    		Neoplasms2993 (6.3)

    		Factors influencing health status2844 (6)

    		Infectious and parasitic diseases2441 (5.1)

    		Symptoms, signs, and ill-defined conditions2328 (4.9)

    		Diseases of the genitourinary system2030 (4.3)

    		Diseases of the nervous system and sense organs1985 (4.2)

    		Diseases of the musculoskeletal system and connective tissue1851 (3.9)

    		Mental disorders1075 (2.3)

    		Other1792 (3.8)

    		Total47,416 (100)
    Due to suspected COVID-19, n=9953 (3.1%)

    		Pneumonia or acute respiratory distress syndrome8277 (83.2)
     Other acute respiratory infections in suspected COVID-19678 (6.8)
     Need for isolation998 (10)
     Total9953 (100)
    Associated with suspected COVID-19 (non–COVID-19 disease), n=29,237 (9%)29,327 (100)
    Total325,810 (100)

    “Due to COVID-19” Patients Analysis

    Table 3 shows the baseline demographic and clinical characteristics of patients hospitalized “due to COVID-19.”

    The cohort included 205,048 patients, with a median age of 72 years. Of these patients, 147,266 (71.8%) were older than 60 years, and there was a higher prevalence of male patients (n=124,181, 60.6%). Additionally, 6.1% (n=12,364) were HMPC citizens coming from Central and Eastern Europe (n=3721, 1.8%), Asia (n=2842, 1.4%), Central and South America (n=2637, 1.3%), Northern Africa (n=1786, 0.9%), and other African countries (n=1378, 0.7%). More than two-thirds of patients resided in the northern Italian regions, while 31.5% (n=64,491) were in Central and Southern Italy. Over the 3 periods of analysis, 39.4% (n=80,843) of patients were hospitalized in the first wave, 53.9% (n=110,446) were hospitalized in the second wave, and only 6.7% (n=13,759) were hospitalized in the postlockdown period. A total of 199,283 (97.2%) patients reported a Charlson Comorbidity Index score lower than or equal to 2.

    The overall 30-day death rate in patients hospitalized “due to COVID-19” was 9.9 per 1000 person-days. After adjustment, mortality was lower for women (HR=0.83; P<.001) and for patients coming from HMPC. Among the latter, North African patients (HR=0.65; P<.001) and patients from Central and Eastern Europe (HR=0.66; P<.001) were around 35% less likely to die within 30 days, compared to patients coming from Italy and other high-income countries; Asian patients as well as Central and Southern American patients were 17% less likely to die within 30 days.

    Table 3. The baseline demographic and clinical characteristics of patients hospitalized “due to COVID-19.”
    Characteristics and categoriesPatients (N=205,048), n (%)The 30-day death rate per 1000 person-days (overall=9.9 per 1000)Adjusted hazard ratio95% CI
    Gender

    		Male124,181 (60.6)10.21N/Aa

    		Female80,867 (39.4)9.40.830.82-0.85
    Age (years)b

    		0-4010,164 (5)0.51N/A

    		41-5015,421 (7.5)1.12.181.80-2.63

    		51-6032,197 (15.7)2.23.773.17-4.49

    		61-7039,080 (19.1)5.88.997.590-10.66

    		>70108,186 (52.7)17.528.9824.49-34.30
    Nationality

    		Italy or other high-income countries192,684 (93.9)10.51N/A

    		North Africa1786 (0.9)2.30.650.54-0.78

    		Other African Countries1378 (0.7)1.40.770.59-1.00

    		Central and South American countries2637 (1.3)2.30.830.71-0.96

    		Asia2842 (1.4)1.60.830.70-0.99

    		Central and Eastern Europe3721 (1.8)2.50.660.59-0.75
    Italian geographic areac

    		North140,557 (68.5)10.21N/A

    		Center34,195 (16.7)8.40.900.88-0.93

    		South and Islands30,296 (14.8)10.31.171.14-1.20
    Periods

    		January 1 to May 5, 202080,843 (39.4)111N/A

    		May 6 to October 7, 202013,759 (6.7)4.40.450.43-0.47

    		October 8 to Dec 31, 2020110,446 (53.9)9.90.850.83-0.86
    Charlson Comorbidity Index

    		096,225 (46.9)6.71

    		1 or 2103,058 (50.3)12.71.371.35-1.40

    		≥35765 (2.8)19.41.681.61-1.76
    Ward of discharge

    		Ordinary care176,801 (86.2)8.51N/A
     Intensive care28,247 (13.8)19.52.552.50-2.60

    aN/A: not applicable.

    bMedian age was 72 years.

    cThe North includes Valle d\'Aosta, Piemonte, Liguria, Lombardia, Trentino-Alto Adige, Friuli-Venezia Giulia, Veneto, and Emilia-Romagna regions. The Center includes Toscana, Umbria, Marche, and Lazio regions. The South and Islands include Abruzzo, Molise, Campania, Puglia, Basilicata, Calabria, Sicilia, and Sardegna regions.

    In addition, we found a higher 30-day survival in the central regions compared with the northern regions, which were the most affected by COVID-19 (HR=0.90; P<.001), while the southern regions and the Islands showed a higher risk of 30-day mortality (HR=1.17; P<.001).

    Transfers to ICUs involved 13.8% (28,247/205,048) of hospitalized patients.

    Patients referred to ICUs were two and a half times more likely to die within 30 days, compared to patients hospitalized in ordinary wards (HR=2.55; P<.001).

    Table 4 shows the survival analysis by geographic area and wards during 2 pandemic waves.

    Over the first wave, among patients in ordinary wards, we found a higher 30-day death rate in the northern regions (10.9 per 1000) compared to the central (7.7 per 1000) and southern islands regions (5.6 per 1000). On the contrary, among patients with a transfer to ICUs, the southern regions and the Islands showed a higher 30-day death rate (28.4 per 1000) compared to the northern regions (15.2 per 1000).

    Throughout the second wave, among patients in ordinary wards, there was a significant decrease in 30-day death rates in the northern regions (8.7 per 1000), whereas such rates increased up to 7.4 per 1000 person-days in the southern regions and the Islands. On the other hand, among patients with a transfer to ICUs, the 30-day death rate increased in all geographic areas, particularly in the southern regions and the Islands, where the number of ICU patients more than tripled from the first to the second wave; this resulted in a significantly higher mortality risk in the southern regions and the Islands compared to the northern ones (HR=2.52; P<.001).

    Table 4. Survival analysis of patients in hospitals due to COVID-19 during the first and the second waves of the pandemic by geographic area and wards of discharge in Italy in 2020.
    Wards and geographic
    area    		First wave (N=80,843)Second wave (N=110,446)
    Patients, n30-day death rate (per 1000 person-days)Adjusted hazard ratioa95% CIPatients, n30-day death rate (per 1000 person-days)Adjusted hazard ratio95% CI
    Ordinary care (first wave: n=68,880; second wave: n=95,915)
     North56,84910.91N/Ab59,4168.71N/A
     Center75917.70.740.70-0.7818,77770.920.88-0.95
     South and Islands44405.60.650.60-0.7117,7227.41.071.03-1.11
    Intensive care unit (first wave: n=11,963; second wave: n=14,531)
     North914415.21N/A7990171N/A
     Center172615.50.920.85-1.01279224.51.501.41-1.60
     South and Islands109328.41.661.52-1.80374942.52.522.39-2.65

    aAdjusted for gender, age class, nationality, and Charlson Comorbidity Index.

    bN/A: not applicable.


    Principal Findings

    In Italian hospitals, due to the lack of COVID-19 codes in the ICD-9-CM system, multiple codes were used for COVID-19 patient admissions throughout the pandemic outbreak; this made it difficult to identify COVID-19 hospitalizations. This study enabled the definition of a specific combination of codes to identify COVID-19 hospitalizations at the national level within the HIS. Additionally, it facilitated the investigation of risk factors associated with mortality due to COVID-19 among patients admitted to Italian hospitals in 2020. The selection algorithm made it possible to identify 239,114 hospitalizations “due to COVID-19” from HDR. This number is very similar to the data on hospital admissions gathered from the COVID-19 Integrated Surveillance System, which amounted to 240,542 [11]. The monthly trend (data are not shown in Tables) is also very similar in the 2 data sources. Only in a few regions, the number of cases derived from HDRs is lower than the number obtained from the Surveillance System, especially in the first wave of the pandemic, likely due to a miscoding of data at the beginning of the epidemiological crisis. Nevertheless, the use of HDR data for COVID-19 impact analyses remains preferable to Surveillance System data for several reasons. The Surveillance System is based on data provided by the regions to monitor the evolution of the pandemic with timely data. Nevertheless, HDR data, although less timely, have higher quality, as they are collected following criteria and rules established by national legislation and guidelines, which guarantee greater homogeneity and comparability across different periods and territories. Moreover, HDR contains more information than the Surveillance System, especially on comorbidities, and allows comparisons with non–COVID-19 hospitalizations as well.

    Overall, we observed a decrease in the 30-day death rate among patients hospitalized “due to COVID-19” between the 2 waves of the pandemic. This cannot be explained by the introduction of new pharmacological therapies [21]; however, it could be due to a better use of existing therapies associated with gained experience in managing the disease by health care providers and a more effective preparedness of the health care system, including reorganization of hospitals and increasing number of ICU beds for patients with COVID-19 [22].

    In particular, in the second wave, mortality was reduced in the northern regions among patients treated in ordinary wards, with a lower proportion of cases referred to ICUs (from 9144/65,993, 13.9.% to 7990/67,406, 11.9%); this may have led to a higher concentration of patients with severe symptoms in ICUs, consistent with a small increase in 30-day death rates (from 15.2 to 17.0 per 1000 person-days). On the contrary, in the southern regions and the Islands, mortality increased between the 2 waves, both in the ordinary wards and in the ICUs, and this might be due to a higher burden on hospital care compared to the first wave, which had affected Southern Italy to a lesser extent. Furthermore, the South reported a decrease in non–COVID-19 hospitalizations, which was greater than the national mean [23], a sign indicating that the hospital system was struggling.

    Our study also highlighted the lower 30-day death rate for women compared to men, regardless of age. The observed disparity may be attributed to sex-based differences in immunological responses, comorbidities prevalence, or differences in behaviours [24,25]. Consistent with several studies, this analysis showed that critically ill older patients were at high risk of disease severity and mortality [26,27]. Some studies also highlighted important post–COVID-19 sequelae in patients with cancer, which adversely affect survival and oncological outcomes after recovery [28].

    Finally, our findings, consistent with the findings of a study carried out in the city of Milan [29], confirmed that mortality among migrants in hospitals due to COVID-19 was lower than that of Italians and people coming from high-income countries. The lower 30-day death rate observed among HMPC migrants can be explained as the effect of specific selection processes (especially in the early and late stages of the migration path), which tend to maintain the overall health status of the foreign population high. One of them is the so-called “healthy migrant effect,” a dynamic for which only people in good health tend to emigrate [30]. Another selection dynamic, known as the “salmon effect,” is due to the habit of older migrants, especially if ill, to go back to their country of origin [31].

    Moreover, an Italian study on data retrieved from the COVID-19 Integrated Surveillance System found an interaction between nationality and the epidemic phases; the analysis stratified by year showed that compared to Italians, mortality and other severe clinical outcomes in non-Italian nationals slightly increased from the last months of 2020 through 2021. The authors hypothesize that reduced access to vaccination by the immigrated population could partly explain such an increase [32], as a high risk of missed vaccination was observed among non–Italian-born people living in Italy, both overall and for individuals ≥50 years of age [33].

    Limitations

    This study has several limitations. First, results may be affected by the information quality of the HDR or by missing or incorrect coding of some diagnoses, leading to possible underestimation of comorbidities. Therefore, to minimize the effect of underreporting bias, comorbidities over the previous 4 years were also retrieved.

    Second, the study is based on the HDR, which collects administrative data and does not contain clinical information of patients, which is useful to better characterize COVID-19 severity. For this reason, the transfer to the ICU has been used in this study as a proxy for the clinical severity of the disease.

    Third, due to the lack of specified COVID-19 codes from the onset of the pandemic, the overall hospitalizations detected through the algorithm might be underestimated.

    Finally, the issue of patients hospitalized due to “suspected COVID-19,” which cannot be defined as COVID-19 cases, remains unsolved.

    Conclusions

    To our knowledge, the algorithm represents the first attempt to define, at the national level, selection criteria for identifying COVID-19 hospitalizations within the HIS.

    While awaiting the adoption of more updated International Classification of Diseases coding systems for hospital diagnoses and procedures in Italy, consistent with the World Health Organization recommendations, the implemented algorithm will be used to monitor the pandemic over time, and the patients selected in 2020 will be followed up in the subsequent years to assess the long-term effects of COVID-19.

    Further analyses will be useful to assess the impact of the anti–COVID-19 vaccination campaign on the severity and mortality of the disease and to investigate possible inequalities between population subgroups.

    Acknowledgments

    We would like to thank Roberto Blaco for his technical support about the COVID-19 coding, Stella Lanzi for the Italian-to-English translation, Simona Gorietti for reference editing, and Elisa Guglielmi for organizational and administrative support. These individuals were not compensated for their role in the study.

    This research received no external funding.

    Data availability

    The data sets generated and analyzed during this study are available from the corresponding author on reasonable request, by contacting direzione.ricerca@agenas.it.

    Authors' Contributions

    All authors discussed the results and contributed to the final manuscript. All authors read and approved the final manuscript.

    Conflicts of Interest

    None declared.

    Multimedia Appendix 1

    Classification of admissions by reason for hospitalisation and clinical presentation of COVID-19.

    DOCX File , 19 KB
    1. Singhal S, Kumar P, Singh S, Saha S, Dey AB. Clinical features and outcomes of COVID-19 in older adults: a systematic review and meta-analysis. BMC Geriatr. May 19, 2021;21(1):321. [FREE Full text] [CrossRef] [Medline]
    2. Wiersinga WJ, Rhodes A, Cheng AC, Peacock SJ, Prescott HC. Pathophysiology, transmission, diagnosis, and treatment of coronavirus disease 2019 (COVID-19): a review. JAMA. Jul 10, 2020.:782-793. [CrossRef] [Medline]
    3. Vetter P, Vu DL, L'Huillier AG, Schibler M, Kaiser L, Jacquerioz F. Clinical features of covid-19. BMJ. Apr 17, 2020;369:m1470. [CrossRef] [Medline]
    4. Grasselli G, Greco M, Zanella A, Albano G, Antonelli M, Bellani G, et al. Risk factors associated with mortality among patients with COVID-19 in intensive care units in Lombardy, Italy. JAMA Intern Med. Oct 01, 2020;180(10):1345-1355. [FREE Full text] [CrossRef] [Medline]
    5. Parohan M, Yaghoubi S, Seraji A, Javanbakht M, Sarraf P, Djalali M. Risk factors for mortality in patients with coronavirus disease 2019 (COVID-19) infection: a systematic review and meta-analysis of observational studies. Aging Male. Dec 2020;23(5):1416-1424. [CrossRef] [Medline]
    6. Docherty AB, Harrison EM, Green CA, Hardwick HE, Pius R, Norman L, et al. Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: prospective observational cohort study. BMJ. May 22, 2020;369:m1985. [FREE Full text] [CrossRef] [Medline]
    7. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. Mar 28, 2020;395(10229):1054-1062. [FREE Full text] [CrossRef] [Medline]
    8. Gao Y, Ding M, Dong X, Zhang J, Kursat Azkur A, Azkur D, et al. Risk factors for severe and critically ill COVID-19 patients: a review. Allergy. Feb 2021;76(2):428-455. [CrossRef] [Medline]
    9. Correale P, Mutti L, Pentimalli F, Baglio G, Saladino RE, Sileri P, et al. HLA-B*44 and C*01 prevalence correlates with Covid19 spreading across Italy. Int J Mol Sci. Jul 23, 2020;21(15):5205. [FREE Full text] [CrossRef] [Medline]
    10. COVID-19 HGI. Mapping the human genetic architecture of COVID-19. Nature. Dec 2021;600(7889):472-477. [CrossRef]
    11. COVID-19 integrated surveillance data in Italy. URL: https://www.epicentro.iss.it/en/coronavirus/sars-cov-2-integrated-surveillance-data [accessed 2022-09-09]
    12. Emergency use ICD codes for COVID-19 disease outbreak. World Health Organization (WHO). URL: https:/​/www.​who.int/​standards/​classifications/​classification-of-diseases/​emergency-use-icd-codes-for-covid-19-disease-outbreak [accessed 2022-08-05]
    13. Linee guida per la codifica della SDO per casi affetti da malattia da Sars-Cov-2 (Covid-19). Nota della Direzione Generale della Programmazione sanitaria del Ministero della Salute. Prot. n. 7648 del 20 marzo 2020. Ministero della Salute. 2020. URL: https:/​/www.​regione.toscana.it/​documents/​10180/​15744217/​Linee+Guida+per+la+codifica+SDO_COVID19.pdf/​8f0b94e7-6751-99d0-4162-c8d31d6e66e2?t=1606819396950
    14. Integrazione dei sistemi di classificazione adottati per la codifica delle informazioni cliniche contenute nella scheda di dimissione ospedaliera e per la remunerazione delle prestazioni ospedaliere in conseguenza della nuova malattia da SARS-CoV-2 (Covid-19). Modifiche al Decreto del 18 dicembre 2008. Ministero della Salute. Decreto 28 ottobre 2020. URL: https://www.gazzettaufficiale.it/eli/id/2021/02/01/21A00441/sg
    15. Valutazione in termini di qualità, sicurezza ed appropriatezza delle attività erogate per l'accreditamento e per gli accordi contrattuali con le strutture sanitarie. Decreto del 19 dicembre 2022. Ministero della Salute. URL: http://tinyurl.com/2eh3jrb7
    16. Istituto Nazionale di Statistica - ISTAT. La presenza straniera in Italia: caratteristiche socio-demografiche. N.10. Roma. Informazioni; 2007.
    17. Charlson M, Szatrowski TP, Peterson J, Gold J. Validation of a combined comorbidity index. J Clin Epidemiol. Nov 1994;47(11):1245-1251. [CrossRef] [Medline]
    18. Deyo RA, Cherkin DC, Ciol MA. Adapting a clinical comorbidity index for use with ICD-9-CM administrative databases. J Clin Epidemiol. Jun 1992;45(6):613-619. [CrossRef] [Medline]
    19. Hall WH, Ramachandran R, Narayan S, Jani AB, Vijayakumar S. An electronic application for rapidly calculating Charlson comorbidity score. BMC Cancer. Dec 20, 2004;4(1):94. [FREE Full text] [CrossRef] [Medline]
    20. Aslam F, Khan NA. Tools for the assessment of comorbidity burden in rheumatoid arthritis. Front Med (Lausanne). Feb 16, 2018;5:39. [FREE Full text] [CrossRef] [Medline]
    21. Sanders JM, Monogue ML, Jodlowski TZ, Cutrell JB. Pharmacologic treatments for coronavirus disease 2019 (COVID-19): a review. JAMA. Apr 13, 2020;18:1824-1836. [CrossRef] [Medline]
    22. Linee di indirizzo organizzative per il potenziamento della rete ospedaliera per emergenza COVID-19 (Art. 2, DL 19 maggio 2020 n. 34). Circolare n. 11254 del Ministero della Salute. 29 maggio 2020. Ministero della Salute. URL: https://reteospedaliera.covid19.difesa.it/
    23. Impatto dell’epidemia COVID-19 sul sistema ospedaliero italiano - anno 2020 [Article in Italian]. Istat. URL: https://www.istat.it/it/archivio/273191 [accessed 2022-07-21]
    24. Mukherjee S, Pahan K. Is COVID-19 gender-sensitive? J Neuroimmune Pharmacol. Mar 06, 2021;16(1):38-47. [FREE Full text] [CrossRef] [Medline]
    25. Jin J, Bai P, He W, Wu F, Liu X, Han D, et al. Gender differences in patients with COVID-19: focus on severity and mortality. Front Public Health. Apr 29, 2020;8:152. [FREE Full text] [CrossRef] [Medline]
    26. Ejaz H, Alsrhani A, Zafar A, Javed H, Junaid K, Abdalla AE, et al. COVID-19 and comorbidities: deleterious impact on infected patients. J Infect Public Health. Dec 2020;13(12):1833-1839. [FREE Full text] [CrossRef] [Medline]
    27. Fang X, Li S, Yu H, Wang P, Zhang Y, Chen Z, et al. Epidemiological, comorbidity factors with severity and prognosis of COVID-19: a systematic review and meta-analysis. Aging (Albany NY). Jul 13, 2020;12(13):12493-12503. [FREE Full text] [CrossRef] [Medline]
    28. Pinato DJ, Tabernero J, Bower M, Scotti L, Patel M, Colomba E, et al. Prevalence and impact of COVID-19 sequelae on treatment and survival of patients with cancer who recovered from SARS-CoV-2 infection: evidence from the OnCovid retrospective, multicentre registry study. Lancet Oncol. Dec 2021;22(12):1669-1680. [FREE Full text] [CrossRef] [Medline]
    29. Giacomelli A, Ridolfo AL, Bonazzetti C, Oreni L, Conti F, Pezzati L, et al. Mortality among Italians and immigrants with COVID-19 hospitalised in Milan, Italy: data from the Luigi Sacco Hospital registry. BMC Infect Dis. Jan 19, 2022;22(1):63. [FREE Full text] [CrossRef] [Medline]
    30. Baglio G, Di Palma R, Eugeni E, Fortino A. Capitolo 8: Gli immigrati irregolari: cosa sappiamo della loro salute? [Undocumented immigrants: what do we know about their health?]. Epidemiol Prev. 2017;41(3-4 (Suppl 1)):57-63. [FREE Full text] [CrossRef] [Medline]
    31. Pablos-Méndez A. Mortality among Hispanics. JAMA. Apr 27, 1994;271(16):1237. [Medline]
    32. Fabiani M, Vescio M, Bressi M, Mateo-Urdiales A, Petrone D, Spuri M, et al. Differences in the incidence and clinical outcomes of SARS-CoV-2 infection between Italian and non-Italian nationals using routine data. Public Health. Oct 2022;211:136-143. [FREE Full text] [CrossRef] [Medline]
    33. Russo A, Tunesi S, Consolazio D, Decarli A, Bergamaschi W. Evaluation of the anti-COVID-19 vaccination campaign in the Metropolitan Area of Milan (Lombardy Region, Northern Italy). Epidemiol Prev. 2021;45(6):568-579. [FREE Full text] [CrossRef] [Medline]

    ‎
    CCI: Charlson Comorbidity Index
    HDR: hospital discharge record
    HIS: Hospital Information System
    HMPC: high migration pressure countries
    HR: hazard ratio
    ICD-9-CM: International Classification of Diseases, Ninth Revision, Clinical Modification
    ICD-10: International Classification of Diseases, 10th Revision
    ICD-11: International Classification of Diseases, 11th Revision
    ICU: intensive care unit
    MoH: Italian Ministry of Health

    Edited by A Mavragani, T Sanchez; submitted 04.11.22; peer-reviewed by F Seccareccia, N Mignolli, H Chen, H Tibble; comments to author 07.07.23; revised version received 27.07.23; accepted 31.10.23; published 23.02.24.

    Copyright

    ©Barbara Giordani, Alessandra Burgio, Francesco Grippo, Alessandra Barone, Erica Eugeni, Giovanni Baglio. Originally published in JMIR Public Health and Surveillance (https://publichealth.jmir.org), 23.02.2024.

    This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Public Health and Surveillance, is properly cited. The complete bibliographic information, a link to the original publication on https://publichealth.jmir.org, as well as this copyright and license information must be included.