Adult-onset degeneration of adipose tissue in mice deficient for the Sox8 transcription factor

Transcriptomic profiling reveals key early response genes during GDF6-mediated differentiation of human adipose-derived stem cells to nucleus pulposus cells

Background

Stem cell-based therapies show promise as a means of repairing the degenerate intervertebral disc, with growth factors often used alongside cells to help direct differentiation toward a nucleus pulposus (NP)-like phenotype. We previously demonstrated adipose-derived stem cell (ASC) differentiation with GDF6 as optimal for generating NP-like cells through evaluating end-stage differentiation parameters. Here we conducted a time-resolved transcriptomic characterization of ASCs response to GDF6 stimulation to understand the early drivers of differentiation to NP-like cells.

Methods

Human ASCs were treated with recombinant human GDF6 for 2, 6, and 12 h. RNA sequencing and detailed bioinformatic analysis were used to assess differential gene expression, gene ontology (GO), and transcription factor involvement during early differentiation. Quantitative polymerase chain reaction (qPCR) was used to validate RNA sequencing findings and inhibitors used to interrogate Smad and Erk signaling pathways, as well as identify primary and secondary response genes.

Results

The transcriptomic response of ASCs to GDF6 stimulation was time-resolved and highly structured, with "cell differentiation" "developmental processes," and "response to stimulus" identified as key biological process GO terms. The transcription factor ERG1 was identified as a key early response gene. Temporal cluster analysis of differentiation genes identified positive regulation NP cell differentiation, as well as inhibition of osteogenesis and adipogenesis. A role for Smad and Erk signaling in the regulation of GDF6-induced early gene expression response was observed and both primary and secondary response genes were identified.

Conclusions

This study identifies a multifactorial early gene response that contributes to lineage commitment, with the identification of a number of potentially useful early markers of differentiation of ASCs to NP cells. This detailed insight into the molecular processes in response to GDF6 stimulation of ASCs is important for the development of an efficient and efficacious cell-based therapy for intervertebral disc degeneration-associated back pain.

Endogenous Sox8 is a critical factor for timely remyelination and oligodendroglial cell repletion in the cuprizone model

Genome-wide association studies identified a single nucleotide polymorphism (SNP) downstream of the transcription factor Sox8, associated with an increased risk of multiple sclerosis (MS). Sox8 is known to influence oligodendrocyte terminal differentiation and is involved in myelin maintenance by mature oligodendrocytes. The possible link of a Sox8 related SNP and MS risk, along with the role of Sox8 in oligodendrocyte physiology prompted us to investigate its relevance during de- and remyelination using the cuprizone model. Sox8-/- mice and wildtype littermates received a cuprizone diet for 5 weeks (wk). Sox8-/- mice showed reduced motor performance and weight compared to wildtype controls. Brains were histologically analysed at the maximum of demyelination (wk 5) and on two time points during remyelination (wk 5.5 and wk 6) for oligodendroglial, astroglial, microglial and myelin markers. We identified reduced proliferation of oligodendrocyte precursor cells at wk 5 as well as reduced numbers of mature oligodendrocytes in Sox8-/- mice at wk 6. Moreover, analysis of myelin markers revealed a delay in remyelination in the Sox8-/- group, demonstrating the potential importance of Sox8 in remyelination processes. Our findings present, for the first time, compelling evidence of a significant role of Sox8 in the context of a disease model.

In mammalian skeletal muscle, phosphorylation of TOMM22 by protein kinase CSNK2/CK2 controls mitophagy

In yeast, Tom22, the central component of the TOMM (translocase of outer mitochondrial membrane) receptor complex, is responsible for the recognition and translocation of synthesized mitochondrial precursor proteins, and its protein kinase CK2-dependent phosphorylation is mandatory for TOMM complex biogenesis and proper mitochondrial protein import. In mammals, the biological function of protein kinase CSNK2/CK2 remains vastly elusive and it is unknown whether CSNK2-dependent phosphorylation of TOMM protein subunits has a similar role as that in yeast. To address this issue, we used a skeletal muscle-specific Csnk2b/Ck2β-conditional knockout (cKO) mouse model. Phenotypically, these skeletal muscle Csnk2b cKO mice showed reduced muscle strength and abnormal metabolic activity of mainly oxidative muscle fibers, which point towards mitochondrial dysfunction. Enzymatically, active muscle lysates from skeletal muscle Csnk2b cKO mice phosphorylate murine TOMM22, the mammalian ortholog of yeast Tom22, to a lower extent than lysates prepared from controls. Mechanistically, CSNK2-mediated phosphorylation of TOMM22 changes its binding affinity for mitochondrial precursor proteins. However, in contrast to yeast, mitochondrial protein import seems not to be affected in vitro using mitochondria isolated from muscles of skeletal muscle Csnk2b cKO mice. PINK1, a mitochondrial health sensor that undergoes constitutive import under physiological conditions, accumulates within skeletal muscle Csnk2b cKO fibers and labels abnormal mitochondria for removal by mitophagy as demonstrated by the appearance of mitochondria-containing autophagosomes through electron microscopy. Mitophagy can be normalized by either introduction of a phosphomimetic TOMM22 mutant in cultured myotubes, or by in vivo electroporation of phosphomimetic Tomm22 into muscles of mice. Importantly, transfection of the phosphomimetic Tomm22 mutant in muscle cells with ablated Csnk2b restored their oxygen consumption rate comparable to wild-type levels. In sum, our data show that mammalian CSNK2-dependent phosphorylation of TOMM22 is a critical switch for mitophagy and reveal CSNK2-dependent physiological implications on metabolism, muscle integrity and behavior.

