DengueSeq: a pan-serotype whole genome amplicon sequencing protocol for dengue virus

BACKGROUND

The increasing burden of dengue virus on public health due to more explosive and frequent outbreaks highlights the need for improved surveillance and control. Genomic surveillance of dengue virus not only provides important insights into the emergence and spread of genetically diverse serotypes and genotypes, but it is also critical to monitor the effectiveness of newly implemented control strategies. Here, we present DengueSeq, an amplicon sequencing protocol, which enables whole-genome sequencing of all four dengue virus serotypes.

RESULTS

We developed primer schemes for the four dengue virus serotypes, which can be combined into a pan-serotype approach. We validated both approaches using genetically diverse virus stocks and clinical specimens that contained a range of virus copies. High genome coverage (>95%) was achieved for all genotypes, except DENV2 (genotype VI) and DENV 4 (genotype IV) sylvatics, with similar performance of the serotype-specific and pan-serotype approaches. The limit of detection to reach 70% coverage was 10-100 RNA copies/μL for all four serotypes, which is similar to other commonly used primer schemes. DengueSeq facilitates the sequencing of samples without known serotypes, allows the detection of multiple serotypes in the same sample, and can be used with a variety of library prep kits and sequencing instruments.

CONCLUSIONS

DengueSeq was systematically evaluated with virus stocks and clinical specimens spanning the genetic diversity within each of the four dengue virus serotypes. The primer schemes can be plugged into existing amplicon sequencing workflows to facilitate the global need for expanded dengue virus genomic surveillance.

Circulation, viral diversity and genomic rearrangement in mpox virus in the Netherlands during the 2022 outbreak and beyond

Mpox is an emerging zoonotic disease which has now spread to over 113 countries as of August 2023, with over 89,500 confirmed human cases. The Netherlands had one of the highest incidence rates in Europe during the peak of the outbreak. In this study, we generated 158 near-complete mpox virus (MPXV) genomes (12.4% of nationwide cases) that were collected throughout the Netherlands from the start of the outbreak in May 2022 to August 2023 to track viral evolution and investigate outbreak dynamics. We detected 14 different viral lineages, suggesting multiple introductions followed by rapid initial spread within the country. The estimated evolutionary rate was relatively high compared to previously described in orthopoxvirus literature, with an estimated 11.58 mutations per year. Genomic rearrangement events occurred at a rate of 0.63% and featured a large deletion event. In addition, based on phylogenetics, we identified multiple potential transmission clusters which could be supported by direct source- and contact tracing data. This led to the identification of at least two main transmission locations at the beginning of the outbreak. We conclude that whole genome sequencing of MPXV is essential to enhance our understanding of outbreak dynamics and evolution of a relatively understudied and emerging zoonotic pathogen.

A comparison of five Illumina, Ion Torrent, and nanopore sequencing technology-based approaches for whole genome sequencing of SARS-CoV-2

Rapid identification of the rise and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern remains critical for monitoring of the efficacy of diagnostics, therapeutics, vaccines, and control strategies. A wide range of SARS-CoV-2 next-generation sequencing (NGS) methods have been developed over the last years, but cross-sequence technology benchmarking studies have been scarce. In the current study, 26 clinical samples were sequenced using five protocols: AmpliSeq SARS-CoV-2 (Illumina), EasySeq RC-PCR SARS-CoV-2 (Illumina/NimaGen), Ion AmpliSeq SARS-CoV-2 (Thermo Fisher), custom primer sets (Oxford Nanopore Technologies (ONT)), and capture probe-based viral metagenomics (Roche/Illumina). Studied parameters included genome coverage, depth of coverage, amplicon distribution, and variant calling. The median SARS-CoV-2 genome coverage of samples with cycle threshold (Ct) values of 30 and lower ranged from 81.6 to 99.8% for, respectively, the ONT protocol and Illumina AmpliSeq protocol. Correlation of coverage with PCR Ct values varied per protocol. Amplicon distribution signatures differed across the methods, with peak differences of up to 4 log₁₀ at disbalanced positions in samples with high viral loads (Ct values ≤ 23). Phylogenetic analyses of consensus sequences showed clustering independent of the workflow used. The proportion of SARS-CoV-2 reads in relation to background sequences, as a (cost-)efficiency metric, was the highest for the EasySeq protocol. The hands-on time was the lowest when using EasySeq and ONT protocols, with the latter additionally having the shortest sequence runtime. In conclusion, the studied protocols differed on a variety of the studied metrics. This study provides data that assist laboratories when selecting protocols for their specific setting.

