Expression and knockdown of zebrafish folliculin suggests requirement for embryonic brain morphogenesis
- Emma J. Kenyon†1, 4Email author,
- Monique N. H. Luijten†2,
- Harmeet Gill†1,
- Nan Li1,
- Matthew Rawlings1,
- James C. Bull3,
- Yavor Hadzhiev1,
- Maurice A. M. van Steensel5, 6,
- Eamonn Maher1 and
- Ferenc Mueller1
© Kenyon et al. 2016
Received: 4 August 2015
Accepted: 15 May 2016
Published: 8 July 2016
Birt-Hogg-Dubé syndrome (BHD) is a dominantly inherited familial cancer syndrome characterised by the development of benign skin fibrofolliculomas, multiple lung and kidney cysts, spontaneous pneumothorax and susceptibility to renal cell carcinoma. BHD is caused by mutations in the gene encoding Folliculin (FLCN). Little is known about what FLCN does in a healthy individual and how best to treat those with BHD. As a first approach to developing a vertebrate model for BHD we aimed to identify the temporal and spatial expression of flcn transcripts in the developing zebrafish embryo. To gain insights into the function of flcn in a whole organism system we generated a loss of function model of flcn by the use of morpholino knockdown in zebrafish.
flcn is expressed broadly and upregulated in the fin bud, somites, eye and proliferative regions of the brain of the Long-pec stage zebrafish embryos. Together with knockdown phenotypes, expression analysis suggest involvement of flcn in zebrafish embryonic brain development. We have utilised the zFucci system, an in vivo, whole organism cell cycle assay to study the potential role of flcn in brain development. We found that at the 18 somite stage there was a significant drop in cells in the S-M phase of the cell cycle in flcn morpholino injected embryos with a corresponding increase of cells in the G1 phase. This was particularly evident in the brain, retina and somites of the embryo. Timelapse analysis of the head region of flcn morpholino injected and mismatch control embryos shows the temporal dynamics of cell cycle misregulation during development.
In conclusion we show that zebrafish flcn is expressed in a non-uniform manner and is likely required for the maintenance of correct cell cycle regulation during embryonic development. We demonstrate the utilisation of the zFucci system in testing the role of flcn in cell proliferation and suggest a function for flcn in regulating cell proliferation in vertebrate embryonic brain development.
KeywordsBirt-Hogg-Dubé syndrome (BHD) Zebrafish flcn Folliculin zFucci Cell cycle Brain
Birt-Hogg-Dubé syndrome (BHD) is a dominantly inherited familial cancer syndrome characterised by the development of benign skin fibrofolliculomas, multiple lung and kidney cysts, spontaneous pneumothorax and susceptibility to renal cell carcinoma [1, 2]. The syndrome results from inactivating mutations in the gene encoding folliculin (FLCN) [3, 4] which has been suggested to be a rab35 specific Guanine Nucleotide Exchange Factor (GEF) . BHD patients only show a mutation in one copy of the FLCN gene , which suggests that several BHD symptoms may be due to abnormal levels of FLCN rather than its complete loss. This could explain why expression of mutant FLCN has been seen in a BHD-associated renal carcinoma . FLCN has been implicated in the regulation of various signalling pathways and cellular processes including cellular metabolism through mTOR, AMPK and HIF1α, transcriptional regulation, JAK-STAT signalling, cell adhesion, ciliogenesis, lysosomal biogenesis and autophagy [7–17]. However, several groups have generated conflicting data on the consequences of FLCN deficiency and how these lead to the clinical manifestations associated with BHD is not clear. This leads to the question of what FLCN does in a healthy individual and how best to treat those with BHD.
Towards answering this question investigations in a whole organism model system provide important insights into the nature and evolutionary conservation of BHD-related effects on target signalling pathways. As homozygous mutations of Flcn is early embryonic lethal in mice, rats and dogs  , the utilisation of alternative animal models may be more informative about the developmental functions of Flcn. The external development of zebrafish make it a valuable vertebrate system, in which to elucidate both the function of flcn and the molecular pathways of oncogenesis .
Here we study the role of flcn in zebrafish development using morpholino oligonucleotides to generate a zebrafish loss of function (LOF) BHD model and we hypothesise that this could provide insights into the biological functions of Folliculin. We aim to identify the temporal and spatial expression of flcn transcripts in the developing embryo and reconcile this with the phenotype associated with the morpholino knock-down of zebrafish flcn to gain an insight into what flcn may do in the developing embryo.
Isolation and expression of the zebrafish flcn gene
Morpholino knockdown indicates embryonic functions of folliculin
Knockdown of the zebrafish flcn gene by microinjection of either 100 μM of splice 2 (Fig. 2f, g and i), 100 μM of splice 1 (Additional file 1: Figure S1C) or 600 μM of ATG- targeting (Additional file 1: Figure S1E) flcn morpholino into wild-type zebrafish embryos resulted in a clear and reproducible phenotype (Additional file 2: Table S1) while control injected embryo appeared unaffected. At 18 somite stage embryos injected with flcn morpholino appeared to show developmental arrest when compared with mismatch controls. At Prim-5 stage experimental embryos showed a 2 h delay in development, with U-shaped somites and non-linear notochord in the tail (Fig. 2f) and cell death in the brain and trunk as implied by an increase in Acridine Orange staining (Fig. 2g), when compared to mismatch control morpholino injected embryos (Fig. 2d and e). At Long-pec stage the phenotype was characterised by brain oedema and problems with tail circulation as well as larger yolk and thinner yolk extension (Fig. 2i, Additional file 1: Figure S1D and S1F) when compared to mismatch controls (Fig. 2h and Additional file 1: Figure S1C and S1E) which were indistinguishable from the uninjected embryos.
To confirm that the morphant phenotype is caused by reduced gene expression of flcn, we aimed to perform rescue experiments by injection of the flcn splice 2 morpholino and full-length flcn RNA. Therefore, the full open reading frame of zebrafish flcn was cloned into the pCS2+ vector and RNA was produced in vitro. As injection of flcn RNA into mismatch morpholino injected embryos resulted in an increased number of embryos showing cyclopia, heart oedema, circulatory defects and reduced axis, rescue of MO knockdown was assessed by abolishment of the hydrocephalus and yolk phenotype. Increasing the amount of flcn RNA in embryos injected with flcn splice 2 morpholino resulted in an increase in the number of embryos that were phenotypically normal and a decrease in the number of embryos showing a severe hydrocephalus and yolk phenotype (Fig. 2j). These results suggest that the flcn morpholino is specifically targeting the flcn transcripts and that flcn may be required for the development of the zebrafish embryo particularly in the growing brain.
flcn morphants do not show kidney defects or defects in motile cilia in early development
As BHD patients have an increased risk of developing renal cysts and tumours  we decided to study the effect of flcn knockdown on kidney development in zebrafish embryos. The Wilms Tumor-1b GFP transgenic line (Tg(wt1b:EGFP)) expresses GFP exclusively in the pronephros and exocrine pancreas . Injection of flcn morpholino into Tg(wt1b:EGFP) did not reveal any obvious defects in pronephric development at Long-pec stage (Additional file 1: Figure S1H) when compared with mismatch control injected (Additional file 1: Figure S1G)
Morpholino knockdown of flcn in zebrafish resulted in a phenotype reminiscent of ciliary gene mutants  and morphants . In addition, BHD in humans was recently described as a novel ciliopathy . However, no clear difference in acetylated alpha tubulin staining (a known ciliary marker) in the pronephros and central canal cilia was observed at this magnification in Prim-5 stage embryos injected with flcn morpholino (Additional file 1: Figure S1J) when compared with mismatch control injected embryos (Additional file 1: Figure S1I).
Cell cycle defects in brain development in flcn morphant embryos
The lack of obvious defects in motile cilia morphology and a lack of defects in the developing kidney prompted us to examine other cellular processes that might be affected by Flcn deficiency, thereby causing the morphant phenotype. In vitro, FLCN has been implicated in cell cycle regulation via late S and G2/M phase and cyclin D1 [12, 24, 25]. In addition, the flcn morphants show developmental arrest and flcn is expressed in the proliferative zone of the brain. Therefore, we decided to examine cell cycle regulation in more detail in these zebrafish embryos.
In this study we show the result of gene expression analysis and morpholino knockdown of flcn during zebrafish early development. Our results suggest that flcn may be required for zebrafish embryonic brain development. Interestingly Birt-Hogg-Dubé (BHD) patients shown no defects in the brain however all known BHD patients only show a mutation in one copy of the FLCN gene . Canine and murine animal models homozygous for FLCN knockout show early embryonic lethality [27, 28]. Chen et al. and Hasumi et al. showed that mice homozygous for flcn mutation die between 3.5 dpc and 8.5 dpc  (E5.5–E6.5 ). Before death mutant mice lacked an organised epiblast and surrounding visceral endoderm. In embryos that did develop visceral endoderm and ectoderm, these layers had a disorganised structure containing cells with an enlarged cytoplasm and disorientated nuclei. In contrast flcn morpholino injected zebrafish, appear to develop normally past gastrulation with the first signs of disruption occurring when the embryos arrest at 18hpf. By Prim-5 stage the embryos exhibit cell death in the brain and trunk and by Long-pec stage the embryos show brain oedema and problems with tail circulation as well as a larger, more fragile yolk and thinner yolk extension.
The morpholino phenotype is supported by whole mount in situ hybridization of flcn showing expression in the tectum and telencephalon in zebrafish. In mice, flcn expression can be seen in the neural tube . At E7.5 Flcn showed a strong signal in the neural ectoderm of the primitive streak region and in the headfold of the ectoderm in mice. At E9.5, Flcn was expressed almost ubiquitously with a strong signal in the neural tube and optic pit. However, at E10.5 Flcn was highly expressed in the branchial arch, forelimb, and hind limb, as well as in the somites. Adult mice showed expression in the heart, pancreas and prostate with moderate expression in adult brain, kidney, liver and lung . These results appear consistent with what we see in zebrafish development where expression of flcn is also seen in the fin bud, somites, eye and regions of the brain. This suggests a critical role for flcn in vertebrate development.
In humans BHD has recently been described as a novel ciliopathy with FLCN being shown to localise to motile and non-motile cilia, centrosomes and the mitotic spindle. Alteration of FLCN levels caused changes to the onset of ciliogenesis . In zebrafish embryos lacking genes important for ciliogenesis one of the reported phenotypes is hydrocephalus , also seen in flcn morphants . It therefore seemed probable that knockdown of flcn in zebrafish embryos would result in defects in the pronephros and/or other ciliated structures. However we have found no evidence for this when we compare the pronephros of Tg(wt1b:EGFP) zebrafish embryos with and without flcn knockdown or when using immunofluorescent staining for acetylated alpha tubulin (a ciliary marker). It is possible that although cilia seem to be present in the flcn knockdown embryos, they may not be functional. Alternatively, mammalian FLCN is present in both motile and non-motile cilia . In this study we primarily focused on the effect of flcn morpholino knockdown on the appearance of motile cilia, as it has been shown that defects in these cilia can cause hydrocephalus and kidney cysts [23, 31, 32]. However aspects of the flcn morpholino phenotype such as hydrocephalus and cell cycle defects could also be explained by defects in the non-motile primary cilia, [33, 34] providing interesting avenues for further investigations. Alternatively, FLCN’s function in ciliogenesis may be specific to mammals.
Given the lack of obvious ciliary defects we sought to understand the developmental arrest phenotype. In particular whether this phenotype was caused by defects in the cell cycle as a recent study by Laviolette et al. (2013) suggested that flcn delays cell cycle progression through late S and G2/M phase . To do this we knocked down flcn in zFucci transgenic embryos where the G1 phase of the cell cycle is labelled red and S, G2 and M are labelled green. We found that at the 18 somite stage there was a significant drop in S-M states (Fig 3j) and a corresponding increase in G1 cells in flcn morpholino injected embryos (Fig 3e, f, h). This was particularly evident in the brain, retina and somites of the embryo, areas that have been shown in both mice  and zebrafish (Fig 1) to express flcn. To investigate at which developmental stage this divergence of zFucci transgene activities arises we undertook timelapse analysis of the head region of flcn morpholino injected and mismatch control embryos which show a gradual and increasing change over time. From this we can conclude that Flcn may play a role in the cell cycle. However, the exact nature of Flcn's control on the cell cycle remains to be determined, as recent in vitro studies have suggested that FLCN delays cell cycle progression through late S and G2/M phase  and that FLCN regulates cyclin D1 expression  which promotes G1 to S phase progression in the cell cycle . Kawai et al. found that RNAi knockdown of FLCN in HeLa cells led to increased cyclin D1 expression and, conversely, that reintroduction of the FLCN gene into FLCN-null Nihon rat kidney tumour cells reduced cyclin D1 expression. They surmised that these results suggested that FLCN inhibits cyclin D1 expression and, by extension, may also regulate G1 to S phase transition of the cell cycle. Conversely in zebrafish we found that knockdown of flcn resulted in no significant change in cyclin D1 (data not shown) suggesting that cell cycle defects in the developing embryo as a whole are likely to be more complex.
Alternatively, the cell cycle defects might be secondary to other consequences of Flcn knockdown. Cells in G1 can only progress if sufficient growth factors such as nutrient supply are present (reviewed ). In vitro, FLCN has been shown to regulate signalling of AMPK, a major sensor of energy status. FLCN directly interacts with AMPK  and loss of FLCN was shown to constitutively activate AMPK leading to mitochondrial biogenesis and HIF-driven aerobic glycolysis [17, 37]. Lui et al. (2013) have shown that Drosophila melanogaster with a genomic deletion of the Flcn gene never survive to adulthood and show characteristics of malnutrition . Zebrafish flcn morpholino injected embryos show large, fragile yolks and thin yolk extensions indicative of defects in nutrient supply and digestive problems [39–41]. It is therefore possible that the phenotype we see is a result of the zebrafish being unable to properly utilise its nutritional reserves.
Recent publications have shown that morpholino analysis may often not be corroborated by genetic mutants generated by reverse genetic approaches [42, 43]. The reason for this is unclear and has been suggested to be due to off target effects of the morpholino or hypomorphic alleles in the mutants. However new data by Rossi et al. may shed some light on this . Rossi et al. generated mutations in the zebrafish egfl7 gene, which showed no obvious phenotype. This was in contradiction to the severe vascular defects seen in wild type embryos injected with an egfl7 morpholino. After comparing the proteome and transcriptome of these mutant and morphant embryos they found a number of proteins and genes which were upregulated in the mutants but not the morphants. Amongst these were genes that could rescue the morphant phenotype. This could argue for a compensatory network that buffers against deleterious mutations .
Morpholino studies – including ours – are supported by controls including mismatch MO injection and mRNA rescue experiments . Since the specific phenotypes we describe are observed by the targeting MO only and the hydrocephalus and yolk phenotypes are rescuable by MO resistant mRNA, it is reasonable to speculate that the observed defects are a combination of a morpholino induced state (such as stress), which generates a yet unexplained embryo environment, in which the targeted flcn gene becomes limiting for normal development. In this morpholino injected environment the defects observed are due to the loss of the targeted flcn gene. The controls indicate that the defects cannot be attributed to flcn-independent off target effects and, therefore, this loss of function approach remains informative about linking flcn function to cell cycle regulation during development. Evidence from cell culture experiments suggests hypomorphic effects of FLCN , while mouse knockouts die in early development prohibiting meaningful phenotype analyses. Thus our morpholino knockdown study provides an informative starting point to unravel the function of flcn in a whole organism vertebrate system. How exactly flcn affects development and the cell cycle will be answered by the future generation of flcn mutant embryos. As flcn is known to have multiple interacting proteins [8, 46–51], we can not rule out that flcn mutants may employ a compensatory network to buffer against loss of flcn, a result that would also provide exciting new avenues of investigation.
In conclusion we show that zebrafish flcn is expressed during embryonic development with elevated expression levels in proliferating tissues of the zebrafish embryo. Using flcn knockdown coupled with a whole organism based in vivo cell proliferation assay we show that flcn is likely required for embryonic brain development and that it affects the normal progression of the cell cycle particularly in the growing brain and suggest a previously undescribed role for flcn in vertebrates.
Zebrafish husbandry and embryo generation
Zebrafish embryos were obtained from sibling crosses from adult AB fish housed at the fish facility at Birmingham University. Zebrafish were raised and bred and embryos staged following standard protocols [52, 53]. Embryos are raised at 28.5 °C. When needed to prevent pigment formation, embryos were raised in 0.003 % phenylthiocarbamide in E3 medium (1 mM NaCl, 0.17 mM KCl, 0.33 mM MgSO4.7H2O, 0.33 mM CaCl2.2H2O) from tailbud stage.
flcn cDNA cloning and mRNA production
The full open reading frame for zebrafish flcn was cloned from total RNA extracted from pooled protruding-mouth stage AB* wild type embryos using TRIZOL reagent according to manufacturers’ instructions. Next, cDNA was produced using M-MLV Reverse Transcriptase, RNase H Minus, Point Mutant (Promega M3681, M3682, M3683) and oligodT primers (Invitrogen). Zebrafish flcn cDNA was amplified using the following forward primer with BamHI restriction site: AATA GGATCC ATGAACGCTTTAGTTGCCCTG and reverse primer with Xba1 restriction site: AATA TCTAGA CCCGCTTTCAGTCTCTCTCAC and cloned into pCS2+ using BamHI (New England Biolabs) and XbaI (New England Biolabs). Plasmid was verified by sequence analysis.
Capped RNA was synthesised using 5 μg (5 μl) of NotI linearised flcn DNA using SP6 mMESSAGE mMACHINE kit (Ambion). The RNA was cleaned using GenElute™ Mammalian Total RNA Miniprep Kit (Sigma RTN70) as per manufacturer’s instructions Appendix 2.
Microinjection of embryos with RNA and morpholino
FLCN ATG MO AGGGCAACTAAAGCGTTCATCTGTG
FLCN splice1 MO ATGACACTCCCCTCTCGCTCACCTC
FLCN splice1 MM ATcACAgTCCCCTgTCGCTgACgTC
FLCN splice2 MO CGTTCATCTGGAGGAAACAAACATA
FLCN splice2 MM CGTTCtTgTGcAGGAtACAAAgATA
Morpholinos were diluted in MO buffer (5 mg/ml phenol red (Sigma), 4 mM HEPES pH 7.2 (Sigma), 160 mM KCl (Sigma)) and 1.4 nl of MO solution was injected into the yolk of the 1 cell stage embryo. RNA was diluted to 100-300 pg/μl in nuclease free water and 1.4 nl was injected into the cell of the 1 cell stage embryo. Embryos were stage matched, anaesthetised using 0.016 % MESAB in E3 embryo medium and imaged with a Nikon SMZ800 stereo microscope and Canon Eos 1100D camera and software.
Quantitative and RT-PCR analysis of splicing morpholino injected embryos
Total RNA was extracted using TRIZOL according to manufactures’ instructions. For RT-PCR RNA was extracted from 50 pooled Long-pec AB* WT embryos either uninjected or injected with either 100 μM control splice 1 mismatch morpholino or 100 μM flcn splice 1 morpholino. For qPCR RNA was extracted from 30 pooled Prim5 Nacre WT embryos either uninjected or injected with either 100 μM of control splice 2 mismatch morpholino or 100 μM of flcn splice 2 morpholino. Following DNAse treatment using TURBO DNA-free™ Kit (Ambion AM1907), 1 μg total RNA was subjected to RT-PCR using M-MLV Reverse Transcriptase, RNase H Minus, Point Mutant (Promega M3681, M3682, M3683) and random hexamers (Fermentas SO142). Specific PCR was performed on 5 μl cDNA using GoTaq polymerase (Promega M300) with the following primers: Splice 1 forward 5’TCCCATCACATGACACACAA3’, Splice 1 reverse 5’GAGACTGCGCACACATGC3’ Exon 1 forward 5’CCCATCACATGACACACAAGA3’, Exon 1 reverse 5’TGCATGTTCATCAGCCTTCC3’, Exon 2 forward 5’TGTTTTGTACCGAAGCCGTC3’, Exon 2 reverse 5’CTCACACATGTCCGCTCTCT3’, Exon 1 to Intron Forward 5’TTGCCATTGTCAACGCTTTG3’ and Exon 1 to Intron Reverse 5’AACTGCGAGATATCCAGCCA3’ and the following parameters 95 °C for 5 min, (95 °C for 30 s → 58 °C for 30 s → 72 °C for 30 s) x34 and 72 °C for 5 min. qPCR was performed using FastStart Universal SYBR Green Master (Roche) and primers Exon 1 Forward, Exon 1 Reverse, Exon 2 Forward and Exon 2 Reverse above and control gene primers β-Actin F CGAGCTGTCTTCCCATCCA, β-Actin R TCACCAACGTAGCTGTCTTTCTG, Rpl13α F TCTGGAGGACTGTAAGAGGTATGC and Rpl13α R AGACGCACAATCTTGAGAGCAG  on a Stratagene Mx3000P qPCR machine.
Whole-mount in situ hybridization, double fluorescent whole-mount in situ hybridization and immunohistochemistry
Whole-mount in situ hybridization was carried out as described by Thisse and Thisse . The flcn gene probe was transcribed directly from cloned cDNA in pCS2+, linearized with NotI (NEB) and transcribed with T7 polymerase (Promega) using Dig labelling mix (Roche). The pcna gene probe was transcribed from a plasmid gifted by F. van Eeden, linearized with Not1 (NEB) and transcribed with SP6 (Promega) using Dig labelling mix (Roche). Pre tail bud embryos were fixed in 4 % paraformaldehyde overnight and then manually dechorionated. Embryos older than Prim-5 stage were manually dechorionated and then fixed in 4 % paraformaldehyde overnight at 4 °C. Embryos were cleared in glycerol and imaged with a Nikon SMZ800 stereo microscope and Canon Eos 1100D camera and software and a Zeiss Axioplan compound microscope and Axiovision software.
Double fluorescent whole mount in situ hybridization was carried out as described by . The flcn gene probe was transcribed as described above and the pcna gene probe was transcribed as described above but labelled with DNP labelling kit (Perkin Elmer). pcna was developed with anti-DNP-POD (1:300) and revealed with Tyramide-Cy3 at a concentration of 1/25. flcn was developed with anti-DIG-POD (1:300) and revealed with Tyramide-FITC at a concentration of 1/150
Whole mount immunofluorescence of acetylated alpha tublin was carried out as follows: embryos were injected with either flcn splice 2 morpholino or mismatch morpholino, stage-matched at Prim-5 stage, manually dechorionated and fixed in 4 % formaldehyde in PBS overnight at 4 °C. Embryos were permeabilised using 10 μg/ml proteinase K (Sigma P4850) in PBS containing 0.1 % tween (PBST) for 20 min at room temperature, followed by re-fixation for 20 min. Non-specific binding sites were blocked by incubation in 10 % newborn calf serum in PBS for 1 h at room temperature, followed by overnight incubation at 4 °C in mouse acetylated alpha tublin antibodies (abcam 24610) diluted 1:200 in block buffer. The following day embryos were washed 4x 30 min in PBST and re-blocked for 1 h. Embryos were incubated in sheep-anti-mouse cy3-conjugated antibodies (Sigma C2181) diluted 1:200 in block buffer overnight at 4 °C, followed by 4x 30 min wash in PBST. Embryos were refixed in 4 % formaldehyde in PBS for 20 min at room temperature. Embryos were mounted in Vectashield (Vector labs) and images captured using a Leica LSI TCS Zoom confocal and Leica software. Images were processed using ImageJ software (NIH).
Acridine orange staining
Twenty four hour post fertilisation Nacre embryos were incubated in 1 μg/ml Acridine Orange (Sigma A6014) in E3 for 30 min. Embryos were washed, anaesthetised and imaged under fluorescent light on a Leica MZ10 Stereo microscope with Jenoptik ProgRes camera.
Embryo preparation and Confocal Imaging
For confocal imaging of zebrafish Fucci embryos  were injected with either 100 μM flcn splice 2 mismatch morpholino or 100 μM flcn splice 2 morpholino. Embryos were maintained at 28 °C in Danieaus solution (58 mM NaCl, 0.7 mM KCl, 0.4 mM MgSO4, 0.6 mM Ca(NO3)2, 5.0 mM HEPES pH 7.6) until bud-stage. A 96-well plate (Greiner Bio-One) was prepared by adding 50 μl of 1 % agarose (Bioline, cat. no. BIO-41025) in Danieaus solution covered by 30 μl of ethanol (VWR) to each well and inserting a brass 96 pin spotter (diameter 0.7 mm) into the plate until the agarose set when the spotter was removed and the plate washed thoroughly with distilled water. Embryos were dechorinated on 1 % agarose in Danieaus solution and transferred with 60 μl of Danieaus in to the 96 well agarose plate. Embryos were oriented using a Microlance 3 needle outer diameter 0.51 mm (BD Biosciences). 40 μl of mineral oil (Sigma, cat. no. M5904) was added over the surface of the Danieaus solution. Leica LSI TCS Zoom confocal and software was used to record the time-lapse video. Time-lapse was conducted in xyzt mode using the 532 nm and 488 nm lasers. 40 slices were taken for each stack for 11 h, z volume 358.395 μm lens 5x, optical zoom 1.362 step size 9.19, red gain 1200 offset -4, frame average 4 green gain 1250, offset -5 frame average 8. Frames were taken at 22 min 36 s intervals from bud stage for flcn splice 2 mismatch morpholino and 2-somite stage for flcn splice 2 morpholino injected embryos. For analysis maximum projection images were prepared, pixel intensities were manually identified in the selected area for all embryos. For still imaging, 3 slices were taken for each image with a frame average 3, 2x lens, z volume 500 μm maximum projections were produced and analysed.
Confocal imaging of Tg(wt1b:EGFP)
Tg(wt1b:EGFP) embryos injected with flcn splice 2 morpholino or flcn splice 2 mismatch morpholino were stage matched anaesthetised and imaged on a Leica TCS LSI confocal. Images were processed using ImageJ.
Analysis of zFucci reporter activity was carried out by statistical evaluation of pixel intensities of fluorescence signals detected for red and green reporters in flcn splice 2 morpholino injected embryos, compared to mismatch controls where a mean of 25 embryos were analysed from each set of injections for 4 independent sets of injections. Mean log pixel intensity was calculated for each imaged embryo. This was performed for red and green channels separately. We tested the hypothesis of differential intensity associated with either flcn MO or control MO using a linear mixed-effects modelling framework . Morpholino treatment was fitted as a categorical fixed effect and other sources of experimental variation (day, parental fish) were modelled as multilevel random effects. Treatment contrasts were tested using t-tests and results were validated by generating parameter confidence distributions using MCMC with 104 iterations. (R package: lme4. )
Time course analysis
Log-transformed intensity data from 3 flcn mismatch MO injected and 4 flcn splice 2 MO injected embryos was modelled using Generalised Additive Models, within a mixed-effects modelling framework. Cubic smoothing splines were used to model mean trajectories for each morpholino treatment. Following initial visual inspection, biological variation was modelled with separate parameters for each morpholino. A first order autoregressive process was used to model temporal autocorrelation within each time course. (R package: mgcv. )
Following replicate (n = 3, mean number embryos per group = 80) sets of injections with different RNA concentrations; phenotypes were categorised as ‘severe’ or ‘mild’ hydrocephalus, or ‘normal’. The proportion of the population showing each of these phenotypes with increasing inject RNA concentration was modelling using ordinal logistic regression (R package: MASS. ).
Quantitative PCR data from morpholino-injected embryos were normalized against mean expression of two housekeeping genes (β-actin and Rpl13α ). Normalized expression data were analysed using a linear mixed effects model. Morpholino (morpholino, MO, vs. mismatch, MM) and exon (1 versus 2) were fitted as interacting fixed effects; with technical replication (n = 3), explicitly nested within biological replication (n = 3), fitted as random effects. Fixed effects were assessed using likelihood ratio tests. Model residuals were found to meet with assumptions of Normality and homogeneity of variance..
All statistical analysis was performed using R, version 3.02 
Animal ethics statement
Animals were raised as per Home Office guidelines and all experiments were in accordance with Home Office Animals (Scientific Procedures) Act 1986.
days post fertilization
hour post fertilization
loss of function
morpholino antisense oligonucleotides
We would like to thank C. Englert for the transgenic line Tg(wt1b:EGFP), Atsushi Miyawaki for zFucci lines and F. van Eeden for the PCNA plasmid.
Funding for this project was provided by the Myrovlytis Trust and a KWF travel grant to Monique Luijten.
EJK drafted the manuscript, participated in experimental design, conducted WMISH, conducted morpholino injections and rescue experiments, produced figures. MNHL helped draft the manuscript, conducted cilia experiments, conducted morpholino injections and rescue experiments, produced figures. HG participated in experimental design, conducted cell cycle experiments, produced figures. NL produced pilot data and contributed to experimental analysis. MR participated in experimental design, produced pilot data. JCB participated in the design of experiments, performed the statistical analysis and reviewed the manuscript. YH conducted bioinformatics analysis. MAMvS and EM helped conceive study and participated in its design and coordination. FM conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Toro JR, Wei MH, Glenn GM, Weinreich M, Toure O, Vocke C, Vocke C, Turner M, Choyke P, Merino MJ, Pinto PA, et al. BHD mutations, clinical and molecular genetic investigations of Birt-Hogg-Dube syndrome: a new series of 50 families and a review of published reports. J Med Genet. 2008;45(6):321–31.
- Menko FH, van Steensel MA, Giraud S, Friis-Hansen L, Richard S, Ungari S, Nordenskjold M, Hansen TV, Solly J, Maher ER, et al. Birt-Hogg-Dube syndrome: diagnosis and management. Lancet Oncol. 2009;10(12):1199–206.View ArticlePubMedGoogle Scholar
- Schmidt LS, Nickerson ML, Warren MB, Glenn GM, Toro JR, Merino MJ, Turner ML, Choyke PL, Sharma N, Peterson J, et al. Germline BHD-mutation spectrum and phenotype analysis of a large cohort of families with Birt-Hogg-Dube syndrome. Am J Hum Genet. 2005;76(6):1023–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Nickerson ML, Warren MB, Toro JR, Matrosova V, Glenn G, Turner ML, Duray P, Merino M, Choyke P, Pavlovich CP, et al. Mutations in a novel gene lead to kidney tumors, lung wall defects, and benign tumors of the hair follicle in patients with the Birt-Hogg-Dube syndrome. Cancer Cell. 2002;2(2):157–64.View ArticlePubMedGoogle Scholar
- Nookala RK, Langemeyer L, Pacitto A, Ochoa-Montano B, Donaldson JC, Blaszczyk BK, Chirgadze DY, Barr FA, Bazan JF, Blundell TL. Crystal structure of folliculin reveals a hidDENN function in genetically inherited renal cancer. Open Biol. 2012;2(8):120071.View ArticlePubMedPubMed CentralGoogle Scholar
- Wei MH, Blake PW, Shevchenko J, Toro JR. The folliculin mutation database: an online database of mutations associated with Birt-Hogg-Dube syndrome. Hum Mutat. 2009;30(9):E880–890.View ArticlePubMedPubMed CentralGoogle Scholar
- Luijten MNH, Basten SG, Claessens T, Vernooij M, Scott CL, Janssen R, Easton JA, Kamps MAF, Vreeburg M, Broers JLV, et al. Birt-Hogg-Dube syndrome is a novel ciliopathy. Hum Mol Genet. 2013;22(21):4383–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Baba M, Hong SB, Sharma N, Warren MB, Nickerson ML, Iwamatsu A, Esposito D, Gillette WK, Hopkins RF, 3rd, Hartley JL, Li Y, Khoo SK, Huang D, Qian CN, et al. Folliculin encoded by the BHD gene interacts with a binding protein, FNIP1, and AMPK, and is involved in AMPK and mTOR signaling. Proc Natl Acad Sci U S A. 2006;103(42):15552–7.
- Chen J, Futami K, Petillo D, Peng J, Wang P, Knol J, et al. Deficiency of FLCN in mouse kidney led to development of polycystic kidneys and renal neoplasia. PLoS One. 2008;3(10):e3581.
- Hartman TR, Nicolas E, Klein-Szanto A, Al-Saleem T, Cash TP, Simon MC, Henske EP. The role of the Birt-Hogg-Dube protein in mTOR activation and renal tumorigenesis. Oncogene. 2009;28(13):1594–604.View ArticlePubMedPubMed CentralGoogle Scholar
- Hudon V, Sabourin S, Dydensborg AB, Kottis V, Ghazi A, Paquet M, Crosby K, Pomerleau V, Uetani N, Pause A. Renal tumour suppressor function of the Birt-Hogg-Dube syndrome gene product folliculin. J Med Genet. 2010;47(3):182–9.View ArticlePubMedGoogle Scholar
- Preston RS, Philp A, Claessens T, Gijezen L, Dydensborg AB, Dunlop EA, Harper KT, Brinkhuizen T, Menko FH, Davies DM, et al. Absence of the Birt-Hogg-Dube gene product is associated with increased hypoxia-inducible factor transcriptional activity and a loss of metabolic flexibility. Oncogene. 2011;30(10):1159–73.View ArticlePubMedGoogle Scholar
- Hong SB, Oh H, Valera VA, Baba M, Schmidt LS, Linehan WM. Inactivation of the FLCN tumor suppressor gene induces TFE3 transcriptional activity by increasing its nuclear localization. PLoS One. 2010;5(12):e15793.View ArticlePubMedPubMed CentralGoogle Scholar
- Betschinger J, Nichols J, Dietmann S, Corrin PD, Paddison PJ, Smith A. Exit from pluripotency is gated by intracellular redistribution of the bHLH transcription factor Tfe3. Cell. 2013;153(2):335–47.View ArticlePubMedPubMed CentralGoogle Scholar
- Petit CS, Roczniak-Ferguson A, Ferguson SM. Recruitment of folliculin to lysosomes supports the amino acid-dependent activation of Rag GTPases. J Cell Biol. 2013;202(7):1107–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsun ZY, Bar-Peled L, Chantranupong L, Zoncu R, Wang T, Kim C, Spooner E, Sabatini DM. The folliculin tumor suppressor is a GAP for the RagC/D GTPases that signal amino acid levels to mTORC1. Mol Cell. 2013;52(4):495–505.View ArticlePubMedGoogle Scholar
- Possik E, Jalali Z, Nouet Y, Yan M, Gingras MC, Schmeisser K, Panaite L, Dupuy F, Kharitidi D, Chotard L, et al. Folliculin regulates ampk-dependent autophagy and metabolic stress survival. PLoS Genet. 2014;10(4):e1004273.View ArticlePubMedPubMed CentralGoogle Scholar
- Hasumi Y, Baba M, Ajima R, Hasumi H, Valera VA, Klein ME, Haines DC, Merino MJ, Hong SB, Yamaguchi TP, et al. Homozygous loss of BHD causes early embryonic lethality and kidney tumor development with activation of mTORC1 and mTORC2. Proc Natl Acad Sci U S A. 2009;106(44):18722–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Berghmans S, Jette C, Langenau D, Hsu K, Stewart R, Look T, Kanki JP. Making waves in cancer research: new models in the zebrafish. Biotechniques. 2005;39(2):227–37.View ArticlePubMedGoogle Scholar
- Nepal C, Hadzhiev Y, Previti C, Haberle V, Li N, Takahashi H, Suzuki AM, Sheng Y, Abdelhamid RF, Anand S, et al. Dynamic regulation of the transcription initiation landscape at single nucleotide resolution during vertebrate embryogenesis. Genome Res. 2013;23(11):1938–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Perner B, Englert C, Bollig F. The Wilms tumor genes wt1a and wt1b control different steps during formation of the zebrafish pronephros. Dev Biol. 2007;309(1):87–96.View ArticlePubMedGoogle Scholar
- Pathak N, Obara T, Mangos S, Liu Y, Drummond IA. The zebrafish fleer gene encodes an essential regulator of cilia tubulin polyglutamylation. Mol Biol Cell. 2007;18(11):4353–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Kramer-Zucker AG, Olale F, Haycraft CJ, Yoder BK, Schier AF, Drummond IA. Cilia-driven fluid flow in the zebrafish pronephros, brain and Kupffer's vesicle is required for normal organogenesis. Development. 2005;132(8):1907–21.View ArticlePubMedGoogle Scholar
- Kawai A, Kobayashi T, Hino O. Folliculin regulates cyclin D1 expression through cis-acting elements in the 3 ' untranslated region of cyclin D1 mRNA. Int J Oncol. 2013;42(5):1597–604.PubMedGoogle Scholar
- Laviolette LA, Wilson J, Koller J, Neil C, Hulick P, Rejtar T, Karger B, Teh BT, Iliopoulos O. Human folliculin delays cell cycle progression through late S and G2/M-phases: effect of phosphorylation and tumor associated mutations. PLoS One. 2013;8(7):e66775.View ArticlePubMedPubMed CentralGoogle Scholar
- Sugiyama M, Sakaue-Sawano A, Iimura T, Fukami K, Kitaguchi T, Kawakami K, Okamoto H, Higashijima SI, Miyawaki A: Illuminating cell-cycle progression in the developing zebrafish embryo. Proc Natl Acad Sci U S A 2009;106(49):20812-7.
- Okimoto K, Sakurai J, Kobayashi T, Mitani H, Hirayama Y, Nickerson ML, Warren MB, Zbar B, Schmidt LS, Hino O. A germ-line insertion in the Birt-Hogg-Dube (BHD) gene gives rise to the Nihon rat model of inherited renal cancer. Proc Natl Acad Sci U S A. 2004;101(7):2023–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Lingaas F, Comstock KE, Kirkness EF, Sorensen A, Aarskaug T, Hitte C, Nickerson ML, Moe L, Schmidt LS, Thomas R, et al. A mutation in the canine BHD gene is associated with hereditary multifocal renal cystadenocarcinoma and nodular dermatofibrosis in the German Shepherd dog. Hum Mol Genet. 2003;12(23):3043–53.View ArticlePubMedGoogle Scholar
- Ferrante MI, Romio L, Castro S, Collins JE, Goulding DA, Stemple DL, Woolf AS, Wilson SW. Convergent extension movements and ciliary function are mediated by ofd1, a zebrafish orthologue of the human oral-facial-digital type 1 syndrome gene. Hum Mol Genet. 2009;18(2):289–303.View ArticlePubMedGoogle Scholar
- Hellman NE, Liu Y, Merkel E, Austin C, Le Corre S, Beier DR, Sun Z, Sharma N, Yoder BK, Drummond IA. The zebrafish foxj1a transcription factor regulates cilia function in response to injury and epithelial stretch. Proc Natl Acad Sci U S A. 2010;107(43):18499–504.View ArticlePubMedPubMed CentralGoogle Scholar
- Brody SL, Yan XH, Wuerffel MK, Song SK, Shapiro SD. Ciliogenesis and left-right axis defects in forkhead factor HFH-4-null mice. Am J Respir Cell Mol Biol. 2000;23(1):45–51.View ArticlePubMedGoogle Scholar
- Hagenlocher C, Walentek P, Ller MC, Thumberger T, Feistel K. Ciliogenesis and cerebrospinal fluid flow in the developing Xenopus brain are regulated by foxj1. Cilia. 2013;2(1):12.View ArticlePubMedPubMed CentralGoogle Scholar
- Guemez-Gamboa A, Coufal NG, Gleeson JG. Primary cilia in the developing and mature brain. Neuron. 2014;82(3):511–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Plotnikova OV, Pugacheva EN, Golemis EA. Primary cilia and the cell cycle. Methods Cell Biol. 2009;94:137–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Fu M, Wang C, Li Z, Sakamaki T, Pestell RG. Minireview: Cyclin D1: normal and abnormal functions. Endocrinology. 2004;145(12):5439–47.View ArticlePubMedGoogle Scholar
- Valk E, Loog M. Multiple Pho85-dependent mechanisms control G1 cyclin abundance in response to nutrient stress. Mol Cell Biol. 2013;33(7):1270–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Yan M, Gingras MC, Dunlop EA, Nouet Y, Dupuy F, Jalali Z, Possik E, Coull BJ, Kharitidi D, Dydensborg AB, et al. The tumor suppressor folliculin regulates AMPK-dependent metabolic transformation. J Clin Invest. 2014;124(6):2640–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu W, Chen Z, Ma Y, Wu X, Jin Y, Hou S. Genetic characterization of the Drosophila birt-hogg-dube syndrome gene. PLoS One. 2013;8(6):e65869.View ArticlePubMedPubMed CentralGoogle Scholar
- Landsverk ML, Weiser DC, Hannibal MC, Kimelman D. Alternative splicing of sept9a and sept9b in zebrafish produces multiple mRNA transcripts expressed throughout development. PLoS One. 2010;5(5):e10712.View ArticlePubMedPubMed CentralGoogle Scholar
- Lumsden AL, Henshall TL, Dayan S, Lardelli MT, Richards RI. Huntingtin-deficient zebrafish exhibit defects in iron utilization and development. Hum Mol Genet. 2007;16(16):1905–20.View ArticlePubMedGoogle Scholar
- Raldua D, Andre M, Babin PJ. Clofibrate and gemfibrozil induce an embryonic malabsorption syndrome in zebrafish. Toxicol Appl Pharmacol. 2008;228(3):301–14.View ArticlePubMedGoogle Scholar
- Kok FO, Shin M, Ni CW, Gupta A, Grosse AS, van Impel A, Kirchmaier BC, Peterson-Maduro J, Kourkoulis G, Male I, et al. Reverse genetic screening reveals poor correlation between morpholino-induced and mutant phenotypes in zebrafish. Dev Cell. 2015;32(1):97–108.View ArticlePubMedGoogle Scholar
- Law SH, Sargent TD. The serine-threonine protein kinase PAK4 is dispensable in zebrafish: identification of a morpholino-generated pseudophenotype. PLoS One. 2014;9(6):e100268.View ArticlePubMedPubMed CentralGoogle Scholar
- Rossi A, Kontarakis Z, Gerri C, Nolte H, Holper S, Kruger M, Stainier DY: Genetic compensation induced by deleterious mutations but not gene knockdowns. Nature 2015;524:230–233
- Eisen JS, Smith JC. Controlling morpholino experiments: don't stop making antisense. Development. 2008;135(10):1735–43.View ArticlePubMedGoogle Scholar
- Gaur K, Li J, Wang D, Dutta P, Yan SJ, Tsurumi A, et al. The Birt-Hogg-Dube tumor suppressor Folliculin negatively regulates ribosomal RNA synthesis. Hum Mol Genet. 2013;22(2):284–99.
- Hasumi H, Baba M, Hong SB, Hasumi Y, Huang Y, Yao M, Valera VA, Linehan WM, Schmidt LS. Identification and characterization of a novel folliculin-interacting protein FNIP2. Gene. 2008;415(1-2):60–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Medvetz DA, Khabibullin D, Hariharan V, Ongusaha PP, Goncharova EA, Schlechter T, Darling TN, Hofmann I, Krymskaya VP, Liao JK, et al. Folliculin, the product of the Birt-Hogg-Dube tumor suppressor gene, interacts with the adherens junction protein p0071 to regulate cell-cell adhesion. PLoS One. 2012;7(11):e47842.View ArticlePubMedPubMed CentralGoogle Scholar
- Nahorski MS, Seabra L, Straatman-Iwanowska A, Wingenfeld A, Reiman A, Lu X, Klomp JA, Teh BT, Hatzfeld M, Gissen P, et al. Folliculin interacts with p0071 (plakophilin-4) and deficiency is associated with disordered RhoA signalling, epithelial polarization and cytokinesis. Hum Mol Genet. 2012;21(24):5268–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Singh SR, Zhen W, Zheng Z, Wang H, Oh SW, Liu W, Zbar B, Schmidt LS, Hou SX. The Drosophila homolog of the human tumor suppressor gene BHD interacts with the JAK-STAT and Dpp signaling pathways in regulating male germline stem cell maintenance. Oncogene. 2006;25(44):5933–41.View ArticlePubMedGoogle Scholar
- Takagi Y, Kobayashi T, Shiono M, Wang L, Piao X, Sun G, Zhang D, Abe M,Hagiwara Y, Takahashi K, et al. Interaction of folliculin (Birt-Hogg-Dube gene product) with a novel Fnip1-like (FnipL/Fnip2) protein. Oncogene. 2008;27(40):5339–47.View ArticlePubMedGoogle Scholar
- Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF. Stages of embryonic development of the zebrafish. Dev Dyn. 1995;203(3):253–310.View ArticlePubMedGoogle Scholar
- Westerfield M. The Zebrafish Book. A Guide for the Laboratory Use of Zebrafish (Danio rerio), 4th edn. Eugene: University of Oregon Press; 2000.Google Scholar
- Tang R, Dodd A, Lai D, McNabb WC, Love DR. Validation of zebrafish (Danio rerio) reference genes for quantitative real-time RT-PCR normalization. Acta Biochim Biophys Sin. 2007;39(5):384–90.View ArticlePubMedGoogle Scholar
- Thisse C, Thisse B. High-resolution in situ hybridization to whole-mount zebrafish embryos. Nat Protoc. 2008;3(1):59–69.View ArticlePubMedGoogle Scholar
- Manfroid I, Delporte F, Baudhuin A, Motte P, Neumann CJ, Voz ML, Martial JA, Peers B. Reciprocal endoderm-mesoderm interactions mediated by fgf24 and fgf10 govern pancreas development. Development. 2007;134(22):4011–21.View ArticlePubMedGoogle Scholar
- Sugiyama M, Sakaue-Sawano A, Iimura T, Fukami K, Kitaguchi T, Kawakami K, Okamoto H, Higashijima S, Miyawaki A. Illuminating cell-cycle progression in the developing zebrafish embryo. Proc Natl Acad Sci U S A. 2009;106(49):20812–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Pinheiro JC, Bates DM: Mixed-effects models in S and S-PLUS: New York ; [London] : Springer, c2000 (2009 printing).
- Bates DM, M.; Bolker, B.; Walker, S.: lme4: Linear mixed-effects models using Eigen and S4. Journal of Statistical Software 2014.
- Wood SN. Stable and efficient multiple smoothing parameter estimation for generalized additive models. J Am Stat Assoc. 2004;99(467):673–86.View ArticleGoogle Scholar
- Venables WN, Ripley BD, Venables WN, Masw SP. Modern applied statistics with S. 4th ed. New York: Springer; 2002.View ArticleGoogle Scholar
- R: A language and environment for statistical computing. R Foundation for Statistical Computing [http://www.R-project.org/]