Sox21 deletion in mice causes postnatal growth deficiency without physiological disruption of hypothalamic-pituitary endocrine axes

The hypothalamic-pituitary axes are the coordinating centers for multiple endocrine gland functions and physiological processes. Defects in the hypothalamus or pituitary gland can cause reduced growth and severe short stature, affecting approximately 1 in 4000 children, and a large percentage of cases of pituitary hormone deficiencies do not have an identified genetic cause. SOX21 is a protein that regulates hair, neural, and trophoblast stem cell differentiation. Mice lacking Sox21 have reduced growth, but the etiology of this growth defect has not been described. We studied the expression of Sox21 in hypothalamic-pituitary development and examined multiple endocrine axes in these mice. We find no evidence of reduced intrauterine growth, food intake, or physical activity, but there is evidence for increased energy expenditure in mutants. In addition, despite changes in pituitary hormone expression, hypothalamic-pituitary axes appear to be functional. Therefore, SOX21 variants may be a cause of non-endocrine short stature in humans.

Gene and Network Analysis of Common Variants Reveals Novel Associations in Multiple Complex Diseases

Genome-wide association (GWA) studies typically lack power to detect genotypes significantly associated with complex diseases, where different causal mutations of small effect may be present across cases. A common, tractable approach for identifying genomic elements associated with complex traits is to evaluate combinations of variants in known pathways or gene sets with shared biological function. Such gene-set analyses require the computation of gene-level P-values or gene scores; these gene scores are also useful when generating hypotheses for experimental validation. However, commonly used methods for generating GWA gene scores are computationally inefficient, biased by gene length, imprecise, or have low true positive rate (TPR) at low false positive rates (FPR), leading to erroneous hypotheses for functional validation. Here we introduce a new method, PEGASUS, for analytically calculating gene scores. PEGASUS produces gene scores with as much as 10 orders of magnitude higher numerical precision than competing methods. In simulation, PEGASUS outperforms existing methods, achieving up to 30% higher TPR when the FPR is fixed at 1%. We use gene scores from PEGASUS as input to HotNet2 to identify networks of interacting genes associated with multiple complex diseases and traits; this is the first application of HotNet2 to common variation. In ulcerative colitis and waist-hip ratio, we discover networks that include genes previously associated with these phenotypes, as well as novel candidate genes. In contrast, existing methods fail to identify these networks. We also identify networks for attention-deficit/hyperactivity disorder, in which GWA studies have yet to identify any significant SNPs.

Elevated Fra-1 expression causes severe lipodystrophy

A shift from osteoblastogenesis to adipogenesis is one of the underlying mechanisms of decreased bone mass and increased fat during aging. We now uncover a new role for the transcription factor Fra-1 in suppressing adipogenesis. Indeed, Fra1 (Fosl1) transgenic (Fra1tg) mice, which developed progressive osteosclerosis as a result of accelerated osteoblast differentiation, also developed a severe general lipodystrophy. The residual fat of these mice appeared immature and expressed lower levels of adipogenic markers, including the fatty acid transporter Cd36 and the CCAAT/enhancer binding protein Cebpa. Consequently accumulation of triglycerides and free fatty acids were detected in the serum of fasting Fra1tg mice. Fra-1 acts cell autonomously because the adipogenic differentiation of Fra1 transgenic primary osteoblasts was drastically reduced, and overexpression of Fra-1 in an adipogenic cell line blocked their differentiation into adipocytes. Strikingly, Cebpa was downregulated in the Fra-1-overexpressing cells and Fra-1 could bind to the Cebpa promoter and directly suppress its activity. Thus, our data add to the known common systemic control of fat and bone mass, a new cell-autonomous level of control of cell fate decision by which the osteogenic transcription factor Fra-1 opposes adipocyte differentiation by inhibiting C/EBPα.

SOX9 controls epithelial branching by activating RET effector genes during kidney development

Congenital abnormalities of the kidney and urinary tract are some of the most common defects detected in the unborn child. Kidney growth is controlled by the GDNF/RET signalling pathway, but the molecular events required for the activation of RET downstream targets are still poorly understood. Here we show that SOX9, a gene involved in campomelic dysplasia (CD) in humans, together with its close homologue SOX8, plays an essential role in RET signalling. Expression of SOX9 can be found from the earliest stages of renal development within the ureteric tip, the ureter mesenchyme and in a segment-specific manner during nephrogenesis. Using a tissue-specific knockout approach, we show that, in the ureteric tip, SOX8 and SOX9 are required for ureter branching, and double-knockout mutants exhibit severe kidney defects ranging from hypoplastic kidneys to renal agenesis. Further genetic analysis shows that SOX8/9 are required downstream of GDNF signalling for the activation of RET effector genes such as Sprouty1 and Etv5. At later stages of development, SOX9 is required to maintain ureteric tip identity and SOX9 ablation induces ectopic nephron formation. Taken together, our study shows that SOX9 acts at multiple steps during kidney organogenesis and identifies SOX8 and SOX9 as key factors within the RET signalling pathway. Our results also explain the aetiology of kidney hypoplasia found in a proportion of CD patients.