Incursion of SARS-CoV-2 BA.2.86.1 variant into Israel: National-scale wastewater surveillance using a novel quantitative real-time PCR assay

The emergence of the SARS-CoV-2 variant BA.2.86.1 raised a considerable concern, due to the large number of potentially virulent mutations. In this study, we developed a novel assay that specifically detects variant BA.2.86.1, and used it to screen environmental samples from wastewater treatment sites across Israel. By using a multiplex assay that included a general SARS-CoV-2 reaction, together with the BA.2.86.1-specific reaction and a control reaction, we quantified the absolute number of viral copies in each sample and its relative abundance, compared with the total copy number of circulating SARS-CoV-2. Evaluation of the new reactions showed that they are both sensitive and specific, detecting down to four copies per reaction, and maintain specificity in the presence of Omicron variants BA.1, 2 and 4 RNA. Examination of 279 samples from 30 wastewater collection sites during August-September 2023 showed that 35 samples (12.5 %) were positive, from 18 sites. Quantitative analysis of the samples showed that the relative abundance of variant BA.2.86.1 with respect to the total viral load of SARS-CoV-2 was very low and consisted between 0.01 % and 0.6 % of the total SARS-CoV-2 circulation. This study demonstrates the importance of combining wastewater surveillance with the development of specialized diagnostic assays, when clinical testing is insufficient. This approach may be useful for timely response by public health authorities in future outbreaks.

Rapid molecular epidemiology investigations into two recent measles outbreaks in Israel detected from October 2023 to January 2024

Between late 2023 and early 2024, two measles outbreaks occurred in Israel, each caused by importation of measles virus strains of respective B3 and D8 genotypes. In this study, we validate transmission pathways uncovered by epidemiological investigations using a rapid molecular approach, based on complete measles virus genomes. The presented findings support this rapid molecular approach in complementing conventional contact tracing and highlight its potential for informing public health interventions.

Nirmatrelvir Resistance-de Novo E166V/L50V Mutations in an Immunocompromised Patient Treated With Prolonged Nirmatrelvir/Ritonavir Monotherapy Leading to Clinical and Virological Treatment Failure-a Case Report

Resistance of SARS-CoV-2 to antivirals was shown to develop in immunocompromised individuals receiving remdesivir. We describe an immunocompromised patient who was treated with repeated and prolonged courses of nirmatrelvir and developed de-novo E166V/L50F mutations in the Mpro region. These mutations were associated with clinical and virological treatment failure.

Factors Associated with Virological Failure in First-Line Antiretroviral Therapy in Patients Diagnosed with HIV-1 between 2010 and 2018 in Israel

Despite the progress in contemporary antiretroviral therapy (ART) and the continuous changes in treatment guidelines, virological failure (VF) is still an ongoing concern. The goal of this study was to assess factors related to VF after first-line ART. A longitudinal cohort retrospective study of individuals on first-line ART diagnosed with HIV-1 in 2010-2018 and followed-up for a median of two years was conducted. Demographics, baseline and longitudinal CD4 counts, treatment regimens, adherence and VF were recorded. The Cox proportional hazards regression and mixed models were used. A cohort of 1130 patients were included. Overall, 80% were males and 62% were Israeli-born individuals. Compared to individuals diagnosed in 2010-2014, when treatment was initiated according to CD4 levels, those diagnosed in 2015-2018 were older and had lower baseline CD4 counts. VF was recorded in 66 (5.8%) patients. Diagnosis with CD4 <200 cells/mmᶟ with AIDS-defining conditions (HR = 2.75, 95%CI:1.52-4.97, p < 0.001) and non-integrase strand transfer inhibitor regimens (non-INSTI, HR = 1.80, 95%CI:1.01-3.24, p = 0.047) increased VF risk. No impact of baseline resistance was observed. We concluded that the early detection of HIV-1 infection and usage of INSTI-based regimens are recommended to reduce VF.

Dengue Types 1 and 3 Identified in Travelers Returning from Kathmandu, Nepal, during the October 2022 Outbreak Are Related to Strains Recently Identified in India

Phylogenetic analysis of dengue serotypes 1 and 3, which were diagnosed in travelers and Nepalese infected in Kathmandu during the October 2022 outbreak, revealed that both serotypes were clustered closest to the sequences sampled in India. This suggests both serotypes may have originated in India.

Environmental surveillance detected type 3 vaccine-derived polioviruses in increasing frequency at multiple sites prior to detection of a poliomyelitis case

Israel conducts routine environmental (15 sites) and acute flaccid paralysis (AFP) surveillance for poliovirus. During September 2021, increasing numbers of wastewater samples collected from more than one site in the Jerusalem region proved positive for ambiguous type 3 vaccine-derived poliovirus (aVDPV3), while environmental samples from remaining sampling sites were negative. In late February 2022, a VDPV3, genetically related to the Jerusalem environmental surveillance samples, was isolated from a stool sample collected from a non-immunodeficient, non-immunized child from Jerusalem who developed AFP, indicating that the aVDPV3s were circulating (cVDPV3s) rather than immunodeficiency-related VDPV3s (iVDPVs). In response to these isolations, the Israel Ministry of Health launched a catch-up immunization program.

Global disparities in SARS-CoV-2 genomic surveillance

