CpG-island-based annotation and analysis of human housekeeping genes

A Review of the Application of Spatial Transcriptomics in Neuroscience

Since spatial transcriptomics can locate and distinguish the gene expression of functional genes in special regions and tissue, it is important for us to investigate the brain development, the development mechanism of brain diseases, and the relationship between brain structure and function in Neuroscience (or Brain science). While previous studies have introduced the crucial spatial transcriptomic techniques and data analysis methods, there are few studies to comprehensively overview the key methods, data resources, and technological applications of spatial transcriptomics in Neuroscience. For these reasons, we first investigate several common spatial transcriptomic data analysis approaches and data resources. Second, we introduce the applications of the spatial transcriptomic data analysis approaches in Neuroscience. Third, we summarize the integrating spatial transcriptomics with other technologies in Neuroscience. Finally, we discuss the challenges and future research directions of spatial transcriptomics in Neuroscience.

Inference of Coalescence Times and Variant Ages Using Convolutional Neural Networks

Accurate inference of the time to the most recent common ancestor (TMRCA) between pairs of individuals and of the age of genomic variants is key in several population genetic analyses. We developed a likelihood-free approach, called CoalNN, which uses a convolutional neural network to predict pairwise TMRCAs and allele ages from sequencing or SNP array data. CoalNN is trained through simulation and can be adapted to varying parameters, such as demographic history, using transfer learning. Across several simulated scenarios, CoalNN matched or outperformed the accuracy of model-based approaches for pairwise TMRCA and allele age prediction. We applied CoalNN to settings for which model-based approaches are under-developed and performed analyses to gain insights into the set of features it uses to perform TMRCA prediction. We next used CoalNN to analyze 2,504 samples from 26 populations in the 1,000 Genome Project data set, inferring the age of ∼80 million variants. We observed substantial variation across populations and for variants predicted to be pathogenic, reflecting heterogeneous demographic histories and the action of negative selection. We used CoalNN's predicted allele ages to construct genome-wide annotations capturing the signature of past negative selection. We performed LD-score regression analysis of heritability using summary association statistics from 63 independent complex traits and diseases (average N=314k), observing increased annotation-specific effects on heritability compared to a previous allele age annotation. These results highlight the effectiveness of using likelihood-free, simulation-trained models to infer properties of gene genealogies in large genomic data sets.

Discovering hematoma-stimulated circuits for secondary brain injury after intraventricular hemorrhage by spatial transcriptome analysis

Introduction

Central nervous system (CNS) diseases, such as neurodegenerative disorders and brain diseases caused by acute injuries, are important, yet challenging to study due to disease lesion locations and other complexities.

Methods

Utilizing the powerful method of spatial transcriptome analysis together with novel algorithms we developed for the study, we report here for the first time a 3D trajectory map of gene expression changes in the brain following acute neural injury using a mouse model of intraventricular hemorrhage (IVH). IVH is a common and representative complication after various acute brain injuries with severe mortality and mobility implications.

Results

Our data identified three main 3D global pseudospace-time trajectory bundles that represent the main neural circuits from the lateral ventricle to the hippocampus and primary cortex affected by experimental IVH stimulation. Further analysis indicated a rapid response in the primary cortex, as well as a direct and integrated effect on the hippocampus after IVH stimulation.

Discussion

These results are informative for understanding the pathophysiological changes, including the spatial and temporal patterns of gene expression changes, in IVH patients after acute brain injury, strategizing more effective clinical management regimens, and developing novel bioinformatics strategies for the study of other CNS diseases. The algorithm strategies used in this study are searchable via a web service (www.combio-lezhang.online/3dstivh/home).

Locus-Specific Detection of DNA Methylation: The Advance, Challenge, and Perspective of CRISPR-Cas Assisted Biosensors

