Polygenic risk score predicting susceptibility and outcome of benign prostatic hyperplasia in the Han Chinese.

BACKGROUND: Given the high prevalence of BPH among elderly men, pinpointing those at elevated risk can aid in early intervention and effective management. This study aimed to explore that polygenic risk score (PRS) is effective in predicting benign prostatic hyperplasia (BPH) incidence, prognosis and risk of operation in Han Chinese. METHODS: A retrospective cohort study included 12,474 male participants (6,237 with BPH and 6,237 non-BPH controls) from the Taiwan Precision Medicine Initiative (TPMI). Genotyping was performed using the Affymetrix Genome-Wide TWB 2.0 SNP Array. PRS was calculated using PGS001865, comprising 1,712 single nucleotide polymorphisms. Logistic regression models assessed the association between PRS and BPH incidence, adjusting for age and prostate-specific antigen (PSA) levels. The study also examined the relationship between PSA, prostate volume, and response to 5-α-reductase inhibitor (5ARI) treatment, as well as the association between PRS and the risk of TURP. RESULTS: Individuals in the highest PRS quartile (Q4) had a significantly higher risk of BPH compared to the lowest quartile (Q1) (OR = 1.51, 95% CI = 1.274-1.783, p < 0.0001), after adjusting for PSA level. The Q4 group exhibited larger prostate volumes and a smaller volume reduction after 5ARI treatment. The Q1 group had a lower cumulative TURP probability at 3, 5, and 10 years compared to the Q4 group. PRS Q4 was an independent risk factor for TURP. CONCLUSIONS: In this Han Chinese cohort, higher PRS was associated with an increased susceptibility to BPH, larger prostate volumes, poorer response to 5ARI treatment, and a higher risk of TURP. Larger prospective studies with longer follow-up are warranted to further validate these findings.

Whole transcriptome analysis to explore the impaired immunological features in critically ill elderly patients with sepsis.

BACKGROUND: Sepsis is a frequent complication in critically ill patients, is highly heterogeneous and is associated with high morbidity and mortality rates, especially in the elderly population. Utilizing RNA sequencing (RNA-Seq) to analyze biological pathways is widely used in clinical and molecular genetic studies, but studies in elderly patients with sepsis are still lacking. Hence, we investigated the mortality-relevant biological features and transcriptomic features in elderly patients who were admitted to the intensive care unit (ICU) for sepsis. METHODS: We enrolled 37 elderly patients with sepsis from the ICU at Taichung Veterans General Hospital. On day-1 and day-8, clinical and laboratory data, as well as blood samples, were collected for RNA-Seq analysis. We identified the dynamic transcriptome and enriched pathways of differentially expressed genes between day-8 and day-1 through DVID enrichment analysis and Gene Set Enrichment Analysis. Then, the diversity of the T cell repertoire was analyzed with MiXCR. RESULTS: Overall, 37 patients had sepsis, and responders and non-responders were grouped through principal component analysis. Significantly higher SOFA scores at day-7, longer ventilator days, ICU lengths of stay and hospital mortality were found in the non-responder group, than in the responder group. On day-8 in elderly ICU patients with sepsis, genes related to innate immunity and inflammation, such as ZDHCC19, ALOX15, FCER1A, HDC, PRSS33, and PCSK9, were upregulated. The differentially expressed genes (DEGs) were enriched in the regulation of transcription, adaptive immune response, immunoglobulin production, negative regulation of transcription, and immune response. Moreover, there was a higher diversity of T-cell receptors on day-8 in the responder group, than on day-1, indicating that they had better regulated recovery from sepsis compared with the non-response patients. CONCLUSION: Sepsis mortality and incidence were both high in elderly individuals. We identified mortality-relevant biological features and transcriptomic features with functional pathway and MiXCR analyses based on RNA-Seq data; and found that the responder group had upregulated innate immunity and increased T cell diversity; compared with the non-responder group. RNA-Seq may be able to offer additional complementary information for the accurate and early prediction of treatment outcome.

Chemoresistance in acute myeloid leukemia: An alternative single-cell RNA sequencing approach.

