Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

COVID-19 Bimodal Clinical and Pathological Phenotypes

View ORCID ProfileSabrina S Batah, Maíra N Benatti, Li Siyuan, Wagner M Telini, Jamile Barbosa, Marcelo B Menezes, View ORCID ProfileTales R Nadai, View ORCID ProfileKeyla S G Sá, Chirag M. Vaswani, Sahil Gupta, View ORCID ProfileDario S Zamboni, View ORCID ProfileDanilo T Wada, View ORCID ProfileRodrigo T Calado, Renê D R Oliveira, Paulo Louzada-Junior, Maria Auxiliadora-Martins, View ORCID ProfileFlávio P Veras, Larissa D Cunha, View ORCID ProfileThiago M Cunha, Rodrigo Luppino-Assad, Marcelo L Balancin, Sirlei S Morais, Ronaldo B Martins, Eurico Arruda, View ORCID ProfileFernando Chahud, View ORCID ProfileMarcel Koenigkam-Santos, Andrea A Cetlin, Fernando Q Cunha, Claudia dos Santos, View ORCID ProfileVera L Capelozzi, Junya Fukuoka, View ORCID ProfileRosane Duarte-Achcar, View ORCID ProfileAlexandre T Fabro
doi: https://doi.org/10.1101/2021.09.03.21262841
Sabrina S Batah
1Department of Pathology and Legal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sabrina S Batah
Maíra N Benatti
2Pulmonary Division, Department of Internal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Li Siyuan
3Department of Surgery, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wagner M Telini
1Department of Pathology and Legal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jamile Barbosa
1Department of Pathology and Legal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marcelo B Menezes
2Pulmonary Division, Department of Internal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tales R Nadai
4Hospital Estadual de Bauru, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Tales R Nadai
Keyla S G Sá
5Department of Cell and Molecular Biology and Pathogenic Bioagents, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Keyla S G Sá
Chirag M. Vaswani
6Department of Physiology, Temerty Faculty of Medicine, University of Toronto, Toronto, ON, CA
7Keenan Research Centre for Biomedical Science, St. Michael’s Hospital, Toronto, ON, CA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sahil Gupta
7Keenan Research Centre for Biomedical Science, St. Michael’s Hospital, Toronto, ON, CA
8Institute of Medical Science, Temerty Faculty of Medicine, University of Toronto, King’s College Circle, Toronto, ON, CA
9Department of Critical Care Medicine, St. Michael’s Hospital, Toronto, ON, CA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dario S Zamboni
5Department of Cell and Molecular Biology and Pathogenic Bioagents, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Dario S Zamboni
Danilo T Wada
10Department of Medical Images, Hematology and Oncology, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Danilo T Wada
Rodrigo T Calado
10Department of Medical Images, Hematology and Oncology, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Rodrigo T Calado
Renê D R Oliveira
11Division of Clinical Immunology, Emergency, Infectious Diseases and Intensive Care Unit, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paulo Louzada-Junior
11Division of Clinical Immunology, Emergency, Infectious Diseases and Intensive Care Unit, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Maria Auxiliadora-Martins
12Division of Intensive Care Medicine, Department of Surgery and Anatomy, Ribeirão Preto Medical School, University of São Paulo; Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Flávio P Veras
13Department of Pharmacology, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Flávio P Veras
Larissa D Cunha
5Department of Cell and Molecular Biology and Pathogenic Bioagents, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Thiago M Cunha
13Department of Pharmacology, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Thiago M Cunha
Rodrigo Luppino-Assad
14Department of Internal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marcelo L Balancin
15Department of Pathology, Faculty of Medicine, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sirlei S Morais
1Department of Pathology and Legal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ronaldo B Martins
5Department of Cell and Molecular Biology and Pathogenic Bioagents, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eurico Arruda
5Department of Cell and Molecular Biology and Pathogenic Bioagents, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fernando Chahud
1Department of Pathology and Legal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Fernando Chahud
Marcel Koenigkam-Santos
10Department of Medical Images, Hematology and Oncology, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Marcel Koenigkam-Santos
Andrea A Cetlin
2Pulmonary Division, Department of Internal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fernando Q Cunha
13Department of Pharmacology, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Claudia dos Santos
7Keenan Research Centre for Biomedical Science, St. Michael’s Hospital, Toronto, ON, CA
16Interdepartmental Division of Critical Care Medicine, University of Toronto, Toronto, ON, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Vera L Capelozzi
15Department of Pathology, Faculty of Medicine, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Vera L Capelozzi
Junya Fukuoka
17Department of Pathology, Nagasaki University Graduate School of Biomedical Sciences, Nagasaki, Japan
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rosane Duarte-Achcar
18National Jewish Health, Department of Medicine, Pathology Division, Denver, Colorado, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Rosane Duarte-Achcar
Alexandre T Fabro
1Department of Pathology and Legal Medicine, Ribeirão Preto Medical School, University of São Paulo, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Alexandre T Fabro
  • For correspondence: fabro{at}fmrp.usp.br
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

ABSTRACT

Background Patients with coronavirus disease-2019 (COVID-19) present varying clinical complications. Different viral load and host response related to genetic and immune background are probably the reasons for these differences. We aimed to sought clinical and pathological correlation that justifies the different clinical outcomes among COVID-19 autopsies cases.

Methods Minimally invasive autopsy was performed on forty-seven confirmed COVID-19 patients from May-July, 2020. Electronic medical record of all patients was collected and a comprehensive histopathological evaluation was performed. Immunohistochemistry, immunofluorescence, special stain, western blotting and post-mortem real-time reverse transcriptase polymerase chain reaction on fresh lung tissue were performed.