Deoxyribonucleic acid (DNA) methylation is one of the epigenetic characteristics that result in heritable and revisable phenotype changes but without sequence changes in DNA. Aberrant methylation occurring at a specific locus was reported to be associated with cancers, insulin resistance, obesity, Alzheimer's disease, Parkinson's disease, etc. Therefore, locus-specific DNA methylation can serve as a valuable biomarker for disease diagnosis and therapy. Recently, Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems are applied to develop biosensors for DNA, ribonucleic acid, proteins, and small molecules detection. Because of their highly specific binding ability and signal amplification capacity, CRISPR-Cas assisted biosensor also serve as a potential tool for locus-specific detection of DNA methylation. In this perspective, based on the detection principle, a detailed classification and comprehensive discussion of recent works about the latest advances in locus-specific detection of DNA methylation using CRISPR-Cas systems are provided. Furthermore, current challenges and future perspectives of CRISPR-based locus-specific detection of DNA methylation are outlined.

The homeodomain of Oct4 is a dimeric binder of methylated CpG elements

Oct4 is essential to maintain pluripotency and has a pivotal role in establishing the germline. Its DNA-binding POU domain was recently found to bind motifs with methylated CpG elements normally associated with epigenetic silencing. However, the mode of binding and the consequences of this capability has remained unclear. Here, we show that Oct4 binds to a compact palindromic DNA element with a methylated CpG core (CpGpal) in alternative states of pluripotency and during cellular reprogramming towards induced pluripotent stem cells (iPSCs). During cellular reprogramming, typical Oct4 bound enhancers are uniformly demethylated, with the prominent exception of the CpGpal sites where DNA methylation is often maintained. We demonstrate that Oct4 cooperatively binds the CpGpal element as a homodimer, which contrasts with the ectoderm-expressed POU factor Brn2. Indeed, binding to CpGpal is Oct4-specific as other POU factors expressed in somatic cells avoid this element. Binding assays combined with structural analyses and molecular dynamic simulations show that dimeric Oct4-binding to CpGpal is driven by the POU-homeodomain whilst the POU-specific domain is detached from DNA. Collectively, we report that Oct4 exerts parts of its regulatory function in the context of methylated DNA through a DNA recognition mechanism that solely relies on its homeodomain.

A Computer Simulation of SARS-CoV-2 Mutation Spectra for Empirical Data Characterization and Analysis

It is very important to compute the mutation spectra, and simulate the intra-host mutation processes by sequencing data, which is not only for the understanding of SARS-CoV-2 genetic mechanism, but also for epidemic prediction, vaccine, and drug design. However, the current intra-host mutation analysis algorithms are not only inaccurate, but also the simulation methods are unable to quickly and precisely predict new SARS-CoV-2 variants generated from the accumulation of mutations. Therefore, this study proposes a novel accurate strand-specific SARS-CoV-2 intra-host mutation spectra computation method, develops an efficient and fast SARS-CoV-2 intra-host mutation simulation method based on mutation spectra, and establishes an online analysis and visualization platform. Our main results include: (1) There is a significant variability in the SARS-CoV-2 intra-host mutation spectra across different lineages, with the major mutations from G- > A, G- > C, G- > U on the positive-sense strand and C- > U, C- > G, C- > A on the negative-sense strand; (2) our mutation simulation reveals the simulation sequence starts to deviate from the base content percentage of Alpha-CoV/Delta-CoV after approximately 620 mutation steps; (3) 2019-NCSS provides an easy-to-use and visualized online platform for SARS-Cov-2 online analysis and mutation simulation.

Spatiotemporal Transformer Neural Network for Time-Series Forecasting

Predicting high-dimensional short-term time-series is a difficult task due to the lack of sufficient information and the curse of dimensionality. To overcome these problems, this study proposes a novel spatiotemporal transformer neural network (STNN) for efficient prediction of short-term time-series with three major features. Firstly, the STNN can accurately and robustly predict a high-dimensional short-term time-series in a multi-step-ahead manner by exploiting high-dimensional/spatial information based on the spatiotemporal information (STI) transformation equation. Secondly, the continuous attention mechanism makes the prediction results more accurate than those of previous studies. Thirdly, we developed continuous spatial self-attention, temporal self-attention, and transformation attention mechanisms to create a bridge between effective spatial information and future temporal evolution information. Fourthly, we show that the STNN model can reconstruct the phase space of the dynamical system, which is explored in the time-series prediction. The experimental results demonstrate that the STNN significantly outperforms the existing methods on various benchmarks and real-world systems in the multi-step-ahead prediction of a short-term time-series.

