Resolving Cell Fate Decisions during Somatic Cell Reprogramming by Single-Cell RNA-Seq.

Somatic cells can be reprogrammed into induced pluripotent stem cells (iPSCs), which is a highly heterogeneous process. Here we report the cell fate continuum during somatic cell reprogramming at single-cell resolution. We first develop SOT to analyze cell fate continuum from Oct4/Sox2/Klf4- or OSK-mediated reprogramming and show that cells bifurcate into two categories, reprogramming potential (RP) or non-reprogramming (NR). We further show that Klf4 contributes to Cd34+/Fxyd5+/Psca+ keratinocyte-like NR fate and that IFN-γ impedes the final transition to chimera-competent pluripotency along the RP cells. We analyze more than 150,000 single cells from both OSK and chemical reprograming and identify additional NR/RP bifurcation points. Our work reveals a generic bifurcation model for cell fate decisions during somatic cell reprogramming that may be applicable to other systems and inspire further improvements for reprogramming.

Regeneration of the human segmentation clock in somitoids in vitro.

Each vertebrate species appears to have a unique timing mechanism for forming somites along the vertebral column, and the process in human remains poorly understood at the molecular level due to technical and ethical limitations. Here, we report the reconstitution of human segmentation clock by direct reprogramming. We first reprogrammed human urine epithelial cells to a presomitic mesoderm (PSM) state capable of long-term self-renewal and formation of somitoids with an anterior-to-posterior axis. By inserting the RNA reporter Pepper into HES7 and MESP2 loci of these iPSM cells, we show that both transcripts oscillate in the resulting somitoids at ~5 h/cycle. GFP-tagged endogenous HES7 protein moves along the anterior-to-posterior axis during somitoid formation. The geo-sequencing analysis further confirmed anterior-to-posterior polarity and revealed the localized expression of WNT, BMP, FGF, and RA signaling molecules and HOXA-D family members. Our study demonstrates the direct reconstitution of human segmentation clock from somatic cells, which may allow future dissection of the mechanism and components of such a clock and aid regenerative medicine.

A survey of generative AI for de novo drug design: new frontiers in molecule and protein generation.

Artificial intelligence (AI)-driven methods can vastly improve the historically costly drug design process, with various generative models already in widespread use. Generative models for de novo drug design, in particular, focus on the creation of novel biological compounds entirely from scratch, representing a promising future direction. Rapid development in the field, combined with the inherent complexity of the drug design process, creates a difficult landscape for new researchers to enter. In this survey, we organize de novo drug design into two overarching themes: small molecule and protein generation. Within each theme, we identify a variety of subtasks and applications, highlighting important datasets, benchmarks, and model architectures and comparing the performance of top models. We take a broad approach to AI-driven drug design, allowing for both micro-level comparisons of various methods within each subtask and macro-level observations across different fields. We discuss parallel challenges and approaches between the two applications and highlight future directions for AI-driven de novo drug design as a whole. An organized repository of all covered sources is available at https://github.com/gersteinlab/GenAI4Drug.

pathMap: a path-based mapping tool for long noisy reads with high sensitivity.

With the rapid development of single-molecule sequencing (SMS) technologies, the output read length is continuously increasing. Mapping such reads onto a reference genome is one of the most fundamental tasks in sequence analysis. Mapping sensitivity is becoming a major concern since high sensitivity can detect more aligned regions on the reference and obtain more aligned bases, which are useful for downstream analysis. In this study, we present pathMap, a novel k-mer graph-based mapper that is specifically designed for mapping SMS reads with high sensitivity. By viewing the alignment chain as a path containing as many anchors as possible in the matched k-mer graph, pathMap treats chaining as a path selection problem in the directed graph. pathMap iteratively searches the longest path in the remaining nodes; more candidate chains with high quality can be effectively detected and aligned. Compared to other state-of-the-art mapping methods such as minimap2 and Winnowmap2, experiment results on simulated and real-life datasets demonstrate that pathMap obtains the number of mapped chains at least 11.50% more than its closest competitor and increases the mapping sensitivity by 17.28% and 13.84% of bases over the next-best mapper for Pacific Biosciences and Oxford Nanopore sequencing data, respectively. In addition, pathMap is more robust to sequence errors and more sensitive to species- and strain-specific identification of pathogens using MinION reads.

Unraveling the assembly mechanism of SADS-CoV virus nucleocapsid protein: insights from RNA binding, dimerization, and epitope diversity profiling.

