An Organoid Biobank of Neuroendocrine Neoplasms Enables Genotype-Phenotype Mapping.

Gastroenteropancreatic (GEP) neuroendocrine neoplasm (NEN) that consists of neuroendocrine tumor and neuroendocrine carcinoma (NEC) is a lethal but under-investigated disease owing to its rarity. To fill the scarcity of clinically relevant models of GEP-NEN, we here established 25 lines of NEN organoids and performed their comprehensive molecular characterization. GEP-NEN organoids recapitulated pathohistological and functional phenotypes of the original tumors. Whole-genome sequencing revealed frequent genetic alterations in TP53 and RB1 in GEP-NECs, and characteristic chromosome-wide loss of heterozygosity in GEP-NENs. Transcriptome analysis identified molecular subtypes that are distinguished by the expression of distinct transcription factors. GEP-NEN organoids gained independence from the stem cell niche irrespective of genetic mutations. Compound knockout of TP53 and RB1, together with overexpression of key transcription factors, conferred on the normal colonic epithelium phenotypes that are compatible with GEP-NEN biology. Altogether, our study not only provides genetic understanding of GEP-NEN, but also connects its genetics and biological phenotypes.

Genomic features of renal cell carcinoma developed during end-stage renal disease and dialysis.

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Variant spectrum of von Hippel-Lindau disease and its genomic heterogeneity in Japan.

Von Hippel-Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80-90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5-22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

Postoperative recurrence detection using individualized circulating tumor DNA in upper tract urothelial carcinoma.

Biomarkers that could detect the postoperative recurrence of upper tract urothelial carcinoma (UTUC) have not been established. In this prospective study, we aim to evaluate the utility of individualized circulating tumor DNA (ctDNA) monitoring using digital PCR (dPCR) as a tumor recurrence biomarker for UTUC in the perioperative period. Twenty-three patients who underwent radical nephroureterectomy (RNU) were included. In each patient, whole exome sequencing by next-generation sequencing and TERT promoter sequencing of tumor DNA were carried out. Case-specific gene mutations were selected from sequencing analysis to examine ctDNA by dPCR analysis. We also prospectively collected plasma and urine ctDNA from each patient. The longitudinal variant allele frequencies of ctDNA during the perioperative period were plotted. Case-specific gene mutations were detected in 22 cases (96%) from ctDNA in the preoperative samples. Frequently detected genes were TERT (39%), FGFR3 (26%), TP53 (22%), and HRAS (13%). In all cases, we obtained plasma and urine samples for 241 time points and undertook individualized ctDNA monitoring for 2 years after RNU. Ten patients with intravesical recurrence had case-specific ctDNA detected in urine at the time of recurrence. The mean lead time of urinary ctDNA in intravesical recurrence was 60 days (range, 0-202 days). Two patients with distal metastasis had case-specific ctDNA in plasma at the time of metastasis. In UTUC, tumor-specific gene mutations can be monitored postoperatively as ctDNA in plasma and urine. Individualized ctDNA might be a minimally invasive biomarker for the early detection of postoperative recurrence.

Hereditary cancer variants and homologous recombination deficiency in biliary tract cancer.

