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.

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.

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.

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 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.

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.

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.

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.

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.

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.

miRNAs control gene expression in the single-cell alga Chlamydomonas reinhardtii.

MicroRNAs (miRNAs) in eukaryotes guide post-transcriptional regulation by means of targeted RNA degradation and translational arrest. They are released by a Dicer nuclease as a 21-24-nucleotide RNA duplex from a precursor in which an imperfectly matched inverted repeat forms a partly double-stranded region. One of the two strands is then recruited by an Argonaute nuclease that is the effector protein of the silencing mechanism. Short interfering RNAs (siRNAs), which are similar to miRNAs, are also produced by Dicer but the precursors are perfectly double-stranded RNA. These siRNAs guide post-transcriptional regulation, as with miRNAs, and epigenetic genome modification. Diverse eukaryotes including fungi, plants, protozoans and metazoans produce siRNAs but, until now, miRNAs have not been described in unicellular organisms and it has been suggested that they evolved together with multicellularity in separate plant and animal lineages. Here we show that the unicellular alga Chlamydomonas reinhardtii contains miRNAs, putative evolutionary precursors of miRNAs and species of siRNAs resembling those in higher plants. The common features of miRNAs and siRNAs in an alga and in higher plants indicate that complex RNA-silencing systems evolved before multicellularity and were a feature of primitive eukaryotic cells.

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.

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.

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.

PolIVb influences RNA-directed DNA methylation independently of its role in siRNA biogenesis.

DNA-dependent RNA polymerase (Pol)IV in Arabidopsis exists in two isoforms (PolIVa and PolIVb), with NRPD1a and NRPD1b as their respective largest subunits. Both isoforms are implicated in production and activity of siRNAs and in RNA-directed DNA methylation (RdDM). Deep sequence analysis of siRNAs in WT Arabidopsis flowers and in nrpd1a and nrpd1b mutants identified >4,200 loci producing siRNAs in a PolIV-dependent manner, with PolIVb reinforcing siRNA production by PolIVa. Transposable element identity and pericentromeric localization are both features that predispose a locus for siRNA production via PolIV proteins and determine the extent to which siRNA production relies on PolIVb. Detailed analysis of DNA methylation at PolIV-dependent loci revealed unexpected deviations from the previously noted association of PolIVb-dependent siRNA production and RdDM. Notably, PolIVb functions independently in DNA methylation and siRNA generation. Additionally, we have uncovered siRNA-directed loss of DNA methylation, a process requiring both PolIV isoforms. From these findings, we infer that the role of PolIVb in siRNA production is secondary to a role in chromatin modification and is influenced by chromatin context.

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.

Highly specific gene silencing by artificial microRNAs in the unicellular alga Chlamydomonas reinhardtii.

MicroRNAs (miRNAs) are small RNAs, 21 to 22 nucleotides long, with important regulatory roles. They are processed from longer RNA molecules with imperfectly matched foldback regions and they function in modulating the stability and translation of mRNA. Recently, we and others have demonstrated that the unicellular alga Chlamydomonas reinhardtii, like diverse multicellular organisms, contains miRNAs. These RNAs resemble the miRNAs of land plants in that they direct site-specific cleavage of target mRNA with miRNA-complementary motifs and, presumably, act as regulatory molecules in growth and development. Utilizing these findings we have developed a novel artificial miRNA system based on ligation of DNA oligonucleotides that can be used for specific high-throughput gene silencing in green algae.

Experimental identification of microRNA-140 targets by silencing and overexpressing miR-140.

MicroRNAs (miRNAs) are short noncoding RNA molecules regulating the expression of mRNAs. Target identification of miRNAs is computationally difficult due to the relatively low homology between miRNAs and their targets. We present here an experimental approach to target identification where the cartilage-specific miR-140 was overexpressed and silenced in cells it is normally expressed in separate experiments. Expression of mRNAs was profiled in both experiments and the intersection of mRNAs repressed by miR-140 overexpression and derepressed by silencing of miR-140 was identified. The intersection contained only 49 genes, although both treatments affected the accumulation of hundreds of mRNAs. These 49 genes showed a very strong enrichment for the miR-140 seed sequence implying that the approach is efficient and specific. Twenty-one of these 49 genes were predicted to be direct targets based on the presence of the seed sequence. Interestingly, none of these were predicted by the published target prediction methods we used. One of the potential target mRNAs, Cxcl12, was experimentally validated by Northern blot analysis and a luciferase reporter assay.