Insights from genomic analysis of a novel Coxiella burnetii strain isolated in Israel

The diagnosis of Q fever is challenging due to nonspecific symptoms and negative standard blood culture results. Serological testing through immunofluorescence assay (IFA) is the most commonly used method for diagnosing this disease. Polymerase chain reaction (PCR) tests can also be used to detect bacterial DNA if taken at an appropriate time. Once the presence of bacteria is confirmed in a sample, an enrichment step is required before characterizing it through sequencing. Cultivating C. burnetii is challenging as it can only be isolated by inoculation into cell culture, embryonated eggs, or animals. In this article, we describe the isolation of C. burnetii from a valve specimen in Vero cells. We conducted genome sequencing and taxonomy profiling of this isolate and were able to determine its taxonomic affiliation. Furthermore, Multispacer sequence typing (MST) analysis suggests that the infection originated from a local strain of C. burnetii found around northern Israel and Lebanon. This novel strain belongs to a previously described genotype MST6, harboring the QpRS plasmid, never reported in Israel.

Enhanced production yields of rVSV-SARS-CoV-2 vaccine using Fibra-Cel® macrocarriers

The COVID-19 pandemic has led to high global demand for vaccines to safeguard public health. To that end, our institute has developed a recombinant viral vector vaccine utilizing a modified vesicular stomatitis virus (VSV) construct, wherein the G protein of VSV is replaced with the spike protein of SARS-CoV-2 (rVSV-ΔG-spike). Previous studies have demonstrated the production of a VSV-based vaccine in Vero cells adsorbed on Cytodex 1 microcarriers or in suspension. However, the titers were limited by both the carrier surface area and shear forces. Here, we describe the development of a bioprocess for rVSV-ΔG-spike production in serum-free Vero cells using porous Fibra-Cel® macrocarriers in fixed-bed BioBLU®320 5p bioreactors, leading to high-end titers. We identified core factors that significantly improved virus production, such as the kinetics of virus production, the use of macrospargers for oxygen supply, and medium replenishment. Implementing these parameters, among others, in a series of GMP production processes improved the titer yields by at least two orders of magnitude (2e9 PFU/mL) over previously reported values. The developed process was highly effective, repeatable, and robust, creating potent and genetically stable vaccine viruses and introducing new opportunities for application in other viral vaccine platforms.

VSV-ΔG-Spike Candidate Vaccine Induces Protective Immunity and Protects K18-hACE2 Mice against SARS-CoV-2 Variants

Since the emergence of the original SARS-CoV-2, several variants were described, raising questions as to the ability of recently developed vaccine platforms to induce immunity and provide protection against these variants. Here, we utilized the K18-hACE2 mouse model to show that VSV-ΔG-spike vaccination provides protection against several SARS-CoV-2 variants: alpha, beta, gamma, and delta. We show an overall robust immune response, regardless of variant identity, leading to reduction in viral load in target organs, prevention of morbidity and mortality, as well as prevention of severe brain immune response, which follows infection with various variants. Additionally, we provide a comprehensive comparison of the brain transcriptomic profile in response to infection with different variants of SARS-CoV-2 and show how vaccination prevents these disease manifestations. Taken together, these results highlight the robust VSV-ΔG-spike protective response against diverse SARS-CoV-2 variants, as well as its promising potential against newly arising variants.

Global transcriptomic analysis of Francisella tularensis SchuS4 differentially expressed genes in response to doxycycline or ciprofloxacin exposure

OBJECTIVE

As part of a research aiming at presenting an alternative approach for rapid determination of antimicrobial susceptibility by quantification of changes in expression levels of specific marker genes and gene sets, cultures of the virulent bacterial strain Francisella tularensis SchuS4 were grown in the presence of inhibitory/sub-inhibitory concentrations of either ciprofloxacin or doxycycline and their transcriptomic profiles were elucidated using differential expression analysis followed by functional annotation.

DATA DESCRIPTION

