The association between high mobility group box 1 (HMGB1) and Interleukin-18 (IL-18) serum concentrations in COVID-19 inpatients.

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.

Sequence analysis of the Spike, RNA-dependent RNA polymerase, and protease genes reveals a distinct evolutionary pattern of SARS-CoV-2 variants circulating in Yogyakarta and Central Java provinces, Indonesia.

During the Covid-19 pandemic, the resurgence of SARS-CoV-2 was due to the development of novel variants of concern (VOC). Thus, genomic surveillance is essential to monitor continuing evolution of SARS-CoV-2 and to track the emergence of novel variants. In this study, we performed phylogenetic, mutation, and selection pressure analyses of the Spike, nsp12, nsp3, and nsp5 genes of SARS-CoV-2 isolates circulating in Yogyakarta and Central Java provinces, Indonesia from May 2021 to February 2022. Various bioinformatics tools were employed to investigate the evolutionary dynamics of distinct SARS-CoV-2 isolates. During the study period, 213 and 139 isolates of Omicron and Delta variants were identified, respectively. Particularly in the Spike gene, mutations were significantly more abundant in Omicron than in Delta variants. Consistently, in all of four genes studied, the substitution rates of Omicron were higher than that of Delta variants, especially in the Spike and nsp12 genes. In addition, selective pressure analysis revealed several sites that were positively selected in particular genes, implying that these sites were functionally essential for virus evolution. In conclusion, our study demonstrated a distinct evolutionary pattern of SARS-CoV-2 variants circulating in Yogyakarta and Central Java provinces, Indonesia.

Molecular Docking Study of the C-10 Massoia Lactone Compound as an Antimicrobial and Antibiofilm Agent against Candida tropicalis.

Antimicrobial resistance is now considered a global health problem because it reduces the effectiveness of antimicrobial drugs. According to the World Health Organization (WHO), the highest mortality rate is associated with infections caused by multidrug-resistant microorganisms, with approximately 700,000 deaths worldwide each year. The aim of this study was to determine the potential of C-10 massoia lactone to inhibit the growth of fungi and C. tropicalis biofilm, and molecular docking studies were performed to determine the nature of the inhibition. The study was conducted using the microdilution method for antifungal and antibiofilm testing and designed with a molecular docking approach. Furthermore, an analysis using the scanning electron microscope (SEM) was performed to evaluate the mechanism of effect. The results obtained showed that C-10 massoia lactone can inhibit the growth of fungi by 84.21% w/v. Meanwhile, the growth of C. tropicalis biofilm in the intermediate phase was 80.23% w/v and in the mature phase was 74.23% w/v. SEM results showed that C-10 massoia lactone damaged the EPS matrix of C. tropicalis so that hyphal formation was hindered due to damage to fungal cells, resulting in a decrease in attachment, density, and lysis of C. tropicalis fungal cells. Based on molecular docking tests, C-10 massoia lactone was able to inhibit biofilm formation without affecting microbial growth, while docking C-10 massoia lactone showed a significant binding and has the potential as an antifungal agent. In conclusion, the C-10 massoia lactone compound has the potential as an antibiofilm against C. tropicalis, so it can become a new antibiofilm agent.

Prognostic factors for the outcomes of COVID-19 patients infected with SARS-CoV-2 Omicron and Delta variants.

