Discovery of antimicrobial peptides in the global microbiome with machine learning.

Novel antibiotics are urgently needed to combat the antibiotic-resistance crisis. We present a machine-learning-based approach to predict antimicrobial peptides (AMPs) within the global microbiome and leverage a vast dataset of 63,410 metagenomes and 87,920 prokaryotic genomes from environmental and host-associated habitats to create the AMPSphere, a comprehensive catalog comprising 863,498 non-redundant peptides, few of which match existing databases. AMPSphere provides insights into the evolutionary origins of peptides, including by duplication or gene truncation of longer sequences, and we observed that AMP production varies by habitat. To validate our predictions, we synthesized and tested 100 AMPs against clinically relevant drug-resistant pathogens and human gut commensals both in vitro and in vivo. A total of 79 peptides were active, with 63 targeting pathogens. These active AMPs exhibited antibacterial activity by disrupting bacterial membranes. In conclusion, our approach identified nearly one million prokaryotic AMP sequences, an open-access resource for antibiotic discovery.

A pangenome reference of 36 Chinese populations.

Human genomics is witnessing an ongoing paradigm shift from a single reference sequence to a pangenome form, but populations of Asian ancestry are underrepresented. Here we present data from the first phase of the Chinese Pangenome Consortium, including a collection of 116 high-quality and haplotype-phased de novo assemblies based on 58 core samples representing 36 minority Chinese ethnic groups. With an average 30.65× high-fidelity long-read sequence coverage, an average contiguity N50 of more than 35.63 megabases and an average total size of 3.01 gigabases, the CPC core assemblies add 189 million base pairs of euchromatic polymorphic sequences and 1,367 protein-coding gene duplications to GRCh38. We identified 15.9 million small variants and 78,072 structural variants, of which 5.9 million small variants and 34,223 structural variants were not reported in a recently released pangenome reference1. The Chinese Pangenome Consortium data demonstrate a remarkable increase in the discovery of novel and missing sequences when individuals are included from underrepresented minority ethnic groups. The missing reference sequences were enriched with archaic-derived alleles and genes that confer essential functions related to keratinization, response to ultraviolet radiation, DNA repair, immunological responses and lifespan, implying great potential for shedding new light on human evolution and recovering missing heritability in complex disease mapping.

Towards the biogeography of prokaryotic genes.

Microbial genes encode the majority of the functional repertoire of life on earth. However, despite increasing efforts in metagenomic sequencing of various habitats1-3, little is known about the distribution of genes across the global biosphere, with implications for human and planetary health. Here we constructed a non-redundant gene catalogue of 303 million species-level genes (clustered at 95% nucleotide identity) from 13,174 publicly available metagenomes across 14 major habitats and use it to show that most genes are specific to a single habitat. The small fraction of genes found in multiple habitats is enriched in antibiotic-resistance genes and markers for mobile genetic elements. By further clustering these species-level genes into 32 million protein families, we observed that a small fraction of these families contain the majority of the genes (0.6% of families account for 50% of the genes). The majority of species-level genes and protein families are rare. Furthermore, species-level genes, and in particular the rare ones, show low rates of positive (adaptive) selection, supporting a model in which most genetic variability observed within each protein family is neutral or nearly neutral.

STAB2: an updated spatio-temporal cell atlas of the human and mouse brain.

The brain is constituted of heterogeneous types of neuronal and non-neuronal cells, which are organized into distinct anatomical regions, and show precise regulation of gene expression during development, aging and function. In the current database release, STAB2 provides a systematic cellular map of the human and mouse brain by integrating recently published large-scale single-cell and single-nucleus RNA-sequencing datasets from diverse regions and across lifespan. We applied a hierarchical strategy of unsupervised clustering on the integrated single-cell transcriptomic datasets to precisely annotate the cell types and subtypes in the human and mouse brain. Currently, STAB2 includes 71 and 61 different cell subtypes defined in the human and mouse brain, respectively. It covers 63 subregions and 15 developmental stages of human brain, and 38 subregions and 30 developmental stages of mouse brain, generating a comprehensive atlas for exploring spatiotemporal transcriptomic dynamics in the mammalian brain. We also augmented web interfaces for querying and visualizing the gene expression in specific cell types. STAB2 is freely available at https://mai.fudan.edu.cn/stab2.

An augmented Mendelian randomization approach provides causality of brain imaging features on complex traits in a single biobank-scale dataset.