Learn from antibody–drug conjugates: consideration in the future construction of peptide-drug conjugates for cancer therapy

Cancer is one of the leading causes of death worldwide due to high heterogeneity. Although chemotherapy remains the mainstay of cancer therapy, non-selective toxicity and drug resistance of mono-chemotherapy incur broad criticisms. Subsequently, various combination strategies have been developed to improve clinical efficacy, also known as cocktail therapy. However, conventional “cocktail administration” is just passable, due to the potential toxicities to normal tissues and unsatisfactory synergistic effects, especially for the combined drugs with different pharmacokinetic properties. The drug conjugates through coupling the conventional chemotherapeutics to a carrier (such as antibody and peptide) provide an alternative strategy to improve therapeutic efficacy and simultaneously reduce the unspecific toxicities, by virtue of the advantages of highly specific targeting ability and potent killing effect. Although 14 antibody–drug conjugates (ADCs) have been approved worldwide and more are being investigated in clinical trials so far, several limitations have been disclosed during clinical application. Compared with ADCs, peptide-drug conjugates (PDCs) possess several advantages, including easy industrial synthesis, low cost, high tissue penetration and fast clearance. So far, only a handful of PDCs have been approved, highlighting tremendous development potential. Herein, we discuss the progress and pitfalls in the development of ADCs and underline what can learn from ADCs for the better construction of PDCs in the future.

An integrated platform for Brucella with knowledge graph technology: From genomic analysis to epidemiological projection

Motivation:Brucella, the causative agent of brucellosis, is a global zoonotic pathogen that threatens both veterinary and human health. The main sources of brucellosis are farm animals. Importantly, the bacteria can be used for biological warfare purposes, requiring source tracking and routine surveillance in an integrated manner. Additionally, brucellosis is classified among group B infectious diseases in China and has been reported in 31 Chinese provinces to varying degrees in urban areas. From a national biosecurity perspective, research on brucellosis surveillance has garnered considerable attention and requires an integrated platform to provide researchers with easy access to genomic analysis and provide policymakers with an improved understanding of both reported patients and detected cases for the purpose of precision public health interventions.

Results: For the first time in China, we have developed a comprehensive information platform for Brucella based on dynamic visualization of the incidence (reported patients) and prevalence (detected cases) of brucellosis in mainland China. Especially, our study establishes a knowledge graph for the literature sources of Brucella data so that it can be expanded, queried, and analyzed. When similar “epidemiological comprehensive platforms” are established in the distant future, we can use knowledge graph to share its information. Additionally, we propose a software package for genomic sequence analysis. This platform provides a specialized, dynamic, and visual point-and-click interface for studying brucellosis in mainland China and improving the exploration of Brucella in the fields of bioinformatics and disease prevention for both human and veterinary medicine.

Artificial intelligence in cancer target identification and drug discovery

Artificial intelligence is an advanced method to identify novel anticancer targets and discover novel drugs from biology networks because the networks can effectively preserve and quantify the interaction between components of cell systems underlying human diseases such as cancer. Here, we review and discuss how to employ artificial intelligence approaches to identify novel anticancer targets and discover drugs. First, we describe the scope of artificial intelligence biology analysis for novel anticancer target investigations. Second, we review and discuss the basic principles and theory of commonly used network-based and machine learning-based artificial intelligence algorithms. Finally, we showcase the applications of artificial intelligence approaches in cancer target identification and drug discovery. Taken together, the artificial intelligence models have provided us with a quantitative framework to study the relationship between network characteristics and cancer, thereby leading to the identification of potential anticancer targets and the discovery of novel drug candidates.

Developing an Embedding, Koopman and Autoencoder Technologies-Based Multi-Omics Time Series Predictive Model (EKATP) for Systems Biology research

