Biomarker RIPK3 Is Silenced by Hypermethylation in Melanoma and Epigenetic Editing Reestablishes Its Tumor Suppressor Function.

For several decades, cancers have demonstrably been one of the most frequent causes of death worldwide. In addition to genetic causes, cancer can also be caused by epigenetic gene modifications. Frequently, tumor suppressor genes are epigenetically inactivated due to hypermethylation of their CpG islands, actively contributing to tumorigenesis. Since CpG islands are usually localized near promoters, hypermethylation of the promoter can have a major impact on gene expression. In this study, the potential tumor suppressor gene Receptor Interacting Serine/Threonine Protein Kinase 3 (RIPK3) was examined for an epigenetic regulation and its gene inactivation in melanomas. A hypermethylation of the RIPK3 CpG island was detected by bisulfite pyrosequencing and was accompanied by a correlated loss of its expression. In addition, an increasing RIPK3 methylation rate was observed with increasing tumor stage of melanomas. For further epigenetic characterization of RIPK3, epigenetic modulation was performed using a modified CRISPR/dCas9 (CRISPRa activation) system targeting its DNA hypermethylation. We observed a reduced fitness of melanoma cells by (re-)expression and demethylation of the RIPK3 gene using the epigenetic editing-based method. The tumor suppressive function of RIPK3 was evident by phenotypic determination using fluorescence microscopy, flow cytometry and wound healing assay. Our data highlight the function of RIPK3 as an epigenetically regulated tumor suppressor in melanoma, allowing it to be classified as a biomarker.

Multi-OMIC analysis of Huntington disease reveals a neuroprotective astrocyte state.

Huntington disease (HD) is an incurable neurodegenerative disease characterized by neuronal loss and astrogliosis. One hallmark of HD is the selective neuronal vulnerability of striatal medium spiny neurons. To date, the underlying mechanisms of this selective vulnerability have not been fully defined. Here, we employed a multi-omic approach including single nucleus RNAseq (snRNAseq), bulk RNAseq, lipidomics, HTT gene CAG repeat length measurements, and multiplexed immunofluorescence on post-mortem brain tissue from multiple brain regions of HD and control donors. We defined a signature of genes that is driven by CAG repeat length and found it enriched in astrocytic and microglial genes. Moreover, weighted gene correlation network analysis showed loss of connectivity of astrocytic and microglial modules in HD and identified modules that correlated with CAG-repeat length which further implicated inflammatory pathways and metabolism. We performed lipidomic analysis of HD and control brains and identified several lipid species that correlate with HD grade, including ceramides and very long chain fatty acids. Integration of lipidomics and bulk transcriptomics identified a consensus gene signature that correlates with HD grade and HD lipidomic abnormalities and implicated the unfolded protein response pathway. Because astrocytes are critical for brain lipid metabolism and play important roles in regulating inflammation, we analyzed our snRNAseq dataset with an emphasis on astrocyte pathology. We found two main astrocyte types that spanned multiple brain regions; these types correspond to protoplasmic astrocytes, and fibrous-like - CD44-positive, astrocytes. HD pathology was differentially associated with these cell types in a region-specific manner. One protoplasmic astrocyte cluster showed high expression of metallothionein genes, the depletion of this cluster positively correlated with the depletion of vulnerable medium spiny neurons in the caudate nucleus. We confirmed that metallothioneins were increased in cingulate HD astrocytes but were unchanged or even decreased in caudate astrocytes. We combined existing genome-wide association studies (GWAS) with a GWA study conducted on HD patients from the original Venezuelan cohort and identified a single-nucleotide polymorphism in the metallothionein gene locus associated with delayed age of onset. Functional studies found that metallothionein overexpressing astrocytes are better able to buffer glutamate and were neuroprotective of patient-derived directly reprogrammed HD MSNs as well as against rotenone-induced neuronal death in vitro. Finally, we found that metallothionein-overexpressing astrocytes increased the phagocytic activity of microglia in vitro and increased the expression of genes involved in fatty acid binding. Together, we identified an astrocytic phenotype that is regionally-enriched in less vulnerable brain regions that can be leveraged to protect neurons in HD.

