Using artificial intelligence to document the hidden RNA virosphere.

Current metagenomic tools can fail to identify highly divergent RNA viruses. We developed a deep learning algorithm, termed LucaProt, to discover highly divergent RNA-dependent RNA polymerase (RdRP) sequences in 10,487 metatranscriptomes generated from diverse global ecosystems. LucaProt integrates both sequence and predicted structural information, enabling the accurate detection of RdRP sequences. Using this approach, we identified 161,979 potential RNA virus species and 180 RNA virus supergroups, including many previously poorly studied groups, as well as RNA virus genomes of exceptional length (up to 47,250 nucleotides) and genomic complexity. A subset of these novel RNA viruses was confirmed by RT-PCR and RNA/DNA sequencing. Newly discovered RNA viruses were present in diverse environments, including air, hot springs, and hydrothermal vents, with virus diversity and abundance varying substantially among ecosystems. This study advances virus discovery, highlights the scale of the virosphere, and provides computational tools to better document the global RNA virome.

Disease-associated astrocyte epigenetic memory promotes CNS pathology.

Disease-associated astrocyte subsets contribute to the pathology of neurologic diseases, including multiple sclerosis and experimental autoimmune encephalomyelitis1-8 (EAE), an experimental model for multiple sclerosis. However, little is known about the stability of these astrocyte subsets and their ability to integrate past stimulation events. Here we report the identification of an epigenetically controlled memory astrocyte subset that exhibits exacerbated pro-inflammatory responses upon rechallenge. Specifically, using a combination of single-cell RNA sequencing, assay for transposase-accessible chromatin with sequencing, chromatin immunoprecipitation with sequencing, focused interrogation of cells by nucleic acid detection and sequencing, and cell-specific in vivo CRISPR-Cas9-based genetic perturbation studies we established that astrocyte memory is controlled by the metabolic enzyme ATP-citrate lyase (ACLY), which produces acetyl coenzyme A (acetyl-CoA) that is used by histone acetyltransferase p300 to control chromatin accessibility. The number of ACLY+p300+ memory astrocytes is increased in acute and chronic EAE models, and their genetic inactivation ameliorated EAE. We also detected the pro-inflammatory memory phenotype in human astrocytes in vitro; single-cell RNA sequencing and immunohistochemistry studies detected increased numbers of ACLY+p300+ astrocytes in chronic multiple sclerosis lesions. In summary, these studies define an epigenetically controlled memory astrocyte subset that promotes CNS pathology in EAE and, potentially, multiple sclerosis. These findings may guide novel therapeutic approaches for multiple sclerosis and other neurologic diseases.

Compendium of human transcription factor effector domains.

Transcription factors (TFs) regulate gene expression by binding to DNA sequences and modulating transcriptional activity through their effector domains. Despite the central role of effector domains in TF function, there is a current lack of a comprehensive resource and characterization of effector domains. Here, we provide a catalog of 924 effector domains across 594 human TFs. Using this catalog, we characterized the amino acid composition of effector domains, their conservation across species and across the human population, and their roles in human diseases. Furthermore, we provide a classification system for effector domains that constitutes a valuable resource and a blueprint for future experimental studies of TF effector domain function.

Identification of astrocyte regulators by nucleic acid cytometry.

