Syntaxin 6-mediated exosome secretion regulates enzalutamide resistance in prostate cancer.

Prostate cancer (PCa) deaths are typically the result of metastatic castration-resistant PCa (mCRPC). Recently, enzalutamide (Enz), an oral androgen receptor inhibitor, was approved for treating patients with mCRPC. Invariably, all PCa patients eventually develop resistance against Enz. Therefore, novel strategies aimed at overcoming Enz resistance are needed to improve the survival of PCa patients. The role of exosomes in drug resistance has not been fully elucidated in PCa. Therefore, we set out to better understand the exosome's role in the mechanism underlying Enz-resistant PCa. Results showed that Enz-resistant PCa cells (C4-2B, CWR-R1, and LNCaP) secreted significantly higher amounts of exosomes (2-4 folds) compared to Enz-sensitive counterparts. Inhibition of exosome biogenesis in resistant cells by GW4869 and dimethyl amiloride strongly decreased their cell viability. Mechanistic studies revealed upregulation of syntaxin 6 as well as its increased colocalization with CD63 in Enz-resistant PCa cells compared to Enz-sensitive cells. Syntaxin 6 knockdown by specific small interfering RNAs in Enz-resistant PCa cells (C4-2B and CWR-R1) resulted in reduced cell number and increased cell death in the presence of Enz. Furthermore, syntaxin 6 knockdown significantly reduced the exosome secretion in both Enz-resistant C4-2B and CWR-R1 cells. The Cancer Genome Atlas analysis showed increased syntaxin 6 expressions associated with higher Gleason score and decreased progression-free survival in PCa patients. Importantly, IHC analysis showed higher syntaxin 6 expression in cancer tissues from Enz-treated patients compared to Enz naïve patients. Overall, syntaxin 6 plays an important role in the secretion of exosomes and increased survival of Enz-resistant PCa cells.

Identification of Hub Genes Associated with Hypertension and Their Interaction with miRNA Based on Weighted Gene Coexpression Network Analysis (WGCNA) Analysis.

BACKGROUND Hypertension is one of the most widespread health conditions in the world, and the molecular mechanism of it is still unclear. In this study, we identified the hub genes (hub miRNA genes) associated with hypertension and explored the relationship between hypertension miRNA-gene by constructing a mRNA co-expression network and a miRNA co-expression network, which can help to reveal the mechanism and predict the prognosis of hypertension progression. MATERIAL AND METHODS Based on gene expression profile data of hypertensive samples from the Gene Expression Omnibus database, WGCNA was used to detect hypertension-related biomarkers and key mRNA and miRNA modules. Then, DAVID was used to perform gene-annotation enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) and miRPath were used for pathway analysis of mRNA and miRNAs genes. RESULTS We identified 3 key modules relating to hypertension, 2 mRNA modules named Msaddlebrown and Mgreenyellow and 1 miRNA module named Msalmon. In addition, 12 hub genes (RPL21, RPS28, LOC442727/PTGAP10, LOC100129599/RPS29P14, TBXAS1, FCER1G, CFP, FURIN, PECAM1, IGSF6, NCF1C, and LOC285296/UNC93B3) and 7 hub miRNAs (hsa-miR-1268a/b, hsa-miR-513c-3p, hsa-miR-4799-5p, hsa-miR-296-3p, hsa-miR-5195-5p, hsa-miR-219-2-3p, and hsa-miR-548d-5p) relating to hypertension were identified. HIF-1 signaling pathway and insulin signaling pathway were closely related to the 3 key modules. We also discovered 4 miRNAs (hsa-miR-548am-3p, hsa-miR-513c-3p, hsa-miR-182-5p, and hsa-miR-548d-5p) and 6 genes (IGF1R, GSK3B, FOXO1, PRKAR2B, HIF1A, and PIK3R1) were the core nodes in the hypertension-related miRNA-gene network, and hsa-miR-548am-3p was at the center of the network. CONCLUSIONS These findings will help improve the understanding of the pathogenesis of hypertension, and the discovered genes can serve as signatures for early diagnosis of hypertension.

Immune signaling-based Cascade Propagation approach re-stratifies HNSCC patients.

