The COVID-19 pandemic in East Java, Indonesia has profoundly affected the region, resulting in a substantial number of cases and elevated mortality rates. ¹ The first confirmed case of COVID-19 in East Java was documented on March 11, 2020, involving a 50-year-old woman who had recently returned from Saudi Arabia and tested positive for the virus. ² Since that time, the incidence of cases has consistently risen, resulting in pressure on the hospital system and public health infrastructure. ³ The COVID-19 pandemic, instigated by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has necessitated the urgent discovery of effective vaccines to control the virus’s transmission and alleviate its effects on public health. ⁴ Indonesia has previously succeeded in developing an inactivated platform SARS-CoV-2 vaccine that has been clinically tested and widely used for the Indonesian population. ^{5– 9} However, there are concerns about the emergence of various new variants that could reduce the efficacy of the current vaccines, necessitating the development of a new vaccine with an mRNA cocktail platform derived from various SARS-CoV-2 variants, which is expected to help combat future COVID-19 outbreaks. Exploration of potential isolates to be developed as vaccine candidates is a crucial and important initial step.

Whole genome sequencing and phylogenetic analysis of SARS-CoV-2 isolates have become essential methodologies for elucidating the genetic diversity, evolution, and prospective vaccine targets of the virus. ^{10– 12} Phylogenetic analysis enhances whole genome sequencing by examining the evolutionary links among various viral strains. ¹³ Researchers can reconstruct the phylogenetic tree by comparing the genetic sequences of SARS-CoV-2 isolates from East Java with those from other places, yielding information about the virus’s genesis, dissemination, and transmission patterns. This information is essential for comprehending the virus’s evolution and dissemination throughout the local population. ^{14, 15} Additionally, whole genome sequencing and phylogenetic analysis can identify conserved sections of the viral genome that are less prone to mutation and may serve as targets for vaccine development. These areas may encompass highly immunogenic epitopes or critical viral proteins implicated in viral replication and host interactions. By focusing on these conserved areas, novel platform vaccines can offer extensive protection against several viral types. ¹⁶

This study aims to conduct whole genome sequencing, phylogenetic analysis, and mutational assessment of SARS-CoV-2 potential isolates from East Java, Indonesia, to find possible next-generation vaccine candidates. The results will enhance comprehension of the genetic variety, evolution, and prospective vaccine targets of the virus in this area. This research is essential for directing the creation of innovative and efficacious vaccines against SARS-CoV-2, especially in regions significantly affected by the disease, such as East Java. Utilizing whole genome sequencing and phylogenetic analysis, we can facilitate the creation of novel vaccine platforms capable of managing the COVID-19 pandemic and future outbreaks of related coronaviruses.

The Delta variant evolved from a variant under investigation (VUI) to a variant of concern (VOC) upon its classification within the B.1.617.2 lineage. The reclassification was predicated on an assessment of its transmissibility, which was demonstrated to be at least comparable to B.1.1.7 (Alpha variant). The mutation is thought to be partially responsible for the commencement of India’s deadly second wave of the pandemic, which began in February 2021. By the end of July, it had also caused an increase in daily infections across various parts of Asia, ²⁰ the United States, ²¹ Australia, and New Zealand. During the Delta epidemic, the initial prevalence of the Delta variant remained below 5% of the total sequenced genomes, as it was outpaced by the Delta sub-lineages, AY.23 and AY.24. Two isolates exhibiting the Delta variant from lineage AY.23 were acquired from the isolation results. These isolates demonstrated vigorous development. The AY.23 mutation emerged in Indonesia and exhibits parallels to the Delta Plus variant identified in Britain. ²² The AY.1 sub-lineage of the delta variant, distinguished by an additional K417N mutation in its spike protein, has been reported in England. The AY.1 variant, also known as the Delta Plus variant, is thought to exhibit the highest mortality rate. The Delta Plus form is thought to exhibit improved antibody evasion due to the K417N mutation. The Beta variant has been previously documented to possess a K417N mutation. Furthermore, relative to other Variant of Concerns (VoCs), the Delta plus variant has an increased propensity for transmission and a pronounced affinity for lung epithelial cells. ²³

