Clinical and Molecular Characteristics of PRKACA L206R Mutant Cortisol-Producing Adenomas in Korean Patients.

BACKGROUND: An activating mutation (c.617A>C/p.Lys206Arg, L206R) in protein kinase cAMP-activated catalytic subunit alpha (PRKACA) has been reported in 35% to 65% of cases of cortisol-producing adenomas (CPAs). We aimed to compare the clinical characteristics and transcriptome analysis between PRKACA L206R mutants and wild-type CPAs in Korea. METHODS: We included 57 subjects with CPAs who underwent adrenalectomy at Seoul National University Hospital. Sanger sequencing for PRKACA was conducted in 57 CPA tumor tissues. RNA sequencing was performed in 13 fresh-frozen tumor tissues. RESULTS: The prevalence of the PRKACA L206R mutation was 51% (29/57). The mean age of the study subjects was 42±12 years, and 87.7% (50/57) of the patients were female. Subjects with PRKACA L206R mutant CPAs showed smaller adenoma size (3.3±0.7 cm vs. 3.8±1.2 cm, P=0.059) and lower dehydroepiandrosterone sulfate levels (218±180 ng/mL vs. 1,511±3,307 ng/mL, P=0.001) than those with PRKACA wild-type CPAs. Transcriptome profiling identified 244 differentially expressed genes (DEGs) between PRKACA L206R mutant (n=8) and wild-type CPAs (n=5), including five upregulated and 239 downregulated genes in PRKACA L206R mutant CPAs (|fold change| ≥2, P<0.05). Among the upstream regulators of DEGs, CTNNB1 was the most significant transcription regulator. In several pathway analyses, the Wnt signaling pathway was downregulated and the steroid biosynthesis pathway was upregulated in PRKACA mutants. Protein-protein interaction analysis also showed that PRKACA downregulates Wnt signaling and upregulates steroid biosynthesis. CONCLUSION: The PRKACA L206R mutation in CPAs causes high hormonal activity with a limited proliferative capacity, as supported by transcriptome profiling.

Genetic profiles of subcutaneous panniculitis-like T-cell lymphoma and clinicopathological impact of HAVCR2 mutations.

Recent studies identified germline mutations in HAVCR2 (encoding T-cell immunoglobulin mucin 3) as a genetic factor that predisposes to subcutaneous panniculitis-like T-cell lymphoma (SPTCL). However, the differences between HAVCR2-mutated (HAVCR2MUT) and HAVCR2 wild-type (HAVCR2WT) SPTCLs remain unclear. A nationwide cohort of 53 patients with SPTCL diagnosed at 8 Korean institutions was established. Whole-exome sequencing and RNA-sequencing were performed on 8 patients in the discovery set. In the validation set, targeted gene sequencing or direct sequencing of HAVCR2 was performed. Of 49 patients with available HAVCR2 status, 25 (51.0%) were HAVCR2Y82C. HAVCR2Y82C was associated with younger age (P = .001), development of hemophagocytic lymphohistiocytosis or hemophagocytic lymphohistiocytosis-like systemic illness (P < .001), and short relapse-free survival (RFS) (P = .023). Most mutated genes in SPTCLs were involved in immune responses, epigenetic modifications, and cell signaling. Mutations in UNC13D, PIAS3, and KMT2D were more frequent in HAVCR2WT SPTCLs. At the gene expression level, HAVCR2Y82C SPTCLs were enriched in genes involved in IL6-JAK-STAT3 signaling and in tumor necrosis factor-α signaling via NF-κB. CCR4 was significantly upregulated in HAVCR2WT SPTCLs both at the messenger RNA level and at the protein level. We established a risk stratification system for SPTCL by integrating clinical and histopathological features, including age and HAVCR2 mutation status. This risk stratification system was strongly associated with RFS (P = .031). In conclusion, the HAVCR2Y82C mutation was common in Korean patients with SPTCL and was associated with unique clinicopathological and genetic features. Combining clinicopathological parameters could aid in predicting prognosis for patients with SPTCL.

