GWAS and meta-analysis identifies 49 genetic variants underlying critical COVID-19.

Critical illness in COVID-19 is an extreme and clinically homogeneous disease phenotype that we have previously shown1 to be highly efficient for discovery of genetic associations2. Despite the advanced stage of illness at presentation, we have shown that host genetics in patients who are critically ill with COVID-19 can identify immunomodulatory therapies with strong beneficial effects in this group3. Here we analyse 24,202 cases of COVID-19 with critical illness comprising a combination of microarray genotype and whole-genome sequencing data from cases of critical illness in the international GenOMICC (11,440 cases) study, combined with other studies recruiting hospitalized patients with a strong focus on severe and critical disease: ISARIC4C (676 cases) and the SCOURGE consortium (5,934 cases). To put these results in the context of existing work, we conduct a meta-analysis of the new GenOMICC genome-wide association study (GWAS) results with previously published data. We find 49 genome-wide significant associations, of which 16 have not been reported previously. To investigate the therapeutic implications of these findings, we infer the structural consequences of protein-coding variants, and combine our GWAS results with gene expression data using a monocyte transcriptome-wide association study (TWAS) model, as well as gene and protein expression using Mendelian randomization. We identify potentially druggable targets in multiple systems, including inflammatory signalling (JAK1), monocyte-macrophage activation and endothelial permeability (PDE4A), immunometabolism (SLC2A5 and AK5), and host factors required for viral entry and replication (TMPRSS2 and RAB2A).

Noradrenergic Signaling and Neuroinflammation Crosstalk Regulate Toxoplasma gondii-Induced Behavioral Changes.

Infections of the nervous system elicit neuroimmune responses and alter neurotransmission, affecting host neurological functions. Chronic infection with the apicomplexan parasite Toxoplasma correlates with certain neurological disorders in humans and alters behavior in rodents. Here, we propose that the crosstalk between neurotransmission and neuroinflammation may underlie some of these cognitive changes. We discuss how T. gondii infection suppresses noradrenergic signaling and how the restoration of this pathway improves behavioral aberrations, suggesting that altered neurotransmission and neuroimmune responses may act in concert to perturb behavior. This interaction might apply to other infectious agents, such as viruses, that elicit cognitive changes. We hypothesize that neurotransmitter signaling in immune cells can contribute to behavioral changes associated with brain infection, offering opportunities for potential therapeutic targeting.

Downregulation of the Central Noradrenergic System by Toxoplasma gondii Infection.

Toxoplasma gondii is associated with physiological effects in the host. Dysregulation of catecholamines in the central nervous system has previously been observed in chronically infected animals. In the study described here, the noradrenergic system was found to be suppressed with decreased levels of norepinephrine (NE) in brains of infected animals and in infected human and rat neural cells in vitro The mechanism responsible for the NE suppression was found to be downregulation of dopamine ß-hydroxylase (DBH) gene expression, encoding the enzyme that synthesizes norepinephrine from dopamine, with downregulation observed in vitro and in infected brain tissue, particularly in the dorsal locus coeruleus/pons region. The downregulation was sex specific, with males expressing reduced DBH mRNA levels whereas females were unchanged. Rather, DBH expression correlated with estrogen receptor in the female rat brains for this estrogen-regulated gene. DBH silencing was not a general response of neurons to infection, as human cytomegalovirus did not downregulate DBH expression. The noradrenergic-linked behaviors of sociability and arousal were altered in chronically infected animals, with a high correlation between DBH expression and infection intensity. A decrease in DBH expression in noradrenergic neurons can elevate dopamine levels, which provides a possible explanation for mixed observations of changes in this neurotransmitter with infection. Decreased NE is consistent with the loss of coordination and motor impairments associated with toxoplasmosis. Further, the altered norepinephrine synthesis observed here may, in part, explain behavioral effects of infection and associations with mental illness.

Interconvertible geometric isomers of Plasmodium falciparum dihydroorotate dehydrogenase inhibitors exhibit multiple binding modes.

Two new tricyclic ß-aminoacrylate derivatives (2e and 3e) have been found to be inhibitors of Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH) with Ki 0.037 and 0.15µM respectively. 1H and 13C NMR spectroscopic data show that these compounds undergo ready cis-trans isomerisation at room temperature in polar solvents. In silico docking studies indicate that for both molecules there is neither conformation nor double bond configuration which bind preferentially to PfDHODH. This flexibility is favourable for inhibitors of this channel that require extensive positioning to reach their binding site.