Endogenous recapitulation of Alzheimer's disease neuropathology through human 3D direct neuronal reprogramming.

Alzheimer's disease (AD) is a neurodegenerative disorder that primarily affects elderly individuals, and is characterized by hallmark neuronal pathologies including extracellular amyloid-ß (Aß) plaque deposition, intracellular tau tangles, and neuronal death. However, recapitulating these age-associated neuronal pathologies in patient-derived neurons has remained a significant challenge, especially for late-onset AD (LOAD), the most common form of the disorder. Here, we applied the high efficiency microRNA-mediated direct neuronal reprogramming of fibroblasts from AD patients to generate cortical neurons in three-dimensional (3D) Matrigel and self-assembled neuronal spheroids. Our findings indicate that neurons and spheroids reprogrammed from both autosomal dominant AD (ADAD) and LOAD patients exhibited AD-like phenotypes linked to neurons, including extracellular Aß deposition, dystrophic neurites with hyperphosphorylated, K63-ubiquitin-positive, seed-competent tau, and spontaneous neuronal death in culture. Moreover, treatment with ß- or γ-secretase inhibitors in LOAD patient-derived neurons and spheroids before Aß deposit formation significantly lowered Aß deposition, as well as tauopathy and neurodegeneration. However, the same treatment after the cells already formed Aß deposits only had a mild effect. Additionally, inhibiting the synthesis of age-associated retrotransposable elements (RTEs) by treating LOAD neurons and spheroids with the reverse transcriptase inhibitor, lamivudine, alleviated AD neuropathology. Overall, our results demonstrate that direct neuronal reprogramming of AD patient fibroblasts in a 3D environment can capture age-related neuropathology and reflect the interplay between Aß accumulation, tau dysregulation, and neuronal death. Moreover, miRNA-based 3D neuronal conversion provides a human-relevant AD model that can be used to identify compounds that can potentially ameliorate AD-associated pathologies and neurodegeneration.

Recapitulation of endogenous 4R tau expression and formation of insoluble tau in directly reprogrammed human neurons.

Tau is a microtubule-binding protein expressed in neurons, and the equal ratios between 4-repeat (4R) and 3-repeat (3R) isoforms are maintained in normal adult brain function. Dysregulation of 3R:4R ratio causes tauopathy, and human neurons that recapitulate tau isoforms in health and disease will provide a platform for elucidating pathogenic processes involving tau pathology. We carried out extensive characterizations of tau isoforms expressed in human neurons derived by microRNA-induced neuronal reprogramming of adult fibroblasts. Transcript and protein analyses showed that miR neurons expressed all six isoforms with the 3R:4R isoform ratio equivalent to that detected in human adult brains. Also, miR neurons derived from familial tauopathy patients with a 3R:4R ratio altering mutation showed increased 4R tau and the formation of insoluble tau with seeding activity. Our results collectively demonstrate the utility of miRNA-induced neuronal reprogramming to recapitulate endogenous tau regulation comparable with the adult brain in health and disease.

Deconstructing Stepwise Fate Conversion of Human Fibroblasts to Neurons by MicroRNAs.

Cell-fate conversion generally requires reprogramming effectors to both introduce fate programs of the target cell type and erase the identity of starting cell population. Here, we reveal insights into the activity of microRNAs miR-9/9∗ and miR-124 (miR-9/9∗-124) as reprogramming agents that orchestrate direct conversion of human fibroblasts into motor neurons by first eradicating fibroblast identity and promoting uniform transition to a neuronal state in sequence. We identify KLF-family transcription factors as direct target genes for miR-9/9∗-124 and show their repression is critical for erasing fibroblast fate. Subsequent gain of neuronal identity requires upregulation of a small nuclear RNA, RN7SK, which induces accessibilities of chromatin regions and neuronal gene activation to push cells to a neuronal state. Our study defines deterministic components in the microRNA-mediated reprogramming cascade.

Quantitative mass spectrometric analysis of the mouse cerebral cortex after ischemic stroke.