RNA sequencing was performed to identify differentially expressed genes (DEGs) in response to exposure of F. tularensis SchuS4 to either ciprofloxacin or doxycycline, the antibiotics of choice for Tularemia therapy. Accordingly, RNA samples were collected 2 h post antibiotic exposure and subjected to RNA sequence analysis. Transcriptomic quantification of RNA representing duplicated samples generated highly similar gene expression data. Exposure to sub-inhibitory concentration [0.5 x MIC (minimal inhibitory concentration)] of doxycycline or ciprofloxacin modulated the expression of 237 or 8 genes, respectively, while exposure to an inhibitory concentration (1 x MIC) resulted in the modulation of 583 or 234 genes, respectively. Amongst the genes modulated upon doxycycline exposure upregulation of 31 genes encoding for translation-functions could be distinguished, as well as downregulation of 14 genes encoding for functions involved in DNA transcription and repair. Ciprofloxacin exposure impacted differently the RNA sequence profile of the pathogen, resulting in upregulation of 27 genes encoding mainly DNA replication and repair functions, transmembrane transporters and molecular chaperons. In addition, 15 downregulated genes were involved in translation processes.

Genome-wide CRISPR screens identify GATA6 as a proviral host factor for SARS-CoV-2 via modulation of ACE2

The global spread of SARS-CoV-2 led to major economic and health challenges worldwide. Revealing host genes essential for infection by multiple variants of SARS-CoV-2 can provide insights into the virus pathogenesis, and facilitate the development of novel therapeutics. Here, employing a genome-scale CRISPR screen, we provide a comprehensive data-set of cellular factors that are exploited by wild type SARS-CoV-2 as well as two additional recently emerged variants of concerns (VOCs), Alpha and Beta. We identified several host factors critical for SARS-CoV-2 infection, including various components belonging to the Clathrin-dependent transport pathway, ubiquitination, Heparan sulfate biogenesis and host phosphatidylglycerol biosynthesis. Comparative analysis of the different VOCs revealed the host factors KREMEN2 and SETDB1 as potential unique candidates required only to the Alpha variant. Furthermore, the analysis identified GATA6, a zinc finger transcription factor, as an essential proviral gene for all variants inspected. We show that GATA6 directly regulates ACE2 transcription and accordingly, is critical for SARS-CoV-2 cell entry. Analysis of clinical samples collected from SARS-CoV-2 infected individuals shows elevated levels of GATA6, suggesting a role in COVID-19 pathogenesis. Finally, pharmacological inhibition of GATA6 resulted in down-modulation of ACE2 and inhibition of viral infectivity. Overall, we show GATA6 may represent a target for the development of anti-SARS-CoV-2 therapeutic strategies and reaffirm the value of the CRISPR loss-of-function screens in providing a list of potential new targets for therapeutic interventions.

Whole genome sequencing and taxonomic profiling of two Pantoea sp. isolated from environmental samples in Israel

OBJECTIVE

As part of a research aiming at the isolation of bacteria secreting growth inhibiting compounds, cultures of Francisella tularensis were implanted in environmental samples and monitored for inhibition zones on agar. Two antibiotic-like secreting bacteria were isolated, their genomic sequence was deciphered and taxonomic profiling analysis classified them as belonging to the Pantoea genus.

DATA DESCRIPTION

Two bacterial isolates exhibiting growth inhibition zones to F. tularensis (LVS) were analyzed using the Oxford Nanopore Technology (ONT). Preliminary de novo assembly of the reads was performed, followed by taxonomic profiling based on Multi Locus Sequence Analysis (MLSA) and implementation of the Average Nucleotide Identity (ANI) measure. The genomic sequences resulted in the identification of two different Pantoea species, denoted EnvD and EnvH. Subsequent de novo genome assembly generated 5 and 10 contigs for EnvD and EnvH, respectively. The largest contig (4,008,183 bps and 3,740,753 bps for EnvD and EnvH, respectively), overlaps to a major extent to the chromosome of closely related Pantoea species. ANI values calculated for both isolates revealed two apparently new species of the Pantoea genus. Our study deciphered the identity of two bacteria producing antibiotic-like compounds, and the genomic sequence revealed they represent distinct Pantoea species.