Multiple sclerosis is a chronic inflammatory disease of the central nervous system1. Astrocytes are heterogeneous glial cells that are resident in the central nervous system and participate in the pathogenesis of multiple sclerosis and its model experimental autoimmune encephalomyelitis2,3. However, few unique surface markers are available for the isolation of astrocyte subsets, preventing their analysis and the identification of candidate therapeutic targets; these limitations are further amplified by the rarity of pathogenic astrocytes. Here, to address these challenges, we developed focused interrogation of cells by nucleic acid detection and sequencing (FIND-seq), a high-throughput microfluidic cytometry method that combines encapsulation of cells in droplets, PCR-based detection of target nucleic acids and droplet sorting to enable in-depth transcriptomic analyses of cells of interest at single-cell resolution. We applied FIND-seq to study the regulation of astrocytes characterized by the splicing-driven activation of the transcription factor XBP1, which promotes disease pathology in multiple sclerosis and experimental autoimmune encephalomyelitis4. Using FIND-seq in combination with conditional-knockout mice, in vivo CRISPR-Cas9-driven genetic perturbation studies and bulk and single-cell RNA sequencing analyses of samples from mouse experimental autoimmune encephalomyelitis and humans with multiple sclerosis, we identified a new role for the nuclear receptor NR3C2 and its corepressor NCOR2 in limiting XBP1-driven pathogenic astrocyte responses. In summary, we used FIND-seq to identify a therapeutically targetable mechanism that limits XBP1-driven pathogenic astrocyte responses. FIND-seq enables the investigation of previously inaccessible cells, including rare cell subsets defined by unique gene expression signatures or other nucleic acid markers.

B-cell-specific checkpoint molecules that regulate anti-tumour immunity.

The role of B cells in anti-tumour immunity is still debated and, accordingly, immunotherapies have focused on targeting T and natural killer cells to inhibit tumour growth1,2. Here, using high-throughput flow cytometry as well as bulk and single-cell RNA-sequencing and B-cell-receptor-sequencing analysis of B cells temporally during B16F10 melanoma growth, we identified a subset of B cells that expands specifically in the draining lymph node over time in tumour-bearing mice. The expanding B cell subset expresses the cell surface molecule T cell immunoglobulin and mucin domain 1 (TIM-1, encoded by Havcr1) and a unique transcriptional signature, including multiple co-inhibitory molecules such as PD-1, TIM-3, TIGIT and LAG-3. Although conditional deletion of these co-inhibitory molecules on B cells had little or no effect on tumour burden, selective deletion of Havcr1 in B cells both substantially inhibited tumour growth and enhanced effector T cell responses. Loss of TIM-1 enhanced the type 1 interferon response in B cells, which augmented B cell activation and increased antigen presentation and co-stimulation, resulting in increased expansion of tumour-specific effector T cells. Our results demonstrate that manipulation of TIM-1-expressing B cells enables engagement of the second arm of adaptive immunity to promote anti-tumour immunity and inhibit tumour growth.

Lactate limits CNS autoimmunity by stabilizing HIF-1α in dendritic cells.

Dendritic cells (DCs) have a role in the development and activation of self-reactive pathogenic T cells1,2. Genetic variants that are associated with the function of DCs have been linked to autoimmune disorders3,4, and DCs are therefore attractive therapeutic targets for such diseases. However, developing DC-targeted therapies for autoimmunity requires identification of the mechanisms that regulate DC function. Here, using single-cell and bulk transcriptional and metabolic analyses in combination with cell-specific gene perturbation studies, we identify a regulatory loop of negative feedback that operates in DCs to limit immunopathology. Specifically, we find that lactate, produced by activated DCs and other immune cells, boosts the expression of NDUFA4L2 through a mechanism mediated by hypoxia-inducible factor 1α (HIF-1α). NDUFA4L2 limits the production of mitochondrial reactive oxygen species that activate XBP1-driven transcriptional modules in DCs that are involved in the control of pathogenic autoimmune T cells. We also engineer a probiotic that produces lactate and suppresses T cell autoimmunity through the activation of HIF-1α-NDUFA4L2 signalling in DCs. In summary, we identify an immunometabolic pathway that regulates DC function, and develop a synthetic probiotic for its therapeutic activation.

Identification of environmental factors that promote intestinal inflammation.