Our previous study demonstrated that myc, mitochondrial oxidative phosphorylation, mTOR, and stemness are independently responsible for chemoresistance in acute myeloid leukemia (AML) cells. This study aimed to identify potential mechanisms of chemoresistance of the "7 + 3" induction in AML by using a single-cell RNA sequencing (scRNA-seq) approach. In the present study, 13 untreated patients with de novo AML were enrolled and stratified into two groups: complete remission (CR; n = 8) and non-CR (n = 5). Single-cell RNA sequencing was used to analyze genetic profiles of 28,950 AML cells from these patients; results were validated using a previously published bulk RNA-seq dataset. Our study results showed chemoresistant AML cells had premature accumulation during early hematopoiesis. Hematopoietic stem cell-like cells from the non-CR group expressed more leukemic stem cell markers (CD9, CD82, IL3RA, and IL1RAP) than those from the CR group. Chemoresistant progenitor cells had impaired myeloid differentiation owing to early arrest of hematopoiesis. Notably, AML cells analyzed by scRNA-seq and bulk RNA-seq harbored a comparable myeloid lineage cell fraction, which internally validated our results. Using the TCGA database, our analysis demonstrated that patients with AML with higher expression of chemoresistant genetic markers (IL3RA and IL1RAP) had a worse overall survival (p < 0.01 for IL3RA; p < 0.05 for IL1RAP). In conclusion, AML cells responsive and resistant to the "7 + 3" induction were derived from a diverse cancerous hematopoietic stem cell population, as indicated by the specific genetic biomarkers obtained using scRNA-seq approach. Furthermore, arrest of hematopoiesis was shown to occur earlier in chemoresistant AML cells, furthering the current understanding of chemoresistance in AML.

Multi-ethnic Imputation System (MI-System): A genotype imputation server for high-dimensional data.

