MondoA regulates gene expression in cholesterol biosynthesis-associated pathways required for zebrafish epiboly
Abstract
The glucose-sensing Mondo pathway regulates expression of metabolic genes in mammals. Here, we characterized its function in the zebrafish and revealed an unexpected role of this pathway in vertebrate embryonic development. We showed that knockdown of mondoa impaired the early morphogenetic movement of epiboly in zebrafish embryos and caused microtubule defects. Expression of genes in the terpenoid backbone and sterol biosynthesis pathways upstream of pregnenolone synthesis was coordinately downregulated in these embryos, including the most downregulated gene nsdhl. Loss of Nsdhl function likewise impaired epiboly, similar to MondoA loss of function. Both epiboly and microtubule defects were partially restored by pregnenolone treatment. Maternal-zygotic mutants of mondoa showed perturbed epiboly with low penetrance and compensatory changes in the expression of terpenoid/sterol/steroid metabolism genes. Collectively, our results show a novel role for MondoA in the regulation of early vertebrate development, connecting glucose, cholesterol and steroid hormone metabolism with early embryonic cell movements.
eLife digest
In most animals, a protein called MondoA closely monitors the amount of glucose in the body, as this type of sugar is the fuel required for many life processes. Glucose levels also act as a proxy for the availability of other important nutrients. Once MondoA has detected glucose molecules, it turns genetic programmes on and off depending on the needs of the cell. So far, these mechanisms have mainly been studied in adult cells.
However, recent studies have shown that proteins that monitor nutrient availability, and their associated pathways, can control early development. MondoA had not been studied in this context before, so Weger et al. decided to investigate its role in embryonic development. The experiments used embryos from zebrafish, a small freshwater fish whose early development is easily monitored and manipulated in the laboratory.
Inhibiting production of the MondoA protein in zebrafish embryos prevented them from maturing any further, stopping their development at an early key stage. This block was caused by defects in microtubules, the tubular molecules that act like a microscopic skeleton to provide structural support for cells and guide transport of cell components.
In addition, the pathway involved in the production of cholesterol and cholesterol-based hormones was far less active in embryos lacking MondoA. Treating MondoA-deficient embryos with one of these hormones corrected the microtubule defects and let the embryos progress to more advanced stages of development.
These results reveal that, during development, the glucose sensor MondoA also controls pathways involved in the creation of cholesterol and associated hormones. These new insights into the metabolic regulation of development could help to understand certain human conditions; for example, certain patients with defective cholesterol pathway genes also show developmental perturbations. In addition, the work highlights a biological link between cholesterol production and cellular responses to glucose, which Weger et al. hope could one day help to identify new cholesterol-lowering drugs.
Introduction
The glucose-sensing Mondo pathway is a key regulator of energy metabolism (Abdul-Wahed et al., 2017; Richards et al., 2017). It consists of three basic helix-loop-helix/leucine zipper (bHLH/Zip) family transcription factors (TFs): the two Mondo family factors ‘Mlx interacting protein’ (Mlxip; MondoA) and ‘Mlx interacting protein-like’ (Mlxipl; or ChREBP for ‘carbohydrate response element binding protein’) and their binding partner ‘Max-like protein X’ (Mlx) (Havula and Hietakangas, 2018; Richards et al., 2017). After their activation by glucose, likely via glucose-6-phosphate (G6P) (Dentin et al., 2012; Li et al., 2010b; Stoltzman et al., 2008), the Mondo family factors regulate gene expression as MondoA:Mlx or ChREBP:Mlx heterodimers by binding E-box type enhancer elements (‘carbohydrate response elements’, ‘ChoREs’) in the regulatory regions of their target genes (Billin et al., 2000; Jeong et al., 2011; Koo and Towle, 2000; Ma et al., 2006; Poungvarin et al., 2015; Yamashita et al., 2001).
While recent studies have highlighted a crucial role for metabolism and metabolic signaling in embryonic development (Miyazawa and Aulehla, 2018), research on Mondo pathway function has concentrated on adult animals (Hunt et al., 2015; Iizuka et al., 2004; Imamura et al., 2014; Song et al., 2019). As it is currently unknown if and how MondoA, ChREBP or Mlx contribute to metabolic regulation of embryogenesis, we investigated a potential developmental role of the Mondo pathway in zebrafish, a model for early vertebrate development (Nüsslein-Volhard, 2002). After fertilization, a cap of cells (the blastoderm) forms on top of the yolk cell, which it subsequently engulfs in a morphogenetic movement called epiboly (Bruce, 2016). Mechanistically, epiboly involves processes both in the enveloping layer (EVL) as well as the deep cell layer (DEL) of the blastoderm and in a teleost specific structure of the yolk cell, the yolk syncytial layer (YSL). By the end of epiboly, germ layers and body axes have been formed. Most developmental pathways as well as many facets of metabolism and its regulation are highly conserved between mammals and zebrafish (Gut et al., 2017; Schier and Talbot, 2005; Solnica-Krezel and Sepich, 2012). However, studies addressing Mondo pathway function in the zebrafish are still lacking.
Herein, we uncover a role of the Mondo pathway in vertebrate development. Mondo pathway gene expression and function was characterized in zebrafish embryos and cultured cells. Morpholino oligonucleotide (MO) mediated loss-of-function of MondoA severely impaired epiboly movements. Subsequent transcriptome analysis revealed genes deregulated upon loss of MondoA function and highlighted an enzyme involved in cholesterol biosynthesis, Nsdhl (NAD(P) dependent steroid dehydrogenase-like), as a potential key mediator of MondoA function in epiboly. Functional analysis of Nsdhl suggests that Nsdhl-mediated cholesterol synthesis downstream of MondoA is required for the synthesis of sufficient levels of the steroid hormone pregnenolone, which stabilizes YSL microtubules necessary for epiboly. Maternal-zygotic (MZ) mutant embryos homozygous for a small deletion allele of mondoa leading to a premature stop codon showed a severe aberrant epiboly phenotype, albeit with low penetrance that our transcriptome analysis indicates to result from compensatory changes in the expression of cholesterol and steroid biosynthesis genes.
Results
Glucose signaling by the Mondo pathway is conserved in zebrafish
For each of the Mondo pathway members MondoA, ChREBP and Mlx one single orthologue is present in the zebrafish genome (GRCz11/danRer11). We cloned the full cDNAs for zebrafish MondoA and ChREBP (GenBank Accession KF713493 [mondoa], KF713494 [chrebp]) based on the (partially) predicted sequences. Phylogenetic analysis showed that zebrafish mondoa, chrebp and mlx cluster with their mammalian and chicken homologs (Figure 1A) and revealed a high level of protein sequence conservation between zebrafish and human orthologs (Figure 1B), especially of the glucose-sensing module (GSM) specific to the Mondo family and of the DNA binding bHLH/Zip domains. This finding suggests that the functions of these proteins are conserved.
To examine whether the zebrafish Mondo pathway factors function similarly in the regulation of gene transcription as their mammalian orthologs, we studied the pathway in zebrafish PAC2 cells (Lin et al., 1994). All three Mondo pathway factors are expressed in these cells (Figure 1—figure supplement 1A). To monitor Mondo signaling, we generated a luciferase reporter gene construct driven by two ChoREs fitting the mammalian consensus (Ma et al., 2006) and a TATA box minimal promoter (Figure 1C). This construct was active in HepG2 cells, a mammalian cell culture model commonly used for Mondo pathway studies (Kim et al., 1996; Yu and Luo, 2009; Figure 1—figure supplement 1B). In PAC2 cells, bioluminescence equally increased in a dose-dependent manner upon glucose treatment (Figure 1C). No significant changes in bioluminescence levels were detected upon glucose treatment in cells expressing a constitutively active luciferase reporter (pGL3-Control; Figure 1D), excluding a general unspecific increase in transcriptional activity by glucose treatment.
We next tested whether overexpression of Mondo pathway members enhances glucose induced pathway activity. Transient overexpression of MondoA and Mlx activated transcription from the reporter also under low glucose conditions (Figure 1E), as shown by the increased bioluminescence compared with transfection of the 2xChoRE reporter alone (p≤0.001). Overexpression of the MondoA and Mlx factors together led to an even more pronounced effect on bioluminescence (p≤0.01), revealing synergistic effects. High glucose levels caused significant reporter gene induction in control cells (p≤0.01; Figure 1E). Importantly, strong glucose induction of the reporter was also shown by cells overexpressing either Mlx (p≤0.01) or MondoA (p≤0.001) alone as well as both factors together (p≤0.001; Figure 1E). An unrelated control reporter construct (pGRE-Luc, Weger et al., 2012) was not responsive to any of the treatments (Figure 1—figure supplement 1C). Together, the data indicate that limited availability of endogenous pathway components restrains pathway activity, and that both basal activity and the response to higher glucose levels are potentiated when more sensor proteins are available.
In addition, we explored the effect of morpholino oligonucleotide (MO) mediated loss-of-function of mondoa to confirm that Mondo pathway function is required for glucose induction of reporter expression. Cells transfected with a control 5 bp mismatch MO (mondoa-mis) showed a 2.0-fold induction of bioluminescence from the 2xChoRE reporter plasmid by high glucose levels, while in cells transfected with a MO directed against the translation start site (mondoa-mo) this glucose response was abolished (Figure 1F). Taken together, our data demonstrate that ChoRE mediated glucose induction of gene expression is regulated by the Mondo pathway also in zebrafish.
The zebrafish Mondo pathway is present in early embryos
After confirming the similarity of zebrafish Mondo pathway function to mammals in cultured cells, we studied its role during development. We reanalyzed a previously published zebrafish developmental transcriptome dataset (White et al., 2017, ENA accession number ERP014517) and detected maternal transcripts of mondoa, chrebp and mlx as well as expression of the three genes throughout development, following distinct temporal patterns (Figure 1G–I). Whole mount in situ hybridization analysis (WISH) for mondoa revealed a ubiquitous expression pattern (Figure 1J). Similar results were observed for chrebp and mlx (Figure 1—figure supplement 1D,E). At the 50% epiboly stage, all three genes are expressed throughout the embryo (Figure 1K; Figure 1—figure supplement 1F).
Next, we tested whether glucose treatment increases the expression of mammalian Mondo pathway target gene homologs in early zebrafish embryos, namely hexokinase 2 (hk2) (Sans et al., 2006), fatty acid synthase (fasn) (Ma et al., 2006) and thioredoxin-interacting protein a (txnipa) (Stoltzman et al., 2008). We injected the glucose analogue 2-deoxy-glucose (2-DG), which is metabolized to the G6P analogue 2-DG6P but not further (Chi et al., 1987; Stoltzman et al., 2008), thereby avoiding activation of other pathways relying on downstream metabolization of glucose. This treatment significantly increased expression of hk2 (p<0.001), fasn (p<0.01) and txnipa (p<0.001; Figure 1L–N), but not of the control gene eef1a1l1/ef1a (eukaryotic translation elongation factor 1 alpha 1 (Figure 1O), ruling out a nonspecific general increase in gene expression by 2-DG treatment. Furthermore, txnipa expression was not induced by 2-DG in embryos injected with mondoa-mo, showing that MondoA function is required for its induction (Figure 1—figure supplement 1G). To directly examine if 2-DG treatment regulates transcription by ChoRE enhancer elements in the early embryo, we injected the 2xChoRE luciferase reporter construct together with 2-DG. Embryos treated with 2-DG (n = 72) showed a 2.3-fold increase (p≤0.01) in bioluminescence over the control (n = 78; Figure 1P). Together, these results suggest that Mondo pathway mediated glucose signaling is present in embryos as early as at the sphere stage.
Loss of mondoa function leads to severe epiboly defects
To examine Mondo pathway function during development, we injected embryos with MOs targeting mondoa and chrebp function. While injections of a MO directed against chrebp did not result in an aberrant phenotype (Figure 2—figure supplement 1A,B), mondoa-mo injected embryos showed a striking delay in epiboly movements when compared with uninjected or mondoa-mis injected embryos (Figure 2A–D; Video 1). This phenotype was dose-dependent, with the strongest delay occurring at the highest concentration (Figure 2A).
To test for efficiency and specificity of the MO-mediated knockdown, a number of controls were employed. GFP expression from a reporter construct carrying the target sequence was efficiently abolished by co-injection of mondoa-mo, but not mondoa-mis (Figure 2—figure supplement 1C). An involvement of P53-mediated apoptosis in the observed phenotypes was ruled out by co-injection of mondoa-mo and a MO directed against p53 (p53-mo), leading to the same phenotype as observed upon injection with mondoa-mo alone (Figure 2—figure supplement 1D). Knockdown of the mondoa cofactor mlx with a MO directed against the translation start site of mlx (mlx-mo) showed a dose-dependent impairment of epiboly similar to the one observed with mondoa-mo, albeit somewhat less severe (Figure 2C,D). These results are consistent with the model that MondoA forms heterodimers with Mlx to exert its function also in epiboly. Finally, to obtain additional evidence that the observed effects are indeed dependent on MondoA function, we attempted to rescue epiboly movements in mondoa-mo injected embryos by co-injecting a mondoa mRNA in which the sequence targeted by mondoa-mo was mutated. Epiboly progression was monitored with WISH against the blastoderm margin marker no-tail (ntl). Embryos injected with only mondoa-mo reached 38% epiboly, whereas embryos co-injected with mondoa mRNA arrived at 47% epiboly (Figure 2E). These values are consistent with those reported in rescue experiments addressing other genes implicated in epiboly (Hsu et al., 2006; Schwend et al., 2011) and might be explained by decreased mRNA stability and/or decreased translation efficiency of the rescue construct in vivo. In conclusion, these observations indicate that mondoa-mo specifically targets mondoa gene function, and reveal that loss of MondoA function severely impacts on epiboly during zebrafish development.
Loss of mondoa function specifically affects YSL epiboly
Next, we aimed at clarifying which epiboly-related processes are perturbed by loss of mondoa function. Since the Mondo pathway has been implicated in the regulation of energy metabolism in adult mammalian tissues, we tested if the aberrant epiboly phenotype of mondoa morphants is a consequence of generally impaired energy metabolism. Our results showed that energy charge was unaltered in the morphants (Figure 3—figure supplement 1A). Alternatively, abnormal cell proliferation or defects in migratory behavior might limit the spreading capacity of the blastoderm. To count cell numbers in the embryos by automated quantification, we imaged transgenic embryos carrying a H2A transgene which labels cell nuclei (Tg(h2afva:h2afva-GFP)) (Pauls et al., 2001). We employed digital scanned laser light sheet microscopy (DSLM) (Kobitski et al., 2015), allowing the optical sectioning of the embryo with a high spatial and temporal resolution. As shown by this analysis, mondoa morphants had similar numbers of cells as control embryos (Figure 3A). Accordingly, the cell density was higher in mondoa morphants compared to mismatch controls, as movement to the vegetal pole was impaired while cells continued to divide (Figure 3B,C; Videos 2–4). These observations strongly argue against deficient cell proliferation. In addition, we observed germ ring-like thickenings (Figure 2A,C; arrowheads) as well as internalization of cells at the blastoderm margin in mondoa morphants (Figure 3D; Figure 3—figure supplement 1B–D; Videos 5 and 6). Furthermore, the average migration speed of deep cells was identical for morphants and controls, demonstrating normal individual cell migration capacity in morphants (Figure 3—figure supplement 1E).
EVL, DEL and YSL have all been implicated in mechanisms of epiboly (Bruce, 2016). Thus, we turned to examining how the three different embryonic cell layers were affected in mondoa morphants. Imaging of Tg(h2afva:h2afva-GFP) embryos allowed us to assign nuclei to the EVL or to the DEL and to distinguish the YSL nuclei. Our observations revealed that EVL cells continued to progress and could reach a position vegetally to the YSL nuclei in the morphants (Figure 4—figure supplement 1A–F, Video 7). Consistently, EVL integrity as revealed by rhodamine-phalloidin staining (Lepage and Bruce, 2010; Lepage et al., 2014) was intact in the morphants (Figure 4—figure supplement 1G–H). These findings indicate that EVL epiboly still occurs in the morphants.
However, deep cells of mondoa morphants lagged behind both EVL and YSL layers, consistent with models suggesting rather passive movement of these cells into spaces left open by the other two layers (Bensch et al., 2013; Song et al., 2013; Figure 4A,B; Figure 4—figure supplement 1A–F). Nevertheless, in mondoa morphants a few deep cells were detected more vegetally than the YSL nuclei. Also, YSL nuclei showed a much broader and more disorganized distribution than in control injected or wild-type embryos (Figure 4—figure supplement 1A–F). The abnormal shape of the YSL was already evident at a slightly earlier stage (compare Figure 4A with Figure 4B). YSL nuclei appeared to cluster and to break away from the blastoderm margin (Figure 4B, white arrow), in contrast to the more even distribution in a relatively narrow band around the margin in wild-type embryos. These findings pointed to a defect in YSL structure or function as a potential reason for the impaired epiboly in mondoa morphants.
To test whether MondoA function in the YSL itself is important for epiboly, we injected mondoa-mo at the 1 k-cell stage into the yolk cell only (Figure 4C). At this late time of YSL injection, additional mondoa transcript and protein have been produced after mid-blastula transition (MBT), which may limit knockdown efficiency compared with injection into the zygote. Strikingly, mondoa knockdown limited to the YSL still led to epiboly delay and developmental arrest compared with uninjected or control injected embryos (Figure 4D). Thus, MondoA function in the YSL indeed is important for correct epiboly movements. Impaired YSL structure or function may lead to impaired vesicle trafficking at the YSL-EVL border (Lepage and Bruce, 2010; Solnica-Krezel et al., 1995; but see also Lepage and Bruce, 2014) or to defective patterning of the overlying blastoderm (Carvalho and Heisenberg, 2010), both of which can perturb epiboly. However, we did not find any evidence for major perturbations of marginal endocytosis or blastoderm patterning in the morphants (Figure 4—figure supplement 1J–Q).
Transcriptome analysis points to the cholesterol synthesis enzyme nsdhl as a main target gene of MondoA in the early embryo
To obtain cues how MondoA function might impact on epiboly we carried out an unbiased examination of differential gene expression upon mondoa knockdown by applying RNA-seq. Embryos were injected either with mondoa-mo or mondoa-mis and left to develop until the control embryos reached sphere stage, just before the morphological phenotype becomes apparent (Figure 5—figure supplement 1A–C; Supplementary file 1). We validated a random subset of the genes expressed differentially between the two conditions by RT-qPCR, confirming the RNA-seq data (Figure 5—figure supplement 1D,E). Interestingly, the gene most down-regulated by mondoa-knockdown was nsdhl (NAD(P) dependent steroid dehydrogenase-like; Figure 5A). Nsdhl is part of the cholesterol biogenesis pathway, catalyzing steps in the conversion of the cholesterol precursor lanosterol to zymosterol (Figure 5B). Consistently, a global metabolic pathway gene set enrichment analysis revealed that both the ‘sterol biosynthesis’ pathway, of which Nsdhl is part, and the upstream ‘terpenoid backbone biosynthesis’ pathway are coordinately downregulated in the morphants (Supplementary file 2, Figure 5B, Figure 7—figure supplement 1A,B).
nsdhl expression is regulated by the Mondo pathway in a glucose-dependent manner and loss of its function perturbs epiboly
Cholesterol has an important role as substrate for the synthesis of pregnenolone, the first steroid hormone within the steroid hormone biosynthesis pathway. Pregnenolone and other steroids have been implicated in the regulation of zebrafish epiboly (Eckerle et al., 2018; Hsu et al., 2006; Schwend et al., 2011; Weng et al., 2013). Therefore, we hypothesized that mondoa knockdown might affect epiboly by interfering with Nsdhl-dependent cholesterol and steroid hormone biogenesis. To test this hypothesis, we first validated regulation of nsdhl expression by MondoA mediated glucose sensing. We observed a significant increase of basal nsdhl transcript levels upon 2-DG injection (p≤0.05) and upon overexpression of both MondoA and Mlx (p≤0.001) in early embryos, which was augmented further when both treatments were combined (p≤0.001; Figure 5C). No changes in expression were observed for ef1a under these conditions (Figure 5D), confirming the target gene specific effects of the treatments. Together, these results strongly suggest that the Mondo pathway mediates glucose regulation of nsdhl expression, thereby linking glucose sensing with cholesterol synthesis.
We next examined whether nsdhl knockdown can affect epiboly using a MO targeted against a splice site of the nsdhl gene (nsdhl-mo). We observed a highly similar phenotype to the mondoa knockdown phenotype, with nsdhl-mo injected embryos arrested at 50% epiboly when control embryos had completed epiboly (Figure 5E,F). As mondoa, nsdhl is ubiquitously expressed in the embryo, including the YSL (Figure 5—figure supplement 1F,G), and injection of nsdhl-mo into the yolk cell at the 1 k-cell stage resulted in an epiboly delay and developmental arrest as observed for YSL specific mondoa knockdown (Figure 5G). Thus, the function of Nsdhl in the biosynthesis pathway to cholesterol in the YSL appears to be essential for epiboly movements.
Pregnenolone partially rescues epiboly in mondoa morphants and normalizes microtubule patterns
Altered cholesterol biosynthesis in the YSL by loss of MondoA/Nsdhl function might affect steroid hormone-dependent epiboly mechanisms. Therefore, we next tested whether pregnenolone treatment can rescue mondoa morphants. Treatment with pregnenolone indeed partially rescued epiboly in mondoa-mo injected embryos (Figure 5H) at a level comparable to the rescue with mondoa mRNA. Accordingly, the aberrant epiboly phenotype of nsdhl morphants was also partially rescued by pregnenolone treatment (Figure 5I), confirming its role downstream of the MondoA/Nsdhl pathway.
Pregnenolone was shown to regulate epiboly by mediating microtubule stability (Hsu et al., 2006; Schwend et al., 2011; Weng et al., 2013). To examine whether yolk cell microtubule integrity is affected in mondoa morphants, we stained embryos against α-tubulin. In uninjected embryos (n = 9/9), the YSL forms a regular band below the blastoderm, with yolk syncytial layer nuclei (YSN) surrounded by microtubule organizing centers (MTOCs) (Solnica-Krezel and Driever, 1994; Strähle and Jesuthasan, 1993; Figure 6A). Below this area, dense arrays of microtubules spread from the MTOCs along the animal-vegetal axis into the yolk cytoplasmic layer (YCL), forming root-like bundles. Compared with the control embryos (Figure 6C,D), mondoa-mo injected embryos showed a rather unstructured YSL with a less well defined border (Figure 6E,F). YSN were less regularly arranged, and their associated MTOCs formed stellar structures (n = 14/16; Figure 6B), apparently reflecting shortened microtubules that did not form the arrays of root-like bundles along the animal-vegetal axis seen in uninjected controls. Importantly, the YSL of morphants treated with pregnenolone appeared to be better organized compared to the untreated morphants because the YSN and their corresponding MTOCs were arranged as in the uninjected control (n = 6/9; Figure 6B,G,H). Furthermore, in pregnenolone treated embryos, the microtubule arrays were longer and formed array-like structures along the animal-vegetal axis of the embryos similar to the control. Taken together, these results strongly suggest that MondoA functions in epiboly via the stabilization of yolk microtubules by pregnenolone formed downstream of cholesterol biosynthesis.
Maternal-zygotic mutants of mondoa show partially penetrant aberrant epiboly phenotypes
Finally, we wanted to validate the knockdown results with a genetic loss-of-function model of mondoa using CRISPR-Cas9. We identified a mutant allele transmitted through the germline that carried a 5 bp deletion in exon 3 (Figure 7A), leading to a frameshift in the coding sequence and a predicted premature stop eliminating both the DNA binding domain and most of the glucose-sensing domain. Embryos homozygous for this allele did not show an aberrant embryonic phenotype and were raised to sexual maturity (Figure 7B,C). By contrast, we found a few maternal-zygotic (MZ) mutants from incrosses of homozygous mutant parents that showed a severe epiboly delay and an arrest in mid-epiboly (Figure 7B,C).
The low penetrance of the epiboly phenotype in the MZ mutants is indicative of compensatory mechanisms that counteract the loss of function of mondoa. To test this assumption, we injected mondoa-mo into MZmondoa mutants. Indeed, we observed a resistance of the mutants to the perturbance of epiboly caused by this MO in wild-types (Figure 7D). This observation strongly suggests the presence of strong compensatory or buffering mechanisms that enable embryonic development even when MondoA function is genetically perturbed. It also further confirms the specificity of the morphant phenotype.
MZmondoa mutants show compensatory changes in expression of cholesterol/steroid biosynthesis genes
To begin to explore potential compensatory mechanisms allowing epiboly progression in MZmondoa mutants, we performed total RNA sequencing of MZmondoa embryos with severe and no epiboly anomalies and of wild-type embryos to determine differential gene expression (Supplementary file 1). To compare the gene expression signatures of MZmondoa mutants with and without aberrant epiboly phenotype to those of mondoa morphants and their respective controls, we employed the rank-rank hypergeometric overlap (RRHO) algorithm (Plaisier et al., 2010). The analysis revealed that differential gene expression significantly overlapped between mutants and morphants, demonstrating that both conditions reflect a lack of MondoA function. Specifically, mutant embryos showed an overlap in differential gene expression with respect to the morphants mainly in the upregulated gene fraction. There were also some differences apparent between the mutants with and without an aberrant epiboly phenotype. We observed that the overlap in downregulated genes was weaker between morphants and mutants without aberrant epiboly phenotype (compare lower left quadrant in Figure 7—figure supplement 2B and C). This observation indicated that mainly the downregulated genes in mondoa morphants were transcriptionally compensated in MZmondoa mutants with unaffected epiboly. In contrast, the affected mutants showed a weaker overlap in the upregulated gene signatures, suggesting that only a fraction of the upregulated genes in MZmondoa mutants was actually important for phenotype compensation. We next explored these differential gene expression patterns in more detail.
As paralogues of genes have been implicated in compensation (El-Brolosy et al., 2019; Rossi et al., 2015), we looked at expression of other Mondo pathway components. Neither chrebp, the paralogue of mondoa, nor the MondoA partner mlx were significantly upregulated in MZmondoa mutants, even though trends (p=0.06) for elevated mRNA levels were seen for mlx in both mutant phenotypes and for chrebp in the affected phenotype only (Figure 7E).
We next turned our attention to genes in the cholesterol biosynthesis pathway based on our previous observations in morphant embryos. Our analysis showed that nsdhl gene expression was comparable to wild-type embryos in both MZmondoa mutant phenotypes (affected and unaffected epiboly) and, thus, was compensated in both (Supplementary file 1 and Figure 7—figure supplement 2E–F). Therefore, differences in nsdhl expression between MZmondoa mutants cannot explain the differences in the epiboly phenotype. This observation indicates that the restoration of expression to wild-type levels of this most downregulated gene upon mondoa knockdown by itself was not sufficient for compensation. Consistently, in both mutant phenotypes the entire sterol biosynthesis pathway of which nsdhl forms part (Figure 5B) lacked the coordinated downregulation that was seen in the morphants (unaffected: p=0.87, affected: p=0.79; Figure 7F and Figure 7—figure supplement 1A).
But what is causing the epiboly arrest if MZmondoa mutants with affected epiboly share such compensatory transcriptional signatures with MZmondoa mutants without aberrant epiboly phenotype? mondoa morphants showed a significant downregulation in terpenoid backbone biosynthesis (p=0.016) as well as in sterol biosynthesis pathways (p=0.022; Supplementary file 2). By contrast, MZmondoa mutants with unaffected epiboly showed a coordinated upregulation of terpenoid pathway genes (i.e., they showed a statistically significant enrichment for upregulation in a directional gene set enrichment test; p=0.012; Figure 7F–G, Figure 7—figure supplement 1B and Supplementary file 2). Intriguingly, mutants with affected epiboly lacked this significant coordinated upregulation in the terpenoid backbone biosynthesis pathway that was present in unaffected mutants (p=0.178; Figure 7F–G and Figure 7—figure supplement 1B). This key difference highlights the upregulation of gene expression in the terpenoid backbone biosynthesis pathway as a compensation response specifically required to restore the normal epiboly phenotype.
We additionally observed that gene expression of steroid hormone biosynthesis pathway genes showed a coordinated upregulation in MZmondoa mutants with unaffected epiboly (p=0.033; Figure 7F and Figure 7—figure supplement 1C), probably to compensate for reduced flux through the biosynthesis pathways upstream. This upregulation was also shown by MZmondoa mutants with aberrant epiboly phenotype (p=0.017; Figure 7F and Figure 7—figure supplement 1C). Strikingly, it was even increased when compared to the mutants with unaffected epiboly (p=0.049), and thereby led to strongly increased expression levels compared to wild-type. This is illustrated by the expression of cyp11a1, the enzyme that is immediately upstream of pregnenolone synthesis (Supplementary file 2, Figure 7—figure supplement 2G). However, in the absence of coordinated regulation of genes in the terpenoid pathway, this regulation was apparently not sufficient to rescue epiboly progression.
Taken together, despite some successful compensatory gene expression regulations, MZmondoa mutants with affected epiboly fail to properly regulate the entire set of genes necessary for epiboly progression. They are not lacking compensation per se, but their compensation is inadequate and fails.
In summary, results from both MO mediated and genetic loss-of-function of MondoA revealed MondoA as a key developmental regulator of progression through epiboly by regulating expression of cholesterol/steroidogenesis pathway genes important for yolk cell microtubule function.
Discussion
Here, we have characterized the Mondo pathway in zebrafish and revealed a developmental function for Mondo signaling, implying a novel role for this glucose-sensing pathway in addition to its functions in adult animals (Song et al., 2019). Previously, studies in non-vertebrates have implicated homologs of Mondo pathway factors in later stages of embryonic or larval development: Loss of function of the Mlx orthologue in Drosophila is lethal at the late pupal stage (Havula et al., 2013), while knockdown of a distant Mlx homolog in C. elegans affected movements of ray one precursor cells (Pickett et al., 2007), which will form a subset of male-specific sense organs in the tail. To the best of our knowledge, ours is the first study demonstrating conserved Mondo signaling mechanisms between zebrafish and mammals and linking its function to early embryonic development of vertebrates.
Our findings demonstrate that MondoA acts upstream of the synthesis of pregnenolone, an important regulator of yolk cell microtubule stability and epiboly movements (Hsu et al., 2006; Schwend et al., 2011; Weng et al., 2013), by regulating the expression of a whole set of genes in the terpenoid/sterol biosynthesis pathway, as exemplified by the cholesterol synthesis enzyme nsdhl. It is tempting to speculate that the role of MondoA in cholesterol/steroid metabolism we here uncover in the zebrafish is a feature conserved across vertebrates, albeit not much is known in mammals. One in vitro study described ‘sterol biosynthetic process’ as one category of genes upregulated in MondoA deficient HeLa cells under conditions of acidosis (Wilde et al., 2019). For the MondoA paralogue Chrebp, a study in vivo reported that ‘cholesterol biosynthesis’ genes are upregulated in the liver of Chrebp knockout mice which were fed with a high fructose diet (Zhang et al., 2017). While these studies were performed under specific disease conditions (i.e., a human cancer cell line under acidosis conditions or animals fed with a high fructose diet), they support our observations that the Mondo pathway plays an important role in the regulation of cholesterol/steroid metabolism also under non-disease conditions. Furthermore, these observations indicate gene expression changes due to compensatory mechanisms upon loss of Mondo pathway function, as observed in our study and discussed in the following section.
The low penetrance of epiboly phenotypes in MZmondoa mutant zebrafish embryos indicates the presence of efficient compensatory pathways, triggered by the genetic lesion but not by antisense-mediated knockdown, that may be operating in mammals as well (El-Brolosy and Stainier, 2017). Nonsense mediated decay products of the mutated transcript are proposed to trigger transcriptional upregulation of homologous genes that compensate for the mutant gene’s function (El-Brolosy et al., 2019). However, chrebp and mlx mRNA levels were not significantly changed compared to wild-type in mutants with unaffected epiboly, implicating the presence of alternative compensatory mechanisms. Interestingly, many terpenoid, sterol and steroid biogenesis genes showed lower than wild-type levels of expression in the morphants and equal or higher levels in the mutants with unaffected epiboly, suggesting that increased steroidogenesis activity rescues epiboly in a majority of MZ mutant embryos. In the remaining mutant embryos, the lack of coordinated upregulation of terpenoid backbone biosynthesis genes is indicative of an incomplete compensation attempt. Apparently, concerted changes across the entire cholesterol/steroid biosynthesis pathway are needed to rescue epiboly progression in the absence of MondoA function. The mere restoration of nsdhl expression levels, for example, appears not to be sufficient for such a rescue, perhaps because a few genes in the sterol biosynthesis pathway still show deregulation of expression in MZ mutant embryos of both phenotypes (Figure 7—figure supplement 1A). Therefore, even though globally there is no coherent up- or downregulation in sterol biosynthesis gene expression of mutants compared to wild-type embryos based on gene set enrichment analysis, pathway function in mutants may still be suboptimal and require the compensatory upregulation in the upstream terpenoid pathway to achieve sufficient synthesis flux. Furthermore, random variation in the amount of material deposited by the mother or in protein expression of key genes may limit compensation in the mutants with affected epiboly. A more detailed understanding of these processes awaits the determination of metabolite fluxes in morphants, mutants and wild-type embryos, specifically in the YSL. Importantly, incomplete penetrance is a frequent phenomenon in human genetic diseases and has been attributed to a wealth of mechanisms, which have been suggested as targets for therapy (Cooper et al., 2013). Future studies how disturbed MondoA function is buffered during development might open leads to manipulating dysregulated Mondo pathway function in adult metabolic regulation.
We observed that mondoa knockdown in the YSL was sufficient to cause an epiboly phenotype, pointing towards a crucial role of MondoA in this tissue. Even though the YSL and epiboly movements are specific for teleosts, there are several indications that the mechanisms detected in the zebrafish are relevant in mammals as well. For example, steroidogenic enzymes are expressed in extraembryonic tissues in both fish and mice, the fish YSL (Hsu et al., 2006; Hsu et al., 2009) and the giant trophoblast cells of the murine placenta (Arensburg et al., 1999). Indeed, it has been proposed that zebrafish epiboly is analogous to the spreading of the trophectoderm during implantation (Kane and Adams, 2002). Consistent with this idea, epiboly movements are severely impaired upon knockdown of the zebrafish homolog of solute carrier family 3 member 2 (Slc3a2) (Takesono et al., 2012), a factor crucial for trophoblast cell adhesion, migration and fusion in mammals (Kabir-Salmani et al., 2008; Kudo et al., 2003). Thus, it is conceivable that developmental functions of MondoA are conserved among vertebrates. Testing this hypothesis awaits a close examination of early embryonic and placental phenotypes in MondoA mutant mice, the analysis of which has so far been limited to adult phenotypes (Ahn et al., 2019; Imamura et al., 2014), and of potential functions of maternal MondoA mRNA or protein in mammals (Li et al., 2010a). Given the small litter size of rodents, phenotypes with low penetrance may easily be overlooked and are difficult to study in these animals.
Importantly, mutations in the NSDHL gene in humans have been linked to an X-linked dominant condition called CHILD syndrome (‘congenital hemidysplasia with ichthyosiform erythroderma and limb defects’, reviewed in Herman, 2000) that is usually lethal in males (Bornholdt et al., 2005). Murine male embryos carrying Nsdhl mutations die shortly after implantation, and those carrying moderate or mild alleles showed reduced placental thickness with a poorly vascularized fetal labyrinth (Caldas et al., 2005; Liu et al., 1999 and references therein). Interestingly, also in female mouse embryos carrying a maternal null allele placental area was severely reduced (Cunningham et al., 2010). In these embryos, paternal X chromosome inactivation in yolk sac endoderm and trophoblast-derived lineages causes complete loss of function only in these tissues. Given the potential analogous functions of placenta and YSL (Carvalho and Heisenberg, 2010), our results regarding Nsdhl function in epiboly make it tempting to speculate that Nsdhl contributes to the development of the mammalian placenta by regulating microtubule stability. Remarkably, stability of microtubules apparently plays a role in trophoblast differentiation and implantation based on in vitro studies (Bates and Kidder, 1984; Douglas and King, 1993).
In summary, our results identify a novel role for the Mondo pathway in vertebrate development and link glucose sensing to cholesterol/steroid biogenesis. These results contribute to the growing recognition of metabolic regulatory functions in development (Miyazawa and Aulehla, 2018) and broaden our understanding of Mondo signaling.
Materials and methods
Fish stocks and maintenance
Request a detailed protocolAll zebrafish husbandry was performed in accordance with the German animal protection standards and approved by the Government of Baden-Württemberg, Regierungspräsidium Karlsruhe, Germany (Aktenzeichen35-9185.64/BH KIT), by the Home Office, United Kingdom (Scientific Procedures, Act 1986), and by the Service de la consommation et des affaires vétérinaires du canton de Vaud, Switzerland. Wild-type fish were descendants of the AB strain (University of Oregon, Eugene) or the golden strain and have been reared for several years in the laboratory. Fish were bred and raised in E3 medium (Nüsslein-Volhard, 2002) (5 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl2, 0.33 mM MgSO4, 0.1% methylene blue). Staging of embryos was performed on a dissecting microscope according to Kimmel et al., 1995.
Cloning of cDNAs and other plasmid constructs
Cloning of cDNAs
Request a detailed protocolIn order to obtain a full length cDNA sequence for chrebp, we performed a PCR with the primer pairs PCR_chrebp_1 (fw: 5’-aaccgcagaaccgaatcag-3’, rv: 5’-tagctgtgcgatgaatgtcc-3’) on a commercially available cDNA template (Life Technologies, 1 cell stage). 5’ and 3’ ends of the cDNA were determined with 5’ and 3’ RACE using the FirstChoice RLM-RACE Kit as recommended by the manufacturer (Ambion). A 5’ RACE was also performed to confirm the translation start site of mondoa.
The open reading frames (ORFs) of mondoa and mlx were cloned by nested PCR using Platinum Taq DNA polymerase as recommended by the manufacturer (Life Technologies). The first PCR was carried out with the primers PCR_mondoa_1 (fw: 5’-ggacgaaggaattgacgaaa-3’, rv: 5’-ctaagaccctgaggaatgtt-3’) and PCR_mlx_1 (fw: 5’-attgtttattccaaaatgac-3’, rv: 5’-gtgctgagccctgctgtatt-3’) for mondoa and mlx, respectively. The cDNA template used was a commercially available cDNA (Life Technologies, 1 cell and 18–30-somite stages). PCR products obtained were used as templates for the second PCR round with the primers PCR_mondoa_2 (fw: 5’-atggctactcgacaataccgga-3’, rv: 5’-ccctgaggaatgttttttgggg-3’) and PCR_mlx_2 (fw: 5’-atgacggaaaacagcgcgt-3’, rv: 5’-cagttgaccattcacctgctgc-3’) for mondoa and mlx. The final PCR products were subcloned into the pGEMT-Easy vector and sequenced.
Expression vectors for mondoa and mlx
Request a detailed protocolFor subcloning into an expression vector to perform overexpression studies, a further PCR step was carried out with the Platinum Taq DNA polymerase kit, using pGEMT-Easy:MondoA and pGEMT-Easy:Mlx as templates and the primers PCR_mondoa_3 (fw: 5’-taggtacccatggctactcgacaataccgga-3’, rv: 5’-tagctcgagccctgaggaatgttttttgggg-3’) and PCR_mlx_3 (fw: 5’-taggtacccatgacggaaaacagcgcgt-3’, rv: 5’-tagctcgagcagttgaccattcacctgctgc-3’), which include XhoI and KpnI overhangs. The PCR products were cloned into the XhoI and KpnI restriction sites of a pcDNA3.1(–) vector (Life Technologies) which lacked the neomycin cassette (removed with an NaeI digest), resulting in the pcDNA3.1(–)ΔNeo:MondoA and pcDNA3.1(–)ΔNeo:Mlx vectors, respectively. For the mondoa mRNA rescue experiments, silent mutations were created in the pcDNA3.1(–) ΔNeo:MondoA vector so that the protein sequence was retained, while mondoa-mo binding was abolished. Mutagenesis was performed with the QuikChange site-directed mutagenesis kit as described by the manufacturer (Stratagene). The primer used was PCR_mondoa_res (5’-ggctagttaagcttggtacccatggcaacgaggcagtataggagggccacaatgatgatcaagca-3’). RNA was transcribed using the mMESSAGE mMACHINE Kit (T7) (Ambion) on templates linearized with DraIII.
Luciferase reporter constructs
Request a detailed protocolThe constitutively expressed pGL3-Control luciferase vector was obtained from Promega (#E1741). The PathDetect pGRE-Luc cis-Reporter plasmid was obtained from Agilent Technologies (#240133). The reporter constructs pLucMCS:2xChoRE and pT2Luci:2xChoRE were generated by cloning the following annealed 5’ phosphate modified oligonucleotides into the XhoI site of the pLucMCS and pT2Luci:MCS recipient plasmids, respectively, as previously described (Weger et al., 2011; Weger et al., 2012; Weger et al., 2013b): 2xChoRE oligos: fw: 5’-tcgatcacgagtgggccgcgtgccggccacgagtgggccgcgtga-3’;rv: 5’-tcgatcacgcggcccactcgtggccggcacgcggcccactcgtga-3’.
pCS2:GFP:MondoA and pCS2:GFP:ChREBP knockdown efficiency constructs
Request a detailed protocolTo test the knockdown efficiency of the mondoa-mo, the ATG region of mondoa was fused to a GFP reporter gene. The ATG region of mondoa (+36/–20) was generated by annealing of the following oligonucleotides containing NcoI and EcoRI overhangs: MondoA_1: fw: 5’-aattatacttatagtaaaagtcgaagcg-3’, rv: 5’-gacgcgcttcgacttttactataagtat-3’; MondoA_2: fw: 5’-cgtatttttaacatggctactcgacaataccg-3’, rv: 5’-catgcggtattgtcgagtagccatgttaaaaa-3’. These oligonucleotides were then inserted in frame in front of the GFP of a pCS2:GFP vector (Ertzer et al., 2007) linearized with NcoI and EcoRI, yielding the pCS2:GFP:MondoA construct. This construct was linearized with NotI and transcribed with the mMESSAGE mMACHINE (Sp6) kit as recommended by the manufacturer (Ambion). The ATG region of chrebp (+26/–34) was generated similarly by annealing the following oligos: chrebp_1: fw: 5'-aattgcagtttgtgaagagtgaataacactag-3', rv: 5'-agaactagtgttattcactcttcacaaactgc-3'; chrebp_2: fw: 5'-ttctcttgcgatggcaaaagcagagtattcagaaac-3', rv: 5'-catggtttctgaatactctgcttttgccatcgcaag-3'.
Templates for digoxigenin labeled probes
Request a detailed protocolIn order to generate the templates for WISH probes directed against the chrebp, mondoa, mlx and nsdhl, PCR was carried out with the following oligos: chrebp probe 1: fw: 5’-ttcatctgcgctcaacactc-3’, rv: 5’-ttgagatctggacgagagca-3’; chrebp probe 2: 5’-ctacagccagagggaagcac-3’, rv: 5’-gatccgtcgtgtgtctgttc-3’; mondoa probe 1: fw: 5’-actaccgtcagccacaggtc-3’, rv: 5’-cactctgcgactctgactgtt-3’; mondoa probe 2: fw: 5’-ttccttctctgcctctggac-3’, rv: 5’-ccagcagaaccacaaggaat-3’; mlx probe 1: fw: 5’-gatggcaccttcagtgacaa-3’, rv: 5’-tggcttgcagtgttcctcta-3’; mlx probe 2: fw: 5’-tccacacttaggaaggaggtg-3’, rv: 5’-tgattggaatcactgccata-3’; nsdhl probe 1: fw: 5’-agcgcttagtacaccgagga-3’, rv: 5’-cgagcgcagctgtagtaatg-3’; nsdhl probe 2: fw: 5’-tcatacaggcttgccatgag-3’, rv: 5’-tcctaggctcgtcttgcagt-3’. PCR products were subcloned into the pGEMT-Easy vector (Promega). The no tail (ntl), goosecoid (gsc), ventral homeobox (vox), ventral expressed homeobox (vent) and casanova/SRY box containing gene 32 (sox32) probes were isolated in a large scale transcription factor screen (Armant et al., 2013). Digoxigenin labeled probe synthesis was done following the manufacturer’s instructions (DIG RNA labeling mix, Roche).
RNA extraction, cDNA Synthesis, RT-qPCR
Request a detailed protocolMZmondoa and wild-type RNA was extracted using RNA Clean and Concentrator-25 columns with an included DNase I treatment (ZYMO RESEARCH) for the corresponding RNA-seq experiment, as described in the manufacturer’s protocol. RNA that was used for cDNA synthesis was extracted from cells and whole embryos with the TRIzol reagent (Life Technologies), as well as RNA that was prepared for the RNA-seq experiment of MondoA morphants. TRIzol extraction was performed as described in the manufacturer’s protocol, with minor modifications. If not otherwise stated, 30 embryos were harvested in 1.5 ml TRIzol. Embryos were disrupted in the TRIzol solution with micropestles (Eppendorf), and 5 × 105 cells in 25 cm2 flasks were harvested with 1.5 ml TRIzol solution. Samples were then stored overnight at –80°C. After thawing the samples on ice, they were passed through a syringe four times (Braun Sterican, 0.45 × 25 mm). The integrity of the total RNA was checked on an agarose gel. A NanoDrop spectrometer was used to determine concentrations and A260/A280 ratios of the samples. Only RNA with an A260/A280 ratio above 1.6 was used. cDNA was synthesized using random primers and Superscript III reverse transcriptase according to the manufacturer’s protocol (Life Technologies). Transcript levels of genes of interest were determined with the StepOnePlus Real-Time PCR System (Applied Biosystems) using the QuantiTect SYBR Green PCR Kit (Qiagen). Copy numbers were normalized using β-actin transcript levels. The melting curves of the amplicons were determined after amplification. Ct values were determined using StepOne Software 2.1 (Applied Biosystems) and data were analyzed by the 2-ΔΔCt method as previously reported (Livak and Schmittgen, 2001). Primer sequences: β-actin: fw: 5'-gcctgacggacaggtcat-3', rv: 5'-accgcaagattccataccc-3'; eef1a: fw: 5'-ccttcgtccaatttcagg-3', rv: 5'-ccttgaaccaccccatgt-3'; mlxip: fw: 5'-ctgccagacgtatcacttcg-3', rv: 5'-ggcgaatcttgtctctccac-3'; mlxipl: fw: 5'-ccacagccagtgttgatgat-3', rv: 5'-tagctgtgcgatgaatgtcc-3'; mlx: fw: 5'-agacttctccatggccacac-3', rv: 5'-ctgatgggccttcacaatg-3'; pck1: fw: 5'-tgacgtcctggaagaacca-3', rv: 5'-gcgtacagaagcgggagtt-3'; hk2: fw: 5'-gaaattggcctcattgtcg-3', rv: 5'-catccaccagctccaagtg-3'; fasn: fw: 5'-aggccatcgttgatggag-3', rv: 5'-tgtagacgccagttttgctg-3'; txnipa: fw: 5'-ccaactagacgaacatccaaca-3', rv: 5'-agacaccagctgcccttg-3'; sox3: fw: 5'-cttagcgcacaactttgcag-3', rv: 5'-caccagtcccgtgtgtctc-3'; tmprss4b: fw: 5'-tcaaagtttcctctcagccagt-3', rv: 5'-ccaccacctcaccacagtc-3'; adh5: fw: 5'-cacgctcctctggataaagtgt-3', rv: 5'-gtagagcccgcttcaacct-3'; cfl1: fw: 5'-atgatctacgccagctccaa-3', rv: 5'-tcacttgccactcgtgctta-3'; vent: fw: 5'-aggagaaatgcagcacagc-3', rv: 5'-tcactctccacatcggtgtatc-3'; lmnl3: fw: 5'-cctggccaactacatcgag-3', rv: 5'-gatctgcatggaggatttgtc-3'; dynlrb1: fw: 5'-aacagcggtttgtcttcaca-3', rv: 5'-atggggattccttctgcatt-3'; gabarapl2: fw: 5'-cgacaaaactgtccctcagtc-3', rv: 5'-ccgtcctcgtcttgttcttt-3'; aurkaip1: fw: 5'-cgagagtggtttcacgactg-3', rv: 5'-ctggctctgtgtctgcaatg-3'; nsdhl: fw: 5'-ggggacctctgtgacaaaca-3', rv: 5'-aggtgaggcacaatgaaaca-3'.
Microinjection
Request a detailed protocolIf not otherwise stated, all injections were performed as follows: Injection mixes were prepared with the desired concentrations of the compound of interest (see below) and 0.1% (w/v) phenol red as an injection control. Each injection mix was injected with borosilicate glass capillary tubes (~0.6 mm diameter) into the yolk of 1 cell stage zebrafish eggs with a gas-driven microinjector (Eppendorf Femtojet express) as described (Müller et al., 1999). After the injection, unfertilized or damaged eggs were removed around the high to sphere stages and were not considered for the subsequent analysis.
Whole-mount in situ hybridization (WISH) and Whole-mount immunohistochemistry (WIHC)
Request a detailed protocolWISH was performed on embryos (n ≥ 20) as previously described (Oxtoby and Jowett, 1993), with minor modifications: The Proteinase K digest step was omitted for early embryonic stages and replaced with four PTw (1xPBS, 0.1% Tween-20) washes, while it was reduced to 6 min for 24 hr old embryos. Epon embedding of embryos was carried out as described (Westerfield, 2007). 5 µm thick sections were imaged with a compound microscope (Leica DM5000). Percentage of epiboly progression of the embryos was determined by performing a ntl in situ hybridization to label the blastoderm margin and using the following formula as previously described (Hsu et al., 2006): Epiboly (%) = (The length between animal pole and blastoderm margin)/(The length between animal pole and vegetal pole).
Microtubule staining of early embryos by WIHC was carried out as recommended by McMenamin et al., 2003. Uninjected embryos and mondoa-mo injected embryos were fixed at late blastula stage just before the onset of epiboly and stained for α-tubulin to visualize the microtubules. The primary antibody used was an anti-α-tubulin antibody (Sigma-Aldrich, 1:500), the secondary antibody was anti-mouse Alexa Fluor 488 (Life technologies, 1:1000). Phalloidin staining of actin filaments was performed according to Nakajima and Burke, 1996. Briefly, embryos were fixed in 4% PFA over night at 4°C. The next day, fixed embryos were washed three times in PTw and then incubated with Rhodamine-Phalloidin (1:500, R415, Molecular Probes (Life Technologies)) for 15 min. Embryos were mounted in 0.5% (w/v) low melting agarose in E3 and imaged with a confocal microscope (Leica TCS2 SP5). As observed for epiboly movement disruption, the severity of the YSL disorganization appears to be dose-dependent, with lower doses of MO leading to a less severe phenotype.
Cell culture
Cell lines and maintenance
Request a detailed protocolZebrafish PAC2 cells (Lin et al., 1994) were maintained essentially as described (Whitmore et al., 2000). Briefly, cells were kept in PAC2 maintenance medium (L-15 Medium, 17% (v/v) FBS, 100 U/ml penicillin/streptomycin, 50 mg/ml gentamycin) in 75 cm2 cell culture flasks at 28° C and passaged every week by trypsinization at a ratio of 1:7 to 1:10. HepG2 cells (ATCC, #HB-8065) were kept in HepG2 maintenance medium (DMEM, 10% (v/v) FBS, 100 U/ml penicillin/streptomycin) at 37°C and 5% CO2 in a humidified incubator and were passaged twice a week at the same ratio. Cells were not authenticated by STR profiling. Tests for mycoplasma contamination were performed every 8 weeks, showing negative results.
Cell transfection
Request a detailed protocolTransient transfections were carried out with FuGENE HD reagent (Promega) following the manufacturer’s protocol. Briefly, one day before transfection, cells were washed once with 1x PBS without Ca2+ or Mg2+, then trypsinized at room temperature to detach cells. To obtain 80% cell confluency on the day of transfection, cell number was determined with a Neubauer chamber using the Trypan blue exclusion method. Cells were diluted with PAC2 maintenance medium and transferred into 6-well plates (800,000 cells per six well), then incubated overnight at 28°C. The next day, maintenance medium was replaced with maintenance medium without antibiotics. The DNA complex was prepared using 100 µl L-15 medium, 1 µg total DNA of interest and 4 µl FuGENE HD (4:1). The mixture was incubated for 20 min at room temperature to ensure complex formation before the entire reagent mixture was added to the cells. Cells were incubated with the reagent overnight at 28°C. The following day, cells were transferred into opaque 96-well plates with PAC2 maintenance medium and were allowed to attach overnight before measurement of bioluminescence (see below). This protocol was modified for the HepG2 cells as follows: 200,000 cells per six wells, DNA complex consisting of 100 µl L-15 medium, 2 µg total DNA of interest and 8 µl FuGENE HD (8:2).
To measure glucose induction of transcription from a ChoRE reporter, cells were transfected either with the 2xChoRE reporter construct or with pGL3-Control (Promega) as a negative control. Treatment of cells with different concentrations of glucose was performed after starvation of HepG2 cells for 1 hr and of PAC2 cells overnight. The starvation medium (i.e. maintenance medium without serum and phenol red) was supplemented with 0.5 mM luciferin. Bioluminescence was determined after 24 hr of glucose incubation (see below).
For overexpression studies with Mondo pathway factors, PAC2 cells were transfected with the 2x ChoRE reporter and co-transfected with full-length mondoa, mlx, or both. The DNA complex was prepared using 600 ng 2xChoRE reporter DNA, 200 ng pRL-CMV vector (to control for transfection efficiency; Promega #E2261), and a total of 200 ng of transcription factor and carrier DNA (pBluescript SK). For the overexpression of one transcription factor, 100 ng transcription factor and 100 ng carrier DNA were used. For the overexpression of two transcription factors, 100 ng of each transcription factor were employed. The zebrafish nr3c1 (GRα) construct (Schaaf et al., 2008) was kindly provided by M. Schaaf (Institute of Biology Leiden, Leiden University, Netherlands). After transfection, cells were starved and then treated for 24 hr with either low (0.3 mM) or high (12 mM) glucose concentrations before bioluminescence was determined. Bioluminescence was measured in cell lysates and normalized to Renilla luciferase activity to control for transfection efficiency.
Morpholino oligonucleotide mediated knockdown in PAC2 cells
Request a detailed protocolTransient transfection of MOs was done by electroporation essentially as previously described (Weger et al., 2011). Briefly, cells (80% confluent 1/2 flask 175 cm2 per transfection) were electroporated with a MO mix containing 18 µg 2xChoRE reporter construct, 6 µg transfection control (pRL-CMV vector), 6 µg carrier DNA (pBluescript SK), 10 µM MO (mondoa-mis or mondoa-mo), in 50 µl dH2O at 0.29 kV, 0 W, 960 mF using a GenePulser II apparatus. Immediately after transfection cells were transferred into 4 ml PAC2 maintenance medium and resuspended. 150 µl of the cell suspension were transferred into 96 well plates containing 100 µl of maintenance medium. Cells were kept overnight at 28°C to attach. The next day cells were washed and then starved for 24 hr in starvation medium supplemented with 0.5 mM luciferin. The following day cells were treated for 24 hr with low (0.3 mM) or high glucose (12 mM) concentrations.
Bioluminescence measurements
Request a detailed protocolIn vivo luciferase measurements of transiently transfected cells was performed at 28°C with an EnVision multilabel plate reader (PerkinElmer) as previously described (Weger et al., 2011; Weger et al., 2012), using L-15 medium without phenol red and supplemented with 0.5 mM luciferin. For HepG2 cells, DMEM medium instead of L-15 medium was used to prepare an analogous luciferin-containing medium. For in vitro luciferase measurements (overexpression and MO knockdown experiments), we employed the Dual-Luciferase Reporter Assay System (Promega, #E1910) as suggested by the manufacturer.
2-Deoxy-D-glucose (2-DG) treatment of embryos
Request a detailed protocolEmbryos were treated with the glucose analog 2-deoxy-glucose (2-DG, Sigma) to activate the Mondo pathway. As Glucose-6-phosphate (G6P) is thought to be the signal activating the Mondo pathway, and because 2-DG can only be metabolized to 2-DG-6-phosphate, but not further, this treatment should avoid activation of other pathways relying on further metabolization of glucose (Li et al., 2010b; Stoltzman et al., 2008).
For recording of in vivo bioluminescence from a ChoRE reporter, a mixture composed of the 2xChoRE luciferase reporter construct (100 ng/µl; see above, ‘cell culture’), phenol red (0.1%) and either 150 mM 2-DG or water was injected into zygotes. The injected embryos were transferred into 96-well plates containing E3 medium with 0.5 mM luciferin (E3L) at high/oblong stage for bioluminescence measurement as described (Weger et al., 2013a; Weger et al., 2012; Weger et al., 2013b). Bioluminescence was measured at sphere stage and subtracted against background luminescence from wells without embryos.
To determine whether nsdhl is a downstream target of MondoA and is regulated by glucose, 150 ng/µl each of mondoa and mlx transcript were injected either with water as a control or with 150 mM 2-DG. RNA was extracted at sphere stage and processed for RT-qPCR against nsdhl. To determine if glucose induction of txnipa expression was affected by mondoa knockdown, embryos were injected with 150 mM 2-DG or with water along with either mondoa-mis or mondoa-mo.
CRISPR/Cas9 mediated zebrafish gene editing
Request a detailed protocolCRISPR/Cas9 mediated gene editing was performed essentially as described (Gagnon et al., 2014), with the following modifications. Guide RNA target regions in exon 3 of mlxip/mondoa (5’-gagtccatattgccatccaa-3’) were determined using the CHOPCHOP web tool (http://chopchop.cbu.uib.no/) (Labun et al., 2019). The guide RNA templates that include an SP6 promotor and Cas9 scaffold were prepared by GeneArt Strings DNA fragment gene synthesis (ThermoFisher) (5’-tataagcttccatggatttaggtgacactatagagtccatattgccatccaagttttagagctagaaatagcaagttaaaataaggctagtccgttatcaacttgaaaaagtggcaccgagtcggtgcttttctagacgatgcccttgagagccttcaacccagtcctatagtgagtcgtattaggatcctac-3’). The templates were amplified by PCR and the amplification product was transcribed using the MEGAscript SP6 Transcription Kit (ThermoFisher). 60 ng/µl guide RNA were injected into 1 cell stage wild-type zebrafish embryos together with 0.9 µg/µl Cas9 protein (GeneArt Platinum Cas9 Nuclease; ThermoFisher), 0.2 M KCl and 0.1% phenol red (SIGMA-Aldrich). Guide RNA efficiency was determined using High Resolution Melting (HRM) analysis (fw: 5’-tgcctttctttcctgaagatgt-3’; rv: 5’-gctttttccattcaaaaccagt-3’). Injected (F0) zebrafish were outcrossed with wild-type zebrafish. F1 mutant carriers were identified by HRM and PCR targeting the mutated region (PCR: fw: 5’-cgcatactcctttatcttgc-3’; rv: 5’-catgctttttccattcaaaacc-3’) on genomic DNA from fin biopsies followed by sequence analysis. Further matings of the F1 fish generated homozygous F2 fish, which were raised to adulthood. Maternal zygotic (MZ) mutants were derived from incrosses of homozygous F2 fish.
Morpholino oligonucleotide knockdown in zebrafish embryos
Request a detailed protocolMorpholino oligonucleotides (MOs) were obtained from GeneTools Inc (http://www.gene-tools.com/). MO sequences were: chrebp-mo (translation blocking MO): 5’-ttctgaatactctgcttttgccatc-3’; mondoa-mo (translation blocking MO): 5’-ggtattgtcgagtagccatgttaaa-3’; mondoa-mis (5 bp mismatch MO): 5’-ggtaatctccagtacccatcttaaa-3’; mlx-mo (translation blocking MO): 5’-cgcgctgttttccgtcattttggaa-3’; nsdhl-mo (splice blocking MO): 5’- ctgaaacaattcacacctgtttgtc-3’; p53-mo (Robu et al., 2007) (translation blocking MO): 5’-gcgccattgctttgcaagaattg-3’; tagged-mo (3’ lissamine tagged MO of unrelated sequence): 5’-gtagcttgtactcgcattcctatct-3’. To target the transcript of interest specifically in the YSL, embryos were injected with MO into the yolk at the 1 k cell stage, as previously described (Sakaguchi et al., 2001). In order to label the area where the MO is localized, either lis-mondoa-mo or lis-mondoa-mis was used, or 0.08 mM of tagged-mo was added to the injection mixture. Injection of tagged-mo alone did not cause an aberrant phenotype.
Knockdown efficiency of the mondoa-mo was tested with 30 ng/µl of the mondoa target sequence-GFP fusion transcript encoded by the pCS2:GFP:MondoA vector. For mRNA transcription, this construct was linearized with NotI and transcribed with the mMESSAGE mMACHINE kit as recommended by the manufacturer (Ambion). Rescue experiments with mondoa mRNA were performed using 250 ng/µl of transcript. mondoa mRNA was transcribed using the mMESSAGE mMACHINE Kit (Ambion) (T7) on templates linearized with DraIII. Pregnenolone (20 µM) rescue experiments were performed as previously reported (Hsu et al., 2006). Knockdown efficiency of the chrebp-mo was tested with 100 ng/µl pCS2:GFP:ChREBP. Resistance of MZmondoa mutant embryos to mondoa knockdown was tested with 0.5 mM lis-mondoa-mo and 0.5 mM lis-mondoa-mis injections into MZmondoa mutants and wild-type controls.
Confocal microscopy
Request a detailed protocolFor live imaging of embryos by confocal microscopy, embryos were dechorionated and transferred with a glass pipette into a chamber containing 0.1% (w/v) low melting agarose in E3. Images were taken on a Leica TCS SP5 upright microscope using either a HCX PL APO 20x/0.70 lambda blue IMM CORR or a HCX APO L 20x/0.50 W U-V-I objective at 26°C. Spatial and sensitivity resolution were strictly kept constant between samples of comparative analysis. Reflection images using 458/488/496 nm incident light were used to assess whether a nucleus is in the EVL or YSL. The reflected light images require no labeling and provide a confocal contrast at the boundary of different tissues (Jester et al., 1991). We used the GFP fluorescence channel to identify the nucleus and the label-free reflection channel to visualize the extracellular matrix. Whether a nucleus belongs to EVL or YSL was individually determined in the XZ or YZ plane using the reflection channel. Images were analyzed and processed in Fiji/ImageJ (Schindelin et al., 2012).
To visualize endocytotic vesicles, embryos were dechorionated at five hpf and placed in 1.5% (w/v) Lucifer Yellow CH (Sigma-Aldrich, L0259) dissolved in E3 medium for 15 min (Solnica-Krezel and Driever, 1994). After the incubation, embryos were washed briefly with excess E3 medium and mounted in 0.5% low-melting temperature agarose for confocal imaging.
Digital scanned laser light sheet microscopy (DSLM)
Request a detailed protocolTime-lapse DSLM images of developing zebrafish embryos were acquired on a home-built microscope described in Kobitski et al., 2015. In the DSLM (Keller et al., 2008; Schmid et al., 2013), the sample is irradiated with two opposing beams of 488 nm laser light (LuxX 488–60 laser, Omicron-Laserage Laserprodukte GmbH, Rodgau, Germany) that were scanned up and down so as to excite fluorescence in the samples in a slab (‘light sheet’) with an average thickness of 6.6 µm across the field of view (1038 × 875 µm2). Fluorescence emission from the illuminated region was collected perpendicular to the light sheet by using a water dipping objective (16x/0.8 w, Nikon) and projected on a 5.5 megapixel sCMOS camera (Neo, Andor, Belfast, UK), yielding 1.5 µm lateral image resolution. Image stacks of 500 frames at 2 µm displacement were collected at a rate of 20 frames per second by moving the sample through the light sheet. To achieve optimal image quality across the entire sample, images from two opposing views were acquired for each light sheet position by rotating the samples by 180°. Consequently, a complete 3D image stack was collected every 60 s.
For DSLM imaging, dechorionated embryos were embedded in 0.1% (w/v) low melting agarose in 1x E3 medium, mounted in fluorinated ethylene propylene tubes as previously described (Kaufmann et al., 2012), and placed inside a sample chamber. During image acquisition over 15 hr, the temperature in the sample chamber was kept constant at 26°C.
Automated image processing and analysis
Fluorescently labeled nuclei were automatically extracted from the 3D image stacks using the fast segmentation approach described previously (Stegmaier et al., 2014). In a post-processing step, the extracted features from complementary rotation images were registered using a bead-based registration technique (Preibisch et al., 2010). Moreover, redundant nuclei were fused on the basis of spatial overlap of segments; a feature based object rejection was performed to minimize the amount of false positive detections. Nucleus counts reflect the number of remaining objects after information fusion and object rejection.
To track the detected nuclei over multiple time points, we made use of a nearest neighbor tracking algorithm implemented in the open-source MATLAB toolbox SciXMiner (Mikut et al., 2017; Stegmaier et al., 2012). To attain more comprehensive tracks, we used an additional heuristic to reconstruct nuclei that were missing in a single frame and we fused extracted sub-tracks in the spatio-temporal domain by applying an additional nearest neighbor matching on the start and end points of sub-tracks. It should be noted, though, that these tracks did not necessarily reflect perfect tracks on the single cell level over the whole duration, but merely are used to identify and trace the qualitative movements of nuclei in certain regions.
To compute the cell density distribution, we counted the number of nuclei located within a bounding sphere with a radius of 40 µm centered on the nucleus of interest. For visualization, we aligned the virtual embryo by manually centering it on the origin and rotating the vegetal pole such that it was centered on the positive y-axis and finally converting the Cartesian coordinates into a spherical projection using azimuth, elevation and radius (see schematic in Figure 3B). This representation allows us to unwrap the embryo by plotting the azimuth (range: [-pi, pi]) and the elevation (range: [-pi/2, pi/2]) in a 2D graph.
For speed determination of ingressing deep cells, we manually selected a minimum number of ten cells for each condition and tracked their movement over a time window of 100 frames, that is, for 200 min. We included only internalizing deep cells which could be tracked over the entire time frame in both conditions. The time point of ingression onset was visually determined using slice-based maximum intensity projection videos. The estimated speed reflects the arithmetic mean of approximately 1000 displacement steps of the selected cells during the time interval.
Implementation details
Request a detailed protocolMaximum intensity projections, seed detection and segmentation were performed using a custom-built C++ application based on the Insight Toolkit (ITK). To speed up calculations, we processed multiple image stacks in parallel on a Hadoop cluster containing 25 nodes (each equipped with two quad-core Intel Xeon E5520 CPUs @ 2.27 GHz and 36 GB of memory) of the Large Scale Data Facility (LSDF) of the KIT (Garcia et al., 2011). Tracking, data analysis and visualizations were performed using MATLAB with the open-source toolbox SciXMiner (Mikut et al., 2017; Stegmaier et al., 2012) and Fiji (Schindelin et al., 2012).
Metabolite quantification
Request a detailed protocolATP, ADP and AMP levels in mondoa-mo and mondoa-mis injected embryos were determined at sphere stage, just before the morphological phenotype becomes apparent, using a mass spectrometry approach. Metabolites were extracted from three pools of 30 zebrafish embryos of each condition with 0.3 ml of 0.1 M HCl in an ultrasonic ice bath for 10 min. The resulting homogenates were centrifuged twice for 10 min at 4°C and 16.400 g to remove cell debris. AMP, ADP and ATP were derivatized with chloroacetaldehyde as described (Bürstenbinder et al., 2007) and separated by reversed phase chromatography on an Acquity BEH C18 column (150 mm x 2.1 mm, 1.7 µm, Waters) connected to an Acquity H-class UPLC system. Prior to separation, the column was heated to 42°C and equilibrated with five column volumes of buffer A (5.7 mM TBAS, 30.5 mM KH2PO4 pH 5.8) at a flow rate of 0.45 ml min−1. Separation of adenosine derivates was achieved by increasing the concentration of buffer B (2/3 acetonitrile in 1/3 buffer A, by volume) in buffer A as follows: 1 min 1% B, 1.6 min 2% B, 3 min 4.5% B, 3.7 min 11% B, 10 min 50% B, and return to 1% B in 2 min. The separated derivatives were detected by fluorescence (Acquity FLR detector, Waters, excitation: 280 nm, emission: 410 nm, gain: 100) and quantified using ultrapure standards (Sigma). Data were acquired and processed with the Empower3 software suite (Waters). Energy charge values were calculated for each batch of embryos separately according to ([ATP] + 0.5 [ADP])/([ATP] + [ADP] + [AMP]) and were analyzed for differences between conditions by One-Way ANOVA.
Next-generation sequencing
Generation and sequencing of mRNA and total RNA libraries
Request a detailed protocolRNA-seq analysis was performed on embryos injected with mondoa-mo or mondoa-mis as well as on MZmondoa mutants affected or unaffected in epiboly and wild-type embryos. MO injected embryos were sampled when control embryos were at sphere stage (age-matched), mutant and wild-type embryos when affected mutants could be clearly distinguished (90% epiboly of controls). Triplicate total RNA samples of each condition were used to generate cDNA libraries. The quality and integrity of the sampled total RNA were assessed with a Bioanalyzer 2100 (Agilent). Sample concentrations were measured using a NanoDrop spectrometer (ThermoScientific). Sequencing libraries were prepared from total RNA with a RNA integrity number (RIN) greater than 8.0. Indexed mRNAseq libraries were generated from morphant samples following the TruSeq mRNA kit protocol (Illumina). Libraries for mutant samples were prepared using the TruSeq Stranded Total RNA Sample Prep Kit with Ribo-Zero Gold depletion (Illumina). Libraries of the morphant experiments were sequenced on a HiSeq 1000 system (Illumina) in 3 lanes of a paired-end flow cell. For each end, 50 nucleotides were sequenced. Cluster and sequencing reactions were performed with the TruSeq SBS kit v3-HS and the paired-end cluster v3-cBot-HS kit. The libraries of the MZmondoa mutants were sequenced on Illumina Hiseq 4000 sequencer as Paired-End 50 base reads following Illumina’s instructions. Image analysis and base calling were performed using RTA 2.7.7 and bcl2fastq2.17.1.14.
Mapping and differential gene expression analysis
Request a detailed protocolPaired-end reads were mapped onto the Danio rerio genome assembly (GRCz11/danRer11) using STAR 2.3.8 (Dobin et al., 2013). Uniquely mapped reads were counted for each gene locus defined in the Ensembl 96 release. Briefly, we applied SAMtools view (Li et al., 2009) to assign uniquely mapped paired-end reads. For downstream applications such as data visualization, we normalized the count data by size factor and applied a variance stabilizing transformation (Huber et al., 2002). Moderated log fold changes were calculated as proposed by Love et al., 2014. DESEq2 (Love et al., 2014) was used to assess differential gene expression between two conditions. To compare the gene expression signatures of MZmondoa mutants with and without epiboly phenotype to those of mondoa morphants and their respective controls, we employed the rank-rank hypergeometric overlap (RRHO) algorithm (Plaisier et al., 2010). Competitive gene set enrichment analyses were performed using CAMERA (Wu and Smyth, 2012). GO term annotations were included in the analysis. In addition, a manually curated set of zebrafish genes linked with important metabolic pathways assembled from KEGG and Ensembl annotations (Weger et al., 2016) complemented the approach. We visualized gene set enrichment using the barcodeplot function from the limma package (Ritchie et al., 2015). Raw files and technical details about the RNA-Seq data have been deposited in NCBI's Gene Expression Omnibus (Edgar et al., 2002) and are accessible through GEO Series accession number GSE144350 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE144350).
Statistical analysis
Request a detailed protocolBioluminescence measurements in embryos were analyzed using Mann–Whitney U test, as the data were not normally distributed. Data sets with a normal distribution, including RT-qPCR data, velocity comparison of ingressing deep cells, epiboly progression, cell based bioluminescence measurements and energy charge of embryos, were analyzed with a Student’s t-test (for two groups) or ANOVA (for more than two groups). Multiple testing was corrected by the Benjamini-Hochberg procedure (Benjamini and Hochberg, 1995) or an adequate ANOVA posttest. Bonferroni’s posttest was chosen if all groups had to be compared, while Dunnett’s posttest was chosen if the groups were compared only to the control. Data were analyzed and visualized by GraphPad software (GraphPad Software, Inc), Microsoft Excel or R (https://www.r-project.org/). ‘n’ indicates number of biological repeats. Asterisks indicate p-values of *p≤0.05, **p≤0.01 and ***p≤0.001.
Data availability
Raw files and technical details about the RNA-Seq data have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE144350 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE144350).
-
NCBI Gene Expression OmnibusID GSE144350. Regulation of cholesterol biogenesis by the glucose-sensing transcription factor MondoA is required for zebrafish epiboly.
-
ENAID accession-number-ERP014517. A high-resolution mRNA expression time course of embryonic development in zebrafish.
References
-
Sweet sixteenth for ChREBP: established roles and future goalsCell Metabolism 26:324–341.https://doi.org/10.1016/j.cmet.2017.07.004
-
Synthesis of putative microtubule-associated proteins by mouse blastocysts during early outgrowth in vitroCanadian Journal of Biochemistry and Cell Biology 62:885–893.https://doi.org/10.1139/o84-113
-
Controlling the false discovery rate: a practical and powerful approach to multiple testingJournal of the Royal Statistical Society: Series B 57:289–300.https://doi.org/10.1111/j.2517-6161.1995.tb02031.x
-
Mutational spectrum of NSDHL in CHILD syndromeJournal of Medical Genetics 42:e17.https://doi.org/10.1136/jmg.2004.024448
-
Zebrafish epiboly: spreading thin over the yolkDevelopmental Dynamics 245:244–258.https://doi.org/10.1002/dvdy.24353
-
The yolk syncytial layer in early zebrafish developmentTrends in Cell Biology 20:586–592.https://doi.org/10.1016/j.tcb.2010.06.009
-
Enzymatic assays for 2-deoxyglucose and 2-deoxyglucose 6-phosphateAnalytical Biochemistry 161:508–513.https://doi.org/10.1016/0003-2697(87)90481-7
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Gene expression omnibus: ncbi gene expression and hybridization array data repositoryNucleic Acids Research 30:207–210.https://doi.org/10.1093/nar/30.1.207
-
Cooperation of sonic hedgehog enhancers in midline expressionDevelopmental Biology 301:578–589.https://doi.org/10.1016/j.ydbio.2006.11.004
-
ConferenceData management and analysis at the large scale data facilityPaper Presented at the 6th International Conference on Digital Information Manag1—Source Data 1ement, Melbourne.
-
Little fish, big data: zebrafish as a model for cardiovascular and metabolic diseasePhysiological Reviews 97:889–938.https://doi.org/10.1152/physrev.00038.2016
-
X-Linked dominant disorders of cholesterol biosynthesis in man and mouseBiochimica Et Biophysica Acta (BBA) - Molecular and Cell Biology of Lipids 1529:357–373.https://doi.org/10.1016/S1388-1981(00)00160-8
-
Zebrafish cyp11a1 and hsd3b genes: structure, expression and steroidogenic development during embryogenesisMolecular and Cellular Endocrinology 312:31–34.https://doi.org/10.1016/j.mce.2009.07.030
-
The glucose-sensing transcription factor MLX promotes myogenesis via myokine signalingGenes & Development 29:2475–2489.https://doi.org/10.1101/gad.267419.115
-
The homeobox genes vox and vent are redundant repressors of dorsal fates in zebrafishDevelopment 128:2407–2420.
-
MondoA deficiency enhances sprint performance in miceBiochemical Journal 464:35–48.https://doi.org/10.1042/BJ20140530
-
In vivo, real-time confocal imagingJournal of Electron Microscopy Technique 18:50–60.https://doi.org/10.1002/jemt.1060180108
-
Life at the edge: epiboly and involution in the zebrafishResults and Problems in Cell Differentiation 40:117–135.https://doi.org/10.1007/978-3-540-46041-1_7
-
Regulation of ATP-citrate lyase gene transcriptionYonsei Medical Journal 37:214–224.https://doi.org/10.3349/ymj.1996.37.3.214
-
Stages of embryonic development of the zebrafishDevelopmental Dynamics 203:253–310.https://doi.org/10.1002/aja.1002030302
-
Glucose regulation of mouse S(14) gene expression in hepatocytes involvement of a novel transcription factor complexThe Journal of Biological Chemistry 275:5200–5207.https://doi.org/10.1074/jbc.275.7.5200
-
CHOPCHOP v3: expanding the CRISPR web toolbox beyond genome editingNucleic Acids Research 47:W171–W174.https://doi.org/10.1093/nar/gkz365
-
Zebrafish epiboly: mechanics and mechanismsThe International Journal of Developmental Biology 54:1213–1228.https://doi.org/10.1387/ijdb.093028sl
-
The sequence alignment/Map format and SAMtoolsBioinformatics 25:2078–2079.https://doi.org/10.1093/bioinformatics/btp352
-
Maternal control of early mouse developmentDevelopment 137:859–870.https://doi.org/10.1242/dev.039487
-
Glucose-6-phosphate mediates activation of the carbohydrate responsive binding protein (ChREBP)Biochemical and Biophysical Research Communications 395:395–400.https://doi.org/10.1016/j.bbrc.2010.04.028
-
ChREBP*mlx is the principal mediator of glucose-induced gene expression in the liverJournal of Biological Chemistry 281:28721–28730.https://doi.org/10.1074/jbc.M601576200
-
Patterning the early zebrafish by the opposing actions of bozozok and vox/ventDevelopmental Biology 224:275–285.https://doi.org/10.1006/dbio.2000.9780
-
Revisiting the role of metabolism during developmentDevelopment 145:dev131110.https://doi.org/10.1242/dev.131110
-
Intronic enhancers control expression of zebrafish sonic hedgehog in floor plate and notochordDevelopment 126:2103–2116.
-
Cloning of the zebrafish krox-20 gene (krx-20) and its expression during hindbrain developmentNucleic Acids Research 21:1087–1095.https://doi.org/10.1093/nar/21.5.1087
-
A zebrafish histone variant H2A.F/Z and a transgenic H2A.F/Z:gfp fusion protein for in vivo studies of embryonic developmentDevelopment Genes and Evolution 211:603–610.https://doi.org/10.1007/s00427-001-0196-x
-
MondoA-Mlx heterodimers are candidate sensors of cellular energy status: mitochondrial localization and direct regulation of glycolysisMolecular and Cellular Biology 26:4863–4871.https://doi.org/10.1128/MCB.00657-05
-
Molecular genetics of Axis formation in zebrafishAnnual Review of Genetics 39:561–613.https://doi.org/10.1146/annurev.genet.37.110801.143752
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
No tail (ntl) is the zebrafish homologue of the mouse T (Brachyury) geneDevelopment 120:1009–1015.
-
Requirement of Npc1 and availability of cholesterol for early embryonic cell movements in zebrafishJournal of Lipid Research 52:1328–1344.https://doi.org/10.1194/jlr.M012377
-
Microtubule arrays of the zebrafish yolk cell: organization and function during epibolyDevelopment 120:2443–2455.
-
Gastrulation: making and shaping germ layersAnnual Review of Cell and Developmental Biology 28:687–717.https://doi.org/10.1146/annurev-cellbio-092910-154043
-
Glucose-Sensing transcription factor MondoA/ChREBP as targets for type 2 diabetes: opportunities and challengesInternational Journal of Molecular Sciences 20:5132.https://doi.org/10.3390/ijms20205132
-
Lithium perturbation and goosecoid expression identify a dorsal specification pathway in the pregastrula zebrafishDevelopment 117:1261–1274.
-
Information fusion of image analysis, video object tracking, and data mining of biological images using the open source MATLAB toolbox Gait-CADBiomedical Engineering / Biomedizinische Technik 57:4073.https://doi.org/10.1515/bmt-2012-4073
-
Ultraviolet irradiation impairs epiboly in zebrafish embryos: evidence for a microtubule-dependent mechanism of epibolyDevelopment 119:909–919.
-
A chemical screening system for glucocorticoid stress hormone signaling in an intact vertebrateACS Chemical Biology 7:1178–1183.https://doi.org/10.1021/cb3000474
-
A chemical screening procedure for glucocorticoid signaling with a zebrafish larva luciferase reporter systemJournal of Visualized Experiments 10:e50439.https://doi.org/10.3791/50439
-
Pregnenolone activates CLIP-170 to promote microtubule growth and cell migrationNature Chemical Biology 9:636–642.https://doi.org/10.1038/nchembio.1321
-
BookThe Zebrafish Book: A Guide for the Laboratory Use of Zebrafish Danio RerioEugene: University of Oregon Press.
-
Camera: a competitive gene set test accounting for inter-gene correlationNucleic Acids Research 40:e133.https://doi.org/10.1093/nar/gks461
-
Lipogenic transcription factor ChREBP mediates fructose-induced metabolic adaptations to prevent hepatotoxicityJournal of Clinical Investigation 127:2855–2867.https://doi.org/10.1172/JCI89934
Article and author information
Author details
Funding
Studienstiftung des Deutschen Volkes
- Meltem Weger
Marie Curie Intra-European Fellowship (PIEF-GA-2013-625827)
- Meltem Weger
Excellence cluster "CellNetworks"
- Thomas Dickmeis
Effect-Net research network in water research, Ministerium für Wissenschaft, Forschung und Kunst Baden-Württemberg
- Thomas Dickmeis
Association Française contre les Myopathies
- Thomas Dickmeis
EU FP7-ZFHEALTH-2007-B2
- Uwe Strähle
Interreg Network for Synthetic Biology in the Upper Rhine Valley
- Uwe Strähle
Helmholtz-Gemeinschaft (BioInterfaces in Technology and Medicine (BIFTM))
- Uwe Strähle
- Gerd Ulrich Nienhaus
- Ralf Mikut
- Thomas Dickmeis
Helmholtz-Gemeinschaft (Science and Technology of Nanostructures (STN))
- Gerd Ulrich Nienhaus
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Ferenc Müller for critically reading the manuscript and Sepand Rastegar, Nicholas S Foulkes, Marco Ferg and Gerald Brenner-Weiß for helpful discussions. Sepand Rastegar and Marcel Schaaf kindly provided plasmids. We acknowledge Christin Lederer and Giuseppina Pace for excellent technical assistance. Olivier Armant, Victor Gourain and colleagues of the NGS core facility of the BIFTM program and Bernard Jost and colleagues of the GenomEast platform at the IGBMC performed library preparations and sequencing of the in the mRNA-seq study and total RNA-seq study, respectively. Gernot Poschet of the Metabolomics Core Technology Platform of the Excellence cluster ‘CellNetworks’ (University of Heidelberg) ran the HPLC-based metabolite quantification, supported by the Deutsche Forschungsgemeinschaft (grant ZUK 40/2010–3009262). Ariel Garcia and Serguei Bourov (Steinbuch Centre for Computing at KIT, Helmholtz Program Supercomputing) as well as Volker Hartmann (Institute for Data Processing and Electronics, KIT) provided support for data handling and Hadoop processing. Funding was provided by the Studienstiftung des deutschen Volkes (to MW), by a Marie Curie Intra-European fellowship (grant PIEF-GA-2013–625827, to MW), by the Excellence Cluster ‘CellNetworks’ (University of Heidelberg, to TD), by the Effect-Net research network in water research of the Ministerium für Wissenschaft, Forschung und Kunst Baden-Württemberg (to TD), by the Deutsche Forschungsgemeinschaft (DFG grant DI913/4-1, to TD), by the Association française contre les myopathies (AFM, to TD), by EU FP7-ZFHEALTH-2007-B2 (to US), by the Interreg Network for Synthetic Biology in the Upper Rhine valley (NSB-Upper Rhine), and by the Helmholtz Programs BioInterfaces in Technology and Medicine (BIFTM) and Science and Technology of Nanostructures (STN).
Ethics
Animal experimentation: All zebrafish husbandry was performed in accordance with the German animal protection standards and approved by the Government of Baden-Württemberg, Regierungspräsidium Karlsruhe, Germany (Aktenzeichen35-9185.64/BH KIT), by the Home Office, United Kingdom (Scientific Procedures, Act 1986), and by the Service de la consommation et des affaires vétérinaires du canton de Vaud, Switzerland.
Copyright
© 2020, Weger et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,874
- views
-
- 131
- downloads
-
- 8
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Neuronal stem cells generate a limited and consistent number of neuronal progenies, each possessing distinct morphologies and functions, which are crucial for optimal brain function. Our study focused on a neuroblast (NB) lineage in Drosophila known as Lin A/15, which generates motoneurons (MNs) and glia. Intriguingly, Lin A/15 NB dedicates 40% of its time to producing immature MNs (iMNs) that are subsequently eliminated through apoptosis. Two RNA-binding proteins, Imp and Syp, play crucial roles in this process. Imp+ MNs survive, while Imp−, Syp+ MNs undergo apoptosis. Genetic experiments show that Imp promotes survival, whereas Syp promotes cell death in iMNs. Late-born MNs, which fail to express a functional code of transcription factors (mTFs) that control their morphological fate, are subject to elimination. Manipulating the expression of Imp and Syp in Lin A/15 NB and progeny leads to a shift of TF code in late-born MNs toward that of early-born MNs, and their survival. Additionally, introducing the TF code of early-born MNs into late-born MNs also promoted their survival. These findings demonstrate that the differential expression of Imp and Syp in iMNs links precise neuronal generation and distinct identities through the regulation of mTFs. Both Imp and Syp are conserved in vertebrates, suggesting that they play a fundamental role in precise neurogenesis across species.
-
- Computational and Systems Biology
- Developmental Biology
The initially homogeneous epithelium of the early Drosophila embryo differentiates into regional subpopulations with different behaviours and physical properties that are needed for morphogenesis. The factors at top of the genetic hierarchy that control these behaviours are known, but many of their targets are not. To understand how proteins work together to mediate differential cellular activities, we studied in an unbiased manner the proteomes and phosphoproteomes of the three main cell populations along the dorso-ventral axis during gastrulation using mutant embryos that represent the different populations. We detected 6111 protein groups and 6259 phosphosites of which 3398 and 3433 were differentially regulated, respectively. The changes in phosphosite abundance did not correlate with changes in host protein abundance, showing phosphorylation to be a regulatory step during gastrulation. Hierarchical clustering of protein groups and phosphosites identified clusters that contain known fate determinants such as Doc1, Sog, Snail, and Twist. The recovery of the appropriate known marker proteins in each of the different mutants we used validated the approach, but also revealed that two mutations that both interfere with the dorsal fate pathway, Toll10B and serpin27aex do this in very different manners. Diffused network analyses within each cluster point to microtubule components as one of the main groups of regulated proteins. Functional studies on the role of microtubules provide the proof of principle that microtubules have different functions in different domains along the DV axis of the embryo.