Epistasis facilitates functional evolution in an ancient transcription factor.

A protein's genetic architecture - the set of causal rules by which its sequence produces its functions - also determines its possible evolutionary trajectories. Prior research has proposed that the genetic architecture of proteins is very complex, with pervasive epistatic interactions that constrain evolution and make function difficult to predict from sequence. Most of this work has analyzed only the direct paths between two proteins of interest - excluding the vast majority of possible genotypes and evolutionary trajectories - and has considered only a single protein function, leaving unaddressed the genetic architecture of functional specificity and its impact on the evolution of new functions. Here, we develop a new method based on ordinal logistic regression to directly characterize the global genetic determinants of multiple protein functions from 20-state combinatorial deep mutational scanning (DMS) experiments. We use it to dissect the genetic architecture and evolution of a transcription factor's specificity for DNA, using data from a combinatorial DMS of an ancient steroid hormone receptor's capacity to activate transcription from two biologically relevant DNA elements. We show that the genetic architecture of DNA recognition consists of a dense set of main and pairwise effects that involve virtually every possible amino acid state in the protein-DNA interface, but higher-order epistasis plays only a tiny role. Pairwise interactions enlarge the set of functional sequences and are the primary determinants of specificity for different DNA elements. They also massively expand the number of opportunities for single-residue mutations to switch specificity from one DNA target to another. By bringing variants with different functions close together in sequence space, pairwise epistasis therefore facilitates rather than constrains the evolution of new functions.

Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae.

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants.

The promise and pitfalls of synteny in phylogenomics.

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Installation of LYRM proteins in early eukaryotes to regulate the metabolic capacity of the emerging mitochondrion.

Core mitochondrial processes such as the electron transport chain, protein translation and the formation of Fe-S clusters (ISC) are of prokaryotic origin and were present in the bacterial ancestor of mitochondria. In animal and fungal models, a family of small Leu-Tyr-Arg motif-containing proteins (LYRMs) uniformly regulates the function of mitochondrial complexes involved in these processes. The action of LYRMs is contingent upon their binding to the acylated form of acyl carrier protein (ACP). This study demonstrates that LYRMs are structurally and evolutionarily related proteins characterized by a core triplet of α-helices. Their widespread distribution across eukaryotes suggests that 12 specialized LYRMs were likely present in the last eukaryotic common ancestor to regulate the assembly and folding of the subunits that are conserved in bacteria but that lack LYRM homologues. The secondary reduction of mitochondria to anoxic environments has rendered the function of LYRMs and their interaction with acylated ACP dispensable. Consequently, these findings strongly suggest that early eukaryotes installed LYRMs in aerobic mitochondria as orchestrated switches, essential for regulating core metabolism and ATP production.

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Dynamics of karyotype evolution.

In the evolution of species, the karyotype changes with a timescale of tens to hundreds of thousand years. In the development of cancer, the karyotype often is modified in cancerous cells over the lifetime of an individual. Characterizing these changes and understanding the mechanisms leading to them has been of interest in a broad range of disciplines including evolution, cytogenetics, and cancer genetics. A central issue relates to the relative roles of random vs deterministic mechanisms in shaping the changes. Although it is possible that all changes result from random events followed by selection, many results point to other non-random factors that play a role in karyotype evolution. In cancer, chromosomal instability leads to characteristic changes in the karyotype, in which different individuals with a specific type of cancer display similar changes in karyotype structure over time. Statistical analyses of chromosome lengths in different species indicate that the length distribution of chromosomes is not consistent with models in which the lengths of chromosomes are random or evolve solely by simple random processes. A better understanding of the mechanisms underlying karyotype evolution should enable the development of quantitative theoretical models that combine the random and deterministic processes that can be compared to experimental determinations of the karyotype in diverse settings.

Characterizing the antigenic evolution of pandemic influenza A (H1N1) pdm09 from 2009 to 2023.