Results We show that 5/47 (10,6%) patients present a progressive decline in oxygenation index for acute respiratory distress syndrome (PaO2/FiO2 ratio), low compliance levels, interstitial fibrosis, high α-SMA+ cells/protein expression, high collagens I/III deposition and NETs(P<0.05), named as fibrotic phenotype (N=5). Conversely, 10/47 (21,2%) patients demonstrated progressive increase in PaO2/FiO2 ratio, high pulmonary compliance levels, preserved elastic framework, increase thrombus formation and high platelets and D-dimer levels at admission (P<0.05), named as thrombotic phenotype. While 32/47 (68,1%) had a mixed phenotypes between both ones.

Conclusions We believe that categorization of patients based on these two phenotypes can be used to develop prognostic tools and potential therapies since the PaO2/FiO2 ratio variation and D-dimer levels correlate with the underlying fibrotic or thrombotic pathologic process, respectively; which may indicate possible clinical outcome of the patient.

  • COVID-19
  • COVID-19 phenotypes
  • clinical outcome
  • fibrosis
  • thrombosis

INTRODUCTION

The Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) was first reported in Wuhan in China in 20191. Infected patients can experience the less severe typical flu-like symptoms with fever, dry cough, and dyspnea2. However, in more severe cases, patients require intensive care unit (ICU) admission and mechanical ventilatory support. Acute respiratory distress syndrome (ARDS) is a major complication for severely infected patients and can lead to multi-organ failure such as acute kidney injury, cardiac complications and thromboembolic events 3,4. Differences in these outcomes can be attributed to a differential host response to infection with variable viral load, age, gender, comorbidities, genetic and immune background. While different studies have stratified heterogenous COVID-19 patients based on ABO blood grouping5, cancer subtypes6, hyper-inflammatory 7,8, the characterization of histopathological phenotypes in COVID-19 patients could reveal the underlying pathophysiological process, but it remains largely unreported.

Studies describing the clinical progression of patients who died from COVID-19 correlated to their histopathological pattern are still lacking. We hypothesized that SARS-CoV-2 infection can lead to different mechanisms of lung injury/repair with different clinical and ventilatory manifestations, even though they all lead to patient death. Therefore, we aimed to assess whether there is a clinical and pathological basis that could explain the different clinical outcomes among COVID-19 patients with commonly associated long-term lung dysfunction. Here we present a unique histopathologic analysis of COVID-19 autopsy cases with clinical and pathological correlation and demonstrate a distinct bimodal phenotypic characterization.

METHODS

Study Design

Forty-seven consecutive COVID-19 patients with positive nasopharyngeal swab for SARS-CoV-2 by reverse transcription polymerase chain reaction (RT-PCR) were considered eligible for this study. Modified minimally invasive autopsy (MIA)9 was performed at University Hospital of Ribeirão Preto Medical School, University of São Paulo, Ribeirão Preto, SP, Brazil – HCFMRP/USP from May to July, 2020. Briefly, all MIAs were performed at bedside through post-mortem surgical lung biopsy within 1 hour of death by a 3 cm incision on the anterior side of the chest between the fourth and fifth ribs. A matching 14-gauge cutting needle (Magnum Needles, Bard) and a biopsy gun (Magnum, Bard) also were used. Two post-COVID-19 biopsies cases were included. This study was approved by the local Research Ethics Committee and written informed consent was waived.

Data Collection

Electronic medical record of all patients (N=47) enrolled in this study was collected. Demographic, clinical data, symptoms, drugs, treatments prescribed and laboratory tests were recorded thoroughly. Laboratory tests were collected and identified as: 1) Admission tests: up to the first 48 hours of admission; and 2) Daily tests: on the fifth, third and first day before death.

In order to better understand patient’s outcome, mechanical ventilation parameters of each individual patient were recorded. All daily change of these variables (from the first to the last day of mechanical ventilation until death) were registered and the means recorded for the last day before death of each individual. In addition, PaO2/FiO2 linear regression was performed with all available values. Then, to determine which patients had changes of lung function, we calculated the larger and smaller angles from the lines of PaO2/FiO2 linear regression. Finally, chest X-ray (CXR) and high-resolution computed tomography (HRCT) images of all patients were reviewed and evaluated by two specialized pulmonary radiologists (MKS and DTW). Unfortunately, only 15 patients had HRCT scan, preventing further analyses.

Histological Evaluation

Lung tissue samples were fixed in 10% buffered formalin, for at least 48 h. Paraffin-embedded sections of 3-μm thickness were stained with Haematoxylin and Eosin, Picrosirius red (Abcam, ab150681) for collagen fibers identification and Verhoeff (Abcam, ab150667) for elastic fibers identification, according to local protocol. Immunohistochemistry for anti-alpha smooth muscle actin (α-SMA) (Abcam, ab5694) and anti-SARS-CoV-2 polyclonal antibody, developed by our group for in situ detection of SARS-CoV-2, were performed in paraffin-embedded sections of 3-μm thickness, following our lab protocol10. Histological evaluation was performed by specialized pulmonary pathologists (MLB, VLC, ATF) blinded to clinical history. For all patients, histopathological features were assessed as present or absent. Cut-offs of area involvement were determined for fibrotic phenotype, OP and AFOP, while thrombus formation was considered present or absent.

Histomorphometry