BACKGROUND & AIMS: The heritability and actionability of variants in homologous recombination-related genes in biliary tract cancers (BTCs) are uncertain. Although associations between BTC and BRCA germline variants have been reported, homologous recombination deficiency has not been investigated in BTCs. METHODS: We sequenced germline variants in 27 cancer-predisposing genes in 1,292 BTC cases and 37,583 controls without a personal nor family history of cancer. We compared pathogenic germline variant frequencies between cases and controls and documented the demographic and clinical characteristics of carriers. In addition, whole-genome sequencing of 45 BTC tissues was performed to evaluate homologous recombination deficiency status. RESULTS: Targeted sequencing identified 5,018 germline variants, which were classified into 317 pathogenic, 3,611 variants of uncertain significance, and 1,090 benign variants. Seventy-one BTC cases (5.5%) had at least one pathogenic variant among 27 cancer-predisposing genes. Pathogenic germline variants enriched in BTCs were present in BRCA1, BRCA2, APC, and MSH6 (p <0.00185). PALB2 variants were marginally associated with BTC (p = 0.01). APC variants were predominantly found in ampulla of Vater carcinomas. Whole-genome sequencing demonstrated that three BTCs with pathogenic germline variants in BRCA2 and PALB2, accompanied by loss of heterozygosity, displayed homologous recombination deficiency. Conversely, pathogenic germline variants without a second hit or variants of other homologous recombination-related genes such as ATM and BRIP1 showed homologous recombination-proficient phenotypes. CONCLUSIONS: In this study, we describe the heritability and actionability of variants in homologous recombination-related genes, which could be used to guide screening and therapeutic strategies for BTCs. IMPACT AND IMPLICATIONS: We found that 5.5% of biliary tract cancers (BTCs) in a Japanese population possessed hereditary cancer-predisposing gene alterations, including in BRCA and genes associated with colorectal cancer. Two hits in homologous recombination-related genes were required to confer a homologous recombination-deficient phenotype. PARP inhibitors and DNA-damaging regimens may be effective strategies against BTCs exhibiting homologous recombination deficiency. Hence, in this study, genome-wide sequencing has revealed a potential new therapeutic strategy that could be applied to a subset of BTCs.

Mitochondrial DNA Variants at Low-Level Heteroplasmy and Decreased Copy Numbers in Chronic Kidney Disease (CKD) Tissues with Kidney Cancer.

The human mitochondrial genome (mtDNA) is a circular DNA molecule with a length of 16.6 kb, which contains a total of 37 genes. Somatic mtDNA mutations accumulate with age and environmental exposure, and some types of mtDNA variants may play a role in carcinogenesis. Recent studies observed mtDNA variants not only in kidney tumors but also in adjacent kidney tissues, and mtDNA dysfunction results in kidney injury, including chronic kidney disease (CKD). To investigate whether a relationship exists between heteroplasmic mtDNA variants and kidney function, we performed ultra-deep sequencing (30,000×) based on long-range PCR of DNA from 77 non-tumor kidney tissues of kidney cancer patients with CKD (stages G1 to G5). In total, this analysis detected 697 single-nucleotide variants (SNVs) and 504 indels as heteroplasmic (0.5% ≤ variant allele frequency (VAF) < 95%), and the total number of detected SNVs/indels did not differ between CKD stages. However, the number of deleterious low-level heteroplasmic variants (pathogenic missense, nonsense, frameshift and tRNA) significantly increased with CKD progression (p < 0.01). In addition, mtDNA copy numbers (mtDNA-CNs) decreased with CKD progression (p < 0.001). This study demonstrates that mtDNA damage, which affects mitochondrial genes, may be involved in reductions in mitochondrial mass and associated with CKD progression and kidney dysfunction.

Immuno-genomic analysis reveals eosinophilic feature and favorable prognosis of female non-smoking esophageal squamous cell carcinomas.

Most of esophageal squamous cell carcinoma (ESCC) develop in smoking males in Japan, but the genomic etiology and immunological characteristics of rare non-smoking female ECSS remain unclear. To elucidate the genomic and immunological features of ESCC in non-smoking females, we analyzed whole-genome or transcriptome sequencing data from 94 ESCCs, including 20 rare non-smoking female cases. In addition, 31,611 immune cells were extracted from four ESCC tissues and subject to single-cell RNA-seq. We compared their immuno-genomic and microbiome profiles between non-smoking female and smoking ESCCs. Non-smoking females showed much better prognosis. Whole-genome sequencing analysis showed no significant differences in driver genes or copy number alterations depending on smoking status. The mutational signatures specifically observed in non-smoking females ESCC could be attributed to aging. Immune profiling from RNA-seq revealed that ESCC in non-smoking females had high tumor microenvironment signatures and a high abundance of eosinophils with a favorable prognosis. Single-cell RNA-sequencing of intratumor immune cells revealed gender differences of eosinophils and their activation in female cases. ESCCs in non-smoking females have age-related mutational signatures and gender-specific tumor immune environment with eosinophils, which is likely to contribute to their favorable prognosis.