The H1N1pdm09 virus has been a persistent threat to public health since the 2009 pandemic. Particularly, since the relaxation of COVID-19 pandemic mitigation measures, the influenza virus and SARS-CoV-2 have been concurrently prevalent worldwide. To determine the antigenic evolution pattern of H1N1pdm09 and develop preventive countermeasures, we collected influenza sequence data and immunological data to establish a new antigenic evolution analysis framework. A machine learning model (XGBoost, accuracy = 0.86, area under the receiver operating characteristic curve = 0.89) was constructed using epitopes, physicochemical properties, receptor binding sites, and glycosylation sites as features to predict the antigenic similarity relationships between influenza strains. An antigenic correlation network was constructed, and the Markov clustering algorithm was used to identify antigenic clusters. Subsequently, the antigenic evolution pattern of H1N1pdm09 was analyzed at the global and regional scales across three continents. We found that H1N1pdm09 evolved into around five antigenic clusters between 2009 and 2023 and that their antigenic evolution trajectories were characterized by cocirculation of multiple clusters, low-level persistence of former dominant clusters, and local heterogeneity of cluster circulations. Furthermore, compared with the seasonal H1N1 virus, the potential cluster-transition determining sites of H1N1pdm09 were restricted to epitopes Sa and Sb. This study demonstrated the effectiveness of machine learning methods for characterizing antigenic evolution of viruses, developed a specific model to rapidly identify H1N1pdm09 antigenic variants, and elucidated their evolutionary patterns. Our findings may provide valuable support for the implementation of effective surveillance strategies and targeted prevention efforts to mitigate the impact of H1N1pdm09.

Molecular characterization of emerging Echovirus 11 (E11) shed light on the recombinant origin of a variant associated with severe hepatitis in neonates.

Echovirus 11 (E11) has gained attention owing to its association with severe neonatal infections. Due to the limited data available, the World Health Organization (WHO) considers public health risk to the general population to be low. The present study investigated the genetic variation and molecular evolution of E11 genomes collected from May to December 2023. Whole genome sequencing (WGS) was performed for 16 E11 strains. Phylogenetic analysis on WG showed how all Italian strains belonged to genogroup D5, similarly to other E11 strains recently reported in France and Germany all together aggregated into separate clusters. A cluster-specific recombination pattern was also identified using phylogenetic analysis of different genome regions. Echovirus 6 was identified as the major recombinant virus in 3Cpro and 3Dpol regions. The molecular clock analysis revealed that the recombination event probably occurred in June 2018 (95% HPD interval: Jan 2016-Jan 2020). Shannon entropy analyses, within P1 region, showed how 11 amino acids exhibited relatively high entropy. Five of them were exposed on the canyon region which is responsible for receptor binding with the neonatal Fc receptor. The present study showed the recombinant origin of a new lineage of E11 associated with severe neonatal infections.

The emergence and circulation of human immunodeficiency virus (HIV)-1 subtype C.

Introduction. Human immunodeficiency virus (HIV)-1 subtype C is the most prevalent globally and is thought to have originated in non-human primates in the Democratic Republic of Congo.Hypothesis/Gap Statement. Although the global dominance of HIV-1 subtype C is well established, a thorough understanding of its evolutionary history and transmission dynamics across various risk populations remains elusive. The current knowledge is insufficient to fully capture the global diversification and dissemination of this subtype.Aim. We for the first time sought to investigate the global evolutionary history and spatiotemporal dynamics of HIV-1 subtype C using a selection of maximum-likelihood-based phylodynamic approaches on a total of 1210 near full-length genomic sequences sampled from 32 countries, collected in 4 continents, with sampling dates between 1986-2019 among various risk groups were analysed.Methodology. We subsampled the HIV-1 subtype C genomic datasets based on continent and risk group traits, and performed nucleotide substitution model selection analysis, maximum likelihood (ML) phylogenetic reconstruction, phylogenetic tree topology similarity analysis, temporal signal analysis and traced the timings of viral spread both geographically and by risk group.Results. Based on the phylodynamic analyses of four datasets (full1210, locrisk626, loc562 and risk393), we inferred the time to the most recent common ancestor (TMRCA) in the 1930s and an evolutionary rate of 0.0023 substitutions per site per year. The total number of introduction events of HIV-1 subtype C between continents and between risk groups is estimated to be 71 and 115, respectively. The largest number of introductions occurred from Africa to Europe (n=32), from not-recorded to heterosexual (n=40) and from heterosexual to not-recorded (n=51) risk groups.Conclusion. Our results emphasize that HIV subtype C has mainly spread from Africa to Europe, likely through heterosexual transmission.