Histochemistry and immunohistochemistry stains were quantified by morphometry. The images were captured with a digital camera on microscope (Novel L3000 LED) and analyzed using Image Pro Plus 7 software. Quantification followed the morphometric standards established by American Thoracic Society and European Thoracic Society (ATS/ERS)11. Positive α-sma cells were evaluated in ten different, randomly selected high-power fields of the lungs. At 400x magnification, the number of positive cells in each field was calculated according to the number of points hitting positive cells as a proportion of the total grid area. The density of collagen and elastic fibers was measured in the lung parenchyma in ten randomly selected microscopic fields at a magnification of ×200. The threshold for collagen fibers was established for all slides after the contrast was enhanced to the point at which the fibers were easily identified as green or orange bands. The density of the collagen and elastic fibers was expressed as the ratio between the measured fibers divided by the total area studied ×100.

Immunofluorescence Analysis

Neutrophils and neutrophil-derived extracellular traps (NETs) were identified and quantified by immunofluorescence with anti-myeloperoxidase (MPO) (Abcam, ab25989), anti-citrullinated histone H3 (H3Cit) (Abcam, ab5103) and DAPI (Life Technologies, D1306). The staining and quantification were performed as our lab protocol12.

Western Blotting Analysis

Protein extraction were performed for all fresh lung samples using QIAzol Lysis Reagent (Qiagen), according to our standered protocol. Protein levels were measured with the Bio-Rad Protein Assay (cat. 5000006) and transferred into a nitrocellulose membrane (GE Healthcare Biosciences, cat. 10600002). The membranes were incubated with anti-α-SMA (Abcam, ab7817), anti-matrix metalloproteinase-2 (anti-MMP-2) (Abcam, ab97779) and anti-glyceraldehyde-3-phosphate dehydrogenase (anti-GAPDH) (Abcam, ab9485). The proteins were detected and quantified with ChemiDoc Imaging Systems (Bio-Rad) and the results were normalized by GAPDH expression.

Real-Time Polymerase Chain Reaction for Viral RNA

Total RNA from fresh lung tissue of all patients was obtained using QIAzol Lysis Reagent (Qiagen). Reverse transcriptase was performed using SuperScript IV reverse transcriptase (Invitrogen, Carlsbad, USA) according to manufacturer’s instructions and our standard protocol. RT-PCR for detection of SARS-CoV-2 nucleocapsid (N2) and envelope (E) viral proteins was performed on Step One Plus Real-Time PCR System (Applied Biosystems,Grand Island, NY) with TaqMan probes (Applied Biosystems). Both primers were designed by ThermoFisher Scientific from N2 and E genes sequences developed by Centers for Disease Control and Prevention (CDC) of Atlanta (USA)13 and Charité Germany14, respectively. Samples were tested in triplicate and cycle threshold values less than 44 were considered positive.

Statistical Analysis

Statistical analysis was performed with SPSS v.13.0 0 software (SPSS, Inc., Chicago, IL, 2004). The data was evaluated using T test. Furthermore, nonparametric correlation (Spearman) was also performed with all collected data. Linear regression was also performed for these patients with all available PaO2/FiO2 values. Data are expressed as mean ± standard deviation and the p value less than 0.05 was be considered statistically significant.

RESULTS

COVID-19 Autopsy Cases: Demographic, Clinical, and Imaging Data

47 COVID-19 cases were confirmed by RT-PCR from nasopharyngeal swab. Following diagnosis, we conducted minimally invasive autopsies (MIA) on all positive patients at our institution (HCFMRP/USP). Demographic and clinical data for the 47 COVID-19 patients are summarized in Table 1.

View this table:
  • View inline
  • View popup
Table 1 – COVID-19 Clinical Data.

The mean time from the onset of symptoms to death was 18.9 days with a mean hospital length of stay of 14.4 days. Patients were equally distributed by gender and the mean age was 67.8. Preexisting medical comorbidities were present in all patients, with systemic arterial hypertension as the most common (55%) followed by obesity (36%). The most frequent initial symptom was cough (55%). Dyspnea was present before the admission (66%), at the time of admission (79%), and during the hospital stay (40%). Clinical complications post hospital admission includes septic shock (62%), acute renal failure (51%) and acute respiratory distress syndrome (ARDS) (45%).

Retrospective review of chest X-ray (CXR) and high-resolution computed tomography (HRCT) of the chest was performed for all patients. On hospital admission, ground-glass opacities (GGO) were found in all 47 cases, and interstitial opacities were noted in 23 cases (48.9%). Consolidations were present in 13 of 15 cases (87%). Additional findings include traction bronchiectasis (N=6; 40%) and architectural distortion (N=8; 53.3%) (Supplementary Table 1).

Retrospective review of in-hospital treatments includes the use of steroids (81%) and prophylactic anticoagulation therapy (100%). Initial laboratory tests, collected in the first 48 hours post admission, demonstrated increased values of lactate dehydrogenase, C-reactive protein and D-dimer (Table 1). Additional laboratory tests are described in Supplementary Table 2. Mechanically ventilated patients (N=36, 76.6%), received neuromuscular blockade and prone positioning treatment (Table 1). PaO2/FiO2 ratio statistically correlated with D-dimer levels (r=0.37; P<0.05) (Supplementary Fig. 1A).

Covid-19 Bimodal Clinical and Pathological Phenotypes