The swine acute diarrhea syndrome coronavirus (SADS-CoV) has caused significant disruptions in porcine breeding and raised concerns about potential human infection. The nucleocapsid (N) protein of SADS-CoV plays a vital role in viral assembly and replication, but its structure and functions remain poorly understood. This study utilized biochemistry, X-ray crystallography, and immunization techniques to investigate the N protein's structure and function in SADS-CoV. Our findings revealed distinct domains within the N protein, including an RNA-binding domain, two disordered domains, and a dimerization domain. Through biochemical assays, we confirmed that the N-terminal domain functions as an RNA-binding domain, and the C-terminal domain is involved in dimerization, with the crystal structure analysis providing visual evidence of dimer formation. Immunization experiments demonstrated that the disordered domain 2 elicited a significant antibody response. These identified domains and their interactions are crucial for viral assembly. This comprehensive understanding of the N protein in SADS-CoV enhances our knowledge of its assembly and replication mechanisms, enabling the development of targeted interventions and therapeutic strategies. IMPORTANCE: SADS-CoV is a porcine coronavirus that originated from a bat HKU2-related coronavirus. It causes devastating swine diseases and poses a high risk of spillover to humans. The coronavirus N protein, as the most abundant viral protein in infected cells, likely plays a key role in viral assembly and replication. However, the structure and function of this protein remain unclear. Therefore, this study employed a combination of biochemistry and X-ray crystallography to uncover distinct structural domains in the N protein, including RNA-binding domains, two disordered domains, and dimerization domains. Additionally, we made the novel discovery that the disordered domain elicited a significant antibody response. These findings provide new insights into the structure and functions of the SADS-CoV N protein, which have important implications for future studies on SADS-CoV diagnosis, as well as the development of vaccines and anti-viral drugs.

Proteomic Dynamics of Breast Cancer Cell Lines Identifies Potential Therapeutic Protein Targets.

Treatment and relevant targets for breast cancer (BC) remain limited, especially for triple-negative BC (TNBC). We identified 6091 proteins of 76 human BC cell lines using data-independent acquisition (DIA). Integrating our proteomic findings with prior multi-omics datasets, we found that including proteomics data improved drug sensitivity predictions and provided insights into the mechanisms of action. We subsequently profiled the proteomic changes in nine cell lines (five TNBC and four non-TNBC) treated with EGFR/AKT/mTOR inhibitors. In TNBC, metabolism pathways were dysregulated after EGFR/mTOR inhibitor treatment, while RNA modification and cell cycle pathways were affected by AKT inhibitor. This systematic multi-omics and in-depth analysis of the proteome of BC cells can help prioritize potential therapeutic targets and provide insights into adaptive resistance in TNBC.

The oxidized phospholipid PGPC impairs endothelial function by promoting endothelial cell ferroptosis via FABP3.

Ferroptosis is a novel cell death mechanism that is mediated by iron-dependent lipid peroxidation. It may be involved in atherosclerosis development. Products of phospholipid oxidation play a key role in atherosclerosis. 1-palmitoyl-2-glutaroyl-sn-glycero-3-phosphocholine (PGPC) is a phospholipid oxidation product present in atherosclerotic lesions. It remains unclear whether PGPC causes atherosclerosis by inducing endothelial cell ferroptosis. In this study, human umbilical vein endothelial cells (HUVECs) were treated with PGPC. Intracellular levels of ferrous iron, lipid peroxidation, superoxide anions (O2•-), and glutathione were detected, and expression of fatty acid binding protein-3 (FABP3), glutathione peroxidase 4 (GPX4), and CD36 were measured. Additionally, the mitochondrial membrane potential (MMP) was determined. Aortas from C57BL6 mice were isolated for vasodilation testing. Results showed that PGPC increased ferrous iron levels, the production of lipid peroxidation and O2•-, and FABP3 expression. However, PGPC inhibited the expression of GPX4 and glutathione production and destroyed normal MMP. These effects were also blocked by ferrostatin-1, an inhibitor of ferroptosis. FABP3 silencing significantly reversed the effect of PGPC. Furthermore, PGPC stimulated CD36 expression. Conversely, CD36 silencing reversed the effects of PGPC, including PGPC-induced FABP3 expression. Importantly, E06, a direct inhibitor of the oxidized 1-palmitoyl-2-arachidonoyl-phosphatidylcholine IgM natural antibody, inhibited the effects of PGPC. Finally, PGPC impaired endothelium-dependent vasodilation, ferrostatin-1 or FABP3 inhibitors inhibited this impairment. Our data demonstrate that PGPC impairs endothelial function by inducing endothelial cell ferroptosis through the CD36 receptor to increase FABP3 expression. Our findings provide new insights into the mechanisms of atherosclerosis and a therapeutic target for atherosclerosis.