BACKGROUND: The SARS-CoV-2 Omicron variant has replaced the previously dominant Delta variant because of high transmissibility. However, studies on the impact of the Omicron variant on the severity of COVID-19 are still limited in developing countries. Our study aimed to determine the prognostic factors for the outcomes of patients infected with SARS-CoV-2 Omicron and Delta variants, including age, sex, comorbidities, and smoking. METHODS: In this retrospective cross-sectional study, we involved 352 patients with COVID-19 from Yogyakarta and Central Java provinces, Indonesia, from May 2021 to February 2022, consisting of 164 males and 188 females. We included all patients with the PCR's Ct value of less than 30 for further whole-genome sequencing. RESULTS: Ct value and mean age of COVID-19 patients were not significantly different between both groups (p = 0.146 and 0.273, respectively). Patients infected with Omicron (n = 139) and Delta (n = 213) variants showed similar hospitalization (p = 0.396) and mortality rates (p = 0.565). Multivariate analysis of both groups showed that older age (≥ 65 years) had a higher risk for hospitalization (OR = 3.86 [95% CI = 1.29-11.5]; p = 0.015) and fatalities (OR = 3.91 [95% CI = 1.35-11.42]; p = 0.012). In both groups, patients with cardiovascular disease had a higher risk for hospitalization (OR = 5.36 [95% CI = 1.08-26.52]; p = 0.039), whereas patients with diabetes revealed a higher risk for fatalities (OR = 9.47 [95% CI = 3.23-27.01]; p = < 0.001). CONCLUSIONS: Our study shows that patients infected with Omicron and Delta variants reveal similar clinical outcomes, including hospitalization and mortality. Our findings further confirm that older age, cardiovascular disease, and diabetes are substantial prognostic factors for the outcomes of COVID-19 patients. Our findings imply that COVID-19 patients with older age, cardiovascular disease, or diabetes should be treated comprehensively and cautiously to prevent further morbidity and mortality. Furthermore, incomplete data on vaccination status hampered us from analyzing further its impact on hospitalization and mortality in our patients.

Screening antibiofilm activity of invasive plants growing at the Slope Merapi Mountain, Central Java, against Candida albicans.

BACKGROUND: Candida albicans causes high-mortality candidiasis. Antifungal drug resistance demands the development of virulence factor-targeting drugs, particularly antibiofilm. This study screened the effects of five invasive plants growing in Indonesia (Mimosa pudica, Lantana camara, Acacia mangium, Ageratina riparia, and Mikania micrantha) against C. albicans biofilms. Antifungal activity, antiphospholipase activity, biofilm morphology of C. albicans, and cytotoxic capacity were also evaluated. METHODS: Maceration was used to extract the plants, and the most active extract inhibiting the biofilms was fractionated using liquid-liquid fractionation. Antibiofilm activity was determined by a colorimetric assay, MTT. Antifungal activity was tested using the broth microdilution method. A phospholipase assay was performed using the egg-yolk agar method. Influence on the C. albicans morphology was assessed using scanning electron microscopy (SEM). The cytotoxic effect was carried out against Vero and HeLa cell lines. RESULTS: M. pudica extracts showed the most potent antifungal efficacy with minimum inhibitory concentration (MIC) of 15.62 µg/mL and 7.81 µg/mL for aerial parts and roots, respectively. At high concentrations (500 µg/mL and 250 µg/mL), ethanol extract of M. pudica aerial parts strongly inhibited the phospholipase activity. Ethyl-acetate fraction of M. pudica aerial parts demonstrated the most potent antibiofilm activity against 24 h old biofilm of C. albicans with an inhibitory concentration (53.89%) of 62.5 µg/mL showed no cytotoxicity in both Vero and HeLa cells. This fraction affected the morphology of C. albicans and contained promising compounds for inhibiting the 24 h old biofilm of C. albicans. CONCLUSIONS: Invasive M. pudica plant inhibited the growth of planktonic C. albicans cells and its ethyl acetate fraction decreased the metabolic activity of C. albicans biofilms. This result demonstrates the potential of invasive M. pudica plant to reduce biofilm-associated candida infection.

Antibiofilm Activity of Invasive Plants against Candida albicans: Focus on Baccharis halimifolia Essential Oil and Its Compounds.

The extracts of five invasive plants were investigated for antifungal and antibiofilm activities against Candida albicans, C. glabrata, C. krusei, and C. parapsilosis. The antifungal activity was evaluated using the microdilution assay and the antibiofilm effect by measurement of the metabolic activity. Ethanol and ethanol-water extracts of Reynoutria japonica leaves inhibited 50 % of planktonic cells at 250 µg mL-1 and 15.6 µg mL-1 , respectively. Ethanol and ethanol-water extracts of Baccharis halimifolia inhibited >75 % of the mature biofilm of C. albicans at 500 µg mL-1 . The essential oil (EO) of B. halimifolia leaves was the most active (50 % inhibition (IC50 ) at 4 and 74 µg mL-1 against the maturation phase and 24 h old-biofilms of C. albicans, respectively). Oxygenated sesquiterpenes were the primary contents in this EO (62.02 %), with ß-caryophyllene oxide as the major component (37 %). Aromadendrene oxide-(2), ß-caryophyllene oxide, and (±)-ß-pinene displayed significant activities against the maturation phase (IC50 =9-310 µ mol l-1 ) and preformed 24 h-biofilm (IC50 =38-630 µ mol l-1 ) of C. albicans with very low cytotoxicity for the first two compounds. C. albicans remained the most susceptible species to this EO and its components. This study highlighted for the first time the antibiofilm potential of B. halimifolia, its EO and some of its components.

