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.

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.

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.

Single-Cell RNA Sequencing Reveals mRNA Splice Isoform Switching during Kidney Development.

BACKGROUND: During mammalian kidney development, nephron progenitors undergo a mesenchymal-to-epithelial transition and eventually differentiate into the various tubular segments of the nephron. Recently, Drop-seq single-cell RNA sequencing technology for measuring gene expression from thousands of individual cells identified the different cell types in the developing kidney. However, that analysis did not include the additional layer of heterogeneity that alternative mRNA splicing creates. METHODS: Full transcript length single-cell RNA sequencing characterized the transcriptomes of 544 individual cells from mouse embryonic kidneys. RESULTS: Gene expression levels measured with full transcript length single-cell RNA sequencing identified each cell type. Further analysis comprehensively characterized splice isoform switching during the transition between mesenchymal and epithelial cellular states, which is a key transitional process in kidney development. The study also identified several putative splicing regulators, including the genes Esrp1/2 and Rbfox1/2. CONCLUSIONS: Discovery of the sets of genes that are alternatively spliced as the fetal kidney mesenchyme differentiates into tubular epithelium will improve our understanding of the molecular mechanisms that drive kidney development.

The SARS-CoV-2 Lambda variant and its neutralisation efficiency following vaccination with Comirnaty, Israel, April to June 2021.

The SARS-CoV-2 Lambda (Pango lineage designation C.37) variant of interest, initially identified in Peru, has spread to additional countries. First detected in Israel in April 2021 following importations from Argentina and several European countries, the Lambda variant infected 18 individuals belonging to two main transmission chains without further spread. Micro-neutralisation assays following Comirnaty (BNT162b2 mRNA, BioNTech-Pfizer) vaccination demonstrated a significant 1.6-fold reduction in neutralising titres compared with the wild type virus, suggesting increased susceptibility of vaccinated individuals to infection.

Characterization of hepatitis B and delta coinfection in Israel.

BACKGROUND: Characteristics of hepatitis B (HBV) and delta (HDV) coinfection in various geographical regions, including Israel, remain unclear. Here we studied HDV seroprevalence in Israel, assessed HDV/HBV viral loads, circulating genotypes and hepatitis delta antigen (HDAg) conservation. METHODS: Serological anti HDV IgG results from 8969 HBsAg positive individuals tested in 2010-2015 were retrospectively analyzed to determine HDV seroprevalence. In a cohort of HBV/HDV coinfected (n=58) and HBV monoinfected (n=27) patients, quantitative real-time PCR (qRT-PCR) and sequencing were performed to determine viral loads, genotypes and hepatitis delta antigen (HDAg) protein sequence. RESULTS: 6.5% (587/8969) of the HBsAg positive patients were positive for anti HDV antibodies. HDV viral load was >2 log copies/ml higher than HBV viral load in most of the coinfected patients with detectable HDV RNA (86%, 50/58). HDV genotype 1 was identified in all patients, most of whom did not express HBV. While 66.6% (4/6) of the HBV/HDV co-expressing patients carried HBV-D2 only 18.5% (5/27) of the HBV monoinfections had HBV-D2 (p=0.03). Higher genetic variability in the HDAg protein sequence was associated with higher HDV viral load. CONCLUSIONS: The overall significant prevalence of HDV (6.5%) mandates HDV RNA testing for all coinfected patients. Patients positive for HDV RNA (characterized by low HBV DNA blood levels) carried HDV genotype 1. Taken together, the significant HDV seroprevalence and the lack of effective anti-HDV therapy, necessitates strict clinical surveillance especially in patients with higher HDV viral loads and increased viral evolution.

Measles outbreak associated with a preschool setting among partially vaccinated children in the Tel Aviv District, Israel, October 2023.

In October 2023, the Tel Aviv District was notified of ten cases of measles. The outbreak initiated in a preschool with high vaccination coverage with one dose of MMR vaccine. Serological testing was available for eight patients (six children and two adults). Among the six children vaccinated with one dose of MMR vaccine, primary vaccine failure was demonstrated. Among the adults, secondary vaccine failure was confirmed. The outbreak was successfully contained due to a combination of factors, notably its occurrence within a population characterized by high vaccination coverage in Tel Aviv, during a period of restricted public interactions due to the prevailing state of war in the country. Despite challenging wartime conditions, effective prophylactic measures were promptly executed, encompassing a 2-dose MMR vaccination schedule for close contacts and the broader community of children in the TA district, successfully curbing the outbreak and preventing widespread infections.

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.

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.

Human metapneumovirus prevalence during 2019-2021 in Israel is influenced by the COVID-19 pandemic.

OBJECTIVES: To compare infection rates and circulating subtypes of human metapneumovirus (hMPV) before (2019-2020) and after the emergence of coronavirus disease 2019 (COVID-19) (2021) in Israel. METHODS: In total, 12,718 respiratory samples were collected from hospitalized patients of all ages during the years 2019 to 2021 at the Sheba Medical Center in Israel and subjected to reverse transcription-polymerase chain reaction analysis. In addition, whole-genome sequencing was performed to characterize the subtypes of hMPV circulating in Israel between 2019 and 2021. RESULTS: A total of 481 samples were found positive for hMPV. Before the emergence of COVID-19, hMPV peaked in winter months and declined thereafter. In sharp contrast, during the COVID-19 pandemic, we observed a delayed peak in hMPV infection cases and higher infection of young children. Viral sequencing showed a shift in the most prevalent circulating hMPV strain from A2b to B1 during the years 2019, 2020, and 2021. CONCLUSION: Compared with the years before the COVID-19 pandemic, in 2021, hMPV mostly affected young children, and the most prevalent circulating subtype shifted from A2b in 2019 to B1.