Genome-wide association studies have identified risk loci linked to inflammatory bowel disease (IBD)1-a complex chronic inflammatory disorder of the gastrointestinal tract. The increasing prevalence of IBD in industrialized countries and the augmented disease risk observed in migrants who move into areas of higher disease prevalence suggest that environmental factors are also important determinants of IBD susceptibility and severity2. However, the identification of environmental factors relevant to IBD and the mechanisms by which they influence disease has been hampered by the lack of platforms for their systematic investigation. Here we describe an integrated systems approach, combining publicly available databases, zebrafish chemical screens, machine learning and mouse preclinical models to identify environmental factors that control intestinal inflammation. This approach established that the herbicide propyzamide increases inflammation in the small and large intestine. Moreover, we show that an AHR-NF-κB-C/EBPß signalling axis operates in T cells and dendritic cells to promote intestinal inflammation, and is targeted by propyzamide. In conclusion, we developed a pipeline for the identification of environmental factors and mechanisms of pathogenesis in IBD and, potentially, other inflammatory diseases.

Gut-licensed IFNγ+ NK cells drive LAMP1+TRAIL+ anti-inflammatory astrocytes.

Astrocytes are glial cells that are abundant in the central nervous system (CNS) and that have important homeostatic and disease-promoting functions1. However, little is known about the homeostatic anti-inflammatory activities of astrocytes and their regulation. Here, using high-throughput flow cytometry screening, single-cell RNA sequencing and CRISPR-Cas9-based cell-specific in vivo genetic perturbations in mice, we identify a subset of astrocytes that expresses the lysosomal protein LAMP12 and the death receptor ligand TRAIL3. LAMP1+TRAIL+ astrocytes limit inflammation in the CNS by inducing T cell apoptosis through TRAIL-DR5 signalling. In homeostatic conditions, the expression of TRAIL in astrocytes is driven by interferon-γ (IFNγ) produced by meningeal natural killer (NK) cells, in which IFNγ expression is modulated by the gut microbiome. TRAIL expression in astrocytes is repressed by molecules produced by T cells and microglia in the context of inflammation. Altogether, we show that LAMP1+TRAIL+ astrocytes limit CNS inflammation by inducing T cell apoptosis, and that this astrocyte subset is maintained by meningeal IFNγ+ NK cells that are licensed by the microbiome.

MoDAFold: a strategy for predicting the structure of missense mutant protein based on AlphaFold2 and molecular dynamics.

Protein structure prediction is a longstanding issue crucial for identifying new drug targets and providing a mechanistic understanding of protein functions. To enhance the progress in this field, a spectrum of computational methodologies has been cultivated. AlphaFold2 has exhibited exceptional precision in predicting wild-type protein structures, with performance exceeding that of other methods. However, predicting the structures of missense mutant proteins using AlphaFold2 remains challenging due to the intricate and substantial structural alterations caused by minor sequence variations in the mutant proteins. Molecular dynamics (MD) has been validated for precisely capturing changes in amino acid interactions attributed to protein mutations. Therefore, for the first time, a strategy entitled 'MoDAFold' was proposed to improve the accuracy and reliability of missense mutant protein structure prediction by combining AlphaFold2 with MD. Multiple case studies have confirmed the superior performance of MoDAFold compared to other methods, particularly AlphaFold2.

MAFG-driven astrocytes promote CNS inflammation.

Multiple sclerosis is a chronic inflammatory disease of the CNS1. Astrocytes contribute to the pathogenesis of multiple sclerosis2, but little is known about the heterogeneity of astrocytes and its regulation. Here we report the analysis of astrocytes in multiple sclerosis and its preclinical model experimental autoimmune encephalomyelitis (EAE) by single-cell RNA sequencing in combination with cell-specific Ribotag RNA profiling, assay for transposase-accessible chromatin with sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing (ChIP-seq), genome-wide analysis of DNA methylation and in vivo CRISPR-Cas9-based genetic perturbations. We identified astrocytes in EAE and multiple sclerosis that were characterized by decreased expression of NRF2 and increased expression of MAFG, which cooperates with MAT2α to promote DNA methylation and represses antioxidant and anti-inflammatory transcriptional programs. Granulocyte-macrophage colony-stimulating factor (GM-CSF) signalling in astrocytes drives the expression of MAFG and MAT2α and pro-inflammatory transcriptional modules, contributing to CNS pathology in EAE and, potentially, multiple sclerosis. Our results identify candidate therapeutic targets in multiple sclerosis.