Ischemic strokes result in the death of brain tissue and a wave of downstream effects, often leading to lifelong disabilities or death. However, the underlying mechanisms of ischemic damage and repair systems remain largely unknown. In order to better understand these mechanisms, TMT-isobaric mass tagging and mass spectrometry were conducted on brain cortex extracts from mice subjected to one hour of middle cerebral artery occlusion (MCAO) and after one hour of reperfusion. In total, 2,690 proteins were identified and quantified, out of which 65% of the top 5% of up- and down-regulated proteins were found to be significant (p < 0.05). Network-based gene ontology analysis was then utilized to cluster all identified proteins by protein functional groups and cellular roles. Although three different cellular functions were identified-organelle outer membrane proteins, cytosolic ribosome proteins, and spliceosome complex proteins-several functional domains were found to be common. Of these, organelle outer membrane proteins were downregulated whereas cytosolic ribosome and spliceosome complex proteins were upregulated, indicating that major molecular events post-stroke were translation-associated and subsequent signaling pathways (e.g., poly (ADP-ribose) (PAR) dependent cell death). By approaching stroke analyses via TMT-isobaric mass tagging, the work herein presents a grand scope of protein-based molecular mechanisms involved with ischemic stroke recovery.

FOXG1 Orchestrates Neocortical Organization and Cortico-Cortical Connections.

The hallmarks of FOXG1 syndrome, which results from mutations in a single FOXG1 allele, include cortical atrophy and corpus callosum agenesis. However, the etiology for these structural deficits and the role of FOXG1 in cortical projection neurons remain unclear. Here we demonstrate that Foxg1 in pyramidal neurons plays essential roles in establishing cortical layers and the identity and axon trajectory of callosal projection neurons. The neuron-specific actions of Foxg1 are achieved by forming a transcription complex with Rp58. The Foxg1-Rp58 complex directly binds and represses Robo1, Slit3, and Reelin genes, the key regulators of callosal axon guidance and neuronal migration. We also found that inactivation of one Foxg1 allele specifically in cortical neurons was sufficient to cause cerebral cortical hypoplasia and corpus callosum agenesis. Together, this study reveals a novel gene regulatory pathway that specifies neuronal characteristics during cerebral cortex development and sheds light on the etiology of FOXG1 syndrome. VIDEO ABSTRACT.

High-Content Genome-Wide RNAi Screen Reveals CCR3 as a Key Mediator of Neuronal Cell Death.

Neuronal loss caused by ischemic injury, trauma, or disease can lead to devastating consequences for the individual. With the goal of limiting neuronal loss, a number of cell death pathways have been studied, but there may be additional contributors to neuronal death that are yet unknown. To identify previously unknown cell death mediators, we performed a high-content genome-wide screening of short, interfering RNA (siRNA) with an siRNA library in murine neural stem cells after exposure to N-methyl-N-nitroso-N'-nitroguanidine (MNNG), which leads to DNA damage and cell death. Eighty genes were identified as key mediators for cell death. Among them, 14 are known cell death mediators and 66 have not previously been linked to cell death pathways. Using an integrated approach with functional and bioinformatics analysis, we provide possible molecular networks, interconnected pathways, and/or protein complexes that may participate in cell death. Of the 66 genes, we selected CCR3 for further evaluation and found that CCR3 is a mediator of neuronal injury. CCR3 inhibition or deletion protects murine cortical cultures from oxygen-glucose deprivation-induced cell death, and CCR3 deletion in mice provides protection from ischemia in vivo. Taken together, our findings suggest that CCR3 is a previously unknown mediator of cell death. Future identification of the neural cell death network in which CCR3 participates will enhance our understanding of the molecular mechanisms of neural cell death.

Critical Roles of the Histone Methyltransferase MLL4/KMT2D in Murine Hepatic Steatosis Directed by ABL1 and PPARγ2.