Solvent constraints for biopolymer folding and evolution in extraterrestrial environments.

We propose that spontaneous folding and molecular evolution of biopolymers are two universal aspects that must concur for life to happen. These aspects are fundamentally related to the chemical composition of biopolymers and crucially depend on the solvent in which they are embedded. We show that molecular information theory and energy landscape theory allow us to explore the limits that solvents impose on biopolymer existence. We consider 54 solvents, including water, alcohols, hydrocarbons, halogenated solvents, aromatic solvents, and low molecular weight substances made up of elements abundant in the universe, which may potentially take part in alternative biochemistries. We find that along with water, there are many solvents for which the liquid regime is compatible with biopolymer folding and evolution. We present a ranking of the solvents in terms of biopolymer compatibility. Many of these solvents have been found in molecular clouds or may be expected to occur in extrasolar planets.

The physical and evolutionary energy landscapes of devolved protein sequences corresponding to pseudogenes.

Protein evolution is guided by structural, functional, and dynamical constraints ensuring organismal viability. Pseudogenes are genomic sequences identified in many eukaryotes that lack translational activity due to sequence degradation and thus over time have undergone "devolution." Previously pseudogenized genes sometimes regain their protein-coding function, suggesting they may still encode robust folding energy landscapes despite multiple mutations. We study both the physical folding landscapes of protein sequences corresponding to human pseudogenes using the Associative Memory, Water Mediated, Structure and Energy Model, and the evolutionary energy landscapes obtained using direct coupling analysis (DCA) on their parent protein families. We found that generally mutations that have occurred in pseudogene sequences have disrupted their native global network of stabilizing residue interactions, making it harder for them to fold if they were translated. In some cases, however, energetic frustration has apparently decreased when the functional constraints were removed. We analyzed this unexpected situation for Cyclophilin A, Profilin-1, and Small Ubiquitin-like Modifier 2 Protein. Our analysis reveals that when such mutations in the pseudogene ultimately stabilize folding, at the same time, they likely alter the pseudogenes' former biological activity, as estimated by DCA. We localize most of these stabilizing mutations generally to normally frustrated regions required for binding to other partners.

Interactions between chaperone and energy storage networks during the evolution of Legionella pneumophila under heat shock.

Waterborne transmission of the bacterium Legionella pneumophila has emerged as a major cause of severe nosocomial infections of major public health impact. The major route of transmission involves the uptake of aerosolized bacteria, often from the contaminated hot water systems of large buildings. Public health regulations aimed at controlling the mesophilic pathogen are generally concerned with acute pasteurization and maintaining high temperatures at the heating systems and throughout the plumbing of hot water systems, but L. pneumophila is often able to survive these treatments due to both bacterium-intrinsic and environmental factors. Previous work has established an experimental evolution system to model the observations of increased heat resistance in repeatedly but unsuccessfully pasteurized L. pneumophila populations. Here, we show rapid fixation of novel alleles in lineages selected for resistance to heat shock and shifts in mutational profile related to increases in the temperature of selection. Gene-level and nucleotide-level parallelisms between independently-evolving lineages show the centrality of the DnaJ/DnaK chaperone system in the heat resistance of L. pneumophila. Inference of epistatic interactions through reverse genetics shows an unexpected interaction between DnaJ/DnaK and the polyhydroxybutyrate-accumulation energy storage mechanism used by the species to survive long-term starvation in low-nutrient environments.

Evolutionary genomics of socially polymorphic populations of Pogonomyrmex californicus.