Comparative analyses define differences between BHD-associated renal tumour and sporadic chromophobe renal cell carcinoma.

BACKGROUND: Birt-Hogg-Dubé (BHD) syndrome, caused by germline alteration of folliculin (FLCN) gene, develops hybrid oncocytic/chromophobe tumour (HOCT) and chromophobe renal cell carcinoma (ChRCC), whereas sporadic ChRCC does not harbor FLCN alteration. To date, molecular characteristics of these similar histological types of tumours have been incompletely elucidated. METHODS: To elucidate renal tumourigenesis of BHD-associated renal tumours and sporadic renal tumours, we conducted whole genome sequencing (WGS) and RNA-sequencing (RNA-seq) of sixteen BHD-associated renal tumours from nine unrelated BHD patients, twenty-one sporadic ChRCCs and seven sporadic oncocytomas. We then compared somatic mutation profiles with FLCN variants and RNA expression profiles between BHD-associated renal tumours and sporadic renal tumours. FINDINGS: RNA-seq analysis revealed that BHD-associated renal tumours and sporadic renal tumours have totally different expression profiles. Sporadic ChRCCs were clustered into two distinct clusters characterized by L1CAM and FOXI1 expressions, molecular markers for renal tubule subclasses. Increased mitochondrial DNA (mtDNA) copy number with fewer variants was observed in BHD-associated renal tumours compared to sporadic ChRCCs. Cell-of-origin analysis using WGS data demonstrated that BHD-associated renal tumours and sporadic ChRCCs may arise from different cells of origin and second hit FLCN alterations may occur in early third decade of life in BHD patients. INTERPRETATION: These data further our understanding of renal tumourigenesis of these two different types of renal tumours with similar histology. FUNDING: This study was supported by JSPS KAKENHI Grants, RIKEN internal grant, and the Intramural Research Program of the National Institutes of Health (NIH), National Cancer Institute (NCI), Center for Cancer Research.

Increased PDGFRB and NF-κB signaling caused by highly prevalent somatic mutations in intracranial aneurysms.

Intracranial aneurysms (IAs) are a high-risk factor for life-threatening subarachnoid hemorrhage. Their etiology, however, remains mostly unknown at present. We conducted screening for sporadic somatic mutations in 65 IA tissues (54 saccular and 11 fusiform aneurysms) and paired blood samples by whole-exome and targeted deep sequencing. We identified sporadic mutations in multiple signaling genes and examined their impact on downstream signaling pathways and gene expression in vitro and an arterial dilatation model in mice in vivo. We identified 16 genes that were mutated in at least one IA case and found that these mutations were highly prevalent (92%: 60 of 65 IAs) among all IA cases examined. In particular, mutations in six genes (PDGFRB, AHNAK, OBSCN, RBM10, CACNA1E, and OR5P3), many of which are linked to NF-κB signaling, were found in both fusiform and saccular IAs at a high prevalence (43% of all IA cases examined). We found that mutant PDGFRBs constitutively activated ERK and NF-κB signaling, enhanced cell motility, and induced inflammation-related gene expression in vitro. Spatial transcriptomics also detected similar changes in vessels from patients with IA. Furthermore, virus-mediated overexpression of a mutant PDGFRB induced a fusiform-like dilatation of the basilar artery in mice, which was blocked by systemic administration of the tyrosine kinase inhibitor sunitinib. Collectively, this study reveals a high prevalence of somatic mutations in NF-κB signaling pathway-related genes in both fusiform and saccular IAs and opens a new avenue of research for developing pharmacological interventions.