Mendelian randomization (MR) is an effective approach for revealing causal risk factors that underpin complex traits and diseases. While MR has been more widely applied under two-sample settings, it is more promising to be used in one single large cohort given the rise of biobank-scale datasets that simultaneously contain genotype data, brain imaging data, and matched complex traits from the same individual. However, most existing multivariable MR methods have been developed for two-sample setting or a small number of exposures. In this study, we introduce a one-sample multivariable MR method based on partial least squares and Lasso regression (MR-PL). MR-PL is capable of considering the correlation among exposures (e.g., brain imaging features) when the number of exposures is extremely upscaled, while also correcting for winner's curse bias. We performed extensive and systematic simulations, and demonstrated the robustness and reliability of our method. Comprehensive simulations confirmed that MR-PL can generate more precise causal estimates with lower false positive rates than alternative approaches. Finally, we applied MR-PL to the datasets from UK Biobank to reveal the causal effects of 36 white matter tracts on 180 complex traits, and showed putative white matter tracts that are implicated in smoking, blood vascular function-related traits, and eating behaviors.

Deciphering the genetic architecture of human brain structure and function: a brief survey on recent advances of neuroimaging genomics.

Brain imaging genomics is an emerging interdisciplinary field, where integrated analysis of multimodal medical image-derived phenotypes (IDPs) and multi-omics data, bridging the gap between macroscopic brain phenotypes and their cellular and molecular characteristics. This approach aims to better interpret the genetic architecture and molecular mechanisms associated with brain structure, function and clinical outcomes. More recently, the availability of large-scale imaging and multi-omics datasets from the human brain has afforded the opportunity to the discovering of common genetic variants contributing to the structural and functional IDPs of the human brain. By integrative analyses with functional multi-omics data from the human brain, a set of critical genes, functional genomic regions and neuronal cell types have been identified as significantly associated with brain IDPs. Here, we review the recent advances in the methods and applications of multi-omics integration in brain imaging analysis. We highlight the importance of functional genomic datasets in understanding the biological functions of the identified genes and cell types that are associated with brain IDPs. Moreover, we summarize well-known neuroimaging genetics datasets and discuss challenges and future directions in this field.

A survey on computational strategies for genome-resolved gut metagenomics.

Recovering high-quality metagenome-assembled genomes (HQ-MAGs) is critical for exploring microbial compositions and microbe-phenotype associations. However, multiple sequencing platforms and computational tools for this purpose may confuse researchers and thus call for extensive evaluation. Here, we systematically evaluated a total of 40 combinations of popular computational tools and sequencing platforms (i.e. strategies), involving eight assemblers, eight metagenomic binners and four sequencing technologies, including short-, long-read and metaHiC sequencing. We identified the best tools for the individual tasks (e.g. the assembly and binning) and combinations (e.g. generating more HQ-MAGs) depending on the availability of the sequencing data. We found that the combination of the hybrid assemblies and metaHiC-based binning performed best, followed by the hybrid and long-read assemblies. More importantly, both long-read and metaHiC sequencings link more mobile elements and antibiotic resistance genes to bacterial hosts and improve the quality of public human gut reference genomes with 32% (34/105) HQ-MAGs that were either of better quality than those in the Unified Human Gastrointestinal Genome catalog version 2 or novel.

Engineering of cytosine base editors with DNA damage minimization and editing scope diversification.

Cytosine base editors (CBEs), which enable precise C-to-T substitutions, have been restricted by potential safety risks, including DNA off-target edits, RNA off-target edits and additional genotoxicity such as DNA damages induced by double-strand breaks (DSBs). Though DNA and RNA off-target edits have been ameliorated via various strategies, evaluation and minimization of DSB-associated DNA damage risks for most CBEs remain to be resolved. Here we demonstrate that YE1, an engineered CBE variant with minimized DNA and RNA off-target edits, could induce prominent DSB-associated DNA damage risks, manifested as γH2AX accumulation in human cells. We then perform deaminase engineering for two deaminases lamprey LjCDA1 and human APOBEC3A, and generate divergent CBE variants with eliminated DSB-associated DNA damage risks, in addition to minimized DNA/RNA off-target edits. Furthermore, the editing scopes and sequence preferences of APOBEC3A-derived CBEs could be further diversified by internal fusion strategy. Taken together, this study provides updated evaluation platform for DSB-associated DNA damage risks of CBEs and further generates a series of safer toolkits with diversified editing signatures to expand their applications.