Genomic sequencing is essential to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments, vaccines, and guide public health responses. To investigate the global SARS-CoV-2 genomic surveillance, we used sequences shared via GISAID to estimate the impact of sequencing intensity and turnaround times on variant detection in 189 countries. In the first two years of the pandemic, 78% of high-income countries sequenced >0.5% of their COVID-19 cases, while 42% of low- and middle-income countries reached that mark. Around 25% of the genomes from high income countries were submitted within 21 days, a pattern observed in 5% of the genomes from low- and middle-income countries. We found that sequencing around 0.5% of the cases, with a turnaround time <21 days, could provide a benchmark for SARS-CoV-2 genomic surveillance. Socioeconomic inequalities undermine the global pandemic preparedness, and efforts must be made to support low- and middle-income countries improve their local sequencing capacity.

Recurrent congenital cytomegalovirus infection in a sequential pregnancy with severe sequelae, and a possible association with prophylactic valacyclovir treatment: a case report

Recurrent congenital cytomegalovirus infections in consecutive pregnancies are rarely reported. Due to the risk of fetal infection from preconception maternal infection, a 6-month interval after primary maternal infection is generally advised before a new conception. Recently, high-dose valacyclovir treatment was shown to prevent fetal infection in first trimester primary infections. We present a case of first trimester primary infection treated with high-dose valacyclovir but resulting in polymerase chain reaction-confirmed fetal infection. Cytomegalovirus-specific immunoglobulin G titers remained very low during treatment and rose only after cessation of antiviral treatment. Six months after primary seroconversion, in a sequential pregnancy, recurrent fetal infection was diagnosed and resulted in severe fetal sequella. Whole genome sequencing of both amniotic fluid isolates proved them to be identical. Both pregnancies were terminated. We hypothesize that valacyclovir treatment, although unsuccessful in preventing fetal infection, had delayed the adaptive maternal immune response and might have contributed to fetal infection during the sequential pregnancy. We suggest that a longer delay might be warranted after valacyclovir treatment and before a new conception.

Emergence of genetically linked vaccine-originated poliovirus type 2 in the absence of oral polio vaccine, Jerusalem, April to July 2022

We report an emergence and increase in poliovirus type 2 detection via routine wastewater surveillance in three non-overlapping regions in the Jerusalem region, Israel, between April and July 2022. Sequencing showed genetic linkage among isolates and accumulation of mutations over time, with two isolates defined as vaccine-derived polioviruses (VDPV). This demonstrates the emergence and potential circulation of type 2 VDPV in a high-income country with high vaccine coverage and underscores the importance of routine wastewater surveillance during the polio eradication.

HIV-1 Circulating Recombinant Forms (CRFs) and Unique Recombinant Forms (URFs) in Israel, 2010-2018

Monitoring HIV-1 circulating recombinant forms (CRFs) and unique recombinant forms (URFs) is important for disease surveillance. Recombination may affect prevention efforts and interfere with the diagnosis and treatment of HIV-1 infection. Here, we characterized the epidemiology of HIV-1 CRFs and URFs in Israel. Partial pol sequences from treatment naïve patients diagnosed in 2010−2018 were assessed using the recombinant identification program (RIP), the recombinant detection program (RDP5), and using the maximum-likelihood phylogenetic method, using 410 reference sequences obtained from the Los Alamos database. CRFs and URFs were identified in 11% (213/1940) of all sequenced cases. The median age at diagnosis was 38 (30−47) years, 61% originated from Israel, and 82% were male. The most common were CRF02_AG (30.5%), CRF01_AE (16.9%), and the more complex forms CRF01_AE/CRF02_AG/A3 (10.8%) and B/F1 (7%). A significant increase in their overall proportion was observed in recent years (8.1% in 2010−2012, 20.3% in 2016−2018, p < 0.001). This increase was most prominent in individuals carrying CRF02_AG (2.5% in 2010−2015, 9.8% in 2016−2018, p < 0.001). Men who have sex with men (MSM) was the most common risk group; however, those infected with the secondary recombinant CRF02_AG/A6 were mainly injecting drug users (IDUs). The most common resistance mutations were K103N (5/213, 2.3%) and E138A (18/213, 8.5%) in the reverse transcriptase. Only E138A was more frequent in the recombinants compared with the classic subtypes and was significantly associated with a specific secondary CRF, CRF02_AG/A4. We concluded that CRFs and URFs were mainly detected in Israeli-born MSM and that an increase in the overall proportion of such HIV-1 sequences could be observed in more recent years.

National Scale Real-Time Surveillance of SARS-CoV-2 Variants Dynamics by Wastewater Monitoring in Israel

In this report, we describe a national-scale monitoring of the SARS-CoV-2 (SC-2) variant dynamics in Israel, using multiple-time sampling of 13 wastewater treatment plants. We used a combination of inclusive and selective quantitative PCR assays that specifically identify variants A19/A20 or B.1.1.7 and tested each sample for the presence and relative viral RNA load of each variant. We show that between December 2020 and March 2021, a complete shift in the SC-2 variant circulation was observed, where the B.1.1.7 replaced the A19 in all examined test points. We further show that the normalized viral load (NVL) values and the average new cases per week reached a peak in January 2021 and then decreased gradually in almost all test points, in parallel with the progression of the national vaccination campaign, during February-March 2021. This study demonstrates the importance of monitoring SC-2 variant by using a combination of inclusive and selective PCR tests on a national scale through wastewater sampling, which is far more amendable for high-throughput monitoring compared with sequencing. This approach may be useful for real-time dynamics surveillance of current and future variants, such as the Omicron (BA.1, BA.2) and other variants.

The impact of SARS-CoV-2 on respiratory syndromic and sentinel surveillance in Israel, 2020: a new perspective on established systems

Background

The COVID-19 pandemic presented new challenges for the existing respiratory surveillance systems, and adaptations were implemented. Systematic assessment of the syndromic and sentinel surveillance platforms during the pandemic is essential for understanding the value of each platform in the context of an emerging pathogen with rapid global spread.