The pathophysiologic continuum of non-alcoholic fatty liver disease begins with steatosis. Despite recent advances in our understanding of the gene regulatory program directing steatosis, how it is orchestrated at the chromatin level is unclear. PPARγ2 is a hepatic steatotic transcription factor induced by overnutrition. Here, we report that the histone H3 lysine 4 methyltransferase MLL4/KMT2D directs overnutrition-induced murine steatosis via its coactivator function for PPARγ2. We demonstrate that overnutrition facilitates the recruitment of MLL4 to steatotic target genes of PPARγ2 and their transactivation via H3 lysine 4 methylation because PPARγ2 phosphorylated by overnutrition-activated ABL1 kinase shows enhanced interaction with MLL4. We further show that Pparg2 (encoding PPARγ2) is also a hepatic target gene of ABL1-PPARγ2-MLL4. Consistently, inhibition of ABL1 improves the fatty liver condition of mice with overnutrition by suppressing the pro-steatotic action of MLL4. Our results uncover a murine hepatic steatosis regulatory axis consisting of ABL1-PPARγ2-MLL4, which may serve as a target of anti-steatosis drug development.

Chx10 Consolidates V2a Interneuron Identity through Two Distinct Gene Repression Modes.

During development, two cell types born from closely related progenitor pools often express identical transcriptional regulators despite their completely distinct characteristics. This phenomenon implies the need for a mechanism that operates to segregate the identities of the two cell types throughout differentiation after initial fate commitment. To understand this mechanism, we investigated the fate specification of spinal V2a interneurons, which share important developmental genes with motor neurons (MNs). We demonstrate that the paired homeodomain factor Chx10 functions as a critical determinant for V2a fate and is required to consolidate V2a identity in postmitotic neurons. Chx10 actively promotes V2a fate, downstream of the LIM-homeodomain factor Lhx3, while concomitantly suppressing the MN developmental program by preventing the MN-specific transcription complex from binding and activating MN genes. This dual activity enables Chx10 to effectively separate the V2a and MN pathways. Our study uncovers a widely applicable gene regulatory principle for segregating related cell fates.

miR-218 is essential to establish motor neuron fate as a downstream effector of Isl1-Lhx3.

While microRNAs have emerged as an important component of gene regulatory networks, it remains unclear how microRNAs collaborate with transcription factors in the gene networks that determines neuronal cell fate. Here we show that in the developing spinal cord, the expression of miR-218 is directly upregulated by the Isl1-Lhx3 complex, which drives motor neuron fate. Inhibition of miR-218 suppresses the generation of motor neurons in both chick neural tube and mouse embryonic stem cells, suggesting that miR-218 plays a crucial role in motor neuron differentiation. Results from unbiased RISC-trap screens, in vivo reporter assays and overexpression studies indicated that miR-218 directly represses transcripts that promote developmental programs for interneurons. In addition, we found that miR-218 activity is required for Isl1-Lhx3 to effectively induce motor neurons and suppress interneuron fates. Together our results reveal an essential role of miR-218 as a downstream effector of the Isl1-Lhx3 complex in establishing motor neuron identity.

Vaccinia-based influenza vaccine overcomes previously induced immunodominance hierarchy for heterosubtypic protection.

Growing concerns about unpredictable influenza pandemics require a broadly protective vaccine against diverse influenza strains. One of the promising approaches was a T cell-based vaccine, but the narrow breadth of T-cell immunity due to the immunodominance hierarchy established by previous influenza infection and efficacy against only mild challenge condition are important hurdles to overcome. To model T-cell immunodominance hierarchy in humans in an experimental setting, influenza-primed C57BL/6 mice were chosen and boosted with a mixture of vaccinia recombinants, individually expressing consensus sequences from avian, swine, and human isolates of influenza internal proteins. As determined by IFN-γ ELISPOT and polyfunctional cytokine secretion, the vaccinia recombinants of influenza expanded the breadth of T-cell responses to include subdominant and even minor epitopes. Vaccine groups were successfully protected against 100 LD50 challenges with PR/8/34 and highly pathogenic avian influenza H5N1, which contained the identical dominant NP366 epitope. Interestingly, in challenge with pandemic A/Cal/04/2009 containing mutations in the dominant epitope, only the group vaccinated with rVV-NP + PA showed improved protection. Taken together, a vaccinia-based influenza vaccine expressing conserved internal proteins improved the breadth of influenza-specific T-cell immunity and provided heterosubtypic protection against immunologically close as well as distant influenza strains.

Gene-set based genome-wide association analysis for the speed of sound in two skeletal sites of Korean women.