MorbidGCN: prediction of multimorbidity with a graph convolutional network based on integration of population phenotypes and disease network.

Exploring multimorbidity relationships among diseases is of great importance for understanding their shared mechanisms, precise diagnosis and treatment. However, the landscape of multimorbidities is still far from complete due to the complex nature of multimorbidity. Although various types of biological data, such as biomolecules and clinical symptoms, have been used to identify multimorbidities, the population phenotype information (e.g. physical activity and diet) remains less explored for multimorbidity. Here, we present a graph convolutional network (GCN) model, named MorbidGCN, for multimorbidity prediction by integrating population phenotypes and disease network. Specifically, MorbidGCN treats the multimorbidity prediction as a missing link prediction problem in the disease network, where a novel feature selection method is embedded to select important phenotypes. Benchmarking results on two large-scale multimorbidity data sets, i.e. the UK Biobank (UKB) and Human Disease Network (HuDiNe) data sets, demonstrate that MorbidGCN outperforms other competitive methods. With MorbidGCN, 9742 and 14 010 novel multimorbidities are identified in the UKB and HuDiNe data sets, respectively. Moreover, we notice that the selected phenotypes that are generally differentially distributed between multimorbidity patients and single-disease patients can help interpret multimorbidities and show potential for prognosis of multimorbidities.

An integrated brain-specific network identifies genes associated with neuropathologic and clinical traits of Alzheimer's disease.

Alzheimer's disease (AD) has a strong genetic predisposition. However, its risk genes remain incompletely identified. We developed an Alzheimer's brain gene network-based approach to predict AD-associated genes by leveraging the functional pattern of known AD-associated genes. Our constructed network outperformed existing networks in predicting AD genes. We then systematically validated the predictions using independent genetic, transcriptomic, proteomic data, neuropathological and clinical data. First, top-ranked genes were enriched in AD-associated pathways. Second, using external gene expression data from the Mount Sinai Brain Bank study, we found that the top-ranked genes were significantly associated with neuropathological and clinical traits, including the Consortium to Establish a Registry for Alzheimer's Disease score, Braak stage score and clinical dementia rating. The analysis of Alzheimer's brain single-cell RNA-seq data revealed cell-type-specific association of predicted genes with early pathology of AD. Third, by interrogating proteomic data in the Religious Orders Study and Memory and Aging Project and Baltimore Longitudinal Study of Aging studies, we observed a significant association of protein expression level with cognitive function and AD clinical severity. The network, method and predictions could become a valuable resource to advance the identification of risk genes for AD.

SemiBin2: self-supervised contrastive learning leads to better MAGs for short- and long-read sequencing.

MOTIVATION: Metagenomic binning methods to reconstruct metagenome-assembled genomes (MAGs) from environmental samples have been widely used in large-scale metagenomic studies. The recently proposed semi-supervised binning method, SemiBin, achieved state-of-the-art binning results in several environments. However, this required annotating contigs, a computationally costly and potentially biased process. RESULTS: We propose SemiBin2, which uses self-supervised learning to learn feature embeddings from the contigs. In simulated and real datasets, we show that self-supervised learning achieves better results than the semi-supervised learning used in SemiBin1 and that SemiBin2 outperforms other state-of-the-art binners. Compared to SemiBin1, SemiBin2 can reconstruct 8.3-21.5% more high-quality bins and requires only 25% of the running time and 11% of peak memory usage in real short-read sequencing samples. To extend SemiBin2 to long-read data, we also propose ensemble-based DBSCAN clustering algorithm, resulting in 13.1-26.3% more high-quality genomes than the second best binner for long-read data. AVAILABILITY AND IMPLEMENTATION: SemiBin2 is available as open source software at https://github.com/BigDataBiology/SemiBin/ and the analysis scripts used in the study can be found at https://github.com/BigDataBiology/SemiBin2_benchmark.

Early-initiated childhood reading for pleasure: associations with better cognitive performance, mental well-being and brain structure in young adolescence.

