Genome-Scale Identification of Essential Metabolic Processes for Targeting the Plasmodium Liver Stage.

Plasmodium gene functions in mosquito and liver stages remain poorly characterized due to limitations in the throughput of phenotyping at these stages. To fill this gap, we followed more than 1,300 barcoded P. berghei mutants through the life cycle. We discover 461 genes required for efficient parasite transmission to mosquitoes through the liver stage and back into the bloodstream of mice. We analyze the screen in the context of genomic, transcriptomic, and metabolomic data by building a thermodynamic model of P. berghei liver-stage metabolism, which shows a major reprogramming of parasite metabolism to achieve rapid growth in the liver. We identify seven metabolic subsystems that become essential at the liver stages compared with asexual blood stages: type II fatty acid synthesis and elongation (FAE), tricarboxylic acid, amino sugar, heme, lipoate, and shikimate metabolism. Selected predictions from the model are individually validated in single mutants to provide future targets for drug development.

Functional Profiling of a Plasmodium Genome Reveals an Abundance of Essential Genes.

The genomes of malaria parasites contain many genes of unknown function. To assist drug development through the identification of essential genes and pathways, we have measured competitive growth rates in mice of 2,578 barcoded Plasmodium berghei knockout mutants, representing >50% of the genome, and created a phenotype database. At a single stage of its complex life cycle, P. berghei requires two-thirds of genes for optimal growth, the highest proportion reported from any organism and a probable consequence of functional optimization necessitated by genomic reductions during the evolution of parasitism. In contrast, extreme functional redundancy has evolved among expanded gene families operating at the parasite-host interface. The level of genetic redundancy in a single-celled organism may thus reflect the degree of environmental variation it experiences. In the case of Plasmodium parasites, this helps rationalize both the relative successes of drugs and the greater difficulty of making an effective vaccine.

A molnupiravir-associated mutational signature in global SARS-CoV-2 genomes.

Molnupiravir, an antiviral medication widely used against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), acts by inducing mutations in the virus genome during replication. Most random mutations are likely to be deleterious to the virus and many will be lethal; thus, molnupiravir-induced elevated mutation rates reduce viral load1,2. However, if some patients treated with molnupiravir do not fully clear the SARS-CoV-2 infections, there could be the potential for onward transmission of molnupiravir-mutated viruses. Here we show that SARS-CoV-2 sequencing databases contain extensive evidence of molnupiravir mutagenesis. Using a systematic approach, we find that a specific class of long phylogenetic branches, distinguished by a high proportion of G-to-A and C-to-T mutations, are found almost exclusively in sequences from 2022, after the introduction of molnupiravir treatment, and in countries and age groups with widespread use of the drug. We identify a mutational spectrum, with preferred nucleotide contexts, from viruses in patients known to have been treated with molnupiravir and show that its signature matches that seen in these long branches, in some cases with onward transmission of molnupiravir-derived lineages. Finally, we analyse treatment records to confirm a direct association between these high G-to-A branches and the use of molnupiravir.

Genomic reconstruction of the SARS-CoV-2 epidemic in England.

The evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus leads to new variants that warrant timely epidemiological characterization. Here we use the dense genomic surveillance data generated by the COVID-19 Genomics UK Consortium to reconstruct the dynamics of 71 different lineages in each of 315 English local authorities between September 2020 and June 2021. This analysis reveals a series of subepidemics that peaked in early autumn 2020, followed by a jump in transmissibility of the B.1.1.7/Alpha lineage. The Alpha variant grew when other lineages declined during the second national lockdown and regionally tiered restrictions between November and December 2020. A third more stringent national lockdown suppressed the Alpha variant and eliminated nearly all other lineages in early 2021. Yet a series of variants (most of which contained the spike E484K mutation) defied these trends and persisted at moderately increasing proportions. However, by accounting for sustained introductions, we found that the transmissibility of these variants is unlikely to have exceeded the transmissibility of the Alpha variant. Finally, B.1.617.2/Delta was repeatedly introduced in England and grew rapidly in early summer 2021, constituting approximately 98% of sampled SARS-CoV-2 genomes on 26 June 2021.

Treenome Browser: co-visualization of enormous phylogenies and millions of genomes.