Aim

We aimed to evaluate systematically the performance of various respiratory syndromic surveillance platforms and the sentinel surveillance system in Israel from 1 January to 31 December 2020.

Methods

We compared the 2020 syndromic surveillance trends to those of the previous 3 years, using Poisson regression adjusted for overdispersion. To assess the performance of the sentinel clinic system as compared with the national SARS-CoV-2 repository, a cubic spline with 7 knots and 95% confidence intervals were applied to the sentinel network's weekly percentage of positive SARS-CoV-2 cases.

Results

Syndromic surveillance trends changed substantially during 2020, with a statistically significant reduction in the rates of visits to physicians and emergency departments to below previous years' levels. Morbidity patterns of the syndromic surveillance platforms were inconsistent with the progress of the pandemic, while the sentinel surveillance platform was found to reflect the national circulation of SARS-CoV-2 in the population.

Conclusion

Our findings reveal the robustness of the sentinel clinics platform for the surveillance of the main respiratory viruses during the pandemic and possibly beyond. The robustness of the sentinel clinics platform during 2020 supports its use in locations with insufficient resources for widespread testing of respiratory viruses.

HiSpike Method for High-Throughput Cost Effective Sequencing of the SARS-CoV-2 Spike Gene

The changing nature of the SARS-CoV-2 pandemic poses unprecedented challenges to the world's health systems. Emerging spike gene variants jeopardize global efforts to produce immunity and reduce morbidity and mortality. These challenges require effective real-time genomic surveillance solutions that the medical community can quickly adopt. The SARS-CoV-2 spike protein mediates host receptor recognition and entry into the cell and is susceptible to generation of variants with increased transmissibility and pathogenicity. The spike protein is the primary target of neutralizing antibodies in COVID-19 patients and the most common antigen for induction of effective vaccine immunity. Tight monitoring of spike protein gene variants is key to mitigating COVID-19 spread and generation of vaccine escape mutants. Currently, SARS-CoV-2 sequencing methods are labor intensive and expensive. When sequence demands are high sequencing resources are quickly exhausted. Consequently, most SARS-CoV-2 strains are sequenced in only a few developed countries and rarely in developing regions. This poses the risk that undetected, dangerous variants will emerge. In this work, we present HiSpike, a method for high-throughput cost effective targeted next generation sequencing of the spike gene. This simple three-step method can be completed in < 30 h, can sequence 10-fold more samples compared to conventional methods and at a fraction of their cost. HiSpike has been validated in Israel, and has identified multiple spike variants from real-time field samples including Alpha, Beta, Delta and the emerging Omicron variants. HiSpike provides affordable sequencing options to help laboratories conserve resources for widespread high-throughput, near real-time monitoring of spike gene variants.

Epidemiology and Transmitted HIV-1 Drug Resistance among Treatment-Naïve Individuals in Israel, 2010-2018

Despite the low prevalence of HIV-1 in Israel, continuous waves of immigration may have impacted the local epidemic. We characterized all people diagnosed with HIV-1 in Israel in 2010–2018. The demographics and clinical data of all individuals (n = 3639) newly diagnosed with HIV-1 were retrieved. Subtypes, transmitted drug-resistance mutations (TDRM), and phylogenetic relations, were determined in >50% of them. In 39.1%, HIV-1 transmission was through heterosexual contact; 34.3% were men who have sex with men (MSM); and 10.4% were people who inject drugs. Many (>65%) were immigrants. Israeli-born individuals were mostly (78.3%) MSM, whereas only 9% of those born in Sub-Saharan Africa (SSA), Eastern Europe and Central Asia (EEU/CA), were MSM. The proportion of individuals from SSA decreased through the years 2010–2018 (21.1% in 2010–2012; 16.8% in 2016–2018) whereas those from EEU/CA increased significantly (21% in 2010–2012; 27.8% in 2016–2018, p < 0.001). TDRM were identified in 12.1%; 3.7, 3.3 and 6.6% had protease inhibitors (PI), nucleotide reverse transcriptase inhibitors (NRTI), and non-nucleoside reverse transcriptase inhibitors (NNRTI) TDRM, respectively, with the overall proportion remaining stable in the studied years. None had integrase TDRM. Subtype B was present in 43.9%, subtype A in 25.2% (A6 in 22.8 and A1 in 2.4%) and subtype C in 17.1% of individuals. Most MSM had subtype B. Subtype C carriers formed small clusters (with one unexpected MSM cluster), A1 formed a cluster mainly of locally-born patients with NNRTI mutations, and A6 formed a looser cluster of individuals mainly from EEU. Israelis, <50 years old, carrying A1, had the highest risk for having TDRM. In conclusion, an increase in immigrants from EEU/CA and a decrease in those from SSA characterized the HIV-1 epidemic in 2010–2018. Baseline resistance testing should still be recommended to identify TDRM, and improve surveillance and care.

Global disparities in SARS-CoV-2 genomic surveillance

Genomic sequencing provides critical information to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments and vaccines, and guide public health responses. To investigate the spatiotemporal heterogeneity in the global SARS-CoV-2 genomic surveillance, we estimated the impact of sequencing intensity and turnaround times (TAT) on variant detection in 167 countries. Most countries submit genomes >21 days after sample collection, and 77% of low and middle income countries sequenced <0.5% of their cases. We found that sequencing at least 0.5% of the cases, with a TAT <21 days, could be a benchmark for SARS-CoV-2 genomic surveillance efforts. Socioeconomic inequalities substantially impact our ability to quickly detect SARS-CoV-2 variants, and undermine the global pandemic preparedness.