The PaO2/FiO2 linear regression, based on daily measure from the last value recorded for PaO2/FiO2 (before death) to the onset of mechanical ventilation, showed a bimodal clinical-pathological phenotype (Table 2):

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2 - Bimodal Clinical and Pathological Phenotypes: Mechanical Ventilation Data.
  1. Fibrotic Phenotype (N=5) – characterized by progressive decline in PaO2/FiO2 ratio and low compliance levels during hospitalization (Figure 1A); and alveolar septal thickening with myxoid fibrosis typical to an organizing-phase of diffuse alveolar damage (DAD) (Figure 1C; Table 3). Compliance correlated negatively with drive pressure for fibrotic phenotype (r=-0.67; P<0.05) (Supplementary Fig. 1C).

  2. Thrombotic Phenotype (N=10) – characterized by a progressive increase in PaO2/FiO2 ratio and high pulmonary compliance levels during hospitalization (Figure 1B); and resolution of the acute/subacute lung injury to near to normal or slight distortion of the underlying lung parenchyma architecture (Figure 1D; Table 3), and with high frequency of thrombosis (Figure 1D).

Figure 1 –
  • Download figure
  • Open in new tab
Figure 1 – Bimodal Clinical and Histopathological Phenotypes.

PaO2/FiO2 linear regression is based on daily measurement from the last value recorded (before death) to the onset of mechanical ventilation. It correlates with distinct histopathological features showing a bimodal phenotype: 1) Fibrotic Phenotype (N=5) showing progressive decline in PaO2/FiO2 ratio (A), and significant alveolar septal thickening by fibrosis (C – double blue arrows, E) without thrombus formation (C); and 2) Thrombotic Phenotype (N=10) showing progressive increase in PaO2/FiO2 ratio (B), and resolution of acute/subacute lung injury to normal or near to normal parenchyma architecture and thrombus formation on vessels (D, E). A population showing a mixed phenotype of lung injury (N=31) with morphologic features from both fibrotic and thrombotic phenotypes, including but not limited to organizing pneumonia, acute and fibrinous organizing pneumonia, and variable amounts of interstitial scarring, was present. This population seems to represent an intermediate non-bimodal phenotype, showing involvement by different stages of acute, organizing, and fibrotic DAD, that appears to fall between two opposite ends of lung injury manifesting as fibrotic and thrombotic respectively. Interestingly, a moderate to low frequency of thrombosis was noticed in these patients (E).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 3 - Bimodal Clinical and Pathological Phenotypes: Histopathological Findings in Minimally Invasive Autopsy

Clinical and radiological data from both phenotypes are described in Supplementary Tables 4 to 6. D-dimer and platelets count was significantly higher in patients with thrombotic phenotype compared to fibrotic phenotype (Supplementary Fig. 2A, B and Supplementary Table 5; P <0.05). Findings from CRX image study showed similar distribution between bimodal phenotypes (Supplementary Table 6).

COVID-19 Mixed Phenotype

A detailed histomorphological review of 47 MIA was performed to better understand COVID-19 induced lung injury (Supplementary Table 3). Variable amounts of alveolar septal thickening by fibrosis (interstitial fibrosis) were identified as the most common feature affecting the majority of the patients (N=41; 87.2%) (Supplementary Table 3). These findings significantly correlated with PaO2/FiO2 ratios (r=-0.42; P<0.05) (Supplementary Fig. 1B). In our cohort of 47 patients, we observe a unique bimodal phenotype in 15 patients (32%). However, 32 patients (68%) revealed a mixed phenotype of lung injury showing morphologic features from both fibrotic and thrombotic phenotypes. This population seems to represent an intermediate non-bimodal phenotype that appears to fall between two opposite ends manifesting as fibrotic and thrombotic respectively (Figure 1E). They may represent an intermediate phenotype showing different stages of involvement acute/subacute lung injury and fibrotic DAD morphologic features including organizing pneumonia (OP), acute fibrinous and organizing pneumonia (AFOP) with intra-alveolar fibrin “balls”, and variable amounts of interstitial scarring. Interestingly, a moderate to low frequency of thrombosis was noticed in these patients (Figure 1E).

Pathological features of Covid-19 Fibrotic Phenotype

The fibrotic phenotype demonstrated in our cohort appears to be triggered by an imbalance between degradation and production of extracellular matrix (ECM) by myofibroblastic activation. Severe decrease in ECM elastic fibers was found in lungs from patients presenting the fibrotic phenotype, contrasting with the preserved elastic framework seen in lungs from patients with the thrombotic phenotype, as highlighted by Verhöeff staining (P<0.05) (Figure 2A, B). This is further supported by the presence of a direct association between the area fraction of elastic fibers in ECM and the compliance levels during mechanical ventilation (r=0.57; P<0.05) (Supplementary Fig. 1D). The fibrotic phenotype patients showed prominent accumulation of α-SMA in ECM, a marker of myofibroblast activity and neovasculogenesis, compared to thrombotic phenotype patients (P <0.05) (Figure 2A, C). There was also an inverse association between the area fraction of elastic fibers and myofibroblasts in ECM (r=-0.57; P<0.05) (Supplementary Fig. 1E). α-SMA and anti-MMP-2 protein expression by Western blot were increased in the fibrotic phenotype (P<0.05) (Figure 2D), reinforcing the fibrotic nature of this process.

Figure 2 –
  • Download figure
  • Open in new tab
Figure 2 – Bimodal Clinical and Pathological Phenotypes: Elastic Fibers Framework and Myofibroblastic Activation.