SUMMARY: Treenome Browser is a web browser tool to interactively visualize millions of genomes alongside huge phylogenetic trees. AVAILABILITY AND IMPLEMENTATION: Treenome Browser for SARS-CoV-2 can be accessed at cov2tree.org, or at taxonium.org for user-provided trees. Source code and documentation are available at github.com/theosanderson/taxonium and docs.taxonium.org/en/latest/treenome.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Using Plasmodium knowlesi as a model for screening Plasmodium vivax blood-stage malaria vaccine targets reveals new candidates.

Plasmodium vivax is responsible for the majority of malaria cases outside Africa. Unlike P. falciparum, the P. vivax life-cycle includes a dormant liver stage, the hypnozoite, which can cause infection in the absence of mosquito transmission. An effective vaccine against P. vivax blood stages would limit symptoms and pathology from such recurrent infections, and therefore could play a critical role in the control of this species. Vaccine development in P. vivax, however, lags considerably behind P. falciparum, which has many identified targets with several having transitioned to Phase II testing. By contrast only one P. vivax blood-stage vaccine candidate based on the Duffy Binding Protein (PvDBP), has reached Phase Ia, in large part because the lack of a continuous in vitro culture system for P. vivax limits systematic screening of new candidates. We used the close phylogenetic relationship between P. vivax and P. knowlesi, for which an in vitro culture system in human erythrocytes exists, to test the scalability of systematic reverse vaccinology to identify and prioritise P. vivax blood-stage targets. A panel of P. vivax proteins predicted to function in erythrocyte invasion were expressed as full-length recombinant ectodomains in a mammalian expression system. Eight of these antigens were used to generate polyclonal antibodies, which were screened for their ability to recognize orthologous proteins in P. knowlesi. These antibodies were then tested for inhibition of growth and invasion of both wild type P. knowlesi and chimeric P. knowlesi lines modified using CRISPR/Cas9 to exchange P. knowlesi genes with their P. vivax orthologues. Candidates that induced antibodies that inhibited invasion to a similar level as PvDBP were identified, confirming the utility of P. knowlesi as a model for P. vivax vaccine development and prioritizing antigens for further follow up.

Genomic and transcriptomic evidence for descent from Plasmodium and loss of blood schizogony in Hepatocystis parasites from naturally infected red colobus monkeys.

Hepatocystis is a genus of single-celled parasites infecting, amongst other hosts, monkeys, bats and squirrels. Although thought to have descended from malaria parasites (Plasmodium spp.), Hepatocystis spp. are thought not to undergo replication in the blood-the part of the Plasmodium life cycle which causes the symptoms of malaria. Furthermore, Hepatocystis is transmitted by biting midges, not mosquitoes. Comparative genomics of Hepatocystis and Plasmodium species therefore presents an opportunity to better understand some of the most important aspects of malaria parasite biology. We were able to generate a draft genome for Hepatocystis sp. using DNA sequencing reads from the blood of a naturally infected red colobus monkey. We provide robust phylogenetic support for Hepatocystis sp. as a sister group to Plasmodium parasites infecting rodents. We show transcriptomic support for a lack of replication in the blood and genomic support for a complete loss of a family of genes involved in red blood cell invasion. Our analyses highlight the rapid evolution of genes involved in parasite vector stages, revealing genes that may be critical for interactions between malaria parasites and mosquitoes.

CRISPR/Cas9 and genetic screens in malaria parasites: small genomes, big impact.

The ∼30 Mb genomes of the Plasmodium parasites that cause malaria each encode ∼5000 genes, but the functions of the majority remain unknown. This is due to a paucity of functional annotation from sequence homology, which is compounded by low genetic tractability compared with many model organisms. In recent years technical breakthroughs have made forward and reverse genome-scale screens in Plasmodium possible. Furthermore, the adaptation of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-Associated protein 9 (CRISPR/Cas9) technology has dramatically improved gene editing efficiency at the single gene level. Here, we review the arrival of genetic screens in malaria parasites to analyse parasite gene function at a genome-scale and their impact on understanding parasite biology. CRISPR/Cas9 screens, which have revolutionised human and model organism research, have not yet been implemented in malaria parasites due to the need for more complex CRISPR/Cas9 gene targeting vector libraries. We therefore introduce the reader to CRISPR-based screens in the related apicomplexan Toxoplasma gondii and discuss how these approaches could be adapted to develop CRISPR/Cas9 based genome-scale genetic screens in malaria parasites. Moreover, since more than half of Plasmodium genes are required for normal asexual blood-stage reproduction, and cannot be targeted using knockout methods, we discuss how CRISPR/Cas9 could be used to scale up conditional gene knockdown approaches to systematically assign function to essential genes.