Characterization of alternative mRNA splicing in cultured cell populations representing progressive stages of human fetal kidney development.

Nephrons are the functional units of the kidney. During kidney development, cells from the cap mesenchyme-a transient kidney-specific progenitor state-undergo a mesenchymal to epithelial transition (MET) and subsequently differentiate into the various epithelial cell types that create the tubular structures of the nephron. Faults in this transition can lead to a pediatric malignancy of the kidney called Wilms' tumor that mimics normal kidney development. While human kidney development has been characterized at the gene expression level, a comprehensive characterization of alternative splicing is lacking. Therefore, in this study, we performed RNA sequencing on cell populations representing early, intermediate, and late developmental stages of the human fetal kidney, as well as three blastemal-predominant Wilms' tumor patient-derived xenografts. Using this newly generated RNAseq data, we identified a set of transcripts that are alternatively spliced between the different developmental stages. Moreover, we found that cells from the earliest developmental stage have a mesenchymal splice-isoform profile that is similar to that of blastemal-predominant Wilms' tumor xenografts. RNA binding motif enrichment analysis suggests that the mRNA binding proteins ESRP1, ESRP2, RBFOX2, and QKI regulate alternative mRNA splicing during human kidney development. These findings illuminate new molecular mechanisms involved in human kidney development and pediatric kidney cancer.

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.

Specific Detection of SARS-CoV-2 Variants B.1.1.7 (Alpha) and B.1.617.2 (Delta) Using a One-Step Quantitative PCR Assay.

In this report, we describe the development of a reverse transcription-quantitative PCR (RT-qPCR) assay, termed Alpha-Delta assay, which can detect all severe acute respiratory syndrome coronavirus 2 (SC-2) variants and distinguish between the Alpha (B.1.1.7) and Delta (B.1.617.2) variants. The Alpha- and Delta-specific reactions in the assay target mutations that are strongly linked to the target variant. The Alpha reaction targets the D3L substitution in the N gene, and the Delta reaction targets the spike gene 156 to 158 mutations. Additionally, we describe a second Delta-specific assay that we use as a confirmatory test for the Alpha-Delta assay that targets the 119 to 120 deletion in the Orf8 gene. Both reactions have similar sensitivities of 15 to 25 copies per reaction, similar to the sensitivity of commercial SC-2 detection tests. The Alpha-Delta assay and the Orf8119del assay were successfully used to classify clinical samples that were subsequently analyzed by whole-genome sequencing. Lastly, the capability of the Alpha-Delta assay and Orf8119del assay to identify correctly the presence of Delta RNA in wastewater samples was demonstrated. This study provides a rapid, sensitive, and cost-effective tool for detecting and classifying two worldwide dominant SC-2 variants. It also highlights the importance of a timely diagnostic response to the emergence of new SC-2 variants with significant consequences on global health. IMPORTANCE The new assays described herein enable rapid, straightforward, and cost-effective detection of severe acute respiratory syndrome coronavirus 2 (SC-2) with immediate classification of the examined sample as Alpha, Delta, non-Alpha, or non-Delta variant. This is highly important for two main reasons: (i) it provides the scientific and medical community with a novel diagnostic tool to rapidly detect and classify any SC-2 sample of interest as Alpha, Delta, or none and can be applied to both clinical and environmental samples, and (ii) it demonstrates how to respond to the emergence of new variants of concern by developing a variant-specific assay. Such assays should improve our preparedness and adjust the diagnostic capacity to serve clinical, epidemiological, and research needs.

Direct sequencing of measles virus complete genomes in the midst of a large-scale outbreak.

Measles outbreaks escalated globally despite worldwide elimination efforts. Molecular epidemiological investigations utilizing partial measles virus (MeV) genomes are challenged by reduction in global genotypes and low evolutionary rates. Greater resolution was reached using MeV complete genomes, however time and costs limit the application to numerous samples. We developed an approach to unbiasedly sequence complete MeV genomes directly from patient urine samples. Samples were enriched for MeV using filtration or nucleases and the minimal number of sequence reads to allocate per sample based on its MeV content was assessed using in-silico reduction of sequencing depth. Application of limited-resource sequencing to treated MeV-positive samples demonstrated that 1-5 million sequences for samples with high/medium MeV quantities and 10-15 million sequences for samples with lower MeV quantities are sufficient to obtain >98% MeV genome coverage and over X50 average depth. This approach enables real-time high-resolution molecular epidemiological investigations of large-scale MeV outbreaks.

Genomic variation and epidemiology of SARS-CoV-2 importation and early circulation in Israel.

Severe acute respiratory disease coronavirus 2 (SARS-CoV-2) which causes corona virus disease (COVID-19) was first identified in Wuhan, China in December 2019 and has since led to a global pandemic. Importations of SARS-CoV-2 to Israel in late February from multiple countries initiated a rapid outbreak across the country. In this study, SARS-CoV-2 whole genomes were sequenced from 59 imported samples with a recorded country of importation and 101 early circulating samples in February to mid-March 2020 and analyzed to infer clades and mutational patterns with additional sequences identified Israel available in public databases. Recorded importations in February to mid-March, mostly from Europe, led to multiple transmissions in all districts in Israel. Although all SARS-CoV-2 defined clades were imported, clade 20C became the dominating clade in the circulating samples. Identification of novel, frequently altered mutated positions correlating with clade-defining positions provide data for surveillance of this evolving pandemic and spread of specific clades of this virus. SARS-CoV-2 continues to spread and mutate in Israel and across the globe. With economy and travel resuming, surveillance of clades and accumulating mutations is crucial for understanding its evolution and spread patterns and may aid in decision making concerning public health issues.

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.

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.

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.

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.

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.