The high mutation rate at the D614G hotspot-furin cleavage site region increases the priming efficiency of the Spike protein by furin protease: analysis of Indonesian SARS-CoV-2 G614 variants obtained during the early COVID-19 pandemic.

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.

Characteristics of clinical isolates of nontuberculous mycobacteria in Java-Indonesia: A multicenter study.

BACKGROUND: Nontuberculous mycobacterial (NTM) lung infections are a major public health concern. Diagnosis of NTM-pulmonary disease (NTM-PD) is difficult because its clinical, microbiological, and radiological features resemble to those of pulmonary tuberculosis (TB), leading to misdiagnosis. Identification at the species level is essential for diagnosis and determination of therapy, which is currently not performed routinely in Indonesian laboratories. METHODOLOGY AND PRINCIPAL FINDINGS: From January 2020 to May 2021, 94 NTM isolates were collected from three TB referral centers in Java Province. Species were identified using matrix-assisted laser desorption-ionization time-of-flight mass spectrometry (MALDI-TOF MS). Tests were performed to determine antibiotic susceptibility, biofilm formation ability, sliding motility characteristics, and the ability to adhere to and invade pneumocytes. After identifying the species of all the isolates, we found nine groups of NTMs: M. fortuitum group 51% (48/94), M. abscessus 38.3% (36/94), M. intracellulare 3.1% (3/94), M. neoaurum 2.1% (2/94), M. chelonae 1.1% (1/94), M. gordonae 1.1% (1/94), M. szulgai 1.1% (1/94), M. mucogenicum 1.1% (1/94), and M. arupense 1.1% (1/94). Amikacin was the most effective antibiotic against M. fortuitum group and M. abscessus. The M. fortuitum group was significantly better at forming biofilms than M. abscessus, but both had the same sliding motility capability. The ability of the M. fortuitum group to adhere to and invade pneumocytes was better than that of M. abscessus, with the number isolates of the M. fortuitum group capable of superior adhesion and invasion to that of M. abscessus. CONCLUSIONS/SIGNIFICANCE: This study shows that M. fortuitum group and M. abscessus were the most common NTM found in Java, Indonesia. The M. fortuitum group and M. abscessus were the most susceptible to amikacin; therefore, this was the empirical treatment of choice. The ability to form biofilms is directly proportional to the ability to adhere to and invade pneumocytes but not to the susceptibility profile or sliding motility characteristics.

The feasibility of SARS-CoV-2 surveillance using wastewater and environmental sampling in Indonesia.

BACKGROUND: Wastewater-based epidemiology (WBE) surveillance as an early warning system (EWS) for monitoring community transmission of SARS-CoV-2 in low- and middle-income country (LMIC) settings, where diagnostic testing capacity is limited, needs further exploration. We explored the feasibility to conduct a WBE surveillance in Indonesia, one of the global epicenters of the COVID-19 pandemic in the middle of 2021, with the fourth largest population in the world where sewer and non-sewered sewage systems are implemented. The feasibility and resource capacity to collect samples on a weekly or fortnightly basis with grab and/or passive sampling methods, as well as to conduct qualitative and quantitative identification of SARS-CoV-2 ribonucleic acid (RNA) using real-time RT-PCR (RT-qPCR) testing of environmental samples were explored. MATERIALS AND METHODS: We initiated a routine surveillance of wastewater and environmental sampling at three predetermined districts in Special Region of Yogyakarta Province. Water samples were collected from central and community wastewater treatment plants (WWTPs), including manholes flowing to the central WWTP, and additional soil samples were collected for the near source tracking (NST) locations (i.e., public spaces where people congregate). RESULTS: We began collecting samples in the Delta wave of the COVID-19 pandemic in Indonesia in July 2021. From a 10-week period, 54% (296/544) of wastewater and environmental samples were positive for SARS-CoV-2 RNA. The sample positivity rate decreased in proportion with the reported incidence of COVID-19 clinical cases in the community. The highest positivity rate of 77% in week 1, was obtained for samples collected in July 2021 and decreased to 25% in week 10 by the end of September 2021. CONCLUSION: A WBE surveillance system for SARS-CoV-2 in Indonesia is feasible to monitor the community burden of infections. Future studies testing the potential of WBE and EWS for signaling early outbreaks of SARS-CoV-2 transmissions in this setting are required.

