Clinical and Genetic Characteristics of Coronaviruses with Particular Emphasis on SARS-CoV-2 Virus

Coronaviruses (CoVs) are unsegmented single-stranded RNA viruses, belonging to the subfamily Coronavirinae, family Coronaviridae, and order Nidovirales (Weiss and Leibowitz 2011). The structural differences of their domains became the basis for the identification of four groups: α, β, γ, and δ (Li 2016). For a long time, coronaviruses were recognized as pathogens not posing a significant threat to humans. In most cases, they caused mild infections, similar to the common cold. However, the occurrence of the etiological factor of a severe acute respiratory syndrome (SARS-CoV) in 2003 and the Middle East respiratory syndrome (MERS-CoV) in 2012, drew the world’s attention to their species barrier-crossing capacity and the destructive potential of human CoVs (Song et al. 2019). It is now known that infections caused by coronaviruses can cause severe dysfunctions of the respiratory, gastrointestinal, or central nervous systems (Perlman and Netland 2009). In late December 2019, the first cases of illness caused by SARS-CoV-2 were reported in Wuhan in China’s Hubei province. The rapidly evolving Coronavirus Disease 2019 (COVID-19) has led to a global healthcare crisis. Initially, the adverse effects of COVID-19 were thought to occur mainly in the respiratory tract, causing pneumonia and acute respiratory distress syndrome (ARDS). It is well known that SARS-CoV-2 infection can lead to extensive thrombotic lesions with microangiopathy, multi-organ failure, and death (Ackermann et al. 2020).

SARS-CoV-2 belongs to the group of zoonotic viruses that can infect both humans and animals. Genetic analysis showed that the new pathogen belongs to the subgenus Sarbecovirus. It is more similar to two bat-derived strains (SL-CoVZC45 and SL-CoVZXC21) than to known coronaviruses, including the virus from the SARS outbreak in 2003 (Lu et al. 2020). The recent publications demonstrated the importance of individual variability in response to SARS-CoV-2 infection. Individual body predisposition, modulated by genetic background, may determine the course of COVID-19 and further prognosis. Despite performing analyses, knowledge of the influence of genetic factors on the interindividual variability of patients with COVID-19 is insufficient to establish fully effective diagnostic and therapeutic interventions. Therefore, a better understanding of this topic is a priority for many researchers worldwide.

The clinical part of the manuscript (risk factors, clinical course, diagnosis, treatment, and prevention) was written in September 2021, with an update in February 2022. Demographic risk factors for disease progression include male gender, age over 65, and smoking. Moreover, comorbidities such as diabetes, cardiovascular diseases, arterial hypertension, malignant neoplasms, respiratory diseases (mainly chronic obstructive pulmonary disease), chronic kidney disease, and liver diseases increase the risk of severe COVID-19 and death (Zheng et al. 2020; Zhou et al. 2020). Laboratory indicators of possible severe disease are thrombocytopenia, elevated IL-6 levels, Ast > 40 U/l, creatinine ≥ 133 mol/l, troponin I (hs-cTnI) > 28 pg/ ml, procalcitonin (PCT) > 0.5 ng/ ml, lactate dehydrogenase (LDH) > 245 U/l, and D-dimers > 0.5 mg/l. The manifestation of symptoms such as dyspnea, weakness, and sputum production has also been associated with the probable fatal course of COVID-19 (Zhou et al. 2020).

In mechanically ventilated patients and those connected to ECMO, the prognosis is poor, and high mortality is observed in these groups. Hospital mortality of patients connected to ECMO due to COVID-19 is 37% and is similar to ARDS caused by diseases other than COVID-19 (Barbaro et al. 2020; Ramanathan et al. 2021). In the USA, the 28-day mortality after the ICU (intensive care unit) admission was defined as 35%. In such conditions, the following factors were considered independent risk factors for death: age over 80, male sex, BMI over 40, ischemic heart disease, active neoplastic disease, PaO₂ : FiO₂ below 100, and liver and kidney failure (Gupta et al. 2020).

SARS-CoV-2 is transmitted mainly by droplets during face-to-face exposure (by sneezing or coughing); however, transmission is possible through direct contact with the infected person or with contaminated surfaces, although it is a marginal transmission route (Wiersinga et al. 2020; Ravindra et al. 2022). In addition, the possibility of fecal-oral transmission has been suggested, which may be evidenced by the presence of SARS-CoV-2 RNA in rectal swabs or stools (assayed by RT-PCR) (Cheung et al. 2020; Bwire et al. 2021).

The disease can take a diverse clinical course. It is often asymptomatic, or it runs as a mild respiratory infection. However, some patients will develop severe pneumonia with cytokine release syndrome (CRS) and progression to respiratory failure (ARDS). Mortality in Poland was estimated at 2.106% on the 4^th of February 2022 (The Johns Hopkins Coronavirus Resource Center 2022).

The incubation period ranges from two to 14 days. The median incubation period is 5.1 days, while 97.5% of all infections will develop within 11.5 days. The incubation period of SARS-CoV-2 is similar to that of other highly pathogenic coronaviruses: SARS (median five days, range 2–14) and MERS (median 5–7 days, range 2–14 days) (Lauer et al. 2020). The SARS-CoV-2 infection has a frequent oligosymptomatic course, including symptoms such as fever, chills, cough, shortness of breath, difficulty breathing, fatigue, headache, throat, muscle or body aches, runny nose, symptoms of conjunctivitis, anorexia, nausea, vomiting, diarrhea, abdominal pain, and loss of sense of smell (anosmia) or taste (ageusia) (Wiersinga et al. 2020). The last two symptoms are frequent in women, children, adolescents, and young adults (Lechien et al. 2020). In the elderly pulmonary manifestations may be preceded by atypical symptoms such as diarrhea with dehydration, hypotonia, cognitive dysfunction, delirium, falls, and body temperature may fluctuate daily with periods of hypothermia (Blain et al. 2020).

As the disease progresses, symptoms of pneumonia and increasing respiratory failure become crucial. Initially, the inflammation is interstitial, but it changes to a mixed form with bacterial superinfections. Dyspnea became the leading symptom, and the chest CT showed several changes: ground-glass opacity (GGO), crazy-paving pattern, and consolidation. The changes mainly affect the lower lobes of the lungs (Alsharif and Quarashi 2021).

There may also be diffuse alveolar damage, edema of the inter-alveolar spaces, protein-rich exudate into the lumen of the alveoli, inflammation, pneumocytes necrosis, formation of hyaline membranes, and fibrosis. Vascular changes may occur as microangiopathy and macroangiopathy (blood clots and hemorrhagic foci) (Englisch et al. 2020). Clinical prognosis depends mainly on the patient’s age and comorbidities (Zheng et al. 2020), while radiographic prognostic criteria assess to what extent the lungs are involved. A generalized inflammatory response of the organism, the so-called cytokine storm, may develop, and respiratory failure may progress to a multi-organ failure. Septic complications may occur (Koçak Tufan et al. 2021). Also, thromboembolic complications can be found, especially in those requiring hospitalization. Venous thromboembolism (VTE) occurs in 26%, pulmonary embolism (PE) in 12%, and deep vein thrombosis (DVT) alone in 14% of the patients. Among those requiring treatment in ICUs, the risk of complications is as follows: VTE 24%, PE 19%, and DVT 7%, and of note is the significantly higher incidence of pulmonary embolism. Independently of the hospital ward, the following vessels may be involved in pulmonary embolism: main trunk and lobar arteries in 37.8%, segmental arteries in 37.9%, and subsegmental arteries in 19% of people (Porfidia et al. 2020).