Reptile-associated Salmonella urinary tract infection: a case report

We present an 18-year-old woman with a urinary tract infection caused by Salmonella Oranienburg. S. Oranienburg was isolated from her pet snake and confirmed as the source of infection using whole genome sequencing. Our case demonstrates the risk of acquiring reptile-associated salmonellosis and stretches the need for awareness to prevent these infections.

Optimizing Antiviral Dosing for HSV and CMV Treatment in Immunocompromised Patients

Herpes simplex virus (HSV) and cytomegalovirus (CMV) are DNA viruses that are common among humans. Severely immunocompromised patients are at increased risk of developing HSV or CMV disease due to a weakened immune system. Antiviral therapy can be challenging because these drugs have a narrow therapeutic window and show significant pharmacokinetic variability. Above that, immunocompromised patients have various comorbidities like impaired renal function and are exposed to polypharmacy. This scoping review discusses the current pharmacokinetic (PK) and pharmacodynamic (PD) knowledge of antiviral drugs for HSV and CMV treatment in immunocompromised patients. HSV and CMV treatment guidelines are discussed, and multiple treatment interventions are proposed: early detection of drug resistance; optimization of dose to target concentration by therapeutic drug monitoring (TDM) of nucleoside analogs; the introduction of new antiviral drugs; alternation between compounds with different toxicity profiles; and combinations of synergistic antiviral drugs. This research will also serve as guidance for future research, which should focus on prospective evaluation of the benefit of each of these interventions in randomized controlled trials.

SARS-CoV-2 outbreaks in secondary school settings in the Netherlands during fall 2020; silent circulation

BACKGROUND

In fall 2020 when schools in the Netherlands operated under a limited set of COVID-19 measures, we conducted outbreaks studies in four secondary schools to gain insight in the level of school transmission and the role of SARS-CoV-2 transmission via air and surfaces.

METHODS

Outbreak studies were performed between 11 November and 15 December 2020 when the wild-type variant of SARS-CoV-2 was dominant. Clusters of SARS-CoV-2 infections within schools were identified through a prospective school surveillance study. All school contacts of cluster cases, irrespective of symptoms, were invited for PCR testing twice within 48 h and 4-7 days later. Combined NTS and saliva samples were collected at each time point along with data on recent exposure and symptoms. Surface and active air samples were collected in the school environment. All samples were PCR-tested and sequenced when possible.

RESULTS

Out of 263 sampled school contacts, 24 tested SARS-CoV-2 positive (secondary attack rate 9.1%), of which 62% remained asymptomatic and 42% had a weakly positive test result. Phylogenetic analysis on 12 subjects from 2 schools indicated a cluster of 8 and 2 secondary cases, respectively, but also other distinct strains within outbreaks. Of 51 collected air and 53 surface samples, none were SARS-CoV-2 positive.

CONCLUSION

Our study confirmed within school SARS-CoV-2 transmission and substantial silent circulation, but also multiple introductions in some cases. Absence of air or surface contamination suggests environmental contamination is not widespread during school outbreaks.

Development of Resistance-Associated Mutations After Sotrovimab Administration in High-risk Individuals Infected With the SARS-CoV-2 Omicron Variant

This study investigates whether resistance-associated mutations develop after treatment with sotrovimab in high-risk patients infected with the SARS-CoV-2 Omicron variant.

[Monkeypox, a new pandemic?]

Monkeypox (MPX) is a disease caused by the monkeypox virus. It is a viral zoonotic disease, endemic in Central and West Africa. Human-to-human spread also occurs and is a feature of the current global outbreak. As far as we know, exponential transmission during this outbreak is not related to changed viral characteristics but due to multiple high-risk contacts in a subset of people that have contracted the virus, so far almost exclusively affecting men who have sex with men (MSM). Appropriate public health measures and increased alertness of all health care providers is needed to increase case-finding and decrease transmission. There is a real chance of MPX to become endemic in large parts of the world.

Hyperimmune Globulin for Severely Immunocompromised Patients Hospitalized With Coronavirus Disease 2019: A Randomized, Controlled Trial