Fibrotic phenotype demonstrates higher levels of α-SMA expression (red arrows), indicating extracellular matrix deposition by myofibroblasts and neovasculogenesis, resulting in lung tissue scarring and stiffening (P <0.05) (A – red arrow, C). Consequently, lung elastic fibers related to tissue elastic capacity (% area fraction of elastolysis), was significantly reduced in fibrotic phenotype when compared to thrombotic phenotype (P <0.05) (A, B). The cicatricial nature of the process is further supported by Western blot, revealing higher protein expression levels of α-SMA and MMP-2 in fibrotic phenotype (D). GAPDH was used as loading control for gene expression.

Picrosirius Red staining demonstrated a significant increase in ECM deposition by green-birefringence thin fibers (most by types III collagen) and red/orange-birefringence thick fibers (most by type III collagen) in lungs from patients with fibrotic phenotype (P<0.05) (Figure. 3B, C). These findings significantly correlate with PaO2/FiO2 ratios (r=-0.64; P<0.05) (Supplementary Fig. 1F). Additionally, neutrophils with significantly higher length of NETs were also observed in these patients (P<0.05) (Supplementary Fig. 3A, C).

Figure 3 –
  • Download figure
  • Open in new tab
Figure 3 – Bimodal Clinical and Pathological Phenotypes: Collagen Depositions.

Extensive collagen production and interstitial deposition was noted, characterized by the presence of alveolar septal expansion by dense fibrosis (yellow arrows), which was more prominent in the fibrotic phenotype, confirmed by picrosirius red staining (A – yellow arrow). Area fraction percentage of collagen fibers labeled by green and yellow orange-polarized picrosirius red stain was significantly higher in fibrotic phenotype compared with thrombotic one, as demonstrated by P <0.05 (B, C).

Fibrotic Phenotype in Post-COVID-19 Survivors Follow Up

Post COVID-19 infection follow up with transbronchial biopsies were performed at our institution (HCFMRP-USP) in two patients, both in their 60’s and male, ranging from three- to five-months’ post infection. Histopathologic analysis revealed minimal chronic inflammation with remodeling of airway wall by submucosal scarring and alveolar septal thickening by fibrosis (Figure 4). While this may represent scarring secondary to resolving infection, the possibility of future post-infection chronic fibrosing disease needs to be excluded. Secondly, more extensive studies are needed to determine if post COVID-19 complications arise in patients who survive this disease.

Figure 4 –
  • Download figure
  • Open in new tab
Figure 4 – Post COVID 19 Infection Transbronchial Biopsies Follow Up.

Two patients, both in their 60’s and male, ranging from three- to five-months’ post infection: There is minimal chronic inflammation and airway wall and lung parenchyma remodeling, characterized by submucosal scarring and alveolar septal thickening by fibrosis associated with alveolar pneumocyte prominence, respectively. While this may represent scarring secondary to resolving infection, the possibility of an evolving post infection chronic fibrosing disease needs to be excluded and more extensive studies are needed to investigate post COVID-19 complications in patients who survive the disease.

DISCUSSION

We discovered two distinctive histopatholgical bimodal phentoypes and correlated histopathologic findings to corresponding clinical outcomes. Patients presenting with confirmed lung fibrosis predominantly displayed either thrombotic or fibrotic characterization, which correlated to differing clinical outcomes. This is in contrast to the initial pathophysiological events that occur in the acute phase of ARDS following post-viral insult to the lungs where the clinical outcomes diverge as infection progresses.

Upon viral insult to the airways and bronchial epithelial cells, a severe lymphocytic bronchiolitis occurs, promoting epithelial injury and cell death, eventually resulting in structural disarray by airway-centered remodeling (Figure 5). The aggressive airway wall infection then spreads to the lung parenchyma, inducing cellular pneumonitis, and damages the alveolar-capillary barrier (Figure 5, yellow arrow). This promotes the denudation of the alveolar septa, characterized as DAD, bronchiolocentric alveolar hemorrhage or fibroblastic plug, called OP, associated or not with fibrin balls, known as AFOP (Figure 5). Subsequently, two different outcomes of repair appear to occur depending on resolution or progression of the lung injury, eventually manifesting as a thrombotic phenotype with gradual recovering of the acute/subacute lung injury associated with thrombosis; or a fibrotic phenotype with sustained myofibroblastic proliferation with interstitial scarring and parenchymal remodeling (Figure 1, 5).

Figure 5 –
  • Download figure
  • Open in new tab
Figure 5 – COVID-19 Lung Pathophysiology.

Following bronchial wall viral infection there is severe lymphocytic bronchiolitis, promoting striking epithelial injury, cell loss, and mucosal ulceration, eventually resulting in airway remodeling. The aggressive airway wall infection then spreads to lung parenchyma, inducing cellular pneumonitis and damage of alveolar-capillary barrier (immunohistochemical stain decorates SARS-CoV-2 infected cells, yellow arrow). This promotes bronchiolocentric alveolar hemorrhage, acute DAD, and fibrin balls associated with fibroblastic plugs of organizing pneumonia, known as acute fibrinous and organizing pneumonia. Subsequently, two different ultimate ends of repair appear to occur depending on resolution or progression of the lung injury, eventually manifesting as a thrombotic phenotype, with gradual recovering of the acute/subacute lung injury associated with prominent thrombosis; or a fibrotic phenotype with sustained myofibroblastic proliferation, which promotes interstitial scarring and parenchymal remodeling.