The outcome of the disease is a resolution of inflammatory changes with gradual recovery or death of the patient. After a severe course of the disease, the patient may require lengthy rehabilitation and may never return to the function and respiratory capacity before the disease.

In December 2020, the Alpha variant (B.1.1.7) was identified in the United Kingdom as the first SARS-CoV-2 VOC (Variant of Concern) (Rambaut et al. 2020). Like the second, the Beta variant (B.1.351) was detected in South Africa in 2020. The Beta variant has three mutations, including K417N, E484K, and N501Y in the Spike protein receptor-binding domain (Tegally et al. 2020). In January 2021, the Gamma variant (P.1) was reported in Brazil. The Gamma variant carries three mutations: K417T, E484K, and N501Y in the Spike protein receptor-binding domain (Faria et al. 2021). In December 2020, the Delta variant (B.1.617.2) carrying mutations: 452R and 478K was first reported in India. The Omicron variant (B.1.1.529) was first described in South Africa in November 2021 (WHO 2022). All VOCs are more transmissible than the wild-type virus. Alpha, Beta, Gamma, and Delta variants cause more severe illnesses than wild-type viruses regarding hospitalization, ICU admission, and mortality. The Beta and Delta variants are at higher risk to patients than the Alpha and Gamma variants (Lin et al. 2021).

To identify SARS-CoV-2, the real-time RT-PCR (real time Reverse Transcription Polymerase Chain Reaction) technique may be used and is considered a reference method. Real time RT-PCR detects regions of the SARS-CoV-2 genome encoding the following proteins: ORF1a/b (Open Reading Frame 1a/b), ORF1a, E (envelope), N (nucleocapsid), and S (Spike). Detectable RdRp (RNA-dependent RNA polymerase) is a part of ORF1a/b (Ravi et al. 2020; WHO 2020; Yüce et al. 2021). The specimens for the assay may include nasal, nasopharyngeal, pharyngeal, sputum, bronchoalveolar lavage fluid (BALF), blood, stool, and rectal swabs. Real time RT-PCR is more sensitive to virus detection when testing biological material from the lower respiratory tract than the upper respiratory tract.

The overall real time RT-PCR sensitivity is high (89.1%), as well as its specificity (98.9%) (Mustafa Hellou et al. 2021). Nasopharyngeal swabs are the most often tested. Real time RT-PCR from bronchial lavage fluid (BALF) has the highest sensitivity (the positive detection rate of 91.8%). Any SARS-CoV-2 RNA molecules were detected in urine or genital swabs. Although real time RT-PCR of bronchial lavage fluid has the highest sensitivity, the fluid is collected by bronchoscopy. It is an invasive method, almost impossible to perform in a patient with severe dyspnea and, therefore, not applicable in daily practice (Böger et al. 2021; Bwire et al. 2021). The real time RT-PCR from rectal swabs has a high positive detection rate and detects the presence of SARS-CoV-2 RNA in stools (Cheung et al. 2020). The timing of material collection for molecular testing is vital in the diagnosis. Clinical materials should be collected while the virus is replicating in the epithelium of the upper respiratory tract. Typically, viral replication continues until ten days after the onset of symptoms. After this time, the sensitivity of molecular tests decreases significantly. The most crucial technique remains real-time signal detection. PCR assay should be repeated when a specific clinical picture (or radiological changes) occurs, and the initial test brings a negative result. Real time RT-PCR is characterized by higher sensitivity than antigen tests (89.1% vs 73,8%).

Serological assays detect and measure the titers of antibodies against S (Spike), N (Nucleocapsid), and RBD (receptor binding domain) antigens of SARS-CoV-2. Therefore, the S and N proteins are the most immunogenic and useful antigens in serological diagnostics (Makoah et al. 2021; Ong et al. 2021). The production of antibodies is a consequence of COVID-19 or vaccination against it. The are several methods used for antibody detection: enzyme-linked immunoassays (ELISAs), chemiluminescent immunoassays (CLIAs), rapid diagnostic tests (RDTs), and neutralization assays. RDTs include lateral flow tests (LFTs) (Ravi et al. 2020). CLIAs have a higher sensitivity and range compared to ELISAs. RDTs are characterized by the lowest sensitivity and the shortest assay period (the results obtained after approx. 15 minutes). In comparison, the results of ELISAs and CLIAs are provided within 30–100 minutes (Ravi et al. 2020). IgA, IgM, IgG antibodies, or a mixture of the above classes can be detected and presented as the so-called total score (Ong et al. 2021). The earliest IgA seroconversion can be found two days after the first clinical symptoms, while IgM seroconversion – in about three to five days from the first symptoms (Yu et al. 2020). The peak of antibody production is in the second and third weeks of the disease. Higher titers of IgA compared to IgM antibodies have been found (Padoan et al. 2020). IgG antibodies become detectable seven to 14 days after the first symptoms, peak in the third and fourth weeks, and remain at a high level until the sixth week of the disease (Guo et al. 2020; Makoah et al. 2021). Seroconversion may be absent or delayed in the asymptomatic group, and specific antibody titers lower than in symptomatic patients. According to metanalysis, the specificity of serological assays is 95–99% (Lisboa et al. 2021). Despite the ongoing SARS-CoV-2 infection, the results may be negative when congenital and acquired humoral deficiencies coexist (e.g., people with AIDS or agammaglobulinemia). Patients may be in the serological window when, despite being infected with SARS-CoV-2, they have not yet developed an immune response, and antibodies concentration is too low. Serological tests are characterized by low sensitivity in the early phase of COVID-19 (first seven to 10 days); thus, they cannot replace PCR during this infection period (Makoah et al. 2021). Serological assays may be used in epidemiological studies, evaluation of the response to a vaccine, and diagnostics of post-inflammatory syndromes. They may be helpful in the diagnosis of COVID-19 plausible in people with a typical, long-lasting clinical picture and/or radiographic image of the lungs but with negative RT-PCR results (WHO 2020; Flisiak et al. 2021).

Antigen tests detect SARS-CoV-2 proteins (most often nucleocapsid) based on the LFIAs (lateral flow immunoassays) technique. The most frequently tested material is a nasopharyngeal swab (Aguilar-Shea et al. 2021). Tests available in healthcare facilities are characterized by high sensitivity and specificity. According to the European Center for Disease Prevention and Control (ECDC), the diagnostic sensitivity of antigen tests should be ≥ 90%, and specificity ≥ 97% (ECDC 2020). Based on the meta-analysis, their sensitivity and specificity were estimated at 73.8% and 99.7%, respectively (Brümmer et al. 2021). The assay’s practical application is limited to the initial stage of the disease when the number of virions in specimens is at the highest level (up to 5–7 days from the onset of clinical symptoms). These assays have several advantages as lower costs when compared to nucleic acid amplification tests (NAAT), short time of analysis (15–30 minutes); the possible use outside health care facilities, without specialized medical equipment (in workplaces, places of service or in patients’ homes). They are also used in ambulances, emergency rooms, and airports (Flisiak et al. 2021). Due to a small number of SARS-CoV-2 virions in the specimens, the results of antigen tests may be negative. Therefore, positive antigen test results correlate better with the period of COVID-19 infectivity than NAAT’ results (Hirotsu et al. 2020; Yamayoshi et al. 2020).