Effectiveness of BNT162b2 mRNA COVID-19 vaccine against SARS-CoV-2 variant Beta (B.1.351) among persons identified through contact tracing in Israel: A prospective cohort study

BACKGROUND

SARS-CoV-2 variant Beta (B.1.351) was designated as a Variant of Concern (VoC) after becoming the dominant strain in South Africa and spreading internationally. BNT162b2 showed lower levels of neutralizing antibodies against Beta than against other strains raising concerns about effectiveness of vaccines against infections caused by Beta. We estimated BNT162b2 vaccine effectiveness (VE) against Beta infections in Israel, a country with high vaccine uptake.

METHODS

The Ministry of Health (MoH) identified Beta cases through mandatory reporting of SARS-CoV-2 cases and whole genome sequencing (WGS) of specimens from vaccination-breakthrough infections, reinfections, arriving international travelers, and a selection of other infected persons. A cohort analysis was conducted of exposure events of contacts of primary Beta cases. WGS was conducted on available PCR-positive specimens collected from contacts. VE estimates with 95% confidence intervals (CIs) against confirmed and probable Beta infections were determined by comparing infection risk between unvaccinated and fully-vaccinated (≥7 days after the second dose) contacts, and between unvaccinated and partially-vaccinated (<7 days after the second dose) contacts.

FINDINGS

MoH identified 310 Beta cases through Jun 27, 2021. During the study period (Dec 11, 2020 - Mar 25, 2021), 164 non-institutionalized primary Beta cases, with 552 contacts aged ≥16 years, were identified. 343/552 (62%) contacts were interviewed and tested. 71/343 (21%) contacts were PCR-positive. WGS was performed on 7/71 (10%) PCR-positive specimens; all were Beta. Among SARS-CoV-2-infected contacts, 48/71 (68%) were symptomatic, 10/71 (14%) hospitalized, and 2/71 (3%) died. Fully-vaccinated VE against confirmed or probable Beta infections was 72% (95% CI -5 - 97%; p=0·04) and against symptomatic confirmed or probable Beta infections was 100% (95% CI 19 - 100%; p=0·01). There was no evidence of protection in partially-vaccinated contacts.

INTERPRETATION

In a prospective observational study, two doses of BNT162b2 were effective against confirmed and probable Beta infections. Through the end of June 2021, introductions of Beta did not interrupt control of the pandemic in Israel.

FUNDING

Israel Ministry of Health and Pfizer.

Rapid and High-Throughput Reverse Transcriptase Quantitative PCR (RT-qPCR) Assay for Identification and Differentiation between SARS-CoV-2 Variants B.1.1.7 and B.1.351

Emerging SARS-CoV-2 (SC-2) variants with increased infectivity and vaccine resistance are of major concern. Rapid identification of such variants is important for the public health decision making and to provide valuable data for epidemiological and policy decision making. We developed a multiplex reverse transcriptase quantitative PCR (RT-qPCR) assay that can specifically identify and differentiate between the emerging B.1.1.7 and B.1.351 SC-2 variants. In a single assay, we combined four reactions-one that detects SC-2 RNA independently of the strain, one that detects the D3L mutation, which is specific to variant B.1.1.7, one that detects the 242 to 244 deletion, which is specific to variant B.1.351, and the fourth reaction, which identifies the human RNAseP gene, serving as an endogenous control for RNA extraction integrity. We show that the strain-specific reactions target mutations that are strongly associated with the target variants and not with other major known variants. The assay's specificity was tested against a panel of respiratory pathogens (n = 16), showing high specificity toward SC-2 RNA. The assay's sensitivity was assessed using both in vitro transcribed RNA and clinical samples and was determined to be between 20 and 40 viral RNA copies per reaction. The assay performance was corroborated with Sanger and whole-genome sequencing, showing complete agreement with the sequencing results. The new assay is currently implemented in the routine diagnostic work at the Central Virology Laboratory, and may be used in other laboratories to facilitate the diagnosis of these major worldwide-circulating SC-2 variants. IMPORTANCE This study describes the design and utilization of a multiplex reverse transcriptase quantitative PCR (RT-qPCR) to identify SARS-COV-2 (SC2) RNA in general and, specifically, to detect whether it is of lineage B.1.1.7 or B.1.351. Implementation of this method in diagnostic and research laboratories worldwide may help the efforts to contain the COVID-19 pandemic. The method can be easily scaled up and be used in high-throughput laboratories, as well as small ones. In addition to immediate help in diagnostic efforts, this method may also help in epidemiological studies focused on the spread of emerging SC-2 lineages.

Rapid and High-Throughput Reverse Transcriptase Quantitative PCR (RT-qPCR) Assay for Identification and Differentiation between SARS-CoV-2 Variants B.1.1.7 and B.1.351