Exportation of Monkeypox Virus From the African Continent

BACKGROUND

The largest West African monkeypox outbreak began September 2017, in Nigeria. Four individuals traveling from Nigeria to the United Kingdom (n = 2), Israel (n = 1), and Singapore (n = 1) became the first human monkeypox cases exported from Africa, and a related nosocomial transmission event in the United Kingdom became the first confirmed human-to-human monkeypox transmission event outside of Africa.

METHODS

Epidemiological and molecular data for exported and Nigerian cases were analyzed jointly to better understand the exportations in the temporal and geographic context of the outbreak.

RESULTS

Isolates from all travelers and a Bayelsa case shared a most recent common ancestor and traveled to Bayelsa, Delta, or Rivers states. Genetic variation for this cluster was lower than would be expected from a random sampling of genomes from this outbreak, but data did not support direct links between travelers.

CONCLUSIONS

Monophyly of exportation cases and the Bayelsa sample, along with the intermediate levels of genetic variation, suggest a small pool of related isolates is the likely source for the exported infections. This may be the result of the level of genetic variation present in monkeypox isolates circulating within the contiguous region of Bayelsa, Delta, and Rivers states, or another more restricted, yet unidentified source pool.

Neutralization of SARS-CoV-2 Variants by rVSV-ΔG-Spike-Elicited Human Sera

The emergence of rapidly spreading variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a major challenge to the ability of vaccines and therapeutic antibodies to provide immunity. These variants contain mutations of specific amino acids that might impede vaccine efficacy. BriLife^® (rVSV-ΔG-spike) is a newly developed SARS-CoV-2 vaccine candidate currently in phase II clinical trials. It is based on a replication-competent vesicular stomatitis virus (VSV) platform. The rVSV-ΔG-spike contains several spontaneously acquired spike mutations that correspond to SARS-CoV-2 variants’ mutations. We show that human sera from BriLife^® vaccinees preserve comparable neutralization titers towards alpha, gamma, and delta variants and show less than a three-fold reduction in the neutralization capacity of beta and omicron compared to the original virus. Taken together, we show that human sera from BriLife^® vaccinees overall maintain a neutralizing antibody response against all tested variants. We suggest that BriLife^®-acquired mutations may prove advantageous against future SARS-CoV-2 VOCs.

Mice with induced pulmonary morbidities display severe lung inflammation and mortality following exposure to SARS-CoV-2

Mice are normally unaffected by SARS coronavirus 2 (SARS-CoV-2) infection since the virus does not bind effectively to the murine version of the angiotensin-converting enzyme 2 (ACE2) receptor molecule. Here, we report that induced mild pulmonary morbidities rendered SARS-CoV-2–refractive CD-1 mice susceptible to this virus. Specifically, SARS-CoV-2 infection after application of low doses of the acute lung injury stimulants bleomycin or ricin caused severe disease in CD-1 mice, manifested by sustained body weight loss and mortality rates greater than 50%. Further studies revealed markedly higher levels of viral RNA in the lungs, heart, and serum of low-dose ricin–pretreated mice compared with non-pretreated mice. Furthermore, lung extracts prepared 2–3 days after viral infection contained subgenomic mRNA and virus particles capable of replication only when derived from the pretreated mice. The deleterious effects of SARS-CoV-2 infection were effectively alleviated by passive transfer of polyclonal or monoclonal antibodies generated against the SARS-CoV-2 receptor binding domain (RBD). Thus, viral cell entry in the sensitized mice seems to depend on viral RBD binding, albeit by a mechanism other than the canonical ACE2-mediated uptake route. This unique mode of viral entry, observed over a mildly injured tissue background, may contribute to the exacerbation of coronavirus disease 2019 (COVID-19) pathologies in patients with preexisting morbidities.