In the early stages of COVID-19 infection, chest X-rays may not be sensitive enough to detect suggestive changes (Zu et al. 2020). Therefore, computer tomography (CT) is the primary imaging measure (Alsharif and Qurashi 2021). It helps to establish the diagnosis, assesses the extent of inflammatory changes (CT Score), detect complications, and determine the outcome of the disease. Milky glass-like lesions (GGO, ground-glass opacity) predominate in the early phase of infection (up to seven days from the onset of clinical symptoms). After day 7, inflammatory consolidations, cobblestone patterns, and fibrosis predominate. CT score strongly correlates with the severity of the disease, the most strongly with CRP and D-dimer levels. CT score > 18 points is associated with significantly higher mortality (Francone et al. 2020).

The principles of treatment and diagnosis of COVID-19 in Poland have been compiled as recommendations by the Polish Society of Epidemiologists and Infectious Disease Doctors. The disease was divided into four stages: 1) asymptomatic or oligosymptomatic, 2) fully symptomatic, 3) respiratory failure (called cytokine storm stage), and 4) full respiratory distress stage (ARDS).

The oligosymptomatic stage is characterized by normal saturation (SpO₂ > 95%). It does not routinely require hospital treatment, and the patient should self-monitor the saturation at home. In addition, the patients should rest, and take adequate oral hydration, anti-inflammatory drugs, and antipyretics if administered. There are no indications for the routine treatment with antibiotics, anti-influenza drugs, vitamin D₃, low-molecular-weight heparins, or systemically administered corticosteroids (Flisiak et al. 2021).

In October 2021, the FDA (Food and Drug Administration) conditionally approved REGEN-COV (casirivimab, imdevimab), a drug for the treatment and post-exposure prophylaxis of COVID-19. The drug may be used in the oligosymptomatic stage. It should be used to manage COVID-19 in patients at high risk of developing severe disease. The treatment reduces the risk of death and hospitalization and the duration of clinical symptoms (Weinreich et al. 2021). In addition, it may be used in post-exposure prophylaxis after home contact with a sick person to prevent the infection and the symptomatic form of the disease (O’Brien et al. 2021). The exposure is considered to be, among others, a 15-minutes contact within 24 hours, without personal protection equipment (PPE), and within 1,8 meters (CDC 2021). REGEN-COV displays activity against SARS-CoV-2 VOC (Variants of concern) (Weinreich et al. 2021).

Molnupiravir (Lagevrio) is approved to treat mild to moderate COVID-19. It should be administered up to five days after the onset of clinical symptoms. It is used orally, at a dose of 800 mg twice a day for 5 days. It can be administered to patients with impaired liver and kidney functions but is absolutely contraindicated during pregnancy. Although it is a well-tolerated drug, some side effects may include diarrhea, nausea, dizziness, and headache. Molnupiravir is a prodrug that is metabolized to the ribonucleoside analog of N-hydroxycytidine (NHC). In cells, NHC is phosphorylated to triphosphate (NHC-TP). With the participation of RNA polymerase, it binds to viral RNA, leading to the accumulation of errors and suppression of replication (Lagevrio, Summary of Product Characteristics). The drug is considered safe and well-tolerated. It reduces the risk of death and hospitalization due to COVID-19, and lowers the time to SARS-CoV-2 RNA clearance (Singh et al. 2021). Molnupiravir is active in vitro against SARS-CoV-2 VOC (Variants of concern) (Vangeel et al. 2022).

In the full symptomatic stage, the saturation falls below 95%, and the patient usually requires hospitalization. In this phase, causal treatment should be implemented. The first-choice drug is remdesivir (Veklury) administered intravenously (by slow infusion) at 200 mg on the first day and 100 mg on the following four days. During the drug administration, aminotransferase activity should be monitored daily as remdesivir is hepatotoxic. Another possibility of causal treatment is the transfusion of a group-compatible plasma from the recovered patients. Low molecular weight heparins in prophylactic or therapeutic doses may also be implemented. From day two to five of remdesivir therapy, systemic dexamethasone at 4 mg is added and maintained until the end of the second week of the illness. Additionally, passive oxygen therapy could be administered.

Remdesivir is an RNA polymerase inhibitor registered for the treatment of COVID-19 in adolescents (aged 12 years and over and weighing more than 40 kg) and adults with pneumonia when the clinical condition requires oxygen administration. Treatment with remdesivir is contraindicated in renal failure with eGFR below 30 ml/min and liver damage with ALT activity exceeding five times the upper limit of the normal level. It should not be taken during pregnancy, breastfeeding, and hypersensitivity to the active substance and excipients. According to the SmPC (summary of product characteristics), the total treatment duration is five to 10 days. However, a randomized phase 3 clinical trial (Goldman et al. 2020) did not demonstrate an advantage of a 10-day treatment over a 5-day treatment (clinical improvement was obtained in 54 and 64% of subjects, respectively). According to the SmPC, the recommended dose is 200 mg on the first day and 100 mg on subsequent days. The activity of eGFR and ALT should be monitored throughout the treatment period. The most common adverse reactions include liver damage with transaminase elevations, nausea, vomiting, headache, and rash (Veklury, Summary of product characteristics). Remdesivir has in vitro activity against SARS-CoV-2 VOC (Variants of concern) (Vangeel et al. 2022).

Stage 3 disease (respiratory failure, cytokine storm) usually occurs at week 2, SpO₂ falls below 90%, and the patient requires hospitalization. Treatment with tocilizumab can be administered. It is a monoclonal antibody that binds specifically to IL6 receptors (both soluble and membrane-bound). The indication for the tocilizumab administration is IL6 concentration above 100 pg/ml. The drug is administered intravenously (8 mg/ kg, max dose 800 mg) in a one-hour infusion. If there is no improvement after the first dose, a second dose may be given after 8–24 hours. Due to the risk of developing severe bacterial infections, the drug should not be given to patients with a neutrophil count below 2 K/ ml. Liver damage with an ALT value above five times the upper limit of normal level is also a contraindication for therapy (RoActemra, Summary of product characteristics). Simultaneously with tocilizumab, we administer dexamethasone intravenously. The initial dose was 8 mg, and it later was reduced to 4 mg after 5 days of therapy. A group of patients with severe COVID-19-associated pneumonia treated with tocilizumab showed decreased mortality, risk of needing mechanical ventilation, and no increased risk of severe bacterial infections (Lan et al. 2020; Tleyjeh et al. 2021). Also, a study conducted on the Polish population showed a beneficial effect of tocilizumab. A decrease in the systemic inflammatory response (CRP, procalcitonin, fibrinogen concentrations) and a rapid improvement of clinical conditions for most patients were observed (Tomasiewicz et al. 2020).

In the last stage of the disease (ARDS), the patient should be admitted to the Intensive Care Unit, where ventilator therapy or ECMO can be administered.

Lopinavir/ritonavir, chloroquine, and hydroxychloroquine should not be used for the treatment due to a lack of efficacy evidence (Horby et al. 2020; Kashour et al. 2021). In addition, concomitant administration of remdesivir and chloroquine/hydroxychloroquine is not recommended since they exert antagonistic effects on intracellular metabolism and antiviral activity.