Evolutionary dynamics of SARS-CoV-2 circulating in Yogyakarta and Central Java, Indonesia: sequence analysis covering furin cleavage site (FCS) region of the spike protein / Dinámica evolutiva del SARS-CoV-2 que circula en Yogyakarta y Java Central, Indonesia: análisis de secuencia que cubre la región del sitio de escisión de furina (FCS) de la proteína espiga

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.(AU)

Co-infection of SARS-CoV-2 with other viral respiratory pathogens in Yogyakarta, Indonesia: A cross-sectional study.

Background: Growing evidence shows that viral co-infection is found repeatedly in patients with Coronavirus Disease-2019 (COVID-19). This is the first report of SARS-CoV-2 co-infection with viral respiratory pathogens in Indonesia. Methods: Over a one month period of April to May 2020, SARS-CoV-2 positive nasopharyngeal swabs in our COVID-19 referral laboratory in Yogyakarta, Indonesia, were tested for viral respiratory pathogens by real-time, reverse transcription polymerase chain reaction (RT-PCR). Proportion of co-infection reported in percentage. Results: Fifty-nine samples were positive for other viral respiratory pathogens among a total of 125 samples. Influenza A virus was detected in 32 samples, Influenza B in 16 samples, Human metapneumovirus in 1 sample, and adenovirus in 10 samples. We did not detect any co-infection with respiratory syncytial virus. Nine (7.2%) patients had co-infection with more than two viruses. Conclusion: Viral co-infection with SARS-CoV-2 is common. These results will provide a helpful reference for diagnosis and clinical treatment of patients with COVID-19.

Evolutionary dynamics of SARS-CoV-2 circulating in Yogyakarta and Central Java, Indonesia: sequence analysis covering furin cleavage site (FCS) region of the spike protein.

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.

Performance of Xpert MTB/RIF and sputum microscopy compared to sputum culture for diagnosis of tuberculosis in seven hospitals in Indonesia.