Immuno-genomic profiling of biopsy specimens predicts neoadjuvant chemotherapy response in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive cancers and is primarily treated with platinum-based neoadjuvant chemotherapy (NAC). Some ESCCs respond well to NAC. However, biomarkers to predict NAC sensitivity and their response mechanism in ESCC remain unclear. We perform whole-genome sequencing and RNA sequencing analysis of 141 ESCC biopsy specimens before NAC treatment to generate a machine-learning-based diagnostic model to predict NAC reactivity in ESCC and analyzed the association between immunogenomic features and NAC response. Neutrophil infiltration may play an important role in ESCC response to NAC. We also demonstrate that specific copy-number alterations and copy-number signatures in the ESCC genome are significantly associated with NAC response. The interactions between the tumor genome and immune features of ESCC are likely to be a good indicator of therapeutic capability and a therapeutic target for ESCC, and machine learning prediction for NAC response is useful.

Proteo-genomic characterization of virus-associated liver cancers reveals potential subtypes and therapeutic targets.

Primary liver cancer is a heterogeneous disease in terms of its etiology, histology, and therapeutic response. Concurrent proteomic and genomic characterization of a large set of clinical liver cancer samples can help elucidate the molecular basis of heterogeneity and thus serve as a valuable resource for personalized liver cancer treatment. In this study, we perform proteomic profiling of ~300 proteins on 259 primary liver cancer tissues with reverse-phase protein arrays, mutational analysis using whole genome sequencing and transcriptional analysis with RNA-Seq. Patients are of Japanese ethnic background and mainly HBV or HCV positive, providing insight into this important liver cancer subtype. Unsupervised classification of tumors based on protein expression profiles reveal three proteomic subclasses R1, R2, and R3. The R1 subclass is immunologically hot and demonstrated a good prognosis. R2 contains advanced proliferative tumor with TP53 mutations, high expression of VEGF receptor 2 and the worst prognosis. R3 is enriched with CTNNB1 mutations and elevated mTOR signaling pathway activity. Twenty-two proteins, including CDK1 and CDKN2A, are identified as potential prognostic markers. The proteomic classification presented in this study can help guide therapeutic decision making for liver cancer treatment.

Actionability evaluation of biliary tract cancer by genome transcriptome analysis and Asian cancer knowledgebase.

INTRODUCTION: Treatment options for biliary tract cancer (BTC) are very limited. It is necessary to investigate actionable genes and candidate drugs using a sophisticated knowledgebase (KB) and characterize BTCs immunologically for evaluating the actionability of molecular and immune therapies. MATERIALS AND METHODS: The genomic and transcriptome data of 219 patients with BTC who underwent surgery were analyzed. Actionable mutations and candidate drugs were annotated using the largest available KB of the Asian population (CancerSCAN®). Predictive biomarkers of immune checkpoint inhibitors were analyzed using DNA and RNA sequencing data. RESULTS: Twenty-two actionable genes and 43 candidate drugs were annotated in 74 patients (33.8%). The most frequent actionable genes were PTEN (7.3%), CDKN2A (6.8%), KRAS (6.4%). BRCA2, CDKN2A, and FGFR2 mutations were most frequently identified in case of intrahepatic cholangiocarcinoma. PTEN and CDKN2A mutations were associated with significantly shorter overall survival. PD-L1 and PD-1 expression was significantly higher in case of extrahepatic cholangiocarcinoma and T-cell-high expression. In total, 49.7% of cases were evaluated as having actionability for molecular therapy or immune checkpoint inhibitors. CONCLUSIONS: Identifying actionable genes and candidate drugs using the KB contribute to the development of therapeutic drugs and personalized treatment for BTC.

Molecular Classification and Tumor Microenvironment Characterization of Gallbladder Cancer by Comprehensive Genomic and Transcriptomic Analysis.

