RNALocate: a resource for RNA subcellular localizations.

Increasing evidence has revealed that RNA subcellular localization is a very important feature for deeply understanding RNA's biological functions after being transported into intra- or extra-cellular regions. RNALocate is a web-accessible database that aims to provide a high-quality RNA subcellular localization resource and facilitate future researches on RNA function or structure. The current version of RNALocate documents more than 37 700 manually curated RNA subcellular localization entries with experimental evidence, involving more than 21 800 RNAs with 42 subcellular localizations in 65 species, mainly including Homo sapiens, Mus musculus and Saccharomyces cerevisiae etc. Besides, RNA homology, sequence and interaction data have also been integrated into RNALocate. Users can access these data through online search, browse, blast and visualization tools. In conclusion, RNALocate will be of help in elucidating the entirety of RNA subcellular localization, and developing new prediction methods. The database is available at http://www.rna-society.org/rnalocate/.

ViRBase: a resource for virus-host ncRNA-associated interactions.

Increasing evidence reveals that diverse non-coding RNAs (ncRNAs) play critically important roles in viral infection. Viruses can use diverse ncRNAs to manipulate both cellular and viral gene expression to establish a host environment conducive to the completion of the viral life cycle. Many host cellular ncRNAs can also directly or indirectly influence viral replication and even target virus genomes. ViRBase (http://www.rna-society.org/virbase) aims to provide the scientific community with a resource for efficient browsing and visualization of virus-host ncRNA-associated interactions and interaction networks in viral infection. The current version of ViRBase documents more than 12,000 viral and cellular ncRNA-associated virus-virus, virus-host, host-virus and host-host interactions involving more than 460 non-redundant ncRNAs and 4400 protein-coding genes from between more than 60 viruses and 20 hosts. Users can query, browse and manipulate these virus-host ncRNA-associated interactions. ViRBase will be of help in uncovering the generic organizing principles of cellular virus-host ncRNA-associated interaction networks in viral infection.

RAID: a comprehensive resource for human RNA-associated (RNA-RNA/RNA-protein) interaction.

Transcriptomic analyses have revealed an unexpected complexity in the eukaryote transcriptome, which includes not only protein-coding transcripts but also an expanding catalog of noncoding RNAs (ncRNAs). Diverse coding and noncoding RNAs (ncRNAs) perform functions through interaction with each other in various cellular processes. In this project, we have developed RAID (http://www.rna-society.org/raid), an RNA-associated (RNA-RNA/RNA-protein) interaction database. RAID intends to provide the scientific community with all-in-one resources for efficient browsing and extraction of the RNA-associated interactions in human. This version of RAID contains more than 6100 RNA-associated interactions obtained by manually reviewing more than 2100 published papers, including 4493 RNA-RNA interactions and 1619 RNA-protein interactions. Each entry contains detailed information on an RNA-associated interaction, including RAID ID, RNA/protein symbol, RNA/protein categories, validated method, expressing tissue, literature references (Pubmed IDs), and detailed functional description. Users can query, browse, analyze, and manipulate RNA-associated (RNA-RNA/RNA-protein) interaction. RAID provides a comprehensive resource of human RNA-associated (RNA-RNA/RNA-protein) interaction network. Furthermore, this resource will help in uncovering the generic organizing principles of cellular function network.

Mapping enhancer and chromatin accessibility landscapes charts the regulatory network of Alzheimer's disease.

BACKGROUND: Enhancers are regulatory elements that target and modulate gene expression and play a role in human health and disease. However, the roles of enhancer regulatory circuit abnormalities driven by epigenetic alterations in Alzheimer's disease (AD) are unclear. METHODS: In this study, a multiomic integrative analysis was performed to map enhancer and chromatin accessibility landscapes and identify regulatory network abnormalities in AD. We identified differentially methylated enhancers and constructed regulatory networks across brain regions using AD brain tissue samples. Through the integration of snATAC-seq and snRNA-seq datasets, we mapped enhancers with DNA methylation alterations (DMA) and chromatin accessibility landscapes. Core regulatory triplets that contributed to AD neuropathology in specific cell types were further prioritized. RESULTS: We revealed widespread DNA methylation alterations (DMA) in the enhancers of AD patients across different brain regions. In addition, the genome-wide transcription factor (TF) binding profiles showed that enhancers with DMA are pervasively regulated by TFs. The TF-enhancer-gene regulatory network analysis identified core regulatory triplets that are associated with brain and immune cell proportions and play important roles in AD pathogenesis. Enhancer regulatory circuits with DMA exhibited distinct chromatin accessibility patterns, which were further characterized at single-cell resolutions. CONCLUSIONS: Our study comprehensively investigated DNA methylation-mediated regulatory circuit abnormalities and provided novel insights into the potential pathogenesis of AD.

Activation of the hypoxia response in the aging cerebrovasculature protects males against cognitive impairment.

Alzheimer's disease (AD) is a neurodegenerative disorder with a distinct sex bias. Age-related vascular alterations, a hallmark of AD onset and progression, are consistently associated with sexual dimorphism. Here, we conducted an integrative meta-analysis of 335,803 single-nucleus transcriptomes and 667 bulk transcriptomes from the vascular system in AD and normal aging to address the underlying sex-dependent vascular aging in AD. All vascular cell types in male AD patients exhibited an activated hypoxia response and downstream signaling pathways including angiogenesis. The female AD vasculature is characterized by increased antigen presentation and decreased angiogenesis. We further confirmed that these sex-biased alterations in the cerebral vascular emerged and were primarily determined in the early stages of AD. Sex-stratified analysis of normal vascular aging revealed that angiogenesis and various stress-response genes were downregulated concurrently with female aging. Conversely, the hypoxia response increased steadily in males upon aging. An investigation of upstream driver transcription factors (TFs) revealed that altered communication between estrogen receptor alpha (ESR1) and hypoxia induced factors during menopause contributes to the inhibition of angiogenesis during normal female vascular aging. Additionally, inhibition of CREB1, a TF that targets estrogen, is also related to female AD. Overall, our study revealed a distinct cerebral vascular profile in females and males, and revealed novel targets for precision medicine therapy for AD.

Development of an m6A-Related lncRNAs Signature Predicts Tumor Stemness and Prognosis for Low-Grade Glioma Patients.

Background: Growing evidence has revealed that m6A modification of long noncoding RNAs (lncRNAs) dynamically controls tumor stemness and tumorigenesis-related processes. However, the prognostic significance of m6A-related lncRNAs and their associations with stemness in low-grade glioma (LGG) remain to be clarified. Methods: A multicenter transcriptome analysis of lncRNA expression in 1,247 LGG samples was performed in this study. The stemness landscape of LGG tumors was presented and associations with clinical features were revealed. The m6A-related lncRNAs were identified between stemness groups and were further prioritized via least absolute shrinkage and selection operator Cox regression analysis. A risk score model based on m6A-related lncRNAs was constructed and validated in external LGG datasets. Results: Based on the expression of LINC02984, PFKP-DT, and CRNDE, a risk model and nomogram were constructed; they successfully predicted the survival of patients and were extended to external datasets. Significant correlations were observed between the risk score and tumor stemness. Moreover, patients in different risk groups exhibited distinct tumor immune microenvironments and immune signatures. We finally provided several potential compounds suitable for specific risk groups, which may aid in LGG treatment. Conclusions: This novel signature presents noteworthy value in the prediction of prognosis and stemness status for LGG patients and will foster future research on the development of clinical regimens.

Identification and functional analysis of N6-methyladenine (m6 A)-related lncRNA across 33 cancer types.

BACKGROUND: N6-methyladenosine (m6 A) plays an essential role in tumorigenesis and cancer progression. Long noncoding RNAs (lncRNAs) are discovered to be important targets of m6 A modification, and they play fundamental roles in diverse biological processes. However, there is still a lack of knowledge with regards to the association between m6 A and lncRNAs in human tumors. METHODS: The relationship between lncRNAs and 21 m6 A regulators was comprehensively explored, through the integration of multi-omics data from M6A2Target, m6A-Atlas, and TCGA (The Cancer Genome Atlas). In order to explore the potential roles of m6A-related lncRNAs in human tumors, three applicable methods were introduced, which include the construction of ceRNA networks, drug sensitivity estimation, and survival analysis. RESULTS: A substantial number of positive correlation events across 33 cancer types were found. Moreover, cancer-specific lncRNAs were associated with tissue specificity, and cancer-common lncRNAs were conserved in cancer-related biological function. In particular, the m6 A-related lncRNA FGD5-AS1 was found to be associated with cancer treatment, through its influence on cisplatin resistance in breast cancer patients. Finally, a user-friendly interface Lnc2m6A, which is enriched with various browsing sections resource for the exhibition of relationships and putative biogenesis between lncRNAs and m6 A modifications, is offered in http://hainmu-biobigdata.com/Lnc2m6A. CONCLUSIONS: In summary, the results from this paper will provide a valuable resource that guides both mechanistic and therapeutic roles of m6 A-related lncRNAs in human tumors.

Enhanced insulin-regulated phagocytic activities support extreme health span and longevity in multiple populations.

The immune system plays a central role in many processes of age-related disorders and it remains unclear if the innate immune system may play roles in shaping extreme longevity. By an integrated analysis with multiple bulk and single cell transcriptomic, so as DNA methylomic datasets of white blood cells, a previously unappreciated yet commonly activated status of the innate monocyte phagocytic activities is identified. Detailed analyses revealed that the life cycle of these monocytes is enhanced and primed to a M2-like macrophage phenotype. Functional characterization unexpectedly revealed an insulin-driven immunometabolic network which supports multiple aspects of phagocytosis. Such reprogramming is associated to a skewed trend of DNA demethylation at the promoter regions of multiple phagocytic genes, so as a direct transcriptional effect induced by nuclear-localized insulin receptor. Together, these highlighted that preservation of insulin sensitivity is a key to healthy lifespan and extended longevity, via boosting the function of innate immune system in advanced ages.

Epigenetically regulated lncRNAs dissect the intratumoural heterogeneity and facilitate immune evasion of glioblastomas.

Background: Glioblastomas are the most common and malignant central nervous system (CNS) tumors that occupied a highly heterogeneous tumor microenvironment (TIME). Long noncoding RNAs (lncRNAs), whose expression can be modified by DNA methylation, are emerging as critical regulators in the immune system. However, knowledge about the epigenetic changes in lncRNAs and their contribution to the immune heterogeneity of glioma is still lacking. Methods: In this study, we integrated paired methylome and transcriptome datasets of glioblastomas and identified 2 robust immune subtypes based on lncRNA methylation features. The immune characteristics of glioma subtypes were compared. Furthermore, immune-related lncRNAs were identified and their relationships with immune evasion were evaluated. Results: Glioma immunophenotypes exhibited distinct immune-related characteristics as well as clinical and epigenetic features. 149 epigenetically regulated (ER) lncRNAs were recognized that possessed inverse variation in epigenetic and transcriptional levels between glioma subtypes. Immune-related lncRNAs were further identified through the investigation of their correlation with immune cell infiltrations and immune-related pathways. In particular, the 'Hot' glioma subtype with higher immunoactivity while a worse survival outcome was found to character immune evasion features. We finally prioritized candidate ER lncRNAs associated with immune evasion markers and response to glioma immunotherapy. Among them, CD109-AS1 and LINC02447 were validated as novel immunoevasive biomarkers for glioma through in vitro experiments. Conclusion: In summary, our study systematically reveals the crosstalk among DNA methylation, lncRNA, and immune regulation in glioblastomas, and will facilitate the development of epigenetic immunotherapy approaches.

Pan-cancer analyses reveal multi-omic signatures and clinical implementations of the forkhead-box gene family.

BACKGROUND: Forkhead box (FOX) proteins belong to one of the largest transcription factor families and play crucial roles in the initiation and progression of cancer. Prior research has linked several FOX genes, such as FOXA1 and FOXM1, to the crucial process of carcinogenesis. However, the overall picture of FOX gene family across human cancers is far from clear. METHODS: To investigate the broad molecular signatures of the FOX gene family, we conducted study on multi-omics data (including genomics, epigenomics and transcriptomics) from over 11,000 patients with 33 different types of human cancers. RESULTS: Pan-cancer analysis reveals that FOX gene mutations were found in 17.4% of tumor patients with a substantial cancer type-dependent pattern. Additionally, high expression heterogeneity of FOX genes across cancer types was discovered, which can be partially attributed to the genomic or epigenomic alteration. Co-expression network analysis reveals that FOX genes may exert functions by regulating the expression of both their own and target genes. For a clinical standpoint, we provided 103 FOX gene-drug target-drug predictions and found FOX gene expression have potential survival predictive value. All of the results have been included in the FOX2Cancer database, which is freely accessible at http://hainmu-biobigdata.com/FOX2Cancer. CONCLUSION: Our findings may provide a better understanding of roles FOX genes played in the development of tumors, and help to offer new avenues for uncovering tumorigenesis and unprecedented therapeutic targets.

Pediatric Pan-Central Nervous System Tumor Methylome Analyses Reveal Immune-Related LncRNAs.

Pediatric central nervous system (CNS) tumors are the second most common cancer diagnosis among children. Long noncoding RNAs (lncRNAs) emerge as critical regulators of gene expression, and they play fundamental roles in immune regulation. However, knowledge on epigenetic changes in lncRNAs in diverse types of pediatric CNS tumors is lacking. Here, we integrated the DNA methylation profiles of 2,257 pediatric CNS tumors across 61 subtypes with lncRNA annotations and presented the epigenetically regulated landscape of lncRNAs. We revealed the prevalent lncRNA methylation heterogeneity across pediatric pan-CNS tumors. Based on lncRNA methylation profiles, we refined 14 lncRNA methylation clusters with distinct immune microenvironment patterns. Moreover, we found that lncRNA methylations were significantly correlated with immune cell infiltrations in diverse tumor subtypes. Immune-related lncRNAs were further identified by investigating their correlation with immune cell infiltrations and potentially regulated target genes. LncRNA with methylation perturbations potentially regulate the genes in immune-related pathways. We finally identified several candidate immune-related lncRNA biomarkers (i.e., SSTR5-AS1, CNTN4-AS1, and OSTM1-AS1) in pediatric cancer for further functional validation. In summary, our study represents a comprehensive repertoire of epigenetically regulated immune-related lncRNAs in pediatric pan-CNS tumors, and will facilitate the development of immunotherapeutic targets.

Comprehensive characterization of human-virus protein-protein interactions reveals disease comorbidities and potential antiviral drugs.

The protein-protein interactions (PPIs) between human and viruses play important roles in viral infection and host immune responses. Rapid accumulation of experimentally validated human-virus PPIs provides an unprecedented opportunity to investigate the regulatory pattern of viral infection. However, we are still lack of knowledge about the regulatory patterns of human-virus interactions. We collected 27,293 experimentally validated human-virus PPIs, covering 8 virus families, 140 viral proteins and 6059 human proteins. Functional enrichment analysis revealed that the viral interacting proteins were likely to be enriched in cell cycle and immune-related pathways. Moreover, we analysed the topological features of the viral interacting proteins and found that they were likely to locate in central regions of human PPI network. Based on network proximity analyses of diseases genes and human-virus interactions in the human interactome, we revealed the associations between complex diseases and viral infections. Network analysis also implicated potential antiviral drugs that were further validated by text mining. Finally, we presented the Human-Virus Protein-Protein Interaction database (HVPPI, http://bio-bigdata.hrbmu.edu.cn/HVPPI), that provides experimentally validated human-virus PPIs as well as seamlessly integrates online functional analysis tools. In summary, comprehensive understanding the regulatory pattern of human-virus interactome will provide novel insights into fundamental infectious mechanism discovery and new antiviral therapy development.

Prognostic Implications and Immune Infiltration Characteristics of Chromosomal Instability-Related Dysregulated CeRNA in Lung Adenocarcinoma.

An accumulating body of research indicates that long-noncoding RNAs (lncRNAs) regulate the target genes and act as competitive endogenous RNAs (ceRNAs) playing an indispensable role in lung adenocarcinoma (LUAD). LUAD is frequently accompanied by the feature of chromosomal instability (CIN); however, CIN-related ceRNAs have not been investigated yet. We systematically analyzed and integrated CIN-related dysregulated ceRNAs characteristics in LUAD samples for the first time. In TCGA LUAD cohort, CIN in tumor samples was significantly higher than that in those of adjacent, and patients with high CIN risk tended to have worse clinical outcomes. We constructed a double-weighted CIN-related dysregulated ceRNA network, in which edge weight and node weight represented the disorder extent of ceRNA and the correlation of RNA expression level and prognosis, respectively. After module mining and analysis, a potential prognostic biomarker composed of 12 RNAs (8 mRNAs and 4 lncRNAs) named CIN-related dysregulated ceRNAs (CRDC) was obtained. The CRDC risk score had a positive relation with clinical stage and CIN, and patients with high CRDC risk scores exhibited poor prognosis. Moreover, CRDC tended to be an independent risk factor with high robustness to overcome the effect of multicollinearity among other explanatory variables for disease-specific survival (DSS) in TCGA and two GEO cohorts. The result of functional analysis indicated that CRDC was involved in multiple cancer progresses, especially immune-related pathways. The patients with lower CRDC risk had higher B cell, T cell CD4+, T cell CD8+, neutrophil, macrophage, and myeloid dendritic cell infiltration than the patients with higher CRDC risk. Meanwhile, patients with lower CRDC risk could get more benefits from immunological therapy. The results suggested that the CRDC could be a potential prognostic biomarker and an immunotherapy predictor for lung adenocarcinoma.

Comprehensive Analysis of the Carcinogenic Process, Tumor Microenvironment, and Drug Response in HPV-Positive Cancers.

Human papillomavirus (HPV) is a common virus, and about 5% of all cancers worldwide is caused by persistent high-risk HPV infections. Here, we reported a comprehensive analysis of the molecular features for HPV-related cancer types using TCGA (The Cancer Genome Atlas) data with HPV status. We found that the HPV-positive cancer patients had a unique oncogenic process, tumor microenvironment, and drug response compared with HPV-negative patients. In addition, HPV improved overall survival for the four cancer types, namely, cervical squamous cell carcinoma (CESC), head and neck squamous cell carcinoma (HNSC), stomach adenocarcinoma (STAD), and uterine corpus endometrial carcinoma (UCEC). The stronger activity of cell-cycle pathways and lower driver gene mutation rates were observed in HPV-positive patients, which implied the different carcinogenic processes between HPV-positive and HPV-negative groups. The increased activities of immune cells and differences in metabolic pathways helped explain the heterogeneity of prognosis between the two groups. Furthermore, we constructed HPV prediction models for different cancers by the virus infection score (VIS) which was linearly correlated with HPV load and found that VIS was associated with drug response. Altogether, our study reveals that HPV-positive cancer patients have unique molecular characteristics which help the development of precision medicine in HPV-positive cancers.

HCDT: an integrated highly confident drug-target resource.

Drug-target association plays an important role in drug discovery, drug repositioning, drug synergy prediction, etc. Currently, a lot of drug-related databases, such as DrugBank and BindingDB, have emerged. However, these databases are separate, incomplete and non-uniform with different criteria. Here, we integrated eight drug-related databases; collected, filtered and supplemented drugs, target genes and experimentally validated (highly confident) associations and built a highly confident drug-target (HCDT: http://hainmu-biobigdata.com/hcdt) database. HCDT database includes 500 681 HCDT associations between 299 458 drugs and 5618 target genes. Compared to individual databases, HCDT database contains 1.1 to 254.2 times drugs, 1.8-5.5 times target genes and 1.4-27.7 times drug-target associations. It is normative, publicly available and easy for searching, browsing and downloading. Together with multi-omics data, it will be a good resource in analyzing the drug functional mechanism, mining drug-related biological pathways, predicting drug synergy, etc. Database URL: http://hainmu-biobigdata.com/hcdt.

The synergistic interaction landscape of chromatin regulators reveals their epigenetic regulation mechanisms across five cancer cell lines.

Chromatin regulators (CRs) regulate the gene transcription process through combinatorial patterns, which currently remain obscure for pan-cancer. This study identified the interaction of CRs and constructed CR-CR interaction networks across five tumor cell lines. The global interaction analysis revealed that CRs tend to function in synergistically. In addition, common and specific CRs in interaction networks were identified, and the epigenetic processes of these CRs in regulating gene transcription were analyzed. Common CRs have conserved binding sites but cooperate with different partners in multiple tumor cell lines. They also participate in gene transcription regulation, through mediation of different histone modifications (HMs). Specific CRs, ATF2 and PRDM10 were found to distinguish liver cancer samples with different prognosis. PRDM10 participates in gene transcription regulation, by exertion of influence on the DNA methylation level of liver cancer. Through analysis of the edges in the CR-CR interaction networks, it was found EP300-TAF1 has genome-wide distinct signaling patterns, which exhibit different effects on downstream targets. This analysis provides novel insights for the understanding of synergistic mechanism of CRs function, as controllers of gene transcription across cancer types.

DysPIA: A Novel Dysregulated Pathway Identification Analysis Method.

Differential co-expression-based pathway analysis is still limited and not widely used. In most current methods, the pathways were considered as gene sets, but the gene regulation relationships were not considered, and the computational speed was slow. In this article, we proposed a novel Dysregulated Pathway Identification Analysis (DysPIA) method to overcome these shortcomings. We adopted the idea of Correlation by Individual Level Product into analysis and performed a fast enrichment analysis. We constructed a combined gene-pair background which was much more sufficient than the background used in Edge Set Enrichment Analysis. In simulation study, DysPIA was able to identify the causal pathways with high AUC (0.9584 to 0.9896). In p53 mutation data, DysPIA obtained better performance than other methods. It obtained more potential dysregulated pathways that could be literature verified, and it ran much faster (∼1,700-8,000 times faster than other methods when 10,000 permutations). DysPIA was also applied to breast cancer relapse dataset and breast cancer subtype dataset. The results show that DysPIA is effective and has a great biological significance. R packages "DysPIA" and "DysPIAData" are constructed and freely available on R CRAN (https://cran.r-project.org/web/packages/DysPIA/index.html and https://cran.r-project.org/web/packages/DysPIAData/index.html), and on GitHub (https://github.com/lemonwang2020).

The Functional Characterization of Epigenetically Related lncRNAs Involved in Dysregulated CeRNA-CeRNA Networks Across Eight Cancer Types.

Numerous studies have demonstrated that lncRNAs could compete with other RNAs to bind miRNAs, as competing endogenous RNAs (ceRNAs), to regulate each other. On the other hand, ceRNAs were found to be recurrently dysregulated in cancer status. However, limited studies considered the upstream epigenetic regulatory factors that disrupted the normal competing mechanism. In the present study, we constructed the lncRNA-associated dysregulated ceRNA networks across eight cancer types. lncRNAs in the individual dysregulated network and pan-cancer core dysregulated ceRNA subnetwork were found to play more important roles than mRNAs. Integrating lncRNA methylation profiles, we identified 49 epigenetically related (ER) lncRNAs involved in the dysregulated ceRNA networks, including 18 epigenetically activated (EA) lncRNAs, 18 epigenetically silenced (ES) lncRNAs, and 13 rewired ER lncRNAs across eight cancer types. Furthermore, we evaluated the epigenetic regulating patterns of these lncRNAs and screened nine pan-cancer ER lncRNAs (six EA and three ES lncRNAs). The nine lncRNAs were found to regulate the cancer hallmarks by competing with mRNAs. Moreover, we found that integrating the expression and methylation profiles of the nine lncRNAs could predict cancer incidence in eight cancer types robustly and the cancer outcome of several cancer types. These results provide an improved understanding of methylation regulation to ceRNA and offer novel potential molecular therapeutic targets for the diagnosis and prognosis across different cancer types.

Whole Blood Transcriptome Analysis Reveals the Correlation between Specific Immune Cells and Septicemic Melioidosis.

Melioidosis is a serious infectious disease caused by the environmental Gram-negative bacillus Burkholderia pseudomallei. It has been shown that the host immune system, mainly comprising various types of immune cells, fights against the disease. The present study was to specify correlation between septicemic melioidosis and the levels of multiple immune cells. First, the genes with differential expression patterns between patients with septicemic melioidosis (B. pseudomallei) and health donors (control/healthy) were identified. These genes being related to cytokine binding, cell adhesion molecule binding, and MHC relevant proteins may influence immune response. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed 23 enriched immune response pathways. We further leveraged the microarray data to investigate the relationship between immune response and septicemic melioidosis, using the CIBERSORT analysis. Comparison of the percentages of 22 immune cell types in B. pseudomallei vs. control/healthy revealed that those of CD4 memory resting cells, CD8+ T cells, B memory cells, and CD4 memory activated cells were low, whereas those of M0 macrophages, neutrophils, and gamma delta T cells were high. The multivariate logistic regression analysis further revealed that CD8+ T cells, M0 macrophages, neutrophils, and naive CD4+ cells were strongly associated with the onset of septicemic melioidosis, and M2 macrophages and neutrophils were associated with the survival in septicemic melioidosis. Taken together, these data point to a complex role of immune cells on the development and progression of melioidosis.

A Pan-Cancer Analysis of Cystatin E/M Reveals Its Dual Functional Effects and Positive Regulation of Epithelial Cell in Human Tumors.

Cystatin E/M (CST6), a representative cysteine protease inhibitor, plays both tumor-promoting and tumor-suppressing functions and is pursued as an epigenetically therapeutic target in special cancer types. However, a comprehensive and systematic analysis for CST6 in pan-cancer level is still lacking. In the present study, we explored the expression pattern of CST6 in multiple cancer types across ∼10,000 samples from TCGA (The Cancer Genome Atlas) and ∼8,000 samples from MMDs (Merged Microarray-acquired Datasets). We found that the dynamic expression alteration of CST6 was consistent with dual function in different types of cancer. In addition, we observed that the expression of CST6 was globally regulated by the DNA methylation in its promoter region. CST6 expression was positively correlated with the epithelial cell infiltration involved in epithelial-to-mesenchymal transition (EMT) and proliferation. The relationship between CST6 and tumor microenvironment was also explored. In particular, we found that CST6 serves a protective function in the process of melanoma metastasis. Finally, the clinical association analysis further revealed the dual function of CST6 in cancer, and a combination of the epithelial cell infiltration and CST6 expression could predict the prognosis for SKCM patients. In summary, this first CST6 pan-cancer study improves the understanding of the dual functional effects on CST6 in different types of human cancer.