Whole-exome sequencing in 168 Korean patients with inherited retinal degeneration.

BACKGROUND: To date, no genetic analysis of inherited retinal disease (IRD) using whole-exome sequencing (WES) has been conducted in a large-scale Korean cohort. The aim of this study was to characterise the genetic profile of IRD patients in Korea using WES. METHODS: We performed comprehensive molecular testing in 168 unrelated Korean IRD patients using WES. The potential pathogenicity of candidate variants was assessed using the American College of Medical Genetics and Genomics and the Association for Molecular Pathology variant interpretation guidelines, in silico prediction tools, published literature, and compatibility with known phenotypes or inheritance patterns. RESULTS: Causative variants were detected in 86/168 (51.2%) IRD patients, including 58/107 (54.2%) with retinitis pigmentosa, 7/15 (46.7%) with cone and cone-rod dystrophy, 2/3 (66.6%) with Usher syndrome, 1/2 (50.0%) with congenital stationary night blindness, 2/2 (100.0%) with Leber congenital amaurosis, 1/1 (100.0%) with Bietti crystalline dystrophy, 1/1 (100.0%) with Joubert syndrome, 9/10 (90.0%) with Stargardt macular dystrophy, 1/10 (10.0%) with vitelliform macular dystrophy, 1/11 (9.1%) with other forms of macular dystrophy, and 3/4 (75.0%) with choroideraemia. USH2A, ABCA4, and EYS were the most common causative genes associated with IRD. For retinitis pigmentosa, variants of USH2A and EYS were the most common causative gene mutations. CONCLUSIONS: This study demonstrated the distribution of causative genetic mutations in Korean IRD patients. The data will serve as a reference for future genetic screening and development of treatment modalities for Korean IRD patients.

Comprehensive Gene Expression Analyses of Immunohistochemically Defined Subgroups of Muscle-Invasive Urinary Bladder Urothelial Carcinoma.

A number of urinary bladder urothelial carcinoma (UB UC) mRNA-based classification systems have been reported. It also has been observed that treatment response and prognosis are different for each molecular subtype. In this study, cytokeratin (CK)5/6 and CK20 immunohistochemistry (IHC) were performed, and IHC-based subgroup classification was applied. UB UC was classified into CK5/6 single-positive (SP), CK20 SP, double-positive (DP) and double-negative (DN) subgroups, and transcriptional analysis was performed. The results of gene ontology (GO) terms and functional analysis using differentially expressed genes indicate that, CK5/6 SP and DP subgroups were enriched in cell migration, immune activation, interleukin 6-Janus kinase-signal transducer and activator of transcription 3 (IL6-JAK-STAT3) signaling pathway and tumor necrosis factor-α signaling via the nuclear factor-κB (NF-κB) signaling pathway signature gene. In addition, compared with the other subgroups, the DN subgroup showed inhibited cell movement, cell migration, and cell activation. Furthermore, in survival analysis, the CK5/6 SP subgroup was significantly associated with poor progression-free survival (p = 0.008). The results of our study indicate that the CK5/6 positive subgroup exhibited high gene expression signature related to aggressive behavior and exhibited worse clinical outcome.

Non-Muscle-Invasive Bladder Carcinoma with Respect to Basal Versus Luminal Keratin Expression.