Pantothenic acid biosynthesis in the parasite Toxoplasma gondii: a target for chemotherapy.

Toxoplasma gondii is a major food pathogen and neglected parasitic infection that causes eye disease, birth defects, and fetal abortion and plays a role as an opportunistic infection in AIDS. In this study, we investigated pantothenic acid (vitamin B5) biosynthesis in T. gondii. Genes encoding the full repertoire of enzymes for pantothenate synthesis and subsequent metabolism to coenzyme A were identified and are expressed in T. gondii. A panel of inhibitors developed to target Mycobacterium tuberculosis pantothenate synthetase were tested and found to exhibit a range of values for inhibition of T. gondii growth. Two inhibitors exhibited lower effective concentrations than the currently used toxoplasmosis drug pyrimethamine. The inhibition was specific for the pantothenate pathway, as the effect of the pantothenate synthetase inhibitors was abrogated by supplementation with pantothenate. Hence, T. gondii encodes and expresses the enzymes for pantothenate synthesis, and this pathway is essential for parasite growth. These promising findings increase our understanding of growth and metabolism in this important parasite and highlight pantothenate synthetase as a new drug target.

Toxoplasma gondii infection and behaviour - location, location, location?

Parasite location has been proposed as an important factor in the behavioural changes observed in rodents infected with the protozoan Toxoplasma gondii. During the chronic stages of infection, encysted parasites are found in the brain but it remains unclear whether the parasite has tropism for specific brain regions. Parasite tissue cysts are found in all brain areas with some, but not all, prior studies reporting higher numbers located in the amygdala and frontal cortex. A stochastic process of parasite location does not, however, seem to explain the distinct and often subtle changes observed in rodent behaviour. One factor that could contribute to the specific changes is increased dopamine production by T. gondii. Recently, it was found that cells encysted with parasites in the brains of experimentally infected rodents have high levels of dopamine and that the parasite encodes a tyrosine hydroxylase, the rate-limiting enzyme in the synthesis of this neurotransmitter. A mechanism is proposed that could explain the behaviour changes due to parasite regulation of dopamine. This could have important implications for T. gondii infections in humans.

Toxoplasma gondii infection, from predation to schizophrenia: can animal behaviour help us understand human behaviour?

We examine the role of the protozoan Toxoplasma gondii as a manipulatory parasite and question what role study of infections in its natural intermediate rodent hosts and other secondary hosts, including humans, may elucidate in terms of the epidemiology, evolution and clinical applications of infection. In particular, we focus on the potential association between T. gondii and schizophrenia. We introduce the novel term 'T. gondii-rat manipulation-schizophrenia model' and propose how future behavioural research on this model should be performed from a biological, clinical and ethically appropriate perspective.

GABAergic regulation of cell proliferation within the adult mouse spinal cord.

Manipulation of neural stem cell proliferation and differentiation in the postnatal CNS is receiving significant attention due to therapeutic potential. In the spinal cord, such manipulations may promote repair in conditions such as multiple sclerosis or spinal cord injury, but may also limit excessive cell proliferation contributing to tumours such as ependymomas. We show that when ambient γ-aminobutyric acid (GABA) is increased in vigabatrin-treated or decreased by GAD67 allele haplodeficiency in glutamic acid decarboxylase67-green fluorescent protein (GAD67-GFP) mice of either sex, the numbers of proliferating cells respectively decreased or increased. Thus, intrinsic spinal cord GABA levels are correlated with the extent of cell proliferation, providing important evidence for manipulating these levels. Diazepam binding inhibitor, an endogenous protein that interacts with GABA receptors and its breakdown product, octadecaneuropeptide, which preferentially activates central benzodiazepine (CBR) sites, were highly expressed in spinal cord, especially in ependymal cells surrounding the central canal. Furthermore, animals with reduced CBR activation via treatment with flumazenil or Ro15-4513, or with a G2F77I mutation in the CBR binding site had greater numbers of Ethynyl-2'-deoxyuridine positive cells compared to control, which maintained their stem cell status since the proportion of newly proliferated cells becoming oligodendrocytes or astrocytes was significantly lower. Altering endogenous GABA levels or modulating GABAergic signalling through specific sites on GABA receptors therefore influences NSC proliferation in the adult spinal cord. These findings provide a basis for further study into how GABAergic signalling could be manipulated to enable spinal cord self-regeneration and recovery or limit pathological proliferative activity.