Introduction: Tuberculosis (TB) is a major public health concern in Indonesia, where the incidence was 301 cases per 100,000 inhabitants in 2020 and the prevalence of multi-drug resistant (MDR) TB is increasing. Diagnostic testing approaches vary across Indonesia due to resource limitations. Acid-fast bacilli (AFB) smear is widely used, though Xpert MTB/RIF has been the preferred assay for detecting TB and rifampicin resistance since 2012 due to higher sensitivity and ability to rapidly identify rifampicin resistance. However, <1,000 Xpert instruments were available in Indonesia as of 2020 and the Xpert supply chain has suffered interruptions. Methods: We compared the performance of Xpert MTB/RIF and AFB smear to facilitate optimization of TB case identification. We analyzed baseline data from a cohort study of adults with pulmonary TB conducted at seven hospitals across Indonesia. We evaluated sensitivity and specificity of AFB smear and Xpert MTB/RIF using Mycobacterium tuberculosis (Mtb) culture as the gold standard, factors associated with assay results, and consistency of Xpert MTB/RIF with drug susceptibility test (DST) in detecting rifampicin resistance. Results: Sensitivity of AFB smear was significantly lower than Xpert MTB/RIF (86.2 vs. 97.4%, p-value <0.001), but specificity was significantly better (86.7 vs. 73.3%, p-value <0.001). Performance varied by hospital. Positivity rate for AFB smear and Mtb culture was higher in subjects with pulmonary cavities and in morning sputum samples. Consistency of Xpert MTB/RIF with DST was lower in those with rifampicin- sensitive TB by DST. Discussion: Additional evaluation using sputa from primary and secondary Indonesian health centers will increase the generalizability of the assessment of AFB smear and Xpert MTB/RIF performance, and better inform health policy. Clinical trial registration: [https://clinicaltrials.gov/], identifier [NCT027 58236].

Is the Infection of the SARS-CoV-2 Delta Variant Associated With the Outcomes of COVID-19 Patients?

Background: Severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) Delta variant (B.1.617.2) has been responsible for the current increase in Coronavirus disease 2019 (COVID-19) infectivity rate worldwide. We compared the impact of the Delta variant and non-Delta variant on the COVID-19 outcomes in patients from Yogyakarta and Central Java provinces, Indonesia. Methods: In this cross-sectional study, we ascertained 161 patients, 69 with the Delta variant and 92 with the non-Delta variant. The Illumina MiSeq next-generation sequencer was used to perform the whole-genome sequences of SARS-CoV-2. Results: The mean age of patients with the Delta variant and the non-Delta variant was 27.3 ± 20.0 and 43.0 ± 20.9 (p = 3 × 10-6). The patients with Delta variant consisted of 23 males and 46 females, while the patients with the non-Delta variant involved 56 males and 36 females (p = 0.001). The Ct value of the Delta variant (18.4 ± 2.9) was significantly lower than that of the non-Delta variant (19.5 ± 3.8) (p = 0.043). There was no significant difference in the hospitalization and mortality of patients with Delta and non-Delta variants (p = 0.80 and 0.29, respectively). None of the prognostic factors were associated with the hospitalization, except diabetes with an OR of 3.6 (95% CI = 1.02-12.5; p = 0.036). Moreover, the patients with the following factors have been associated with higher mortality rate than the patients without the factors: age ≥65 years, obesity, diabetes, hypertension, and cardiovascular disease with the OR of 11 (95% CI = 3.4-36; p = 8 × 10-5), 27 (95% CI = 6.1-118; p = 1 × 10-5), 15.6 (95% CI = 5.3-46; p = 6 × 10-7), 12 (95% CI = 4-35.3; p = 1.2 × 10-5), and 6.8 (95% CI = 2.1-22.1; p = 0.003), respectively. Multivariate analysis showed that age ≥65 years, obesity, diabetes, and hypertension were the strong prognostic factors for the mortality of COVID-19 patients with the OR of 3.6 (95% CI = 0.58-21.9; p = 0.028), 16.6 (95% CI = 2.5-107.1; p = 0.003), 5.5 (95% CI = 1.3-23.7; p = 0.021), and 5.8 (95% CI = 1.02-32.8; p = 0.047), respectively. Conclusions: We show that the patients infected by the SARS-CoV-2 Delta variant have a lower Ct value than the patients infected by the non-Delta variant, implying that the Delta variant has a higher viral load, which might cause a more transmissible virus among humans. However, the Delta variant does not affect the COVID-19 outcomes in our patients. Our study also confirms that older age and comorbidity increase the mortality rate of patients with COVID-19.

Association between prognostic factors and the outcomes of patients infected with SARS-CoV-2 harboring multiple spike protein mutations.

The outcome of SARS-CoV-2 infection is determined by multiple factors, including the viral, host genetics, age, and comorbidities. This study investigated the association between prognostic factors and disease outcomes of patients infected by SARS-CoV-2 with multiple S protein mutations. Fifty-one COVID-19 patients were recruited in this study. Whole-genome sequencing of 170 full-genomes of SARS-CoV-2 was conducted with the Illumina MiSeq sequencer. Most patients (47%) had mild symptoms of COVID-19 followed by moderate (19.6%), no symptoms (13.7%), severe (4%), and critical (2%). Mortality was found in 13.7% of the COVID-19 patients. There was a significant difference between the age of hospitalized patients (53.4 ± 18 years) and the age of non-hospitalized patients (34.6 ± 19) (p = 0.001). The patients' hospitalization was strongly associated with hypertension, diabetes, and anticoagulant and were strongly significant with the OR of 17 (95% CI 2-144; p = 0.001), 4.47 (95% CI 1.07-18.58; p = 0.039), and 27.97 (95% CI 1.54-507.13; p = 0.02), respectively; while the patients' mortality was significantly correlated with patients' age, anticoagulant, steroid, and diabetes, with OR of 8.44 (95% CI 1.5-47.49; p = 0.016), 46.8 (95% CI 4.63-472.77; p = 0.001), 15.75 (95% CI 2-123.86; p = 0.009), and 8.5 (95% CI 1.43-50.66; p = 0.019), respectively. This study found the clade: L (2%), GH (84.3%), GR (11.7%), and O (2%). Besides the D614G mutation, we found L5F (18.8%), V213A (18.8%), and S689R (8.3%). No significant association between multiple S protein mutations and the patients' hospitalization or mortality. Multivariate analysis revealed that hypertension and anticoagulant were the significant factors influencing the hospitalization and mortality of patients with COVID-19 with an OR of 17.06 (95% CI 2.02-144.36; p = 0.009) and 46.8 (95% CI 4.63-472.77; p = 0.001), respectively. Moreover, the multiple S protein mutations almost reached a strong association with patients' hospitalization (p = 0.07). We concluded that hypertension and anticoagulant therapy have a significant impact on COVID-19 outcomes. This study also suggests that multiple S protein mutations may impact the COVID-19 outcomes. This further emphasized the significance of monitoring SARS-CoV-2 variants through genomic surveillance, particularly those that may impact the COVID-19 outcomes.