The speed of sound (SOS) value is an indicator of bone mineral density (BMD). Previous genome-wide association (GWA) studies have identified a number of genes, whose variations may affect BMD levels. However, their biological implications have been elusive. We re-analyzed the GWA study dataset for the SOS values in skeletal sites of 4,659 Korean women, using a gene-set analysis software, GSA-SNP. We identified 10 common representative GO terms, and 17 candidate genes between these two traits (PGS ＜ 0.05). Implication of these GO terms and genes in the bone mechanism is well supported by the literature survey. Interestingly, the significance levels of some member genes were inversely related, in several gene-sets that were shared between two skeletal sites. This implies that biological process, rather than SNP or gene, is the substantial unit of genetic association for SOS in bone. In conclusion, our findings may provide new insights into the biological mechanisms for BMD.

Gene set analyses of genome-wide association studies on 49 quantitative traits measured in a single genetic epidemiology dataset.

Gene set analysis is a powerful tool for interpreting a genome-wide association study result and is gaining popularity these days. Comparison of the gene sets obtained for a variety of traits measured from a single genetic epidemiology dataset may give insights into the biological mechanisms underlying these traits. Based on the previously published single nucleotide polymorphism (SNP) genotype data on 8,842 individuals enrolled in the Korea Association Resource project, we performed a series of systematic genome-wide association analyses for 49 quantitative traits of basic epidemiological, anthropometric, or blood chemistry parameters. Each analysis result was subjected to subsequent gene set analyses based on Gene Ontology (GO) terms using gene set analysis software, GSA-SNP, identifying a set of GO terms significantly associated to each trait (pcorr < 0.05). Pairwise comparison of the traits in terms of the semantic similarity in their GO sets revealed surprising cases where phenotypically uncorrelated traits showed high similarity in terms of biological pathways. For example, the pH level was related to 7 other traits that showed low phenotypic correlations with it. A literature survey implies that these traits may be regulated partly by common pathways that involve neuronal or nerve systems.

Performance Comparison of Two Gene Set Analysis Methods for Genome-wide Association Study Results: GSA-SNP vs i-GSEA4GWAS.

Gene set analysis (GSA) is useful in interpreting a genome-wide association study (GWAS) result in terms of biological mechanism. We compared the performance of two different GSA implementations that accept GWAS p-values of single nucleotide polymorphisms (SNPs) or gene-by-gene summaries thereof, GSA-SNP and i-GSEA4GWAS, under the same settings of inputs and parameters. GSA runs were made with two sets of p-values from a Korean type 2 diabetes mellitus GWAS study: 259,188 and 1,152,947 SNPs of the original and imputed genotype datasets, respectively. When Gene Ontology terms were used as gene sets, i-GSEA4GWAS produced 283 and 1,070 hits for the unimputed and imputed datasets, respectively. On the other hand, GSA-SNP reported 94 and 38 hits, respectively, for both datasets. Similar, but to a lesser degree, trends were observed with Kyoto Encyclopedia of Genes and Genomes (KEGG) gene sets as well. The huge number of hits by i-GSEA4GWAS for the imputed dataset was probably an artifact due to the scaling step in the algorithm. The decrease in hits by GSA-SNP for the imputed dataset may be due to the fact that it relies on Z-statistics, which is sensitive to variations in the background level of associations. Judicious evaluation of the GSA outcomes, perhaps based on multiple programs, is recommended.

Ethnic variability in the allelic distribution of pharmacogenes between Korean and other populations.