BACKGROUND: Childhood is a crucial neurodevelopmental period. We investigated whether childhood reading for pleasure (RfP) was related to young adolescent assessments of cognition, mental health, and brain structure. METHODS: We conducted a cross-sectional and longitudinal study in a large-scale US national cohort (10 000 + young adolescents), using the well-established linear mixed model and structural equation methods for twin study, longitudinal and mediation analyses. A 2-sample Mendelian randomization (MR) analysis for potential causal inference was also performed. Important factors including socio-economic status were controlled. RESULTS: Early-initiated long-standing childhood RfP (early RfP) was highly positively correlated with performance on cognitive tests and significantly negatively correlated with mental health problem scores of young adolescents. These participants with higher early RfP scores exhibited moderately larger total brain cortical areas and volumes, with increased regions including the temporal, frontal, insula, supramarginal; left angular, para-hippocampal; right middle-occipital, anterior-cingulate, orbital areas; and subcortical ventral-diencephalon and thalamus. These brain structures were significantly related to their cognitive and mental health scores, and displayed significant mediation effects. Early RfP was longitudinally associated with higher crystallized cognition and lower attention symptoms at follow-up. Approximately 12 h/week of youth regular RfP was cognitively optimal. We further observed a moderately significant heritability of early RfP, with considerable contribution from environments. MR analysis revealed beneficial causal associations of early RfP with adult cognitive performance and left superior temporal structure. CONCLUSIONS: These findings, for the first time, revealed the important relationships of early RfP with subsequent brain and cognitive development and mental well-being.

mBodyMap: a curated database for microbes across human body and their associations with health and diseases.

mBodyMap is a curated database for microbes across the human body and their associations with health and diseases. Its primary aim is to promote the reusability of human-associated metagenomic data and assist with the identification of disease-associated microbes by consistently annotating the microbial contents of collected samples using state-of-the-art toolsets and manually curating the meta-data of corresponding human hosts. mBodyMap organizes collected samples based on their association with human diseases and body sites to enable cross-dataset integration and comparison. To help users find microbes of interest and visualize and compare their distributions and abundances/prevalence within different body sites and various diseases, the mBodyMap database is equipped with an intuitive interface and extensive graphical representations of the collected data. So far, it contains a total of 63 148 runs, including 14 401 metagenomes and 48 747 amplicons related to health and 56 human diseases, from within 22 human body sites across 136 projects. Also available in the database are pre-computed abundances and prevalence of 6247 species (belonging to 1645 genera) stratified by body sites and diseases. mBodyMap can be accessed at: https://mbodymap.microbiome.cloud.

GMrepo v2: a curated human gut microbiome database with special focus on disease markers and cross-dataset comparison.

GMrepo (data repository for Gut Microbiota) is a database of curated and consistently annotated human gut metagenomes. Its main purposes are to increase the reusability and accessibility of human gut metagenomic data, and enable cross-project and phenotype comparisons. To achieve these goals, we performed manual curation on the meta-data and organized the datasets in a phenotype-centric manner. GMrepo v2 contains 353 projects and 71,642 runs/samples, which are significantly increased from the previous version. Among these runs/samples, 45,111 and 26,531 were obtained by 16S rRNA amplicon and whole-genome metagenomics sequencing, respectively. We also increased the number of phenotypes from 92 to 133. In addition, we introduced disease-marker identification and cross-project/phenotype comparison. We first identified disease markers between two phenotypes (e.g. health versus diseases) on a per-project basis for selected projects. We then compared the identified markers for each phenotype pair across datasets to facilitate the identification of consistent microbial markers across datasets. Finally, we provided a marker-centric view to allow users to check if a marker has different trends in different diseases. So far, GMrepo includes 592 marker taxa (350 species and 242 genera) for 47 phenotype pairs, identified from 83 selected projects. GMrepo v2 is freely available at: https://gmrepo.humangut.info.

Unnatural Direct Interspecies Electron Transfer Enabled by Living Cell-Cell Click Chemistry.

Direct interspecies electron transfer (DIET) is essential for maintaining the function and stability of anaerobic microbial consortia. However, only limited natural DIET modes have been identified and DIET engineering remains highly challenging. In this study, an unnatural DIET between Shewanella oneidensis MR-1 (SO, electron donating partner) and Rhodopseudomonas palustris (RP, electron accepting partner) was artificially established by a facile living cell-cell click chemistry strategy. By introducing alkyne- or azide-modified monosaccharides onto the cell outer surface of the target species, precise covalent connections between different species in high proximity were realized through a fast click chemistry reaction. Remarkably, upon covalent connection, outer cell surface C-type cytochromes mediated DIET between SO and RP was achieved and identified, although this was never realized naturally. Moreover, this connection directly shifted the natural H2 mediated interspecies electron transfer (MIET) to DIET between SO and RP, which delivered superior interspecies electron exchange efficiency. Therefore, this work demonstrated a naturally unachievable DIET and an unprecedented MIET shift to DIET accomplished by cell-cell distance engineering, offering an efficient and versatile solution for DIET engineering, which extends our understanding of DIET and opens up new avenues for DIET exploration and applications.

