Myocardial infarction (MI) is a leading cause of death in the world and many genes are involved in it. Transcription factor (TFs) and microRNAs (miRNAs) are key regulators of gene expression. We hypothesized that miRNAs and TFs might play combinatory regulatory roles in MI. After collecting MI candidate genes and miRNAs from various resources, we constructed a comprehensive MI-specific miRNA-TF co-regulatory network by integrating predicted and experimentally validated TF and miRNA targets. We found some hub nodes (e.g. miR-16 and miR-26) in this network are important regulators, and the network can be severed as a bridge to interpret the associations of previous results, which is shown by the case of miR-29 in this study. We also constructed a regulatory network for MI recurrence and found several important genes (e.g. DAB2, BMP6, miR-320 and miR-103), the abnormal expressions of which may be potential regulatory mechanisms and markers of MI recurrence. At last we proposed a cellular model to discuss major TF and miRNA regulators with signaling pathways in MI. This study provides more details on gene expression regulation and regulators involved in MI progression and recurrence. It also linked up and interpreted many previous results.

The plastids and mitochondria of the eukaryotic cell are of endosymbiotic origin. These events occurred ~2 billion years ago and produced significant changes in the genomes of the host and the endosymbiont. Previous studies demonstrated that the invasion of land affected plastids and mitochondria differently and that the paths of mitochondrial integration differed between animals and plants. Other studies examined the reasons why a set of proteins remained encoded in the organelles and were not transferred to the nuclear genome. However, our understanding of the functional relations of the transferred genes is insufficient. In this paper, we report a high-throughput phylogenetic analysis to identify genes of cyanobacterial origin for plants of different levels of complexity: Arabidopsis thaliana, Chlamydomonas reinhardtii, Physcomitrella patens, Populus trichocarpa, Selaginella moellendorffii, Sorghum bicolor, Oryza sativa, and Ostreococcus tauri. Thus, a census of cyanobacterial gene recruits and a study of their function are presented to better understand the functional aspects of plastid symbiogenesis. From algae to angiosperms, the GO terms demonstrated a gradual expansion over functionally related genes in the nuclear genome, beginning with genes related to thylakoids and photosynthesis, followed by genes involved in metabolism, and finally with regulation-related genes, primarily in angiosperms. The results demonstrate that DNA is supplied to the nuclear genome on a permanent basis with no regard to function, and only what is needed is kept, which thereby expands on the GO space along the related genes.