Gallbladder cancer (GBC), a rare but lethal disease, is often diagnosed at advanced stages. So far, molecular characterization of GBC is insufficient, and a comprehensive molecular portrait is warranted to uncover new targets and classify GBC. We performed a transcriptome analysis of both coding and non-coding RNAs from 36 GBC fresh-frozen samples. The results were integrated with those of comprehensive mutation profiling based on whole-genome or exome sequencing. The clustering analysis of RNA-seq data facilitated the classification of GBCs into two subclasses, characterized by high or low expression levels of TME (tumor microenvironment) genes. A correlation was observed between gene expression and pathological immunostaining. TME-rich tumors showed significantly poor prognosis and higher recurrence rate than TME-poor tumors. TME-rich tumors showed overexpression of genes involved in epithelial-to-mesenchymal transition (EMT) and inflammation or immune suppression, which was validated by immunostaining. One non-coding RNA, miR125B1, exhibited elevated expression in stroma-rich tumors, and miR125B1 knockout in GBC cell lines decreased its invasion ability and altered the EMT pathway. Mutation profiles revealed TP53 (47%) as the most commonly mutated gene, followed by ELF3 (13%) and ARID1A (11%). Mutations of ARID1A, ERBB3, and the genes related to the TGF-ß signaling pathway were enriched in TME-rich tumors. This comprehensive analysis demonstrated that TME, EMT, and TGF-ß pathway alterations are the main drivers of GBC and provides a new classification of GBCs that may be useful for therapeutic decision-making.

Whole genome sequencing analysis identifies recurrent structural alterations in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) is the predominant type of esophageal cancer in the Asian region, including Japan. A previous study reported mutational landscape of Japanese ESCCs by using exome sequencing. However, somatic structural alterations were yet to be explored. To provide a comprehensive mutational landscape, we performed whole genome sequencing (WGS) analysis of biopsy specimens from 20 ESCC patients in a Japanese population. WGS analysis identified non-silent coding mutations of TP53, ZNF750 and FAT1 in ESCC. We detected six mutational signatures in ESCC, one of which showed significant association with smoking status. Recurrent structural variations, many of which were chromosomal deletions, affected genes such as LRP1B, TTC28, CSMD1, PDE4D, SDK1 and WWOX in 25%-30% of tumors. Somatic copy number amplifications at 11q13.3 (CCND1), 3q26.33 (TP63/SOX2), and 8p11.23 (FGFR1) and deletions at 9p21.3 (CDKN2A) were identified. Overall, these multi-dimensional view of genomic alterations improve the understanding of the ESCC development at molecular level and provides future prognosis and therapeutic implications for ESCC in Japan.

Genomic and Transcriptomic Profiling of Combined Hepatocellular and Intrahepatic Cholangiocarcinoma Reveals Distinct Molecular Subtypes.

We performed genomic and transcriptomic sequencing of 133 combined hepatocellular and intrahepatic cholangiocarcinoma (cHCC-ICC) cases, including separate, combined, and mixed subtypes. Integrative comparison of cHCC-ICC with hepatocellular carcinoma and intrahepatic cholangiocarcinoma revealed that combined and mixed type cHCC-ICCs are distinct subtypes with different clinical and molecular features. Integrating laser microdissection, cancer cell fraction analysis, and single nucleus sequencing, we revealed both mono- and multiclonal origins in the separate type cHCC-ICCs, whereas combined and mixed type cHCC-ICCs were all monoclonal origin. Notably, cHCC-ICCs showed significantly higher expression of Nestin, suggesting Nestin may serve as a biomarker for diagnosing cHCC-ICC. Our results provide important biological and clinical insights into cHCC-ICC.

Potential protective function of the sterol regulatory element binding factor 1-fatty acid desaturase 1/2 axis in early-stage age-related macular degeneration.

Age-related macular degeneration (AMD) is the most common cause of vision loss in elderly individuals throughout the developed world. Inhibitors of vascular endothelial growth factor have been successfully used to treat choroidal neovascularization in late-stage AMD. The pathogenesis of early-stage AMD, however, remains largely unknown, impairing efforts to develop effective therapies that prevent progression to late-stage AMD. To address this, we performed comparative transcriptomics of macular and extramacular retinal pigmented epithelium-choroid (RPE-choroid) tissue from early-stage AMD patients. We found that expression of fatty acid desaturase 1 (FADS1), FADS2, and acetyl-CoA acetyltransferase 2 (ACAT2) is increased in macular but not extramacular tissue, possibly through activation of sterol regulatory element binding factor 1 (SREBF1). Consistent with this, we also found that expression of Fads1 is increased in RPE-choroid in a mouse model of early-stage AMD. In zebrafish, deletion of fads2, which encodes a protein that functions as both Fads1 and Fads2 in other species, enhanced apoptosis in the retina upon exposure to intense light. Similarly, pharmacological inhibition of Srebf1 enhanced apoptosis and reduced fads2 expression in zebrafish exposed to intense light. These results suggest that the SREBF1-FADS1/2 axis may be activated in macular RPE-choroid as a protective response during early-stage AMD and could thus be a therapeutic target for early-stage AMD.

