The BNT162b2 mRNA vaccine is highly effective against SARS-CoV-2. However, apprehension exists that variants of concern (VOCs) may evade vaccine protection, due to evidence of reduced neutralization of the VOCs B.1.1.7 and B.1.351 by vaccine sera in laboratory assays. We performed a matched cohort study to examine the distribution of VOCs in infections of BNT162b2 mRNA vaccinees from Clalit Health Services (Israel) using viral genomic sequencing, and hypothesized that if vaccine effectiveness against a VOC is reduced, its proportion among breakthrough cases would be higher than in unvaccinated controls. Analyzing 813 viral genome sequences from nasopharyngeal swabs, we showed that vaccinees who tested positive at least 7 days after the second dose were disproportionally infected with B.1.351, compared with controls. Those who tested positive between 2 weeks after the first dose and 6 days after the second dose were disproportionally infected by B.1.1.7. These findings suggest reduced vaccine effectiveness against both VOCs within particular time windows. Our results emphasize the importance of rigorously tracking viral variants, and of increasing vaccination to prevent the spread of VOCs.

Full genome sequences are increasingly used to track the geographic spread and transmission dynamics of viral pathogens. Here, with a focus on Israel, we sequenced 212 SARS-CoV-2 sequences and use them to perform a comprehensive analysis to trace the origins and spread of the virus. A phylogenetic analysis including thousands of globally sampled sequences allowed us to infer multiple independent introductions into Israel, followed by local transmission. Returning travelers from the U.S. contributed dramatically more to viral spread relative to their proportion in incoming infected travelers. Using phylodynamic analysis, we estimated that the basic reproduction number of the virus was initially around ~2.0-2.6, dropping by two-thirds following the implementation of social distancing measures. A comparison between reported and model-estimated case numbers indicated high levels of transmission heterogeneity in SARS-CoV-2 spread, with between 1-10% of infected individuals resulting in 80% of secondary infections. Overall, our findings underscore the ability of this virus to efficiently transmit between and within countries, as well as demonstrate the effectiveness of social distancing measures for reducing its spread.

Full genome sequences are increasingly used to track the geographic spread and transmission dynamics of viral pathogens. Here, with a focus on Israel, we sequence 212 SARS-CoV-2 sequences and use them to perform a comprehensive analysis to trace the origins and spread of the virus. We find that travelers returning from the United States of America significantly contributed to viral spread in Israel, more than their proportion in incoming infected travelers. Using phylodynamic analysis, we estimate that the basic reproduction number of the virus was initially around 2.5, dropping by more than two-thirds following the implementation of social distancing measures. We further report high levels of transmission heterogeneity in SARS-CoV-2 spread, with between 2-10% of infected individuals resulting in 80% of secondary infections. Overall, our findings demonstrate the effectiveness of social distancing measures for reducing viral spread.

The SARS-CoV-2 pandemic has been raging for over a year, creating global detrimental impact. The BNT162b2 mRNA vaccine has demonstrated high protection levels, yet apprehension exists that several variants of concerns (VOCs) can surmount the immune defenses generated by the vaccines. Neutralization assays have revealed some reduction in neutralization of VOCs B.1.1.7 and B.1.351, but the relevance of these assays in real life remains unclear. Here, we performed a case-control study that examined whether BNT162b2 vaccinees with documented SARS-CoV-2 infection were more likely to become infected with B.1.1.7 or B.1.351 compared with unvaccinated individuals. Vaccinees infected at least a week after the second dose were disproportionally infected with B.1.351 (odds ratio of 8:1). Those infected between two weeks after the first dose and one week after the second dose, were disproportionally infected by B.1.1.7 (odds ratio of 26:10), suggesting reduced vaccine effectiveness against both VOCs under different dosage/timing conditions. Nevertheless, the B.1.351 incidence in Israel to-date remains low and vaccine effectiveness remains high against B.1.1.7, among those fully vaccinated. These results overall suggest that vaccine breakthrough infection is more frequent with both VOCs, yet a combination of mass-vaccination with two doses coupled with non-pharmaceutical interventions control and contain their spread.

RNA viruses are responsible for some of the worst pandemics known to mankind, including outbreaks of Influenza, Ebola, and the recent COVID-19. One major challenge in tackling RNA viruses is the fact they are extremely genetically diverse. Nevertheless, they share common features that include their dependence on host cells for replication, and high mutation rates. We set out to search for shared evolutionary characteristics that may aid in gaining a broader understanding of RNA virus evolution, and constructed a phylogeny-based dataset spanning thousands of sequences from diverse single-stranded RNA viruses of animals. Strikingly, we found that the vast majority of these viruses have a skewed nucleotide composition, manifested as adenine rich (A-rich) coding sequences. In order to test whether A-richness is driven by selection or by biased mutation processes, we harnessed the effects of incomplete purifying selection at the tips of virus phylogenies. Our results revealed consistent mutational biases towards U rather than A in genomes of all viruses. In +ssRNA viruses we found that this bias is compensated by selection against U and selection for A, which leads to A-rich genomes. In -ssRNA viruses the genomic mutational bias towards U on the negative strand manifests as A-rich coding sequences, on the positive strand. We investigated possible reasons for the advantage of A-rich sequences including weakened RNA secondary structures, codon usage bias, and selection for a particular amino-acid composition, and conclude that host immune pressures may have led to similar biases in coding sequence composition across very divergent RNA viruses.