It is very important for systems biologists to predict the state of the multi-omics time series for disease occurrence and health detection. However, it is difficult to make the prediction due to the high-dimensional, nonlinear and noisy characteristics of the multi-omics time series data. For this reason, this study innovatively proposes an Embedding, Koopman and Autoencoder technologies-based multi-omics time series predictive model (EKATP) to predict the future state of a high-dimensional nonlinear multi-omics time series. We evaluate this EKATP by using a genomics time series with chaotic behavior, a proteomics time series with oscillating behavior and a metabolomics time series with flow behavior. The computational experiments demonstrate that our proposed EKATP can substantially improve the accuracy, robustness and generalizability to predict the future state of a time series for multi-omics data.

MCDB: A comprehensive curated mitotic catastrophe database for retrieval, protein sequence alignment, and target prediction

Mitotic catastrophe (MC) is a form of programmed cell death induced by mitotic process disorders, which is very important in tumor prevention, development, and drug resistance. Because rapidly increased data for MC is vigorously promoting the tumor-related biomedical and clinical study, it is urgent for us to develop a professional and comprehensive database to curate MC-related data. Mitotic Catastrophe Database (MCDB) consists of 1214 genes/proteins and 5014 compounds collected and organized from more than 8000 research articles. Also, MCDB defines the confidence level, classification criteria, and uniform naming rules for MC-related data, which greatly improves data reliability and retrieval convenience. Moreover, MCDB develops protein sequence alignment and target prediction functions. The former can be used to predict new potential MC-related genes and proteins, and the latter can facilitate the identification of potential target proteins of unknown MC-related compounds. In short, MCDB is such a proprietary, standard, and comprehensive database for MC-relate data that will facilitate the exploration of MC from chemists to biologists in the fields of medicinal chemistry, molecular biology, bioinformatics, oncology and so on. The MCDB is distributed on http://www.combio-lezhang.online/MCDB/index_html/.

2019nCoVAS: Developing the Web Service for Epidemic Transmission Prediction, Genome Analysis, and Psychological Stress Assessment for 2019-nCoV

