Consensus Sequence Research Articles

Venomics and Peptidomics of Palearctic Vipers: A Clade-Wide Analysis of Seven Taxa of the Genera Vipera, Montivipera, Macrovipera, and Daboia across Türkiye.

Snake venom variations are a crucial factor to understand the consequences of snakebite envenoming worldwide, and therefore it is important to know about toxin composition alterations between taxa. Palearctic vipers of the genera Vipera, Montivipera, Macrovipera, and Daboia have high medical impacts across the Old World. One hotspot for their occurrence and diversity is Türkiye, located on the border between continents, but many of their venoms remain still understudied. Here, we present the venom compositions of seven Turkish viper taxa. By complementary mass spectrometry-based bottom-up and top-down workflows, the venom profiles were investigated on proteomics and peptidomics level. This study includes the first venom descriptions of Vipera berus barani, Vipera darevskii, Montivipera bulgardaghica albizona, and Montivipera xanthina, as well as the first snake venomics profiles of Turkish Macrovipera lebetinus obtusa, and Daboia palaestinae, including an in-depth reanalysis of M. bulgardaghica bulgardaghica venom. Additionally, we identified the modular consensus sequence pEXW(PZ)1-2P(EI)/(KV)PPLE for bradykinin-potentiating peptides in viper venoms. For better insights into variations and potential impacts of medical significance, the venoms were compared against other Palearctic viper proteomes, including the first genus-wide Montivipera venom comparison. This will help the risk assessment of snakebite envenoming by these vipers and aid in predicting the venoms' pathophysiology and clinical treatments.

Type I Interferon Activates PD-1 Expression through Activation of the STAT1-IRF2 Pathway in Myeloid Cells.

PD-1 (Programmed cell death protein 1) regulates the metabolic reprogramming of myeloid-derived suppressor cells and myeloid cell differentiation, as well as the type I interferon (IFN-I) signaling pathway in myeloid cells in the tumor microenvironment. PD-1, therefore, is a key inhibitory receptor in myeloid cells. However, the regulation of PD-1 expression in myeloid cells is unknown. We report that the expression level of PDCD1, the gene that encodes the PD-1 protein, is positively correlated with the levels of IFNB1 and IFNAR1 in myeloid cells in human colorectal cancer. Treatment of mouse myeloid cell lines with recombinant IFNβ protein elevated PD-1 expression in myeloid cells in vitro. Knocking out IFNAR1, the gene that encodes the IFN-I-specific receptor, diminished the inductive effect of IFNβ on PD-1 expression in myeloid cells in vitro. Treatment of tumor-bearing mice with a lipid nanoparticle-encapsulated IFNβ-encoding plasmid (IFNBCOL01) increased IFNβ expression, resulting in elevated PD-1 expression in tumor-infiltrating myeloid cells. At the molecular level, we determined that IFNβ activates STAT1 (signal transducer and activator of transcription 1) and IRFs (interferon regulatory factors) in myeloid cells. Analysis of the cd279 promoter identified IRF2-binding consensus sequence elements. ChIP (chromatin immunoprecipitation) analysis determined that the pSTAT1 directly binds to the irf2 promoter and that IRF2 directly binds to the cd279 promoter in myeloid cells in vitro and in vivo. In colon cancer patients, the expression levels of STAT1, IRF2 and PDCD1 are positively correlated in tumor-infiltrating myeloid cells. Our findings determine that IFNβ activates PD-1 expression at least in part by an autocrine mechanism via the stimulation of the pSTAT1-IRF2 axis in myeloid cells.

Sources of Human Immunodeficiency Virus Infections Among Men Who Have Sex With Men With a Migration Background: A Viral Phylogenetic Case Study in Amsterdam, The Netherlands.

Men and women with a migration background comprise an increasing proportion of incident human immunodeficiency virus (HIV) cases across Western Europe. To characterize sources of transmission in local transmission chains, we used partial HIV consensus sequences with linked demographic and clinical data from the opt-out AIDS Therapy Evaluation in the Netherlands (ATHENA) cohort of people with HIV in the Netherlands and identified phylogenetically and epidemiologically possible HIV transmission pairs in Amsterdam. We interpreted these in the context of estimated infection dates, and quantified population-level sources of transmission to foreign-born and Dutch-born Amsterdam men who have sex with men (MSM) within Amsterdam transmission chains. We estimate that Dutch-born MSM were the predominant sources of infections among all Amsterdam MSM who acquired their infection locally in 2010-2021, and among almost all foreign-born Amsterdam MSM subpopulations. Stratifying by 2-year intervals indicated time trends in transmission dynamics, with a majority of infections originating from foreign-born MSM since 2016, although uncertainty ranges remained wide. Native-born MSM have predominantly driven HIV transmissions in Amsterdam in 2010-2021. However, in the context of rapidly declining incidence in Amsterdam, the contribution from foreign-born MSM living in Amsterdam is increasing, with some evidence that most local transmissions have been from foreign-born Amsterdam MSM since 2016.