BACKGROUND

The aim of this randomized, controlled trial is to determine whether antisevere acute respiratory syndrome coronavirus 2 hyperimmune globulin (COVIG) protects against severe coronavirus disease 2019 (COVID-19) in severely immunocompromised, hospitalized, COVID-19 patients.

METHODS

Patients were randomly assigned to receive COVIG or intravenous immunoglobulin (IVIG) without SARS-CoV-2 antibodies.

RESULTS

Severe COVID-19 was observed in 2 of 10 (20%) patients treated with COVIG compared to 7 of 8 (88%) in the IVIG control group (P = .015, Fisher's exact test).

CONCLUSIONS

Antisevere acute respiratory syndrome coronavirus 2 hyperimmune globulin may be a valuable treatment in severely immunocompromised, hospitalized, COVID-19 patients and should be considered when no monoclonal antibody therapies are available.

Clinically relevant DNA viruses in pregnancy

Infections by DNA viruses during pregnancy are associated with increased health risks to both mother and fetus. Although not all DNA viruses are related to an increased risk of complications during pregnancy, several can directly infect the fetus and/or cause placental dysfunction. During NIPT analysis, the presence of viral DNA can be detected, theoretically allowing screening early in pregnancy. Although treatment options are currently limited, this might rapidly change in the near future. It is therefore important to be aware of the potential impact of these viruses on feto-maternal health. In this manuscript we provide a brief introduction into the most commonly detected DNA viruses in human cell-free DNA sequencing experiments and their pathogenic potential during pregnancy. This article is protected by copyright. All rights reserved.

Rapid reinfection with SARS-CoV-2 variant-of-concern Alpha detected in a nurse during an outbreak at a non-covid inpatient ward: lessons learned

We describe the lessons learned during a SARS-CoV-2 variant-of-concern Alpha outbreak investigation at a normal care unit in a university hospital in Amsterdam in December 2020. The outbreak consisted of nine nurses and two roomed-in patient family members. (attack rate 18%). One nurse tested positive with a phylogenetically distinct variant, after a documented infection 83 days prior. Three key points were taken from this investigation. First, it was controlled by adherence to existing guidelines, despite increased transmissibility of the variant. Second, viral sequencing can inform transmission cluster inference, but the epidemiological context is essential to draw appropriate conclusions. Third, reinfections with Alpha variants can occur rapidly after primary infection.

Possible host-adaptation of SARS-CoV-2 due to improved ACE2 receptor binding in mink

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections on mink farms are increasingly observed in several countries, leading to the massive culling of animals on affected farms. Recent studies showed multiple (anthropo)zoonotic transmission events between humans and mink on these farms. Mink-derived SARS-CoV-2 sequences from The Netherlands and Denmark contain multiple substitutions in the S protein receptor binding domain (RBD). Molecular modeling showed that these substitutions increase the mean binding energy, suggestive of potential adaptation of the SARS-CoV-2 S protein to the mink angiotensin-converting enzyme 2 (ACE2) receptor. These substitutions could possibly also impact human ACE2 binding affinity as well as humoral immune responses directed to the RBD region of the SARS-CoV-2 S protein in humans. We wish to highlight these observations to raise awareness and urge for the continued surveillance of mink (and other animal)-related infections.

Genetic diversity and host adaptation of avian H5N1 influenza viruses during human infection

The continuing pandemic threat posed by avian influenza A/H5N1 viruses calls for improved insights into their evolution during human infection. We performed whole genome deep sequencing of respiratory specimens from 44 H5N1-infected individuals from Indonesia and found substantial within-host viral diversity. At nearly 30% of genome positions multiple amino acids were observed within or across samples, including positions implicated in aerosol transmission between ferrets. Amino acid variants detected our cohort were often found more frequently in available H5N1 sequences of human than avian isolates. We additionally identified previously unreported amino acid variants and multiple variants that increased in proportion over time in available sequential samples. Given the importance of the polymerase complex for host adaptation, we tested 121 amino acid variants found in the PB2, PB1 and PA subunits for their effects on polymerase activity in human cells. We identified multiple single amino acid variants in all three polymerase subunits that substantially increase polymerase activity including some with effects comparable to that of the widely recognized adaption and virulence marker PB2-E627 K. These results indicate highly dynamic evolutionary processes during human H5N1 virus infection and the potential existence of previously undocumented adaptive pathways.

Deep sequencing identifies hepatitis B virus core protein signatures in chronic hepatitis B patients

BACKGROUND

We aimed to identify HBc amino acid differences between subgroups of chronic hepatitis B (CHB) patients.

METHODS

Deep sequencing of HBc was performed in samples of 89 CHB patients (42 HBeAg positive, 47 HBeAg negative). Amino acid types were compared using Sequence Harmony to identify subgroup specific sites between HBeAg-positive and -negative patients, and between patients with combined response and non-response to peginterferon/adefovir combination therapy.

RESULTS

We identified 54 positions in HBc where the frequency of appearing amino acids was significantly different between HBeAg-positive and -negative patients. In HBeAg negative patients, 22 positions in HBc were identified which differed between patients with treatment response and those with non-response. The fraction non-consensus sequence on selected positions was significantly higher in HBeAg-negative patients, and was negatively correlated with HBV DNA and HBsAg levels.

CONCLUSIONS

Sequence Harmony identified a number of amino acid changes associated with HBeAg-status and response to peginterferon/adefovir combination therapy.

Viral minority variants in the core promoter and precore region identified by deep sequencing are associated with response to peginterferon and adefovir in HBeAg negative chronic hepatitis B patients

BACKGROUND AND AIM

Precore (PC) and basal core promoter (BCP) mutations are associated with responses to interferon-based treatment in HBeAg-positive chronic hepatitis B (CHB) patients. Here, we identify viral minority variants in these regions and assess association with response to peginterferon-alfa (Peg-IFN) and adefovir combination therapy.

PATIENTS AND METHODS

Ultra-deep pyrosequencing analysis of the BCP and PC region was performed for 89 CHB patients (42 HBeAg-positive; 47 HBeAg-negative), at baseline and during treatment. Specifically, associations of individual positions with the HBeAg-negative phenotype were studied, as well as the association of the most prevalent mutations with combined response in HBeAg-positive and -negative patients at week 72 (HBeAg negativity, HBV-DNA <2000 IU/mL and ALT normalization at 24 weeks of treatment-free follow-up).

RESULTS

The mutations most strongly correlated with the HBeAg-negative phenotype were at positions 1762/1764 and 1896/1899 in the BCP and PC region, respectively. No major changes in nucleotide composition of these positions were observed during treatment. In HBeAg-negative patients, a combined presence of 1764A and 1896A was correlated with lower ALT levels (p = 0.004), whereas the presence of 1899A was correlated with higher age (p = 0.030), lower HBV-DNA level (p = 0.036), and previous IFN therapy (p = 0.032). The presence of 1764A/1896A or the absence of 1899A at baseline, was associated with lower response rates, after adjustment for HBV genotype (p = 0.031 and p = 0.017) and HBsAg level (p = 0.035 and p = 0.022).

CONCLUSION

We identified novel correlations between common BCP and PC variants with response to Peg-IFN and adefovir in HBeAg-negative patients. Ultimately, this may guide the selection of those patients most likely to benefit from Peg-IFN-based treatment.

A comparison of 454 sequencing and clonal sequencing for the characterization of hepatitis C virus NS3 variants

We compared 454 amplicon sequencing with clonal sequencing for the characterization of intra-host hepatitis C virus (HCV) NS3 variants. Clonal and 454 sequences were obtained from 12 patients enrolled in a clinical phase I study for telaprevir, an NS3-4a protease inhibitor. Thirty-nine datasets were used to compare the consensus sequence, average pairwise distance, normalized Shannon entropy, phylogenetic tree topology and the number and frequency of variants derived from both sequencing techniques. In general, a good concordance was observed between both techniques for the majority of datasets. Discordant results were observed for 5 out of 39 clonal and 454 datasets, which could be attributed to primer-related selective amplification used for clonal sequencing. Both 454 and clonal datasets consisted of a few major variants and a large number of low-frequency variants. Telaprevir resistance-associated variants were observed in low frequencies and were detected more often by 454. We conclude that performance of 454 and clonal sequencing is comparable for the characterization of intra-host virus populations. Not surprisingly, 454 is superior for the detection of low frequency resistance-associated variants. However, despite the greater coverage, 454 failed to detect some low frequency variants detected by clonal sequencing.

Improved detection of artifactual viral minority variants in high-throughput sequencing data

High-throughput sequencing (HTS) of viral samples provides important information on the presence of viral minority variants. However, detection and accurate quantification is limited by the capacity to distinguish biological from artificial variation. In this study, errors related to the Illumina HiSeq2000 library generation and HTS process were investigated by determining minority variant frequencies in an influenza A/WSN/1933(H1N1) virus reverse-genetics plasmid pool. Errors related to amplification and sequencing were determined using the same plasmid pool, by generation of infectious virus using reverse genetics followed by in duplo reverse-transcriptase PCR (RT-PCR) amplification and HTS in the same sequence run. Results showed that after “best practice” quality control (QC), within the plasmid pool, one minority variant with a frequency >0.5% was identified, while 84 and 139 were identified in the RT-PCR amplified samples, indicating RT-PCR amplification artificially increased variation. Detailed analysis showed that artifactual minority variants could be identified by two major technical characteristics: their predominant presence in a single read orientation and uneven distribution of mismatches over the length of the reads. We demonstrate that by addition of two QC steps 95% of the artifactual minority variants could be identified. When our analysis approach was applied to three clinical samples 68% of the initially identified minority variants were identified as artifacts. Our study clearly demonstrated that, without additional QC steps, overestimation of viral minority variants is very likely to occur, mainly as a consequence of the required RT-PCR amplification step. The improved ability to detect and correct for artifactual minority variants, increases data resolution and could aid both past and future studies incorporating HTS. The source code has been made available through Sourceforge (https://sourceforge.net/projects/mva-ngs).

Pandemic H1N1 virus transmission and shedding dynamics in index case households of a prospective Vietnamese cohort

Objectives

Influenza household transmission studies are required to guide prevention strategies but most passively recruit index cases that seek healthcare. We investigated A(H1N1)pdm09 transmission in a household-based cohort during 2009.

Methods

Health-workers visited 270 households weekly, and collected swabs from influenza-like-illness cases. If A(H1N1)pdm09 was RT-PCR-confirmed, all household members had symptoms assessed and swabs collected daily for 10–15 days. Viral RNA was quantified and sequenced and serology performed on pre-pandemic sera.

Results

Index cases were detected in 20 households containing 81 people. 98.5% lacked A(H1N1)pdm09 neutralizing antibodies in pre-pandemic sera. Eleven (18.6%, 95% CI 10.7–30.4%) of 59 contacts were infected. Virus genetic diversity within households was negligible and less than between households. Index and secondary cases were distributed between mothers, daughters and sons, and had similar virus-RNA shedding and symptom dynamics. Fathers were rarely infected. Five secondary cases (45%) had no apparent symptoms and three shed virus before symptoms. Secondary infection was associated with index case wet cough (OR 1.56, 95% CI 1.22–1.99).

Conclusions

In this cohort of A(H1N1)pdm09 susceptible persons, virus sequencing was capable of discriminating household from community transmission. Household transmission involved mothers and children but rarely fathers. Asymptomatic or pre-symptomatic shedding was common.

Emergence of the virulence-associated PB2 E627K substitution in a fatal human case of highly pathogenic avian influenza virus A(H7N7) infection as determined by Illumina ultra-deep sequencing

Avian influenza viruses are capable of crossing the species barrier and infecting humans. Although evidence of human-to-human transmission of avian influenza viruses to date is limited, evolution of variants toward more-efficient human-to-human transmission could result in a new influenza virus pandemic. In both the avian influenza A(H5N1) and the recently emerging avian influenza A(H7N9) viruses, the polymerase basic 2 protein (PB2) E627K mutation appears to be of key importance for human adaptation. During a large influenza A(H7N7) virus outbreak in the Netherlands in 2003, the A(H7N7) virus isolated from a fatal human case contained the PB2 E627K mutation as well as a hemagglutinin (HA) K416R mutation. In this study, we aimed to investigate whether these mutations occurred in the avian or the human host by Illumina Ultra-Deep sequencing of three previously uninvestigated clinical samples obtained from the fatal case. In addition, we investigated three chicken samples, two of which were obtained from the source farm. Results showed that the PB2 E627K mutation was not present in any of the chicken samples tested. Surprisingly, the avian samples were characterized by the presence of influenza virus defective RNA segments, suggestive for the synthesis of defective interfering viruses during infection in poultry. In the human samples, the PB2 E627K mutation was identified with increasing frequency during infection. Our results strongly suggest that human adaptation marker PB2 E627K has emerged during virus infection of a single human host, emphasizing the importance of reducing human exposure to avian influenza viruses to reduce the likelihood of viral adaptation to humans.

The presence of CXCR4-using HIV-1 prior to start of antiretroviral therapy is an independent predictor of delayed viral suppression

The emergence of CXCR4-using HIV variants (X4-HIV) is associated with accelerated disease progression in the absence of antiretroviral therapy. However, the effect of X4-HIV variants on the treatment response remains unclear. Here we determined whether the presence of X4-HIV variants influenced the time to undetectable viral load and CD4+ T cell reconstitution after initiation of cART in 732 patients. The presence of X4-HIV variants was determined by MT-2 assay prior to cART initiation and viral load and CD4+ T cell counts were analyzed every 3 to 6 months during a three year follow-up period. Kaplan-Meier and Cox proportional hazard analyses were performed to compare time to viral suppression and the absolute CD4+ T cell counts and increases in CD4+ T cell counts during follow-up were compared for patients with and without X4-HIV at start of cART. Patients harboring X4-HIV variants at baseline showed a delay in time to achieve viral suppression below the viral load detection limit. This delay in viral suppression was independently associated with high viral load and the presence of X4-HIV variants. Furthermore, the absolute CD4+ T cell counts were significantly lower in patients harboring X4-HIV variants at all time points during follow-up. However, no differences were observed in the increase in absolute CD4+ T cell numbers after treatment initiation, indicating that the reconstitution of CD4+ T cells is independent of the presence of X4-HIV variants. The emergence of X4-HIV has been associated with an accelerated CD4+ T cell decline during the natural course of infection and therefore, patients who develop X4-HIV variants may benefit from earlier treatment initiation in order to obtain faster reconstitution of the CD4+ T cell population to normal levels.

Viral population analysis and minority-variant detection using short read next-generation sequencing

RNA viruses within infected individuals exist as a population of evolutionary-related variants. Owing to evolutionary change affecting the constitution of this population, the frequency and/or occurrence of individual viral variants can show marked or subtle fluctuations. Since the development of massively parallel sequencing platforms, such viral populations can now be investigated to unprecedented resolution. A critical problem with such analyses is the presence of sequencing-related errors that obscure the identification of true biological variants present at low frequency. Here, we report the development and assessment of the Quality Assessment of Short Read (QUASR) Pipeline (http://sourceforge.net/projects/quasr) specific for virus genome short read analysis that minimizes sequencing errors from multiple deep-sequencing platforms, and enables post-mapping analysis of the minority variants within the viral population. QUASR significantly reduces the error-related noise in deep-sequencing datasets, resulting in increased mapping accuracy and reduction of erroneous mutations. Using QUASR, we have determined influenza virus genome dynamics in sequential samples from an in vitro evolution of 2009 pandemic H1N1 (A/H1N1/09) influenza from samples sequenced on both the Roche 454 GSFLX and Illumina GAIIx platforms. Importantly, concordance between the 454 and Illumina sequencing allowed unambiguous minority-variant detection and accurate determination of virus population turnover in vitro.

The human severe acute respiratory syndrome coronavirus (SARS-CoV) 8b protein is distinct from its counterpart in animal SARS-CoV and down-regulates the expression of the envelope protein in infected cells

The severe acute respiratory syndrome coronavirus (SARS-CoV), isolated from humans infected during the peak of epidemic, encodes two accessory proteins termed as 8a and 8b. Interestingly, the SARS-CoV isolated from animals contains an extra 29-nucleotide in this region such that these proteins are fused to become a single protein, 8ab. Here, we compared the cellular properties of the 8a, 8b and 8ab proteins by examining their cellular localizations and their abilities to interact with other SARS-CoV proteins. These results may suggest that the conformations of 8a and 8b are different from 8ab although nearly all the amino acids in 8a and 8b are found in 8ab. In addition, the expression of the structural protein, envelope (E), was down-regulated by 8b but not 8a or 8ab. Consequently, E was not detectable in SARS-CoV-infected cells that were expressing high levels of 8b. These findings suggest that 8b may modulate viral replication and/or pathogenesis.