MRDPDA: A multi-Laplacian regularized deepFM model for predicting piRNA-disease associations.

PIWI-interacting RNAs (piRNAs) are a typical class of small non-coding RNAs, which are essential for gene regulation, genome stability and so on. Accumulating studies have revealed that piRNAs have significant potential as biomarkers and therapeutic targets for a variety of diseases. However current computational methods face the challenge in effectively capturing piRNA-disease associations (PDAs) from limited data. In this study, we propose a novel method, MRDPDA, for predicting PDAs based on limited data from multiple sources. Specifically, MRDPDA integrates a deep factorization machine (deepFM) model with regularizations derived from multiple yet limited datasets, utilizing separate Laplacians instead of a simple average similarity network. Moreover, a unified objective function to combine embedding loss about similarities is proposed to ensure that the embedding is suitable for the prediction task. In addition, a balanced benchmark dataset based on piRPheno is constructed and a deep autoencoder is applied for creating reliable negative set from the unlabeled dataset. Compared with three latest methods, MRDPDA achieves the best performance on the pirpheno dataset in terms of the five-fold cross validation test and independent test set, and case studies further demonstrate the effectiveness of MRDPDA.

Oncogenic and immunological roles of RACGAP1 in pan-cancer and its potential value in nasopharyngeal carcinoma.

A particular GTPase-activating protein called RACGAP1 is involved in apoptosis, proliferation, invasion, metastasis, and drug resistance in a variety of malignancies. Nevertheless, the role of RACGAP1 in pan-cancer was less studied, and its value of the expression and prognostic of nasopharyngeal carcinoma (NPC) has not been explored. Hence, the goal of this study was to investigate the oncogenic and immunological roles of RACGAP1 in various cancers and its potential value in NPC. We comprehensively analyzed RACGAP1 expression, prognostic value, function, methylation levels, relationship with immune cells, immune infiltration, and immunotherapy response in pan-cancer utilizing multiple databases. The results discovered that RACGAP1 expression was elevated in most cancers and suggested poor prognosis, which could be related to the involvement of RACGAP1 in various cancer-related pathways such as the cell cycle and correlated with RACGAP1 methylation levels, immune cell infiltration and reaction to immunotherapy, and chemoresistance. RACGAP1 could inhibit anti-tumor immunity and immunotherapy responses by fostering immune cell infiltration and cytotoxic T lymphocyte dysfunction. Significantly, we validated that RACGAP1 mRNA and protein were highly expressed in NPC. The Gene Expression Omnibus database revealed that elevated RACGAP1 expression was associated with shorter PFS in patients with NPC, and RACGAP1 potentially influenced cell cycle progression, DNA replication, metabolism, and immune-related pathways, resulting in the recurrence and metastasis of NPC. This study indicated that RACGAP1 could be a potential biomarker in pan-cancer and NPC.

Confining Thiolysis of Dinitrophenyl Ether to a Luminescent Metal-Organic Framework with a Large Stokes Shift for Highly Efficient Detection of Hydrogen Sulfide in Rat Brain.

Hydrogen sulfide (H2S) is a gaseous signaling molecule that regulates various physiological and pathological processes in the central nervous system. It is vital to develop an effective method to detect H2S in vivo to elucidate its critical role. However, current fluorescent probes for accurate quantification of H2S still face big challenges due to complicated fabrication, small Stokes shift, unsatisfactory selectivity, and especially delayed response time. Herein, based on simple postsynthetic modification, we present an innovative strategy by confining H2S-triggered thiolysis of dinitrophenyl (DNP) ether within a luminescent metal-organic framework (MOF) to address those issues. Due to the cleavage of the DNP moiety by H2S, the nanoprobe gives rise to a remarkable fluorescence turn-on signal with a large Stokes shift of 190 nm and also provides high selectivity to H2S against various interferents including competing biothiols. In particular, by virtue of the unique structural property of the MOF, it exhibits an ultrafast sensing ability for H2S (only 5 s). Moreover, the fluorescence enhancement efficiency displays a good linear correlation with H2S concentration in the range of 0-160 µM with a detection limit of 0.29 µM. Importantly, these superior sensing performances enable the nanoprobe to measure the basal value and monitor the change of H2S level in the rat brain.

