Representativeness of whole-genome sequencing approaches in England: the importance for understanding inequalities associated with SARS-CoV-2 infection.

Whole-genome sequencing (WGS) information has played a crucial role in the SARS-CoV-2 (COVID-19) pandemic by providing evidence about variants to inform public health policy. The purpose of this study was to assess the representativeness of sequenced cases compared with all COVID-19 cases in England, between March 2020 and August 2021, by demographic and socio-economic characteristics, to evaluate the representativeness and utility of these data in epidemiological analyses. To achieve this, polymerase chain reaction (PCR)-confirmed COVID-19 cases were extracted from the national laboratory system and linked with WGS data. During the study period, over 10% of COVID-19 cases in England had WGS data available for epidemiological analysis. With sequencing capacity increasing throughout the period, sequencing representativeness compared to all reported COVID-19 cases increased over time, allowing for valuable epidemiological analyses using demographic and socio-economic characteristics, particularly during periods with emerging novel SARS-CoV-2 variants. This study demonstrates the comprehensiveness of England's sequencing throughout the COVID-19 pandemic, rapidly detecting variants of concern, and enabling representative epidemiological analyses to inform policy.

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.

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.

The antimalarial efficacy and mechanism of resistance of the novel chemotype DDD01034957.

New antimalarial therapeutics are needed to ensure that malaria cases continue to be driven down, as both emerging parasite resistance to frontline chemotherapies and mosquito resistance to current insecticides threaten control programmes. Plasmodium, the apicomplexan parasite responsible for malaria, causes disease pathology through repeated cycles of invasion and replication within host erythrocytes (the asexual cycle). Antimalarial drugs primarily target this cycle, seeking to reduce parasite burden within the host as fast as possible and to supress recrudescence for as long as possible. Intense phenotypic drug screening efforts have identified a number of promising new antimalarial molecules. Particularly important is the identification of compounds with new modes of action within the parasite to combat existing drug resistance and suitable for formulation of efficacious combination therapies. Here we detail the antimalarial properties of DDD01034957-a novel antimalarial molecule which is fast-acting and potent against drug resistant strains in vitro, shows activity in vivo, and possesses a resistance mechanism linked to the membrane transporter PfABCI3. These data support further medicinal chemistry lead-optimization of DDD01034957 as a novel antimalarial chemical class and provide new insights to further reduce in vivo metabolic clearance.

An enhanced toolkit for the generation of knockout and marker-free fluorescent Plasmodium chabaudi.

The rodent parasite Plasmodium chabaudi is an important in vivo model of malaria. The ability to produce chronic infections makes it particularly useful for investigating the development of anti- Plasmodium immunity, as well as features associated with parasite virulence during both the acute and chronic phases of infection. P. chabaudi also undergoes asexual maturation (schizogony) and erythrocyte invasion in culture, so offers an experimentally-amenable in vivo to in vitro model for studying gene function and drug activity during parasite replication. To extend the usefulness of this model, we have further optimised transfection protocols and plasmids for P. chabaudi and generated stable, fluorescent lines that are free from drug-selectable marker genes. These mother-lines show the same infection dynamics as wild-type parasites throughout the lifecycle in mice and mosquitoes; furthermore, their virulence can be increased by serial blood passage and reset by mosquito transmission. We have also adapted the large-insert, linear PlasmoGEM vectors that have revolutionised the scale of experimental genetics in another rodent malaria parasite and used these to generate barcoded P. chabaudi gene-deletion and -tagging vectors for transfection in our fluorescent P. chabaudi mother-lines. This produces a tool-kit of P. chabaudi lines, vectors and transfection approaches that will be of broad utility to the research community.

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.

A forward genetic screen reveals a primary role for Plasmodium falciparum Reticulocyte Binding Protein Homologue 2a and 2b in determining alternative erythrocyte invasion pathways.

