Differentially expressed genes in Diffuse Alveolar Damage (DAD) patterns of COVID-19

ABSTRACT

Severe COVID-19 induces DAD, a condition with temporal-spatial heterogeneity. We determined the differentially expressed genes (DEGs) in the histological patterns of DAD. Twelve fatal COVID-19 cases were classified in acute DAD (n=5) and intermediate/advanced (IA) DAD (n=7). Autopsy lung RNA was extracted from COVID-19 and 4 control cases. RNA sequencing was performed on the Illumina NovaSeq 6000. Enrichment analysis was performed with clusterProfiler using Genome-wide annotation for Human R package. GO terms and KEGG pathways were considered enriched if adjusted p<=0.05. Principal component analysis showed that IA-DAD samples were grouped, while acute DAD samples were scattered. The differential expression analysis between these two groups and the control cases revealed 261 DEGs in the acute DAD (143 Up- and 53 Down-regulated), 244 DEGs in the IA- DAD tissues (67 Up- and 116 Down-regulated), and 61 DEGs were shared between them (45 Up- and 16 Downregulated). Patients with acute DAD had up-regulated genes related to oxidative phosphorylation, blood coagulation, megakaryocytes differentiation/regulation, and platelet degranulation/activation. Patients with IA-DAD had DEGs related to immunoglobulins and extracellular matrix. The shared up-regulated DEGs between both patterns are involved in innate and adaptive immune responses. We selected 3 DEGs in each DAD pattern for validation by realtime PCR. There were no differences in acute DAD DEGs, but DEGs overexpressed in intermediate DAD (COL3A1, IGLV3-19, IGHV1-58) were significantly higher. Genes related to thrombotic events occur at the acute stage of DAD, whereas immunoglobulin production and remodeling occur at later stages of DAD.