Non-muscle-invasive bladder cancer (NMIBC) consists of transcriptional subtypes that are distinguishable from those of muscle-invasive cancer. We aimed to identify genetic signatures of NMIBC related to basal (K5/6) and luminal (K20) keratin expression. Based on immunohistochemical staining, papillary high-grade NMIBC was classified into K5/6-only (K5/6High-K20Low), K20-only (K5/6Low-K20High), double-high (K5/6High-K20High), and double-low (K5/6Low-K20Low) groups (n = 4 per group). Differentially expressed genes identified between each group using RNA sequencing were subjected to functional enrichment analyses. A public dataset was used for validation. Machine learning algorithms were implemented to predict our samples against UROMOL subtypes. Transcriptional investigation demonstrated that the K20-only group was enriched in the cell cycle, proliferation, and progression gene sets, and this result was also observed in the public dataset. The K5/6-only group was closely regulated by basal-type gene sets and showed activated invasive or adhesive functions. The double-high group was enriched in cell cycle arrest, macromolecule biosynthesis, and FGFR3 signaling. The double-low group moderately expressed genes related to cell cycle and macromolecule biosynthesis. All K20-only group tumors were classified as UROMOL "class 2" by the machine learning algorithms. K5/6 and K20 expression levels indicate the transcriptional subtypes of NMIBC. The K5/6Low-K20High expression is a marker of high-risk NMIBC.

Discovery of Novel Recurrent Mutations and Clinically Meaningful Subgroups in Nodal Marginal Zone Lymphoma.

Nodal marginal zone lymphoma (NMZL) is a rare B-cell neoplasm, the genetic and transcriptomic landscape of which are unclear. Using high-throughput sequencing for whole-exome and transcriptome, we investigated the genetic characteristics of NMZL in a discovery cohort (n = 8) and validated their features in an extended cohort (n = 30). Novel mutations in NFKBIE and ITPR2 were found in 7.9% (3/38) and 13.9% (5/36), respectively, suggesting roles for the NF-κB pathway and B-cell-receptor-mediated calcium signaling pathway in the pathogenesis of NMZL. RNA-seq showed that NMZLs were characterized by an aberrant marginal zone differentiation, associated with an altered IRF4-NOTCH2 axis and the enrichment of various oncogenic pathways. Based on gene expression profile, two subgroups were identified. Compared with subgroup 1, subgroup 2 showed the following: the significant enrichment of cell cycle-associated and MYC-signaling pathways, a more diverse repertoire of upstream regulators, and higher Ki-67 proliferation indices. We designated two subgroups according to Ki-67 labeling, and subgroup 2 was significantly associated with a shorter progression-free survival (p = 0.014), a greater proportion of large cells (p = 0.009), and higher MYC expression (p = 0.026). We suggest that NMZL has unique features and, in this study, we provide information as to the heterogeneity of this enigmatic entity.

CK14 Expression Identifies a Basal/Squamous-Like Type of Papillary Non-Muscle-Invasive Upper Tract Urothelial Carcinoma.

Object: CK14 expression is an important marker of basal/squamous-like (BASQ)-type muscle-invasive bladder carcinoma, and this molecularly defined subtype has a poor prognosis and a distinct response to chemotherapy. However, CK14 expression and its clinicopathological and molecular significance in papillary non-muscle-invasive upper tract urothelial carcinoma (NMIUTUC) remain unknown. Herein, we investigated the prognostic implications of immunohistochemical (IHC) staining for CK14 and the transcriptional characteristics associated with CK14 expression in papillary NMIUTUC. Materials and Methods: IHC staining for CK14 was conducted in 204 papillary NMIUTUC specimens. Positive CK14 IHC staining was defined as a positive signal in >0% of tumor cells. RNA sequencing data were analyzed from 8 papillary high-grade NMIUTUC specimens consisting of 4 CK14-positive and 4 CK14-negative tumors. Results: CK14 positivity was associated with a high TNM stage (p < 0.001) and a high World Health Organization grade (p = 0.003). Survival analysis showed that CK14 positivity was significantly associated with poor progression-free survival (p = 0.015; hazard ratio [HR] = 2.990; 95% confidence interval [CI] = 1.180-7.580) and was marginally associated with poor cancer-specific survival (p = 0.052; HR = 3.77; 95% CI = 0.900-15.780). Gene set enrichment analysis demonstrated that the CK14-positive tumors were associated with a basal subtype of breast cancer, squamous cell carcinoma development, p40, tumor necrosis factor α-nuclear factor-κB, and p53 pathways, and embryonic stem cells; these characteristics are reminiscent of the BASQ subtype. In addition, with a p < 0.05 and |fold change| ≥2 as the cutoffs, we identified 178 differentially expressed genes when comparing CK14-positive and CK14-negative tumors. Functional analysis of these genes revealed several networks and gene ontology terms related to the positive regulation of cellular proliferation in CK14-positive tumors. Consistent with these results, we demonstrated that the mean Ki-67 proliferative index was higher in CK14-positive tumors than it was in CK14-negative tumors (2.3 vs. 0.8%, respectively, p = 0.002). Conclusion: CK14-positive papillary NMIUTUC is an aggressive subtype with BASQ-like molecular characteristics and dynamic proliferative activity. We propose that CK14 IHC staining can be a useful biomarker of BASQ-type papillary NMIUTUC that can be applied in daily practice with the aim of precision oncology.