The availability of high-throughput genomic assays and rich electronic medical records allows us to identify cancer subtypes with greater accuracy and resolution. The integration of multiplatform, heterogenous, and high dimensional data remains an enormous challenge in using big data in bioinformatics research. Previous methods have been developed for patient stratification, however, these approaches did not incorporate prior knowledge and offer limited biology insight. New computational methods are needed to better utilize multiple types of information to identify clinically meaningful subtypes. Recent studies have shown that many immune functional genes are associated with cancer progression, recurrence and prognosis in head and neck squamous cell carcinoma (HNSCC). Therefore, we developed a novel immune signaling based Cascade Propagation (CasP) subtyping approach to stratify HNSCC patients. Unlike previous stratification methods that use only patient genomic data, our approach makes use of prior biological information such as immune signaling and protein-protein interactions, as well as patient survival information. CasP is a multi-step stratification procedure, composed of a dynamic network tree cutting step followed by a mutational stratification step. Using this approach, HNSCC patients were first stratified into clinically relative subgroups with different survival outcomes and distinct immunogenic features. We found that the good outcome of a subgroup of HNSCC patients was due to an enhanced immune response. The gene sets were characterized by a significant activation of T cell receptor signaling pathways, in addition to other important cancer related pathways such as PI3K and JAK/STAT signaling pathways. Further stratification of patients based on somatic mutation profiles detected three survival-distinct subnetworks. Our newly developed CasP subtyping approach allowed us to integrate multiple data types and identify clinically relevant subtypes of HNSCC patients.

A novel integrated approach to predicting cancer immunotherapy efficacy.

Immunotherapies have revolutionized cancer treatment modalities; however, predicting clinical response accurately and reliably remains challenging. Neoantigen load is considered as a fundamental genetic determinant of therapeutic response. However, only a few predicted neoantigens are highly immunogenic, with little focus on intratumor heterogeneity (ITH) in the neoantigen landscape and its link with different features in the tumor microenvironment. To address this issue, we comprehensively characterized neoantigens arising from nonsynonymous mutations and gene fusions in lung cancer and melanoma. We developed a composite NEO2IS to characterize interplays between cancer and CD8+ T-cell populations. NEO2IS improved prediction accuracy of patient responses to immune-checkpoint blockades (ICBs). We found that TCR repertoire diversity was consistent with the neoantigen heterogeneity under evolutionary selections. Our defined neoantigen ITH score (NEOITHS) reflected infiltration degree of CD8+ T lymphocytes with different differentiation states and manifested the impact of negative selection pressure on CD8+ T-cell lineage heterogeneity or tumor ecosystem plasticity. We classified tumors into distinct immune subtypes and examined how neoantigen-T cells interactions affected disease progression and treatment response. Overall, our integrated framework helps profile neoantigen patterns that elicit T-cell immunoreactivity, enhance the understanding of evolving tumor-immune interplays and improve prediction of ICBs efficacy.

The Potential Regulation of A-to-I RNA Editing on Genes in Parkinson's Disease.

Parkinson's disease (PD) is characterized by dopaminergic neurodegeneration and an abnormal accumulation of α-synuclein aggregates. A number of genetic factors have been shown to increase the risk of PD. Exploring the underlying molecular mechanisms that mediate PD's transcriptomic diversity can help us understand neurodegenerative pathogenesis. In this study, we identified 9897 A-to-I RNA editing events associated with 6286 genes across 372 PD patients. Of them, 72 RNA editing events altered miRNA binding sites and this may directly affect miRNA regulations of their host genes. However, RNA editing effects on the miRNA regulation of genes are more complex. They can (1) abolish existing miRNA binding sites, which allows miRNAs to regulate other genes; (2) create new miRNA binding sites that may sequester miRNAs from regulating other genes; or (3) occur in the miRNA seed regions and change their targets. The first two processes are also referred to as miRNA competitive binding. In our study, we found 8 RNA editing events that may alter the expression of 1146 other genes via miRNA competition. We also found one RNA editing event that modified a miRNA seed region, which was predicted to disturb the regulation of four genes. Considering the PD-related functions of the affected genes, 25 A-to-I RNA editing biomarkers for PD are proposed, including the 3 editing events in the EIF2AK2, APOL6, and miR-4477b seed regions. These biomarkers may alter the miRNA regulation of 133 PD-related genes. All these analyses reveal the potential mechanisms and regulations of RNA editing in PD pathogenesis.

DOTA: Deep Learning Optimal Transport Approach to Advance Drug Repositioning for Alzheimer's Disease.