Invasion of human erythrocytes is essential for Plasmodium falciparum parasite survival and pathogenesis, and is also a complex phenotype. While some later steps in invasion appear to be invariant and essential, the earlier steps of recognition are controlled by a series of redundant, and only partially understood, receptor-ligand interactions. Reverse genetic analysis of laboratory adapted strains has identified multiple genes that when deleted can alter invasion, but how the relative contributions of each gene translate to the phenotypes of clinical isolates is far from clear. We used a forward genetic approach to identify genes responsible for variable erythrocyte invasion by phenotyping the parents and progeny of previously generated experimental genetic crosses. Linkage analysis using whole genome sequencing data revealed a single major locus was responsible for the majority of phenotypic variation in two invasion pathways. This locus contained the PfRh2a and PfRh2b genes, members of one of the major invasion ligand gene families, but not widely thought to play such a prominent role in specifying invasion phenotypes. Variation in invasion pathways was linked to significant differences in PfRh2a and PfRh2b expression between parasite lines, and their role in specifying alternative invasion was confirmed by CRISPR-Cas9-mediated genome editing. Expansion of the analysis to a large set of clinical P. falciparum isolates revealed common deletions, suggesting that variation at this locus is a major cause of invasion phenotypic variation in the endemic setting. This work has implications for blood-stage vaccine development and will help inform the design and location of future large-scale studies of invasion in clinical isolates.

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 genome-scale vector resource enables high-throughput reverse genetic screening in a malaria parasite.

The genome-wide identification of gene functions in malaria parasites is hampered by a lack of reverse genetic screening methods. We present a large-scale resource of barcoded vectors with long homology arms for effective modification of the Plasmodium berghei genome. Cotransfecting dozens of vectors into the haploid blood stages creates complex pools of barcoded mutants, whose competitive fitness can be measured during infection of a single mouse using barcode sequencing (barseq). To validate the utility of this resource, we rescreen the P. berghei kinome, using published kinome screens for comparison. We find that several protein kinases function redundantly in asexual blood stages and confirm the targetability of kinases cdpk1, gsk3, tkl3, and PBANKA_082960 by genotyping cloned mutants. Thus, parallel phenotyping of barcoded mutants unlocks the power of reverse genetic screening for a malaria parasite and will enable the systematic identification of genes essential for in vivo parasite growth and transmission.

PlasmoGEM, a database supporting a community resource for large-scale experimental genetics in malaria parasites.