Activation of Sterol Regulatory Element Binding Factors by Fenofibrate and Gemfibrozil Stimulates Myelination in Zebrafish.

Oligodendrocytes are major myelin-producing cells and play essential roles in the function of a healthy nervous system. However, they are also one of the most vulnerable neural cell types in the central nervous system (CNS), and myelin abnormalities in the CNS are found in a wide variety of neurological disorders, including multiple sclerosis, adrenoleukodystrophy, and schizophrenia. There is an urgent need to identify small molecular weight compounds that can stimulate myelination. In this study, we performed comparative transcriptome analysis to identify pharmacodynamic effects common to miconazole and clobetasol, which have been shown to stimulate myelination by mouse oligodendrocyte progenitor cells (OPCs). Of the genes differentially expressed in both miconazole- and clobetasol-treated mouse OPCs compared with untreated cells, we identified differentially expressed genes (DEGs) common to both drug treatments. Gene ontology analysis revealed that these DEGs are significantly associated with the sterol biosynthetic pathway, and further bioinformatics analysis suggested that sterol regulatory element binding factors (SREBFs) might be key upstream regulators of the DEGs. In silico screening of a public database for chemicals associated with SREBF activation identified fenofibrate, a peroxisome proliferator-activated receptor α (PPARα) agonist, as a drug that increases the expression of known SREBF targets, raising the possibility that fenofibrate may also stimulate myelination. To test this, we performed in vivo imaging of zebrafish expressing a fluorescent reporter protein under the control of the myelin basic protein (mbp) promoter. Treatment of zebrafish with fenofibrate significantly increased expression of the fluorescent reporter compared with untreated zebrafish. This increase was attenuated by co-treatment with fatostatin, a specific inhibitor of SREBFs, confirming that the fenofibrate effect was mediated via SREBFs. Furthermore, incubation of zebrafish with another PPARα agonist, gemfibrozil, also increased expression of the mbp promoter-driven fluorescent reporter in an SREBF-dependent manner. These results suggest that activation of SREBFs by small molecular weight compounds may be a feasible therapeutic approach to stimulate myelination.

Comparative Transcriptome Analysis Identifies CCDC80 as a Novel Gene Associated with Pulmonary Arterial Hypertension.

Pulmonary arterial hypertension (PAH) is a heterogeneous disorder associated with a progressive increase in pulmonary artery resistance and pressure. Although various therapies have been developed, the 5-year survival rate of PAH patients remains low. There is thus an important need to identify novel genes that are commonly dysregulated in PAH of various etiologies and could be used as biomarkers and/or therapeutic targets. In this study, we performed comparative transcriptome analysis of five mammalian PAH datasets downloaded from a public database. We identified 228 differentially expressed genes (DEGs) from a rat PAH model caused by inhibition of vascular endothelial growth factor receptor under hypoxic conditions, 379 DEGs from a mouse PAH model associated with systemic sclerosis, 850 DEGs from a mouse PAH model associated with schistosomiasis, 1598 DEGs from one cohort of human PAH patients, and 4260 DEGs from a second cohort of human PAH patients. Gene-by-gene comparison identified four genes that were differentially upregulated or downregulated in parallel in all five sets of DEGs. Expression of coiled-coil domain containing 80 (CCDC80) and anterior gradient two genes was significantly increased in the five datasets, whereas expression of SMAD family member six and granzyme A was significantly decreased. Weighted gene co-expression network analysis revealed a connection between CCDC80 and collagen type I alpha 1 (COL1A1) expression. To validate the function of CCDC80 in vivo, we knocked out ccdc80 in zebrafish using the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system. In vivo imaging of zebrafish expressing a fluorescent protein in endothelial cells showed that ccdc80 deletion significantly increased the diameter of the ventral artery, a vessel supplying blood to the gills. We also demonstrated that expression of col1a1 and endothelin-1 mRNA was significantly decreased in the ccdc80-knockout zebrafish. Finally, we examined Ccdc80 immunoreactivity in a rat PAHmodel and found increased expression in the hypertrophied media and adventitia of the pre-acinar pulmonary arteries (PAs) and in the thickened intima, media, and adventitia of the obstructed intra-acinar PAs. These results suggest that increased expression of CCDC80 may be involved in the pathogenesis of PAH, potentially by modulating the expression of endothelin-1 and COL1A1.