The fibrotic phenotype seen in some patients is the end are the result of the activation and termination of many different pathways, one of which is collagen fibers deposition (type I and III) to cause alveolar septal thickening (Figure 3). Additionally, the fibrotic phenotype revealed higher active myofibroblasts α-SMA area fraction (Figure 2A, C) and higher α-SMA and MMP-2 expression (Figure 2D) (P <0.05) than the thrombotic phenotype. Similar studies also describes COVID-19 organizing DAD with α-SMA+ myofibroblasts15 and fibrous parenchyma remodeling16,17. Alveolar septal thickening by fibrosis impairs adequate lung function and gas exchange, which is reflected in our study by significant correlation between collagen fibers deposition and PaO2/FiO2 ratio in the fibrotic phenotype (r=-0.64; P <0.05) (Supplementary Fig. 1F). Furthermore, lung elastic capacity, represented by elastic fibers area fraction (Verhoeff stain quatification) was significantly reduced (P<0.05) (Figure 2A, B). Moreover, fibrotic phenotype patients presented NETs in a significantly higher extent throughout lung parenchyma compared to thrombotic phenotype subjects (P<0.05) (Supplementary Fig. 3). As described by our group previously12, the formation of neutrophil extracelular traps (NETs), networks of extracellular fibers, primarily composed of DNA from neutrophils, binded with pathogens 18, is an important response to SARS-CoV-2 during the acute and a potential contributory factor to later fibrotic phase by amplify the chronic reparative phase. This result suggests a more prominent host inflammatory response in the fibrotic phenotype subjects with consequently greater damage to lung parenchyma, trigging a more aggressive reparative fibroblastic process, as observed in these patients. In accordance with the histological analysis, the clinical data also suggested a poor clinical outcome. These patients presented progressive decline in PaO2/FiO2 ratio (Figure 1A) (mean 101.16) and low compliance levels during hospitalization (mean 19.8 ml/cmH2O) (Table 2). In accordance to other studies, low ventilatory parameters were also observed19,20 and patients presented a median value of 129 PaO2/FiO2 ratio and 27 ml/cmH2O compliance, as described by Cummings M et al21,

Conversely, subjects from the thrombotic phenotype showed an unexpected progressive increase in PaO2/FiO2 ratio (Figure 1B) (mean 141) and high pulmonary compliance levels (mean 41.4 ml/cmH2O) (Table 2). Similarly, a nearly normal lung compliance has also surprised Gattinoni et al22 as their patients presented 50.2 ml/cmH2O, despite the severe hypoxaemia. However, better ventilatory parameters in our thrombotic phenotype subjects could be justified by the gradual recovering of the acute/subacute lung injury (Figure 1D) with increased area fraction of elastic fibers (Figure 2A, B). Nevertheless, sudden death occurred, probably due to higher frequency of pulmonary thromboembolism (Figure 1D, Table 3), as also described in other studies23,24. As expected, D-dimer levels and platelets count were significantly higher in these patients than in fibrotic phenotype subjects (Supplementary Fig. 2A, B and Supplementary Table 5). High D-dimer levels in COVID-19 patients has been recently related with thrombosis as an important prognostic biomarker of severity 25,26. In addition, platelets have been reported with a hyperactivated phenotype in COVID-19 patients, possibly related to hypercoagulation8. While it’s well known that thrombus formation is a feature of organizing DAD, it’s unclear why some COVID-19 patients develop a more frequently form thrombi.

In light of this, we observed two distinct bimodal histopathological phenotypes with their correlative clinical outcomes among the 15 COVID-19 (32%) evidenced with fibrosis. This may justify the differing clinical outcomes observed in out cohort, despite the same etiology and the convergence to their endpoints- mortality. Likewise, Gattinoni et al3 have also reported different behaviors throughout COVID-19 patients’ outcome and classified them into two phenotypes: Type L: low elastance; and Type H: high elastance. Thus, COVID-19 patients with different clinical presentations are frequently described. However, our work provides insight into the importance of histopathological screening among COVID19 autopsies correlated to their differing clinical outcomes, with a particular focus in the lung repair resolving phase of injury.

Besides those two crucial phenotypes, it is important to mention that 32 patients (68%) had a mixed pattern of lung injury appearing as a non-bimodal phenotype (Figure 1E). The pathophysiology of these patients probably represents different stages of viral lung injury, with both features of acute/subacute and fibrotic/organizing lung injury with variable degrees of interstitial scarring, OP and AFOP. Despite the importance of this histopathological finding, the correlation with clinical data did not highlight any relevance. On the other hand, the clinical and pathological correlation of our 15 patients evidenced two distinctive bimodal opposite COVID-19 phenotypes, manifesting as fibrotic predominant and thrombotic predominant, being a unique and important finding.

Therefore, the reported sudden deaths of COVID-19 patients in clinical improvement and the description of post-COVID-19 tomographic fibrotic changes may be justified by our findings. In fact, two COVID-19 patients’ follow up with transbronchial biopsies, performed at 3 to 5 months post infection, demonstrated lung parenchyma remodeling with alveolar septal thickening by dense fibrosis (Figure 4). Pathophysiological differences in the virus versus immunity interaction may cause different clinical outcomes in COVID-19 patients, as demonstrated in our study after 47 COVID-19 autopsies. Our described bimodal clinical and pathological phenotypes are critically important in the management of infected patients, given that a progressive drop in PaO2/FiO2 ratio or high platelets/D-Dimer levels are the clinical findings that may raise the clinical suspicion for one of the two outcome phenotypes at the beginning of the hospitalization, allowing a possible better therapeutic management.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request. Source data are provided with this paper.

CONFLICT OF INTEREST

The authors declare no competing financial interests.

STUDY APPROVAL

The procedures followed in the study were approved by the National Ethics Committee – Brazil (CAAE: 32475220.5.0000.5440). The written informed consent was waived.