BACKGROUND: Social insects vary considerably in their social organization both between and within species. In the California harvester ant, Pogonomyrmex californicus (Buckley 1867), colonies are commonly founded and headed by a single queen (haplometrosis, primary monogyny). However, in some populations in California (USA), unrelated queens cooperate not only during founding (pleometrosis) but also throughout the life of the colony (primary polygyny). The genetic architecture and evolutionary dynamics of this complex social niche polymorphism (haplometrosis vs pleometrosis) have remained unknown. RESULTS: We provide a first analysis of its genomic basis and evolutionary history using population genomics comparing individuals from a haplometrotic population to those from a pleometrotic population. We discovered a recently evolved (< 200 k years), 8-Mb non-recombining region segregating with the observed social niche polymorphism. This region shares several characteristics with supergenes underlying social polymorphisms in other socially polymorphic ant species. However, we also find remarkable differences from previously described social supergenes. Particularly, four additional genomic regions not in linkage with the supergene show signatures of a selective sweep in the pleometrotic population. Within these regions, we find for example genes crucial for epigenetic regulation via histone modification (chameau) and DNA methylation (Dnmt1). CONCLUSIONS: Altogether, our results suggest that social morph in this species is a polygenic trait involving a potential young supergene. Further studies targeting haplo- and pleometrotic individuals from a single population are however required to conclusively resolve whether these genetic differences underlie the alternative social phenotypes or have emerged through genetic drift.

Similar enzymatic functions in distinct bioluminescence systems: evolutionary recruitment of sulfotransferases in ostracod light organs.

Genes from ancient families are sometimes involved in the convergent evolutionary origins of similar traits, even across vast phylogenetic distances. Sulfotransferases are an ancient family of enzymes that transfer sulfate from a donor to a wide variety of substrates, including probable roles in some bioluminescence systems. Here, we demonstrate multiple sulfotransferases, highly expressed in light organs of the bioluminescent ostracod Vargula tsujii, transfer sulfate in vitro to the luciferin substrate, vargulin. We find luciferin sulfotransferases (LSTs) of ostracods are not orthologous to known LSTs of fireflies or sea pansies; animals with distinct and convergently evolved bioluminescence systems compared to ostracods. Therefore, distantly related sulfotransferases were independently recruited at least three times, leading to parallel evolution of luciferin metabolism in three highly diverged organisms. Reuse of homologous genes is surprising in these bioluminescence systems because the other components, including luciferins and luciferases, are completely distinct. Whether convergently evolved traits incorporate ancient genes with similar functions or instead use distinct, often newer, genes may be constrained by how many genetic solutions exist for a particular function. When fewer solutions exist, as in genetic sulfation of small molecules, evolution may be more constrained to use the same genes time and again.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Speciation patterns of Aedes mosquitoes in the Scutellaris Group: a mitochondrial perspective.

The Scutellaris Group of Aedes comprises 47 mosquito species, including Aedes albopictus. While Ae. albopictus is widely distributed, the other species are mostly found in the Asia-Pacific region. Evolutionary history researches of Aedes species within the Scutellaris Group have mainly focused on Ae. albopictus, a species that raises significant public health concerns, neglecting the other species. In this study, we aimed to assess genetic diversity and estimate speciation times of several species within the Scutellaris Group. Mosquitoes were therefore collected from various Asia-Pacific countries. Their mitochondrial cytochrome c oxidase subunit 1 (cox1) and subunit 3 (cox3) sequences were analyzed alongside those of other Scutellaris Group species available in the GenBank database. To estimate the divergence time, we analyzed 1849 cox1 gene sequences from 21 species, using three species (Aedes aegypti, Aedes notoscriptus and Aedes vigilax) as outgroups. We found that most of the speciation dates occurred during the Paleogene and the Neogene periods. A separation between the Scutellaris Subgroup and the Albopictus Subgroup occurred approximately 64-61 million years ago (MYA). We also identified a split between species found in Asia/Micronesia and those collected in Melanesia/Polynesia approximately 36-35 MYA. Our findings suggest that the speciation of Aedes species within the Scutellaris Group may be driven by diversity in mammalian hosts, climate and environmental changes, and geological dynamics rather than human migration.

A novel approach to exploring the dark genome and its application to mapping of the vertebrate virus fossil record.

