OGEE v3: Online GEne Essentiality database with increased coverage of organisms and human cell lines

Bacterial genome reduction for optimal chassis of synthetic biology: a review

Bacteria with streamlined genomes, that harbor full functional genes for essential metabolic networks, are able to synthesize the desired products more effectively and thus have advantages as production platforms in industrial applications. To obtain streamlined chassis genomes, a large amount of effort has been made to reduce existing bacterial genomes. This work falls into two categories: rational and random reduction. The identification of essential gene sets and the emergence of various genome-deletion techniques have greatly promoted genome reduction in many bacteria over the past few decades. Some of the constructed genomes possessed desirable properties for industrial applications, such as: increased genome stability, transformation capacity, cell growth, and biomaterial productivity. The decreased growth and perturbations in physiological phenotype of some genome-reduced strains may limit their applications as optimized cell factories. This review presents an assessment of the advancements made to date in bacterial genome reduction to construct optimal chassis for synthetic biology, including: the identification of essential gene sets, the genome-deletion techniques, the properties and industrial applications of artificially streamlined genomes, the obstacles encountered in constructing reduced genomes, and the future perspectives.

DTYMK is an essential gene in mice and heterozygosity does not cause neural tube defects

Regulation of nucleotide biosynthesis is necessary for maintaining cellular processes including DNA replication and repair. A key enzyme in this process is deoxythymidylate kinase (dTYMK), which catalyzes the initial step in the production of dTTP from dTMP. This gene constitutes the first merged step of dTTP synthesis from the de novo and salvage pathways which regulate dTMP biosynthesis. Decreased de novo dTMP biosynthesis causes dysregulated dTTP:dUTP pools, and leads to increased uracil in DNA and neural tube closure defect (NTD) development in mice. The goal of this research was to investigate if dTYMK, the downstream enzyme in dTTP production, is an essential gene in mice and if impairments in dTYMK play a causal role in development including NTD pathology in mice. Dtymk+/- C57BL/6J females were weaned onto either a control, excess folic acid, or folic acid deficient diet and timed breeding was performed after 8 weeks on diet. The offspring were analyzed for NTDs and other reproductive outcomes at embryonic day 12.5 (E12.5). Dtymk-/- mice were confirmed to be embryonic lethal before E12.5, and Dtymk+/- mice on all three experimental diets did not show the presence of open neural tube defects, spina bifida or exencephaly. However, the expression of dTYMK in Dtymk+/- mouse embryos was confirmed to be decreased by approximately 3-fold compared to Dtymk+/+ embryos. Although dTYMK was demonstrated to be an essential gene in mice and is required for the regulation of nucleotide pools in vitro, there was no evidence of increased risk of NTDs because of a reduction in expression of this enzyme during embryonic development. It is possible that a further reduction in expression may be required to see developmental anomalies in C57BL/6J mice.

VarPPUD: Variant post prioritization developed for undiagnosed genetic disorders

Rare and ultra-rare genetic conditions are estimated to impact nearly 1 in 17 people worldwide, yet accurately pinpointing the diagnostic variants underlying each of these conditions remains a formidable challenge. Because comprehensive, in vivo functional assessment of all possible genetic variants is infeasible, clinicians instead consider in silico variant pathogenicity predictions to distinguish plausibly disease-causing from benign variants across the genome. However, in the most difficult undiagnosed cases, such as those accepted to the Undiagnosed Diseases Network (UDN), existing pathogenicity predictions cannot reliably discern true etiological variant(s) from other deleterious candidate variants that were prioritized through N-of-1 efforts. Pinpointing the disease-causing variant from a pool of plausible candidates remains a largely manual effort requiring extensive clinical workups, functional and experimental assays, and eventual identification of genotype- and phenotype-matched individuals. Here, we introduce VarPPUD, a tool trained on prioritized variants from UDN cases, that leverages gene-, amino acid-, and nucleotide-level features to discern pathogenic variants from other deleterious variants that are unlikely to be confirmed as disease relevant. VarPPUD achieves a cross-validated accuracy of 79.3% and precision of 77.5% on a held-out subset of uniquely challenging UDN cases, respectively representing an average 18.6% and 23.4% improvement over nine traditional pathogenicity prediction approaches on this task. We validate VarPPUD's ability to discriminate likely from unlikely pathogenic variants on synthetic, GAN-generated candidate variants as well. Finally, we show how VarPPUD can be probed to evaluate each input feature's importance and contribution toward prediction-an essential step toward understanding the distinct characteristics of newly-uncovered disease-causing variants.

Significance Statement

Patients with chronic, undiagnosed and underdiagnosed genetic conditions often endure expensive and excruciating years-long diagnostic odysseys without clear results. In many instances, clinical genome sequencing of patients and their family members fails to reveal known disease-causing variants, although compelling variants of uncertain significance are frequently encountered. Existing computational tools struggle to reliably differentiate truly disease-causing variants from other plausible candidate variants within these prioritized sets. Consequently, the confirmation of disease-causing variants often necessitates extensive experimental follow-up, including studies in model organisms and identification of other similarly presenting genotype-matched individuals, a process that can extend for several years. Here, we present VarPPUD, a tool trained specifically to distinguish likely from unlikely to be confirmed pathogenic variants that were prioritized across cases in the Undiagnosed Diseases Network. By evaluating the importance and impact of different input feature values on prediction, we gain deeper insights into the distinctive attributes of difficult-to-identify diagnostic variants. For patients who remain undiagnosed following comprehensive whole genome sequencing, our new method VarPPUD may reveal pathogenic variants amid a pool of candidate variants, thereby advancing diagnostic efforts where progress has otherwise stalled.

Machine learning methods for predicting guide RNA effects in CRISPR epigenome editing experiments