The first cases of COVID-19 were diagnosed in Europe in January 2020 (Lescure et al. 2020), and already in December of that year, the EMA (European Medicines Agency) conditionally approved the vaccine from BioNTech and Pfizer. The trials were based on mRNA platforms (BNT162b2, Pfizer, BioNTech; mRNA-1273, Moderna), vector vaccines (chAdOx1-S, AstraZeneca; Ad26.COV2-S, Janssen) or inactivated viruses (e.g., Sinovac). Laboratory evaluation of vaccine efficacy was performed by assessing their ability to induce the production of total IgG or neutralizing antibodies (against S protein) and their ability to induce a response from T lymphocytes (assessed, for example, by the ELISpot assay) (Ong et al. 2021). The highly immunogenic S protein plays a crucial role in infecting the host cells – it mediates the entry of virions into the host cells. Blocking it will prevent the virus from multiplying; hence, this protein has become the primary target of COVID-19 vaccines. Sequencing of SARS-CoV-2 new variants provides valuable information about mutations, especially in S protein, which may lead to the development new and better vaccines that provide long-term protection and neutralize mutated variants of SARS-CoV-2 (Harvey et al. 2021; Malik et al. 2022).

The most effective vaccines are those based on the mRNA platform (BNT162b2; mRNA-1273), with an efficacy of more than 94% (McDonald et al. 2021). A Russian-made vector vaccine (rAd26/Ad5) is also characterized by high (91.6%) efficacy and good tolerability (Logunov et al. 2021).

Vaccines (BNT162b2, mRNA-1273, ChAdOx1) are effective against VOC ((Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), and Delta (B.1.617.2)). After two doses of vaccines, their effectiveness against symptomatic infection on the seventh day and later was 89–92%, 87%, 88%, 82–89%, and 87–95% against Alpha, Beta, Gamma, Beta/ Gamma, and Delta variants, respectively (Nasreen et al. 2022). T cells after vaccination recognize SARS-CoV-2 VOC, including Omicron (B.1.1.529). Receptor-binding domain (RBD) of the memory B cells’ recognition of Omicron is reduced to 42% compared to other variants (Tarke et al. 2022).

By the end of January 2022, almost 10 billion doses of COVID-19 vaccines had been administered, and over 4 billion people were fully vaccinated. In Poland, 21,81 million people, i.e., 57,46% of the population, were fully vaccinated, which is the worst result compared to Western European countries. In France, 76.94% of the population was fully vaccinated; in Germany – 74.29%, Italy – 76.79%, the United Kingdom – 72.56%, Ireland – 78.75%, and Spain – 81.38% (The Johns Hopkins Coronavirus Resource Center 2022).

Comirnata’s mRNA vaccine (BNT162b2) was the first to receive the EMA approval and was authorized for use in Europe. It was registered for the active immunization of adults and children above 12 years of age. Basic vaccination is administered intramuscularly in a two-dose schedule. A booster dose (third dose) is given six months after the second. The most common adverse reactions were injection site pain, fatigue, headache, myalgia, chills, fever, and swelling at the injection site (Comirnaty, Summary of product characteristics).

Moreover, there is an increased risk of myocarditis and pericarditis (more often after the second dose in young males) (Chua et al. 2021). mRNA is encapsulated in lipid nanoparticles, allowing its delivery to host cells. Then, the SARS-CoV-2 antigen (glycoprotein S) is transiently expressed. It is characterized by high efficacy, 95% for the primary endpoint of the study, and 91.3% for the secondary endpoint (Comirnaty, Summary of product characteristics). A dose of ten micrograms is approved for children 5 to 11 years of age.

Another mRNA vaccine is Spikevax, previously COVID-19 Vaccine Moderna (mRNA-1273-SARS-CoV). The vaccine contains mRNA encoding the SARS-CoV-2 S protein encapsulated in lipid nanoparticles. After intramuscular injection, cells at the administration site absorb the lipid nanoparticles, providing mRNA sequences for translation and protein biosynthesis. The delivered mRNA does not penetrate the cell nucleus, does not interact with the genome, is incapable of replication, and its expression is transient. The immune system recognizes the S protein as a foreign antigen, which triggers a response from T and B lymphocytes. The vaccine is licensed for administration to persons 12 years of age and older. The basic vaccination cycle consists of two intramuscular doses (0.5 ml each), with the second dose given 28 days after the first dose. A booster vaccination is given six months after the second dose (0.25 ml, 50 micrograms mRNA) (Spikevax, Summary of product characteristics). Vaccine efficacy has been estimated at 94% and is characterized by the strong humoral response directed against S-glycoprotein (McDonald et al. 2021). The most reported adverse reactions were mild, self-limiting, and included injection site pain, fatigue, headache, muscle pain, joint pain, chills, nausea or vomiting, armpit swelling or tenderness, fever, injection site swelling, and redness. These occurred more frequently in younger age groups (Spikevax, Summary of product characteristics). There is an increased risk of myocarditis and pericarditis following vaccination. It most often develops a few days after vaccination, mainly within 14 days. These site effects are more common in young males. Clinical course is similar as in unvaccinated person (Gargano et al. 2021).

Another vaccine is ChAdOx1-S (COVID-19 Vaccine AstraZeneca, Vaxzevria). It is a replication-deficient recombinant chimpanzee adenovirus with the S (Spike) SARS-CoV-2 glycoprotein coding sequence and is registered for persons over 18 years of age. The vaccination cycle consists of two doses, with the second given between the fourth to 12 weeks after the first. Common adverse reactions are similar to other vaccines, more pronounced after the first dose. Local synthesis of the SARS-CoV-2 glycoprotein S, production of neutralizing antibodies, and stimulation of the cellular immune response occur following the vaccine administration. ChAdOx1-S induces the most potent T cell response in the ELISpot assay (McDonald et al. 2021).

Based on registration studies (COV002, COV003), vaccine efficacy in subjects vaccinated with two doses was estimated at 59.5%. In the group vaccinated with two doses of those who contracted COVID-19, none required hospitalization (COVID-19 Vaccine Astra-Zeneca, Summary of product characteristics).

COVID-19 Vaccine Janssen (Ad26.COV2-S) is a monovalent recombinant vaccine. It is based on adenovirus type 26 (with no replication capability) containing the S-glycoprotein coding sequence. It was registered for use in adults. It has the advantage of being a single-dose scheduled vaccine. Common adverse reactions were similar to other preparations. Sporadic cases of vaccine-associated immune thrombosis and thrombocytopenia (VITT) (Franchini et al. 2021) and capillary leaking syndrome (CLS) have been observed (Choi et al. 2021). Based on the registration studies, the vaccine efficacy was estimated at 66.1–85.4%. Moreover, vaccine efficacy was higher against severe/ critical COVID-19 (Sadoff et al. 2021; COVID-19 Vaccine Janssen, Summary of product characteristics).

A very rare (1:150,000 vaccinations) and severe complication following the use of the vaccines based on adenoviral vectors (ChAdOx1-S; less commonly Ad26. COV2-S) is the VITT syndrome. It is characterized by a venous thrombosis of atypical localization, usually in the central nervous system (CNS), or as a sinus vein thrombosis (CVT). Coexisting thrombocytopenia indicates the immunological background of the pathology. Severe thrombocytopenia (defined as a platelet level < 25,000/μl) was present in 52.6% of cases. The syndrome most commonly affects young women, who take oral hormonal contraception often, with a median age of 40.5 years. Mortality is high, above 40%; furthermore, mortality increases to 60% when CVT occurs. In treatment, novel oral anticoagulants (NOAC’s), intravenous immunoglobulin preparations (IVIG, 1 g/kg body weight), and high-dose steroid therapy (e.g., dexa-methasone 40 mg/ day for 4 days) are administered to the patients (Franchini et al. 2021). Brief characteristics of all vaccines are presented in Table I.

