Huntingtin loss in hepatocytes is associated with altered metabolism, adhesion, and liver zonation.

Huntington's disease arises from a toxic gain of function in the huntingtin (HTT) gene. As a result, many HTT-lowering therapies are being pursued in clinical studies, including those that reduce HTT RNA and protein expression in the liver. To investigate potential impacts, we characterized molecular, cellular, and metabolic impacts of chronic HTT lowering in mouse hepatocytes. Lifelong hepatocyte HTT loss is associated with multiple physiological changes, including increased circulating bile acids, cholesterol and urea, hypoglycemia, and impaired adhesion. HTT loss causes a clear shift in the normal zonal patterns of liver gene expression, such that pericentral gene expression is reduced. These alterations in liver zonation in livers lacking HTT are observed at the transcriptional, histological, and plasma metabolite levels. We have extended these phenotypes physiologically with a metabolic challenge of acetaminophen, for which the HTT loss results in toxicity resistance. Our data reveal an unexpected role for HTT in regulating hepatic zonation, and we find that loss of HTT in hepatocytes mimics the phenotypes caused by impaired hepatic ß-catenin function.

Huntingtin loss in hepatocytes is associated with altered metabolism, adhesion, and liver zonation.

Huntington's disease arises from a toxic gain of function in the huntingtin ( HTT ) gene. As a result, many HTT-lowering therapies are being pursued in clinical studies, including those that reduce HTT RNA and protein expression in the liver. To investigate potential impacts, we characterized molecular, cellular, and metabolic impacts of chronic HTT lowering in mouse hepatocytes. Lifelong hepatocyte HTT loss is associated with multiple physiological changes, including increased circulating bile acids, cholesterol and urea, hypoglycemia, and impaired adhesion. HTT loss causes a clear shift in the normal zonal patterns of liver gene expression, such that pericentral gene expression is reduced. These alterations in liver zonation in livers lacking HTT are observed at the transcriptional, histological and plasma metabolite level. We have extended these phenotypes physiologically with a metabolic challenge of acetaminophen, for which the HTT loss results in toxicity resistance. Our data reveal an unexpected role for HTT in regulating hepatic zonation, and we find that loss of HTT in hepatocytes mimics the phenotypes caused by impaired hepatic ß-catenin function.

Mapping brain gene coexpression in daytime transcriptomes unveils diurnal molecular networks and deciphers perturbation gene signatures.

Brain tissue transcriptomes may be organized into gene coexpression networks, but their underlying biological drivers remain incompletely understood. Here, we undertook a large-scale transcriptomic study using 508 wild-type mouse striatal tissue samples dissected exclusively in the afternoons to define 38 highly reproducible gene coexpression modules. We found that 13 and 11 modules are enriched in cell-type and molecular complex markers, respectively. Importantly, 18 modules are highly enriched in daily rhythmically expressed genes that peak or trough with distinct temporal kinetics, revealing the underlying biology of striatal diurnal gene networks. Moreover, the diurnal coexpression networks are a dominant feature of daytime transcriptomes in the mouse cortex. We next employed the striatal coexpression modules to decipher the striatal transcriptomic signatures from Huntington's disease models and heterozygous null mice for 52 genes, uncovering novel functions for Prkcq and Kdm4b in oligodendrocyte differentiation and bipolar disorder-associated Trank1 in regulating anxiety-like behaviors and nocturnal locomotion.

Developing HDAC4-Selective Protein Degraders To Investigate the Role of HDAC4 in Huntington's Disease Pathology.

Huntington's disease (HD) is a lethal autosomal dominant neurodegenerative disorder resulting from a CAG repeat expansion in the huntingtin (HTT) gene. The product of translation of this gene is a highly aggregation-prone protein containing a polyglutamine tract >35 repeats (mHTT) that has been shown to colocalize with histone deacetylase 4 (HDAC4) in cytoplasmic inclusions in HD mouse models. Genetic reduction of HDAC4 in an HD mouse model resulted in delayed aggregation of mHTT, along with amelioration of neurological phenotypes and extended lifespan. To further investigate the role of HDAC4 in cellular models of HD, we have developed bifunctional degraders of the protein and report the first potent and selective degraders of HDAC4 that show an effect in multiple cell lines, including HD mouse model-derived cortical neurons. These degraders act via the ubiquitin-proteasomal pathway and selectively degrade HDAC4 over other class IIa HDAC isoforms (HDAC5, HDAC7, and HDAC9).