ACKNOWLEDGMENTS

The authors would like to thank the institutional support of the University Hospital of Ribeirão Preto Medical School, University of São Paulo, Ribeirão Preto, SP, Brazil – HCFMRP/USP. This work was supported by São Paulo Research Foundation (FAPESP) [2019/01517-3, 2019/19591-5 and 2020/13370-4].

REFERENCES

  1. 1.↵
    Wu C, Chen X, Cai Y, et al. Risk Factors Associated With Acute Respiratory Distress Syndrome and Death in Patients With Coronavirus Disease 2019 Pneumonia in Wuhan, China. JAMA internal medicine. Jul 1 2020;180(7):934–943.
    OpenUrl
  2. 2.↵
    Du Y, Tu L, Zhu P, et al. Clinical Features of 85 Fatal Cases of COVID-19 from Wuhan. A Retrospective Observational Study. Am J Respir Crit Care Med. Jun 1 2020;201(11):1372–1379.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Gattinoni L, Chiumello D, Caironi P, et al. COVID-19 pneumonia: different respiratory treatments for different phenotypes? Intensive care medicine. Jun 2020;46(6):1099–1102.
    OpenUrlPubMed
  4. 4.↵
    De Michele S, Sun Y, Yilmaz MM, et al. Forty Postmortem Examinations in COVID-19 Patients. American journal of clinical pathology. Nov 4 2020;154(6):748–760.
    OpenUrlPubMed
  5. 5.↵
    Severe Covid GG, Ellinghaus D, Degenhardt F, et al. Genomewide Association Study of Severe Covid-19 with Respiratory Failure. The New England journal of medicine. Oct 15 2020;383(16):1522–1534.
    OpenUrlCrossRefPubMed
  6. 6.↵
    Lee LYW, Cazier JB, Starkey T, et al. COVID-19 prevalence and mortality in patients with cancer and the effect of primary tumour subtype and patient demographics: a prospective cohort study. The Lancet. Oncology. Oct 2020;21(10):1309–1316.
    OpenUrlCrossRefPubMed
  7. 7.↵
    Amendola A, Garoffolo G, Songia P, et al. Human cardiosphere-derived stromal cells exposed to SARS-CoV-2 evolve into hyper-inflammatory/pro-fibrotic phenotype and produce infective viral particles depending on the levels of ACE2 receptor expression. Cardiovascular research. May 25 2021;117(6):1557–1566.
    OpenUrl
  8. 8.↵
    Bongiovanni D, Klug M, Lazareva O, et al. SARS-CoV-2 infection is associated with a pro-thrombotic platelet phenotype. Cell death & disease. Jan 5 2021;12(1):50.
    OpenUrlCrossRef
  9. 9.↵
    Weustink AC, Hunink MM, Van Dijke CF, Renken NS, Krestin GP, Oosterhuis JW. Minimally invasive autopsy: an alternative to conventional autopsy? Radiology. 2009;250(3):897–904.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    Pontelli MC, Castro IA, Martins RB, et al. Infection of human lymphomononuclear cells by SARS-CoV-2. bioRxiv. 2020:2020.2007.2028.225912.
  11. 11.↵
    Hsia CC, Hyde DM, Ochs M, Weibel ER, Structure AEJTFoQAoL. An official research policy statement of the American Thoracic Society/European Respiratory Society: standards for quantitative assessment of lung structure. Am J Respir Crit Care Med. Feb 15 2010;181(4):394–418.
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    Veras FP, Pontelli M, Silva C, et al. SARS-CoV-2 triggered neutrophil extracellular traps (NETs) mediate COVID-19 pathology. medRxiv. 2020:2020.2006.2008.20125823.
  13. 13.↵
    Control CfD, Prevention. CDC 2019-novel coronavirus (2019-nCoV) real-time RT-PCR diagnostic panel. 2020.
  14. 14.↵
    Corman VM, Landt O, Kaiser M, et al. Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin. Jan 2020;25(3).
  15. 15.↵
    Chen J, Wu H, Yu Y, Tang N. Pulmonary alveolar regeneration in adult COVID-19 patients. Cell research. Aug 2020;30(8):708–710.
    OpenUrlCrossRefPubMed
  16. 16.↵
    Bradley BT, Maioli H, Johnston R, et al. Histopathology and ultrastructural findings of fatal COVID-19 infections in Washington State: a case series. Lancet. Aug 1 2020;396(10247):320–332.
    OpenUrlCrossRefPubMed
  17. 17.↵
    Grillo F, Barisione E, Ball L, Mastracci L, Fiocca R. Lung fibrosis: an undervalued finding in COVID-19 pathological series. The Lancet. Infectious diseases. Jul 28 2020.
  18. 18.↵
    Brinkmann V, Reichard U, Goosmann C, et al. Neutrophil extracellular traps kill bacteria. Science. Mar 5 2004;303(5663):1532–1535.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    Schenck EJ, Hoffman K, Goyal P, et al. Respiratory Mechanics and Gas Exchange in COVID-19-associated Respiratory Failure. Ann Am Thorac Soc. Sep 2020;17(9):1158–1161.
    OpenUrlCrossRefPubMed
  20. 20.↵
    Ziehr DR, Alladina J, Petri CR, et al. Respiratory Pathophysiology of Mechanically Ventilated Patients with COVID-19: A Cohort Study. Am J Respir Crit Care Med. Jun 15 2020;201(12):1560–1564.
    OpenUrlCrossRefPubMed
  21. 21.↵
    Cummings MJ, Baldwin MR, Abrams D, et al. Epidemiology, clinical course, and outcomes of critically ill adults with COVID-19 in New York City: a prospective cohort study. medRxiv. Apr 20 2020.
  22. 22.↵
    Gattinoni L, Coppola S, Cressoni M, Busana M, Rossi S, Chiumello D. COVID-19 Does Not Lead to a “Typical” Acute Respiratory Distress Syndrome. Am J Respir Crit Care Med. May 15 2020;201(10):1299–1300.
    OpenUrlCrossRefPubMed
  23. 23.↵
    Del Nonno F, Colombo D, Nardacci R, Falasca L. Fatal pulmonary arterial thrombosis in a COVID-19 patient, with asymptomatic history, occurred after swab negativization. Thrombosis journal. Jan 6 2021;19(1):1.
    OpenUrl
  24. 24.↵
    Polat V, Bostanci GI. Sudden death due to acute pulmonary embolism in a young woman with COVID-19. Journal of thrombosis and thrombolysis. Jul 2020;50(1):239–241.
    OpenUrl
  25. 25.↵
    Zhou F, Yu T, Du R, 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.
    OpenUrlCrossRefPubMed
  26. 26.↵
    Huang C, Wang Y, Li X, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. Feb 15 2020;395(10223):497–506.
    OpenUrlCrossRefPubMed