Emerging SARS-CoV-2 (SC-2) variants with increased infectivity and vaccine resistance are of major concern. Rapid identification of such variants is important for the public health decision making and to provide valuable data for epidemiological and policy decision making. We developed a multiplex reverse transcriptase quantitative PCR (RT-qPCR) assay that can specifically identify and differentiate between the emerging B.1.1.7 and B.1.351 SC-2 variants. In a single assay, we combined four reactions-one that detects SC-2 RNA independently of the strain, one that detects the D3L mutation, which is specific to variant B.1.1.7, one that detects the 242 to 244 deletion, which is specific to variant B.1.351, and the fourth reaction, which identifies the human RNAseP gene, serving as an endogenous control for RNA extraction integrity. We show that the strain-specific reactions target mutations that are strongly associated with the target variants and not with other major known variants. The assay's specificity was tested against a panel of respiratory pathogens (n = 16), showing high specificity toward SC-2 RNA. The assay's sensitivity was assessed using both in vitro transcribed RNA and clinical samples and was determined to be between 20 and 40 viral RNA copies per reaction. The assay performance was corroborated with Sanger and whole-genome sequencing, showing complete agreement with the sequencing results. The new assay is currently implemented in the routine diagnostic work at the Central Virology Laboratory, and may be used in other laboratories to facilitate the diagnosis of these major worldwide-circulating SC-2 variants. IMPORTANCE This study describes the design and utilization of a multiplex reverse transcriptase quantitative PCR (RT-qPCR) to identify SARS-COV-2 (SC2) RNA in general and, specifically, to detect whether it is of lineage B.1.1.7 or B.1.351. Implementation of this method in diagnostic and research laboratories worldwide may help the efforts to contain the COVID-19 pandemic. The method can be easily scaled up and be used in high-throughput laboratories, as well as small ones. In addition to immediate help in diagnostic efforts, this method may also help in epidemiological studies focused on the spread of emerging SC-2 lineages.

BNT162b2 vaccine effectiveness was marginally affected by the SARS-CoV-2 beta variant in fully vaccinated individuals

Objective

To evaluate the effectiveness of the Pfizer BNT162b2 vaccine against the SARS-Cov-2 Beta variant.

Study Design and Setting

Israel's mass vaccination program, using two doses of the Pfizer BNT162b2 vaccine, successfully curtailed the Alpha variant outbreak during winter 2020–2021, However, the virus may mutate and partially evade the immune system. To monitor this, sequencing of selected positive swab samples of interest was initiated. Comparing vaccinated with unvaccinated PCR positive persons, we estimated the odds ratio for a vaccinated case to have the Beta vs. the Alpha variant, using logistic regression, controlling for important confounders.

Results

There were 19 cases of Beta variant (3.2%) among those vaccinated more than 14 days before the positive sample and 79 (3.4%) among the unvaccinated. The estimated odds ratio was 1.26 (95% CI: 0.65–2.46). Assuming the effectiveness against the Alpha variant to be 95%, the estimated effectiveness against the Beta variant was 94% (95% CI: 88%–98%).

Conclusion

Despite concerns over the Beta variant, the BNT162b2 vaccine seemed to provide substantial immunity against both the Beta and the Alpha variants. From 14 days following the second vaccine dose, the effectiveness of BNT162b2 vaccine was at most marginally affected by the Beta variant.

Nosocomial outbreak caused by the SARS-CoV-2 Delta variant in a highly vaccinated population, Israel, July 2021

A nosocomial outbreak of SARS-CoV-2 Delta variant infected 42 patients, staff and family members; 39 were fully vaccinated. The attack rate was 10.6% (16/151) among exposed staff and reached 23.7% (23/97) among exposed patients in a highly vaccinated population, 16–26 weeks after vaccination (median: 25 weeks). All cases were linked and traced to one patient. Several transmissions occurred between individuals wearing face masks. Fourteen of 23 patients became severely sick or died, raising a question about possible waning immunity.

Detection of SARS-CoV-2 variants by genomic analysis of wastewater samples in Israel

Investigation of SARS-CoV-2 spread and identification of variants in sewers has been demonstrated to accurately detect prevalence of viral strains and is advantageous to clinical sampling in population catchment size. Herein, we utilized an established nationwide system of wastewater sampling and viral concentration approaches to perform large-scale surveillance of SARS-CoV-2 variants in nine different locations across Israel that were sampled from August 2020 to February 2021 and sequenced (n = 58). Viral sequences obtained from the wastewater samples had high coverages of the genome, and mutation analyses successfully identified the penetration of the B.1.1.7 variant into Israel in December 2020 in the central and north regions, and its spread into additional regions in January and February 2021, corresponding with clinical sampling results. Moreover, the wastewater analysis identified the B.1.1.7 variant in December 2020 in regions in which non-sufficient clinical sampling was available. Other variants of concern examined, including P.1 (Brazil/Manaus), B.1.429 (USA/California), B.1.526 (USA/New York), A.23.1 (Uganda) and B.1.525 (Unknown origin), did not show consistently elevated frequencies. This study exemplifies that surveillance by sewage is a robust approach which allows to monitor the diversity of SARS-CoV-2 strains circulating in the community. Most importantly, this approach can pre-identify the emergence of epidemiologically or clinically relevant mutations/variants, aiding in public health decision making.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2 with grinch

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

The Rise and Fall of a Local SARS-CoV-2 Variant with the Spike Protein Mutation L452R

Emerging SARS-CoV-2 variants may threaten global vaccination efforts and the awaited reduction in outbreak burden. In this study, we report a novel variant carrying the L452R mutation that emerged from a local B.1.362 lineage, B.1.362+L452R. The L452R mutation is associated with the Delta and Epsilon variants and was shown to cause increased infection and reduction in neutralization in pseudoviruses. Indeed, the B.1.362+L452R variant demonstrated a X4-fold reduction in neutralization capacity of sera from BNT162b2-vaccinated individuals compared to a wild-type strain. The variant infected 270 individuals in Israel between December 2020 and March 2021, until diminishing due to the gain in dominance of the Alpha variant in February 2021. This study demonstrates an independent, local emergence of a variant carrying a critical mutation, L452R, which may have the potential of becoming a variant of concern and emphasizes the importance of routine surveillance and detection of novel variants among efforts undertaken to prevent further disease spread.

The Rise and Fall of a Local SARS-CoV-2 Variant with the Spike Protein Mutation L452R