Long non-coding RNAs (lncRNAs) play key roles in various cellular contexts and diseases by diverse mechanisms. With the rapid growth of identified lncRNAs and disease-associated single nucleotide polymorphisms (SNPs), there is a great demand to study SNPs in lncRNAs. Aiming to provide a useful resource about lncRNA SNPs, we systematically identified SNPs in lncRNAs and analyzed their potential impacts on lncRNA structure and function. In total, we identified 495,729 and 777,095 SNPs in more than 30,000 lncRNA transcripts in human and mouse, respectively. A large number of SNPs were predicted with the potential to impact on the miRNA-lncRNA interaction. The experimental evidence and conservation of miRNA-lncRNA interaction, as well as miRNA expressions from TCGA were also integrated to prioritize the miRNA-lncRNA interactions and SNPs on the binding sites. Furthermore, by mapping SNPs to GWAS results, we found that 142 human lncRNA SNPs are GWAS tagSNPs and 197,827 lncRNA SNPs are in the GWAS linkage disequilibrium regions. All these data for human and mouse lncRNAs were imported into lncRNASNP database (http://bioinfo.life.hust.edu.cn/lncRNASNP/), which includes two sub-databases lncRNASNP-human and lncRNASNP-mouse. The lncRNASNP database has a user-friendly interface for searching and browsing through the SNP, lncRNA and miRNA sections.

Transcription factors (TFs) are proteins that bind to specific DNA sequences, thereby playing crucial roles in gene-expression regulation through controlling the transcription of genetic information from DNA to RNA. Transcription cofactors and chromatin remodeling factors are also essential in the gene transcriptional regulation. Identifying and annotating all the TFs are primary and crucial steps for illustrating their functions and understanding the transcriptional regulation. In this study, based on manual literature reviews, we collected and curated 72 TF families for animals, which is currently the most complete list of TF families in animals. Then, we systematically characterized all the TFs in 50 animal species and constructed a comprehensive animal TF database, AnimalTFDB. To better serve the community, we provided detailed annotations for each TF, including basic information, gene structure, functional domain, 3D structure hit, Gene Ontology, pathway, protein-protein interaction, paralogs, orthologs, potential TF-binding sites and targets. In addition, we collected and annotated transcription cofactors and chromatin remodeling factors. AnimalTFDB has a user-friendly web interface with multiple browse and search functions, as well as data downloading. It is freely available at http://www.bioguo.org/AnimalTFDB/.

The Macleaya spp., including Macleaya cordata and Macleaya microcarpa, are traditional anti-virus, inflammation eliminating, and insecticide herb medicines for their isoquinoline alkaloids. They are also known as the basis of the popular natural animal food addictive in Europe. However, few studies especially at genomics level were conducted on them. Hence, we performed the Macleaya spp. transcriptome and integrated it with iTRAQ proteome analysis in order to identify potential genes involved in alkaloids biosynthesis.

Cancer genes tend to be highly mutated under positive selection. Better understanding the recurrently mutated genes (RMGs) in cancer is critical for explicating the mechanisms of tumorigenesis and providing vital clues for therapy. Although some studies have investigated functional impacts of RMGs in specific cancer types, a comprehensive analysis of RMGs and their mutational impacts across cancers is still needed.

The protein translational system, including transfer RNAs (tRNAs) and several categories of enzymes, plays a key role in regulating cell proliferation. Translation dysregulation also contributes to cancer development, though relatively little is known about the changes that occur to the translational system in cancer. Here, we present global analyses of tRNAs and three categories of enzymes involved in translational regulation in ~10,000 cancer patients across 31 cancer types from The Cancer Genome Atlas. By analyzing the expression levels of tRNAs at the gene, codon, and amino acid levels, we identified unequal alterations in tRNA expression, likely due to the uneven distribution of tRNAs decoding different codons. We find that overexpression of tRNAs recognizing codons with a low observed-over-expected ratio may overcome the translational bottleneck in tumorigenesis. We further observed overall overexpression and amplification of tRNA modification enzymes, aminoacyl-tRNA synthetases, and translation factors, which may play synergistic roles with overexpression of tRNAs to activate the translational systems across multiple cancer types.

Renal impairment (RI) is one of the hallmarks of multiple myeloma (MM) and carries a poor prognosis. Microvesicles (MVs) are membrane vesicles and play an important role in disease progression. Here, we investigated the role of MVs derived from MM cells (MM-MVs) in RI of MM. We found that MM-MVs significantly inhibited viability and induced apoptosis, but not epithelial-mesenchymal transition in human kidney-2 (HK-2), a human renal tubular epithelial cell line. The protein levels of cleaved caspase-3, 8, and 9, and E-cadherin, were increased, but vementin levels were decreased in the HK-2 cells treated with MM-MVs. Through a comparative sequencing and analysis of RNA content between the MVs from RPMI8226 MM cells (RPMI8226-MVs) and K562 leukemia cells, RPMI8226-MVs were enriched with more renal-pathogenic miRNAs, in which the selective miRNAs may participate in the up-regulation of the levels of cleaved caspase-3. Furthermore, the levels of CD138+ circulating MVs (cirMVs) in the peripheral blood were positively correlated with the severity of RI in newly-diagnosed MM. Our study supports MM-MVs representing a previously undescribed factor and playing a potential role in the development of RI of MM patients, and sheds light on the potential application of CD138+ cirMV counts in precise diagnosis of RI in MM and exploring MM-MVs as a therapeutic target.

Half of the chronic myeloid leukemia (CML) patients with sustained deep molecular response suffer from relapse after discontinuation mainly because tyrosine kinase inhibitors (TKIs) cannot eradicate leukemia stem cells (LSCs). In addition, tumor necrosis factor α (TNF-α) is highly detected in CML patients. Our aim was to explore whether TNF-α is a potential target for LSC elimination.

T cells and the T-cell receptor (TCR) repertoire play pivotal roles in immune response and immunotherapy. TCR sequencing (TCR-Seq) technology has enabled accurate profiling TCR repertoire and currently a large number of TCR-Seq data are available in public. Based on the urgent need to effectively re-use these data, we developed TCRdb, a comprehensive human TCR sequences database, by a uniform pipeline to characterize TCR sequences on TCR-Seq data. TCRdb contains more than 277 million highly reliable TCR sequences from over 8265 TCR-Seq samples across hundreds of tissues/clinical conditions/cell types. The unique features of TCRdb include: (i) comprehensive and reliable sequences for TCR repertoire in different samples generated by a strict and uniform pipeline of TCRdb; (ii) powerful search function, allowing users to identify their interested TCR sequences in different conditions; (iii) categorized sample metadata, enabling comparison of TCRs in different sample types; (iv) interactive data visualization charts, describing the TCR repertoire in TCR diversity, length distribution and V-J gene utilization. The TCRdb database is freely available at http://bioinfo.life.hust.edu.cn/TCRdb/ and will be a useful resource in the research and application community of T cell immunology.

DNA methylation plays essential roles in tumor occurrence and stemness maintenance. Tumor-repopulating cells (TRCs) are cancer stem cell (CSC)-like cells with highly tumorigenic and self-renewing abilities, which were selected from tumor cells in soft three-dimensional (3D) fibrin gels.

Although CD19 chimeric antigen receptor-T (CAR-T) cells therapy has achieved unparalleled success in B cell malignancies. The dysfunction of CAR-T cells due to exhaustion is considered as a key factor for treatment failure, and the mechanisms of exhaustion remain elusive. Extracellular vesicles (EVs), important media for communication between tumor and immune cells, may contribute to CAR-T cell exhaustion. Here, we demonstrated that CD19+ tumor cells derived EVs (NALM6-EVs) can carry CD19 antigen and activate CD19 CAR-T cells. The transient activation induced a supraphysiologic inflammatory state with increased release of multiple cytokines. Besides, the sustained activation led CD19 CAR-T cells to enter an exhausted state with upregulated inhibitory receptors, decreased expansion ability, exaggerated effector cell differentiation and impaired antitumor activity. Transcriptomic profiling validated these findings and identified dynamic changes in CD8+ effector T, CD8+ exhausted T, CD8+RRM2+ T and T helper cell subpopulations during activation to exhaustion, as well as changes in many cytokines, inflammatory and immune-related pathways. Our findings identify a credible mechanism of CAR-T cell exhaustion that driven by tumor-derived EVs and provide a novel possible trigger for early cytokine release syndrome.

The success of immune checkpoint blockade (ICB) promotes the immunotherapy to be a new pillar in cancer treatment. However, the low response rate of the ICB therapy limits its application. To increase the response rate and enhance efficacy, the ICB combination therapy has emerged and its clinical trials are increasing. Nevertheless, the gene expression profile and its pattern of ICB combination were not comprehensively studied, which limits the understanding of the ICB combination therapy and the identification of new drugs. Here, we constructed ICBcomb (http://bioinfo.life.hust.edu.cn/ICBcomb/), a comprehensive database, by analyzing the human and mouse expression data of the ICB combination therapy and comparing them between groups treated with ICB, other drugs or their combinations. ICBcomb contains 1399 samples across 29 cancer types involving 52 drugs. It provides a user-friendly web interface for demonstrating the results of the available comparisons in the ICB combination therapy datasets with five functional modules: [1, 2] the 'Dataset/Disease' modules for browsing the expression, enrichment and comparison results in each dataset or disease; [3] the 'Gene' module for inputting a gene symbol and displaying its expression and comparison results across datasets/diseases; [4] the 'Gene Set' module for GSVA/GSEA enrichment analysis on the built-in gene sets and the user-input gene sets in different comparisons; [5] the 'Immune Cell' module for immune cell infiltration comparison between different groups by immune cell abundance analysis. The ICBcomb database provides the first resource for gene expression profile and comparison in ICB combination therapy, which may provide clues for discovering the mechanism of effective combination strategies and new combinatory drugs.

Formalin-fixed paraffin-embedded (FFPE) tissues are widely available specimens for clinical studies. However, RNA degradation in FFPE tissues often restricts their utility. In this study, we determined optimal FFPE preparation conditions, including tissue ischemia at 4°C (<48 h) or 25°C for a short time (0.5 h), 48-h fixation at 25°C and sampling from FFPE scrolls instead of sections. Notably, we observed an increase in intronic reads and a significant change in gene rank based on expression level in the FFPE as opposed to fresh-frozen (FF) samples. Additionally, we found that more reads were mapped to genes associated with chemical stimulus in FFPE samples. Furthermore, we demonstrated that more degraded genes in FFPE samples were enriched in genes with short transcripts and high free energy. Besides, we found 40 housekeeping genes exhibited stable expression in FF and FFPE samples across various tissues. Moreover, our study showed that FFPE samples yielded comparable results to FF samples in dimensionality reduction and pathway analyses between case and control samples. Our study established the optimal conditions for FFPE preparation and identified gene attributes associated with degradation, which would provide useful clues for the utility of FFPE tissues in clinical practice and research.

Type 2 diabetes (T2D) is generally regarded as a metabolic disorder disease with various phenotypic expressions. Traditional Chinese medicine (TCM) has been widely used for preventing and treating diabetes. In our study, we demonstrated that Cyclocarya paliurus formula extractum (CPE), a compound of TCM, can ameliorate diabetes in diabetic rats. Transcriptome profiles were performed to elucidate the anti-diabetic mechanisms of CPE on pancreas and liver. Pancreatic transcriptome analysis showed CPE treatment significantly inhibited gene expressions related to inflammation and apoptosis pathways, among which the transcription factors (TFs) nuclear factor κB (NF-κB), STAT, and miR-9a/148/200 may serve as core regulators contributing to ameliorate diabetes. Biochemical studies also demonstrated CPE treatment decreased pro-inflammatory cytokines (tumor necrosis factor alpha [TNF-α], interleukin [IL]-1β, and IL-6) and reduced β cell apoptosis. In liver tissue, our transcriptome and biochemical experiments showed that CPE treatment reduced lipid accumulation and liver injury, and it promoted glycogen synthesis, which may be regulated by TFs Srebf1, Mlxipl, and miR-122/128/192. Taken together, our findings revealed CPE could be used as a potential therapeutic agent to prevent and treat diabetes. It is the first time to combine transcriptome and regulatory network analyses to study the mechanism of CPE in preventing diabetes, giving a demonstration of exploring the mechanism of TCM on complex diseases.

Mutations in EYS are associated with autosomal recessive retinitis pigmentosa (arRP) and autosomal recessive cone-rod dystrophy (arCRD) however, the function of EYS and the molecular mechanisms of how these mutations cause retinal degeneration are still unclear. Because EYS is absent in mouse and rat, and the structure of the retina differs substantially between humans and Drosophila, we utilised zebrafish as a model organism to study the function of EYS in the retina. We constructed an EYS-knockout zebrafish-line by TALEN technology which showed visual impairment at an early age, while the histological and immunofluorescence assays indicated the presence of progressive retinal degeneration with a cone predominately affected pattern. These phenotypes recapitulate the clinical manifestations of arCRD patients. Furthermore, the EYS-/- zebrafish also showed mislocalisation of certain outer segment proteins (rhodopsin, opn1lw, opn1sw1, GNB3 and PRPH2), and disruption of actin filaments in photoreceptors. Protein mislocalisation may, therefore, disrupt the function of cones and rods in these zebrafish and cause photoreceptor death. Collectively, these results point to a novel role for EYS in maintaining the morphological structure of F-actin and in protein transport, loss of this function might be the trigger for the resultant cellular events that ultimately lead to photoreceptor death.

The pediatric T cell acute lymphoblastic leukemia (T-ALL) still remains a cancer with worst prognosis for high recurrence. Massive studies were conducted for the leukemia relapse based on diagnosis and relapse paired samples. However, the initially diagnostic samples may contain the relapse information and mechanism, which were rarely studied. In this study, we collected mRNA and microRNA (miRNA) data from initially diagnosed pediatric T-ALL samples with their relapse or remission status after treatment. Integrated differential co-expression and miRNA-transcription factor (TF)-gene regulatory network analyses were used to reveal the possible relapse mechanisms for pediatric T-ALL. We detected miR-1246/1248 and NOTCH2 served as key nodes in the relapse network, and they combined with TF WT1/SOX4/REL to form regulatory modules that influence the progress of T-ALL. A regulatory loop miR-429-MYCN-MFHAS1 was found potentially associated with the remission of T-ALL. Furthermore, we proved miR-1246/1248 combined with NOTCH2 could promote cell proliferation in the T-ALL cell line by experiments. Meanwhile, analysis based on the miRNA-drug relationships demonstrated that drugs 5-fluorouracil, ascorbate, and trastuzumab targeting miR-1246 could serve as potential supplements for the standard therapy. In conclusion, our findings revealed the potential molecular mechanisms of T-ALL relapse by the combination of co-expression and regulatory network, and they provide preliminary clues for precise treatment of T-ALL patients.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy. The understanding of its gene expression regulation and molecular mechanisms still remains elusive. Started from experimentally verified T-ALL-related miRNAs and genes, we obtained 120 feed-forward loops (FFLs) among T-ALL-related genes, miRNAs and TFs through combining target prediction. Afterwards, a T-ALL miRNA and TF co-regulatory network was constructed, and its significance was tested by statistical methods. Four miRNAs in the miR-17-92 cluster and four important genes (CYLD, HOXA9, BCL2L11 and RUNX1) were found as hubs in the network. Particularly, we found that miR-19 was highly expressed in T-ALL patients and cell lines. Ectopic expression of miR-19 represses CYLD expression, while miR-19 inhibitor treatment induces CYLD protein expression and decreases NF-κB expression in the downstream signaling pathway. Thus, miR-19, CYLD and NF-κB form a regulatory FFL, which provides new clues for sustained activation of NF-κB in T-ALL. Taken together, we provided the first miRNA-TF co-regulatory network in T-ALL and proposed a model to demonstrate the roles of miR-19 and CYLD in the T-cell leukemogenesis. This study may provide potential therapeutic targets for T-ALL and shed light on combining bioinformatics with experiments in the research of complex diseases.

Cytogenetically normal acute myeloid leukemia (CN-AML) presents with diverse outcomes in different patients and is categorized as an intermediate prognosis group. It is important to identify prognostic factors for CN-AML risk stratification. In this study, using the TCGA CN-AML dataset, we found that the scavenger receptor stabilin-1 (STAB1) is a prognostic factor for poor outcomes and validated it in three other independent CN-AML datasets. The high STAB1 expression (STAB1high) group had shorter event-free survival compared with the low STAB1 expression (STAB1low) group in both the TCGA dataset (n = 79; p = 0.0478) and GEO: GSE6891 dataset (n = 187; p = 0.0354). Differential expression analysis between the STAB1high and STAB1low groups revealed that upregulated genes in the STAB1high group were enriched in pathways related to cell adhesion and migration and immune responses. We confirmed that STAB1 suppression inhibits cell growth in KG1a and NB4 leukemia cells. Expression correlation analyses between STAB1 and cancer drug targets suggested that patients in the STAB1low group are more sensitive to the BCL2 inhibitor venetoclax, and we confirmed it in cell lines. In conclusion, we identified STAB1 as a prognostic factor for CN-AML in multiple datasets, explored its underlying mechanism, and provided potential therapeutic indications.

DNA methylation is an important epigenetic mechanism for regulating gene expression. Aberrant DNA methylation has been observed in various human diseases, including cancer. Single-nucleotide polymorphisms can contribute to tumor initiation, progression and prognosis by influencing DNA methylation, and DNA methylation quantitative trait loci (meQTL) have been identified in physiological and pathological contexts. However, no database has been developed to systematically analyze meQTLs across multiple cancer types. Here, we present Pancan-meQTL, a database to comprehensively provide meQTLs across 23 cancer types from The Cancer Genome Atlas by integrating genome-wide genotype and DNA methylation data. In total, we identified 8 028 964 cis-meQTLs and 965 050 trans-meQTLs. Among these, 23 432 meQTLs are associated with patient overall survival times. Furthermore, we identified 2 214 458 meQTLs that overlap with known loci identified through genome-wide association studies. Pancan-meQTL provides a user-friendly web interface (http://bioinfo.life.hust.edu.cn/Pancan-meQTL/) that is convenient for browsing, searching and downloading data of interest. This database is a valuable resource for investigating the roles of genetics and epigenetics in cancer.