The Plasmodium Genetic Modification (PlasmoGEM) database (http://plasmogem.sanger.ac.uk) provides access to a resource of modular, versatile and adaptable vectors for genome modification of Plasmodium spp. parasites. PlasmoGEM currently consists of >2000 plasmids designed to modify the genome of Plasmodium berghei, a malaria parasite of rodents, which can be requested by non-profit research organisations free of charge. PlasmoGEM vectors are designed with long homology arms for efficient genome integration and carry gene specific barcodes to identify individual mutants. They can be used for a wide array of applications, including protein localisation, gene interaction studies and high-throughput genetic screens. The vector production pipeline is supported by a custom software suite that automates both the vector design process and quality control by full-length sequencing of the finished vectors. The PlasmoGEM web interface allows users to search a database of finished knock-out and gene tagging vectors, view details of their designs, download vector sequence in different formats and view available quality control data as well as suggested genotyping strategies. We also make gDNA library clones and intermediate vectors available for researchers to produce vectors for themselves.

Phosphoinositide metabolism links cGMP-dependent protein kinase G to essential Ca²âº signals at key decision points in the life cycle of malaria parasites.

Many critical events in the Plasmodium life cycle rely on the controlled release of Ca²âº from intracellular stores to activate stage-specific Ca²âº-dependent protein kinases. Using the motility of Plasmodium berghei ookinetes as a signalling paradigm, we show that the cyclic guanosine monophosphate (cGMP)-dependent protein kinase, PKG, maintains the elevated level of cytosolic Ca²âº required for gliding motility. We find that the same PKG-dependent pathway operates upstream of the Ca²âº signals that mediate activation of P. berghei gametocytes in the mosquito and egress of Plasmodium falciparum merozoites from infected human erythrocytes. Perturbations of PKG signalling in gliding ookinetes have a marked impact on the phosphoproteome, with a significant enrichment of in vivo regulated sites in multiple pathways including vesicular trafficking and phosphoinositide metabolism. A global analysis of cellular phospholipids demonstrates that in gliding ookinetes PKG controls phosphoinositide biosynthesis, possibly through the subcellular localisation or activity of lipid kinases. Similarly, phosphoinositide metabolism links PKG to egress of P. falciparum merozoites, where inhibition of PKG blocks hydrolysis of phosphatidylinostitol (4,5)-bisphosphate. In the face of an increasing complexity of signalling through multiple Ca²âº effectors, PKG emerges as a unifying factor to control multiple cellular Ca²âº signals essential for malaria parasite development and transmission.

Comparative genomics in Chlamydomonas and Plasmodium identifies an ancient nuclear envelope protein family essential for sexual reproduction in protists, fungi, plants, and vertebrates.

Fertilization is a crucial yet poorly characterized event in eukaryotes. Our previous discovery that the broadly conserved protein HAP2 (GCS1) functioned in gamete membrane fusion in the unicellular green alga Chlamydomonas and the malaria pathogen Plasmodium led us to exploit the rare biological phenomenon of isogamy in Chlamydomonas in a comparative transcriptomics strategy to uncover additional conserved sexual reproduction genes. All previously identified Chlamydomonas fertilization-essential genes fell into related clusters based on their expression patterns. Out of several conserved genes in a minus gamete cluster, we focused on Cre06.g280600, an ortholog of the fertilization-related Arabidopsis GEX1. Gene disruption, cell biological, and immunolocalization studies show that CrGEX1 functions in nuclear fusion in Chlamydomonas. Moreover, CrGEX1 and its Plasmodium ortholog, PBANKA_113980, are essential for production of viable meiotic progeny in both organisms and thus for mosquito transmission of malaria. Remarkably, we discovered that the genes are members of a large, previously unrecognized family whose first-characterized member, KAR5, is essential for nuclear fusion during yeast sexual reproduction. Our comparative transcriptomics approach provides a new resource for studying sexual development and demonstrates that exploiting the data can lead to the discovery of novel biology that is conserved across distant taxa.

A Plasmodium calcium-dependent protein kinase controls zygote development and transmission by translationally activating repressed mRNAs.

Calcium-dependent protein kinases (CDPKs) play key regulatory roles in the life cycle of the malaria parasite, but in many cases their precise molecular functions are unknown. Using the rodent malaria parasite Plasmodium berghei, we show that CDPK1, which is known to be essential in the asexual blood stage of the parasite, is expressed in all life stages and is indispensable during the sexual mosquito life-cycle stages. Knockdown of CDPK1 in sexual stages resulted in developmentally arrested parasites and prevented mosquito transmission, and these effects were independent of the previously proposed function for CDPK1 in regulating parasite motility. In-depth translational and transcriptional profiling of arrested parasites revealed that CDPK1 translationally activates mRNA species in the developing zygote that in macrogametes remain repressed via their 3' and 5'UTRs. These findings indicate that CDPK1 is a multifunctional protein that translationally regulates mRNAs to ensure timely and stage-specific protein expression.

The UEA sRNA workbench: a suite of tools for analysing and visualizing next generation sequencing microRNA and small RNA datasets.

SUMMARY: RNA silencing is a complex, highly conserved mechanism mediated by small RNAs (sRNAs), such as microRNAs (miRNAs), that is known to be involved in a diverse set of biological functions including development, pathogen control, genome maintenance and response to environmental change. Advances in next generation sequencing technologies are producing increasingly large numbers of sRNA reads per sample at a fraction of the cost of previous methods. However, many bioinformatics tools do not scale accordingly, are cumbersome, or require extensive support from bioinformatics experts. Therefore, researchers need user-friendly, robust tools, capable of not only processing large sRNA datasets in a reasonable time frame but also presenting the results in an intuitive fashion and visualizing sRNA genomic features. Herein, we present the UEA sRNA workbench, a suite of tools that is a successor to the web-based UEA sRNA Toolkit, but in downloadable format and with several enhanced and additional features. AVAILABILITY: The program and help pages are available at http://srna-workbench.cmp.uea.ac.uk. CONTACT: vincent.moulton@cmp.uea.ac.uk.

PfSET10, a Plasmodium falciparum methyltransferase, maintains the active var gene in a poised state during parasite division.

A major virulence factor of the malaria parasite Plasmodium falciparum is erythrocyte membrane protein 1 (PfEMP1), a variant protein expressed on the infected erythrocyte surface. PfEMP1 is responsible for adherence of infected erythrocytes to the endothelium and plays an important role in pathogenesis. Mutually exclusive transcription and switched expression of one of 60 var genes encoding PfEMP1 in each parasite genome provides a mechanism for antigenic variation. We report the identification of a parasite protein, designated PfSET10, which localizes exclusively to the perinuclear active var gene expression site. PfSET10 is a histone 3 lysine 4 methyltransferase required to maintain the active var gene in a poised state during division for reactivation in daughter parasites, and as such is required for P. falciparum antigenic variation. PfSET10 likely maintains the transcriptionally permissive chromatin environment of the active var promoter and thus retains memory for heritable transmission of epigenetic information during parasite division.

A scalable pipeline for highly effective genetic modification of a malaria parasite.

In malaria parasites, the systematic experimental validation of drug and vaccine targets by reverse genetics is constrained by the inefficiency of homologous recombination and by the difficulty of manipulating adenine and thymine (A+T)-rich DNA of most Plasmodium species in Escherichia coli. We overcame these roadblocks by creating a high-integrity library of Plasmodium berghei genomic DNA (>77% A+T content) in a bacteriophage N15-based vector that can be modified efficiently using the lambda Red method of recombineering. We built a pipeline for generating P. berghei genetic modification vectors at genome scale in serial liquid cultures on 96-well plates. Vectors have long homology arms, which increase recombination frequency up to tenfold over conventional designs. The feasibility of efficient genetic modification at scale will stimulate collaborative, genome-wide knockout and tagging programs for P. berghei.

mRNA expression profiling reveals conserved and non-conserved miR-140 targets.

microRNAs are non-coding RNAs that regulate gene expression. A significant proportion of microRNAs is perfectly conserved across the vertebrate clade, including miR-140, which is specifically expressed in cartilage. Although it has been computationally predicted that a large majority of microRNA targets are conserved, experimental evidence for this hypothesis remains scarce. In this work we use mRNA expression profiles obtained after manipulation of miR-140 activity levels in human and chicken primary chondrocytes to explore the extent of miR-140 target conservation. Our data suggest that miR-140 has a large number of targets conserved between human and chicken and we validate one of these, BMP2. However, we also found a significant number of non-conserved targets in the two species. In addition, we found that a commercially available scrambled siRNA, which is regularly used as a negative control, regulate the accumulation of many genes.

Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns.

Plants feature a particularly diverse population of short (s)RNAs, the central component of all RNA silencing pathways. Next generation sequencing techniques enable deeper insights into this complex and highly conserved mechanism and allow identification and quantification of sRNAs. We employed deep sequencing to monitor the sRNAome of developing tomato fruits covering the period between closed flowers and ripened fruits by profiling sRNAs at 10 time-points. It is known that microRNAs (miRNAs) play an important role in development but very little information is available about the majority of sRNAs that are not miRNAs. Here we show distinctive patterns of sRNA expression that often coincide with stages of the developmental process such as flowering, early and late fruit maturation. Moreover, thousands of non-miRNA sRNAs are differentially expressed during fruit development and ripening. Some of these differentially expressed sRNAs derived from transposons but many derive from protein coding genes or regions that show homology to protein coding genes, several of which are known to play a role in flower and fruit development. These findings raise the possibility of a regulative role of these sRNAs during fruit onset and maturation in a crop species. We also identified six new miRNAs and experimentally validated two target mRNAs. These two mRNAs are targeted by the same miRNA but do not belong to the same gene family, which is rare for plant miRNAs. Expression pattern and putative function of these targets indicate a possible role in glutamate accumulation, which contributes to establishing the taste of the fruit.

The Arabidopsis RNA-directed DNA methylation argonautes functionally diverge based on their expression and interaction with target loci.

Argonaute (AGO) effectors of RNA silencing bind small RNA (sRNA) molecules and mediate mRNA cleavage, translational repression, or epigenetic DNA modification. In many organisms, these targeting mechanisms are devolved to different products of AGO multigene families. To investigate the basis of AGO functional diversification, we characterized three closely related Arabidopsis thaliana AGOs (AGO4, AGO6, and AGO9) implicated in RNA-directed DNA methylation. All three AGOs bound 5' adenosine 24-nucleotide sRNAs, but each exhibited different preferences for sRNAs from different heterochromatin-associated loci. This difference was reduced when AGO6 and AGO9 were expressed from the AGO4 promoter, indicating that the functional diversification was partially due to differential expression of the corresponding genes. However, the AGO4-directed pattern of sRNA accumulation and DNA methylation was not fully recapitulated with AGO6 or AGO9 expressed from the AGO4 promoter. Here, we show that sRNA length and 5' nucleotide do not account for the observed functional diversification of these AGOs. Instead, the selectivity of sRNA binding is determined by the coincident expression of the AGO and sRNA-generating loci, and epigenetic modification is influenced by interactions between the AGO protein and the different target loci. These findings highlight the importance of tissue specificity and AGO-associated proteins in influencing epigenetic modifications.