Emerging SARS-CoV-2 variants may threaten global vaccination efforts and the awaited reduction in outbreak burden. In this study, we report a novel variant carrying the L452R mutation that emerged from a local B.1.362 lineage, B.1.362+L452R. The L452R mutation is associated with the Delta and Epsilon variants and was shown to cause increased infection and reduction in neutralization in pseudoviruses. Indeed, the B.1.362+L452R variant demonstrated a X4-fold reduction in neutralization capacity of sera from BNT162b2-vaccinated individuals compared to a wild-type strain. The variant infected 270 individuals in Israel between December 2020 and March 2021, until diminishing due to the gain in dominance of the Alpha variant in February 2021. This study demonstrates an independent, local emergence of a variant carrying a critical mutation, L452R, which may have the potential of becoming a variant of concern and emphasizes the importance of routine surveillance and detection of novel variants among efforts undertaken to prevent further disease spread.

Prolonged detection of complete viral genomes demonstrated by SARS-CoV-2 sequencing of serial respiratory specimens

Accurate and timely diagnosis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is clinically essential, and is required also to monitor confirmed cases aiming to prevent further spread. Positive real-time PCR results at late time points following initial diagnosis may be clinically misleading as this methodology cannot account for the infection capabilities and the existence of whole genome sequences. In this study, 47 serial respiratory samples were tested by Allplex-nCoV test (Seegene), a triplex of three assays targeting the SARS-CoV-2 RdRP, E and N genes and subsequently assessed by next generation sequencing (NGS). COVID19 patients were tested at an early stage of the disease, when all these viral gene targets were positive, and at an advanced stage, when only the N gene target was positive in the Allplex-nCoV test. The corresponding NGS results showed the presence of complete viral genome copies at both early and advanced stages of the disease, although the total number of mapped sequences was lower in samples from advanced disease stages. We conclude that reduced viral transmission at this late disease stage may result from the low quantities of complete viral sequences and not solely from transcription favoring the N gene.

A Unique SARS-CoV-2 Spike Protein P681H Variant Detected in Israel

The routine detection, surveillance, and reporting of novel SARS-CoV-2 variants is crucial, as these threaten to hinder global vaccination efforts. Herein we report a novel local variant with a non-synonymous mutation in the spike (S) protein P681H. This local Israeli variant was not associated with a higher infection rate or higher prevalence. Furthermore, the local variant was successfully neutralized by sera from fully vaccinated individuals at a comparable level to the B.1.1.7 variant and an Israel wild-type strain. While it is not a variant of concern, routine monitoring by sequencing is still required.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Diversity in the Circulation of Influenza A(H3N2) Viruses in the Northern Hemisphere in the 2018-19 Season

While vaccination is considered the most effective means to prevent influenza infection, its seasonal effectiveness varies, depending on the circulating influenza strains. Here, we characterized the circulation of influenza strains in October-2018 and March-2019 around the world. For this, we used nasopharyngeal samples collected from outpatient and hospitalized patients in Israel and data reported in ECDC, CDC, and WHO databases. Influenza A(H3N2) was dominant in Israel, while in Europe, Asia, and USA, A(H1N1)pdm09 virus circulated first, and then the A(H3N2) virus also appeared. Phylogenetic analysis indicated that A(H3N2) viruses circulating in Israel belonged to clade-3C.3a, while in Europe, Asia, and USA, A(H3N2) viruses belonged to subclade-3C.2a1, but were later replaced by clade-3C.3a viruses in USA. The vaccine A(H3N2) components of that year, A/Singapore/INFIMH-16-0019/2016-(H3N2)-like-viruses, belonged to clade-3C.2a1. The circulation of different influenza subtypes and clades of A(H3N2) viruses in a single season highlights the need for universal influenza vaccines.

Comprehensive Analyses of SARS-CoV-2 Transmission in a Public Health Virology Laboratory

SARS-CoV-2 has become a major global concern as of December 2019, particularly affecting healthcare workers. As person-to-person transmission is airborne, crowded closed spaces have high potential for rapid virus spread, especially early in the pandemic when social distancing and mask wearing were not mandatory. This retrospective study thoroughly investigates a small-scale SARS-CoV-2 outbreak in Israel’s central virology laboratory (ICVL) in mid-March 2020, in which six staff members and two related family members were infected. Suspicions regarding infection by contaminated surfaces in ICVL facilities were nullified by SARS-CoV-2 negative real time polymerase chain reaction (PCR) of work surfaces swipe tests. Complete SARS-CoV-2 genomes were sequenced and mutation analyses showed inclusion of all samples to clades 20B and 20C, possessing the spike mutation D614G. Phylogenetic analysis clarified transmission events, confirming S1 as having infected at least three other staff members and refuting the association of a staff member’s infected spouse with the ICVL transmission cluster. Finally, serology tests exhibited IgG and IgA antibodies in all infected individuals and revealed the occurrence of asymptomatic infections in additional staff members. This study demonstrates the advantages of molecular epidemiology in elucidating transmission events and exemplifies the importance of good laboratory practice, distancing and mask wearing in preventing SARS-CoV-2 spread, specifically in healthcare facilities.

The V(H) repertoire and clonal diversification of B cells in inflammatory myopathies

The contribution of antigen-driven B-cell adaptive immune responses within the inflamed muscle of inflammatory myopathies (IMs) is largely unknown. In this study, we investigated the immunoglobulin V(H) gene repertoire, somatic hypermutation, clonal diversification, and selection of infiltrating B cells in muscle biopsies from IM patients (dermatomyositis and polymyositis), to determine whether B cells and/or plasma cells contribute to the associated pathologies of these diseases. The data reveal that Ig V(H) gene repertoires of muscle-infiltrating B cells deviate from the normal V(H) gene repertoire in individual patients, and differ between different types of IMs. Analysis of somatic mutations revealed clonal diversification of muscle-infiltrating B cells and evidence for a chronic B-cell response within the inflamed muscle. We conclude that muscle-infiltrating B cells undergo selection, somatic hypermutation and clonal diversification in situ during antigen-driven immune responses in patients with IMs, providing insight into the contribution of B cells to the pathological mechanisms of these disorders.