Detection of SARS-CoV-2 spike protein D614G mutation by qPCR-HRM analysis.

OBJECTIVES: Monitoring the spread of the G614 in specific locations is critical as this variant is highly transmissible and can trigger the emergence of other mutations. Therefore, a rapid and accurate method that can reliably detect the D614G mutation will be beneficial. This study aims to analyze the potential use of the two-step Reverse Transcriptase quantitative polymerase chain reaction - high resolution melting analysis (RT-qPCR-HRM) to detect a specific mutation in the SARS-CoV-2 genome. METHODS: Six SARS-CoV-2 RNA samples were synthesized into cDNA and analyzed with the qPCR-HRM method in order to detect the D614G mutation in Spike protein of SARS-CoV-2. The primers are designed to target the specific Spike region containing the D614G mutation. The qPCR-HRM analysis was conducted simultaneously, and the identification of the SARS-CoV-2 variant was confirmed by conventional PCR and Sanger sequencing methods. RESULTS: The results showed that the melting temperature (Tm) of the D614 variant was 79.39 ± 0.03 °C, which was slightly lower than the Tm of the G614 variant (79.62 ± 0.015 °C). The results of the HRM analysis, visualized by the normalized melting curve and the difference curve were able to discriminate the D614 and G614 variant samples. All samples were identified as G614 variants by qPCR-HRM assay, which was subsequently confirmed by Sanger sequencing. CONCLUSIONS: This study demonstrated a sensitive method that can identify the D614G mutation by a simple two-step RT-qPCR-HRM assay procedure analysis, which can be useful for active surveillance of the transmission of a specific mutation.

Molecular epidemiology of SARS-CoV-2 isolated from COVID-19 family clusters.

BACKGROUND: Transmission within families and multiple spike protein mutations have been associated with the rapid transmission of SARS-CoV-2. We aimed to: (1) describe full genome characterization of SARS-CoV-2 and correlate the sequences with epidemiological data within family clusters, and (2) conduct phylogenetic analysis of all samples from Yogyakarta and Central Java, Indonesia and other countries. METHODS: The study involved 17 patients with COVID-19, including two family clusters. We determined the full-genome sequences of SARS-CoV-2 using the Illumina MiSeq next-generation sequencer. Phylogenetic analysis was performed using a dataset of 142 full-genomes of SARS-CoV-2 from different regions. RESULTS: Ninety-four SNPs were detected throughout the open reading frame (ORF) of SARS-CoV-2 samples with 58% (54/94) of the nucleic acid changes resulting in amino acid mutations. About 94% (16/17) of the virus samples showed D614G on spike protein and 56% of these (9/16) showed other various amino acid mutations on this protein, including L5F, V83L, V213A, W258R, Q677H, and N811I. The virus samples from family cluster-1 (n = 3) belong to the same clade GH, in which two were collected from deceased patients, and the other from the survived patient. All samples from this family cluster revealed a combination of spike protein mutations of D614G and V213A. Virus samples from family cluster-2 (n = 3) also belonged to the clade GH and showed other spike protein mutations of L5F alongside the D614G mutation. CONCLUSIONS: Our study is the first comprehensive report associating the full-genome sequences of SARS-CoV-2 with the epidemiological data within family clusters. Phylogenetic analysis revealed that the three viruses from family cluster-1 formed a monophyletic group, whereas viruses from family cluster-2 formed a polyphyletic group indicating there is the possibility of different sources of infection. This study highlights how the same spike protein mutations among members of the same family might show different disease outcomes.