Reducing Quantitative Uncertainty Caused by Data Processing in Untargeted Metabolomics.

Processing liquid chromatography-mass spectrometry-based metabolomics data using computational programs often introduces additional quantitative uncertainty, termed computational variation in a previous work. This work develops a computational solution to automatically recognize metabolic features with computational variation in a metabolomics data set. This tool, AVIR (short for "Accurate eValuation of alIgnment and integRation"), is a support vector machine-based machine learning strategy (https://github.com/HuanLab/AVIR). The rationale is that metabolic features with computational variation have a poor correlation between chromatographic peak area and peak height-based quantifications across the samples in a study. AVIR was trained on a set of 696 manually curated metabolic features and achieved an accuracy of 94% in a 10-fold cross-validation. When tested on various external data sets from public metabolomics repositories, AVIR demonstrated an accuracy range of 84%-97%. Finally, tested on a large-scale metabolomics study, AVIR clearly indicated features with computational variation and thus guided us to manually correct them. Our results show that 75.3% of the samples with computational variation had a relative intensity difference of over 20% after correction. This demonstrates the critical role of AVIR in reducing computational variation to improve quantitative certainty in untargeted metabolomics analysis.

Association of vaccine-specific regulatory T cells with reduced antibody response to repeated influenza vaccination.

Repeated annual influenza vaccinations have been associated with reduced vaccine-induced antibody responses. This prospective study aimed to explore the role of vaccine antigen-specific regulatory T (Treg) cells in antibody response to repeated annual influenza vaccination. We analyzed pre- and postvaccination hemagglutination inhibition (HI) titers, seroconversion rates, seroprotection rates, vaccine antigen hemagglutinin (HA)-specific Treg cells, and conventional T (Tconv) cells. We compared these parameters between vaccinees with or without vaccine-induced seroconversion. Our multivariate logistic regression revealed that prior vaccination was significantly associated with a decreased likelihood of achieving seroconversion for both H1N1(adjusted OR, 0.03; 95% CI, 0.01-0.13) and H3N2 (adjusted OR, 0.09; 95% CI, 0.03-0.30). Furthermore, individuals who received repeated vaccinations had significantly higher levels of pre-existing HA-specific Treg cells than those who did not. We also found that vaccine-induced fold-increases in HI titers and seroconversion were negatively correlated with pre-existing HA-specific Treg cells and positively correlated with the ratio of Tconv to Treg cells. Overall, our findings suggest that repeated annual influenza vaccination is associated with a lower vaccine-induced antibody response and a higher frequency of vaccine-specific Treg cells. However, a lower frequency of pre-existing Treg cells correlates with a higher postvaccination antibody response.

Single-nucleus transcriptome analysis identifies a novel FKBP5+ endothelial cell subtype involved in endothelial-to-mesenchymal transition in adipose tissue during aging.

Age-associated adipose tissue (AT) dysfunction is multifactorial and often leads to detrimental health consequences. AT is highly vascularized and endothelial cells (ECs) has been recently identified as a key regulator in the homeostasis of AT. However, the alteration of cell composition in AT during aging and the communication between endothelial cells and adipocytes remain poorly understood. In this study, we take advantage of single nucleus RNA sequencing analysis, and discovered a group of FKBP5+ ECs specifically resident in aged AT. Of interest, FKBP5+ ECs exhibited the potential for endothelial-to-mesenchymal transition (EndoMT) and exhibited a critical role in regulating adipocytes. Furthermore, lineage tracing experiments demonstrated that ECs in aged AT tend to express FKBP5 and undergo EndoMT with progressive loss of endothelial marker. This study may provide a basis for a new mechanism of microvascular ECs-induced AT dysfunction during aging.

Disruption of the interaction between caveolae and Piezo1 promotes pressure overload-induced cardiac remodeling.

Piezo1 channels are activated by mechanical stress and play a significant role in cardiac hypertrophy and fibrosis. However, the molecular mechanisms underlying Piezo1 activation on the cell membrane following pressure overload remain unclear. Caveolae are known to mitigate mechanical forces and regulate Piezo1 function. Therefore, this study aimed to investigate the interaction between caveolae and Piezo1 in the development of pressure overload-induced cardiac remodeling. We observed reduced colocalization between Piezo1 and Caveolin-3 in hypertrophic cardiomyocytes following abdominal aortic constriction and Angiotensin-II treatment, accompanied by increased Piezo1 function and expression. Furthermore, enhanced Piezo1 function was also noted upon caveolae disruption using methyl-beta-cyclodextrin (mßCD). Thus, our findings suggested that pressure overload led to Piezo1 translocation from caveolae, thereby augmenting its function and expression, which may contribute to cardiac remodeling.

Effective Bidirectional Mott-Schottky Catalysts Derived from Spent LiFePO4 Cathodes for Robust Lithium-Sulfur Batteries.

It is deemed as a tough yet profound project to comprehensively cope with a range of detrimental problems of lithium-sulfur batteries (LSBs), mainly pertaining to the shuttle effect of lithium polysulfides (LiPSs) and sluggish sulfur conversion. Herein, a Co2P-Fe2P@N-doped carbon (Co2P-Fe2P@NC) Mott-Schottky catalyst is introduced to enable bidirectionally stimulated sulfur conversion. This catalyst is prepared by simple carbothermal reduction of spent LiFePO4 cathode and LiCoO2. The experimental and theoretical calculation results indicate that thanks to unique surface/interface properties derived from the Mott-Schottky effect, full anchoring of LiPSs, mediated Li2S nucleation/dissolution, and bidirectionally expedited "solid⇌liquid⇌solid" kinetics can be harvested. Consequently, the S/Co2P-Fe2P@NC manifests high reversible capacity (1569.9 mAh g-1), superb rate response (808.9 mAh g-1 at 3C), and stable cycling (a low decay rate of 0.06% within 600 cycles at 3C). Moreover, desirable capacity (5.35 mAh cm-2) and cycle stability are still available under high sulfur loadings (4-5 mg cm-2) and lean electrolyte (8 µL mg-1) conditions. Furthermore, the as-proposed universal synthetic route can be extended to the preparation of other catalysts such as Mn2P-Fe2P@NC from spent LiFePO4 and MnO2. This work unlocks the potential of carbothermal reduction phosphating to synthesize bidirectional catalysts for robust LSBs.

MLRDFM: a multi-view Laplacian regularized DeepFM model for predicting miRNA-disease associations.

MOTIVATION: MicroRNAs (miRNAs), as critical regulators, are involved in various fundamental and vital biological processes, and their abnormalities are closely related to human diseases. Predicting disease-related miRNAs is beneficial to uncovering new biomarkers for the prevention, detection, prognosis, diagnosis and treatment of complex diseases. RESULTS: In this study, we propose a multi-view Laplacian regularized deep factorization machine (DeepFM) model, MLRDFM, to predict novel miRNA-disease associations while improving the standard DeepFM. Specifically, MLRDFM improves DeepFM from two aspects: first, MLRDFM takes the relationships among items into consideration by regularizing their embedding features via their similarity-based Laplacians. In this study, miRNA Laplacian regularization integrates four types of miRNA similarity, while disease Laplacian regularization integrates two types of disease similarity. Second, to judiciously train our model, Laplacian eigenmaps are utilized to initialize the weights in the dense embedding layer. The experimental results on the latest HMDD v3.2 dataset show that MLRDFM improves the performance and reduces the overfitting phenomenon of DeepFM. Besides, MLRDFM is greatly superior to the state-of-the-art models in miRNA-disease association prediction in terms of different evaluation metrics with the 5-fold cross-validation. Furthermore, case studies further demonstrate the effectiveness of MLRDFM.

Predicting drug-drug interactions by graph convolutional network with multi-kernel.

Drug repositioning is proposed to find novel usages for existing drugs. Among many types of drug repositioning approaches, predicting drug-drug interactions (DDIs) helps explore the pharmacological functions of drugs and achieves potential drugs for novel treatments. A number of models have been applied to predict DDIs. The DDI network, which is constructed from the known DDIs, is a common part in many of the existing methods. However, the functions of DDIs are different, and thus integrating them in a single DDI graph may overlook some useful information. We propose a graph convolutional network with multi-kernel (GCNMK) to predict potential DDIs. GCNMK adopts two DDI graph kernels for the graph convolutional layers, namely, increased DDI graph consisting of 'increase'-related DDIs and decreased DDI graph consisting of 'decrease'-related DDIs. The learned drug features are fed into a block with three fully connected layers for the DDI prediction. We compare various types of drug features, whereas the target feature of drugs outperforms all other types of features and their concatenated features. In comparison with three different DDI prediction methods, our proposed GCNMK achieves the best performance in terms of area under receiver operating characteristic curve and area under precision-recall curve. In case studies, we identify the top 20 potential DDIs from all unknown DDIs, and the top 10 potential DDIs from the unknown DDIs among breast, colorectal and lung neoplasms-related drugs. Most of them have evidence to support the existence of their interactions. fangxiang.wu@usask.ca.

DeepLncLoc: a deep learning framework for long non-coding RNA subcellular localization prediction based on subsequence embedding.

Long non-coding RNAs (lncRNAs) are a class of RNA molecules with more than 200 nucleotides. A growing amount of evidence reveals that subcellular localization of lncRNAs can provide valuable insights into their biological functions. Existing computational methods for predicting lncRNA subcellular localization use k-mer features to encode lncRNA sequences. However, the sequence order information is lost by using only k-mer features. We proposed a deep learning framework, DeepLncLoc, to predict lncRNA subcellular localization. In DeepLncLoc, we introduced a new subsequence embedding method that keeps the order information of lncRNA sequences. The subsequence embedding method first divides a sequence into some consecutive subsequences and then extracts the patterns of each subsequence, last combines these patterns to obtain a complete representation of the lncRNA sequence. After that, a text convolutional neural network is employed to learn high-level features and perform the prediction task. Compared with traditional machine learning models, popular representation methods and existing predictors, DeepLncLoc achieved better performance, which shows that DeepLncLoc could effectively predict lncRNA subcellular localization. Our study not only presented a novel computational model for predicting lncRNA subcellular localization but also introduced a new subsequence embedding method which is expected to be applied in other sequence-based prediction tasks. The DeepLncLoc web server is freely accessible at http://bioinformatics.csu.edu.cn/DeepLncLoc/, and source code and datasets can be downloaded from https://github.com/CSUBioGroup/DeepLncLoc.

HyMM: hybrid method for disease-gene prediction by integrating multiscale module structure.

MOTIVATION: Identifying disease-related genes is an important issue in computational biology. Module structure widely exists in biomolecule networks, and complex diseases are usually thought to be caused by perturbations of local neighborhoods in the networks, which can provide useful insights for the study of disease-related genes. However, the mining and effective utilization of the module structure is still challenging in such issues as a disease gene prediction. RESULTS: We propose a hybrid disease-gene prediction method integrating multiscale module structure (HyMM), which can utilize multiscale information from local to global structure to more effectively predict disease-related genes. HyMM extracts module partitions from local to global scales by multiscale modularity optimization with exponential sampling, and estimates the disease relatedness of genes in partitions by the abundance of disease-related genes within modules. Then, a probabilistic model for integration of gene rankings is designed in order to integrate multiple predictions derived from multiscale module partitions and network propagation, and a parameter estimation strategy based on functional information is proposed to further enhance HyMM's predictive power. By a series of experiments, we reveal the importance of module partitions at different scales, and verify the stable and good performance of HyMM compared with eight other state-of-the-arts and its further performance improvement derived from the parameter estimation. CONCLUSIONS: The results confirm that HyMM is an effective framework for integrating multiscale module structure to enhance the ability to predict disease-related genes, which may provide useful insights for the study of the multiscale module structure and its application in such issues as a disease-gene prediction.

Biomedical data, computational methods and tools for evaluating disease-disease associations.

In recent decades, exploring potential relationships between diseases has been an active research field. With the rapid accumulation of disease-related biomedical data, a lot of computational methods and tools/platforms have been developed to reveal intrinsic relationship between diseases, which can provide useful insights to the study of complex diseases, e.g. understanding molecular mechanisms of diseases and discovering new treatment of diseases. Human complex diseases involve both external phenotypic abnormalities and complex internal molecular mechanisms in organisms. Computational methods with different types of biomedical data from phenotype to genotype can evaluate disease-disease associations at different levels, providing a comprehensive perspective for understanding diseases. In this review, available biomedical data and databases for evaluating disease-disease associations are first summarized. Then, existing computational methods for disease-disease associations are reviewed and classified into five groups in terms of the usages of biomedical data, including disease semantic-based, phenotype-based, function-based, representation learning-based and text mining-based methods. Further, we summarize software tools/platforms for computation and analysis of disease-disease associations. Finally, we give a discussion and summary on the research of disease-disease associations. This review provides a systematic overview for current disease association research, which could promote the development and applications of computational methods and tools/platforms for disease-disease associations.