CRISPR epigenomic editing technologies enable functional interrogation of non-coding elements. However, current computational methods for guide RNA (gRNA) design do not effectively predict the power potential, molecular and cellular impact to optimize for efficient gRNAs, which are crucial for successful applications of these technologies. We present "launch-dCas9" (machine LeArning based UNified CompreHensive framework for CRISPR-dCas9) to predict gRNA impact from multiple perspectives, including cell fitness, wildtype abundance (gauging power potential), and gene expression in single cells. Our launchdCas9, built and evaluated using experiments involving >1 million gRNAs targeted across the human genome, demonstrates relatively high prediction accuracy (AUC up to 0.81) and generalizes across cell lines. Method-prioritized top gRNA(s) are 4.6-fold more likely to exert effects, compared to other gRNAs in the same cis-regulatory region. Furthermore, launchdCas9 identifies the most critical sequence-related features and functional annotations from >40 features considered. Our results establish launch-dCas9 as a promising approach to design gRNAs for CRISPR epigenomic experiments.

Methylation entropy landscape of Chinese long-lived individuals reveals lower epigenetic noise related to human healthy aging

The transition from ordered to noisy is a significant epigenetic signature of aging and age-related disease. As a paradigm of healthy human aging and longevity, long-lived individuals (LLI, >90 years old) may possess characteristic strategies in coping with the disordered epigenetic regulation. In this study, we constructed high-resolution blood epigenetic noise landscapes for this cohort by a methylation entropy (ME) method using whole genome bisulfite sequencing (WGBS). Although a universal increase in global ME occurred with chronological age in general control samples, this trend was suppressed in LLIs. Importantly, we identified 38,923 genomic regions with LLI-specific lower ME (LLI-specific lower entropy regions, for short, LLI-specific LERs). These regions were overrepresented in promoters, which likely function in transcriptional noise suppression. Genes associated with LLI-specific LERs have a considerable impact on SNP-based heritability of some aging-related disorders (e.g., asthma and stroke). Furthermore, neutrophil was identified as the primary cell type sustaining LLI-specific LERs. Our results highlight the stability of epigenetic order in promoters of genes involved with aging and age-related disorders within LLI epigenomes. This unique epigenetic feature reveals a previously unknown role of epigenetic order maintenance in specific genomic regions of LLIs, which helps open a new avenue on the epigenetic regulation mechanism in human healthy aging and longevity.

Comparative structure activity and target exploration of 1,2-diphenylethynes in Haemonchus contortus and Caenorhabditis elegans

Infections and diseases caused by parasitic nematodes have a major adverse impact on the health and productivity of animals and humans worldwide. The control of these parasites often relies heavily on the treatment with commercially available chemical compounds (anthelmintics). However, the excessive or uncontrolled use of these compounds in livestock animals has led to major challenges linked to drug resistance in nematodes. Therefore, there is a need to develop new anthelmintics with novel mechanism(s) of action. Recently, we identified a small molecule, designated UMW-9729, with nematocidal activity against the free-living model organism Caenorhabditis elegans. Here, we evaluated UMW-9729's potential as an anthelmintic in a structure-activity relationship (SAR) study in C. elegans and the highly pathogenic, blood-feeding Haemonchus contortus (barber's pole worm), and explored the compound-target relationship using thermal proteome profiling (TPP). First, we synthesised and tested 25 analogues of UMW-9729 for their nematocidal activity in both H. contortus (larvae and adults) and C. elegans (young adults), establishing a preliminary nematocidal pharmacophore for both species. We identified several compounds with marked activity against either H. contortus or C. elegans which had greater efficacy than UMW-9729, and found a significant divergence in compound bioactivity between these two nematode species. We also identified a UMW-9729 analogue, designated 25, that moderately inhibited the motility of adult female H. contortus in vitro. Subsequently, we inferred three H. contortus proteins (HCON_00134350, HCON_00021470 and HCON_00099760) and five C. elegans proteins (F30A10.9, F15B9.8, B0361.6, DNC-4 and UNC-11) that interacted directly with UMW-9729; however, no conserved protein target was shared between the two nematode species. Future work aims to extend the SAR investigation in these and other parasitic nematode species, and validate individual proteins identified here as possible targets of UMW-9729. Overall, the present study evaluates this anthelmintic candidate and highlights some challenges associated with early anthelmintic investigation.

ECDEP: identifying essential proteins based on evolutionary community discovery and subcellular localization

BACKGROUND

In cellular activities, essential proteins play a vital role and are instrumental in comprehending fundamental biological necessities and identifying pathogenic genes. Current deep learning approaches for predicting essential proteins underutilize the potential of gene expression data and are inadequate for the exploration of dynamic networks with limited evaluation across diverse species.

RESULTS

We introduce ECDEP, an essential protein identification model based on evolutionary community discovery. ECDEP integrates temporal gene expression data with a protein-protein interaction (PPI) network and employs the 3-Sigma rule to eliminate outliers at each time point, constructing a dynamic network. Next, we utilize edge birth and death information to establish an interaction streaming source to feed into the evolutionary community discovery algorithm and then identify overlapping communities during the evolution of the dynamic network. SVM recursive feature elimination (RFE) is applied to extract the most informative communities, which are combined with subcellular localization data for classification predictions. We assess the performance of ECDEP by comparing it against ten centrality methods, four shallow machine learning methods with RFE, and two deep learning methods that incorporate multiple biological data sources on Saccharomyces. Cerevisiae (S. cerevisiae), Homo sapiens (H. sapiens), Mus musculus, and Caenorhabditis elegans. ECDEP achieves an AP value of 0.86 on the H. sapiens dataset and the contribution ratio of community features in classification reaches 0.54 on the S. cerevisiae (Krogan) dataset.