Mutation Profile of SARS-CoV-2 Genome Sequences Originating from Eight Israeli Patient Isolates

We report the genome sequences and the identification of genetic variations in eight clinical samples of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Samples were collected from nasopharyngeal swabs of symptomatic and asymptomatic individuals from five care homes for elderly and infirm persons in Israel. The sequences obtained are valuable, as they carry a newly reported nonsynonymous substitution located within the nucleoprotein open reading frame.

Comparative Analysis of the Global Transcriptomic Response to Oxidative Stress of Bacillus anthracis htrA-Disrupted and Parental Wild Type Strains

We previously demonstrated that the HtrA (High Temperature Requirement A) protease/chaperone active in the quality control of protein synthesis, represents an important virulence determinant of Bacillus anthracis. Virulence attenuation of htrA-disrupted Bacillus anthracis strains was attributed to susceptibility of ΔhtrA strains to stress insults, as evidenced by affected growth under various stress conditions. Here, we report a comparative RNA-seq transcriptomic study generating a database of differentially expressed genes in the B. anthracishtrA-disrupted and wild type parental strains under oxidative stress. The study demonstrates that, apart from protease and chaperone activities, HtrA exerts a regulatory role influencing expression of more than 1000 genes under stress. Functional analysis of groups or individual genes exhibiting strain-specific modulation, evidenced (i) massive downregulation in the ΔhtrA and upregulation in the WT strains of various transcriptional regulators, (ii) downregulation of translation processes in the WT strain, and (iii) downregulation of metal ion binding functions and upregulation of sporulation-associated functions in the ΔhtrA strain. These modulated functions are extensively discussed. Fifteen genes uniquely upregulated in the wild type strain were further interrogated for their modulation in response to other stress regimens. Overexpression of one of these genes, encoding for MazG (a nucleoside triphosphate pyrophosphohydrolase involved in various stress responses in other bacteria), in the ΔhtrA strain resulted in partial alleviation of the H₂O₂-sensitive phenotype.

Coding-Complete Genome Sequences of Two SARS-CoV-2 Isolates from Early Manifestations of COVID-19 in Israel

We announce the genome sequences of two strains of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) isolated in Israel, one imported by a traveler who returned from Japan and the second strain collected from a patient infected by a traveler returning from Italy. The sequences obtained are valuable as early manifestations for future follow-up of the local spread of the virus in Israel.

A rapid high-throughput sequencing-based approach for the identification of unknown bacterial pathogens in whole blood

High-throughput DNA sequencing (HTS) of pathogens in whole blood samples is hampered by the high host/pathogen nucleic acids ratio. We describe a novel and rapid bacterial enrichment procedure whose implementation is exemplified in simulated bacteremic human blood samples. The procedure involves depletion of the host DNA, rapid HTS and bioinformatic analyses. Following this procedure, Y. pestis, F. tularensis and B. anthracis spiked-in samples displayed an improved host/pathogen DNA ratio of 2.5-5.9 orders of magnitude, in samples with bacteria spiked-in at 103-105 CFU/ml. The procedure described in this study enables rapid and detailed metagenomic profiling of pathogens within 8-9 h, circumventing the challenges imposed by the high background present in the bacteremic blood and by the unknown nature of the sample.

Case Report: Imported Melioidosis from Goa, India to Israel, 2018

A previously healthy young man presented with a chronic cavitary pulmonary infection that began while in Goa, India. Burkholderia pseudomallei was cultured from sputum samples. The infection fully resolved after prolonged antibiotic treatment. Other than traveling during the monsoon season, extensive use of well-water for water-pipe smoking of cannabis was identified as a possible risk factor for infection. This is one of the first reports of travel-associated melioidosis from India. Genomic and immunological characterization suggested that the B. pseudomallei isolate collected from the reported case exhibited limited similarity to other B. pseudomallei strains.

