The transcriptional response of human blood leukocytes to SARS-CoV-2 infection was investigated focusing on the differences between mild and severe cases and between age subgroups. Weighted gene co-expression network analysis and comparative gene expression analysis were used. Three transcriptional modules positively associated with the traits of interest and their respective high hierarchy genes were identified. Enrichment analyses showed that the yellow module, associated with severe cases and older patients, had an overrepresentation of genes involved in inflammatory and innate immune responses, and neutrophil activation. The magenta and black modules, associated with disease severity and younger patients, contained genes related to innate immunity and inflammation and genes that regulate those responses. Subnetworks for these modules were constructed using genes enriched for innate immunity, inflammation, immunoregulation and differentially expressed genes (severe vs. mild). Their analysis evidenced that immunoregulatory functions are more activated in the modules associated with younger patients, what may help to explain the better disease course and faster recovery observed in younger COVID-19 patients. Comparative gene expression analysis between severe and mild groups, followed by gene enrichment and normalized gene expression analyses, revelated a set of 23 potential biomarkers for COVID-19 severity, of which 13 are newly described.