CONCLUSIONS

Our proposed method adeptly integrates network dynamics and yields outstanding results across various datasets. Furthermore, the incorporation of evolutionary community discovery algorithms amplifies the capacity of gene expression data in classification.

Lethal phenotypes in Mendelian disorders

Essential genes are those whose function is required for cell proliferation and/or organism survival. A gene's intolerance to loss-of-function can be allocated within a spectrum, as opposed to being considered a binary feature, since this function might be essential at different stages of development, genetic backgrounds or other contexts. Existing resources that collect and characterise the essentiality status of genes are based on either proliferation assessment in human cell lines, embryonic and postnatal viability evaluation in different model organisms, and gene metrics such as intolerance to variation scores derived from human population sequencing studies. There are also several repositories available that document phenotypic annotations for rare disorders in humans such as the Online Mendelian Inheritance in Man (OMIM) and the Human Phenotype Ontology (HPO) knowledgebases. This raises the prospect of being able to use clinical data, including lethality as the most severe phenotypic manifestation, to further our characterisation of gene essentiality. Here we queried OMIM for terms related to lethality and classified all Mendelian genes into categories, according to the earliest age of death recorded for the associated disorders, from prenatal death to no reports of premature death. To showcase this curated catalogue of human essential genes, we developed the Lethal Phenotypes Portal (https://lethalphenotypes.research.its.qmul.ac.uk), where we also explore the relationships between these lethality categories, constraint metrics and viability in cell lines and mouse. Further analysis of the genes in these categories reveals differences in the mode of inheritance of the associated disorders, physiological systems affected and disease class. We highlight how the phenotypic similarity between genes in the same lethality category combined with gene family/group information can be used for novel disease gene discovery. Finally, we explore the overlaps and discrepancies between the lethal phenotypes observed in mouse and human and discuss potential explanations that include differences in transcriptional regulation, functional compensation and molecular disease mechanisms. We anticipate that this resource will aid clinicians in the diagnosis of early lethal conditions and assist researchers in investigating the properties that make these genes essential for human development.

Untangling the Context-Specificity of Essential Genes by Means of Machine Learning: A Constructive Experience

Gene essentiality is a genetic concept crucial for a comprehensive understanding of life and evolution. In the last decade, many essential genes (EGs) have been determined using different experimental and computational approaches, and this information has been used to reduce the genomes of model organisms. A growing amount of evidence highlights that essentiality is a property that depends on the context. Because of their importance in vital biological processes, recognising context-specific EGs (csEGs) could help for identifying new potential pharmacological targets and to improve precision therapeutics. Since most of the computational procedures proposed to identify and predict EGs neglect their context-specificity, we focused on this aspect, providing a theoretical and experimental overview of the literature, data and computational methods dedicated to recognising csEGs. To this end, we adapted existing computational methods to exploit a specific context (the kidney tissue) and experimented with four different prediction methods using the labels provided by four different identification approaches. The considerations derived from the analysis of the obtained results, confirmed and validated also by further experiments for a different tissue context, provide the reader with guidance on exploiting existing tools for achieving csEGs identification and prediction.

A comprehensive atlas of epigenetic regulators reveals tissue-specific epigenetic regulation patterns

Epigenetic machinery contributes to gene regulation in eukaryotic species. However, the machinery including more than 600 epigenetic regulator (ER) genes responsible for reading, writing, and erasing histone modifications and DNA modifications remains largely uncharacterized across species. We compile a comprehensive list of ERs based on an evolutionary analysis across 23 species, which is the most comprehensive ER list in various species until recently. We further perform comparative transcriptomic analyses across different tissues in humans, mice, as well as other amniote species. We observe a consistent tissue-of-origin expression specificity pattern of duplicated ER genes across species and suggest links between expression specificity and ER gene evolution as well as ER function. Additional analyses further suggest that ER duplication can generate tissue-specific ER genes with the same epigenetic substrates, which may be closely related to their regulatory specificity in tissue development. Our work can serve as a foundation to better comprehend the tissue-specific expression patterns of ER genes from an evolutionary perspective and also the functional implications of ERs in tissue-specific epigenetic regulation.

'Bingo'-a large language model- and graph neural network-based workflow for the prediction of essential genes from protein data

The identification and characterization of essential genes are central to our understanding of the core biological functions in eukaryotic organisms, and has important implications for the treatment of diseases caused by, for example, cancers and pathogens. Given the major constraints in testing the functions of genes of many organisms in the laboratory, due to the absence of in vitro cultures and/or gene perturbation assays for most metazoan species, there has been a need to develop in silico tools for the accurate prediction or inference of essential genes to underpin systems biological investigations. Major advances in machine learning approaches provide unprecedented opportunities to overcome these limitations and accelerate the discovery of essential genes on a genome-wide scale. Here, we developed and evaluated a large language model- and graph neural network (LLM-GNN)-based approach, called 'Bingo', to predict essential protein-coding genes in the metazoan model organisms Caenorhabditis elegans and Drosophila melanogaster as well as in Mus musculus and Homo sapiens (a HepG2 cell line) by integrating LLM and GNNs with adversarial training. Bingo predicts essential genes under two 'zero-shot' scenarios with transfer learning, showing promise to compensate for a lack of high-quality genomic and proteomic data for non-model organisms. In addition, the attention mechanisms and GNNExplainer were employed to manifest the functional sites and structural domain with most contribution to essentiality. In conclusion, Bingo provides the prospect of being able to accurately infer the essential genes of little- or under-studied organisms of interest, and provides a biological explanation for gene essentiality.

A genome-scale metabolic model of parasitic whipworm