MetNetMaker: a free and open-source tool for the creation of novel metabolic networks in SBML format.

UNLABELLED: An application has been developed to help with the creation and editing of Systems Biology Markup Language (SBML) format metabolic networks up to the organism scale. Networks are defined as a collection of Kyoto Encyclopedia of Genes and Genomes (KEGG) LIGAND reactions with an optional associated Enzyme Classification (EC) number for each reaction. Additional custom reactions can be defined by the user. Reactions within the network can be assigned flux constraints and compartmentalization is supported for each reaction in addition to the support for reactions that occur across compartment boundaries. Exported networks are fully SBML L2V4 compatible with an optional L2V1 export for compatibility with old versions of the COBRA toolbox. AVAILABILITY AND IMPLEMENTATION: The software runs in the free Microsoft Access 2007 Runtime (Microsoft Inc.), which is included with the installer and works on Windows XP SP2 or better. Full source code is viewable in the full version of Access 2007 or 2010. Users must have a license to use the KEGG LIGAND database (free academic licensing is available). Please go to www.bioinformatics.leeds.ac.uk/~pytf/metnetmaker for software download, help and tutorials.

Gene function prediction using semantic similarity clustering and enrichment analysis in the malaria parasite Plasmodium falciparum.

MOTIVATION: Functional genomics data provides a rich source of information that can be used in the annotation of the thousands of genes of unknown function found in most sequenced genomes. However, previous gene function prediction programs are mostly produced for relatively well-annotated organisms that often have a large amount of functional genomics data. Here, we present a novel method for predicting gene function that uses clustering of genes by semantic similarity, a naïve Bayes classifier and 'enrichment analysis' to predict gene function for a genome that is less well annotated but does has a severe effect on human health, that of the malaria parasite Plasmodium falciparum. RESULTS: Predictions for the molecular function, biological process and cellular component of P.falciparum genes were created from eight different datasets with a combined prediction also being produced. The high-confidence predictions produced by the combined prediction were compared to those produced by a simple K-nearest neighbour classifier approach and were shown to improve accuracy and coverage. Finally, two case studies are described, which investigate two biological processes in more detail, that of translation initiation and invasion of the host cell. AVAILABILITY: Predictions produced are available at http://www.bioinformatics.leeds.ac.uk/∼bio5pmrt/PAGODA.

metaTIGER: a metabolic evolution resource.

Metabolic networks are a subject that has received much attention, but existing web resources do not include extensive phylogenetic information. Phylogenomic approaches (phylogenetics on a genomic scale) have been shown to be effective in the study of evolution and processes like horizontal gene transfer (HGT). To address the lack of phylogenomic information relating to eukaryotic metabolism, metaTIGER (www.bioinformatics.leeds.ac.uk/metatiger) has been created, using genomic information from 121 eukaryotes and 404 prokaryotes and sensitive sequence search techniques to predict the presence of metabolic enzymes. These enzyme sequences were used to create a comprehensive database of 2257 maximum-likelihood phylogenetic trees, some containing over 500 organisms. The trees can be viewed using iTOL, an advanced interactive tree viewer, enabling straightforward interpretation of large trees. Complex high-throughput tree analysis is also available through user-defined queries, allowing the rapid identification of trees of interest, e.g. containing putative HGT events. metaTIGER also provides novel and easy-to-use facilities for viewing and comparing the metabolic networks in different organisms via highlighted pathway images and tables. metaTIGER is demonstrated through evolutionary analysis of Plasmodium, including identification of genes horizontally transferred from chlamydia.

A study of the effects of substituents on the selectivity of the binding of N-arylaminomethylene malonate inhibitors to DHODH.

A series of mono- and di-substituted N-arylaminomethylene malonates have been used to probe the potential of utilizing additional H-bonding contacts in the ubiquinone binding channel, for selective inhibition between either human or Plasmodium DHODH. Altered 'head' group functionalities have been utilized in order to probe the role of specific functionalities within the inhibitors in terms of enzyme affinity and selectivity.

Toxoplasma gondii-altered host behaviour: clues as to mechanism of action.

A convincing body of evidence now exists, from both human and animal studies, and encompassing epidemiological to experimental, to indicate that the common protozoan Toxoplasma gondii can cause specific behavioural changes in its host. Such behavioural alterations are likely to be the product of strong selective pressures for the parasite to enhance transmission from its intermediate host reservoir, primarily rodent, to its feline definitive host, wherein sexual reproduction can occur and the parasite's life cycle completed. Here we consider what the available data to date may reveal about the potential mechanisms involved, the future research that needs to be performed, and the subsequent implications for animal and human health.