Alzheimer's disease (AD) is the leading cause of age-related dementia, affecting over 5 million people in the United States and incurring a substantial global healthcare cost. Unfortunately, current treatments are only palliative and do not cure AD. There is an urgent need to develop novel anti-AD therapies; however, drug discovery is a time-consuming, expensive, and high-risk process. Drug repositioning, on the other hand, is an attractive approach to identify drugs for AD treatment. Thus, we developed a novel deep learning method called DOTA (Drug repositioning approach using Optimal Transport for Alzheimer's disease) to repurpose effective FDA-approved drugs for AD. Specifically, DOTA consists of two major autoencoders: (1) a multi-modal autoencoder to integrate heterogeneous drug information and (2) a Wasserstein variational autoencoder to identify effective AD drugs. Using our approach, we predict that antipsychotic drugs with circadian effects, such as quetiapine, aripiprazole, risperidone, suvorexant, brexpiprazole, olanzapine, and trazadone, will have efficacious effects in AD patients. These drugs target important brain receptors involved in memory, learning, and cognition, including serotonin 5-HT2A, dopamine D2, and orexin receptors. In summary, DOTA repositions promising drugs that target important biological pathways and are predicted to improve patient cognition, circadian rhythms, and AD pathogenesis.

Phenotype-Genotype analysis of caucasian patients with high risk of osteoarthritis.

Background: Osteoarthritis (OA) is a common cause of disability and pain around the world. Epidemiologic studies of family history have revealed evidence of genetic influence on OA. Although many efforts have been devoted to exploring genetic biomarkers, the mechanism behind this complex disease remains unclear. The identified genetic risk variants only explain a small proportion of the disease phenotype. Traditional genome-wide association study (GWAS) focuses on radiographic evidence of OA and excludes sex chromosome information in the analysis. However, gender differences in OA are multifactorial, with a higher frequency in women, indicating that the chromosome X plays an essential role in OA pathology. Furthermore, the prevalence of comorbidities among patients with OA is high, indicating multiple diseases share a similar genetic susceptibility to OA. Methods: In this study, we performed GWAS of OA and OA-associated key comorbidities on 3366 OA patient data obtained from the Osteoarthritis Initiative (OAI). We performed Mendelian randomization to identify the possible causal relationship between OA and OA-related clinical features. Results: One significant OA-associated locus rs2305570 was identified through sex-specific genome-wide association. By calculating the LD score, we found OA is positively correlated with heart disease and stroke. A strong genetic correlation was observed between knee OA and inflammatory disease, including eczema, multiple sclerosis, and Crohn's disease. Our study also found that knee alignment is one of the major risk factors in OA development, and we surprisingly found knee pain is not a causative factor of OA, although it was the most common symptom of OA. Conclusion: We investigated several significant positive/negative genetic correlations between OA and common chronic diseases, suggesting substantial genetic overlaps between OA and these traits. The sex-specific association analysis supports the critical role of chromosome X in OA development in females.

APAview: A web-based platform for alternative polyadenylation analyses in hematological cancers.

Background: Hematologic malignancies, such as acute promyelocytic leukemia (APL) and acute myeloid leukemia (AML), are cancers that start in blood-forming tissues and can affect the blood, bone marrow, and lymph nodes. They are often caused by genetic and molecular alterations such as mutations and gene expression changes. Alternative polyadenylation (APA) is a post-transcriptional process that regulates gene expression, and dysregulation of APA contributes to hematological malignancies. RNA-sequencing-based bioinformatic methods can identify APA sites and quantify APA usages as molecular indexes to study APA roles in disease development, diagnosis, and treatment. Unfortunately, APA data pre-processing, analysis, and visualization are time-consuming, inconsistent, and laborious. A comprehensive, user-friendly tool will greatly simplify processes for APA feature screening and mining. Results: Here, we present APAview, a web-based platform to explore APA features in hematological cancers and perform APA statistical analysis. APAview server runs on Python3 with a Flask framework and a Jinja2 templating engine. For visualization, APAview client is built on Bootstrap and Plotly. Multimodal data, such as APA quantified by QAPA/DaPars, gene expression data, and clinical information, can be uploaded to APAview and analyzed interactively. Correlation, survival, and differential analyses among user-defined groups can be performed via the web interface. Using APAview, we explored APA features in two hematological cancers, APL and AML. APAview can also be applied to other diseases by uploading different experimental data.

PredTAD: A machine learning framework that models 3D chromatin organization alterations leading to oncogene dysregulation in breast cancer cell lines.