The Omicron variant (B.1.1.529) was first reported in South Africa and recognized as VOC by the WHO on November 26, 2021 (WHO 2021). The Omicron SARS-CoV-2 has over 30 mutations in the Spike protein sequence, 15 of which are located within the RBD (Receptor Binding Domain). Omicron has a higher affinity for ACE2 than the Delta variant due to many mutations in RBD (Kumar et al. 2022). Omicron may be associated with a higher risk of reinfection with SARS-CoV-2 (Pulliam et al. 2021). Diagnostics of this variant do not differ from other variants. Omicron does not substantially impair NAATs performance except for S gene target failure (SGTF). It is a deletion at the amino acid positions 69 and 70 of the Spike protein sequence (it was also noticed in the Alpha variant). So, detecting only one genomic region in PCR may lead to a false-negative result, but NAATs commonly detect two, three, or four genomic regions (Ferré et al. 2022; Kumar et al. 2022). Treatment of Omicron-induced COVID-19 is the same as the disease caused by other variants (Araf et al. 2022).

The coronavirus genome (SARS-CoV) size is between 26 and 32 kilobases and has a variable number of open reading frames (ORFs) (Song et al. 2019; Kim et al. 2020). SARS-CoV-2 is an enveloped virus with a single-stranded RNA genome with positive polarity (Song et al. 2019). The genome at the 3ʹ end contained ORFs encoding structural and accessory proteins: surface glycoprotein (S), a small envelope protein (E), matrix protein (M), nucleocapsid (N), and eight accessory proteins: 3a, 3b, p6, 7a, 7b, 8b, 9b, and ORF 14 (Wu et al. 2020). The 20 kb region at the 5ʹ end contains two OFRs (ORF1a/1ab) that encode 16 non-structural proteins (NSP11-NSP16) essential for virus replication (Davies et al. 2020). In addition, the genome contains non-coding sequences and sequences that regulate the transcription process (i.e., hexanucleotide sequence ACGAAC and CUAAAC) (Yoshimoto 2020).