DrugMAP: molecular atlas and pharma-information of all drugs.

The efficacy and safety of drugs are widely known to be determined by their interactions with multiple molecules of pharmacological importance, and it is therefore essential to systematically depict the molecular atlas and pharma-information of studied drugs. However, our understanding of such information is neither comprehensive nor precise, which necessitates the construction of a new database providing a network containing a large number of drugs and their interacting molecules. Here, a new database describing the molecular atlas and pharma-information of drugs (DrugMAP) was therefore constructed. It provides a comprehensive list of interacting molecules for >30 000 drugs/drug candidates, gives the differential expression patterns for >5000 interacting molecules among different disease sites, ADME (absorption, distribution, metabolism and excretion)-relevant organs and physiological tissues, and weaves a comprehensive and precise network containing >200 000 interactions among drugs and molecules. With the great efforts made to clarify the complex mechanism underlying drug pharmacokinetics and pharmacodynamics and rapidly emerging interests in artificial intelligence (AI)-based network analyses, DrugMAP is expected to become an indispensable supplement to existing databases to facilitate drug discovery. It is now fully and freely accessible at: https://idrblab.org/drugmap/.

A task-specific encoding algorithm for RNAs and RNA-associated interactions based on convolutional autoencoder.

RNAs play essential roles in diverse physiological and pathological processes by interacting with other molecules (RNA/protein/compound), and various computational methods are available for identifying these interactions. However, the encoding features provided by existing methods are limited and the existing tools does not offer an effective way to integrate the interacting partners. In this study, a task-specific encoding algorithm for RNAs and RNA-associated interactions was therefore developed. This new algorithm was unique in (a) realizing comprehensive RNA feature encoding by introducing a great many of novel features and (b) enabling task-specific integration of interacting partners using convolutional autoencoder-directed feature embedding. Compared with existing methods/tools, this novel algorithm demonstrated superior performances in diverse benchmark testing studies. This algorithm together with its source code could be readily accessed by all user at: https://idrblab.org/corain/ and https://github.com/idrblab/corain/.

ncRNAInter: a novel strategy based on graph neural network to discover interactions between lncRNA and miRNA.

In recent years, many studies have illustrated the significant role that non-coding RNA (ncRNA) plays in biological activities, in which lncRNA, miRNA and especially their interactions have been proved to affect many biological processes. Some in silico methods have been proposed and applied to identify novel lncRNA-miRNA interactions (LMIs), but there are still imperfections in their RNA representation and information extraction approaches, which imply there is still room for further improving their performances. Meanwhile, only a few of them are accessible at present, which limits their practical applications. The construction of a new tool for LMI prediction is thus imperative for the better understanding of their relevant biological mechanisms. This study proposed a novel method, ncRNAInter, for LMI prediction. A comprehensive strategy for RNA representation and an optimized deep learning algorithm of graph neural network were utilized in this study. ncRNAInter was robust and showed better performance of 26.7% higher Matthews correlation coefficient than existing reputable methods for human LMI prediction. In addition, ncRNAInter proved its universal applicability in dealing with LMIs from various species and successfully identified novel LMIs associated with various diseases, which further verified its effectiveness and usability. All source code and datasets are freely available at https://github.com/idrblab/ncRNAInter.

ConSIG: consistent discovery of molecular signature from OMIC data.