Since the COVID-19 epidemic is still expanding around the world and poses a serious threat to human life and health, it is necessary for us to carry out epidemic transmission prediction, whole genome sequence analysis, and public psychological stress assessment for 2019-nCoV. However, transmission prediction models are insufficiently accurate and genome sequence characteristics are not clear, and it is difficult to dynamically assess the public psychological stress state under the 2019-nCoV epidemic. Therefore, this study develops a 2019nCoVAS web service (http://www.combio-lezhang.online/2019ncov/home.html) that not only offers online epidemic transmission prediction and lineage-associated underrepresented permutation (LAUP) analysis services to investigate the spreading trends and genome sequence characteristics, but also provides psychological stress assessments based on such an emotional dictionary that we built for 2019-nCoV. Finally, we discuss the shortcomings and further study of the 2019nCoVAS web service.

The "Genomic Code": DNA Pervasively Moulds Chromatin Structures Leaving no Room for "Junk"

The chromatin of the human genome was analyzed at three DNA size levels. At the first, compartment level, two "gene spaces" were found many years ago: A GC-rich, gene-rich "genome core" and a GC-poor, gene-poor "genome desert", the former corresponding to open chromatin centrally located in the interphase nucleus, the latter to closed chromatin located peripherally. This bimodality was later confirmed and extended by the discoveries (1) of LADs, the Lamina-Associated Domains, and InterLADs; (2) of two "spatial compartments", A and B, identified on the basis of chromatin interactions; and (3) of "forests and prairies" characterized by high and low CpG islands densities. Chromatin compartments were shown to be associated with the compositionally different, flat and single- or multi-peak DNA structures of the two, GC-poor and GC-rich, "super-families" of isochores. At the second, sub-compartment, level, chromatin corresponds to flat isochores and to isochore loops (due to compositional DNA gradients) that are susceptible to extrusion. Finally, at the short-sequence level, two sets of sequences, GC-poor and GC-rich, define two different nucleosome spacings, a short one and a long one. In conclusion, chromatin structures are moulded according to a "genomic code" by DNA sequences that pervade the genome and leave no room for "junk".

Review of Artificial Intelligence Applications and Algorithms for Brain Organoid Research

The human brain organoid is a miniature three-dimensional tissue culture that can simulate the structure and function of the brain in an in vitro culture environment. Although we consider that human brain organoids could be used to understand brain development and diseases, experimental models of human brain organoids are so highly variable that we apply artificial intelligence (AI) techniques to investigate the development mechanism of the human brain. Therefore, this study briefly reviewed commonly used AI applications for human brain organoid-magnetic resonance imaging, electroencephalography, and gene editing techniques, as well as related AI algorithms. Finally, we discussed the limitations, challenges, and future study direction of AI-based technology for human brain organoids.

Transcriptome-based selection and validation of optimal house-keeping genes for skin research in goats (Capra hircus)

BACKGROUND

In quantitative real-time polymerase chain reaction (qRT-PCR) experiments, accurate and reliable target gene expression results are dependent on optimal amplification of house-keeping genes (HKGs). RNA-seq technology offers a novel approach to detect new HKGs with improved stability. Goat (Capra hircus) is an economically important livestock species and plays an indispensable role in the world animal fiber and meat industry. Unfortunately, uniform and reliable HKGs for skin research have not been identified in goat. Therefore, this study seeks to identify a set of stable HKGs for the skin tissue of C. hircus using high-throughput sequencing technology.

RESULTS

Based on the transcriptome dataset of 39 goat skin tissue samples, 8 genes (SRP68, NCBP3, RRAGA, EIF4H, CTBP2, PTPRA, CNBP, and EEF2) with relatively stable expression levels were identified and selected as new candidate HKGs. Commonly used HKGs including SDHA and YWHAZ from a previous study, and 2 conventional genes (ACTB and GAPDH) were also examined. Four different experimental variables: (1) different development stages, (2) hair follicle cycle stages, (3) breeds, and (4) sampling sites were used for determination and validation. Four algorithms (geNorm, NormFinder, BestKeeper, and ΔCt method) and a comprehensive algorithm (ComprFinder, developed in-house) were used to assess the stability of each HKG. It was shown that NCBP3 + SDHA + PTPRA were more stably expressed than previously used genes in all conditions analysis, and that this combination was effective at normalizing target gene expression. Moreover, a new algorithm for comprehensive analysis, ComprFinder, was developed and released.

CONCLUSION

This study presents the first list of candidate HKGs for C. hircus skin tissues based on an RNA-seq dataset. We propose that the NCBP3 + SDHA + PTPRA combination could be regarded as a triplet set of HKGs in skin molecular biology experiments in C. hircus and other closely related species. In addition, we also encourage researchers who perform candidate HKG evaluations and who require comprehensive analysis to adopt our new algorithm, ComprFinder.

Development of an Early Prediction Model for Subarachnoid Hemorrhage With Genetic and Signaling Pathway Analysis

Subarachnoid hemorrhage (SAH) is devastating disease with high mortality, high disability rate, and poor clinical prognosis. It has drawn great attentions in both basic and clinical medicine. Therefore, it is necessary to explore the therapeutic drugs and effective targets for early prediction of SAH. Firstly, we demonstrate that LCN2 can effectively intervene or treat SAH from the perspective of cell signaling pathway. Next, three potential genes that we explored have been validated by manually reviewed experimental evidences. Finally, we turn out that the SAH early ensemble learning predictive model performs better than the classical LR, SVM, and Naïve-Bayes models.

Exploring the computational methods for protein-ligand binding site prediction

Proteins participate in various essential processes in vivo via interactions with other molecules. Identifying the residues participating in these interactions not only provides biological insights for protein function studies but also has great significance for drug discoveries. Therefore, predicting protein-ligand binding sites has long been under intense research in the fields of bioinformatics and computer aided drug discovery. In this review, we first introduce the research background of predicting protein-ligand binding sites and then classify the methods into four categories, namely, 3D structure-based, template similarity-based, traditional machine learning-based and deep learning-based methods. We describe representative algorithms in each category and elaborate on machine learning and deep learning-based prediction methods in more detail. Finally, we discuss the trends and challenges of the current research such as molecular dynamics simulation based cryptic binding sites prediction, and highlight prospective directions for the near future.