OBJECTIVE: We examined the differences in allele frequencies for pharmacogenes among the Korean (KOR), Chinese (CHB), Japanese (JPT), Caucasian (CEU), and Nigerian (YRI) populations. METHODS: Fifty-seven pharmacogenes were selected from the imputed Korean Association REsource and HapMap databases. Minor allele frequencies were analyzed using the sample size-modified single nucleotide polymorphism-specific fixation index (FST) and the χ-test with Bonferroni's correction. Geneset analysis was also carried out to identify pharmacogenes that have significantly different allele frequencies among the various populations tested. RESULTS: The KOR population was the most divergent group from the YRI population (FST: 0.079) but very similar to the CHB and JPT populations (FST: 0.003). VKORC1 showed a large population divergence in the KOR-YRI (0.439) comparison. CYP3A4 was also highly divergent in the KOR-YRI (FST: 0.361) comparison. The calcium signaling pathway gene set was divergent in all pairwise population comparisons. CONCLUSION: In terms of the 57 pharmacogenes studied, there were no significant differences among the KOR, CHB, and JPT populations. However, the YRI and CEU populations were significantly differentiated from the three Eastern Asian groups. Future pharmacogenomics studies can utilize the polymorphisms identified in this study, as these variants may have important implications for the selection of highly informative single nucleotide polymorphisms for future clinical trials.

Generation of reassortant influenza viruses within the non-industrial poultry system.

We compared the genetic and biologic characteristics of 35 influenza viruses of different epidemiological backgrounds in Korea, including H3N2 canine influenza virus (CIV). Phylogenetic analysis revealed that chicken adapted H9N2 viruses (A/chicken/Korea/96006/96 [CK/Kor/96006-like]) have acquired aquatic avian gene segments through reassortment, and these reassorted H9N2 viruses were more frequently detected from minor poultry species than from industrial poultry. Conversely, gene segments from CK/Kor/96006-like viruses were also detected in most of the viruses from domestic ducks. Interestingly, domestic ducks, rather than wild aquatic birds, harbored close relatives of all eight gene segments of H3N2 CIV, which preferred binding to avian receptors. Therefore, bidirectional virus transmission events are assumed to have occurred between land-based poultry and aquatic poultry, in particular within the non-industrial poultry system. These events have contributed to the generation of a novel reassortant, H3N2 CIV. To prevent generating other reassortants capable of interspecies transmission, gene movements in the non-industrial poultry systems should be clarified and managed.

Posterior circulation stroke after bronchial artery embolization.

We report on a 66-year-old woman with a posterior circulation stroke that occurred after bronchial artery embolization (BAE) due to post-tuberculous bronchiectasis. Stroke is a rare complication of BAE and is usually thought to be caused by inadvertent embolization via a bronchial artery-pulmonary vein shunt. However, the possibility of thromboembolic stroke should be considered, because of the patient's possible underlying anatomical variations or atherothrombosis.

Field application of the H9M2e enzyme-linked immunosorbent assay for differentiation of H9N2 avian influenza virus-infected chickens from vaccinated chickens.

Vaccination for control of H9N2 low-pathogenicity avian influenza (LPAI) in chickens began in 2007 in South Korea where the H9N2 virus is prevalent. Recently, an enzyme-linked immunosorbent assay (ELISA) using the extracellular domain of the M2 protein (M2e ELISA) was developed as another strategy to differentiate between vaccinated and infected chickens. Here, an ELISA using the extracellular domain of the M2 protein of H9N2 LPAI virus (H9M2e ELISA) was applied to differentiate infected from vaccinated chickens using the H9N2 LPAI virus M2 peptide. The specificity and sensitivity of the optimized H9M2e ELISA were 96.1% and 83.8% (the absorbance of the sample to the absorbance for the positive control [S/P ratio] ≥ 0.6), respectively, with the cutoff value (S/P ratio = 0.6), and the criterion of avian influenza (AI) infection in a chicken house was established as >20% reactivity of anti-M2e antibody per house with this cutoff value. After infection in naïve chickens and once-vaccinated chickens with a hemagglutination inhibition (HI) assay titer of 9.25 ± 0.75 log(2) units, the sera from infected chickens were confirmed as AI infected when the chickens were 1 week old in both groups, and AI infection lasted for 24 weeks and 9 weeks in naïve and once-vaccinated chickens, respectively, although in twice-vaccinated chickens with a higher HI titer of 11.17 ± 0.37 log(2) units, anti-M2e antibody in infected sera did not reach a level indicating AI infection. In field application, anti-M2e antibody produced in infected chickens after vaccination or in reinfected chickens could be identified as AI infection, although HI test could not distinguish infected from vaccinated sera. These results indicate the utility of H9M2e ELISA as a surveillance tool in control of H9N2 LPAI infections.