Back to top
PreviousNext
Posted September 14, 2021.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
COVID-19 Bimodal Clinical and Pathological Phenotypes
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
COVID-19 Bimodal Clinical and Pathological Phenotypes
Sabrina S Batah, Maíra N Benatti, Li Siyuan, Wagner M Telini, Jamile Barbosa, Marcelo B Menezes, Tales R Nadai, Keyla S G Sá, Chirag M. Vaswani, Sahil Gupta, Dario S Zamboni, Danilo T Wada, Rodrigo T Calado, Renê D R Oliveira, Paulo Louzada-Junior, Maria Auxiliadora-Martins, Flávio P Veras, Larissa D Cunha, Thiago M Cunha, Rodrigo Luppino-Assad, Marcelo L Balancin, Sirlei S Morais, Ronaldo B Martins, Eurico Arruda, Fernando Chahud, Marcel Koenigkam-Santos, Andrea A Cetlin, Fernando Q Cunha, Claudia dos Santos, Vera L Capelozzi, Junya Fukuoka, Rosane Duarte-Achcar, Alexandre T Fabro
medRxiv 2021.09.03.21262841; doi: https://doi.org/10.1101/2021.09.03.21262841
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
COVID-19 Bimodal Clinical and Pathological Phenotypes
Sabrina S Batah, Maíra N Benatti, Li Siyuan, Wagner M Telini, Jamile Barbosa, Marcelo B Menezes, Tales R Nadai, Keyla S G Sá, Chirag M. Vaswani, Sahil Gupta, Dario S Zamboni, Danilo T Wada, Rodrigo T Calado, Renê D R Oliveira, Paulo Louzada-Junior, Maria Auxiliadora-Martins, Flávio P Veras, Larissa D Cunha, Thiago M Cunha, Rodrigo Luppino-Assad, Marcelo L Balancin, Sirlei S Morais, Ronaldo B Martins, Eurico Arruda, Fernando Chahud, Marcel Koenigkam-Santos, Andrea A Cetlin, Fernando Q Cunha, Claudia dos Santos, Vera L Capelozzi, Junya Fukuoka, Rosane Duarte-Achcar, Alexandre T Fabro
medRxiv 2021.09.03.21262841; doi: https://doi.org/10.1101/2021.09.03.21262841

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Respiratory Medicine
Subject Areas
All Articles
  • Addiction Medicine (349)
  • Allergy and Immunology (668)
  • Allergy and Immunology (668)
  • Anesthesia (181)
  • Cardiovascular Medicine (2648)
  • Dentistry and Oral Medicine (316)
  • Dermatology (223)
  • Emergency Medicine (399)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (942)
  • Epidemiology (12228)
  • Forensic Medicine (10)
  • Gastroenterology (759)
  • Genetic and Genomic Medicine (4103)
  • Geriatric Medicine (387)
  • Health Economics (680)
  • Health Informatics (2657)
  • Health Policy (1005)
  • Health Systems and Quality Improvement (985)
  • Hematology (363)
  • HIV/AIDS (851)
  • Infectious Diseases (except HIV/AIDS) (13695)
  • Intensive Care and Critical Care Medicine (797)
  • Medical Education (399)
  • Medical Ethics (109)
  • Nephrology (436)
  • Neurology (3882)
  • Nursing (209)
  • Nutrition (577)
  • Obstetrics and Gynecology (739)
  • Occupational and Environmental Health (695)
  • Oncology (2030)
  • Ophthalmology (585)
  • Orthopedics (240)
  • Otolaryngology (306)
  • Pain Medicine (250)
  • Palliative Medicine (75)
  • Pathology (473)
  • Pediatrics (1115)
  • Pharmacology and Therapeutics (466)
  • Primary Care Research (452)
  • Psychiatry and Clinical Psychology (3432)
  • Public and Global Health (6527)
  • Radiology and Imaging (1403)
  • Rehabilitation Medicine and Physical Therapy (814)
  • Respiratory Medicine (871)
  • Rheumatology (409)
  • Sexual and Reproductive Health (410)
  • Sports Medicine (342)
  • Surgery (448)
  • Toxicology (53)
  • Transplantation (185)
  • Urology (165)