An apicoplast-resident folate transporter is essential for sporogony of malaria parasites.

Malaria parasites are fast replicating unicellular organisms and require substantial amounts of folate for DNA synthesis. Despite the central role of this critical co-factor for parasite survival, only little is known about intraparasitic folate trafficking in Plasmodium. Here, we report on the expression, subcellular localisation and function of the parasite's folate transporter 2 (FT2) during life cycle progression in the murine malaria parasite Plasmodium berghei. Using live fluorescence microscopy of genetically engineered parasites, we demonstrate that FT2 localises to the apicoplast. In invasive P. berghei stages, a fraction of FT2 is also observed at the apical end. Upon genetic disruption of FT2, blood and liver infection, gametocyte production and mosquito colonisation remain unaltered. But in the Anopheles vector, FT2-deficient parasites develop inflated oocysts with unusual pulp formation consisting of numerous single-membrane vesicles, which ultimately fuse to form large cavities. Ultrastructural analysis suggests that this defect reflects aberrant sporoblast formation caused by abnormal vesicular traffic. Complete sporogony in FT2-deficient oocysts is very rare, and mutant sporozoites fail to establish hepatocyte infection, resulting in a complete block of parasite transmission. Our findings reveal a previously unrecognised organellar folate transporter that exerts critical roles for pathogen maturation in the arthropod vector.

Synergistic malaria vaccine combinations identified by systematic antigen screening.

A highly effective vaccine would be a valuable weapon in the drive toward malaria elimination. No such vaccine currently exists, and only a handful of the hundreds of potential candidates in the parasite genome have been evaluated. In this study, we systematically evaluated 29 antigens likely to be involved in erythrocyte invasion, an essential developmental stage during which the malaria parasite is vulnerable to antibody-mediated inhibition. Testing antigens alone and in combination identified several strain-transcending targets that had synergistic combinatorial effects in vitro, while studies in an endemic population revealed that combinations of the same antigens were associated with protection from febrile malaria. Video microscopy established that the most effective combinations targeted multiple discrete stages of invasion, suggesting a mechanistic explanation for synergy. Overall, this study both identifies specific antigen combinations for high-priority clinical testing and establishes a generalizable approach that is more likely to produce effective vaccines.

Randomized controlled trial of molnupiravir SARS-CoV-2 viral and antibody response in at-risk adult outpatients.

Viral clearance, antibody response and the mutagenic effect of molnupiravir has not been elucidated in at-risk populations. Non-hospitalised participants within 5 days of SARS-CoV-2 symptoms randomised to receive molnupiravir (n = 253) or Usual Care (n = 324) were recruited to study viral and antibody dynamics and the effect of molnupiravir on viral whole genome sequence from 1437 viral genomes. Molnupiravir accelerates viral load decline, but virus is detectable by Day 5 in most cases. At Day 14 (9 days post-treatment), molnupiravir is associated with significantly higher viral persistence and significantly lower anti-SARS-CoV-2 spike antibody titres compared to Usual Care. Serial sequencing reveals increased mutagenesis with molnupiravir treatment. Persistence of detectable viral RNA at Day 14 in the molnupiravir group is associated with higher transition mutations following treatment cessation. Viral viability at Day 14 is similar in both groups with post-molnupiravir treated samples cultured up to 9 days post cessation of treatment. The current 5-day molnupiravir course is too short. Longer courses should be tested to reduce the risk of potentially transmissible molnupiravir-mutated variants being generated. Trial registration: ISRCTN30448031.

Addressing pandemic-wide systematic errors in the SARS-CoV-2 phylogeny.