Sequence analysis of 103 SARS-CoV-2 genomes allowed the identification of two distinct lineages (L and S), characterized by different polymorphic variants at positions 8782 (ORF1ab: T8517C) and 28144 (ORF8: C251T, S84L). The L line was more frequent than the S line in the trial studies. The T8517C mutation in ORF1ab does not affect the sequence of the encoded protein (changes the codon AGT (Ser) to AGC (Ser). However, it may be necessary in the translation of ORF1ab because the AGT codon is preferentially chosen in contrast to AGC (Tang et al. 2020). ORF8 promotes the expression of ATF6, which, being an activating transcription factor, plays an essential role in the transcription of endoplasmic reticulum chaperone proteins (Hu et al. 2017).

The genome sequence of SARS-CoV-2 from infected patients and the GenBank database showed discrepancies indicating the significant variability of the pathogen. For coronaviruses, the average rate of evolutionary change is approximately ~ 3 × 10^–4 substitutions during each replication cycle (Pyrc et al. 2006). These mutations may affect the transmission and infectivity of SARS-CoV-2. Therefore, further studies considering the unique genomic features, epidemiological information, and medical records of the clinical course of COVID-19 are essential.

After entry into host cells, viral genomic RNA is translated. Viral non-structural proteins (NSPs) provide conditions favorable for viral infection and viral mRNA synthesis (Lim et al. 2016). Among those, non-structural protein 1 (NSP1), known as a potent inhibitor of host gene expression, is crucial. NSP1 binds to ribosomes reducing the pool of those that can participate in the translation process. More efficient translation based on viral mRNA will be favored over host mRNA in such conditions. The most recent studies show that NSP1 inhibits the translation of both native and reporter mRNAs containing the 5ʹUTR region of the virus. SARS-CoV-2 can adjust the cellular concentration of NSP1 to values that will prevent the reduction of viral mRNA translation but will be sufficient to inhibit translation initiation from less efficiently recruited cellular mRNAs. This mechanism helps to explain how infected host cells target their protein synthesis machinery to produce foreign pathogen proteins.

The above mechanism ensures efficient viral gene expression (Schubert et al. 2020). NSP2 is a highly conserved protein specific to SARS-CoV viruses. It binds to two host proteins, prohibitin 1 and prohibitin 2 (PHB1, prohibitin 1; PHB2, prohibitin 2). These proteins play essential roles in regulating cell cycle progression, cell migration, cell differentiation, apoptosis, and mitochondrial biogenesis. The binding of NSP2 to PHB1 and PHB2 proteins suggests that NSP2 disrupts the host cell environment and interferes with its intracellular signaling (Cornillez-Ty et al. 2009; Yoshimoto 2020). The 200 kDa multi-domain NSP3 protein is the largest protein encoded by SARS-CoVs (not including ORF 1a and ORF1ab). It has a variable structure. Individual coronaviruses can have between 10 and 16 NSP3 domains, of which eight and two transmembrane regions of the protein are highly conserved. NSP3 is an essential component of the replication and transcription complex. This protein interacts with the NSP4 in the host cell membrane rearrangement, which is necessary to initiate the replication and transcription of the viral genome (Sakai et al. 2017; Lei et al. 2018).

The nucleocapsid (N) protein is a multifunctional structural protein. Its domain structure consists of three distinct, highly conserved components: an N-terminal RNA-binding domain (NTD domain), a C-terminal dimerization domain (CTD domain), and a region rich in serine/arginine (SRD domain) (Zeng et al. 2020). The N protein allows the viral genomic RNA to be packaged into long, helical ribonucleoprotein complexes called nucleocapsids. The nucleocapsid protects the genome and ensures its replication and transmission of the pathogen. The N protein interacts with viral membrane proteins during viral assembly and assists in RNA synthesis and folding (McBride et al. 2014). Previous studies also indicate its role in regulating host-pathogen interactions by affecting the cell cycle progression and host cell apoptosis.

Furthermore, the N protein is highly immunogenic and induces host immune responses (Surjit et al. 2006; Du et al. 2008). The analyses confirmed the presence of IgA, IgM, and IgG class antibodies against the N antigen in the sera of patients with COVID-19, indicating this protein’s role in modulating immunity and its important diagnostic value (Zeng, 2020). S protein is another protein of structural importance. Each protein monomer, with a molecular mass of approximately 180 kDa, consists of two functional subunits responsible for binding to the host cell receptor (S1 subunit) and viral and cell membrane fusion (S2 subunit) (Ou et al. 2020). The S1 subunit is characterized by the presence of an N-terminal and ACE2 receptor-binding domain (RBD). The S2 subunit is composed of a fusion peptide (FP) domain, heptapeptide 1 (HR) repeat sequence, heptapeptide 2 (HR2) repeat sequence, transmembrane domain (TM), and cytoplasmic domain (CT) (Huang et al. 2020). The S protein is involved in two key events: binding to the cell receptor and inducing the fusion between the viral and cell membranes (Hussain et al. 2020).

Completion of these two processes leads to the entry of the viral RNA genome into the host cell and the subsequent start of the viral replication cycle (Belouzard et al. 2012). Replication and virions release lead to cell pyroptosis, a highly inflammatory form of cell death commonly found in many viruses. Additionally, pattern recognition receptors (PRRs), including Toll-like receptors (TLRs) detect molecular patterns of the pathogen, which initiates an increased inflammatory response and triggers the release of pro-inflammatory cytokines from neighboring cells (McFadyen et al. 2020). It is also worth noting that SARS-CoV-2 significant transmission and infectivity may be due to the unique, compared to other CoVs, a furin cleavage site in protein S. The significant difference lies in the insertion of four amino acids (PRRA) that form a cleavage site recognized by the ubiquitous protease furin. Probably due to this additional amino acid insertion, SARS-CoV-2 gains the ability to infect organs or tissues that are not susceptible to other coronaviruses, thus providing ample opportunities for the virus to be released into the environment and significantly increasing its transmission (Wang et al. 2020). The M protein is a type III glycoprotein, the smallest among other structural proteins. According to the studies performed on SARS-CoV and MERS-CoV viruses, the M protein consists of 230 amino acids, and its molecular mass is between 25–35 kDa (Asghari et al. 2020). It has a short ectodomain at the amino end, three transmembrane domains, and a long domain at the carboxyl end of the protein.

The M protein gives the virus shape and is crucial in coordinating virus folding and forming mature viral envelopes. Homotypic interactions between M proteins are the main factor for viral envelope formation but are not sufficient by themselves for efficient viral assembly (de Haan et al. 2000). Presumably, only the interactions of the NSP3, M, and N proteins stabilize the nucleocapsid and help bind the viral genome to the replicase-transcriptase complex (RTC), which promotes the completion of viral particle assembly (Fehr and Perlman 2015).

Among all structural proteins, protein E is less understood. The E protein is hydrophobic and consists of a sequence of 74–109 amino acids with a molecular weight of 8.4–10.9 kDa. This protein plays an important role in the pathogenesis of coronaviruses. SARS-CoV with abolished E protein channel activity was significantly less infectious (Nieto-Torres et al. 2014; Asghari et al. 2020). The analyses confirmed the cation-conducting ability of the SARS-CoV-2 E-protein, which is considered a potential ion channel and an important target for the anti-COVID-19 drug development strategy (Singh Tomar and Arkin 2020). The comparative analysis of the SARS-CoV-2 genome with other coronaviruses revealed significant differences in the sequences of the multi-domain NSP3 protein and the surface protein S. The NSP3 protein contains a new large sequence insertion between two structurally and probably functionally independent domains. This newly added extended peptide sequence may act as a long inter-domain junction, extending the conformational flexibility of the protein. Furthermore, the presence of four new insertions and one highly variable surface region of the S protein has been demonstrated, which may have a structural impact on the homo-trimeric form of this protein and the functions performed by the N-terminal domain of the NTD (Srinivasan et al. 2020).

The receptor recognition is the first step of viral infection and a key determinant of host cell and tissue tropism (Walls et al. 2020). It is also an important target in a host immune surveillance and the intervention strategies (Shang et al. 2020). The SARS-CoV-2 virus uses human angiotensin-converting enzyme 2 (ACE2) as a receptor (Zou et al. 2020). ACE2 is a type I surface glycoprotein of approximately 100 kDa, consisting of 805 amino acid residues. In its structure, one can identify two functional domains: the N-terminal M2 peptidase domain and the C-terminal cofilin domain (Hussain et al. 2020). The receptor is encoded by the ACE2 gene, which locus is on chromosome X (Xp22), spreads over 39.98 kb of genomic DNA, and contains 18 exons and 20 introns. The mechanism of S protein binding to ACE2 may depend on several factors. Upon binding, the RBD receptor-binding domain of the S1 subunit undergoes the conformational change. The dynamic structure of RBD can occur in two forms: “standing”, which binds the receptor, and “lying”, which responds to a state of unavailability to the host cell receptor and allows it to evade immune surveillance (Yuan et al. 2017).

Furthermore, the S protein of SARS-CoV-2 acquired mutations that enhance its affinity for the human receptor by 10–20 fold compared to the SARS-CoV, making it highly infectious (Luan et al. 2020). The analysis of glycoprotein S structure using cryo-electron microscopy (cryo-EM) revealed that the RBD domain of SARS-CoV-2 is mainly in a “lying” position, making it unavailable to ACE2 (Wrapp et al. 2020). Therefore, it is hypothesized that despite its high affinity, the RBD domain may similarly, or to a lesser degree, interact with the ACE2 receptor compared to SARS-CoV (Shang et al. 2020).

After engaging the human receptor (ACE2), S protein is modified by the membrane-bound trans-mem-brane serine protease type II (TMPRSS2). The S-glycoprotein is cleaved into S1 and S2 subunits, allowing membrane fusion, the release of the virus into host cells, and its further spread (Hoffmann et al. 2020). Recent analyses of intestinal epithelial cell lines (HEK293) confirmed the importance of TMPRSS2 coexpression with ACE2 in enhancing virus infectivity. However, not only the serine protease TMPRSS2 plays an important role in S-protein cleavage and enhancement of the membrane fusion. The expression of TMPRSS4 also significantly increases viral RNA in the presence of ACE2. In addition, coexpression of TMPRSS4 and TMPRSS2 shows an additive effect favoring maximum infectivity in the cell culture (Zang et al. 2020).

On the other hand, the role of host cell proteases in SARS-CoV infection may not be limited to protein S cleaving. Previous studies suggested that ACE2 might be proteolytically modified by host cell proteases, which influenced the entry and pathogenesis of SARS-CoV. Furthermore, the arginine and lysine residues of ACE2 at positions 697 and 716 are essential for receptor modification by these proteases. In summary, TMPRSS2 and potentially related proteases promote SARS-CoV entry into the host cell through two different mechanisms: cleavage of ACE2, which may promote viral uptake, and cutting of S protein, activating it for membrane fusion process (Heurich et al. 2014).

The tissue expression and distribution of the ACE2 receptor determine the tropism of viral infection, which is crucial for understanding pathogenesis and designing therapeutic strategies. Gene expression profile analysis revealed that ACE2 is mainly expressed in alveolar epithelial type II (AECII) cells, suggesting that these cells may serve as a reservoir for SARS-CoV-2. In addition, many other genes that positively regulate viral entry, multiplication, and transmission were found to be expressed in the ACE2-expressing cells (Zhao et al. 2020). ACE2 is expressed in the lungs and in the blood vessels, heart, kidney, and intestines. Therefore, it is now hypothesized that viral transmission via the endothelial ACE2 receptor represents a mechanism by which the virus can efficiently spread in the body, infecting different host tissues (McFadyen et al. 2020). Viral RNA has been found in urine and stool samples from COVID-19 patients, which enlarges the possibility of pathogen spread (Chen et al. 2020; Peng et al. 2020).

Understanding the genetic background that promotes SARS-CoV-2 transmission and infectivity is crucial for identifying the phenotype of patients predisposed to COVID-19 infection and its severity. Despite the increased efforts worldwide and the growing knowledge on the genetic determinants of COVID-19, many questions remain unanswered.

Recent autopsy analysis of the microvascular architecture of the lungs in patients who died from respiratory failure due to COVID-19 and AH1N1 influenza provided some interesting observations. More endothelial ACE2 receptors and significant changes in endothelial morphology were observed in the group with COVID-19 compared to the influenza-infected group. In addition, patients in this group were characterized by more enhanced and sprouting angiogenesis compared to post-influenza patients. The above study was supported by multiplex expression analysis of 323 angiogenesis-related genes. Among these, the expression of 69 genes was different in the COVID-19 group, while the other 26 were differentially expressed in the post-influenza group. A common expression pattern was found for 45 genes in both groups of patients (Ackermann et al. 2020). Recently, an analysis of 1,700 polymorphic variants of ACE2 was also performed. Polymorphisms of this gene may determine changes in its expression and modulate plasma levels of the encoded protein, which increase is considered a molecular marker of bad prognosis of COVID-19 (Zipeto et al. 2020). The obtained results indicated a significant population homogeneity, which did not allow the clear selection of polymorphisms that could limit the binding of protein S to the ACE2 receptor and thus protect against viral transmission. However, it should be noted that 15 polymorphisms (rs112171234, rs12010448, rs143695310, rs1996225, rs200781818, rs2158082, rs4060, rs4646127, rs4830974, rs4830983, rs5936011, rs5936029 rs6629110, rs6632704, and rs75979613) can upregulate ACE2 expression in various host tissues and they are more frequent in the East Asian than the European population. The differences observed between the populations indicate that there may be significant differences in susceptibility and/ or response to SARS-CoV-2 infection, despite similar conditions (Cao et al. 2020). Also noteworthy is the higher mortality rate among men infected with SARS-CoV-2 compared to women. It was initially hypothesized that this might be related to the fact that ACE2 shows an unusual heterogeneous male-female expression pattern, with higher expression in men (Tukiainen et al. 2017). However, closer examination of the expression of the mentioned gene within the lung tissue did not confirm significant differences between patients of different sexes (Cai 2020). Looking for genetic prognostic markers for early identification of COVID-19 high-risk individuals, the rs2285666 (G8790A) polymorphism of ACE2, known as a potential risk factor for hypertension, type 2 diabetes, and coronary artery disease, may be important. Recent studies in the Italian population have shown that replacing the G allele by the A allele increases the strength of the splicing site, reflected by an increase in serum ACE2 protein levels in patients (Asselta et al. 2020). The analyses of the probable correlation between ACE2 expression and the immune response of a patient infected with SARS-CoV-2 also seem interesting. It appeared that lung tissue characterized by a high level of ACE2 expression in women or young people (aged ≤ 49 years) determines a weaker immune response to infection. The high expression of ACE2 in the lung observed in men or older individuals (aged > 49 years) determines stronger immune signatures. It suggests that men and elderly individuals infected with SARS-CoV-2 may be particularly predisposed to an enhanced immune response manifested by a cytokine storm or immunopathological damage (Li et al. 2020a).

The ACE2 locus is not the only one considered therapeutic target. Genome-wide association studies (GWAS) have identified a likely genetic susceptibility locus for developing chronic heart failure (CHF) associated with COVID-19. More than eight million single nucleotide polymorphisms have been analyzed, and statistically significant correlations were shown for loci 3p21.31 and 9q34.2. Within this locus, the frequency of risk alleles of selected polymorphic variants was higher in mechanically ventilated patients than those receiving oxygen supplementation alone. Furthermore, the 3p21.31 locus is characterized by a cluster of genes (SLC6A20, LZTFL1, CCR9, FYCO1, CXCR6, and XCR1) that may shape susceptibility to COVID-19 and be also relevant to the AB0 group system. Among these, of note is the SLC6A20 gene, which encodes sodium-imino acid transporter 1 (SIT1) and interacts with ACE2. In addition, the 3p21.31 locus also contains genes encoding chemokine receptors, including C-C motif chemokine receptor 9 (CCR9) and C-X-C motif chemokine receptor 6 (CXCR6).

The last one regulates the specific localization of tissue-resident memory T cells (TRMs) in the lung, which is the first line of defense against respiratory system pathogens (Wein et al. 2019). This meta-analysis also indicates the association of polymorphic variants: rs11385942 locus 3p21.31 and rs657152 locus 9q34.2 with SARS-CoV-2-induced respiratory failure in Italian and Spanish populations. In addition, higher susceptibility to SARS-CoV-2 infection was observed among individuals with blood group A and a protective effect of group 0 against individuals with other blood groups (Severe COVID-19 GWAS Group et al. 2020; Zhao et al. 2021).

On the other hand, individual susceptibility to SARS-CoV-2 infection and the course of COVID-19 depends on the efficiency of the natural defense mechanisms of the immune system. The ZAP protein, known as ZC3HAV1 in mammals and hZAP in humans, is crucial in the antiviral immune response. The zinc-finger of the antiviral ZAP protein binds to CpG dinucleotides in viral RNA genomes (Meagher et al. 2019). Thus, it inhibits virus replication and simultaneously mediates the degradation of the virus genome (Gao et al. 2002; Mao et al. 2013). The experimental evidence is that CpG deficiency in RNA viruses has evolved in response to specific antiviral activities of the host organism. SARS-CoV-2 is characterized by extreme CpG deficiency, allowing the virus to evolve in host tissue with high expression rates (Meagher et al. 2019). Considering the human body’s defense mechanisms, the endogenous enzyme APOBEC3G, found in innate immune cells, also exhibits antiviral activity. Both ZAP and APOBEC3G show tissue-specific expression patterns in humans. The enzyme modifies viral genetic material by deaminating cytosine to uracil (Sharma et al. 2015). CpG modification on UpG in non-functional regions may reduce the virus’s susceptibility to the ZAP protein activity. In contrast, deamination induced by APOBEC3G may contribute to CpG deficiency and thereby reduce the antiviral effect of ZAP protein directed against CpG sequences (Xia 2020).

Another gene involved in the regulation of the immune system is IFITM3, which encodes interferon-induced transmembrane protein 3. Such trans-mem-brane proteins are involved in the innate immune response to viral infections. It regulates the fusion of invading endocytic vesicles and directs them to lysosomes. Furthermore, IFITM3 can alter cell membrane stiffness and curvature to inhibit viral membrane fusion and further pathogen spread (Li et al. 2013; Suddala et al. 2019). Previous reports indicate that the C allele of the rs12252 IFITM3 polymorphism is associated with an increased probability of influenza virus infection as well as an increased risk of severe influenza and death (Xuan et al. 2015). Therefore, it should be investigated within the Chinese population, where the prevalence of the CC genotype is 26.5%. Analyses in Chinese patients have shown an association of the mentioned genotypic variant with an increased predisposition to severe SARS-CoV-2 infection (Thevarajan et al. 2020; Zhang et al. 2020). Presumably, the severe and complicated course of COVID-19 may be conditioned by selected polymorphic variants of IFITM3 affecting gene expression or the lack of functionality of the encoded IFITM3 protein, which favors SARS-CoV-2 tropism (Hachim et al. 2020).

Another hypothesis for human susceptibility to COVID-19 is based on the human leukocyte antigen (HLA) system, responsible for presenting viral antigens to T lymphocytes. Different HLA alleles may determine individual susceptibility to COVID-19 infection and spread, as confirmed for SARS and MERS virus (Li et al. 2020b). Analysis of HLA variation and effects on cellular immune response in patients with confirmed COVID- 19 revealed that HLA-B*46:01 contains the lowest number of SARS-CoV-2 binding sites, suggesting that individuals with this allele may be particularly predisposed to SARS-CoV-2 infection. In contrast, HLA-B*15:03 shows the most remarkable ability to present highly conserved SARS-CoV-2 peptides, which may enable T-cell-based defensive immunity (Nguyen et al. 2020). Studies performed within the Italian population have revealed that healthy individuals carrying the HLA-B*44 and/ or C*01 alleles, and the HLA-A*25, HLA-B*08 alleles may be more susceptible to SARS-CoV-2 infection. It may be caused by the limited ability to present viral epitopes and, consequently, the inability to initiate a rapid and effective antiviral immune response. Based on these data, it can be hypothesized that the virus will efficiently replicate and spread from these patients’ oral and pharyngeal mucosal areas. In addition, previous studies indicate the role of the HLA-B*44 alleles in question, as well as C*01, in the progression of inflammatory autoimmune diseases (Grams et al. 2002; Jung et al. 2016), highlighting their involvement in the process of inducing imperfect and often undesirable immune responses (Correale et al. 2020). Noteworthy, the HLA-C*01 allele represents a specific ligand for killer cell immunoglobulin-like receptors (KIRs). These receptors can inhibit the activity of cells that represent the first line of host defense against the onset of a specific T cell response. In the Chinese population, the alleles particularly predisposing to COVID-19 infection may be HLA-C*07:29 and B*15:27, but further analyses are needed due to the small size of the study group (Wang et al. 2020).

The severe and complicated course of COVID- 19 may also be conditioned by activating the transmembrane serine protease TMPRSS2. TMPRSS2 promotes SARS-CoV-2 entry into the host cell via two distinct mechanisms: by cutting SARS-S, which activates the S protein to fuse with the membrane, or cleaving ACE2, which may promote viral uptake via a cathepsin L-dependent pathway (Heurich et al. 2014). Furthermore, among the hypotheses explored regarding gender-dependent differences in susceptibility to infection, the analysis of the expression and genetic variability of TMPRSS2 seems interesting. Data extracted from a large Italian cohort compared the frequency of the selected polymorphisms between Italian, European, and East Asian populations. The study identified four polymorphic variants (rs2298659, rs17854725, rs12329760, and rs3787950) with a significantly different frequency between the mentioned populations. Moreover, a haplotype consisting of polymorphisms: rs2070788, rs9974589, rs7364083 and rs463727, rs34624090, rs55964536, rs734056, rs4290734, rs34783969, rs11702475, rs35899679, and rs35041537, typical of the European population and not found in the Asian population, may regulate androgen-dependent TMPRSS2 expression. It could explain why men are more susceptible to severe COVID-19 and die more often (Asselta et al. 2020). Also noteworthy is the haplotype of the variants: rs2070788, rs9974589, and rs7364083, whose minor allele frequency (MAF) is significantly more frequent in European than in East Asian populations. Previous studies indicate a role for the rs2070788 polymorphism in the severity of influenza A(H7N9) and A(H1N1) (Cheng et al. 2015).

Further data concerned the expression profiles of the TMPRSS2 gene. Polymorphisms rs464397, rs469390, rs2070788, and rs383510 were found to determine TMPRSS2 expression in lung tissue. The allele frequency of each polymorphism was then estimated in African, American, European, and three Asian populations (Chinese, Japanese and Taiwanese).

Interestingly, variants that increase TMPRSS2 expression occur much more frequently in European and American populations than in Asian populations, which may indicate an increased susceptibility of these populations to SARS-CoV-2 infection. Recent genome-wide association studies (GWAS) involving patients from intensive care units in the U.K. identified new loci predisposing to severe and complicated COVID-19. At least two biological mechanisms determine the severe course of the disease. The first one is the innate antiviral defense mechanism, which is particularly important in the early stages of the disease (IFNAR2 and OAS genes) and the life-threatening inflammatory lung damage, which is a consequence of the later pathomechanism (DPP9, TYK2, and CCR2 genes). A new significant genome-wide association rs10735079 (chr12q24.13), located in the oligoadenylate interferon synthetase gene cluster (OAS1, OAS2, and OAS3) and encoding antiviral restriction enzyme activators were identified, and replicated. The mentioned associations were described for the rs2109069 (chr19p13.2) polymorphism in the gene encoding tyrosine kinase 2 (TYK2), the rs2109069 (ch19p13.3) polymorphism in the gene encoding dipeptidyl peptidase 9 (DPP9), and the rs2236757 (chr21q22.1) polymorphism in the IFNAR2 gene, which encodes the receptor for interferons Alpha and Beta variants.

The rapid spread of the SARS-CoV-2 virus favors its molecular evolution. In September 2020, Alpha (B.1.1.7), also known as VUI 202012/01 (Variant Under Investigation, the year 2020, month 12, variant 01), the first new variant of the SARS-CoV-2 virus was discovered. Phylogenetic analyses have shown that SARS-CoV-2 variants accumulate nucleotide mutations at approximately 1–2 mutations per month (Duchene et al. 2020). The VUI 202012/01 has a higher number of genetic changes acquired during global transmission. The genetic structure of the Alpha variant contains eight mutations within the S protein (deletions 69-70HV, 144Y, substitutions N501Y, A570D, P681H, T716I, S982A, D1118H). Attention was paid to three of these, including N501Y within the S protein RBD domain, which the virus uses to bind to the human ACE2 receptor. The mutation involves one of the six key amino acid residues of the domain and can affect the increased affinity of the virus for host cell receptors, increasing its infectivity and virulence (Starr et al. 2020; Santos and Passos 2021). The 69-70del mutation plays a vital role in the viral mechanism of human immune response evasion. The P681H mutation is directly adjacent to the furin S1/S2 cleavage site and may be biologically relevant (Rambaut et al. 2020). Additionally, newly detected mutations affect ORF8. This small protein contains 121 amino acids. A stop codon at position 27 due to the mutation Q27stop, results in the mentioned protein loss of function. In variant B.1.1.7, there were also six synonymous mutations in ORF1ab (C913T, C5986T, C14676T, C15279T, C16176T) and one in the M gene (T26801C). The new VUI 202012/01 strain has many genetic alterations and has rapidly led to a significant increase in COVID-19 cases in the U.K. So far, other variants of SARS-CoV-2 have also been identified, such as Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and the latest Omicron (B 1.1.529). Of concern is the rapid spread of the Omicron SARS-CoV-2 variant, which was first recorded in South Africa. This variant was characterized by approximately 50 different mutations, including more than 30 in the spike protein. Some identified mutations, e.g., 69–70del, T95I, G142D/143–145del, K417N, T478K, N501Y, N655Y, N679K, and P681H, overlap with those in the Alpha, Beta, Gamma, or Delta variants. Omicron may be associated with increased infectivity and the ability to avoid the host’s antibodies blocking the infection, primarily due to mutations in the furin cleavage site (Karim and Karim 2021; Callaway 2021; Torjesen 2021).

Phenotypic consequences of the SARS-CoV-2 mutations and the unknown effects of their co-occurrence indicate that further research and enhanced genomic surveillance worldwide are needed.

Coronaviruses are commonly found in our surroundings, constantly pressured by anthropogenic influences. It is also worth noting that, in recent years, there have been re-peated warning signals indicating the possibility of the emergence of a virus with pandemic potential. However, scientists’ suggestion on the risk of such a crisis has been largely downplayed. Therefore, it is so essential to draw the correct conclusions from the recent experience, especially as COVID-19 is probably not the last threat of viral and etiology that will confront us.