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.

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.

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.

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.

Prevalence and Distribution of Dengue Virus in Aedes aegypti in Yogyakarta City before Deployment of Wolbachia Infected Aedes aegypti.

Indonesia is one of the countries where dengue infection is prevalent. In this study we measure the prevalence and distribution of dengue virus (DENV) DENV-infected Aedes aegypti in Yogyakarta City, Indonesia, during the wet season when high dengue transmission period occurred, as baseline data before implementation of a Wolbachia-infected Aedes aegypti trial for dengue control. We applied One-Step Multiplex Real Time PCR (RT-PCR) for the type-specific-detection of dengue viruses in field-caught adult Aedes aegypti mosquitoes. In a prospective field study conducted from December 2015 to May 2016, adult female Aedes aegypti were caught from selected areas in Yogyakarta City, and then screened by using RT-PCR. During the survey period, 36 (0.12%) mosquitoes from amongst 29,252 female mosquitoes were positive for a DENV type. In total, 22.20% of dengue-positive mosquitoes were DENV-1, 25% were DENV-2, 17% were DENV-3, but none were positive for DENV-4. This study has provided dengue virus infection prevalence in field-caught Aedes aegypti and its circulating serotype in Yogyakarta City before deployment of Wolbachia-infected Aedes aegypti.

Anti-SIRPα antibodies as a potential new tool for cancer immunotherapy.

Tumor cells are thought to evade immune surveillance through interaction with immune cells. Much recent attention has focused on the modification of immune responses as a basis for new cancer treatments. SIRPα is an Ig superfamily protein that inhibits phagocytosis in macrophages upon interaction with its ligand CD47 expressed on the surface of target cells. Here, we show that SIRPα is highly expressed in human renal cell carcinoma and melanoma. Furthermore, an anti-SIRPα Ab that blocks the interaction with CD47 markedly suppressed tumor formation by renal cell carcinoma or melanoma cells in immunocompetent syngeneic mice. This inhibitory effect of the Ab appeared to be mediated by dual mechanisms: direct induction of Ab-dependent cellular phagocytosis of tumor cells by macrophages and blockade of CD47-SIRPα signaling that negatively regulates such phagocytosis. The antitumor effect of the Ab was greatly attenuated by selective depletion not only of macrophages but also of NK cells or CD8+ T cells. In addition, the anti-SIRPα Ab also enhances the inhibitory effects of Abs against CD20 and programmed cell death 1 (PD-1) on tumor formation in mice injected with SIRPα-nonexpressing tumor cells. Anti-SIRPα Abs thus warrant further study as a potential new therapy for a broad range of cancers.

Protein tyrosine phosphatase SAP-1 protects against colitis through regulation of CEACAM20 in the intestinal epithelium.

Intestinal epithelial cells contribute to regulation of intestinal immunity in mammals, but the detailed molecular mechanisms of such regulation have remained largely unknown. Stomach-cancer-associated protein tyrosine phosphatase 1 (SAP-1, also known as PTPRH) is a receptor-type protein tyrosine phosphatase that is localized specifically at microvilli of the brush border in gastrointestinal epithelial cells. Here we show that SAP-1 ablation in interleukin (IL)-10-deficient mice, a model of inflammatory bowel disease, resulted in a marked increase in the severity of colitis in association with up-regulation of mRNAs for various cytokines and chemokines in the colon. Tyrosine phosphorylation of carcinoembryonic antigen-related cell adhesion molecule (CEACAM) 20, an intestinal microvillus-specific transmembrane protein of the Ig superfamily, was greatly increased in the intestinal epithelium of the SAP-1-deficient animals, suggesting that this protein is a substrate for SAP-1. Tyrosine phosphorylation of CEACAM20 by the protein tyrosine kinase c-Src and the consequent association of CEACAM20 with spleen tyrosine kinase (Syk) promoted the production of IL-8 in cultured cells through the activation of nuclear factor-κB (NF-κB). In addition, SAP-1 and CEACAM20 were found to form a complex through interaction of their ectodomains. SAP-1 and CEACAM20 thus constitute a regulatory system through which the intestinal epithelium contributes to intestinal immunity.

Tyrosine Phosphorylation of Carcinoembryonic Antigen-related Cell Adhesion Molecule 20 and Its Functional Role.

Carcinoembryonic antigen-related cell adhesion molecule (CEACAM) 20 is an immunoglobulin-superfamily transmembrane protein that contains an immunoreceptor tyrosine-based activation motif (ITAM) in its cytoplasmic region. However, the mechanism for tyrosine phosphorylation of, or the physiological function of, this protein remains largely unknown. Here we have shown that CEACAM20 is indeed tyrosine-phosphorylated by either treatment with pervanadate or forced expression of c-Src. In addition, Tyr522, Tyr559 or Tyr570, the latter two of which are within the ITAM, is likely important for such tyrosine phosphorylation. Forced expression of Myc-tagged wild-type CEACAM20 promoted the phagocytic activity of cultured cells for microbeads coupled with anti-Myc antibodies. By contrast, such phagocytic activity was markedly reduced when a mutant form of CEACAM20, in which Tyr559 and Tyr570 were substituted with phenylalanine, was expressed. Furthermore, the CEACAM20-mediated phagocytic activity was markedly prevented by the treatment with an inhibitor for either Src family kinases (SFKs), Syk, phosphoinositide 3-kinase (PI3K) or phospholipase C-γ (PLCγ). Inhibition of actin polymerization by Cytochalasin D significantly inhibited the CEACAM20-mediated phagocytosis. These results thus suggest that tyrosine phosphorylation of CEACAM20 likely promotes phagocytic activity of the cells. The CEACAM20-mediated phagocytic activity requires the activation of SFKs, Syk, PI3K or PLCγ.