The SARS-CoV-2 genome occupies a unique place in infection biology - it is the most highly sequenced genome on earth (making up over 20% of public sequencing datasets) with fine scale information on sampling date and geography, and has been subject to unprecedented intense analysis. As a result, these phylogenetic data are an incredibly valuable resource for science and public health. However, the vast majority of the data was sequenced by tiling amplicons across the full genome, with amplicon schemes that changed over the pandemic as mutations in the viral genome interacted with primer binding sites. In combination with the disparate set of genome assembly workflows and lack of consistent quality control (QC) processes, the current genomes have many systematic errors that have evolved with the virus and amplicon schemes. These errors have significant impacts on the phylogeny, and therefore over the last few years, many thousands of hours of researchers time has been spent in "eyeballing" trees, looking for artefacts, and then patching the tree. Given the huge value of this dataset, we therefore set out to reprocess the complete set of public raw sequence data in a rigorous amplicon-aware manner, and build a cleaner phylogeny. Here we provide a global tree of 3,960,704 samples, built from a consistently assembled set of high quality consensus sequences from all available public data as of March 2023, viewable at https://viridian.taxonium.org. Each genome was constructed using a novel assembly tool called Viridian (https://github.com/iqbal-lab-org/viridian), developed specifically to process amplicon sequence data, eliminating artefactual errors and mask the genome at low quality positions. We provide simulation and empirical validation of the methodology, and quantify the improvement in the phylogeny. Phase 2 of our project will address the fact that the data in the public archives is heavily geographically biased towards the Global North. We therefore have contributed new raw data to ENA/SRA from many countries including Ghana, Thailand, Laos, Sri Lanka, India, Argentina and Singapore. We will incorporate these, along with all public raw data submitted between March 2023 and the current day, into an updated set of assemblies, and phylogeny. We hope the tree, consensus sequences and Viridian will be a valuable resource for researchers.

ProteInfer, deep neural networks for protein functional inference.

Predicting the function of a protein from its amino acid sequence is a long-standing challenge in bioinformatics. Traditional approaches use sequence alignment to compare a query sequence either to thousands of models of protein families or to large databases of individual protein sequences. Here we introduce ProteInfer, which instead employs deep convolutional neural networks to directly predict a variety of protein functions - Enzyme Commission (EC) numbers and Gene Ontology (GO) terms - directly from an unaligned amino acid sequence. This approach provides precise predictions which complement alignment-based methods, and the computational efficiency of a single neural network permits novel and lightweight software interfaces, which we demonstrate with an in-browser graphical interface for protein function prediction in which all computation is performed on the user's personal computer with no data uploaded to remote servers. Moreover, these models place full-length amino acid sequences into a generalised functional space, facilitating downstream analysis and interpretation. To read the interactive version of this paper, please visit https://google-research.github.io/proteinfer/.

Regulators of male and female sexual development are critical for the transmission of a malaria parasite.

Malaria transmission to mosquitoes requires a developmental switch in asexually dividing blood-stage parasites to sexual reproduction. In Plasmodium berghei, the transcription factor AP2-G is required and sufficient for this switch, but how a particular sex is determined in a haploid parasite remains unknown. Using a global screen of barcoded mutants, we here identify genes essential for the formation of either male or female sexual forms and validate their importance for transmission. High-resolution single-cell transcriptomics of ten mutant parasites portrays the developmental bifurcation and reveals a regulatory cascade of putative gene functions in the determination and subsequent differentiation of each sex. A male-determining gene with a LOTUS/OST-HTH domain as well as the protein interactors of a female-determining zinc-finger protein indicate that germ-granule-like ribonucleoprotein complexes complement transcriptional processes in the regulation of both male and female development of a malaria parasite.

Taxonium, a web-based tool for exploring large phylogenetic trees.

The COVID-19 pandemic has resulted in a step change in the scale of sequencing data, with more genomes of SARS-CoV-2 having been sequenced than any other organism on earth. These sequences reveal key insights when represented as a phylogenetic tree, which captures the evolutionary history of the virus, and allows the identification of transmission events and the emergence of new variants. However, existing web-based tools for exploring phylogenies do not scale to the size of datasets now available for SARS-CoV-2. We have developed Taxonium, a new tool that uses WebGL to allow the exploration of trees with tens of millions of nodes in the browser for the first time. Taxonium links each node to associated metadata and supports mutation-annotated trees, which are able to capture all known genetic variation in a dataset. It can either be run entirely locally in the browser, from a server-based backend, or as a desktop application. We describe insights that analysing a tree of five million sequences can provide into SARS-CoV-2 evolution, and provide a tool at cov2tree.org for exploring a public tree of more than five million SARS-CoV-2 sequences. Taxonium can be applied to any tree, and is available at taxonium.org, with source code at github.com/theosanderson/taxonium.