High frequency of azole resistant Candida spp. colonization among presumptive multidrug resistant tuberculosis (MDR-TB) patients.

BACKGROUND: Tuberculosis is one of the major causes of death globally. The problems become even more complicated with the rise in prevalence of multidrug resistant tuberculosis (MDR-TB). Many diseases have been reported to occur with tuberculosis making it more difficult to manage. Candida spp., which are yeast-like fungi and a constituent of normal flora in humans, are notoriously reported to be one of the most common opportunistic nosocomial infections. This study aimed to measure the proportion of presumptive MDR-TB patients colonized with Candida spp. and to characterize its susceptibility against azole group antifungal agents. METHODS: Sputum from presumptive MDR-TB patients were collected and examined for the presence of Mycobacterium tuberculosis and its rifampicin resistant status using GeneXpert. It was further cultured on Sabouroud's Dextrose Agar (SDA) to isolate the Candida spp. The Candida species were determined using HiCrome™ Candidal Differential Agar. Antifungal susceptibility was tested using microbroth dilution methods. Checkerboard microdilution assays were performed to measure the interaction between rifampicin and fluconazole to C. albicans. RESULTS: There were 355 presumptive MDR-TB patients enrolled. A total of 101 (28.4%) patients were confirmed to have M. tuberculosis. There were 113 (31.8%) sputum positive for Candida spp., which corresponded to 149 Candida spp. isolates. Candida albicans was the most frequent (53.7%) species isolated from all patients. The susceptibility of Candida spp. against fluconazole, itraconazole, and ketoconazole were 38.3%, 1.3%, and 10.7% respectively. There was significant association between rifampicin exposure history and susceptibility of Candida albicans against fluconazole (Odds Ratio: 9.96; 95% CI: 1.83-54.19; p <0.01), but not for ketoconazole and itraconazole. The checkerboard microdilution assays showed that rifampicin decreased the fungicidal activity of fluconazole to C. albicans in a dose-dependent manner. CONCLUSION: There was high frequency of azole resistant Candida spp. isolates colonizing the respiratory tract of presumptive MDR-TB patients. This presence might indicate the association of chronic exposure to rifampicin, the main drug for tuberculosis therapy, with the induction of azole resistance.

Full-length genome characterization and phylogenetic analysis of SARS-CoV-2 virus strains from Yogyakarta and Central Java, Indonesia.

BACKGROUND: Recently, SARS-CoV-2 virus with the D614G mutation has become a public concern due to rapid dissemination of this variant across many countries. Our study aims were (1) to report full-length genome sequences of SARS-CoV-2 collected from four COVID-19 patients in the Special Region of Yogyakarta and Central Java provinces, Indonesia; (2) to compare the clade distribution of full-length genome sequences from Indonesia (n = 60) from March to September 2020 and (3) to perform phylogenetic analysis of SARS-CoV-2 complete genomes from different countries, including Indonesia. METHODS: Whole genome sequencing (WGS) was performed using next-generation sequencing (NGS) applied in the Illumina MiSeq instrument. Full-length virus genomes were annotated using the reference genome of hCoV-19/Wuhan/Hu-1/2019 (NC_045512.2) and then visualized in UGENE v. 1.30. For phylogenetic analysis, a dataset of 88 available SARS-CoV-2 complete genomes from different countries, including Indonesia, was retrieved from GISAID. RESULTS: All patients were hospitalized with various severities of COVID-19. Phylogenetic analysis revealed that one and three virus samples belong to clade L and GH. These three clade GH virus samples (EPI_ISL_525492, EPI_ISL_516800 and EPI_ISL_516829) were not only located in a cluster with SARS-CoV-2 genomes from Asia but also those from Europe, whereas the clade L virus sample (EPI_ISL_516806) was located amongst SARS-CoV-2 genomes from Asia. Using full-length sequences available in the GISAID EpiCoV Database, 39 of 60 SARS-CoV-2 (65%) from Indonesia harbor the D614G mutation. CONCLUSION: These findings indicate that SARS-CoV-2 with the D614G mutation appears to become the major circulating virus in Indonesia, concurrent with the COVID-19 situation worldwide.