Bioinformatics services for analyzing massive genomic datasets.

The explosive growth of next-generation sequencing data has resulted in ultra-large-scale datasets and ensuing computational problems. In Korea, the amount of genomic data has been increasing rapidly in the recent years. Leveraging these big data requires researchers to use large-scale computational resources and analysis pipelines. A promising solution for addressing this computational challenge is cloud computing, where CPUs, memory, storage, and programs are accessible in the form of virtual machines. Here, we present a cloud computing-based system, Bio-Express, that provides user-friendly, cost-effective analysis of massive genomic datasets. Bio-Express is loaded with predefined multi-omics data analysis pipelines, which are divided into genome, transcriptome, epigenome, and metagenome pipelines. Users can employ predefined pipelines or create a new pipeline for analyzing their own omics data. We also developed several web-based services for facilitating downstream analysis of genome data. Bio-Express web service is freely available at https://www.bioexpress.re.kr/.

Comparative analysis of commonly used peak calling programs for ChIP-Seq analysis.

Chromatin immunoprecipitation coupled with high-throughput DNA sequencing (ChIP-Seq) is a powerful technology to profile the location of proteins of interest on a whole-genome scale. To identify the enrichment location of proteins, many programs and algorithms have been proposed. However, none of the commonly used peak calling programs could accurately explain the binding features of target proteins detected by ChIP-Seq. Here, publicly available data on 12 histone modifications, including H3K4ac/me1/me2/me3, H3K9ac/me3, H3K27ac/me3, H3K36me3, H3K56ac, and H3K79me1/me2, generated from a human embryonic stem cell line (H1), were profiled with five peak callers (CisGenome, MACS1, MACS2, PeakSeq, and SISSRs). The performance of the peak calling programs was compared in terms of reproducibility between replicates, examination of enriched regions to variable sequencing depths, the specificity-to-noise signal, and sensitivity of peak prediction. There were no major differences among peak callers when analyzing point source histone modifications. The peak calling results from histone modifications with low fidelity, such as H3K4ac, H3K56ac, and H3K79me1/me2, showed low performance in all parameters, which indicates that their peak positions might not be located accurately. Our comparative results could provide a helpful guide to choose a suitable peak calling program for specific histone modifications.

Transcription-dependent targeting of Hda1C to hyperactive genes mediates H4-specific deacetylation in yeast.

In yeast, Hda1 histone deacetylase complex (Hda1C) preferentially deacetylates histones H3 and H2B, and functionally interacts with Tup1 to repress transcription. However, previous studies identified global increases in histone H4 acetylation in cells lacking Hda1, a component of Hda1C. Here, we find that Hda1C binds to hyperactive genes, likely via the interaction between the Arb2 domain of Hda1 and RNA polymerase II. Additionally, we report that Hda1C specifically deacetylates H4, but not H3, at hyperactive genes to partially inhibit elongation. This role is contrast to that of the Set2-Rpd3S pathway deacetylating histones at infrequently transcribed genes. We also find that Hda1C deacetylates H3 at inactive genes to delay the kinetics of gene induction. Therefore, in addition to fine-tuning of transcriptional response via H3-specific deacetylation, Hda1C may modulate elongation by specifically deacetylating H4 at highly transcribed regions.