GSCA: an integrated platform for gene set cancer analysis at genomic, pharmacogenomic and immunogenomic levels.

Cancer initiation and progression are likely caused by the dysregulation of biological pathways. Gene set analysis (GSA) could improve the signal-to-noise ratio and identify potential biological insights on the gene set level. However, platforms exploring cancer multi-omics data using GSA methods are lacking. In this study, we upgraded our GSCALite to GSCA (gene set cancer analysis, http://bioinfo.life.hust.edu.cn/GSCA) for cancer GSA at genomic, pharmacogenomic and immunogenomic levels. In this improved GSCA, we integrated expression, mutation, drug sensitivity and clinical data from four public data sources for 33 cancer types. We introduced useful features to GSCA, including associations between immune infiltration with gene expression and genomic variations, and associations between gene set expression/mutation and clinical outcomes. GSCA has four main functional modules for cancer GSA to explore, analyze and visualize expression, genomic variations, tumor immune infiltration, drug sensitivity and their associations with clinical outcomes. We used case studies of three gene sets: (i) seven cell cycle genes, (ii) tumor suppressor genes of PI3K pathway and (iii) oncogenes of PI3K pathway to prove the advantage of GSCA over single gene analysis. We found novel associations of gene set expression and mutation with clinical outcomes in different cancer types on gene set level, while on single gene analysis level, they are not significant associations. In conclusion, GSCA is a user-friendly web server and a useful resource for conducting hypothesis tests by using GSA methods at genomic, pharmacogenomic and immunogenomic levels.

Population structure and selection signal analysis of indigenous sheep from the southern edge of the Taklamakan Desert.

Analyzing the genetic diversity and selection characteristics of sheep (Ovis aries) holds significant value in understanding their environmental adaptability, enhancing breeding efficiency, and achieving effective conservation and rational utilization of genetic resources. In this study, we utilized Illumina Ovine SNP 50 K BeadChip data from four indigenous sheep breeds from the southern margin of the Taklamakan Desert (Duolang sheep: n = 36, Hetian sheep: n = 74, Kunlun sheep: n = 27, Qira black sheep: n = 178) and three foreign meat sheep breeds (Poll Dorset sheep: n = 105, Suffolk sheep: n = 153, Texel sheep: n = 150) to investigate the population structure, genetic diversity, and genomic signals of positive selection within the indigenous sheep. According to the Principal component analysis (PCA), the Neighbor-Joining tree (NJ tree), and Admixture, we revealed distinct clustering patterns of these seven sheep breeds based on their geographical distribution. Then used Cross Population Extended Haplotype Homozygosity (XP-EHH), Fixation Index (FST), and Integrated Haplotype Score (iHS), we identified a collective set of 32 overlapping genes under positive selection across four indigenous sheep breeds. These genes are associated with wool follicle development and wool traits, desert environmental adaptability, disease resistance, reproduction, and high-altitude adaptability. This study reveals the population structure and genomic selection characteristics in the extreme desert environments of native sheep breeds from the southern edge of the Taklimakan Desert, providing new insights into the conservation and sustainable use of indigenous sheep genetic resources in extreme environments. Additionally, these findings offer valuable genetic resources for sheep and other mammals to adapt to global climate change.

EVAtool: an optimized reads assignment tool for small ncRNA quantification and its application in extracellular vesicle datasets.

Extracellular vesicles (EVs) carrying various small non-coding RNAs (sncRNAs) play a vital roles in cell communication and diseases. Correct quantification of multiple sncRNA biotypes simultaneously in EVs is a challenge due to the short reads (<30 bp) could be mapped to multiple sncRNA types. To address this question, we developed an optimized reads assignment algorithm (ORAA) to dynamically map multi-mapping reads to the sncRNA type with a higher proportion. We integrated ORAA with reads processing steps into EVAtool Python-package (http://bioinfo.life.hust.edu.cn/EVAtool) to quantify sncRNAs, especially for sncRNA-seq from EV samples. EVAtool allows users to specify interested sncRNA types in advanced mode or use default seven sncRNAs (microRNA, small nucleolar RNA, PIWI-interacting RNAs, small nuclear RNA, ribosomal RNA, transfer RNA and Y RNA). To prove the utilities of EVAtool, we quantified the sncRNA expression profiles for 200 samples from cognitive decline and multiple sclerosis. We found that more than 20% of short reads on average were mapped to multiple sncRNA biotypes in multiple sclerosis. In cognitive decline, the proportion of Y RNA is significantly higher than other sncRNA types. EVAtool is a flexible and extensible tool that would benefit to mine potential biomarkers and functional molecules in EVs.

EVAtlas: a comprehensive database for ncRNA expression in human extracellular vesicles.

Extracellular vesicles (EVs) packing various molecules play vital roles in intercellular communication. Non-coding RNAs (ncRNAs) are important functional molecules and biomarkers in EVs. A comprehensive investigation of ncRNAs expression in EVs under different conditions is a fundamental step for functional discovery and application of EVs. Here, we curated 2030 small RNA-seq datasets for human EVs (1506 sEV and 524 lEV) in 24 conditions and over 40 diseases. We performed a unified reads dynamic assignment algorithm (RDAA) considering mismatch and multi-mapping reads to quantify the expression profiles of seven ncRNA types (miRNA, snoRNA, piRNA, snRNA, rRNA, tRNA and Y RNA). We constructed EVAtlas (http://bioinfo.life.hust.edu.cn/EVAtlas), a comprehensive database for ncRNA expression in EVs with four functional modules: (i) browse and compare the distribution of ncRNAs in EVs from 24 conditions and eight sources (plasma, serum, saliva, urine, sperm, breast milk, primary cell and cell line); (ii) prioritize candidate ncRNAs in condition related tissues based on their expression; (iii) explore the specifically expressed ncRNAs in EVs from 24 conditions; (iv) investigate ncRNA functions, related drugs, target genes and EVs isolation methods. EVAtlas contains the most comprehensive ncRNA expression in EVs and will be a key resource in this field.

A comprehensive survey for human transcription factors on expression, regulation, interaction, phenotype and cancer survival.

Transcription factors (TFs) act as key regulators in biological processes through controlling gene expression. Here, we conducted a systematic study for all human TFs on the expression, regulation, interaction, mutation, phenotype and cancer survival. We revealed that the average expression levels of TFs in normal tissues were lower than 50% expression of non-TFs, whereas TF expression was increased in cancers. TFs that are specifically expressed in an individual tissue or cancer may be potential marker genes. For instance, TGIF2LX/Y were preferentially expressed in testis and NEUROG1, PRDM14, SRY, ZNF705A and ZNF716 were specifically highly expressed in germ cell tumors. We found different distributions of target genes and TF co-regulations in different TF families. Some small TF families have huge protein interaction pairs, suggesting their central roles in transcriptional regulation. The bZIP family is a small family involving many signaling pathways. Survival analysis indicated that most TFs significantly affect survival of one or more cancers. Some survival-related TFs were also specifically highly expressed in the corresponding cancer types, which may be potential targets for cancer therapy. Finally, we identified 43 TFs whose mutations were closely correlated to survival, suggesting their cancer-driven roles. The systematic analysis of TFs provides useful clues for further investigation of TF regulatory mechanisms and the role of TFs in diseases.

CCLA: an accurate method and web server for cancer cell line authentication using gene expression profiles.

Cancer cell lines (CCLs) as important model systems play critical roles in cancer research. The misidentification and contamination of CCLs are serious problems, leading to unreliable results and waste of resources. Current methods for CCL authentication are mainly based on the CCL-specific genetic polymorphism, whereas no method is available for CCL authentication using gene expression profiles. Here, we developed a novel method and homonymic web server (CCLA, Cancer Cell Line Authentication, http://bioinfo.life.hust.edu.cn/web/CCLA/) to authenticate 1291 human CCLs of 28 tissues using gene expression profiles. CCLA showed an excellent speed advantage and high accuracy for CCL authentication, a top 1 accuracy of 96.58 or 92.15% (top 3 accuracy of 100 or 95.11%) for microarray or RNA-Seq validation data (719 samples, 461 CCLs), respectively. To the best of our knowledge, CCLA is the first approach to authenticate CCLs using gene expression data. Users can freely and conveniently authenticate CCLs using gene expression profiles or NCBI GEO accession on CCLA website.

Expression profile of immune checkpoint genes and their roles in predicting immunotherapy response.

Immune checkpoint genes (ICGs) play critical roles in circumventing self-reactivity and represent a novel target to develop treatments for cancers. However, a comprehensive analysis for the expression profile of ICGs at a pan-cancer level and their correlation with patient response to immune checkpoint blockade (ICB) based therapy is still lacking. In this study, we defined three expression patterns of ICGs using a comprehensive survey of RNA-seq data of tumor and immune cells from the functional annotation of the mammalian genome (FANTOM5) project. The correlation between the expression patterns of ICGs and patients survival and response to ICB therapy was investigated. The expression patterns of ICGs were robust across cancers, and upregulation of ICGs was positively correlated with high lymphocyte infiltration and good prognosis. Furthermore, we built a model (ICGe) to predict the response of patients to ICB therapy using five features of ICG expression. A validation scenario of six independent datasets containing data of 261 patients with CTLA-4 and PD-1 blockade immunotherapies demonstrated that ICGe achieved area under the curves of 0.64-0.82 and showed a robust performance and outperformed other mRNA-based predictors. In conclusion, this work revealed expression patterns of ICGs and underlying correlations between ICGs and response to ICB, which helps to understand the mechanisms of ICGs in ICB signal pathways and other anticancer treatments.

miRNASNP-v3: a comprehensive database for SNPs and disease-related variations in miRNAs and miRNA targets.

MicroRNAs (miRNAs) related single-nucleotide variations (SNVs), including single-nucleotide polymorphisms (SNPs) and disease-related variations (DRVs) in miRNAs and miRNA-target binding sites, can affect miRNA functions and/or biogenesis, thus to impact on phenotypes. miRNASNP is a widely used database for miRNA-related SNPs and their effects. Here, we updated it to miRNASNP-v3 (http://bioinfo.life.hust.edu.cn/miRNASNP/) with tremendous number of SNVs and new features, especially the DRVs data. We analyzed the effects of 7 161 741 SNPs and 505 417 DRVs on 1897 pre-miRNAs (2630 mature miRNAs) and 3'UTRs of 18 152 genes. miRNASNP-v3 provides a one-stop resource for miRNA-related SNVs research with the following functions: (i) explore associations between miRNA-related SNPs/DRVs and diseases; (ii) browse the effects of SNPs/DRVs on miRNA-target binding; (iii) functional enrichment analysis of miRNA target gain/loss caused by SNPs/DRVs; (iv) investigate correlations between drug sensitivity and miRNA expression; (v) inquire expression profiles of miRNAs and their targets in cancers; (vi) browse the effects of SNPs/DRVs on pre-miRNA secondary structure changes; and (vii) predict the effects of user-defined variations on miRNA-target binding or pre-miRNA secondary structure. miRNASNP-v3 is a valuable and long-term supported resource in functional variation screening and miRNA function studies.

An ultra-sensitive T-cell receptor detection method for TCR-Seq and RNA-Seq data.

MOTIVATION: T-cell receptors (TCRs) function to recognize antigens and play vital roles in T-cell immunology. Surveying TCR repertoires by characterizing complementarity-determining region 3 (CDR3) is a key issue. Due to the high diversity of CDR3 and technological limitation, accurate characterization of CDR3 repertoires remains a great challenge. RESULTS: We propose a computational method named CATT for ultra-sensitive and precise TCR CDR3 sequences detection. CATT can be applied on TCR sequencing, RNA-Seq and single-cell TCR(RNA)-Seq data to characterize CDR3 repertoires. CATT integrated de Bruijn graph-based micro-assembly algorithm, data-driven error correction model and Bayesian inference algorithm, to self-adaptively and ultra-sensitively characterize CDR3 repertoires with high performance. Benchmark results of datasets from in silico and experimental data demonstrated that CATT showed superior recall and precision compared with existing tools, especially for data with short read length and small size and single-cell sequencing data. Thus, CATT will be a useful tool for TCR analysis in researches of cancer and immunology. AVAILABILITY AND IMPLEMENTATION: http://bioinfo.life.hust.edu.cn/CATT or https://github.com/GuoBioinfoLab/CATT. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SAGD: a comprehensive sex-associated gene database from transcriptomes.

Many animal species present sex differences. Sex-associated genes (SAGs), which have female-biased or male-biased expression, have major influences on the remarkable sex differences in important traits such as growth, reproduction, disease resistance and behaviors. However, the SAGs resulting in the vast majority of phenotypic sex differences are still unknown. To provide a useful resource for the functional study of SAGs, we manually curated public RNA-seq datasets with paired female and male biological replicates from the same condition and systematically re-analyzed the datasets using standardized methods. We identified 27,793 female-biased SAGs and 64,043 male-biased SAGs from 2,828 samples of 21 species, including human, chimpanzee, macaque, mouse, rat, cow, horse, chicken, zebrafish, seven fly species and five worm species. All these data were cataloged into SAGD, a user-friendly database of SAGs (http://bioinfo.life.hust.edu.cn/SAGD) where users can browse SAGs by gene, species, drug and dataset. In SAGD, the expression, annotation, targeting drugs, homologs, ontology and related RNA-seq datasets of SAGs are provided to help researchers to explore their functions and potential applications in agriculture and human health.

tRic: a user-friendly data portal to explore the expression landscape of tRNAs in human cancers.

Transfer RNAs (tRNAs) play critical roles in human cancer. Currently, no database provides the expression landscape and clinical relevance of tRNAs across a variety of human cancers. Utilizing miRNA-seq data from The Cancer Genome Atlas, we quantified the relative expression of tRNA genes and merged them into the codon level and amino level across 31 cancer types. The expression of tRNAs is associated with clinical features of patient smoking history and overall survival, and disease stage, subtype, and grade. We further analysed codon frequency and amino acid frequency for each protein coding gene and linked alterations of tRNA expression with protein translational efficiency. We include these data resources in a user-friendly data portal, tRic (tRNA in cancer, https://hanlab.uth.edu/tRic/ or http://bioinfo.life.hust.edu.cn/tRic/), which can be of significant interest to the research community.

GSCALite: a web server for gene set cancer analysis.

Summary: The availability of cancer genomic data makes it possible to analyze genes related to cancer. Cancer is usually the result of a set of genes and the signal of a single gene could be covered by background noise. Here, we present a web server named Gene Set Cancer Analysis (GSCALite) to analyze a set of genes in cancers with the following functional modules. (i) Differential expression in tumor versus normal, and the survival analysis; (ii) Genomic variations and their survival analysis; (iii) Gene expression associated cancer pathway activity; (iv) miRNA regulatory network for genes; (v) Drug sensitivity for genes; (vi) Normal tissue expression and eQTL for genes. GSCALite is a user-friendly web server for dynamic analysis and visualization of gene set in cancer and drug sensitivity correlation, which will be of broad utilities to cancer researchers. Availability and implementation: GSCALite is available on http://bioinfo.life.hust.edu.cn/web/GSCALite/. Supplementary information: Supplementary data are available at Bioinformatics online.

Unexpected Neuroprotective Effects of Loganin on 1-Methyl-4-Phenyl-1,2,3,6-Tetrahydropyridine-Induced Neurotoxicity and Cell Death in Zebrafish.

1-methyl-4-phenyl-1,2,3,6 tetrahydropyridine (MPTP), which induces the pathological characteristics of Parkinson's disease in rodents, also specifically targets dopaminergic neurons in zebrafish embryos and larvae. Loganin, a traditional Chinese drug, was reported to regulate immune function and possess anti-inflammatory and anti-shock effects. Here, we investigate the role of loganin in MPTP-induced Parkinson-like abnormalities in zebrafish. MPTP treatment-induced abnormal development, in larvae, such as pericardium edema, increased yolk color, yolk sac edema, and retarded yolk sac resorption, as well as defects in brain development. Loganin could block MPTP-induced defects, with little toxicity to the eggs. Results of whole mount in situ hybridization showed loganin prevented the loss of both dopaminergic neurons and locomotor activity, exhibited by larvae treated with MPTP. In addition, loganin significantly rescued MPTP-induced neurotoxicity on PC12 cells, possibly through the suppression of PI3K/Akt/mTOR axis and JNK signaling pathways. In conclusion, loganin blocks MPTP-induced neurotoxicity and abnormal development in zebrafish. J. Cell. Biochem. 118: 615-628, 2017. © 2016 Wiley Periodicals, Inc.

Activating transcription factor 3 promotes spinal cord regeneration of adult zebrafish.

Zebrafish is an excellent model to study the mechanisms underlying successful central nervous system (CNS) regeneration. Previous study shows that activating transcription factor 3 (ATF3) promotes neurite outgrowth and is involved in optic nerve regeneration in zebrafish. Here, we used zebrafish model to investigate the role of ATF3 in regeneration following spinal cord injury (SCI). Quantitative polymerase chain reaction (qPCR) and in situ hybridization revealed that ATF3 mRNA levels increased at 12 h and 6 d following SCI. Double labeled immunofluorescence showed that ATF3 expressed in motoneurons. Treatment of anti-sense ATF3 morpholino (MO) inhibited locomotor recovery and decreased axon regeneration of spinal cord injured zebrafish. Further, inhibition of ATF3 up-regulated the expression of inflammatory factors tumor necrosis factor-α (TNF-α) and interleukin-1ß (IL-1ß). These data suggest that ATF3 could promote locomotor recovery and axon regrowth in zebrafish SCI model possibly by regulating inflammatory response.

AnimalTFDB 2.0: a resource for expression, prediction and functional study of animal transcription factors.

Transcription factors (TFs) are key regulators for gene expression. Here we updated the animal TF database AnimalTFDB to version 2.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB/). Using the improved prediction pipeline, we identified 72 336 TF genes, 21 053 transcription co-factor genes and 6502 chromatin remodeling factor genes from 65 species covering main animal lineages. Besides the abundant annotations (basic information, gene model, protein functional domain, gene ontology, pathway, protein interaction, ortholog and paralog, etc.) in the previous version, we made several new features and functions in the updated version. These new features are: (i) gene expression from RNA-Seq for nine model species, (ii) gene phenotype information, (iii) multiple sequence alignment of TF DNA-binding domains, and the weblogo and phylogenetic tree based on the alignment, (iv) a TF prediction server to identify new TFs from input sequences and (v) a BLAST server to search against TFs in AnimalTFDB. A new nice web interface was designed for AnimalTFDB 2.0 allowing users to browse and search all data in the database. We aim to maintain the AnimalTFDB as a solid resource for TF identification and studies of transcription regulation and comparative genomics.

Tumor Cell-Derived Microvesicles Induced Not Epithelial-Mesenchymal Transition but Apoptosis in Human Proximal Tubular (HK-2) Cells: Implications for Renal Impairment in Multiple Myeloma.

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.

Beneficial roles of melanoma cell adhesion molecule in spinal cord transection recovery in adult zebrafish.

Melanoma cell adhesion molecule (MCAM) is a multifunctional protein involved in miscellaneous processes, including development and tumor angiogenesis. Here, spinal cord transection in adult zebrafish was used to investigate the effects of MCAM on spinal cord injury (SCI) and subsequent recovery. Expression of MCAM mRNA increased and co-localized with motoneurons in the spinal cord after SCI. With MCAM morpholino treatment, inhibition of MCAM retarded both axon regrowth and locomotor recovery in the spinal cord injured zebrafish. Furthermore, MCAM mRNA expression was also observed in fli1a:EGFP transgenic zebrafish, which specifically show labeled blood vessels. Inhibition of MCAM down-regulated the expression of angiogenesis-related factors, such as VEGFR-2, p-p38 and p-AKT, and the inflammatory factors TNF-α, IL-1ß and IL-8. Taken together, these data suggest that MCAM may have a beneficial role in the recovery from SCI, via the promotion of neurogenesis and angiogenesis. In the context of adult zebrafish spinal cord injury, we proved that Melanoma cell adhesion molecule (MCAM) is beneficial to the recovery, possibly via mechanisms of angiogenensis and inflammation. MCAM promotes angiogenesis by adjusting VEGFR-2, p-p38 and p-AKT. MCAM affects inflammatory factors such as TNF-α, IL-1ß and IL-8. Our results extend the beneficial role of MCAM in the regeneration of central nervous system.

The Upregulation of TRAF1 Induced by Helicobacter pylori Plays an Antiapoptotic Effect on the Infected Cells.

BACKGROUND: Tumor necrosis factor receptor-associated factor 1 (TRAF1) is a member of the TRAF family and is dysregulated in diseases, such as atheroma, lymphoma, and solid tumors, but the role of TRAF1 in gastric cancer remains unknown. This study was aimed to investigate the role of TRAF1 in Helicobacter pylori (H. pylori)-related cell apoptosis and gastric carcinogenesis. MATERIALS AND METHODS: The mRNA and protein expression levels of TRAF1 were measured in a panel of gastric cancer cell lines and in H. pylori -infected mice by quantitative real-time PCR (qPCR) and Western blotting. The transcription factor that mainly affects transcription of TRAF1 during H. pylori infection was identified. The roles of H. pylori virulence factors that regulate TRAF1 expression were analyzed using isogenic cagA-, vacA-, and cagE-null mutants. The effects of TRAF1 on gastric cell viability and apoptosis during H. pylori infection were detected using the standard MTS (cell viability) assay and flow cytometry, respectively. RESULTS: H. pylori infection induced TRAF1 overexpression in both gastric epithelial cells and mice. The expression of TRAF1 in response to H. pylori infection was majorly regulated by the activation of NF-κB and was strongly related to H. pylori virulence factor CagA. The upregulation of TRAF1 inhibited cell apoptosis and increased the viability of infected cells. CONCLUSIONS: H. pylori infection induces the overexpression of TRAF1 in gastric epithelial cells. The upregulation of TRAF1 plays an antiapoptotic role in H. pylori -infected gastric cells and may contribute to the gastric carcinogenesis.

Production and cell surface display of recombinant anthrax protective antigen on the surface layer of attenuated Bacillus anthracis.

To investigate the surface display of the anthrax protective antigen (PA) on attenuated Bacillus anthracis, a recombinant B. anthracis strain, named AP429 was constructed by integrating into the chromosome a translational fusion harboring the DNA fragments encoding the cell wall-targeting domain of the S-layer protein EA1 and the anthrax PA. Crerecombinase action at the loxP sites excised the antibiotic marker. Western blot analysis, fluorescence-activated cell sorting and immunofluorescence analysis confirmed that PA was successfully expressed on the S-layer of the recombinant antibiotic marker-free strain. Notwithstanding extensive proteolytic degradation of the hybrid protein SLHs-PA, quantitative ELISA revealed that approximately 8.1 × 10(6) molecules of SLHs-PA were gained from each Bacillus cell. Moreover, electron microscopy assay indicated that the typical S-layer structures could be clearly observed from the recombinant strain micrographs.