Genome-scale metabolic models are widely used to enhance our understanding of metabolic features of organisms, host-pathogen interactions and to identify therapeutics for diseases. Here we present iTMU798, the genome-scale metabolic model of the mouse whipworm Trichuris muris. The model demonstrates the metabolic features of T. muris and allows the prediction of metabolic steps essential for its survival. Specifically, that Thioredoxin Reductase (TrxR) enzyme is essential, a prediction we validate in vitro with the drug auranofin. Furthermore, our observation that the T. muris genome lacks gsr-1 encoding Glutathione Reductase (GR) but has GR activity that can be inhibited by auranofin indicates a mechanism for the reduction of glutathione by the TrxR enzyme in T. muris. In addition, iTMU798 predicts seven essential amino acids that cannot be synthesised by T. muris, a prediction we validate for the amino acid tryptophan. Overall, iTMU798 is as a powerful tool to study not only the T. muris metabolism but also other Trichuris spp. in understanding host parasite interactions and the rationale design of new intervention strategies.

Essential genes: a cross-species perspective

Protein coding genes exhibit different degrees of intolerance to loss-of-function variation. The most intolerant genes, whose function is essential for cell or/and organism survival, inform on fundamental biological processes related to cell proliferation and organism development and provide a window on the molecular mechanisms of human disease. Here we present a brief overview of the resources and knowledge gathered around gene essentiality, from cancer cell lines to model organisms to human development. We outline the implications of using different sources of evidence and definitions to determine which genes are essential and highlight how information on the essentiality status of a gene can inform novel disease gene discovery and therapeutic target identification.

Machine learning model for anti-cancer drug combinations: Analysis, prediction, and validation

Drug combination therapy is a highly effective approach for enhancing the therapeutic efficacy of anti-cancer drugs and overcoming drug resistance. However, the innumerable possible drug combinations make it impractical to screen all synergistic drug pairs. Moreover, biological insights into synergistic drug pairs are still lacking. To address this challenge, we systematically analyzed drug combination datasets curated from multiple databases to identify drug pairs more likely to show synergy. We classified drug pairs based on their MoA and discovered that 110 MoA pairs were significantly enriched in synergy in at least one type of cancer. To improve the accuracy of predicting synergistic effects of drug pairs, we developed a suite of machine learning models that achieve better predictive performance. Unlike most previous methods that were rarely validated by wet-lab experiments, our models were validated using two-dimensional cell lines and three-dimensional tumor slice culture (3D-TSC) models, implying their practical utility. Our prediction and validation results indicated that the combination of the RTK inhibitors Lapatinib and Pazopanib exhibited a strong therapeutic effect in breast cancer by blocking the downstream PI3K/AKT/mTOR signaling pathway. Furthermore, we incorporated molecular features to identify potential biomarkers for synergistic drug pairs, and almost all potential biomarkers found connections between drug targets and corresponding molecular features using protein-protein interaction network. Overall, this study provides valuable insights to complement and guide rational efforts to develop drug combination treatments.

A new metabolic model of Drosophila melanogaster and the integrative analysis of Parkinson's disease

