ABSTRACT
Background: High mobility group box 1 (HMGB1) and interleukin-18 (IL-18) are involved in various non-coronavirus disease pathogenesis and are reported as potential biomarkers for coronavirus disease (COVID-19). However, their association with COVID-19 pathogenesis has not yet been explored. Aim: This study aimed to investigate the association between HMGB1 and IL-18 concentrations in the sera of COVID-19 patients versus non-COVID-19 patients. Material and methods: We used stored serum samples obtained from 30 COVID-19 patients and 30 non-COVID-19 patients. We collected data on age, gender, treatment status, principal diagnosis, and comorbidity from patient medical records. HMGB1 and IL-18 concentrations were analyzed in the serum by enzyme-linked immunosorbent assay (ELISA). The swab samples' RT-PCR cycle threshold (CT) values were obtained from the laboratory database. Results: HMGB1 concentrations were increased in the COVID-19 inpatients and non-COVID-19 inpatients compared to non-COVID-19 outpatients (COVID-19 inpatients vs. non-COVID-19 outpatients: 151.33 (90.27-192.38) vs. 80.75 (54.16-128.72) ng/ml; p = 0.0316; non-COVID-19 inpatients vs. non-COVID-19 outpatients: 152.66 (104.04-288.51) vs. 80.75 (54.16-128.72) ng/ml; p = 0.0199). IL-18 concentrations were also higher in the COVID-19 inpatients and non-COVID-19 inpatients compared to non-COVID-19 outpatients (COVID-19 inpatients vs. non-COVID-19 outpatients: 620.00 (461.50-849.6) vs. 403.10 (372.70-556.90) pg/ml; p = 0.0376; non-COVID-19 inpatients vs. non-COVID-19 outpatients: 835.70 (558.30-1602.00) vs. 403.10 (372.70-556.90) pg/ml; p = 0.0026). Moreover, HMGB1 was associated with IL-18 concentrations in the sera of COVID-19 inpatients (p = 0.0337; r = 0.5500). Conclusion: The association of HMGB1 and IL-18 in COVID-19 might indicate the potential for a dangerous cycle leading to a cytokine storm to occur.
ABSTRACT
D614G mutation plays a significant role in the transmissibility of SARS-CoV-2. Identification of other mutations related to D614G mutation within the Spike protein is pivotal as they might contribute to the pathogenicity of SARS-CoV-2. This study aims to analyze the mutation rate of furin cleavage site (FCS) region of Indonesian origin SARS-CoV-2 and to predict the effect of mutation against Spike priming efficiency by furin. A total of 375 sequences of Indonesian isolates obtained during the early pandemic were used for mutation analysis. Mutation analysis includes mutation pattern, variability, frequency of mutation, amino acid conservation, and mutation rate. The effect of mutation against Spike priming efficiency by furin protease from eight sequences with mutation in the FCS region was analyzed by protein-protein docking. We showed that mutations related to the G614 variant were increasing through time, in contrast to the D614 variant. The FCS region at the position 675-692 contained the most variable (66.67%) as well as the highest mutation frequency (85.92%) and has been observed to be the hotspot mutations linked to the D614G mutation. The D614G hotspot-FCS region (residue 600-700) had the highest amino acid change per site (20.8%) as well as the highest mutation rate as 1.34 × 10-2 substitution per site per year (95% CI 1.79 × 10-3-2.74 × 10-2), compared with other Spike protein regions. Mutations in the FCS region were the most common mutation found after the D614G mutation. These mutations were predicted to increase the Spike priming efficiency by furin. Thus, this study elucidates the importance of D614G mutation to other mutations located in the FCS region and their significance to Spike priming efficiency by furin. Supplementary Information: The online version contains supplementary material available at 10.1007/s13337-023-00827-w.
ABSTRACT
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a new virus responsible for the COVID-19 pandemic. The emergence of the new SARS-CoV-2 has been attributed to the possibility of evolutionary dynamics in the furin cleavage site (FCS) region. This study aimed to analyze the sequence of the FCS region in the spike protein of SARS-CoV-2 isolates that circulated in the Special Region of Yogyakarta and Central Java provinces in Indonesia. The RNA solution extracted from nasopharyngeal swab samples of confirmed COVID-19 patients were used and subjected to cDNA synthesis, PCR amplification, sequencing, and analysis of the FCS region. The sequence data from GISAID were also retrieved for further genome analysis. This study included 52 FCS region sequences. Several mutations were identified in the FCS region, i.e., D614G, Q675H, Q677H, S680P, and silent mutation in 235.57 C > T. The most important mutation in the FCS region is D614G. This finding indicated the G614 variant was circulating from May 2020 in those two provinces. Eventually, the G614 variant totally replaced the D614 variant from September 2020. All Indonesian SARS-CoV-2 isolates during this study and those deposited in GISAID showed the formation of five clade clusters from the FCS region, in which the D614 variant is in one specific cluster, and the G614 variant is dispersed into four clusters. The data indicated there is evolutionary advantage of the D614G mutation in the FCS region of the spike protein of SARS-CoV-2 circulating in the Special Region of Yogyakarta and Central Java provinces in Indonesia.