webTWAS: a resource for disease candidate susceptibility genes identified by transcriptome-wide association study

TWAS Atlas: a curated knowledgebase of transcriptome-wide association studies

Transcriptome-wide association studies (TWASs), as a practical and prevalent approach for detecting the associations between genetically regulated genes and traits, are now leading to a better understanding of the complex mechanisms of genetic variants in regulating various diseases and traits. Despite the ever-increasing TWAS outputs, there is still a lack of databases curating massive public TWAS information and knowledge. To fill this gap, here we present TWAS Atlas (https://ngdc.cncb.ac.cn/twas/), an integrated knowledgebase of TWAS findings manually curated from extensive literature. In the current implementation, TWAS Atlas collects 401,266 high-quality human gene-trait associations from 200 publications, covering 22,247 genes and 257 traits across 135 tissue types. In particular, an interactive knowledge graph of the collected gene-trait associations is constructed together with single nucleotide polymorphism (SNP)-gene associations to build up comprehensive regulatory networks at multi-omics levels. In addition, TWAS Atlas, as a user-friendly web interface, efficiently enables users to browse, search and download all association information, relevant research metadata and annotation information of interest. Taken together, TWAS Atlas is of great value for promoting the utility and availability of TWAS results in explaining the complex genetic basis as well as providing new insights for human health and disease research.

Machine learning based identification of hub genes in renal clear cell carcinoma using multi-omics data

Renal cell carcinoma is one of the most universal urinary system cancers in the world. The most common renal cell carcinoma subtype is renal clear cell carcinoma. It is usually associated with high rates of metastasis and mortality. Therefore, finding effective therapeutic targets and prognostic molecular markers is of great significance to improve the early diagnosis rate and prognostic accuracy of renal clear cell carcinoma. In this work, we successfully identified six hub genes that are closely related to the occurrence, development and prognosis of renal clear cell carcinoma and proposed three new potential prognostic markers, namely ATP4B, AC144831.1 and Tfcp2l1 through differentially expressed genes (DEGs) analysis, GO functional enrichment and KEGG pathway analysis, WGCNA analysis, and survival analysis. In addition, we established machine learning models to predict the occurrence of tumors through the gene expression data of patients. It is expected that the results of this study can provide reference value for the treatment of renal clear cell carcinoma.

Dynamic transcriptome analysis suggests the key genes regulating seed development and filling in Tartary buckwheat (Fagopyrum tataricum Garetn.)

Tartary buckwheat is highly attractive for the richness of nutrients and quality, yet post-embryonic seed abortion greatly halts the yield. Seed development is crucial for determining grain yield, whereas the molecular basis and regulatory network of Tartary buckwheat seed development and filling is not well understood at present. Here, we assessed the transcriptional dynamics of filling stage Tartary buckwheat seeds at three developmental stages by RNA sequencing. Among the 4249 differentially expressed genes (DEGs), genes related to seed development were identified. Specifically, 88 phytohormone biosynthesis signaling genes, 309 TFs, and 16 expansin genes participating in cell enlargement, 37 structural genes involved in starch biosynthesis represented significant variation and were candidate key seed development genes. Cis-element enrichment analysis indicated that the promoters of differentially expressed expansin genes and starch biosynthesis genes are rich of hormone-responsive (ABA-, AUX-, ET-, and JA-), and seed growth-related (MYB, MYC and WRKY) binding sites. The expansin DEGs showed strong correlations with DEGs in phytohormone pathways and transcription factors (TFs). In total, phytohormone ABA, AUX, ET, BR and CTK, and related TFs could substantially regulate seed development in Tartary buckwheat through targeting downstream expansin genes and structural starch biosynthetic genes. This transcriptome data could provide a theoretical basis for improving yield of Tartary buckwheat.

Identify Bitter Peptides by Using Deep Representation Learning Features

A bitter taste often identifies hazardous compounds and it is generally avoided by most animals and humans. Bitterness of hydrolyzed proteins is caused by the presence of bitter peptides. To improve palatability, bitter peptides need to be identified experimentally in a time-consuming and expensive process, before they can be removed or degraded. Here, we report the development of a machine learning prediction method, iBitter-DRLF, which is based on a deep learning pre-trained neural network feature extraction method. It uses three sequence embedding techniques, soft symmetric alignment (SSA), unified representation (UniRep), and bidirectional long short-term memory (BiLSTM). These were initially combined into various machine learning algorithms to build several models. After optimization, the combined features of UniRep and BiLSTM were finally selected, and the model was built in combination with a light gradient boosting machine (LGBM). The results showed that the use of deep representation learning greatly improves the ability of the model to identify bitter peptides, achieving accurate prediction based on peptide sequence data alone. By helping to identify bitter peptides, iBitter-DRLF can help research into improving the palatability of peptide therapeutics and dietary supplements in the future. A webserver is available, too.

Heterozygosity and homozygosity regions affect reproductive success and the loss of reproduction: a case study with litter traits in pigs

Runs of heterozygosity (ROHet) and homozygosity (ROH) harbor useful information related to traits of interest. There is a lack of investigating the effect of ROHet and ROH on reproductive success and the loss of reproduction in mammals. Here, we detected and characterized the ROHet and ROH patterns in the genomes of Chinese indigenous pigs (i.e., Jinhua, Chun’an, Longyou Black, and Shengxian Spotted pigs), revealing the similar genetic characteristics of indigenous pigs. Later, we highlighted the underlying litter traits-related ROHet and ROH using association analysis with linear model in these four indigenous pig breeds. To pinpoint the promising candidate genes associated with litter traits, we further in-depth explore the selection patterns of other five pig breeds (i.e., Erhualian, Meishan, Minzhu, Rongchang, and Diqing pigs) with different levels of reproduction performance at the underlying litter traits-related ROHet and ROH using F_ST and genetic diversity ratio. Then, we identified a set of known and novel candidate genes associated with reproductive performance in pigs. For the novel candidate genes (i.e., CCDC91, SASH1, SAMD5, MACF1, MFSD2A, EPC2, and MBD5), we obtained public available datasets and performed multi-omics analyses integrating transcriptome-wide association studies and comparative single-cell RNA-seq analyses to uncover the roles of them in mammalian reproductive performance. The genes have not been widely reported to be fertility-related genes and can be complementally considered as prior biological information to modify genomic selections models that benefits pig genetic improvement of litter traits. Besides, our findings provide new insights into the function of ROHet and ROH in mammals.

A survey on protein–DNA-binding sites in computational biology

Transcription factors are important cellular components of the process of gene expression control. Transcription factor binding sites are locations where transcription factors specifically recognize DNA sequences, targeting gene-specific regions and recruiting transcription factors or chromatin regulators to fine-tune spatiotemporal gene regulation. As the common proteins, transcription factors play a meaningful role in life-related activities. In the face of the increase in the protein sequence, it is urgent how to predict the structure and function of the protein effectively. At present, protein–DNA-binding site prediction methods are based on traditional machine learning algorithms and deep learning algorithms. In the early stage, we usually used the development method based on traditional machine learning algorithm to predict protein–DNA-binding sites. In recent years, methods based on deep learning to predict protein–DNA-binding sites from sequence data have achieved remarkable success. Various statistical and machine learning methods used to predict the function of DNA-binding proteins have been proposed and continuously improved. Existing deep learning methods for predicting protein–DNA-binding sites can be roughly divided into three categories: convolutional neural network (CNN), recursive neural network (RNN) and hybrid neural network based on CNN–RNN. The purpose of this review is to provide an overview of the computational and experimental methods applied in the field of protein–DNA-binding site prediction today. This paper introduces the methods of traditional machine learning and deep learning in protein–DNA-binding site prediction from the aspects of data processing characteristics of existing learning frameworks and differences between basic learning model frameworks. Our existing methods are relatively simple compared with natural language processing, computational vision, computer graphics and other fields. Therefore, the summary of existing protein–DNA-binding site prediction methods will help researchers better understand this field.

Mineralocorticoid receptor and glucocorticoid receptor work alone and together in cell-type-specific manner: Implications for resilience prediction and targeted therapy

‘You can't roll the clock back and reverse the effects of experiences' Bruce McEwen used to say when explaining how allostasis labels the adaptive process. Here we will for once roll the clock back to the times that the science of the glucocorticoid hormone was honored with a Nobel prize and highlight the discovery of their receptors in the hippocampus as inroad to its current status as master regulator in control of stress coping and adaptation. Glucocorticoids operate in concert with numerous neurotransmitters, neuropeptides, and other hormones with the aim to facilitate processing of information in the neurocircuitry of stress, from anticipation and perception of a novel experience to behavioral adaptation and memory storage. This action, exerted by the glucocorticoids, is guided by two complementary receptor systems, mineralocorticoid receptors (MR) and glucocorticoid receptors (GR), that need to be balanced for a healthy stress response pattern. Here we discuss the cellular, neuroendocrine, and behavioral studies underlying the MR:GR balance concept, highlight the relevance of hypothalamic-pituitary-adrenal (HPA) -axis patterns and note the limited understanding yet of sexual dimorphism in glucocorticoid actions. We conclude with the prospect that (i) genetically and epigenetically regulated receptor variants dictate cell-type-specific transcriptome signatures of stress-related neuropsychiatric symptoms and (ii) selective receptor modulators are becoming available for more targeted treatment. These two new developments may help to ‘restart the clock’ with the prospect to support resilience.

Identification of Diagnostic Markers for Breast Cancer Based on Differential Gene Expression and Pathway Network

Background: Breast cancer is the second largest cancer in the world, the incidence of breast cancer continues to rise worldwide, and women's health is seriously threatened. Therefore, it is very important to explore the characteristic changes of breast cancer from the gene level, including the screening of differentially expressed genes and the identification of diagnostic markers. Methods: The gene expression profiles of breast cancer were obtained from the TCGA database. The edgeR R software package was used to screen the differentially expressed genes between breast cancer patients and normal samples. The function and pathway enrichment analysis of these genes revealed significant enrichment of functions and pathways. Next, download these pathways from KEGG website, extract the gene interaction relations, construct the KEGG pathway gene interaction network. The potential diagnostic markers of breast cancer were obtained by combining the differentially expressed genes with the key genes in the network. Finally, these markers were used to construct the diagnostic prediction model of breast cancer, and the predictive ability of the model and the diagnostic ability of the markers were verified by internal and external data. Results: 1060 differentially expressed genes were identified between breast cancer patients and normal controls. Enrichment analysis revealed 28 significantly enriched pathways (p < 0.05). They were downloaded from KEGG website, and the gene interaction relations were extracted to construct the gene interaction network of KEGG pathway, which contained 1277 nodes and 7345 edges. The key nodes with a degree greater than 30 were extracted from the network, containing 154 genes. These 154 key genes shared 23 genes with differentially expressed genes, which serve as potential diagnostic markers for breast cancer. The 23 genes were used as features to construct the SVM classification model, and the model had good predictive ability in both the training dataset and the validation dataset (AUC = 0.960 and 0.907, respectively). Conclusion: This study showed that the difference of gene expression level is important for the diagnosis of breast cancer, and identified 23 breast cancer diagnostic markers, which provides valuable information for clinical diagnosis and basic treatment experiments.

Bioinformatics Analysis of Long Non-coding RNA and Related Diseases: An Overview

Long non-coding RNAs (lncRNAs) are usually located in the nucleus and cytoplasm of cells. The transcripts of lncRNAs are >200 nucleotides in length and do not encode proteins. Compared with small RNAs, lncRNAs have longer sequences, more complex spatial structures, and more diverse and complex mechanisms involved in the regulation of gene expression. LncRNAs are widely involved in the biological processes of cells, and in the occurrence and development of many human diseases. Many studies have shown that lncRNAs can induce the occurrence of diseases, and some lncRNAs undergo specific changes in tumor cells. Research into the roles of lncRNAs has covered the diagnosis of, for example, cardiovascular, cerebrovascular, and central nervous system diseases. The bioinformatics of lncRNAs has gradually become a research hotspot and has led to the discovery of a large number of lncRNAs and associated biological functions, and lncRNA databases and recognition models have been developed. In this review, the research progress of lncRNAs is discussed, and lncRNA-related databases and the mechanisms and modes of action of lncRNAs are described. In addition, disease-related lncRNA methods and the relationships between lncRNAs and human lung adenocarcinoma, rectal cancer, colon cancer, heart disease, and diabetes are discussed. Finally, the significance and existing problems of lncRNA research are considered.