The development and progression of cancer arise due to the accumulation of mutations in driver genes. Correctly identifying the driver genes that lead to cancer development can significantly assist the drug design, cancer diagnosis and treatment. Most computer methods detect cancer drivers based on gene-gene networks by assuming that driver genes tend to work together, form protein complexes and enrich pathways. However, they ignore that microribonucleic acid (RNAs; miRNAs) regulate the expressions of their targeted genes and are related to human diseases. In this work, we propose a graph convolution network (GCN) approach called GM-GCN to identify the cancer driver genes based on a gene-miRNA network. First, we constructed a gene-miRNA network, where the nodes are miRNAs and their targeted genes. The edges connecting miRNA and genes indicate the regulatory relationship between miRNAs and genes. We prepared initial attributes for miRNA and genes according to their biological properties and used a GCN model to learn the gene feature representations in the network by aggregating the features of their neighboring miRNA nodes. And then, the learned features were passed through a 1D convolution module for feature dimensionality change. We employed the learned and original gene features to optimize model parameters. Finally, the gene features learned from the network and the initial input gene features were fed into a logistic regression model to predict whether a gene is a driver gene. We applied our model and state-of-the-art methods to predict cancer drivers for pan-cancer and individual cancer types. Experimental results show that our model performs well in terms of the area under the receiver operating characteristic curve and the area under the precision-recall curve compared to state-of-the-art methods that work on gene networks. The GM-GCN is freely available via https://github.com/weiba/GM-GCN.

Analysis of cell-cell communication (CCC) in the tumor micro-environment helps decipher the underlying mechanism of cancer progression and drug tolerance. Currently, single-cell RNA-Seq data are available on a large scale, providing an unprecedented opportunity to predict cellular communications. There have been many achievements and applications in inferring cell-cell communication based on the known interactions between molecules, such as ligands, receptors and extracellular matrix. However, the prior information is not quite adequate and only involves a fraction of cellular communications, producing many false-positive or false-negative results. To this end, we propose an improved hierarchical variational autoencoder (HiVAE) based model to fully use single-cell RNA-seq data for automatically estimating CCC. Specifically, the HiVAE model is used to learn the potential representation of cells on known ligand-receptor genes and all genes in single-cell RNA-seq data, respectively, which are then utilized for cascade integration. Subsequently, transfer entropy is employed to measure the transmission of information flow between two cells based on the learned representations, which are regarded as directed communication relationships. Experiments are conducted on single-cell RNA-seq data of the human skin disease dataset and the melanoma dataset, respectively. Results show that the HiVAE model is effective in learning cell representations, and transfer entropy could be used to estimate the communication scores between cell types.

Long-range ribonucleic acid (RNA)-RNA interactions (RRI) are prevalent in positive-strand RNA viruses, including Beta-coronaviruses, and these take part in regulatory roles, including the regulation of sub-genomic RNA production rates. Crosslinking of interacting RNAs and short read-based deep sequencing of resulting RNA-RNA hybrids have shown that these long-range structures exist in severe acute respiratory syndrome coronavirus (SARS-CoV)-2 on both genomic and sub-genomic levels and in dynamic topologies. Furthermore, co-evolution of coronaviruses with their hosts is navigated by genetic variations made possible by its large genome, high recombination frequency and a high mutation rate. SARS-CoV-2's mutations are known to occur spontaneously during replication, and thousands of aggregate mutations have been reported since the emergence of the virus. Although many long-range RRIs have been experimentally identified using high-throughput methods for the wild-type SARS-CoV-2 strain, evolutionary trajectory of these RRIs across variants, impact of mutations on RRIs and interaction of SARS-CoV-2 RNAs with the host have been largely open questions in the field. In this review, we summarize recent computational tools and experimental methods that have been enabling the mapping of RRIs in viral genomes, with a specific focus on SARS-CoV-2. We also present available informatics resources to navigate the RRI maps and shed light on the impact of mutations on the RRI space in viral genomes. Investigating the evolution of long-range RNA interactions and that of virus-host interactions can contribute to the understanding of new and emerging variants as well as aid in developing improved RNA therapeutics critical for combating future outbreaks.

As a dynamical system, complex disease always has a sudden state transition at the tipping point, which is the result of the long-term accumulation of abnormal regulations. This paper proposes a novel approach to detect the early-warning signals of influenza A (H3N2 and H1N1) outbreaks by dysregulated dynamic network biomarkers (dysregulated DNBs) for individuals. The results of cross-validation show that our approach can detect early-warning signals before the symptom appears successfully. Unlike the traditional DNBs, our dysregulated DNBs are anchored and very few, which is essential for disease early diagnosis in clinical practice. Moreover, the genes of dysregulated DNBs are significantly enriched in the influenza-related pathways. The source code of this paper can be freely downloaded from https://github.com/YanhaoHuo/dysregulated-DNBs.git.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first detected in December 2019. As of mid-2021, the delta variant was the primary type; however, in January 2022, the omicron (BA.1) variant rapidly spread and became the dominant type in the United States. In June 2022, its subvariants surpassed previous variants in different temporal and spatial situations. To investigate the high transmissibility of omicron variants, we assessed the complex of spike protein 1 receptor-binding domain (S1RBD) and human angiotensin-converting enzyme 2 (hACE2) from the Protein Data Bank (6m0j, 7a91, 7mjn, 7v80, 7v84, 7v8b, 7wbl and 7xo9) and directly mutated specific amino acids to simulate several variants, including variants of concern (alpha, beta, gamma, delta), variants of interest (delta plus, epsilon, lambda, mu, mu without R346K) and omicron variants (BA.1, BA.2, BA.2.12.1, BA.4, BA.5). Molecular dynamics (MD) simulations for 100 ns under physiological conditions were then performed. We found that the omicron S1RBD-hACE2 complexes become more compact with increases in hydrogen-bond interactions at the interface, which is related to the transmissibility of SARS-CoV-2. Moreover, the relaxation time of hydrogen bonds is relatively short among the omicron variants, which implies that the interface conformation alterations are fast. From the molecular perspective, PHE486 and TYR501 in omicron S1RBDs need to involve hydrogen bonds and hydrophobic interactions on the interface. Our study provides structural features of the dominant variants that explain the evolution trend and their increased contagiousness and could thus also shed light on future variant changes.