High conservation of the disease-associated genes between flies and humans facilitates the common use of Drosophila melanogaster to study metabolic disorders under controlled laboratory conditions. However, metabolic modeling studies are highly limited for this organism. We here report a comprehensively curated genome-scale metabolic network model of Drosophila using an orthology-based approach. The gene coverage and metabolic information of the draft model derived from a reference human model were expanded via Drosophila-specific KEGG and MetaCyc databases, with several curation steps to avoid metabolic redundancy and stoichiometric inconsistency. Furthermore, we performed literature-based curations to improve gene-reaction associations, subcellular metabolite locations, and various metabolic pathways. The performance of the resulting Drosophila model (8,230 reactions, 6,990 metabolites, and 2,388 genes), iDrosophila1 (https://github.com/SysBioGTU/iDrosophila), was assessed using flux balance analysis in comparison with the other currently available fly models leading to superior or comparable results. We also evaluated the transcriptome-based prediction capacity of iDrosophila1, where differential metabolic pathways during Parkinson's disease could be successfully elucidated. Overall, iDrosophila1 is promising to investigate system-level metabolic alterations in response to genetic and environmental perturbations.

De novo mutations disturb early brain development more frequently than common variants in schizophrenia

Investigating functional, temporal, and cell-type expression features of mutations is important for understanding a complex disease. Here, we collected and analyzed common variants and de novo mutations (DNMs) in schizophrenia (SCZ). We collected 2,636 missense and loss-of-function (LoF) DNMs in 2,263 genes across 3,477 SCZ patients (SCZ-DNMs). We curated three gene lists: (a) SCZ-neuroGenes (159 genes), which are intolerant to LoF and missense DNMs and are neurologically important, (b) SCZ-moduleGenes (52 genes), which were derived from network analyses of SCZ-DNMs, and (c) SCZ-commonGenes (120 genes) from a recent GWAS as reference. To compare temporal gene expression, we used the BrainSpan dataset. We defined a fetal effect score (FES) to quantify the involvement of each gene in prenatal brain development. We further employed the specificity indexes (SIs) to evaluate cell-type expression specificity from single-cell expression data in cerebral cortices of humans and mice. Compared with SCZ-commonGenes, SCZ-neuroGenes and SCZ-moduleGenes were highly expressed in the prenatal stage, had higher FESs, and had higher SIs in fetal replicating cells and undifferentiated cell types. Our results suggested that gene expression patterns in specific cell types in early fetal stages might have impacts on the risk of SCZ during adulthood.

Methylation of nonessential genes in cutaneous melanoma - Rule Out hypothesis

Differential methylation plays an important role in melanoma development and is associated with survival, progression and response to treatment. However, the mechanisms by which methylation promotes melanoma development are poorly understood. The traditional explanation of selective advantage provided by differential methylation postulates that hypermethylation of regulatory 5'-cytosine-phosphate-guanine-3' dinucleotides (CpGs) downregulates the expression of tumor suppressor genes and therefore promotes tumorigenesis. We believe that other (not necessarily alternative) explanations of the selective advantages of methylation are also possible. Here, we hypothesize that melanoma cells use methylation to shut down transcription of nonessential genes - those not required for cell survival and proliferation. Suppression of nonessential genes allows tumor cells to be more efficient in terms of energy and resource usage, providing them with a selective advantage over the tumor cells that transcribe and subsequently translate genes they do not need. We named the hypothesis the Rule Out (RO) hypothesis. The RO hypothesis predicts higher methylation of CpGs located in regulatory regions (CpG islands) of nonessential genes. It also predicts the higher methylation of regulatory CpGs linked to nonessential genes in melanomas compared to nevi and lower expression of nonessential genes in malignant (derived from melanoma) versus normal (derived from nonaffected skin) melanocytes. The analyses conducted using in-house and publicly available data found that all predictions derived from the RO hypothesis hold, providing observational support for the hypothesis.

Recent advances in genetic tools for engineering probiotic lactic acid bacteria

Synthetic biology has grown exponentially in the last few years, with a variety of biological applications. One of the emerging applications of synthetic biology is to exploit the link between microorganisms, biologics, and human health. To exploit this link, it is critical to select effective synthetic biology tools for use in appropriate microorganisms that would address unmet needs in human health through the development of new game-changing applications and by complementing existing technological capabilities. Lactic acid bacteria (LAB) are considered appropriate chassis organisms that can be genetically engineered for therapeutic and industrial applications. Here, we have reviewed comprehensively various synthetic biology techniques for engineering probiotic LAB strains, such as clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 mediated genome editing, homologous recombination, and recombineering. In addition, we also discussed heterologous protein expression systems used in engineering probiotic LAB. By combining computational biology with genetic engineering, there is a lot of potential to develop next-generation synthetic LAB with capabilities to address bottlenecks in industrial scale-up and complex biologics production. Recently, we started working on Lactochassis project where we aim to develop next generation synthetic LAB for biomedical application.

The occurrence, characteristics, and adaptation of A-to-I RNA editing in bacteria: A review

A-to-I RNA editing is a very important post-transcriptional modification or co-transcriptional modification that creates isoforms and increases the diversity of proteins. In this process, adenosine (A) in RNA molecules is hydrolyzed and deaminated into inosine (I). It is well known that ADAR (adenosine deaminase acting on RNA)-dependent A-to-I mRNA editing is widespread in animals. Next, the discovery of A-to-I mRNA editing was mediated by TadA (tRNA-specific adenosine deaminase) in Escherichia coli which is ADAR-independent event. Previously, the editing event S128P on the flagellar structural protein FliC enhanced the bacterial tolerance to oxidative stress in Xoc. In addition, the editing events T408A on the enterobactin iron receptor protein XfeA act as switches by controlling the uptake of Fe³⁺ in response to the concentration of iron in the environment. Even though bacteria have fewer editing events, the great majority of those that are currently preserved have adaptive benefits. Interestingly, it was found that a TadA-independent A-to-I RNA editing event T408A occurred on xfeA, indicating that there may be other new enzymes that perform a function like TadA. Here, we review recent advances in the characteristics, functions, and adaptations of editing in bacteria.

Gene co-expression analyses of health(span) across multiple species

Health(span)-related gene clusters/modules were recently identified based on knowledge about the cross-species genetic basis of health, to interpret transcriptomic datasets describing health-related interventions. However, the cross-species comparison of health-related observations reveals a lot of heterogeneity, not least due to widely varying health(span) definitions and study designs, posing a challenge for the exploration of conserved healthspan modules and, specifically, their transfer across species. To improve the identification and exploration of conserved/transferable healthspan modules, here we apply an established workflow based on gene co-expression network analyses employing GEO/ArrayExpress data for human and animal models, and perform a comprehensive meta-study of the resulting modules related to health(span), yielding a small set of literature backed health(span) candidate genes. For each experiment, WGCNA (weighted gene correlation network analysis) was used to infer modules of genes which correlate in their expression with a 'health phenotype score' and to determine the most-connected (hub) genes (and their interactions) for each such module. After mapping these hub genes to their human orthologs, 12 health(span) genes were identified in at least two species (ACTN3, ANK1, MRPL18, MYL1, PAXIP1, PPP1CA, SCN3B, SDCBP, SKIV2L, TUBG1, TYROBP, WIPF1), for which enrichment analysis by g:profiler found an association with actin filament-based movement and associated organelles, as well as muscular structures. We conclude that a meta-study of hub genes from co-expression network analyses for the complex phenotype health(span), across multiple species, can yield molecular-mechanistic insights and can direct experimentalists to further investigate the contribution of individual genes and their interactions to health(span).

Rule-Based Pruning and In Silico Identification of Essential Proteins in Yeast PPIN

Proteins are vital for the significant cellular activities of living organisms. However, not all of them are essential. Identifying essential proteins through different biological experiments is relatively more laborious and time-consuming than the computational approaches used in recent times. However, practical implementation of conventional scientific methods sometimes becomes challenging due to poor performance impact in specific scenarios. Thus, more developed and efficient computational prediction models are required for essential protein identification. An effective methodology is proposed in this research, capable of predicting essential proteins in a refined yeast protein–protein interaction network (PPIN). The rule-based refinement is done using protein complex and local interaction density information derived from the neighborhood properties of proteins in the network. Identification and pruning of non-essential proteins are equally crucial here. In the initial phase, careful assessment is performed by applying node and edge weights to identify and discard the non-essential proteins from the interaction network. Three cut-off levels are considered for each node and edge weight for pruning the non-essential proteins. Once the PPIN has been filtered out, the second phase starts with two centralities-based approaches: (1) local interaction density (LID) and (2) local interaction density with protein complex (LIDC), which are successively implemented to identify the essential proteins in the yeast PPIN. Our proposed methodology achieves better performance in comparison to the existing state-of-the-art techniques.

A deep learning framework for identifying essential proteins based on multiple biological information

Background

Essential Proteins are demonstrated to exert vital functions on cellular processes and are indispensable for the survival and reproduction of the organism. Traditional centrality methods perform poorly on complex protein–protein interaction (PPI) networks. Machine learning approaches based on high-throughput data lack the exploitation of the temporal and spatial dimensions of biological information.

Results

We put forward a deep learning framework to predict essential proteins by integrating features obtained from the PPI network, subcellular localization, and gene expression profiles. In our model, the node2vec method is applied to learn continuous feature representations for proteins in the PPI network, which capture the diversity of connectivity patterns in the network. The concept of depthwise separable convolution is employed on gene expression profiles to extract properties and observe the trends of gene expression over time under different experimental conditions. Subcellular localization information is mapped into a long one-dimensional vector to capture its characteristics. Additionally, we use a sampling method to mitigate the impact of imbalanced learning when training the model. With experiments carried out on the data of Saccharomyces cerevisiae, results show that our model outperforms traditional centrality methods and machine learning methods. Likewise, the comparative experiments have manifested that our process of various biological information is preferable.

Conclusions

Our proposed deep learning framework effectively identifies essential proteins by integrating multiple biological data, proving a broader selection of subcellular localization information significantly improves the results of prediction and depthwise separable convolution implemented on gene expression profiles enhances the performance.

Computational estimation of quality and clinical relevance of cancer cell lines

Immortal cancer cell lines (CCLs) are the most widely used system for investigating cancer biology and for the preclinical development of oncology therapies. Pharmacogenomic and genome‐wide editing screenings have facilitated the discovery of clinically relevant gene–drug interactions and novel therapeutic targets via large panels of extensively characterised CCLs. However, tailoring pharmacological strategies in a precision medicine context requires bridging the existing gaps between tumours and in vitro models. Indeed, intrinsic limitations of CCLs such as misidentification, the absence of tumour microenvironment and genetic drift have highlighted the need to identify the most faithful CCLs for each primary tumour while addressing their heterogeneity, with the development of new models where necessary. Here, we discuss the most significant limitations of CCLs in representing patient features, and we review computational methods aiming at systematically evaluating the suitability of CCLs as tumour proxies and identifying the best patient representative in vitro models. Additionally, we provide an overview of the applications of these methods to more complex models and discuss future machine‐learning‐based directions that could resolve some of the arising discrepancies.

dbEssLnc: A manually curated database of human and mouse essential lncRNA genes

Long non-coding RNAs (lncRNAs) play important roles in many biological processes. Knocking out or knocking down some lncRNAs will lead to lethality or infertility. These lncRNAs are called essential lncRNAs. Knowledges of essential lncRNAs are important in establishing minimal genomes of living cells, developing drug therapies and early diagnostic approaches for complex diseases. However, existing databases focus on collecting essential coding genes. Essential non-coding gene records are rare in existing databases. A comprehensive collection of essential non-coding genes, particularly essential lncRNA genes, is demanded. We manually curated 207 essential lncRNAs from literatures for establishing a database on essential lncRNAs, which is named as dbEssLnc (Database of essential lncRNAs). The dbEssLnc database has a web-based user-friendly interface for the users to browse, to search, to visualize and to blast search records in the database. The dbEssLnc database is freely accessible at https://esslnc.pufengdu.org. All data and source codes for mirroring the dbEssLnc database have been deposited in GitHub (https://github.com/yyZhang14/dbEssLnc).

Network cartographs for interpretable visualizations

Networks offer an intuitive visual representation of complex systems. Important network characteristics can often be recognized by eye and, in turn, patterns that stand out visually often have a meaningful interpretation. In conventional network layout algorithms, however, the precise determinants of a node's position within a layout are difficult to decipher and to control. Here we propose an approach for directly encoding arbitrary structural or functional network characteristics into node positions. We introduce a series of two- and three-dimensional layouts, benchmark their efficiency for model networks, and demonstrate their power for elucidating structure-to-function relationships in large-scale biological networks.

RBM22, a Key Player of Pre-mRNA Splicing and Gene Expression Regulation, Is Altered in Cancer

RNA-Binding Proteins (RBP) are very diverse and cover a large number of functions in the cells. This review focuses on RBM22, a gene encoding an RBP and belonging to the RNA-Binding Motif (RBM) family of genes. RBM22 presents a Zinc Finger like and a Zinc Finger domain, an RNA-Recognition Motif (RRM), and a Proline-Rich domain with a general structure suggesting a fusion of two yeast genes during evolution: Cwc2 and Ecm2. RBM22 is mainly involved in pre-mRNA splicing, playing the essential role of maintaining the conformation of the catalytic core of the spliceosome and acting as a bridge between the catalytic core and other essential protein components of the spliceosome. RBM22 is also involved in gene regulation, and is able to bind DNA, acting as a bona fide transcription factor on a large number of target genes. Undoubtedly due to its wide scope in the regulation of gene expression, RBM22 has been associated with several pathologies and, notably, with the aggressiveness of cancer cells and with the phenotype of a myelodysplastic syndrome. Mutations, enforced expression level, and haploinsufficiency of RBM22 gene are observed in those diseases. RBM22 could represent a potential therapeutic target in specific diseases, and, notably, in cancer.

Molecular Biology in Treatment Decision Processes-Neuro-Oncology Edition

Computational approaches including machine learning, deep learning, and artificial intelligence are growing in importance in all medical specialties as large data repositories are increasingly being optimised. Radiation oncology as a discipline is at the forefront of large-scale data acquisition and well positioned towards both the production and analysis of large-scale oncologic data with the potential for clinically driven endpoints and advancement of patient outcomes. Neuro-oncology is comprised of malignancies that often carry poor prognosis and significant neurological sequelae. The analysis of radiation therapy mediated treatment and the potential for computationally mediated analyses may lead to more precise therapy by employing large scale data. We analysed the state of the literature pertaining to large scale data, computational analysis, and the advancement of molecular biomarkers in neuro-oncology with emphasis on radiation oncology. We aimed to connect existing and evolving approaches to realistic avenues for clinical implementation focusing on low grade gliomas (LGG), high grade gliomas (HGG), management of the elderly patient with HGG, rare central nervous system tumors, craniospinal irradiation, and re-irradiation to examine how computational analysis and molecular science may synergistically drive advances in personalised radiation therapy (RT) and optimise patient outcomes.

Identifying essential genes across eukaryotes by machine learning

Identifying essential genes on a genome scale is resource intensive and has been performed for only a few eukaryotes. For less studied organisms essentiality might be predicted by gene homology. However, this approach cannot be applied to non-conserved genes. Additionally, divergent essentiality information is obtained from studying single cells or whole, multi-cellular organisms, and particularly when derived from human cell line screens and human population studies. We employed machine learning across six model eukaryotes and 60 381 genes, using 41 635 features derived from the sequence, gene function information and network topology. Within a leave-one-organism-out cross-validation, the classifiers showed high generalizability with an average accuracy close to 80% in the left-out species. As a case study, we applied the method to Tribolium castaneum and Bombyx mori and validated predictions experimentally yielding similar performances. Finally, using the classifier based on the studied model organisms enabled linking the essentiality information of human cell line screens and population studies.

Selective destabilization of polypeptides synthesized from NMD-targeted transcripts

The translation of mRNAs that contain a premature termination codon (PTC) generates truncated proteins that may have toxic dominant negative effects. Nonsense-mediated decay (NMD) is an mRNA surveillance pathway that degrades PTC-containing mRNAs to limit the production of truncated proteins. NMD activation requires a ribosome terminating translation at a PTC, but what happens to the polypeptides synthesized during the translation cycle needed to activate NMD is incompletely understood. Here, by establishing reporter systems that encode the same polypeptide sequence before a normal termination codon or PTC, we show that termination of protein synthesis at a PTC is sufficient to selectively destabilize polypeptides in mammalian cells. Proteasome inhibition specifically rescues the levels of nascent polypeptides produced from PTC-containing mRNAs within an hour, but also disrupts mRNA homeostasis within a few hours. PTC-terminated polypeptide destabilization is also alleviated by depleting the central NMD factor UPF1 or SMG1, the kinase that phosphorylates UPF1 to activate NMD, but not by inhibiting SMG1 kinase activity. Our results suggest that polypeptide degradation is linked to PTC recognition in mammalian cells and clarify a framework to investigate these mechanisms.

Mouse Genome Informatics (MGI): latest news from MGD and GXD

The Mouse Genome Informatics (MGI) database system combines multiple expertly curated community data resources into a shared knowledge management ecosystem united by common metadata annotation standards. MGI's mission is to facilitate the use of the mouse as an experimental model for understanding the genetic and genomic basis of human health and disease. MGI is the authoritative source for mouse gene, allele, and strain nomenclature and is the primary source of mouse phenotype annotations, functional annotations, developmental gene expression information, and annotations of mouse models with human diseases. MGI maintains mouse anatomy and phenotype ontologies and contributes to the development of the Gene Ontology and Disease Ontology and uses these ontologies as standard terminologies for annotation. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are MGI's two major knowledgebases. Here, we highlight some of the recent changes and enhancements to MGD and GXD that have been implemented in response to changing needs of the biomedical research community and to improve the efficiency of expert curation. MGI can be accessed freely at http://www.informatics.jax.org .

Harnessing model organism genomics to underpin the machine learning-based prediction of essential genes in eukaryotes - Biotechnological implications

The availability of high-quality genomes and advances in functional genomics have enabled large-scale studies of essential genes in model eukaryotes, including the 'elegant worm' (Caenorhabditis elegans; Nematoda) and the 'vinegar fly' (Drosophila melanogaster; Arthropoda). However, this is not the case for other, much less-studied organisms, such as socioeconomically important parasites, for which functional genomic platforms usually do not exist. Thus, there is a need to develop innovative techniques or approaches for the prediction, identification and investigation of essential genes. A key approach that could enable the prediction of such genes is machine learning (ML). Here, we undertake an historical review of experimental and computational approaches employed for the characterisation of essential genes in eukaryotes, with a particular focus on model ecdysozoans (C. elegans and D. melanogaster), and discuss the possible applicability of ML-approaches to organisms such as socioeconomically important parasites. We highlight some recent results showing that high-performance ML, combined with feature engineering, allows a reliable prediction of essential genes from extensive, publicly available 'omic data sets, with major potential to prioritise such genes (with statistical confidence) for subsequent functional genomic validation. These findings could 'open the door' to fundamental and applied research areas. Evidence of some commonality in the essential gene-complement between these two organisms indicates that an ML-engineering approach could find broader applicability to ecdysozoans such as parasitic nematodes or arthropods, provided that suitably large and informative data sets become/are available for proper feature engineering, and for the robust training and validation of algorithms. This area warrants detailed exploration to, for example, facilitate the identification and characterisation of essential molecules as novel targets for drugs and vaccines against parasitic diseases. This focus is particularly important, given the substantial impact that such diseases have worldwide, and the current challenges associated with their prevention and control and with drug resistance in parasite populations.

Protein Integrated Network Analysis to Reveal Potential Drug Targets Against Extended Drug-Resistant Mycobacterium tuberculosis XDR1219

The reconstruction and analysis of the protein-protein interaction (PPI) network is a powerful approach to understand the complex biological and molecular functions in normal and disease states of the cell. The interactome of most organisms is largely unidentified except some model organisms. The current study focused on the construction of PPI network for the human pathogen Mycobacterium tuberculosis (MTB)-resistant strain XDR1219 using computational methods. In this work, a bioinformatics approach was employed to reveal potential drug targets. The pipeline adopted the combination of an extensive integrated network analysis that led to identify 22 key proteins involved in drug resistance, resistant metabolic pathways, virulence, pathogenesis and persistency of the infection. The MTB XDR1219 interactome consists of 11,383 non-redundant PPIs among 1499 proteins covering 38% of the entire MTB XDR1219 proteome. The overall quality of the network was assessed and topological parameters of the PPI were calculated. The predicted interactions were functionally annotated and their relevance was assessed with the functional similarity. The study attempts to present the interactome of previously unidentified MTB XDR1219 and revealed potential drug targets that can be further explored by scientific community.

Thousands of induced germline mutations affecting immune cells identified by automated meiotic mapping coupled with machine learning

Forward genetic studies use meiotic mapping to adduce evidence that a particular mutation, normally induced by a germline mutagen, is causative of a particular phenotype. Particularly in small pedigrees, cosegregation of multiple mutations, occasional unawareness of mutations, and paucity of homozygotes may lead to erroneous declarations of cause and effect. We sought to improve the identification of mutations causing immune phenotypes in mice by creating Candidate Explorer (CE), a machine-learning software program that integrates 67 features of genetic mapping data into a single numeric score, mathematically convertible to the probability of verification of any putative mutation-phenotype association. At this time, CE has evaluated putative mutation-phenotype associations arising from screening damaging mutations in ∼55% of mouse genes for effects on flow cytometry measurements of immune cells in the blood. CE has therefore identified more than half of genes within which mutations can be causative of flow cytometric phenovariation in Mus musculus The majority of these genes were not previously known to support immune function or homeostasis. Mouse geneticists will find CE data informative in identifying causative mutations within quantitative trait loci, while clinical geneticists may use CE to help connect causative variants with rare heritable diseases of immunity, even in the absence of linkage information. CE displays integrated mutation, phenotype, and linkage data, and is freely available for query online.

Semi-Supervised Learning Using Hierarchical Mixture Models: Gene Essentiality Case Study

Integrating gene-level data is useful for predicting the role of genes in biological processes. This problem has typically focused on supervised classification, which requires large training sets of positive and negative examples. However, training data sets that are too small for supervised approaches can still provide valuable information. We describe a hierarchical mixture model that uses limited positively labeled gene training data for semi-supervised learning. We focus on the problem of predicting essential genes, where a gene is required for the survival of an organism under particular conditions. We applied cross-validation and found that the inclusion of positively labeled samples in a semi-supervised learning framework with the hierarchical mixture model improves the detection of essential genes compared to unsupervised, supervised, and other semi-supervised approaches. There was also improved prediction performance when genes are incorrectly assumed to be non-essential. Our comparisons indicate that the incorporation of even small amounts of existing knowledge improves the accuracy of prediction and decreases variability in predictions. Although we focused on gene essentiality, the hierarchical mixture model and semi-supervised framework is standard for problems focused on prediction of genes or other features, with multiple data types characterizing the feature, and a small set of positive labels.

Non-essential ribosomal proteins in bacteria and archaea identified using COGs

Ribosomal proteins (RPs) are highly conserved across the bacterial and archaeal domains. Although many RPs are essential for survival, genome analysis demonstrates the absence of some RP genes in many bacterial and archaeal genomes. Furthermore, global transposon mutagenesis and/or targeted deletion showed that elimination of some RP genes had only a moderate effect on the bacterial growth rate. Here, we systematically analyze the evolutionary conservation of RPs in prokaryotes by compiling the list of the ribosomal genes that are missing from one or more genomes in the recently updated version of the Clusters of Orthologous Genes (COG) database. Some of these absences occurred because the respective genes carried frameshifts, presumably, resulting from sequencing errors, while others were overlooked and not translated during genome annotation. Apart from these annotation errors, we identified multiple genuine losses of RP genes in a variety of bacteria and archaea. Some of these losses are clade-specific, whereas others occur in symbionts and parasites with dramatically reduced genomes. The lists of computationally and experimentally defined non-essential ribosomal genes show a substantial overlap, revealing a common trend in prokaryote ribosome evolution that could be linked to the architecture and assembly of the ribosomes. Thus, RPs that are located at the surface of the ribosome and/or are incorporated at a late stage of ribosome assembly are more likely to be non-essential and to be lost during microbial evolution, particularly, in the course of genome compaction.IMPORTANCEIn many prokaryote genomes, one or more ribosomal protein (RP) genes are missing. Analysis of 1,309 prokaryote genomes included in the COG database shows that only about half of the RPs are universally conserved in bacteria and archaea. In contrast, up to 16 other RPs are missing in some genomes, primarily, tiny (<1 Mb) genomes of host-associated bacteria and archaea. Ten universal and nine archaea-specific ribosomal proteins show clear patterns of lineage-specific gene loss. Most of the RPs that are frequently lost from bacterial genomes are located on the ribosome periphery and are non-essential in Escherichia coli and Bacillus subtilis These results reveal general trends and common constraints in the architecture and evolution of ribosomes in prokaryotes.