Prediction of horizontal gene transfers in eukaryotes: approaches and challenges.

HGT (horizontal gene transfer) is recognized as an important force in bacterial evolution. Now that many eukaryotic genomes have been sequenced, it has become possible to carry out studies of HGT in eukaryotes. The present review compares the different approaches that exist for identifying HGT genes and assess them in the context of studying eukaryotic evolution. The metabolic evolution resource metaTIGER is then described, with discussion of its application in identification of HGT in eukaryotes.

metaSHARK: a WWW platform for interactive exploration of metabolic networks.

The metaSHARK (metabolic search and reconstruction kit) web server offers users an intuitive, fully interactive way to explore the KEGG metabolic network via a WWW browser. Metabolic reconstruction information for specific organisms, produced by our automated SHARKhunt tool or from other programs or genome annotations, may be uploaded to the website and overlaid on the generic network. Additional data from gene expression experiments can also be incorporated, allowing the visualization of differential gene expression in the context of the predicted metabolic network. metaSHARK is available at http://bioinformatics.leeds.ac.uk/shark/.

Metabolic reconstruction and analysis for parasite genomes.

With the completion of sequencing projects for several parasite genomes, efforts are ongoing to make sense of this mass of information in terms of the gene products encoded and their interactions in the growth, development and survival of parasites. The emerging science of systems biology aims to explain the complex relationship between genotype and phenotype by using network models. One area in which this approach has been particularly successful is in the modeling of metabolism. With an accurate picture of the set of metabolic reactions encoded in a genome, it is now possible to identify enzymes or transporters that might be viable targets for new drugs. Because these predictions greatly depend on the quality and completeness of the genome annotation, there are substantial efforts in the scientific community to increase the numbers of metabolic enzymes identified. In this review, we discuss the opportunities for using metabolic reconstruction and analysis tools in parasitology research, and their applications to protozoan parasites.

Design and synthesis of potent inhibitors of the malaria parasite dihydroorotate dehydrogenase.

Pyrimidine biosynthesis presents an attractive drug target in malaria parasites due to the absence of a pyrimidine salvage pathway. A set of compounds designed to inhibit the Plasmodium falciparum pyrimidine biosynthetic enzyme dihydroorotate dehydrogenase (PfDHODH) was synthesized. PfDHODH-specific inhibitors with low nanomolar binding affinities were identified that bind in the N-terminal hydrophobic channel of dihydroorotate dehydrogenase, the presumed site of ubiquinone binding during oxidation of dihydroorotate to orotate. These compounds also prevented growth of cultured parasites at low micromolar concentrations. Models that suggest the mode of inhibitor binding is based on shape complementarity, matching hydrophobic regions of inhibitor and enzyme, and interaction of inhibitors with amino acid residues F188, H185, and R265 are supported by mutagenesis data. These results further highlight PfDHODH as a promising new target for chemotherapeutic intervention in prevention of malaria and provide better understanding of the factors that determine specificity over human dihydroorotate dehydrogenase.

metaSHARK: software for automated metabolic network prediction from DNA sequence and its application to the genomes of Plasmodium falciparum and Eimeria tenella.

The metabolic SearcH And Reconstruction Kit (metaSHARK) is a new fully automated software package for the detection of enzyme-encoding genes within unannotated genome data and their visualization in the context of the surrounding metabolic network. The gene detection package (SHARKhunt) runs on a Linux system and requires only a set of raw DNA sequences (genomic, expressed sequence tag and/or genome survey sequence) as input. Its output may be uploaded to our web-based visualization tool (SHARKview) for exploring and comparing data from different organisms. We first demonstrate the utility of the software by comparing its results for the raw Plasmodium falciparum genome with the manual annotations available at the PlasmoDB and PlasmoCyc websites. We then apply SHARKhunt to the unannotated genome sequences of the coccidian parasite Eimeria tenella and observe that, at an E-value cut-off of 10(-20), our software makes 142 additional assertions of enzymatic function compared with a recent annotation package working with translated open reading frame sequences. The ability of the software to cope with low levels of sequence coverage is investigated by analyzing assemblies of the E.tenella genome at estimated coverages from 0.5x to 7.5x. Lastly, as an example of how metaSHARK can be used to evaluate the genomic evidence for specific metabolic pathways, we present a study of coenzyme A biosynthesis in P.falciparum and E.tenella.