PRC2/EED-EZH2 complex is up-regulated in breast cancer lymph node metastasis compared to primary tumor and correlates with tumor proliferation in situ

BACKGROUND

Lymph node metastasis is a key event in the progression of breast cancer. Therefore it is important to understand the underlying mechanisms which facilitate regional lymph node metastatic progression.

METHODOLOGY/PRINCIPAL FINDINGS

We performed gene expression profiling of purified tumor cells from human breast tumor and lymph node metastasis. By microarray network analysis, we found an increased expression of polycomb repression complex 2 (PRC2) core subunits EED and EZH2 in lymph node metastatic tumor cells over primary tumor cells which were validated through real-time PCR. Additionally, immunohistochemical (IHC) staining and quantitative image analysis of whole tissue sections showed a significant increase of EZH2 expressing tumor cells in lymph nodes over paired primary breast tumors, which strongly correlated with tumor cell proliferation in situ. We further explored the mechanisms of PRC2 gene up-regulation in metastatic tumor cells and found up-regulation of E2F genes, MYC targets and down-regulation of tumor suppressor gene E-cadherin targets in lymph node metastasis through GSEA analyses. Using IHC, the expression of potential EZH2 target, E-cadherin was examined in paired primary/lymph node samples and was found to be significantly decreased in lymph node metastases over paired primary tumors.

CONCLUSIONS/SIGNIFICANCE

This study identified an over expression of the epigenetic silencing complex PRC2/EED-EZH2 in breast cancer lymph node metastasis as compared to primary tumor and its positive association with tumor cell proliferation in situ. Concurrently, PRC2 target protein E-cadherin was significant decreased in lymph node metastases, suggesting PRC2 promotes epithelial mesenchymal transition (EMT) in lymph node metastatic process through repression of E-cadherin. These results indicate that epigenetic regulation mediated by PRC2 proteins may provide additional advantage for the outgrowth of metastatic tumor cells in lymph nodes. This opens up epigenetic drug development possibilities for the treatment and prevention of lymph node metastasis in breast cancer.

Altered local and systemic immune profiles underlie lymph node metastasis in breast cancer patients

Cancer-mediated immune dysfunction contributes to tumor progression and correlates with patient outcome. Metastasis to tumor draining lymph nodes (TDLNs) is an important step in breast cancer progression and is used to predict patient outcome and survival. Although lymph nodes are important immune organs, the role of immune cells in TDLNs has not been thoroughly investigated. We hypothesized that the host immune response in node negative (NN) patients is more intact and thereby can resist tumor invasion compared to node positive (NP) patients. As such, lymph node metastasis requires breakdown of the host immune response in addition to escape of cancer cells from the tumor. To investigate the immunological differences between NN and NP breast cancer patients, we purified and profiled immune cells from the three major compartments where cancer and immune cells interact: tumor, TDLNs and peripheral blood. Significant down-regulation of genes associated with immune-related pathways and up-regulation of genes associated with tumor-promoting pathways was consistently observed in NP patients' TDLNs compared to NN patients. Importantly, these signatures were seen even in NP patients' tumor-free TDLNs, suggesting that such immune changes are not driven solely by local tumor invasion. Furthermore, similar patterns were also observed in NP patients' tumor and blood immune cells, suggesting that immunological differences between NN and NP patients are systemic. Together, these findings suggest that alterations in overall immune function may underlie risk for LN metastasis in breast cancer patients.

Ig gene diversification and selection in follicular lymphoma, diffuse large B cell lymphoma and primary central nervous system lymphoma revealed by lineage tree and mutation analyses

Follicular lymphoma (FL), diffuse large B cell lymphoma (DLBCL) and primary central nervous system lymphoma are B cell malignancies. FL and DLBCL have a germinal center origin. We have applied mutational analyses and a novel algorithm for quantifying shape properties of mutational lineage trees to investigate the nature of the diversification, somatic hypermutation and selection processes that affect B cell clones in these malignancies and reveal whether they differ from normal responses. Lineage tree analysis demonstrated higher diversification and mutations per cell in the lymphoma clones. This was caused solely by the longer diversification times of the malignant clones, as their recent diversification processes were similar to those of normal responses, implying similar mutation frequencies. Since previous analyses of antigen-driven selection were shown to yield false positives, we performed a corrected analysis of replacement and silent mutation patterns, which revealed selection against replacement mutations in the framework regions, responsible for the structural integrity of the B cell receptor, but not for positive selection for replacements in the complementary determining regions. Most replacements, however, were neutral or conservative, suggesting that if at all selection operates in these malignancies it is for structural B cell receptor integrity but not for antigen binding.

Somatic hypermutation and antigen-driven selection of B cells are altered in autoimmune diseases