Multiple types of nuclear localization signals in Entamoeba histolytica

Entamoeba histolytica is a protozoan parasite that belongs to the Amoebozoa supergroup whose study related to the nucleocytoplasmic transport of proteins through the nucleus is poorly studied. In this work, we have performed in silico predictions of the potential nuclear localization signals (NLS) corresponding to the proteome of 8201 proteins from Entamoeba histolytica annotated in the AmoebaDB database. We have found the presence of monopartite nuclear localization signals (MNLSs), bipartite nuclear localization signals (BNLSs), and non-canonical monopartite NLSs with lengths exceeding 20 amino acid residues. Additionally, we detected a new type of NLS consisting of multiple juxtaposed bipartite NLSs (JNLSs) that have not been described in any eukaryotic organism. Also, we have generated consensus sequences for the nuclear import of proteins with the NLSs obtained. Docking experiments between EhImportin α and an MNLS, BNLS, and JNLS outlined the interacting residues between the Importin and cargo proteins, emphasizing their putative roles in nuclear import. By transfecting HA-tagged protein constructs, we assessed the nuclear localization of MNLS (U1A and U2AF1), JMNLS (U2AF2), and non-canonical NLS (N-terminus of Pol ll) in vivo. Our data provide the basis for understanding the nuclear transport process in E. histolytica.

Latent Dirichlet allocation mixture models for nucleotide sequence analysis.

Strings of nucleotides carrying biological information are typically described as sequence motifs represented by weight matrices or consensus sequences. However, many signals in DNA or RNA are recognized by multiple factors in temporal sequence, consist of distinct alternative motifs, or are best described by base composition. Here we apply the latent Dirichlet allocation (LDA) mixture model to nucleotide sequences. Using positions in an alignment of human or Drosophila splice sites as samples, we show that LDA readily identifies motifs, including such elusive cases as the intron branch site. Using whole sequences with positional k-mers as features, LDA can identify sequence subtypes enriched in long vs. short introns. LDA with bulk k-mers can reliably distinguish reading frame and species of origin in coding sequences from humans and Drosophila. We find that LDA is a useful model for describing heterogeneous signals, for assigning individual sequences to subtypes, and for identifying and characterizing sequences that do not fit recognized subtypes. Because LDA topic models are interpretable, they also aid the discovery of new motifs, even those present in a small fraction of samples. In summary, LDA can identify and characterize signals in nucleotide sequences, including candidate regulatory factors involved in biological processes.

Towards a better future for DNA barcoding: Evaluating monophyly- and distance-based species identification using COI gene fragments of Dacini fruit flies.

The utility of a universal DNA 'barcode' fragment (658 base pairs of the Cytochrome C Oxidase I [COI] gene) has been established as a useful tool for species identification, and widely criticized as one for understanding the evolutionary history of a group. Large amounts of COI sequence data have been produced that hold promise for rapid species identification, for example, for biosecurity. The fruit fly tribe Dacini holds about a thousand species, of which 80 are pests of economic concern. We generated a COI reference library for 265 species of Dacini containing 5601 sequences that span most of the COI gene using circular consensus sequencing. We compared distance metrics versus monophyly assessments for species identification and although we found a 'soft' barcode gap around 2% pairwise distance, the exceptions to this rule dictate that a monophyly assessment is the only reliable method for species identification. We found that all fragments regularly used for Dacini fruit fly identification >450 base pairs long provide similar resolution. 11.3% of the species in our dataset were non-monophyletic in a COI tree, which is mostly due to species complexes. We conclude with recommendations for the future generation and use of COI libraries. We revise the generic assignment of Dacus transversus stat. rev. Hardy 1982, and Dacus perpusillus stat. rev. Drew 1971 and we establish Dacus maculipterus White 1998 syn. nov. as a junior synonym of Dacus satanas Liang etal. 1993.

Long-read sequencing identifies an SVA_D retrotransposon insertion deep within the intron of ATP7A as a novel cause of occipital horn syndrome