Topologically associating domains, or TADs, play important roles in genome organization and gene regulation; however, they are often altered in diseases. High-throughput chromatin conformation capturing assays, such as Hi-C, can capture domains of increased interactions, and TADs and boundaries can be identified using well-established analytical tools. However, generating Hi-C data is expensive. In our study, we addressed the relationship between multi-omics data and higher-order chromatin structures using a newly developed machine-learning model called PredTAD. Our tool uses already-available and cost-effective datatypes such as transcription factor and histone modification ChIPseq data. Specifically, PredTAD utilizes both epigenetic and genetic features as well as neighboring information to classify the entire human genome as boundary or non-boundary regions. Our tool can predict boundary changes between normal and breast cancer genomes. Among the most important features for predicting boundary alterations were CTCF, subunits of cohesin (RAD21 and SMC3), and chromosome number, suggesting their roles in conserved and dynamic boundaries formation. Upon further analysis, we observed that genes near altered TAD boundaries were found to be involved in several important breast cancer signaling pathways such as Ras, Jak-STAT, and estrogen signaling pathways. We also discovered a TAD boundary alteration that contributes to RET oncogene overexpression. PredTAD can also successfully predict TAD boundary changes in other conditions and diseases. In conclusion, our newly developed machine learning tool allowed for a more complete understanding of the dynamic 3D chromatin structures involved in signaling pathway activation, altered gene expression, and disease state in breast cancer cells.

The Functional Impact of Alternative Splicing on the Survival Prognosis of Triple-Negative Breast Cancer.

PURPOSE: Triple-negative breast cancer (TNBC) is a type of breast cancer (BC) showing a high recurrence ratio and a low survival probability, which requires novel actionable molecular targets. The involvement of alternative splicing (AS) in TNBC promoted us to study the potential roles of AS events in the survival prognosis of TNBC patients. METHODS: A total of 150 TNBC patients from The Cancer Genome Atlas (TCGA) were involved in this work. To study the effects of AS in the recurrence-free survival (RFS) prognosis of TNBC, we performed the analyses as follows. First, univariate Cox regression model was applied to identify RFS-related AS events. Their host genes were analyzed by Metascape to discover the potential functions and involved pathways. Next, least absolute shrinkage and selection operator (LASSO) method was used to select the most informative RFS-related AS events to constitute an AS risk factor for RFS prognosis, which was evaluated by Kaplan-Meier (KM) and receiver operating characteristic (ROC) curves in all the data and also in different clinical subgroups. Furthermore, we analyzed the relationships between splicing factors (SFs) and these RFS-related AS events to seek the possibility that SFs regulated AS events to influence RFS. Then, we evaluated the potential of these RFS-related AS events in the overall survival (OS) prognosis from all the above aspects. RESULTS: We identified a total of 546 RFS-related AS events, which were enriched in some splicing and TNBC-associated pathways. Among them, seven RFS-related events were integrated into a risk factor, exhibiting satisfactory RFS prognosis alone and even better performance when combined with clinical tumor-node-metastasis stages. Furthermore, the correlation analysis between SFs and the seven AS events revealed the hypotheses that SRPK3 might upregulate PCYT2_44231_AA to have an effect on RFS prognosis and that three other SFs may work together to downregulate FLAD1_7874_RI to influence RFS prognosis. In addition, the seven RFS-related AS events were validated to be promising in the OS prognosis of TNBC as well. CONCLUSION: The abnormal AS events regulated by SFs may act as a kind of biomarker for the survival prognosis of TNBC.

Causal Discovery in Radiographic Markers of Knee Osteoarthritis and Prediction for Knee Osteoarthritis Severity With Attention-Long Short-Term Memory.

The goal of this study is to build a prognostic model to predict the severity of radiographic knee osteoarthritis (KOA) and to identify long-term disease progression risk factors for early intervention and treatment. We designed a long short-term memory (LSTM) model with an attention mechanism to predict Kellgren/Lawrence (KL) grade for knee osteoarthritis patients. The attention scores reveal a time-associated impact of different variables on KL grades. We also employed a fast causal inference (FCI) algorithm to estimate the causal relation of key variables, which will aid in clinical interpretability. Based on the clinical information of current visits, we accurately predicted the KL grade of the patient's next visits with 90% accuracy. We found that joint space narrowing was a major contributor to KOA progression. Furthermore, our causal structure model indicated that knee alignments may lead to joint space narrowing, while symptoms (swelling, grinding, catching, and limited mobility) have little impact on KOA progression. This study evaluated a broad spectrum of potential risk factors from clinical data, questionnaires, and radiographic markers that are rarely considered in previous studies. Using our statistical model, providers are able to predict the risk of the future progression of KOA, which will provide a basis for selecting proper interventions, such as proceeding to joint arthroplasty for patients. Our causal model suggests that knee alignment should be considered in the primary treatment and KOA progression was independent of clinical symptoms.