Compared to histamine-2 receptor antagonist, proton pump inhibitor induces stronger oral-to-gut microbial transmission and gut microbiome alterations: a randomised controlled trial.

OBJECTIVE: We aim to compare the effects of proton pump inhibitors (PPIs) and histamine-2 receptor antagonists (H2RAs) on the gut microbiota through longitudinal analysis. DESIGN: Healthy volunteers were randomly assigned to receive either PPI (n=23) or H2RA (n=26) daily for seven consecutive days. We collected oral (saliva) and faecal samples before and after the intervention for metagenomic next-generation sequencing. We analysed intervention-induced alterations in the oral and gut microbiome including microbial abundance and growth rates, oral-to-gut transmissions, and compared differences between the PPI and H2RA groups. RESULTS: Both interventions disrupted the gut microbiota, with PPIs demonstrating more pronounced effects. PPI usage led to a significantly higher extent of oral-to-gut transmission and promoted the growth of specific oral microbes in the gut. This led to a significant increase in both the number and total abundance of oral species present in the gut, including the identification of known disease-associated species like Fusobacterium nucleatum and Streptococcus anginosus. Overall, gut microbiome-based machine learning classifiers could accurately distinguish PPI from non-PPI users, achieving an area under the receiver operating characteristic curve (AUROC) of 0.924, in contrast to an AUROC of 0.509 for H2RA versus non-H2RA users. CONCLUSION: Our study provides evidence that PPIs have a greater impact on the gut microbiome and oral-to-gut transmission than H2RAs, shedding light on the mechanism underlying the higher risk of certain diseases associated with prolonged PPI use. TRIAL REGISTRATION NUMBER: ChiCTR2300072310.

Neurodevelopmental risk and adaptation as a model for comorbidity among internalizing and externalizing disorders: genomics and cell-specific expression enriched morphometric study.

BACKGROUND: Comorbidity is the rule rather than the exception for childhood and adolescent onset mental disorders, but we cannot predict its occurrence and do not know the neural mechanisms underlying comorbidity. We investigate if the effects of comorbid internalizing and externalizing disorders on anatomical differences represent a simple aggregate of the effects on each disorder and if these comorbidity-associated cortical surface differences relate to a distinct genetic underpinning. METHODS: We studied the cortical surface area (SA) and thickness (CT) of 11,878 preadolescents (9-10 years) from the Adolescent Brain and Cognitive Development Study. Linear mixed models were implemented in comparative and association analyses among internalizing (dysthymia, major depressive disorder, disruptive mood dysregulation disorder, agoraphobia, panic disorder, specific phobia, separation anxiety disorder, social anxiety disorder, generalized anxiety disorder, post-traumatic stress disorder), externalizing (attention-deficit/hyperactivity disorder, oppositional defiant disorder, conduct disorder) diagnostic groups, a group with comorbidity of the two and a healthy control group. Genome-wide association analysis (GWAS) and cell type specificity analysis were performed on 4468 unrelated European participants from this cohort. RESULTS: Smaller cortical surface area but higher thickness was noted across patient groups when compared to controls. Children with comorbid internalizing and externalizing disorders had more pronounced areal reduction than those without comorbidity, indicating an additive burden. In contrast, cortical thickness had a non-linear effect with comorbidity: the comorbid group had no significant CT differences, while those patient groups without comorbidity had significantly higher thickness compare to healthy controls. Distinct biological pathways were implicated in regional SA and CT differences. Specifically, CT differences were associated with immune-related processes implicating astrocytes and oligodendrocytes, while SA-related differences related mainly to inhibitory neurons. CONCLUSION: The emergence of comorbidity across distinct clusters of psychopathology is unlikely to be due to a simple additive neurobiological effect alone. Distinct developmental risk moderated by immune-related adaptation processes, with unique genetic and cell-specific factors, may contribute to underlying SA and CT differences. Children with the highest risk but lowest resilience, both captured in their developmental morphometry, may develop a comorbid illness pattern.

nMAGMA: a network-enhanced method for inferring risk genes from GWAS summary statistics and its application to schizophrenia.