Rapid identification of unknown pathogens in environmental samples using a high-throughput sequencing-based approach

In the event of a bioterror attack, a prompt, sensitive and definite identification of the agents involved is of major concern for confirmation of the event and for mitigation of countermeasures. Whether the information from intelligence forces is limited concerning the biothreat identity or one suspects the presence of a novel or engineered agent, the genetic identification of microorganisms in an unknown sample is challenging. High-throughput sequencing (HTS) technologies can sequence a heterogeneous mixture of genetic materials with high sensitivity and speed; nevertheless, despite the enormous advantages of HTS, all previous reports have analyzed unknown samples in a timeframe of a few days to a few weeks. This timeframe might not be relevant to an emergency scenario. Here, we present an HTS-based approach for deciphering the genetic composition of unknown samples within a working day. This outcome is accomplished by a rapid library preparation procedure, short-length sequencing and a prompt bioinformatics comparison against all available microbial genomic sequences. Using this approach, as a proof of concept, we were able to detect two spiked-in biothreat agents, B. anthracis and Y. pestis, in a variety of environmental samples at relevant concentrations and within a short timeframe of eight hours.

Processes of fungal proteome evolution and gain of function: gene duplication and domain rearrangement

During evolution, organisms have gained functional complexity mainly by modifying and improving existing functioning systems rather than creating new ones ab initio. Here we explore the interplay between two processes which during evolution have had major roles in the acquisition of new functions: gene duplication and protein domain rearrangements. We consider four possible evolutionary scenarios: gene families that have undergone none of these event types; only gene duplication; only domain rearrangement, or both events. We characterize each of the four evolutionary scenarios by functional attributes. Our analysis of ten fungal genomes indicates that at least for the fungi clade, species significantly appear to gain complexity by gene duplication accompanied by the expansion of existing domain architectures via rearrangements. We show that paralogs gaining new domain architectures via duplication tend to adopt new functions compared to paralogs that preserve their domain architectures. We conclude that evolution of protein families through gene duplication and domain rearrangement is correlated with their functional properties. We suggest that in general, new functions are acquired via the integration of gene duplication and domain rearrangements rather than each process acting independently.

Evolution of domain promiscuity in eukaryotic genomes--a perspective from the inferred ancestral domain architectures

Most eukaryotic proteins are composed of two or more domains. These assemble in a modular manner to create new proteins usually by the acquisition of one or more domains to an existing protein. Promiscuous domains which are found embedded in a variety of proteins and co-exist with many other domains are of particular interest and were shown to have roles in signaling pathways and mediating network communication. The evolution of domain promiscuity is still an open problem, mostly due to the lack of sequenced ancestral genomes. Here we use inferred domain architectures of ancestral genomes to trace the evolution of domain promiscuity in eukaryotic genomes. We find an increase in average promiscuity along many branches of the eukaryotic tree. Moreover, domain promiscuity can proceed at almost a steady rate over long evolutionary time or exhibit lineage-specific acceleration. We also observe that many signaling and regulatory domains gained domain promiscuity around the Bilateria divergence. In addition we show that those domains that played a role in the creation of two body axes and existed before the divergence of the bilaterians from fungi/metazoan achieve a boost in their promiscuities during the bilaterian evolution.

Modular genes with metazoan-specific domains have increased tissue specificity

We have systematically examined the domain composition across a comprehensive set of tissue-specific, midrange and housekeeping genes as defined by their mode of expression in 52 normal mouse tissues. We show a definite correlation between the number of domains and the degree of tissue specificity. This trend is further supported by a novel analysis involving the time of origin of each domain. Genes containing metazoan-specific domains are more prevalent in signal transduction and cell-communication pathways, and are depleted in primary metabolism. Our analyses suggest that highly modular gene products have been recruited for tissue-specific functions that are required in complex organisms.