B cells have been found to play a critical role in the pathogenesis of several autoimmune (AI) diseases. A common feature amongst many AI diseases is the formation of ectopic germinal centers (GC) within the afflicted tissue or organ, in which activated B cells expand and undergo somatic hypermutation (SHM) and antigen-driven selection on their immunoglobulin variable region (IgV) genes. However, it is not yet clear whether these processes occurring in ectopic GCs are identical to those in normal GCs. The analysis of IgV mutations has aided in revealing many aspects concerning B cell expansion, mutation and selection in GC reactions. We have applied several mutation analysis methods, based on lineage tree construction, to a large set of data, containing IgV productive and non-productive heavy and light chain sequences from several different tissues, to examine three of the most profoundly studied AI diseases - Rheumatoid Arthritis (RA), Multiple Sclerosis (MS) and Sjögren's Syndrome (SS). We have found that RA and MS sequences exhibited normal mutation spectra and targeting motifs, but a stricter selection compared to normal controls, which was more apparent in RA. SS sequence analysis results deviated from normal controls in both mutation spectra and indications of selection, also showing differences between light and heavy chain IgV and between different tissues. The differences revealed between AI diseases and normal control mutation patterns may result from the different microenvironmental influences to which ectopic GCs are exposed, relative to those in normal secondary lymphoid tissues.

IgTree: creating Immunoglobulin variable region gene lineage trees

Lineage trees describe the microevolution of cells within an organism. They have been useful in the study of B cell affinity maturation, which is based on somatic hypermutation of immunoglobulin genes in germinal centers and selection of the resulting mutants. Our aim was to create and implement an algorithm that can generate lineage trees from immunoglobulin variable region gene sequences. The IgTree program implements the algorithm we developed, and generates lineage trees. Original sequences found in experiments are assigned to either leaves or internal nodes of the tree. Each tree node represents a single mutation separating the sequences. The mutations that separate the sequences from each other can be point mutations, deletions or insertions. The program can deal with gaps and find potential reversion mutations. The program also enumerates mutation frequencies and sequence motifs around each mutation, on a per-tree basis. The algorithm has proven useful in several studies of immunoglobulin variable region gene mutations.

Antigen-driven selection in germinal centers as reflected by the shape characteristics of immunoglobulin gene lineage trees: a large-scale simulation study

During the immune response, the generation of memory B lymphocytes in germinal centers involves affinity maturation of the cells' antigen receptors, based on somatic hypermutation of receptor genes and antigen-driven selection of the resulting mutants. Affinity maturation is vital for immune protection, and is the basis of humoral immune learning and memory. Lineage trees of somatically hypermutated immunoglobulin genes often serve to qualitatively illustrate claims concerning the dynamics of affinity maturation in germinal centers. Here, we derive the quantitative relationships between parameters characterizing affinity maturation dynamics (proliferation, differentiation and mutation rates, initial affinity of the Ig to the antigen, and selection thresholds) and the mathematical properties of lineage trees, using a computer simulation which combines mathematical models for all mature B cell populations, stochastic models of hypermutation and selection, lineage tree generation and measurement of graphical tree characteristics. We identified seven key lineage tree properties, and found correlations of these with initial clone affinity and with the selection threshold. These two parameters were found to be the main factors affecting lineage tree shapes in both primary and secondary response trees. The results also confirm that recycling from centrocytes back to centroblasts is highly likely.

Lineage tree analysis of immunoglobulin variable-region gene mutations in autoimmune diseases: chronic activation, normal selection

Autoimmune diseases show high diversity in the affected organs, clinical manifestations and disease dynamics. Yet they all share common features, such as the ectopic germinal centers found in many affected tissues. Lineage trees depict the diversification, via somatic hypermutation (SHM), of immunoglobulin variable-region (IGV) genes. We previously developed an algorithm for quantifying the graphical properties of IGV gene lineage trees, allowing evaluation of the dynamical interplay between SHM and antigen-driven selection in different lymphoid tissues, species, and disease situations. Here, we apply this method to ectopic GC B cell clones from patients with Myasthenia Gravis, Rheumatoid Arthritis, and Sjögren's Syndrome, using data scaling to minimize the effects of the large variability due to methodological differences between groups. Autoimmune trees were found to be significantly larger relative to normal controls. In contrast, comparison of the measurements for tree branching indicated that similar selection pressure operates on autoimmune and normal control clones.

Novel Analysis of Clonal Diversification in Blood B Cell and Bone Marrow Plasma Cell Clones in Immunoglobulin Light Chain Amyloidosis

Immunoglobulin light chain amyloidosis (AL) is characterized by a limited clonal expansion of plasma cells and amyloid formation. Here, we report restriction in the diversity of VL gene usage with a dominance of clonally related B cells in the peripheral blood (PB) isotype-specific repertoire of AL patients. A rigorous quantification of lineage trees reveals presence of intraclonal variations in the PB clones compared to the bone marrow (BM) clones, which suggests a common precursor that is still subject to somatic mutation. When compared to normal BM and PB B cells, AL clones showed significant but incomplete impairment of antigenic selection, which could not be detected by conventional R and S mutation analysis. Therefore, graphical analysis of B cell lineage trees and mathematical quantification of tree properties provide novel insights into the process of B cell clonal evolution in AL.

Quantitative analysis of clonal bone marrow CD19+ B cells: use of B cell lineage trees to delineate their role in the pathogenesis of light chain amyloidosis

Light chain amyloidosis (AL) is a bone marrow (BM) plasma cell neoplasia with systemic deposition of Ig light chain amyloid fibrils. Here, we report the identification of clonal CD19 B cells in the BM and the use of a novel mathematical algorithm to generate B cell lineage trees of the clonal CD19 B cells and CD138 plasma cells from the BM of AL patients to delineate the relationship between these two clonal populations. The CD19+ clonal B cells in the BM of AL patients related to the clonal plasma cells represent a pre-plasma cell precursor population. The B cell lineage trees from AL patients also show significant differences in clonal diversification and antigenic selection compared to clones from normal, healthy controls. These data provide a robust example of the use of graphical quantification methods in delineating the role of neoplastic precursors in the pathogenesis of hematopoietic malignancies.