BACKGROUND: Genomic regions that remain poorly understood, often referred to as the dark genome, contain a variety of functionally relevant and biologically informative features. These include endogenous viral elements (EVEs)-virus-derived sequences that can dramatically impact host biology and serve as a virus fossil record. In this study, we introduce a database-integrated genome screening (DIGS) approach to investigate the dark genome in silico, focusing on EVEs found within vertebrate genomes. RESULTS: Using DIGS on 874 vertebrate genomes, we uncover approximately 1.1 million EVE sequences, with over 99% originating from endogenous retroviruses or transposable elements that contain EVE DNA. We show that the remaining 6038 sequences represent over a thousand distinct horizontal gene transfer events across 10 virus families, including some that have not previously been reported as EVEs. We explore the genomic and phylogenetic characteristics of non-retroviral EVEs and determine their rates of acquisition during vertebrate evolution. Our study uncovers novel virus diversity, broadens knowledge of virus distribution among vertebrate hosts, and provides new insights into the ecology and evolution of vertebrate viruses. CONCLUSIONS: We comprehensively catalog and analyze EVEs within 874 vertebrate genomes, shedding light on the distribution, diversity, and long-term evolution of viruses and reveal their extensive impact on vertebrate genome evolution. Our results demonstrate the power of linking a relational database management system to a similarity search-based screening pipeline for in silico exploration of the dark genome.

Experimental Evaluation of a Direct Fitness Effect of the De Novo Evolved Mouse Gene Pldi.

De novo evolved genes emerge from random parts of noncoding sequences and have, therefore, no homologs from which a function could be inferred. While expression analysis and knockout experiments can provide insights into the function, they do not directly test whether the gene is beneficial for its carrier. Here, we have used a seminatural environment experiment to test the fitness of the previously identified de novo evolved mouse gene Pldi, which has been implicated to have a role in sperm differentiation. We used a knockout mouse strain for this gene and competed it against its parental wildtype strain for several generations of free reproduction. We found that the knockout (ko) allele frequency decreased consistently across three replicates of the experiment. Using an approximate Bayesian computation framework that simulated the data under a demographic scenario mimicking the experiment's demography, we could estimate a selection coefficient ranging between 0.21 and 0.61 for the wildtype allele compared to the ko allele in males, under various models. This implies a relatively strong selective advantage, which would fix the new gene in less than hundred generations after its emergence.

Mitochondrial genome fragmentation is correlated with increased rates of molecular evolution.

While mitochondrial genome content and organization is quite diverse across all Eukaryotes, most bilaterian animal mitochondrial genomes (mitogenomes) exhibit highly conserved gene content and organisation, with genes typically encoded on a single circular chromosome. However, many species of parasitic lice (Insecta: Phthiraptera) are among the notable exceptions, having mitogenomes fragmented into multiple circular chromosomes. To better understand the process of mitogenome fragmentation, we conducted a large-scale genomic study of a major group of lice, Amblycera, with extensive taxon sampling. Analyses of the evolution of mitogenome structure across a phylogenomic tree of 90 samples from 53 genera revealed evidence for multiple independent origins of mitogenome fragmentation, some inferred to have occurred less than five million years ago. We leveraged these many independent origins of fragmentation to compare the rates of DNA substitution and gene rearrangement, specifically contrasting branches with fragmented and non-fragmented mitogenomes. We found that lineages with fragmented mitochondrial genomes had significantly higher rates of mitochondrial sequence evolution. In addition, lineages with fragmented mitochondrial genomes were more likely to have mitogenome gene rearrangements than those with single-chromosome mitochondrial genomes. By combining phylogenomics and mitochondrial genomics we provide a detailed portrait of mitogenome evolution across this group of insects with a remarkably unstable mitogenome structure, identifying processes of molecular evolution that are correlated with mitogenome fragmentation.

A stochastic approach for co-evolution process of virus and human immune system.

Infectious diseases have long been a shaping force in human history, necessitating a comprehensive understanding of their dynamics. This study introduces a co-evolution model that integrates both epidemiological and evolutionary dynamics. Utilizing a system of differential equations, the model represents the interactions among susceptible, infected, and recovered populations for both ancestral and evolved viral strains. Methodologically rigorous, the model's existence and uniqueness have been verified, and it accommodates both deterministic and stochastic cases. A myriad of graphical techniques have been employed to elucidate the model's dynamics. Beyond its theoretical contributions, this model serves as a critical instrument for public health strategy, particularly predicting future outbreaks in scenarios where viral mutations compromise existing interventions.