Tables 2 and 3 present the nucleotide and amino acid mutations of the structural proteins of selected SARS-CoV-2 isolates. T19R increases the affinity to the ACE2 protein and facilitates its evasion of monoclonal antibodies, including regdanvimab and bamlanivimab, by reducing their binding to the virus. ²⁴ The Delta variant’s unique triple mutation, E156/F157/R158G, was precisely detected in our consensus nucleotide sequences. Nonetheless, the COVID-19 genome annotator program inadequately identified this non-codon-aligned triple mutation. ²⁵ The L452R mutation increases transmission, significantly intensifies severity, and facilitates the evasion of neutralizing antibodies generated by vaccines. ^{26– 28} The T478K mutation increases transmissibility and results in heightened disease severity. Furthermore, it diminishes the virus’s susceptibility to neutralizing antibodies generated by vaccines. ^{29, 30} The P681R mutation increases transmission, does not significantly impact disease severity, and partially diminishes vaccination efficacy by evading neutralizing antibodies. ^{29, 30} The P681R mutation increases transmission, does not significantly impact disease severity, and partially diminishes vaccination efficacy by evading neutralizing antibodies. ³¹ The D614G mutation substantially enhances transmission rates, exerts no significant impact on disease severity, and markedly diminishes vaccine efficacy. The vaccine resulted in a 3-6-fold reduction in the titer of serum neutralizing antibodies (NABs) against the pseudovirus. ^{32, 33}

The D950N mutation in the spike protein of Delta variant was identified in a region external of the receptor-binding domain (RBD), which is recognized for its stability and resistance to recurrent mutations. ³⁴ The V1264L variant in the spike protein introduces an acidic dileucine motif. This motif may affect the endocytosis process for several receptors associated with the resulting protein. ^{35– 37} V1264L possesses the capability to enhance the functionality of the spike protein. The presence of D63G was associated with a notable increase in viral load and may enable the virus to evade the immune system during replication. ³⁸

The mutation in the M gene confers a greater biological fitness benefit, potentially associated with glucose absorption during viral replication. Consequently, it is essential to integrate this lineage into current genomic surveillance efforts and meticulously evaluate its possible effects on heightened pathogenicity and treatment considerations. ³⁹ The mutations D63G, R203M, G215C, and D377Y in the N gene were associated with the highest incidence of intensive care unit admissions attributed to the Delta variation. Alterations were detected in an instance of breakthrough reinfection, in which the patient had hypoxia and necessitated hospitalization. These mutations were significantly associated with mortality, as detailed in a study. ⁴⁰ These modifications may facilitate the propagation of the virus.

Numerous Omicron lineages have been identified following the classification of B.1.1.529 as a Variant of Concern (VOC) on November 26, 2021. The Omicron variant has been categorized into several sub-lineages, one of which is BA.11. The Omicron sub-variants demonstrate an increased transmissibility relative to the Omicron variant (BA.1) and Delta. Nonetheless, an Omicron infection is less severe than a Delta infection. The BA.5 subvariant of Omicron is presently the most transmissible subvariant and has emerged as the predominant variant in circulation in the United States and numerous other regions globally. BA.5 is gradually replacing the original BA.1 and BA.2 Omicron subvariants of SARS-CoV-2 in several countries. ⁴⁰ The Omicron BQ.1 subvariant is swiftly becoming as the predominant variant in numerous countries globally, despite its recent emergence. Initial research indicates that BQ.1 is a subvariant of BA.5, with many progeny including BQ.1.1, BQ.1.2, BQ.1.3, and BQ.1.4.2. Alongside BQ.1, its progeny BQ.1.1 is also being monitored with increasing prevalence globally. The distribution of BQ1.23 is as follows: the United States of America comprises 26.0%, Australia 12.0%, Indonesia 9.0%, Japan 7.0%, and Germany 5.0%. ⁴¹

We observed that the isolated SARS-CoV-2 omicron variant exhibited a mutation at T91 in the E protein. The T9I mutation results in the formation of a nonselective ion channel with reduced sensitivity to acid. T9I also decreased cytokine production and mitigated cell death. The alterations in channel properties may account for the Omicron variant’s diminished efficacy and lower extent of induced cellular damage. The severity of ionic dyshomeostasis, membrane disruption, and lysosome-related cell death would be mitigated by the enhanced chloride influx and reduced potassium efflux facilitated by T9I channels. ⁴²

Alterations Q19E and A63T were identified in all major Omicron subvariants. The closeness of the A63T mutation to a key location indicates that it may influence the stabilization of the M protein dimer. ⁴² The N-terminal D3G and D3N mutations were solely identified in BA.1 and BA.5, respectively, possibly resulting in the establishment of the N-myristoylation site at positions 3–8. ⁴³ The P13L mutation may enable the virus to evade cellular immunity. ⁴⁴ R203K mutations augment the transmission and severity of particular SARS-CoV-2 variants. ⁴⁵ The Alpha and Gamma VOCs display modifications identical to the R203K and G204R mutations in the N protein, which may result in increased viral loads and augmented expression of subgenomic RNA. ^{46, 47}

It was determined that several isolates are available for the development of next-generation vaccines targeting distinct SARS-CoV-2 variants, including Alpha, Delta, and Omicron. Sequencing test results for nine isolates indicated that four samples were identified as the Alpha (2) and Delta (2) variants, lineage B.1 and AY.23, respectively, while five samples were identified as the Omicron variant, with lineages BA.1.1, BA.5.2, and BQ.1.23.

Ethics and consent were not required.