Syntaxin 6: A novel predictive and prognostic biomarker in papillary renal cell carcinoma.

Syntaxin 6 is a SNARE family protein known to play an important role in intracellular trafficking. Here, we examined the tumorogenic role of syntaxin 6 in renal cell carcinoma (RCC). The Cancer Genome Atlas (TCGA) was queried for clinicopathologic data and syntaxin 6 expression. We found a significant difference in overall survival (OS) between groups, with high syntaxin 6 expression correlating with decreased survival. When stratifying the data based on histological subtype, the papillary RCC subtype exhibited a significant correlation between syntaxin 6 expression and survival. Using ROC curve, we calculated the area under the curve (AUC) to determine the ability of syntaxin 6 to predict 3-year overall survival. The AUC for syntaxin 6 was 0.73, significantly higher compared to 0.52 for T stage. Next, syntaxin 6 expression was evaluated in clear cell (786-O and Caki-1) and papillary (Caki-2 and ACHN) RCC cells. Syntaxin 6 expression was higher in Caki-1 and ACHN RCC cells. Silencing of syntaxin 6 in ACHN cells significantly decreased the cell viability (p < 0.001). Overall, syntaxin 6 could be a prognostic biomarker for patients with papillary RCC and syntaxin 6 inhibitors hold promise as a novel therapy against RCC.

LSCC SNP variant regulates SOX2 modulation of VDAC3.

Lung squamous cell carcinoma (LSCC) is a genomically complex malignancy with no effective treatments. Recent studies have found a large number of DNA alterations such as SOX2 amplification in LSCC patients. As a stem cell transcription factor, SOX2 is important for the maintenance of pluripotent cells and may play a role in cancer. To study the downstream mechanisms of SOX2, we employed expression quantitative trait loci (eQTLs) technology to investigate how the presence of SOX2 affects the expression of target genes. We discovered unique eQTLs, such as rs798827-VDAC3 (FDR p-value = 0.0034), that are only found in SOX2-active patients but not in SOX2-inactive patients. SNP rs798827 is within strong linkage disequilibrium (r2 = 1) to rs58163073, where rs58163073 [T] allele increases the binding affinity of SOX2 and allele [TA] decreases it. In our analysis, SOX2 silencing downregulates VDAC3 in two LSCC cell lines. Chromatin conformation capturing data indicates that this SNP is located within the same Topologically Associating Domain (TAD) of VDAC3, further suggesting SOX2's role in the regulation of VDAC3 through the binding of rs58163073. By first subgrouping patients based on SOX2 activity, we made more relevant eQTL discoveries and our analysis can be applied to other diseases.

Genome-Wide Association and Mechanistic Studies Indicate That Immune Response Contributes to Alzheimer's Disease Development.

Alzheimer's disease (AD) is the most common cause of dementia. Although genome-wide association study (GWAS) have reported hundreds of single-nucleotide polymorphisms (SNPs) and genes linked to AD, the mechanisms about how these SNPs modulate the development of AD remain largely unknown. In this study, we performed GWAS for three traits in cerebrospinal fluid (CSF) and one clinical trait in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. Our analysis identified five most significant AD related SNPs (FDR < 0.05) within or proximal to APOE, APOC1, and TOMM40. One of the SNPs was co-inherited with APOE allele 4, which is the most important genetic risk factor for AD. Three of the five SNPs were located in promoter or enhancer regions, and transcription factor (TF) binding affinity calculations showed dramatic changes (| Log2FC| > 2) of three TFs (PLAG1, RREB1, and ZBTB33) for two motifs containing SNPs rs2075650 and rs157580. In addition, our GWAS showed that both rs2075650 and rs157580 were significantly associated with the poliovirus receptor-related 2 (PVRL2) gene (FDR < 0.25), which is involved in spreading of herpes simplex virus (HSV). The altered regulation of PVRL2 may increase the susceptibility AD patients to HSV and other virus infections of the brain. Our work suggests that AD is a type of immune disorder driven by viral or microbial infections of the brain during aging.