HDinHD: A Rich Data Portal for Huntington's Disease Research.

HDinHD (Huntington's Disease in High Definition; HDinHD.org) is an open online portal for the HD research community that presents a synthesized view of HD-related scientific data. Here, we present a broad overview of HDinHD and highlight the newly launched HDinHD Explorer tool that enables researchers to discover and explore a wide range of diverse yet interconnected HD-related data. We demonstrate the utility of HDinHD Explorer through data mining of a single collection of newly released in vivo therapeutic intervention study reports alongside previously published reports.

Shape deformation analysis reveals the temporal dynamics of cell-type-specific homeostatic and pathogenic responses to mutant huntingtin.

Loss of cellular homeostasis has been implicated in the etiology of several neurodegenerative diseases (NDs). However, the molecular mechanisms that underlie this loss remain poorly understood on a systems level in each case. Here, using a novel computational approach to integrate dimensional RNA-seq and in vivo neuron survival data, we map the temporal dynamics of homeostatic and pathogenic responses in four striatal cell types of Huntington's disease (HD) model mice. This map shows that most pathogenic responses are mitigated and most homeostatic responses are decreased over time, suggesting that neuronal death in HD is primarily driven by the loss of homeostatic responses. Moreover, different cell types may lose similar homeostatic processes, for example, endosome biogenesis and mitochondrial quality control in Drd1-expressing neurons and astrocytes. HD relevance is validated by human stem cell, genome-wide association study, and post-mortem brain data. These findings provide a new paradigm and framework for therapeutic discovery in HD and other NDs.

DNA methylation study of Huntington's disease and motor progression in patients and in animal models.

Although Huntington's disease (HD) is a well studied Mendelian genetic disorder, less is known about its associated epigenetic changes. Here, we characterize DNA methylation levels in six different tissues from 3 species: a mouse huntingtin (Htt) gene knock-in model, a transgenic HTT sheep model, and humans. Our epigenome-wide association study (EWAS) of human blood reveals that HD mutation status is significantly (p < 10-7) associated with 33 CpG sites, including the HTT gene (p = 6.5 × 10-26). These Htt/HTT associations were replicated in the Q175 Htt knock-in mouse model (p = 6.0 × 10-8) and in the transgenic sheep model (p = 2.4 × 10-88). We define a measure of HD motor score progression among manifest HD cases based on multiple clinical assessments. EWAS of motor progression in manifest HD cases exhibits significant (p < 10-7) associations with methylation levels at three loci: near PEX14 (p = 9.3 × 10-9), GRIK4 (p = 3.0 × 10-8), and COX4I2 (p = 6.5 × 10-8). We conclude that HD is accompanied by profound changes of DNA methylation levels in three mammalian species.

Ultrasensitive quantitative measurement of huntingtin phosphorylation at residue S13.

Huntington's disease (HD) is a progressive neurodegenerative disorder caused by an expansion of a CAG triplet repeat (encoding for a polyglutamine tract) within the first exon of the huntingtin gene. Expression of the mutant huntingtin (mHTT) protein can result in the production of N-terminal fragments with a robust propensity to form oligomers and aggregates, which may be causally associated with HD pathology. Several lines of evidence indicate that N17 phosphorylation or pseudophosphorylation at any of the residues T3, S13 or S16, alone or in combination, modulates mHTT aggregation, subcellular localization and toxicity. Consequently, increasing N17 phosphorylation has been proposed as a potential therapeutic approach. However, developing genetic/pharmacological tools to quantify these phosphorylation events is necessary in order to subsequently develop tool modulators, which is difficult given the transient and incompletely penetrant nature of such post-translational modifications. Here we describe the first ultrasensitive sandwich immunoassay that quantifies HTT phosphorylated at residue S13 and demonstrate its utility for specific analyte detection in preclinical models of HD.

MicroRNA signatures of endogenous Huntingtin CAG repeat expansion in mice.

In Huntington's disease (HD) patients and in model organisms, messenger RNA transcriptome has been extensively studied; in contrast, comparatively little is known about expression and potential role of microRNAs. Using RNA-sequencing, we have quantified microRNA expression in four brain regions and liver, at three different ages, from an allelic series of HD model mice with increasing CAG length in the endogenous Huntingtin gene. Our analyses reveal CAG length-dependent microRNA expression changes in brain, with 159 microRNAs selectively altered in striatum, 102 in cerebellum, 51 in hippocampus, and 45 in cortex. In contrast, a progressive CAG length-dependent microRNA dysregulation was not observed in liver. We further identify microRNAs whose transcriptomic response to CAG length expansion differs significantly among the brain regions and validate our findings in data from a second, independent cohort of mice. Using existing mRNA expression data from the same animals, we assess the possible relationships between microRNA and mRNA expression and highlight candidate microRNAs that are negatively correlated with, and whose predicted targets are enriched in, CAG-length dependent mRNA modules. Several of our top microRNAs (Mir212/Mir132, Mir218, Mir128 and others) have been previously associated with aspects of neuronal development and survival. This study provides an extensive resource for CAG length-dependent changes in microRNA expression in disease-vulnerable and -resistant brain regions in HD mice, and provides new insights for further investigation of microRNAs in HD pathogenesis and therapeutics.

Genetic invalidation of Lp-PLA2 as a therapeutic target: Large-scale study of five functional Lp-PLA2-lowering alleles.

Aims Darapladib, a potent inhibitor of lipoprotein-associated phospholipase A2 (Lp-PLA2), has not reduced risk of cardiovascular disease outcomes in recent randomized trials. We aimed to test whether Lp-PLA2 enzyme activity is causally relevant to coronary heart disease. Methods In 72,657 patients with coronary heart disease and 110,218 controls in 23 epidemiological studies, we genotyped five functional variants: four rare loss-of-function mutations (c.109+2T > C (rs142974898), Arg82His (rs144983904), Val279Phe (rs76863441), Gln287Ter (rs140020965)) and one common modest-impact variant (Val379Ala (rs1051931)) in PLA2G7, the gene encoding Lp-PLA2. We supplemented de-novo genotyping with information on a further 45,823 coronary heart disease patients and 88,680 controls in publicly available databases and other previous studies. We conducted a systematic review of randomized trials to compare effects of darapladib treatment on soluble Lp-PLA2 activity, conventional cardiovascular risk factors, and coronary heart disease risk with corresponding effects of Lp-PLA2-lowering alleles. Results Lp-PLA2 activity was decreased by 64% ( p = 2.4 × 10-25) with carriage of any of the four loss-of-function variants, by 45% ( p < 10-300) for every allele inherited at Val279Phe, and by 2.7% ( p = 1.9 × 10-12) for every allele inherited at Val379Ala. Darapladib 160 mg once-daily reduced Lp-PLA2 activity by 65% ( p < 10-300). Causal risk ratios for coronary heart disease per 65% lower Lp-PLA2 activity were: 0.95 (0.88-1.03) with Val279Phe; 0.92 (0.74-1.16) with carriage of any loss-of-function variant; 1.01 (0.68-1.51) with Val379Ala; and 0.95 (0.89-1.02) with darapladib treatment. Conclusions In a large-scale human genetic study, none of a series of Lp-PLA2-lowering alleles was related to coronary heart disease risk, suggesting that Lp-PLA2 is unlikely to be a causal risk factor.

Huntington's disease accelerates epigenetic aging of human brain and disrupts DNA methylation levels.

Age of Huntington's disease (HD) motoric onset is strongly related to the number of CAG trinucleotide repeats in the huntingtin gene, suggesting that biological tissue age plays an important role in disease etiology. Recently, a DNA methylation based biomarker of tissue age has been advanced as an epigenetic aging clock. We sought to inquire if HD is associated with an accelerated epigenetic age. DNA methylation data was generated for 475 brain samples from various brain regions of 26 HD cases and 39 controls. Overall, brain regions from HD cases exhibit a significant epigenetic age acceleration effect (p=0.0012). A multivariate model analysis suggests that HD status increases biological age by 3.2 years. Accelerated epigenetic age can be observed in specific brain regions (frontal lobe, parietal lobe, and cingulate gyrus). After excluding controls, we observe a negative correlation (r=-0.41, p=5.5×10-8) between HD gene CAG repeat length and the epigenetic age of HD brain samples. Using correlation network analysis, we identify 11 co-methylation modules with a significant association with HD status across 3 broad cortical regions. In conclusion, HD is associated with an accelerated epigenetic age of specific brain regions and more broadly with substantial changes in brain methylation levels.

Plasma Levels of Risk-Variant APOL1 Do Not Associate with Renal Disease in a Population-Based Cohort.

Two common missense variants in APOL1 (G1 and G2) have been definitively linked to CKD in black Americans. However, not all individuals with the renal-risk genotype develop CKD, and little is known about how APOL1 variants drive disease. Given the association of APOL1 with HDL particles, which are cleared by the kidney, differences in the level or quality of mutant APOL1­HDL particles could be causal for disease and might serve as a useful risk stratification marker. We measured plasma levels of G0 (low risk), G1, and G2 APOL1 in 3450 individuals in the Dallas Heart Study using a liquid chromatography-MS method that enabled quantitation of the different variants. Additionally, we characterized native APOL1­HDL from donors with no or two APOL1 risk alleles by size-exclusion chromatography and analysis of immunopurified APOL1­HDL particles. Finally, we identified genetic loci associated with plasma APOL1 levels and tested for APOL1-dependent association with renal function. Although we replicated the previous association between APOL1 variant status and renal function in nondiabetic individuals, levels of circulating APOL1 did not associate with microalbuminuria or GFR. Furthermore, the size or known components of APOL1­HDL did not consistently differ in subjects with the renal-risk genotype. Genetic association studies implicated variants in loci harboring haptoglobin-related protein (HPR), APOL1, and ubiquitin D (UBD) in the regulation of plasma APOL1 levels, but these variants did not associate with renal function. Collectively, these data demonstrate that the risk of renal disease associated with APOL1 is probably not related to circulating levels of the mutant protein.

Exome-wide association study identifies a TM6SF2 variant that confers susceptibility to nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is the most common form of liver disease. To elucidate the molecular basis of NAFLD, we performed an exome-wide association study of liver fat content. Three variants were associated with higher liver fat levels at the exome-wide significance level of 3.6 × 10(-7): two in PNPLA3, an established locus for NAFLD, and one (encoding p.Glu167Lys) in TM6SF2, a gene of unknown function. The TM6SF2 variant encoding p.Glu167Lys was also associated with higher circulating levels of alanine transaminase, a marker of liver injury, and with lower levels of low-density lipoprotein-cholesterol (LDL-C), triglycerides and alkaline phosphatase in 3 independent populations (n > 80,000). When recombinant protein was expressed in cultured hepatocytes, 50% less Glu167Lys TM6SF2 protein was produced relative to wild-type TM6SF2. Adeno-associated virus-mediated short hairpin RNA knockdown of Tm6sf2 in mice increased liver triglyceride content by threefold and decreased very-low-density lipoprotein (VLDL) secretion by 50%. Taken together, these data indicate that TM6SF2 activity is required for normal VLDL secretion and that impaired TM6SF2 function causally contributes to NAFLD.

Demonstration of diet-induced decoupling of fatty acid and cholesterol synthesis by combining gene expression array and 2H2O quantification.

The liver is a crossroad for metabolism of lipid and carbohydrates, with acetyl-CoA serving as an important metabolic intermediate and a precursor for fatty acid and cholesterol biosynthesis pathways. A better understanding of the regulation of these pathways requires an experimental approach that provides both quantitative metabolic flux measurements and mechanistic insight. Under conditions of high carbohydrate availability, excess carbon is converted into free fatty acids and triglyceride for storage, but it is not clear how excessive carbohydrate availability affects cholesterol biosynthesis. To address this, C57BL/6J mice were fed either a low-fat, high-carbohydrate diet or a high-fat, carbohydrate-free diet. At the end of the dietary intervention, the two groups received (2)H(2)O to trace de novo fatty acid and cholesterol synthesis, and livers were collected for gene expression analysis. Expression of lipid and glucose metabolism genes was determined using a custom-designed pathway focused PCR-based gene expression array. The expression analysis showed downregulation of cholesterol biosynthesis genes and upregulation of fatty acid synthesis genes in mice receiving the high-carbohydrate diet compared with the carbohydrate-free diet. In support of these findings, (2)H(2)O tracer data showed that fatty acid synthesis was increased 10-fold and cholesterol synthesis was reduced by 1.6-fold in mice fed the respective diets. In conclusion, by applying gene expression analysis and tracer methodology, we show that fatty acid and cholesterol synthesis are differentially regulated when the carbohydrate intake in mice is altered.

Meiotic functions of RAD18.

RAD18 is an ubiquitin ligase that is involved in replication damage bypass and DNA double-strand break (DSB) repair processes in mitotic cells. Here, we investigated the testicular phenotype of Rad18-knockdown mice to determine the function of RAD18 in meiosis, and in particular, in the repair of meiotic DSBs induced by the meiosis-specific topoisomerase-like enzyme SPO11. We found that RAD18 is recruited to a specific subfraction of persistent meiotic DSBs. In addition, RAD18 is recruited to the chromatin of the XY chromosome pair, which forms the transcriptionally silent XY body. At the XY body, RAD18 mediates the chromatin association of its interaction partners, the ubiquitin-conjugating enzymes HR6A and HR6B. Moreover, RAD18 was found to regulate the level of dimethylation of histone H3 at Lys4 and maintain meiotic sex chromosome inactivation, in a manner similar to that previously observed for HR6B. Finally, we show that RAD18 and HR6B have a role in the efficient repair of a small subset of meiotic DSBs.

AAV8-mediated long-term expression of human LCAT significantly improves lipid profiles in hCETP;Ldlr(+/-) mice.

Lecithin:cholesterol acyltransferase (LCAT) is the key circulating enzyme responsible for high-density lipoprotein (HDL) cholesterol esterification, HDL maturation, and potentially reverse cholesterol transport. To further explore LCAT's mechanism of action on lipoprotein metabolism, we employed adeno-associated viral vector (AAV) serotype 8 to achieve long-term (32-week) high level expression of human LCAT in hCETP;Ldlr(+/-) mice, and characterized the lipid profiles in detail. The mice had a marked increase in HDL cholesterol, HDL particle size, and significant reduction in low-density lipoprotein (LDL) cholesterol, plasma triglycerides, and plasma apoB. Plasma LCAT activity significantly increased with humanized substrate specificity. HDL cholesteryl esters increased in a fashion that fits human LCAT specificity. HDL phosphatidylcholines trended toward decrease, with no change observed for HDL lysophosphatidylcholines. Triglycerides reduction appeared to reside in all lipoprotein particles (very low-density lipoprotein (VLDL), LDL, and HDL), with HDL triglycerides composition highly reflective of VLDL, suggesting that changes in HDL triglycerides were primarily driven by the altered triglycerides metabolism in VLDL. In summary, in this human-like model for lipoprotein metabolism, AAV8-mediated overexpression of human LCAT resulted in profound changes in plasma lipid profiles. Detailed lipid analyses in the lipoprotein particles suggest that LCAT's beneficial effect on lipid metabolism includes not only enhanced HDL cholesterol esterification but also improved metabolism of apoB-containing particles and triglycerides. Our findings thus shed new light on LCAT's mechanism of action and lend support to its therapeutic potential in treating dyslipidemia.

PDK1 attenuation fails to prevent tumor formation in PTEN-deficient transgenic mouse models.

PDK1 activates AKT suggesting that PDK1 inhibition might suppress tumor development. However, while PDK1 has been investigated intensively as an oncology target, selective inhibitors suitable for in vivo studies have remained elusive. In this study we present the results of in vivo PDK1 inhibition through a universally applicable RNAi approach for functional drug target validation in oncogenic pathway contexts. This approach, which relies on doxycycline-inducible shRNA expression from the Rosa26 locus, is ideal for functional studies of genes like PDK1 where constitutive mouse models lead to strong developmental phenotypes or embryonic lethality. We achieved more than 90% PDK1 knockdown in vivo, a level sufficient to impact physiological functions resulting in hyperinsulinemia and hyperglycemia. This phenotype was reversible on PDK1 reexpression. Unexpectedly, long-term PDK1 knockdown revealed a lack of potent antitumor efficacy in 3 different mouse models of PTEN-deficient cancer. Thus, despite efficient PDK1 knockdown, inhibition of the PI3K pathway was marginal suggesting that PDK1 was not a rate limiting factor. Ex vivo analysis of pharmacological inhibitors revealed that AKT and mTOR inhibitors undergoing clinical development are more effective than PDK1 inhibitors at blocking activated PI3K pathway signaling. Taken together our findings weaken the widely held expectation that PDK1 represents an appealing oncology target.

In vivo analysis of gene knockdown in tetracycline-inducible shRNA mice.

Expression of small hairpin RNA (shRNA) in mammalian cells can trigger potent RNAi-mediated gene silencing. The dominant-acting RNAi can often result in phenotypes similar to that of a null allele. Moreover, the generation of shRNA knockdown mice and subsequent phenotypic analysis can be achieved in a condensed timeline compared to that of conventional gene-targeting knockout strategies. Here, we discuss methods for the in vivo analysis of gene function in adult mouse tissues following tetracycline-induced RNA knockdown in a single-copy inducible polymerase III promoter-driven shRNA system.

Manic fringe is not required for embryonic development, and fringe family members do not exhibit redundant functions in the axial skeleton, limb, or hindbrain.

Tight regulation of Notch pathway signaling is important in many aspects of embryonic development. Notch signaling can be modulated by expression of fringe genes, encoding glycosyltransferases that modify EGF repeats in the Notch receptor. Although Lunatic fringe (Lfng) has been shown to play important roles in vertebrate segmentation, comparatively little is known regarding the developmental functions of the other vertebrate fringe genes, Radical fringe (Rfng) and Manic fringe (Mfng). Here we report that Mfng expression is not required for embryonic development. Further, we find that despite significant overlap in expression patterns, we detect no obvious synergistic defects in mice in the absence of two, or all three, fringe genes during development of the axial skeleton, limbs, hindbrain, and cranial nerves.

Dual role of Nr2e3 in photoreceptor development and maintenance.

Nr2e3 is a photoreceptor-specific nuclear receptor believed to play a role in photoreceptor development, differentiation, and survival. Much research has focused on the interaction of Nr2e3 with other transcription factors in determining the milieu of target gene expression in photoreceptors of the neonatal and adult retina. To investigate the downstream targets of Nr2e3 and thereby shed light on the functional pathways relevant to photoreceptor development and maintenance, expression profiling was performed on retinas from two different mouse knockout lines, one containing a targeted disruption of the Nr2e3 gene (Nr2e3 -/-), the other containing a spontaneous null allele of the Nr2e3 locus (rd7). Using whole genome microarrays, mRNA expression profiles of retinas from the two mutant strains were compared to those of wildtype C57BL/6 mice over a time course that ranged from postnatal day (p) 2 to 6months of age (p180). Additionally, expression profiling was performed on retinal explants treated with a putative NR2E3 agonist. The molecular profiling of Nr2e3 -/- and rd7/rd7 retinas identified 281 putative Nr2e3-dependent genes that were differentially expressed between wildtype and mutant retinas during at least one time point. Consistent with previous reports that Nr2e3 is necessary for the repression of cone-specific genes, increased expression of cone-specific genes was observed in the mutant samples, thereby providing proof-of-concept for the microarray screen. Further annotation of these data sets revealed ten predominant functional classes involved in the Nr2e3-mediated development and/or maintenance of photoreceptors. Interestingly, differences in the expression of Nr2e3-dependent genes exhibited two distinct temporal patterns. One group of genes showed a sustained difference in expression as compared to wildtype over the entire time course of the study, whereas a second group showed only transient differences which were largest around p10. Comparison of gene expression changes in Nr2e3 -/- and rd7/rd7 retinas with those uncovered by treating retinal explants with a putative NR2E3 agonist revealed four genes that were down-regulated in mutant retinas that lack Nr2e3 function but were up-regulated in agonist-treated explants. These results strongly suggest that the four genes may be direct targets of Nr2e3. Our identification of two sets of Nr2e3-regulated genes provides further evidence of a dual role for Nr2e3 in specification of photoreceptor fate during development as well as photoreceptor maintenance in the adult.