The discovery of proper molecular signature from OMIC data is indispensable for determining biological state, physiological condition, disease etiology, and therapeutic response. However, the identified signature is reported to be highly inconsistent, and there is little overlap among the signatures identified from different biological datasets. Such inconsistency raises doubts about the reliability of reported signatures and significantly hampers its biological and clinical applications. Herein, an online tool, ConSIG, was constructed to realize consistent discovery of gene/protein signature from any uploaded transcriptomic/proteomic data. This tool is unique in a) integrating a novel strategy capable of significantly enhancing the consistency of signature discovery, b) determining the optimal signature by collective assessment, and c) confirming the biological relevance by enriching the disease/gene ontology. With the increasingly accumulated concerns about signature consistency and biological relevance, this online tool is expected to be used as an essential complement to other existing tools for OMIC-based signature discovery. ConSIG is freely accessible to all users without login requirement at https://idrblab.org/consig/.

Biological activities of drug inactive ingredients.

In a drug formulation (DFM), the major components by mass are not Active Pharmaceutical Ingredient (API) but rather Drug Inactive Ingredients (DIGs). DIGs can reach much higher concentrations than that achieved by API, which raises great concerns about their clinical toxicities. Therefore, the biological activities of DIG on physiologically relevant target are widely demanded by both clinical investigation and pharmaceutical industry. However, such activity data are not available in any existing pharmaceutical knowledge base, and their potentials in predicting the DIG-target interaction have not been evaluated yet. In this study, the comprehensive assessment and analysis on the biological activities of DIGs were therefore conducted. First, the largest number of DIGs and DFMs were systematically curated and confirmed based on all drugs approved by US Food and Drug Administration. Second, comprehensive activities for both DIGs and DFMs were provided for the first time to pharmaceutical community. Third, the biological targets of each DIG and formulation were fully referenced to available databases that described their pharmaceutical/biological characteristics. Finally, a variety of popular artificial intelligence techniques were used to assess the predictive potential of DIGs' activity data, which was the first evaluation on the possibility to predict DIG's activity. As the activities of DIGs are critical for current pharmaceutical studies, this work is expected to have significant implications for the future practice of drug discovery and precision medicine.

RNAenrich: a web server for non-coding RNA enrichment.

MOTIVATION: With the rapid advances of RNA sequencing and microarray technologies in non-coding RNA (ncRNA) research, functional tools that perform enrichment analysis for ncRNAs are needed. On the one hand, because of the rapidly growing interest in circRNAs, snoRNAs, and piRNAs, it is essential to develop tools for enrichment analysis for these newly emerged ncRNAs. On the other hand, due to the key role of ncRNAs' interacting target in the determination of their function, the interactions between ncRNA and its corresponding target should be fully considered in functional enrichment. Based on the ncRNA-mRNA/protein-function strategy, some tools have been developed to functionally analyze a single type of ncRNA (the majority focuses on miRNA); in addition, some tools adopt predicted target data and lead to only low-confidence results. RESULTS: Herein, an online tool named RNAenrich was developed to enable the comprehensive and accurate enrichment analysis of ncRNAs. It is unique in (i) realizing the enrichment analysis for various RNA types in humans and mice, such as miRNA, lncRNA, circRNA, snoRNA, piRNA, and mRNA; (ii) extending the analysis by introducing millions of experimentally validated data of RNA-target interactions as a built-in database; and (iii) providing a comprehensive interacting network among various ncRNAs and targets to facilitate the mechanistic study of ncRNA function. Importantly, RNAenrich led to a more comprehensive and accurate enrichment analysis in a COVID-19-related miRNA case, which was largely attributed to its coverage of comprehensive ncRNA-target interactions. AVAILABILITY AND IMPLEMENTATION: RNAenrich is now freely accessible at https://idrblab.org/rnaenr/.

Naringin attenuates Actinobacillus pleuropneumoniae-induced acute lung injury via MAPK/NF-κB and Keap1/Nrf2/HO-1 pathway.