BackgroundSINE-VNTR-Alu (SVA) retrotransposons move from one genomic location to another in a ‘copy-and-paste’ manner. They continue to move actively and cause monogenic diseases through various mechanisms. Currently, disease-causing SVA retrotransposons...

Linguistic and cultural specifics of the sports media narrative (based on the example of sports commentaries in Canada and Russia)

This study is devoted to the comparative study of the specifics of the sports media narrative. The aim of the study is to compare the sports media narrative in the linguistic and cultural paradigm (based on the material of sports commentaries from Canada and Russia). The object of the study is the sports media narrative in the Canadian and Russian linguistic cultures. The subject of the study is the linguistic and cultural aspect of the sports media narrative of Russia and Canada. The relevance of this article lies in the fact that events in the field of sports reflect the problems of globalization and forecasting the future, which determines the importance of both a comprehensive and aspect-based study of communication in this field and, in particular, the sports media narrative. In addition, comparative studies that identify the linguistic and cultural features of the modern media text, its expressive capabilities as a means of organizing narration, correspond to the methodology of modern language science. The main methods used in this work are comparative discourse analysis and the linguistic narrative method. In addition, quantitative data processing methods were used. The novelty of this work is due to the fact that for the first time a comparative analysis of the sports media narrative in English and Russian (based on the material of sports commentary) is carried out. Conclusions: oral texts produced within the framework of both discourses are characterized by internal dynamics, polycode (the commentator's speech is accompanied by a video sequence) and are arranged in accordance with the canonical narrative sequence. When describing the image of an event, commentators, speakers of different linguistic cultures, use an extra-egitic strategy (an omniscient observer-narrator), sometimes using an introdiegetic strategy (the desire to indirectly become a participant in the narrated events). Linguistic and cultural differences are most clearly seen at the level of speech types. Russian commentators are characterized by the use of informative narration, while Canadian commentators use descriptive narration. The second most significant type of speech in the Russian media discourse is reasoning. For Canadian commentators, the description is more typical. The Russian sports media narrative follows a dramatic (plot) way of showing events, while the Canadian one follows a chronotopic one, i.e. it is determined by reference to time and place. The revealed differences may indicate that when narrating an event, commentators create a "second reality" based on their axiological attitudes due to cultural differences between the two countries.

Dimerization is required for the glycosylation of S1-S2 linker of sea urchin voltage-gated proton channel Hv1

Multimerization of ion channels is essential for establishing the ion-selective pathway and tuning the gating regulated by membrane potential, second messengers, and temperature. Voltage-gated proton channel, Hv1, consists of voltage-sensor domain and coiled-coil domain. Hv1 forms dimer, whereas voltage-dependent channel activity is self-contained in monomer unlike many ion channels, which assemble to form ion-conductive pathways among multiple subunits. Dimerization of Hv1 is necessary for cooperative gating, but other roles of dimerization in physiological aspects are still largely unclear. In this study, we show that dimerization of Hv1 takes place in ER. Sea urchin Hv1 (Strongylocentrotus purpuratus Hv1: SpHv1) was glycosylated in the consensus sequence for N-linked glycosylation within the S1-S2 extracellular loop. However, glycosylation was not observed in the monomeric SpHv1 that lacks the coiled-coil domain. A version of mHv1 in which the S1-S2 loop was replaced by that of SpHv1 showed glycosylation and its monomeric form was not glycosylated. Tandem dimer of monomeric SpHv1 underwent glycosylation, suggesting that dimerization of Hv1 is required for glycosylation. Moreover, when monomeric Hv1 has a dilysine motif in the C-terminal end, which is known to act as a retrieval signal from Golgi to ER, prolonging the time of residency in ER, it was glycosylated. Overall, our results suggest that monomeric SpHv1 does not stay long in ER, thereby escaping glycosylation, while the dimerization causes the proteins to stay longer in ER. Thus, the findings highlight the novel significance of dimerization of Hv1: regulation of biogenesis and maturation of the proteins in intracellular compartments.

Interepidemic xenosurveillance of Japanese encephalitis virus and Zika virus in Culex mosquitoes from Ubon Ratchathani province, Thailand.

Some Culex mosquitoes are competent vectors for Japanese encephalitis virus (JEV) and Zika virus (ZIKV), which cause public health problems worldwide, especially in South-east Asia. Xenosurveillance of Culex mosquitoes remains limited compared with other common mosquito-borne diseases. This study aimed to identify JEV and ZIKV in field-caught Culex mosquitoes collected from Ubon Ratchathani province. We investigated the presence of JEV and ZIKV in Culex mosquitoes from two districts in Ubon Ratchathani province, Thailand, and examined their role in viral interepidemic circulation. Female Culex mosquitoes (5,587) were collected using a mechanical aspirator from indoors and outdoors. The consensus sequences of the E and NS1 genes of JEV and the E gene of ZIKV were identified using real-time reverse transcription-polymerase chain reaction. From 335 sample pools that contain a total of 5587 adult female Culex mosquitoes collected from Don Yung, Mueang district (4,406) and Phon Duan, Det Udom district (1,181), none of the collected mosquitoes tested positive for either JEV or ZIKV. This study did not find JEV and ZIKV in Culex mosquitoes collected from the area of collection, which may be due to the low circulating amount of the virus in the vectors in the area, making it undetectable, or it may be because Culex mosquitoes are not suitable vector for the virus being tested. However, further xenosurveillance study of JEV and ZIKV in mosquito is suggested to prepare for the next outbreak.

Taxonomy of introns and the evolution of minor introns.

Classification of introns, which is crucial to understanding their evolution and splicing, has historically been binary and has resulted in the naming of major and minor introns that are spliced by their namesake spliceosome. However, a broad range of intron consensus sequences exist, leading us to here reclassify introns as minor, minor-like, hybrid, major-like, majorand non-canonical introns in 263 species across six eukaryotic supergroups. Through intron orthology analysis, we discovered that minor-like introns are a transitory node for intron conversion across evolution. Despite close resemblance of their consensus sequences to minor introns, these introns possess an AG dinucleotide at the -1 and -2 position of the 5' splice site, a salient feature of major introns. Through combined analysis of CoLa-seq, CLIP-seq for major and minor spliceosome components, and RNAseq from samples in which the minor spliceosome is inhibited we found that minor-like introns are also an intermediate class from a splicing mechanism perspective. Importantly, this analysis has provided insight into the sequence elements that have evolved to make minor-like introns amenable to recognition by both minor and major spliceosome components. We hope that this revised intron classification provides a new framework to study intron evolution and splicing.

Predictive and Experimental Motif Interaction Analysis Identifies Functions of the WNK-OSR1/SPAK Pathway.

The WNK-OSR1/SPAK protein kinase signaling pathway regulates ion homeostasis and cell volume, but its other functions are poorly understood. To uncover undefined signaling functions of the pathway we analyzed the binding specificity of the conserved C-terminal (CCT) domains of OSR1 and SPAK to find all possible interaction motifs in human proteins. These kinases bind the core consensus sequences R-F-x-V/I and R-x-F-x-V/I. Motifs were ranked based on sequence, conservation, cellular localization, and solvent accessibility. Out of nearly 3,700 motifs identified, 90% of previously published motifs were within the top 2% of those predicted. Selected candidates (TSC22D1, CAVIN1, ATG9A, NOS3, ARHGEF5) were tested. Upstream kinases WNKs 1-4 and their close relatives, the pseudokinases NRBP1/2, contain CCT-like domains as well. We identified additional distinct motif variants lacking the conserved arginine previously thought to be required, and found that the NRBP1 CCT-like domain binds TSC22D1 via the same motif as OSR1 and SPAK. Our results further highlight the rich and diverse functionality of CCT and CCT-like domains in connecting WNK signaling to cellular processes.

Development of a transcription factor decoy-nanocarrier system as a successful inhibitor of Enterococcus faecalis virulence in vitro and in vivo

Enterococcus faecalis is a troublesome nosocomial pathogen that acquired resistance to most available antimicrobial agents. Antivirulence agents represent an unconventional treatment approach. Here, transcription factor decoy (TFD)-loaded cationic liposomes (TLL) were developed as an inhibitor of the Fsr quorum-sensing system and its associated virulence traits, in E. faecalis. The consensus sequence of the FsrA binding site was found conserved among 651 E. faecalis annotated genomes. The TFD was synthesized as an 82 bp DNA duplex, containing the conserved binding sequence, and loaded onto cationic liposomes. The optimum loading capacity, mean particle size, and zeta potential of the TLL were characterized. The developed TLL lacked any effect on E. faecalis growth and significantly inhibited the in vitro production of the proteolytic enzymes controlled by the Fsr system; gelatinase and serine protease, in a concentration-dependent manner. This inhibition was accompanied by a significant reduction in the transcription levels of FsrA-regulated genes (fsrB, gelE, and sprE). The developed TLL were safe as evidenced by the nonhemolytic effect on human RBCs and the negligible cytotoxicity on human skin fibroblast cells. Moreover, in the larvae infection model, TLL displayed a significant abolish in the mortality rates of Galleria mellonella larvae infected with E. faecalis. In conclusion, the developed TLL offer a new safe strategy for combating E. faecalis infection through the inhibition of quorum-sensing-mediated virulence; providing a platform for the development of similar agents to combat many other pathogens.

The nuclear localization signal of CPSF6 governs post-nuclear import steps of HIV-1 infection.

The early stages of HIV-1 infection include the trafficking of the viral core into the nucleus of infected cells. However, much remains to be understood about how HIV-1 accomplishes nuclear import and the consequences of the import pathways utilized on nuclear events. The host factor cleavage and polyadenylation specificity factor 6 (CPSF6) assists HIV-1 nuclear localization and post-entry integration targeting. Here, we used a CPSF6 truncation mutant lacking a functional nuclear localization signal (NLS), CPSF6-358, and appended heterologous NLSs to rescue nuclear localization. We show that some, but not all, NLSs drive CPSF6-358 into the nucleus. Interestingly, we found that some nuclear localized CPSF6-NLS chimeras supported inefficient HIV-1 infection. We found that HIV-1 still enters the nucleus in these cell lines but fails to traffic to speckle-associated domains (SPADs). Additionally, we show that HIV-1 fails to efficiently integrate in these cell lines. Collectively, our results demonstrate that the NLS of CPSF6 facilitates steps of HIV-1 infection subsequent to nuclear import and additionally identify the ability of canonical NLS sequences to influence cargo localization in the nucleus following nuclear import.

Replication initiation sites and zones in the mammalian genome: Where are they located and how are they defined?

Eukaryotic DNA replication is a tightly controlled process that occurs in two main steps, i.e., licensing and firing, which take place in the G1 and S phases of the cell cycle, respectively. In Saccharomyces cerevisiae, the budding yeast, replication origins contain consensus sequences that are recognized and bound by the licensing factor Orc1–6, which then recruits the replicative Mcm2–7 helicase. By contrast, mammalian initiation sites lack such consensus sequences, and the mammalian ORC does not exhibit sequence specificity. Studies performed over the past decades have identified replication initiation sites in the mammalian genome using sequencing-based assays, raising the question of whether replication initiation occurs at confined sites or in broad zones across the genome. Although recent reports have shown that the licensed MCMs in mammalian cells are broadly distributed, suggesting that ORC-dependent licensing may not determine the initiation sites/zones, they are predominantly located upstream of actively transcribed genes. This review compares the mechanism of replication initiation in yeast and mammalian cells, summarizes the sequencing-based technologies used for the identification of initiation sites/zones, and proposes a possible mechanism of initiation-site/zone selection in mammalian cells. Future directions and challenges in this field are also discussed.

Pumilio-1 mediated translational control of claudin-5 at the blood-brain barrier

Claudin-5 is one of the most essential tight junction proteins at the blood-brain barrier. A single nucleotide polymorphism rs10314 is located in the 3’-untranslated region of claudin-5 and has been shown to be a risk factor for schizophrenia. Here, we show that the pumilio RNA-binding protein, pumilio-1, is responsible for rs10314-mediated claudin-5 regulation. The RNA sequence surrounding rs10314 is highly homologous to the canonical pumilio-binding sequence and claudin-5 mRNA with rs10314 produces 25% less protein due to its inability to bind to pumilio-1. Pumilio-1 formed cytosolic granules under stress conditions and claudin-5 mRNA appeared to preferentially accumulate in these granules. Added to this, we observed granular pumilio-1 in endothelial cells in human brain tissues from patients with psychiatric disorders or epilepsy with increased/accumulated claudin-5 mRNA levels, suggesting translational claudin-5 suppression may occur in a brain-region specific manner. These findings identify a key regulator of claudin-5 translational processing and how its dysregulation may be associated with neurological and neuropsychiatric disorders.Graphical

Short Notes: First Report of Cucurbite Aphid-Borne Yellows Virus in Iraq

Squash (Cucurbita pepo) fields in Mosul suffered from virus-like symptoms, such as vein clearing, severe mosaic, mottling, chlorosis, necrosis, stunting and fruit discoloration. The symptoms were associated with aphids presence in the infected fields. A metatranscriptomic analysis has been applied to detect the causal agent. The result showed that 559,024 mapped reads were related to Cucurbit aphid-borne yellows virus with full coverage to produce 5,662 bp consensus sequence. The virus sequence was deposited in National Center for Biotechnology Information (NCBI) under accession number OQ865059, and the isolate named Iraq. For the first time in Iraq, the infection of the virus has been reported in squash, highlighting the impact of such virus infection that spread epidemically across the squash and other cucurbits fields.

Pioneer factors for DNA replication initiation in metazoans.

Precise DNA replication is fundamental for genetic inheritance. In eukaryotes, replication initiates at multiple origins that are first "licensed" and subsequently "fired" to activate DNA synthesis. Despite the success in identifying origins with specific DNA motifs in Saccharomyces cerevisiae, no consensus sequence or sequences with a predictive value of replication origins have been recognized in metazoan genomes. Rather, epigenetic rules and chromatin structures are believed to play important roles in governing the selection and activation of replication origins. We propose that replication initiation is facilitated by a group of sequence-specific "replication pioneer factors," which function to increase chromatin accessibility and foster a chromatin environment that is conducive to the loading of the prereplication complex. Dysregulation of the function of these factors may lead to gene duplication, genomic instability, and ultimately the occurrence of pathological conditions such as cancer.

PCNA molecular recognition of different PIP motifs: Role of Tyr211 phosphorylation

The coordination of enzymes and regulatory proteins for eukaryotic DNA replication and repair is largely achieved by Proliferating Cell Nuclear Antigen (PCNA), a toroidal homotrimeric protein that embraces the DNA duplex. Many proteins bind PCNA through a conserved sequence known as the PCNA interacting protein motif (PIP). PCNA is further regulated by different post-translational modifications. Phosphorylation at residue Y211 facilitates unlocking stalled replication forks to bypass DNA damage repair processes but increasing nucleotide misincorporation. We explore here how phosphorylation at Y211 affects PCNA recognition of the canonical PIP sequences of the regulatory proteins p21 and p15, which bind with nM and μM affinity, respectively. For that purpose, we have prepared PCNA with p-carboxymethyl-L-phenylalanine (pCMF, a mimetic of phosphorylated tyrosine) at position 211. We have also characterized PCNA binding to the non-canonical PIP sequence of the catalytic subunit of DNA polymerase δ (p125), and to the canonical PIP sequence of the enzyme ubiquitin specific peptidase 29 (USP29) which deubiquitinates PCNA. Our results show that Tyr211 phosphorylation has little effect on the molecular recognition of p21 and p15, and that the PIP sequences of p125 and USP29 bind to the same site on PCNA as other PIP sequences, but with very low affinity.

The effect of molnupiravir and nirmatrelvir on SARS-CoV-2 genome diversity in severe models of COVID-19.

Immunocompromised individuals are susceptible to severe COVID-19 and potentially contribute to the emergence of variants with altered pathogenicity due to persistent infection. This study investigated the impact of immunosuppression on SARS-CoV-2 infection in k18-hACE2 mice and the effectiveness of antiviral treatments in this context during the first 7 days of infection. Mice were immunosuppressed using cyclophosphamide and infected with a B lineage of SARS-CoV-2. Molnupiravir and nirmatrelvir, alone and in combination, were administered and viral load and viral sequence diversity was assessed. Treatment of infected but immune compromised mice with both compounds either singly or in combination resulted in decreased viral loads and pathological changes compared to untreated animals. Treatment also abrogated infection of neuronal tissue. However, no consistent changes in the viral consensus sequence were observed, except for the emergence of the S:H655Y mutation. Molnupiravir, but not nirmatrelvir or immunosuppression alone, increased the transition/transversion (Ts/Tv) ratio, representative of A>G and C>U mutations and this increase was not altered by the co-administration of nirmatrelvir with molnupiravir.Notably, immunosuppression itself did not appear to promote the emergence of mutational characteristic of variants of concern (VOCs). Further investigations are warranted to fully understand the role of immunocompromised individuals in VOC development, especially by taking persistence into consideration, and to inform optimised public health strategies. It is more likely that immunodeficiency promotes viral persistence but does not necessarily lead to substantial consensus-level changes in the absence of antiviral selection pressure. Consistent with mechanisms of action, molnupiravir showed a stronger mutagenic effect than nirmatrelvir in this model.