MOTIVATION: Annotating genetic variants from summary statistics of genome-wide association studies (GWAS) is crucial for predicting risk genes of various disorders. The multimarker analysis of genomic annotation (MAGMA) is one of the most popular tools for this purpose, where MAGMA aggregates signals of single nucleotide polymorphisms (SNPs) to their nearby genes. In biology, SNPs may also affect genes that are far away in the genome, thus missed by MAGMA. Although different upgrades of MAGMA have been proposed to extend gene-wise variant annotations with more information (e.g. Hi-C or eQTL), the regulatory relationships among genes and the tissue specificity of signals have not been taken into account. RESULTS: We propose a new approach, namely network-enhanced MAGMA (nMAGMA), for gene-wise annotation of variants from GWAS summary statistics. Compared with MAGMA and H-MAGMA, nMAGMA significantly extends the lists of genes that can be annotated to SNPs by integrating local signals, long-range regulation signals (i.e. interactions between distal DNA elements), and tissue-specific gene networks. When applied to schizophrenia (SCZ), nMAGMA is able to detect more risk genes (217% more than MAGMA and 57% more than H-MAGMA) that are involved in SCZ compared with MAGMA and H-MAGMA, and more of nMAGMA results can be validated with known SCZ risk genes. Some disease-related functions (e.g. the ATPase pathway in Cortex) are also uncovered in nMAGMA but not in MAGMA or H-MAGMA. Moreover, nMAGMA provides tissue-specific risk signals, which are useful for understanding disorders with multitissue origins.

Identifying age-specific gene signatures of the human cerebral cortex with joint analysis of transcriptomes and functional connectomes.

The human cerebral cortex undergoes profound structural and functional dynamic variations across the lifespan, whereas the underlying molecular mechanisms remain unclear. Here, with a novel method transcriptome-connectome correlation analysis (TCA), which integrates the brain functional magnetic resonance images and region-specific transcriptomes, we identify age-specific cortex (ASC) gene signatures for adolescence, early adulthood and late adulthood. The ASC gene signatures are significantly correlated with the cortical thickness (P-value <2.00e-3) and myelination (P-value <1.00e-3), two key brain structural features that vary in accordance with brain development. In addition to the molecular underpinning of age-related brain functions, the ASC gene signatures allow delineation of the molecular mechanisms of neuropsychiatric disorders, such as the regulation between ARNT2 and its target gene ETF1 involved in Schizophrenia. We further validate the ASC gene signatures with published gene sets associated with the adult cortex, and confirm the robustness of TCA on other brain image datasets. Availability: All scripts are written in R. Scripts for the TCA method and related statistics result can be freely accessed at https://github.com/Soulnature/TCA. Additional data related to this paper may be requested from the authors.

CITEdb: a manually curated database of cell-cell interactions in human.

MOTIVATION: The interactions among various types of cells play critical roles in cell functions and the maintenance of the entire organism. While cell-cell interactions are traditionally revealed from experimental studies, recent developments in single-cell technologies combined with data mining methods have enabled computational prediction of cell-cell interactions, which have broadened our understanding of how cells work together, and have important implications in therapeutic interventions targeting cell-cell interactions for cancers and other diseases. Despite the importance, to our knowledge, there is no database for systematic documentation of high-quality cell-cell interactions at the cell type level, which hinders the development of computational approaches to identify cell-cell interactions. RESULTS: We develop a publicly accessible database, CITEdb (Cell-cell InTEraction database, https://citedb.cn/), which not only facilitates interactive exploration of cell-cell interactions in specific physiological contexts (e.g. a disease or an organ) but also provides a benchmark dataset to interpret and evaluate computationally derived cell-cell interactions from different tools. CITEdb contains 728 pairs of cell-cell interactions in human that are manually curated. Each interaction is equipped with structured annotations including the physiological context, the ligand-receptor pairs that mediate the interaction, etc. Our database provides a web interface to search, visualize and download cell-cell interactions. Users can search for cell-cell interactions by selecting the physiological context of interest or specific cell types involved. CITEdb is the first attempt to catalogue cell-cell interactions at the cell type level, which is beneficial to both experimental, computational and clinical studies of cell-cell interactions. AVAILABILITY AND IMPLEMENTATION: CITEdb is freely available at https://citedb.cn/ and the R package implementing benchmark is available at https://github.com/shanny01/benchmark. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.