Since 2020, the SARS-CoV-2 virus has infected billions of people and spread to 185 countries. The virus spreads by making new copies of its genome inside human cells and exploits the cells' machinery to synthesise viral proteins it needs to infect further cells. Each time the virus copies its genetic material there's a chance that the replication process introduces an error to the genetic sequence. Over time, these mutations accumulate which can give rise to new variants with different properties. These new variants, originating from a common ancestor, may spread faster or be able to evade immune systems that have learnt to recognise previous variants. To understand where new variants of SARS-CoV-2 come from and how related they are to each other, scientists build family trees called 'phylogenetic trees' based on similarities in the genetic sequences of different variants of the virus. Looking at these trees researchers can track how a variant spreads geographically, and also attempt to identify new worrying variants that might lead to a new wave of infections. The scale of the COVID-19 pandemic together with the global effort by clinicians and researchers to sequence SARS-CoV-2 genetic material means a library of over 13 million SARS-CoV-2 genomes now exists, making it the largest such collection for any organism. Although phylogenetic trees of viruses have been studied for a long time, exploring the SARS-CoV-2 library presents technical and practical challenges due to its sheer size. Sanderson has developed an open-source web tool called Taxonium that allows users to explore phylogenetic trees with millions of sequences. With help from collaborators at the University of California, Santa Cruz, Sanderson built a website called Cov2Tree, that uses the Taxonium platform to allow immediate access to an expansive tree of all publicly available SARS-CoV-2 sequences. Cov2Tree enables users to visualise all SARS-CoV-2 genomes in a birds-eye view akin to a 'Google Earth for virus sequences' where anyone can zoom in on a related family of viruses down to the level of individual sequences. This can be used to compare variants and follow geographic spread. Using Taxonium, scientists can explore how virus sequences are related to each other. They can also see the individual mutations that have occurred at each branch of the tree, and can search for sequences based on mutation, geographical location, or other factors. Interestingly, a trend appearing in the SARS-CoV-2 phylogenetic tree is the emergence of identical mutations at different branches of the tree without a common origin. These mutations may be a result of convergent evolution, a phenomenon that occurs when a mutation appears independently in different variants as it confers an advantage to the virus making such mutations more likely to persist. This means that scientists may be able to expect certain mutations to appear in more distantly related variants if they have appeared independently in several different variants already. Overall, Taxonium is an important tool for monitoring SARS-CoV-2 genomes, but it also has broader applications. The tool can be used to browse phylogenetic trees of other viruses and organisms. Furthermore, the Taxonium website offers a way to browse a tree of life, with images and links to Wikipedia. The SARS-CoV-2 library might be the largest now, but in the future even bigger datasets will likely be available, highlighting the importance of tools like Taxonium.

Using deep learning to annotate the protein universe.

Understanding the relationship between amino acid sequence and protein function is a long-standing challenge with far-reaching scientific and translational implications. State-of-the-art alignment-based techniques cannot predict function for one-third of microbial protein sequences, hampering our ability to exploit data from diverse organisms. Here, we train deep learning models to accurately predict functional annotations for unaligned amino acid sequences across rigorous benchmark assessments built from the 17,929 families of the protein families database Pfam. The models infer known patterns of evolutionary substitutions and learn representations that accurately cluster sequences from unseen families. Combining deep models with existing methods significantly improves remote homology detection, suggesting that the deep models learn complementary information. This approach extends the coverage of Pfam by >9.5%, exceeding additions made over the last decade, and predicts function for 360 human reference proteome proteins with no previous Pfam annotation. These results suggest that deep learning models will be a core component of future protein annotation tools.

Determinants of Spike infectivity, processing, and neutralization in SARS-CoV-2 Omicron subvariants BA.1 and BA.2.

SARS-CoV-2 Omicron rapidly outcompeted other variants and currently dominates the COVID-19 pandemic. Its enhanced transmission and immune evasion are thought to be driven by numerous mutations in the Omicron Spike protein. Here, we systematically introduced BA.1 and/or BA.2 Omicron Spike mutations into the ancestral Spike protein and examined the impacts on Spike function, processing, and susceptibility to neutralization. Individual mutations of S371F/L, S375F, and T376A in the ACE2-receptor-binding domain as well as Q954H and N969K in the hinge region 1 impaired infectivity, while changes to G339D, D614G, N764K, and L981F moderately enhanced it. Most mutations in the N-terminal region and receptor-binding domain reduced the sensitivity of the Spike protein to neutralization by sera from individuals vaccinated with the BNT162b2 vaccine and by therapeutic antibodies. Our results represent a systematic functional analysis of Omicron Spike adaptations that have allowed this SARS-CoV-2 variant to dominate the current pandemic.