OBJECTIVE: Genotype imputation is a commonly used technique that infers un-typed variants into a study's genotype data, allowing better identification of causal variants in disease studies. However, due to overrepresentation of Caucasian studies, there's a lack of understanding of genetic basis of health-outcomes in other ethnic populations. Therefore, facilitating imputation of missing key-predictor-variants that can potentially improve a risk health-outcome prediction model, specifically for Asian ancestry, is of utmost relevance. METHODS: We aimed to construct an imputation and analysis web-platform, that primarily facilitates, but is not limited to genotype imputation on East-Asians. The goal is to provide a collaborative imputation platform for researchers in the public domain towards rapidly and efficiently conducting accurate genotype imputation. RESULTS: We present an online genotype imputation platform, Multi-ethnic Imputation System (MI-System) (https://misystem.cgm.ntu.edu.tw/), that offers users 3 established pipelines, SHAPEIT2-IMPUTE2, SHAPEIT4-IMPUTE5, and Beagle5.1 for conducting imputation analyses. In addition to 1000 Genomes and Hapmap3, a new customized Taiwan Biobank (TWB) reference panel, specifically created for Taiwanese-Chinese ancestry is provided. MI-System further offers functions to create customized reference panels to be used for imputation, conduct quality control, split whole genome data into chromosomes, and convert genome builds. CONCLUSION: Users can upload their genotype data and perform imputation with minimum effort and resources. The utility functions further can be utilized to preprocess user uploaded data with easy clicks. MI-System potentially contributes to Asian-population genetics research, while eliminating the requirement for high performing computational resources and bioinformatics expertise. It will enable an increased pace of research and provide a knowledge-base for genetic carriers of complex diseases, therefore greatly enhancing patient-driven research. STATEMENT OF SIGNIFICANCE: Multi-ethnic Imputation System (MI-System), primarily facilitates, but is not limited to, imputation on East-Asians, through 3 established prephasing-imputation pipelines, SHAPEIT2-IMPUTE2, SHAPEIT4-IMPUTE5, and Beagle5.1, where users can upload their genotype data and perform imputation and other utility functions with minimum effort and resources. A new customized Taiwan Biobank (TWB) reference panel, specifically created for Taiwanese-Chinese ancestry is provided. Utility functions include (a) create customized reference panels, (b) conduct quality control, (c) split whole genome data into chromosomes, and (d) convert genome builds. Users can also combine 2 reference panels using the system and use combined panels as reference to conduct imputation using MI-System.

GJB2 p.V37I Mutation Associated With Moderate Nonsyndromic Hearing Loss in an Adult Taiwanese Population.

BACKGROUND: Gap junction protein beta 2 ( GJB2 ) p.V37I mutations are the most important hereditary cause of sensorineural hearing loss (SNHL) in Taiwan. Hearing outcomes are associated with hearing levels at baseline and the duration of follow-up. However, the audiological features of GJB2 p.V37I mutations in the adult population are unknown. The objectives of the present study were to investigate the audiological features, progression rate, and allele frequency of GJB2 p.V37I mutations among an adult Taiwanese population. METHODS: Subjects of this case-control study were chosen from 13,580 participants of the Taiwan Precision Medicine Initiative. The genetic variations of GJB2 p.V37I were determined by polymerase chain reaction. We analyzed existing pure-tone threshold data from 38 individuals who were homozygous or compound heterozygotes for GJB2 p.V37I, 129 who were heterozygotes, and 602 individuals who were wild-type. Phenome-wide association studies (PheWAS) analysis was also performed to identify phenotypes associated with GJB2 p.V37I. RESULTS: The minor allele frequency of GJB2 p.V37I was 0.92% in our study population. The mean hearing level of participants with a p.V37I mutation indicated moderate to severe hearing loss with 38.2% ± 22.3% binaural hearing impairment. GJB2 p.V37I was associated with an increased risk of hearing disability (odds ratio: 21.46, 95% confidence interval: 8.62 to 53.44, p < 0.001) in an autosomal recessive pattern. In addition, PheWAS discovered a significant association between GJB2 p.V37I and fracture of the humerus. GJB2 p.V37I is a pathogenic and prevalent variant of SNHL among the adult population. CONCLUSIONS: The present study recommends patients with known GJB2 p.V37I mutations receive regular audiometric evaluation and genetic counseling. Early assistive listening device intervention is suggested to improve the quality of hearing.

Association test using Copy Number Profile Curves (CONCUR) enhances power in rare copy number variant analysis.

Copy number variants (CNVs) are the gain or loss of DNA segments in the genome that can vary in dosage and length. CNVs comprise a large proportion of variation in human genomes and impact health conditions. To detect rare CNV associations, kernel-based methods have been shown to be a powerful tool due to their flexibility in modeling the aggregate CNV effects, their ability to capture effects from different CNV features, and their accommodation of effect heterogeneity. To perform a kernel association test, a CNV locus needs to be defined so that locus-specific effects can be retained during aggregation. However, CNV loci are arbitrarily defined and different locus definitions can lead to different performance depending on the underlying effect patterns. In this work, we develop a new kernel-based test called CONCUR (i.e., copy number profile curve-based association test) that is free from a definition of locus and evaluates CNV-phenotype associations by comparing individuals' copy number profiles across the genomic regions. CONCUR is built on the proposed concepts of "copy number profile curves" to describe the CNV profile of an individual, and the "common area under the curve (cAUC) kernel" to model the multi-feature CNV effects. The proposed method captures the effects of CNV dosage and length, accounts for the numerical nature of copy numbers, and accommodates between- and within-locus etiological heterogeneity without the need to define artificial CNV loci as required in current kernel methods. In a variety of simulation settings, CONCUR shows comparable or improved power over existing approaches. Real data analyses suggest that CONCUR is well powered to detect CNV effects in the Swedish Schizophrenia Study and the Taiwan Biobank.

Integrating resistance functions to predict response to induction chemotherapy in de novo acute myeloid leukemia.

OBJECTIVES: This study explored resistance functions and their interactions in de novo AML treated with the "7 + 3" induction regimen. METHODS: We analyzed RNA-sequencing profiles of whole bone marrow samples from 52 de novo AML patients who completed the "7 + 3" regimen and stratified patients into CR (n = 35) and non-CR (n = 17) groups. RESULTS: A systematic gene set analysis revealed significant associations between chemoresistance and mTOR (P < .001), myc (P < .001), mitochondrial oxidative phosphorylation (P < .001), and stemness (P = .002). These functions were independent with regard to gene contents and activity scores. An integration of these four functions showed a prediction of chemoresistance (area under the receiver operating characteristic curve = 0.815) superior to that of each function alone. Moreover, our proposed seven-gene scoring system significantly correlated with the four-function model (r = .97; P < .001) to predict chemoresistance to the "7 + 3" regimen. On multivariate analysis, a seven-gene score of ≥-0.027 (hazard ratio: 11.18; 95% confidence interval: 2.06-60.65; P = .005) was an independent risk factor for induction failure. CONCLUSIONS: Myc, OXPHOS, mTOR, and stemness were responsive for chemoresistance in AML. Treatments other than the "7 + 3" regimen need to be considered for de novo AML patients predicted to be refractory to the "7 + 3" regimen.

Differential correlation analysis of glioblastoma reveals immune ceRNA interactions predictive of patient survival.

BACKGROUND: Recent studies illuminated a novel role of microRNA (miRNA) in the competing endogenous RNA (ceRNA) interaction: two genes (ceRNAs) can achieve coexpression by competing for a pool of common targeting miRNAs. Individual biological investigations implied ceRNA interaction performs crucial oncogenic/tumor suppressive functions in glioblastoma multiforme (GBM). Yet, a systematic analysis has not been conducted to explore the functional landscape and prognostic significance of ceRNA interaction. RESULTS: Incorporating the knowledge that ceRNA interaction is highly condition-specific and modulated by the expressional abundance of miRNAs, we devised a ceRNA inference by differential correlation analysis to identify the miRNA-modulated ceRNA pairs. Analyzing sample-paired miRNA and gene expression profiles of GBM, our data showed that this alternative layer of gene interaction is essential in global information flow. Functional annotation analysis revealed its involvement in activated processes in brain, such as synaptic transmission, as well as critical tumor-associated functions. Notably, a systematic survival analysis suggested the strength of ceRNA-ceRNA interactions, rather than expressional abundance of individual ceRNAs, among three immune response genes (CCL22, IL2RB, and IRF4) is predictive of patient survival. The prognostic value was validated in two independent cohorts. CONCLUSIONS: This work addresses the lack of a comprehensive exploration into the functional and prognostic relevance of ceRNA interaction in GBM. The proposed efficient and reliable method revealed its significance in GBM-related functions and prognosis. The highlighted roles of ceRNA interaction provide a basis for further biological and clinical investigations.

Genome-wide identification of key modulators of gene-gene interaction networks in breast cancer.

BACKGROUND: With the advances in high-throughput gene profiling technologies, a large volume of gene interaction maps has been constructed. A higher-level layer of gene-gene interaction, namely modulate gene interaction, is composed of gene pairs of which interaction strengths are modulated by (i.e., dependent on) the expression level of a key modulator gene. Systematic investigations into the modulation by estrogen receptor (ER), the best-known modulator gene, have revealed the functional and prognostic significance in breast cancer. However, a genome-wide identification of key modulator genes that may further unveil the landscape of modulated gene interaction is still lacking. RESULTS: We proposed a systematic workflow to screen for key modulators based on genome-wide gene expression profiles. We designed four modularity parameters to measure the ability of a putative modulator to perturb gene interaction networks. Applying the method to a dataset of 286 breast tumors, we comprehensively characterized the modularity parameters and identified a total of 973 key modulator genes. The modularity of these modulators was verified in three independent breast cancer datasets. ESR1, the encoding gene of ER, appeared in the list, and abundant novel modulators were illuminated. For instance, a prognostic predictor of breast cancer, SFRP1, was found the second modulator. Functional annotation analysis of the 973 modulators revealed involvements in ER-related cellular processes as well as immune- and tumor-associated functions. CONCLUSIONS: Here we present, as far as we know, the first comprehensive analysis of key modulator genes on a genome-wide scale. The validity of filtering parameters as well as the conservativity of modulators among cohorts were corroborated. Our data bring new insights into the modulated layer of gene-gene interaction and provide candidates for further biological investigations.

Small Molecule T315 Promotes Casitas B-Lineage Lymphoma-Dependent Degradation of Epidermal Growth Factor Receptor via Y1045 Autophosphorylation.

RATIONALE: Despite the fact that tyrosine kinase inhibitors (TKIs) have been found effective in treating patients harboring activating mutations of epidermal growth factor receptor (EGFR), an acquired secondary mutation, T790M, which lowers the affinity to TKIs, can lead to EGFR TKI resistance after this standard treatment. OBJECTIVES: To evaluate the effect of small molecule T315 on EGFR degradation and its therapeutic efficacy in vitro and in vivo. METHODS: Lung adenocarcinoma cells were treated with T315, and cell proliferation and apoptotic proportion were determined by the CellTiter 96 AQueous MTS (3-[4,5-dimethylthiazol-2-yl]-5-[3-carboxymethoxyphenyl]-2-[4-sulfophenyl]-2H-tetrazolium, inner salt) assay and flow cytometry. The effects of T315 on EGFR mRNA and protein levels, autophosphorylation, ubiquitination, and degradation were evaluated by real-time polymerase chain reaction and Western blot, respectively. Direct targeting of T315 to EGFR was confirmed by the in vitro kinase assay and mass spectrometry. Finally, the preclinical effect of T315 was validated in the murine xenograft model in combination with a second-generation TKI, afatinib. MEASUREMENTS AND MAIN RESULTS: We identified T315 as a novel, potent small molecule for suppressing cancer cell proliferation in vitro and in vivo. The therapeutic effect was verified after T315 was combined with a second-generation TKI, afatinib, compared with a single drug administration. We found a new mechanism of action, in that T315 appears to directly bind EGFR and triggers EGFR-Y1045 autophosphorylation, whereby its degradation is triggered through the ubiquitin-proteasome pathway. CONCLUSIONS: Our evidence suggests that T315 is a novel class of anticancer drug that is able to inhibit the growth of EGFR-TKI-resistant lung adenocarcinoma cells by inducing the degradation of EGFR.

Crystal structures of Staphylococcal SaeR reveal possible DNA-binding modes.

Two-component system SaeRS of Staphylococcus regulates virulence factor expression through phosphorylation of the DNA-binding regulator SaeR by the sensor histidine kinase SaeS. Here crystal structures of the DNA-binding domain (DBD) of SaeR from two Staphylococcal species Staphylococcus epidermidis and Staphylococcus aureus were determined and showed similar folds. Analyzing the DNA binding activity of three mutants of SeSaeR, we observed that Thr217 is important in binding to the phosphate group of DNA and Trp219 may interact with the base pairs. Additionally, the tandem arrangement of DBD may represent a possible way for SaeR oligomerization on DNA.

Parameter optimization for constructing competing endogenous RNA regulatory network in glioblastoma multiforme and other cancers.

BACKGROUND: In addition to direct targeting and repressing mRNAs, recent studies reported that microRNAs (miRNAs) can bridge up an alternative layer of post-transcriptional gene regulatory networks. The competing endogenous RNA (ceRNA) regulation depicts the scenario where pairs of genes (ceRNAs) sharing, fully or partially, common binding miRNAs (miRNA program) can establish coexpression through competition for a limited pool of the miRNA program. While the dynamics of ceRNA regulation among cellular conditions have been verified based on in silico and in vitro experiments, comprehensive investigation into the strength of ceRNA regulation in human datasets remains largely unexplored. Furthermore, pan-cancer analysis of ceRNA regulation, to our knowledge, has not been systematically investigated. RESULTS: In the present study we explored optimal conditions for ceRNA regulation, investigated functions governed by ceRNA regulation, and evaluated pan-cancer effects. We started by investigating how essential factors, such as the size of miRNA programs, the number of miRNA program binding sites, and expression levels of miRNA programs and ceRNAs affect the ceRNA regulation capacity in tumors derived from glioblastoma multiforme patients captured by The Cancer Genome Atlas (TCGA). We demonstrated that increased numbers of common targeting miRNAs as well as the abundance of binding sites enhance ceRNA regulation and strengthen coexpression of ceRNA pairs. Also, our investigation revealed that the strength of ceRNA regulation is dependent on expression levels of both miRNA programs and ceRNAs. Through functional annotation analysis, our results indicated that ceRNA regulation is highly associated with essential cellular functions and diseases including cancer. Furthermore, the highly intertwined ceRNA regulatory relationship enables constitutive and effective intra-function regulation of genes in diverse types of cancer. CONCLUSIONS: Using gene and microRNA expression datasets from TCGA, we successfully quantified the optimal conditions for ceRNA regulation, which hinge on four essential parameters of ceRNAs. Our analysis suggests optimized ceRNA regulation is related to disease pathways and essential cellular functions. Furthermore, although the strength of ceRNA regulation is dynamic among cancers, its governing functions are stably maintained. The findings of this report contribute to better understanding of ceRNA dynamics and its crucial roles in cancers.

Co-modulation analysis of gene regulation in breast cancer reveals complex interplay between ESR1 and ERBB2 genes.

BACKGROUND: Gene regulation is dynamic across cellular conditions and disease subtypes. From the aspect of regulation under modulation, regulation strength between a pair of genes can be modulated by (dependent on) expression abundance of another gene (modulator gene). Previous studies have demonstrated the involvement of genes modulated by single modulator genes in cancers, including breast cancer. However, analysis of multi-modulator co-modulation that can further delineate the landscape of complex gene regulation is, to our knowledge, unexplored previously. In the present study we aim to explore the joint effects of multiple modulator genes in modulating global gene regulation and dissect the biological functions in breast cancer. RESULTS: To carry out the analysis, we proposed the Covariability-based Multiple Regression (CoMRe) method. The method is mainly built on a multiple regression model that takes expression levels of multiple modulators as inputs and regulation strength between genes as output. Pairs of genes were divided into groups based on their co-modulation patterns. Analyzing gene expression profiles from 286 breast cancer patients, CoMRe investigated ten candidate modulator genes that interacted and jointly determined global gene regulation. Among the candidate modulators, ESR1, ERBB2, and ADAM12 were found modulating the most numbers of gene pairs. The largest group of gene pairs was composed of ones that were modulated by merely ESR1. Functional annotation revealed that the group was significantly related to tumorigenesis and estrogen signaling in breast cancer. ESR1-ERBB2 co-modulation was the largest group modulated by more than one modulators. Similarly, the group was functionally associated with hormone stimulus, suggesting that functions of the two modulators are performed, at least partially, through modulation. The findings were validated in majorities of patients (> 99%) of two independent breast cancer datasets. CONCLUSIONS: We have showed CoMRe is a robust method to discover critical modulators in gene regulatory networks, and it is capable of achieving reproducible and biologically meaningful results. Our data reveal that gene regulatory networks modulated by single modulator or co-modulated by multiple modulators play important roles in breast cancer. Findings of this report illuminate complex and dynamic gene regulation under modulation and its involvement in breast cancer.

A novel significance score for gene selection and ranking.

MOTIVATION: When identifying differentially expressed (DE) genes from high-throughput gene expression measurements, we would like to take both statistical significance (such as P-value) and biological relevance (such as fold change) into consideration. In gene set enrichment analysis (GSEA), a score that can combine fold change and P-value together is needed for better gene ranking. RESULTS: We defined a gene significance score π-value by combining expression fold change and statistical significance (P-value), and explored its statistical properties. When compared to various existing methods, π-value based approach is more robust in selecting DE genes, with the largest area under curve in its receiver operating characteristic curve. We applied π-value to GSEA and found it comparable to P-value and t-statistic based methods, with added protection against false discovery in certain situations. Finally, in a gene functional study of breast cancer profiles, we showed that using π-value helps elucidating otherwise overlooked important biological functions. AVAILABILITY: http://gccri.uthscsa.edu/Pi_Value_Supplementary.asp CONTACT: xy@ieee.org, cheny8@uthscsa.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

microRNA-449a functions as a tumor suppressor in neuroblastoma through inducing cell differentiation and cell cycle arrest.

microRNA-449a (miR-449a) has been identified to function as a tumor suppressor in several types of cancers. However, the role of miR-449a in neuroblastoma has not been intensively investigated. We recently found that the overexpression of miR-449a significantly induces neuroblastoma cell differentiation, suggesting its potential tumor suppressor function in neuroblastoma. In this study, we further investigated the mechanisms underlying the tumor suppressive function of miR-449a in neuroblastoma. We observed that miR-449a inhibits neuroblastoma cell survival and growth through 2 mechanisms--inducing cell differentiation and cell cycle arrest. Our comprehensive investigations on the dissection of the target genes of miR-449a revealed that 3 novel targets- MFAP4, PKP4 and TSEN15 -play important roles in mediating its differentiation-inducing function. In addition, we further found that its function in inducing cell cycle arrest involves down-regulating its direct targets CDK6 and LEF1. To determine the clinical significance of the miR-449a-mediated tumor suppressive mechanism, we examined the correlation between the expression of these 5 target genes in neuroblastoma tumor specimens and the survival of neuroblastoma patients. Remarkably, we noted that high tumor expression levels of all the 3 miR-449a target genes involved in regulating cell differentiation, but not the target genes involved in regulating cell cycle, are significantly correlated with poor survival of neuroblastoma patients. These results suggest the critical role of the differentiation-inducing function of miR-449a in determining neuroblastoma progression. Overall, our study provides the first comprehensive characterization of the tumor-suppressive function of miR-449a in neuroblastoma, and reveals the potential clinical significance of the miR-449a-mediated tumor suppressive pathway in neuroblastoma prognosis.

Predicting drug response through tumor deconvolution by cancer cell lines.

Large-scale cancer drug sensitivity data have become available for a collection of cancer cell lines, but only limited drug response data from patients are available. Bridging the gap in pharmacogenomics knowledge between in vitro and in vivo datasets remains challenging. In this study, we trained a deep learning model, Scaden-CA, for deconvoluting tumor data into proportions of cancer-type-specific cell lines. Then, we developed a drug response prediction method using the deconvoluted proportions and the drug sensitivity data from cell lines. The Scaden-CA model showed excellent performance in terms of concordance correlation coefficients (>0.9 for model testing) and the correctly deconvoluted rate (>70% across most cancers) for model validation using Cancer Cell Line Encyclopedia (CCLE) bulk RNA data. We applied the model to tumors in The Cancer Genome Atlas (TCGA) dataset and examined associations between predicted cell viability and mutation status or gene expression levels to understand underlying mechanisms of potential value for drug repurposing.

Polygenic risk score-based prediction of breast cancer risk in Taiwanese women with dense breast using a retrospective cohort study.

Mammographic screening has contributed to a significant reduction in breast cancer mortality. Several studies have highlighted the correlation between breast density, as detected through mammography, and a higher likelihood of developing breast cancer. A polygenic risk score (PRS) is a numerical score that is calculated based on an individual's genetic information. This study aims to explore the potential roles of PRS as candidate markers for breast cancer development and investigate the genetic profiles associated with clinical characteristics in Asian females with dense breasts. This is a retrospective cohort study integrated breast cancer screening, population genotyping, and cancer registry database. The PRSs of the study cohort were estimated using genotyping data of 77 single nucleotide polymorphisms based on the PGS000001 Catalog. A subgroup analysis was conducted for females without breast symptoms. Breast cancer patients constituted a higher proportion of individuals in PRS Q4 (37.8% vs. 24.8% in controls). Among dense breast patients with no symptoms, the high PRS group (Q4) consistently showed a significantly elevated breast cancer risk compared to the low PRS group (Q1-Q3) in both univariate (OR = 2.25, 95% CI 1.43-3.50, P < 0.001) and multivariate analyses (OR: 2.23; 95% CI 1.41-3.48, P < 0.001). The study was extended to predict breast cancer risk using common low-penetrance risk variants in a PRS model, which could be integrated into personalized screening strategies for Taiwanese females with dense breasts without prominent symptoms.

The Expression of TP63 as a Biomarker of Early Recurrence in Resected Esophageal Squamous Cell Carcinoma after Neoadjuvant Chemoradiotherapy.

Esophageal cancer ranks among the ten most common cancers worldwide. Despite the adoption of neoadjuvant concurrent chemoradiotherapy (nCCRT) followed by surgery as the standard treatment approach in recent years, the local recurrence rate remains high. In this study, we employed RNA-seq to investigate distinctive gene expression profiles in esophageal squamous cell carcinoma (ESCC) with or without recurrence following a standard treatment course. Our findings indicate that recurrent ESCC exhibits heightened keratinizing and epidermis development activity compared to non-recurrent ESCC. We identified TP63 as a potential candidate for distinguishing clinical outcomes. Furthermore, immunohistochemistry confirmed the trend of TP63 overexpression in ESCC recurrence. Patients with elevated TP63 expression had poorer overall survival and lower 3-year recurrence-free survival. This study underscores the potential of TP63 as a biomarker for detecting cancer recurrence and suggests its role in guiding future treatment options.

Machine learning approaches to identify systemic lupus erythematosus in anti-nuclear antibody-positive patients using genomic data and electronic health records.

BACKGROUND: Although the 2019 EULAR/ACR classification criteria for systemic lupus erythematosus (SLE) has required at least a positive anti-nuclear antibody (ANA) titer (≥ 1:80), it remains challenging for clinicians to identify patients with SLE. This study aimed to develop a machine learning (ML) approach to assist in the detection of SLE patients using genomic data and electronic health records. METHODS: Participants with a positive ANA (≥ 1:80) were enrolled from the Taiwan Precision Medicine Initiative cohort. The Taiwan Biobank version 2 array was used to detect single nucleotide polymorphism (SNP) data. Six ML models, Logistic Regression, Random Forest (RF), Support Vector Machine, Light Gradient Boosting Machine, Gradient Tree Boosting, and Extreme Gradient Boosting (XGB), were used to identify SLE patients. The importance of the clinical and genetic features was determined by Shapley Additive Explanation (SHAP) values. A logistic regression model was applied to identify genetic variations associated with SLE in the subset of patients with an ANA equal to or exceeding 1:640. RESULTS: A total of 946 SLE and 1,892 non-SLE controls were included in this analysis. Among the six ML models, RF and XGB demonstrated superior performance in the differentiation of SLE from non-SLE. The leading features in the SHAP diagram were anti-double strand DNA antibodies, ANA titers, AC4 ANA pattern, polygenic risk scores, complement levels, and SNPs. Additionally, in the subgroup with a high ANA titer (≥ 1:640), six SNPs positively associated with SLE and five SNPs negatively correlated with SLE were discovered. CONCLUSIONS: ML approaches offer the potential to assist in diagnosing SLE and uncovering novel SNPs in a group of patients with autoimmunity.

Polygenic Risk Score Predicts Earlier-Onset Adult Systemic Lupus Erythematosus and First-Year Renal Diseases in a Taiwanese Cohort.

OBJECTIVES: This study aimed to develop a predictive model using polygenic risk score (PRS) to forecast renal outcomes for adult systemic lupus erythematosus (SLE) in a Taiwanese population. METHODS: Patients with SLE (n=2782) and matched non-SLE controls (n=11 128) were genotyped using Genome-Wide TWB 2.0 single-nucleotide polymorphism (SNP) array. PRS models (C+T, LDpred2, Lassosum, PRSice-2, PRS-continuous shrinkage (CS)) were constructed for predicting SLE susceptibility. Logistic regression was assessed for C+T-based PRS association with renal involvement in patients with SLE. RESULTS: In the training set, C+T-based SLE-PRS, only incorporating 27 SNPs, outperformed other models with area under the curve (AUC) values of 0.629, surpassing Lassosum (AUC=0.621), PRSice-2 (AUC=0.615), LDpred2 (AUC=0.609) and PRS-CS (AUC=0.602). Additionally, C+T-based SLE-PRS demonstrated consistent predictive capacity in the testing set (AUC=0.620). Individuals in the highest quartile exhibited earlier SLE onset (39.06 vs 44.22 years, p<0.01), higher Systemic Lupus Erythematosus Disease Activity Index scores (3.00 vs 2.37, p=0.04), elevated risks of renal involvement within the first year of SLE diagnosis, including WHO class III-IV lupus nephritis (OR 2.36, 95% CI 1.47 to 3.80, p<0.01), estimated glomerular filtration rate <60 mL/min/1.73m2 (OR 1.49, 95% CI 1.18 to 1.89, p<0.01) and urine protein-to-creatinine ratio >150 mg/day (OR 2.07, 95% CI 1.49 to 2.89, p<0.01), along with increased seropositivity risks, compared with those in the lowest quartile. Furthermore, among patients with SLE with onset before 50 years, the highest PRS quartile was significantly associated with more serious renal diseases within the first year of SLE diagnosis. CONCLUSIONS: PRS of SLE is associated with earlier onset, renal involvement within the first year of SLE diagnosis and seropositivity in Taiwanese patients. Integrating PRS with clinical decision-making may enhance lupus nephritis screening and early treatment to improve renal outcomes in patients with SLE.