Actinobacillus pleuropneumoniae (APP) causes porcine pleuropneumonia (PCP), which is clinically characterized by acute hemorrhagic, necrotizing pneumonia, and chronic fibrinous pneumonia. Although many measures have been taken to prevent the disease, prevention and control of the disease are becoming increasingly difficult due to the abundance of APP sera, weak vaccine cross-protection, and increasing antibiotic resistance in APP. Therefore, there is an urgent need to develop novel drugs against APP infection to prevent the spread of APP. Naringin (NAR) has been reported to have an excellent therapeutic effect on pulmonary diseases, but its therapeutic effect on lung injury caused by APP is not apparent. Our research has shown that NAR was able to alleviate APP-induced weight loss and quantity of food taken and reduce the number of WBCs and NEs in peripheral blood in mice; pathological tissue sections showed that NAR was able to prevent and control APP-induced pathological lung injury effectively; based on the establishment of an in vivo/in vitro model of APP inflammation, it was found that NAR was able to play an anti-inflammatory role through inhibiting the MAPK/NF-κB signaling pathway and exerting anti-inflammatory effects; additionally, NAR activating the Nrf2 signalling pathway, increasing the secretion of antioxidant enzymes Nqo1, CAT, and SOD1, inhibiting the secretion of oxidative damage factors NOS2 and COX2, and enhancing the antioxidant stress ability, thus playing an antioxidant role. In summary, NAR can relieve severe lung injury caused by APP by reducing excessive inflammatory response and improving antioxidant capacity.

SYNBIP: synthetic binding proteins for research, diagnosis and therapy.

The success of protein engineering and design has extensively expanded the protein space, which presents a promising strategy for creating next-generation proteins of diverse functions. Among these proteins, the synthetic binding proteins (SBPs) are smaller, more stable, less immunogenic, and better of tissue penetration than others, which make the SBP-related data attracting extensive interest from worldwide scientists. However, no database has been developed to systematically provide the valuable information of SBPs yet. In this study, a database named 'Synthetic Binding Proteins for Research, Diagnosis, and Therapy (SYNBIP)' was thus introduced. This database is unique in (a) comprehensively describing thousands of SBPs from the perspectives of scaffolds, biophysical & functional properties, etc.; (b) panoramically illustrating the binding targets & the broad application of each SBP and (c) enabling a similarity search against the sequences of all SBPs and their binding targets. Since SBP is a human-made protein that has not been found in nature, the discovery of novel SBPs relied heavily on experimental protein engineering and could be greatly facilitated by in-silico studies (such as AI and computational modeling). Thus, the data provided in SYNBIP could lay a solid foundation for the future development of novel SBPs. The SYNBIP is accessible without login requirement at both official (https://idrblab.org/synbip/) and mirror (http://synbip.idrblab.net/) sites.

NPCDR: natural product-based drug combination and its disease-specific molecular regulation.

Natural product (NP) has a long history in promoting modern drug discovery, which has derived or inspired a large number of currently prescribed drugs. Recently, the NPs have emerged as the ideal candidates to combine with other therapeutic strategies to deal with the persistent challenge of conventional therapy, and the molecular regulation mechanism underlying these combinations is crucial for the related communities. Thus, it is urgently demanded to comprehensively provide the disease-specific molecular regulation data for various NP-based drug combinations. However, no database has been developed yet to describe such valuable information. In this study, a newly developed database entitled 'Natural Product-based Drug Combination and Its Disease-specific Molecular Regulation (NPCDR)' was thus introduced. This database was unique in (a) providing the comprehensive information of NP-based drug combinations & describing their clinically or experimentally validated therapeutic effect, (b) giving the disease-specific molecular regulation data for a number of NP-based drug combinations, (c) fully referencing all NPs, drugs, regulated molecules/pathways by cross-linking them to the available databases describing their biological or pharmaceutical characteristics. Therefore, NPCDR is expected to have great implications for the future practice of network pharmacology, medical biochemistry, drug design, and medicinal chemistry. This database is now freely accessible without any login requirement at both official (https://idrblab.org/npcdr/) and mirror (http://npcdr.idrblab.net/) sites.