Downregulation of GSTK1 Is a Common Mechanism Underlying Hypertrophic Cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is characterized by left ventricular hypertrophy and is associated with a number of potential outcomes, including impaired diastolic function, heart failure, and sudden cardiac death. Various etiologies have been described for HCM, including pressure overload and mutations in sarcomeric and non-sarcomeric genes. However, the molecular pathogenesis of HCM remains incompletely understood. In this study, we performed comparative transcriptome analysis to identify dysregulated genes common to five mouse HCM models of differing etiology: (i) mutation of myosin heavy chain 6, (ii) mutation of tropomyosin 1, (iii) expressing human phospholamban on a null background, (iv) knockout of frataxin, and (v) transverse aortic constriction. Gene-by-gene comparison identified five genes dysregulated in all five HCM models. Glutathione S-transferase kappa 1 (Gstk1) was significantly downregulated in the five models, whereas myosin heavy chain 7 (Myh7), connective tissue growth factor (Ctgf), periostin (Postn), and reticulon 4 (Rtn4) were significantly upregulated. Gene ontology comparison revealed that 51 cellular processes were significantly enriched in genes dysregulated in each transcriptome dataset. Among them, six processes (oxidative stress, aging, contraction, developmental process, cell differentiation, and cell proliferation) were related to four of the five genes dysregulated in all HCM models. GSTK1 was related to oxidative stress only, whereas the other four genes were related to all six cell processes except MYH7 for oxidative stress. Gene-gene functional interaction network analysis suggested correlative expression of GSTK1, MYH7, and actin alpha 2 (ACTA2). To investigate the implications of Gstk1 downregulation for cardiac function, we knocked out gstk1 in zebrafish using the clustered regularly interspaced short palindromic repeats/Cas9 system. We found that expression of the zebrafish homologs of MYH7, ACTA2, and actin alpha 1 were increased in the gstk1-knockout zebrafish. In vivo imaging of zebrafish expressing a fluorescent protein in cardiomyocytes showed that gstk1 deletion significantly decreased the end diastolic volume and, to a lesser extent, end systolic volume. These results suggest that downregulation of GSTK1 may be a common mechanism underlying HCM of various etiologies, possibly through increasing oxidative stress and the expression of sarcomere genes.

DNA Damage Response Is Involved in the Developmental Toxicity of Mebendazole in Zebrafish Retina.

Intestinal helminths cause iron-deficiency anemia in pregnant women, associated with premature delivery, low birth weight, maternal ill health, and maternal death. Although benzimidazole compounds such as mebendazole (MBZ) are highly efficacious against helminths, there are limited data on its use during pregnancy. In this study, we performed in vivo imaging of the retinas of zebrafish larvae exposed to MBZ, and found that exposure to MBZ during 2 and 3 days post-fertilization caused malformation of the retinal layers. To identify the molecular mechanism underlying the developmental toxicity of MBZ, we performed transcriptome analysis of zebrafish eyes. The analysis revealed that the DNA damage response was involved in the developmental toxicity of MBZ. We were also able to demonstrate that inhibition of ATM significantly attenuated the apoptosis induced by MBZ in the zebrafish retina. These results suggest that MBZ causes developmental toxicity in the zebrafish retina at least partly by activating the DNA damage response, including ATM signaling, providing a potential adverse outcome pathway in the developmental toxicity of MBZ in mammals.