Predicting gene essentiality in Caenorhabditis elegans by feature engineering and machine-learning

Defining genes that are essential for life has major implications for understanding critical biological processes and mechanisms. Although essential genes have been identified and characterised experimentally using functional genomic tools, it is challenging to predict with confidence such genes from molecular and phenomic data sets using computational methods. Using extensive data sets available for the model organism Caenorhabditis elegans, we constructed here a machine-learning (ML)-based workflow for the prediction of essential genes on a genome-wide scale. We identified strong predictors for such genes and showed that trained ML models consistently achieve highly-accurate classifications. Complementary analyses revealed an association between essential genes and chromosomal location. Our findings reveal that essential genes in C. elegans tend to be located in or near the centre of autosomal chromosomes; are positively correlated with low single nucleotide polymorphim (SNP) densities and epigenetic markers in promoter regions; are involved in protein and nucleotide processing; are transcribed in most cells; are enriched in reproductive tissues or are targets for small RNAs bound to the argonaut CSR-1. Based on these results, we hypothesise an interplay between epigenetic markers and small RNA pathways in the germline, with transcription-based memory; this hypothesis warrants testing. From a technical perspective, further work is needed to evaluate whether the present ML-based approach will be applicable to other metazoans (including Drosophila melanogaster) for which comprehensive data set (i.e. genomic, transcriptomic, proteomic, variomic, epigenetic and phenomic) are available

Identification of G protein-coupled receptors in Schistosoma haematobium and S. mansoni by comparative genomics

BACKGROUND: Schistosomiasis is a parasitic disease affecting ~200 million people worldwide. Schistosoma haematobium and S. mansoni are two relatively closely related schistosomes (blood flukes), and the causative agents of urogenital and hepatointestinal schistosomiasis, respectively. The availability of genomic, transcriptomic and proteomic data sets for these two schistosomes now provides unprecedented opportunities to explore their biology, host interactions and schistosomiasis at the molecular level. A particularly important group of molecules involved in a range of biological and developmental processes in schistosomes and other parasites are the G protein-coupled receptors (GPCRs). Although GPCRs have been studied in schistosomes, there has been no detailed comparison of these receptors between closely related species. Here, using a genomic-bioinformatic approach, we identified and characterised key GPCRs in S. haematobium and S. mansoni (two closely related species of schistosome). METHODS: Using a Hidden Markov Model (HMM) and Support Vector Machine (SVM)-based pipeline, we classified and sub-classified GPCRs of S. haematobium and S. mansoni, combined with phylogenetic and transcription analyses. RESULTS: We identified and classified classes A, B, C and F as well as an unclassified group of GPCRs encoded in the genomes of S. haematobium and S. mansoni. In addition, we characterised ligand-specific subclasses (i.e. amine, peptide, opsin and orphan) within class A (rhodopsin-like). CONCLUSIONS: Most GPCRs shared a high degree of similarity and conservation, except for members of a particular clade (designated SmGPR), which appear to have diverged between S. haematobium and S. mansoni and might explain, to some extent, some of the underlying biological differences between these two schistosomes. The present set of annotated GPCRs provides a basis for future functional genomic studies of cellular GPCR-mediated signal transduction and a resource for future drug discovery efforts in schistosomes

Analysis of Haemonchus embryos at single cell resolution identifies two eukaryotic elongation factors as intervention target candidates

Advances in single cell technologies are allowing investigations of a wide range of biological processes and pathways in animals, such as the multicellular model organism Caenorhabditis elegans – a free-living nematode. However, there has been limited application of such technology to related parasitic nematodes which cause major diseases of humans and animals worldwide. With no vaccines against the vast majority of parasitic nematodes and treatment failures due to drug resistance or inefficacy, new intervention targets are urgently needed, preferably informed by a deep understanding of these nematodes’ cellular and molecular biology – which is presently lacking for most worms. Here, we created the first single cell atlas for an early developmental stage of Haemonchus contortus – a highly pathogenic, C. elegans-related parasitic nematode. We obtained and curated RNA sequence (snRNA-seq) data from single nuclei from embryonating eggs of H. contortus (150,000 droplets), and selected high-quality transcriptomic data for > 14,000 single nuclei for analysis, and identified 19 distinct clusters of cells. Guided by comparative analyses with C. elegans, we were able to reproducibly assign seven cell clusters to body wall muscle, hypodermis, neuronal, intestinal or seam cells, and identified eight genes that were transcribed in all cell clusters/types, three of which were inferred to be essential in H. contortus. Two of these genes (i.e. Hc-eef-1A and Hc-eef1G), coding for eukaryotic elongation factors (called Hc-eEF1A and Hc-eEF1G), were also demonstrated to be transcribed and expressed in all key developmental stages of H. contortus. Together with these findings, sequence- and structure-based comparative analyses indicated the potential of Hc-eEF1A and/or Hc-eEF1G as intervention targets within the protein biosynthesis machinery of H. contortus. Future work will focus on single cell studies of all key developmental stages and tissues of H. contortus, and on evaluating the suitability of the two elongation factor proteins as drug targets in H. contortus and related nematodes, with a view to finding new nematocidal drug candidates