- Research article
- Open Access
Whole organism transcriptome analysis of zebrafish models of Bardet-Biedl Syndrome and Alström Syndrome provides mechanistic insight into shared and divergent phenotypes
© Hostelley et al. 2016
- Received: 15 December 2015
- Accepted: 28 April 2016
- Published: 3 May 2016
Bardet-Biedl Syndrome (BBS) and Alström Syndrome are two pleiotropic ciliopathies with significant phenotypic overlap between them across many tissues. Although BBS and Alström genes are necessary for the proper function of primary cilia, their role in defects across multiple organ systems is unclear.
To provide insight into the pathways underlying BBS and Alström phenotypes, we carried out whole organism transcriptome analysis by RNA sequencing in established zebrafish models of the syndromes.
We analyzed all genes that were significantly differentially expressed and found enrichment of phenotypically significant pathways in both models. These included multiple pathways shared between the two disease models as well as those unique to each model. Notably, we identified significant downregulation of genes in pathways relevant to visual system deficits and obesity in both disorders, consistent with those shared phenotypes. In contrast, neuronal pathways were significantly downregulated only in the BBS model but not in the Alström model. Our observations also suggested an important role for G-protein couple receptor and calcium signaling defects in both models.
Pathway network analyses of both models indicate that visual system defects may be driven by genetic mechanisms independent of other phenotypes whereas the majority of other phenotypes are a result of genetic players that contribute to multiple pathways simultaneously. Additionally, examination of genes differentially expressed in opposing directions between the two models suggest a deficit in pancreatic function in the Alström model, that is not present in the BBS model.
These findings provide important novel insight into shared and divergent phenotypes between two similar but distinct genetic syndromes.
- Primary Cilium
- Downregulated Gene
- Enrich Pathway
- Retinal Degeneration
- Neuronal Pathway
Many Mendelian disorders are pleiotropic, characterized by dysfunction across multiple organ systems. Deciphering the underlying molecular etiology of such diseases can be difficult due to the complexity of interpreting the role of a single gene in multiple disparate cell types. This complexity remains even for syndromes in which the subcellular defect is relatively well understood. For example, the ciliopathies are caused by mutations in genes encoding proteins associated with primary cilia or associated structures. However, the molecular pathways disrupted by loss of function in various ciliopathy genes remain relatively unclear. Two ciliopathies in particular, Bardet-Biedl Syndrome (BBS) and Alström Syndrome represent pleiotropic syndromes with overlapping features, but for which the primary molecular defects underlying those features are not well understood.
Alström Syndrome is a rare autosomal recessive, monogenic disorder caused by mutations in the ALMS1 gene. The ALMS1 gene is located on chromosome 2p13 and consists of 23 exons, which encode for a 461-kDa protein [1, 2]. Over 100 different disease-causing mutations have been identified in the ALMS1 gene, the majority of which are nonsense and frameshift mutations occurring in exons 8, 10, or 16 [1, 2]. The reported prevalence of Alström is less than one per million in the general population, with less than one thousand documented cases. Bardet-Biedl Syndrome (BBS) is also a rare autosomal recessive disorder, but exhibits significant genetic heterogeneity with 20 associated genes identified to date  The most common gene mutated in BBS is BBS1, located on chromosome 11q13, which accounts for approximately 25 % of BBS cases . The prevalence of BBS is 1 in 160,000 in the general population, but as high as 1 in 13,500 in a few isolated communities . Both Alström and BBS are classified as ciliopathies due to the localization of their proteins to the basal body of the primary cilia and the subsequent ciliary dysfunction that occurs when these proteins are disrupted. ALMS1 may be involved in intracellular trafficking and protein transport from the Golgi to the primary cilium . Several of the BBS proteins, including BBS1, form a protein complex called the BBSome, which is involved in the transport of molecules to the cilium and along the ciliary shaft [5, 6].
Given the potential similarities in intracellular localization and function, it is not surprising that many of the hallmark clinical features of Alström and BBS are shared between them. Both are uniquely characterized among the ciliopathies by highly penetrant early-onset obesity . They are also both characterized by progressive neurosensory deficits beginning in the first years of life that typically results in hearing loss in adolescence and blindness by adulthood [1, 2, 4, 7]. They exhibit similar endocrine and reproductive defects, such as hypogonadotropic hypogonadism as well as resulting infertility. Renal disease and respiratory illnesses are also features commonly found in both disorders [1, 4]. Despite these similarities, a number of differences between the disorders are present, particularly in the rates of type 2 diabetes mellitus (T2DM), developmental defects, and cognitive deficits. Alström is characterized by highly penetrant early onset T2DM with 70 % prevalence by the age of 20, while in BBS the prevalence is as low as 2–6 %. Developmental delays are common in BBS. They can also be present in Alström but with much lower penetrance, and are generally global, although it can be specific to certain areas of development. BBS is characterized by polydactyly, which can affect all four limbs or only the upper or lower limbs , but polydactyly is not present in Alström. Cognitive deficits present in BBS can vary from mild to severe with many patients requiring special schooling. Development of autism spectrum disorder and psychosis is also present . In contrast, the majority of Alström patients exhibit normal levels of intelligence although mild developmental delay of fine motor skills has been observed, likely as a result of sensory deficits .
Due to the high degree of phenotypic overlap between these disorders, they are often mistaken for one another, resulting in misdiagnosis. Here, we set out to characterize the global molecular profiles of each disorder in an effort to determine the underlying etiology distinguishing each and to identify common pathways that may be shared between both. We generated zebrafish models of both disorders by targeted knockdown approaches as previously described [8, 9]. Using whole transcriptome analysis by sequencing of RNA isolated from whole zebrafish larvae, we examined gene expression changes in each model and identified pathways altered as a result. Our findings revealed gene expression signatures that offer insight into the mechanisms driving unique and common features across the disorders. These include potentially common mechanisms underlying obesity and retinal degeneration, hallmarks of both disorders, and suggest genetic mechanisms for susceptibility to cognitive delay in BBS but not Alström. Taken together, these findings shed light onto the molecular etiology of two distinct but similar ciliopathies, providing important insight into the extent of similarity as well as important differences between them.
Identification of differentially expressed genes
To identify differentially expressed genes in Alström Syndrome (Alms) and Bardet-Biedl Syndrome (BBS), we generated zebrafish models of either syndrome using splice blocking morpholinos (MOs) targeting either alms1 or bbs1 transcripts that were previously validated to suppress protein production without introduction of off-target effects . MOs were injected into wild-type zebrafish embryos. At 48 h post fertilization (hpf) two replicates of RNA, from pools of 20 embryos per treatment as well as control embryos injected with control MO were collected for whole transcriptome sequencing. To fully capture all organ systems that may be impacted by disruption of either bbs1 or alms1, we isolated total RNA from whole zebrafish at 48 hpf, a time point at which most organ systems have developed and begin to be functional. MO-driven knockdown at 48 hpf was verified by western blot (Additional file 1: Figure S1.) RNA from each sample was sequenced to a depth of ~110 million reads with 91–100 million reads aligned to the zebrafish genome. Of these, between 85 and 88 % were mapped to exonic regions.
Gene expression changes in BBS model suggest deficiencies in neuronal and visual pathways
In total, 92 pathways were downregulated in the BBS model. Aside from the very broad grouping of “disease”, the most highly enriched pathway was the neuronal system, with 61 genes downregulated, accounting for 22.3 % of the genes involved in this pathway (Fig. 2a). This included genes such as grik1 (−3.74 log fold change, LFC), kcnv1 (−3.53 LFC), kcnk3 (−3.42 LFC), and slc17a7 (−3.35 LFC), all of which were significantly downregulated (Additional file 2: Table S1 and Additional file 3: Table S2). Other highly enriched neuronal pathways included transmission across chemical synapses, which contained 37 downregulated genes, neuroactive ligand-receptor interaction with 32 genes, and axon guidance with 24 genes (Fig. 2a). This also included specific neuronal functional pathways such as Glutamatergic and GABAergic pathways, containing 28 and 21 genes, respectively. In addition, an important neuronal pathway that was highly enriched among the downregulated pathways was defective ACTH and pro-opiomelanocortinin deficiency (POMCD), with 31 downregulated genes (Fig. 2a). Pro-opiomelanocortinin (POMC) is directly related to regulation of satiety and hyperphagia and is known to be disrupted in BBS [12–14]. It is processed through several enzymatic steps to produce melanocyte-stimulating hormone (MSH) and corticotrophin (ACTH) both of which are implicated in other forms of obesity such as melanocortin receptor mutations and Cushing’s [15, 16]. In all, pathways related to neuronal function accounted for approximately one-third of the downregulated pathways.
Perhaps unsurprisingly, some of the most highly enriched pathways were the visual pathways, consistent with the centrality of retinal degeneration to the BBS phenotype. These pathways included visual transduction and visual phototransduction, both with 19 genes downregulated. These included genes required for photoreceptor response, such as rhodopsin, rho (−4.67 LFC) and rhol (−3.17 LFC), guca1b (−4.18 LFC), rbp1 (−3.84 LFC) and several opsins, opn1mw2 (−3.74 LFC), opn1lw2 (−3.36 LFC), opn1sw2 (−2.98 LFC), and opn1sw1 (−2.83 LFC).
Several other pathways with general intracellular and signaling roles were also significantly downregulated in the BBS model. For example, the G protein coupled receptor (GPCR) ligand binding pathway was heavily effected with 31 genes downregulated as well as G alpha (i) signaling events with 22 genes, or 10 % of the genes involved in the pathway. This may be consistent with the important link between ciliary localization of GPCRs and the role of BBS proteins in ciliary trafficking or trafficking of GPCRs specifically  Calcium signaling was also enriched among the downregulated pathways. 11.4 % (19) of the genes involved in the calcium-signaling pathway were downregulated. The presence of several neuronal and visual pathways is consistent with observed phenotypes of the syndrome. BBS patients are characterized by retinal degeneration and cognitive impairment [2, 18], potentially as a result of disruption of neuronal signaling pathways . In addition, two pathways directly related to energy homeostasis were among the most enriched pathways in downregulated genes. These were related to insulin secretion and integration of energy metabolism. The former may be directly related to glucose homeostasis and evidence suggests that it may be directly perturbed in BBS [19, 20], but both pathways are likely relevant to overall defects in energy metabolism, consistent with the obesity phenotype that is a hallmark of BBS.
169 GO terms were enriched among the downregulated genes (Fig. 2b). The top downregulated GO terms were highly consistent with the top downregulated pathways. The most highly enriched GO terms included several related to the neuronal and visual systems such as: regulation of neuronal synaptic plasticity and regulation of synapse structure or activity, and phototransduction, detection of light stimulus, and visual perception. Several GPCR GO terms were also found, including termination of GPCR signaling pathway and negative regulation of GPCR protein signaling (Fig. 2b).
Although the number of genes upregulated in the BBS model was considerably lower than the number downregulated, these genes encompassed a much greater number of pathways. This is likely due to the smaller number of differentially expressed genes found in each upregulated pathway. In total, 170 pathways were upregulated in the BBS model. The majority of the upregulated pathways are those involved with the cell cycle, including mitosis and the S phase, and those involved with transcription, including several RNA polymerases (Fig. 2c). 159 GO terms were enriched among the upregulated genes and were consistent with the upregulated pathways (Fig. 2d). The majority were involved in RNA processing such as: exonucleolytic trimming involved in rRNA processing, maturation of 5.8S rRNA, tRNA methylation, rRNA processing, and ribosome biogenesis. These observations suggest a previously underappreciated role for BBS genes in regulation of cell cycle events, transcription and translation. Interestingly, the most highly enriched GO term among the upregulated genes is cilium-dependent cell motility consistent with the known roles of BBS genes in ciliogenesis and cilia function [5, 13].
Analysis of pathways disrupted in Alström model suggests importance of phototransduction signaling
Similar to the BBS model, many of the downregulated pathways involved the visual pathways including, diseases associated with visual transduction and visual phototransduction. Other enriched visual pathways include: inactivation, recovery and regulation of the phototransduction cascade and the phototransduction cascade, each with 7 downregulated genes, as well as, activation of the phototransduction cascade with 5 of the 12 involved genes (45 %) being downregulated. This is consistent with the retinal degeneration phenotype that is likewise characteristic of Alström [1, 2, 21]. Furthermore, the POMC pathway was also downregulated in alms1 as it was in bbs1, consistent with hyperphagia induced obesity observed in the disorder . Pathways involved with GPCR and calcium signaling were found to be downregulated in alms1 as they were in bbs1. However, in each of the cases mentioned, considerably fewer genes were found in the pathways in alms1 than were found in bbs1.
47 GO terms were enriched among the downregulated genes. Consistent with the downregulated pathways, the most highly enriched GO terms associated with downregulated genes were all involved with the visual system, including phototransduction, visual perception, and detection of light stimulus (Fig. 3b). Several GPCR related GO terms were also highly enriched, such as termination of GPCR signaling pathway and regulation of GPCR protein signaling pathway. Although these terms reflect a high degree of overlap between the BBS and Alström models, the neuronal pathways that were prevalent in bbs1-depleted larvae were noticeably absent from both the downregulated pathways and enriched GO terms in alms1-deficient animals. This may be consistent with absent or very mild cognitive impairment in Alström patients , contrasting with BBS.
Overlap of differential gene expression between BBS and Alström models
We then looked at the pathways and GO terms downregulated in both models to assess overlap. 24 pathways were downregulated in both models (Fig. 4b). In almost every case, the number of genes downregulated in the shared pathways was greater in the BBS model. The majority of genes differentially expressed in alms1-deficient animals were also differentially expressed in animals targeted for bbs1. Pathways associated with the visual system were the most commonly downregulated pathways in both models. The POMC pathway was also one of the top downregulated pathways found in both. Other pathways that were downregulated in both disorders include those associated with GPCRs and the calcium-signaling pathway.
47 GO terms were enriched among the downregulated genes found in both models. The 30 terms with the greatest enrichment in the downregulated genes found in both are shown in Fig. 4c. The enriched GO terms were consistent with the pathways, again, with the most highly enriched GO terms associated with the visual system including: phototransduction, detection of light stimulus, visual perception, and sensory perception of light stimulus. Several GPCR terms were also found to be downregulated in both, such as termination of GPCR signaling pathway, negative regulation of GPCR protein signaling, and GPCR signaling pathway.
Our examination of the overlap between these pathways allowed us to identify genes that exert a greater impact on the spectrum of phenotypes observed in these disorders based on their contribution to a larger number of pathways. For example, downregulation of the G protein genes gngt1, gng7, and gng8, as well as the GABA receptor genes gabrb2, gabrb3, gabra1, gabra2, gabra3, and gabrg2 in BBS may be driving the cognitive deficits characteristic of the disorder, as these genes were found in the majority of the downregulated neuronal system pathways. Each of these genes contributed to greater than 20 % of the pathways identified perhaps indicating a broad phenotypic spectrum driven by a relatively small number of genes. Indeed, gngt1 alone contributed to nearly 60 % of pathways, indicating its potential importance in contributing to multiple phenotypes (Additional file 4: Figure S2A). Additionally, downregulation of rhodopsin rho, gnat1, grm1, grm5, pde6b, and opn1mw2 in both Alström and BBS may be major drivers of the retinal degeneration that is characteristic of both disorders, as these genes were found in the majority of the downregulated visual and phototransduction pathways (Additional file 4: Figure S2A-B).
Gene expression changes provide insight into unique molecular features of BBS
117 GO terms were enriched among the unique to bbs1 downregulated genes, the GO terms with the greatest enrichment are shown in Fig. 6b. Again, most of the top GO terms were related to the neuronal system, including regulation of neuronal synaptic plasticity, synapse assembly, and synaptic transmission, supporting the evidence that neuronal system defects may be a characteristic marker of BBS that is not shared with Alström. Also absent in the unique enriched GO terms were those related to the visual system and GPCR signaling, similar to the unique downregulated pathways.
Opposite changes in expression
In this study, we have undertaken a comprehensive transcriptomic approach in models of BBS and Alström syndrome with the goal of identifying genetic pathway mechanisms underlying phenotypes associated with these disorders. Our use of the zebrafish model allowed for assessment of gene expression across the whole organism, an approach that is particularly relevant to multi-organ syndromes. Although our assessments were limited to the window of analysis to embryonic and early larval stages due to the use of transient knockdown to generate the disease models, the nature of these syndromes as developmental disorders suggests that these may be the most relevant time points for onset of most phenotypes associated with disease. Indeed, the majority of the hallmark features of either disorder manifest within the first or second decade of life, consistent with developmental or early childhood onset [1, 2, 4, 7]. Our global analysis revealed important insights into the pathways likely contributing to disease phenotypes. These included known mechanisms, such as the POMC deficiency reported in BBS models to underlie hyperphagia and obesity [13, 14], as well as perhaps novel mechanisms, such as the neuronal pathways identified that may relate to cognitive impairment. Taken together, the results of these analyses offer novel candidate genes and pathways for further exploration in the study of these disorders and related ciliopathies.
Our observations in the BBS model shed light onto the potential role of BBS genes as well as phenotypes driven by their dysfunction. For example, perhaps the most prominent feature of BBS is retinal degeneration characterized by rod-cone dystrophy [4, 19]. Previous reports have attributed this feature to the importance of the photoreceptor primary cilium and trafficking within it . Our findings suggest the specific perturbation of genes in the phototransduction pathway, such as rho (−4.67 LFC), gngt1 (−2.38 LFC), and pde6g (−1.80 LFC), or the visual transduction pathway, such as gnat1 (−3.32 LFC), cnga1 (−3.64 LFC), and pde6b (−1.71 LFC). It is unclear if this disruption is a cause or consequence of other primary factors driving retinal degeneration, but they potentially offer novel factors to examine in models of BBS or in patients. Interestingly, our analysis of the interconnectedness between pathways suggests little overlap between the disrupted vision pathways and neuronal signaling. The exception to this in the BBS model was the presence of the gngt1 gene, which contributed to both the visual and neuronal pathways. It is unlikely, however, that downregulation of this gene is sufficient for neuronal deficits given the absence of neuronal system effects in the Alström model. This potentially indicates the specificity of the retinal phenotype to that structure, rather than a general consequence of neuronal dysfunction. BBS is characterized by deficits in both functions and it is therefore unsurprising that pathways important to both tissues are perturbed with loss of bbs1. In contrast, a higher degree of connectedness was observed between neuronal pathways and energy metabolism regulation, including insulin secretion, potentially implicating neuronal regulation in those defects.
Other pathways disrupted in the BBS model offer novel insight into phenotypes associated with the disorder. For example, a number of basic cellular functions appear to be upregulated in the BBS model. This included cell cycle processes as well as transcriptional and translational regulation. Much evidence supports the role of ciliary proteins in regulation of the cell cycle, given the tight light between ciliogenesis and proliferation , but our observations perhaps suggest a more direct role for the BBS genes in regulation of transcription or translation. This may be consistent with the recently proposed role of BBS7 as a transcription factor .
Perhaps the most intriguing observation relates to the genes that were differentially regulated in each model but with divergent directionality. Only 8 genes fit this criteria, two of which were uncharacterized and therefore of unknown relevance. Interestingly, however, two of the eight genes were pancreatic proteases downregulated in the Alström model and upregulated in the BBS model. We also identified 3 other pancreatic protease genes that were downregulated in the Alström model but unchanged in the BBS model. Although these enzymes relate specifically to exocrine pancreatic function, they may offer a novel mechanism by which differences in either exocrine or endocrine pancreatic function occur between the two disorders. The prevalence of diabetes in Alström is far higher than that in BBS and is at least in part likely due to deficits in endocrine pancreas function . Our data might suggest a novel interaction between exocrine and endocrine pancreas. Moreover, given the importance of pancreatic enzymes in digestion, these observations might suggest a mechanism by which digestive phenotypes, a secondary feature of Alström [1, 7], might arise.
In summary, our analysis of zebrafish models of two related but distinct syndromes offers novel insight into the mechanisms underlying their phenotypes. These data confirm previously proposed mechanisms of dysfunction for the common phenotypes, such as obesity, while at the same time offering new mechanisms for many features of both syndromes. Importantly, these findings distinguish the two disorders from each other and perhaps suggest novel disruptions of cell types and pathways in each. Finally, these results inform the understanding of how BBS1, ALMS1 and cilia may function at the cellular level. These global gene expression data represent a starting point from which detailed mechanistic links can be explored.
Experiments were carried out in wild type Tubingen (WT) zebrafish. Adult zebrafish were housed and naturally mated according to standard protocol. All zebrafish work was conducted in accordance with University of Maryland IACUC guidelines.
Morpholino antisense oligonucleotides (MOs) that block splicing of targeted mRNAs were injected into one- to two- cell stage embryos. We used previously validated MOs to target bbs1 and alms1 transcripts . A control non-specific MO was used (5′-CCTCTTACCTCAGTTACAATTTATA-3′). The embryos were grown at 28.5 °C until 48 h post fertilization (hpf) for analyses.
RNA extraction and sequencing
Groups of 20 embryos were pooled for each condition (control MO, alms1 MO, bbs1 MO) at 48hpf for RNA extraction. RNA was extracted using Isol-RNA Lysis Reagent (5 Prime) and purified using Qiagen RNeasy cleanup kit according to manufacturer’s protocol. The three RNA samples were sent in duplicates for sequencing and primary analysis. 2 ug of total RNA per sample was used. Illumina RNAseq libraries were prepared using the TruSeq RNA Sample Prep kit (Illumina, San Diego, CA). Libraries were barcoded and pooled, and each received 0.33 of a lane of sequencing on an Illumina HiSeq 2500 with a paired-end 100 base configuration. Library construction and sequencing was performed at the Institute for Genome Sciences of the University of Maryland, School of Medicine (IGS).
RNA-sequencing data analysis
Sequenced reads were aligned to the zebrafish genomes (Zv10) using TopHat . Gene expression levels (read counts) were calculated using HTseq  based on ENSEMBL gene annotations . The zebrafish datasets were then normalized using reads per kilobase per million reads (RPKM). An average of 28,700 features were detected as expressed in each sample set. Replicate analysis and differentially expressed genes were determined using DEGseq . Both alms1 and bbs1 knockdowns were compared to the control, cutoffs of false discovery rate (FDR) <0.05 and p value <0.05 were used to identify genes with significant changes in expression. Primary data analysis was performed at the Institute for Genome Sciences of the University of Maryland, School of Medicine (IGS). Sequence data for samples is publicly available through the NCBI Sequence Read Archive (SRA), Submission ID “University of Maryland BBS Alstrom Zebrafish RNA-Seq”.
RNA was extracted from pools of 20 48hpf embryos using Isol-RNA Lysis Reagent (5 PRIME). The RNA was purified by 3 M sodium acetate and isopropanol precipitation, and spun for 10 min at 4° (12,500 rpm) following a 20 min incubation. The RNA pellet was then washed in 70 % ethanol and re-suspended in RNase free water. cDNA was transcribed using Fermentas First Strand cDNA transcription kit (Thermo Scientific) according to manufacturer’s protocol. cDNA was then diluted 1:3. qRT-PCR was performed using target specific primers and LightCycler 480 SybrGreen (Roche) according to manufacturer’s protocol on a LightCycler 480 machine (Roche). All samples were run in duplicate with CT value normalized to β-actin.
Western blotting was performed according to previously described protocol .
Determination of enriched pathways
Pathways enriched among the genes exhibiting changes in expression were determined by uploading sets of the differentially expressed genes to ConsensusPathDB (http://cpdb.molgen.mpg.de). Genes were identified based on their gene symbol (HGNC symbol). Pathways from the KEGG and Reactome databases were used, along with a q-value cutoff of <0.05 and a minimum of 2 genes overlap with the input list.
Generation of pathway networks
Pathway networks were generated using the ConsensusPathDB in the same manner as the enriched pathways. Once the enriched pathways were determined, the top 30 most highly significantly enriched pathways were selected for visualization. Pathways were shown as connected if there were at least 2 genes overlapping and the pathways shared at least 20 % of their genes.
Determination of enriched gene ontologies
Gene ontologies enriched among the differentially expressed genes were determined using the GO Enrichment Analysis tool of the Gene Ontology Consortium (http://geneontology.org/page/go-enrichment-analysis). Gene lists were uploaded to the tool and queried against the biological process Danio rerio gene ontology reference database, a cutoff p-value of <0.05 was applied.
This work was supported by R01DK102001 (N.A.Z.), P30DK090868 Baltimore PKD Center Pilot and Feasibility grant (N.A.Z.), and the University of Maryland Meyerhoff Graduate Fellows program (T.L.H.). The funding body played no role in design, collection or analysis of data, writing of the manuscript, or the decision to submit the manuscript for publication.
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.
- Marshall JD, Maffei P, Collin GB, Naggert JK. Alström syndrome: genetics and clinical overview. Curr Genomics. 2011;12:225–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Girard D, Petrovsky N. Alström syndrome: insights into the pathogenesis of metabolic disorders. Nat Rev Endocrinol. 2011;7:77–88.View ArticlePubMedGoogle Scholar
- Aldahmesh MA, Li Y, Alhashem A, et al. IFT27, encoding a small GTPase component of IFT particles, is mutated in a consanguineous family with Bardet-Biedl syndrome. Hum Mol Genet. 2014;23:3307–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Forsythe E, Beales PL. Bardet-Biedl syndrome. Eur J Hum Genet. 2013;21:8–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Nachury MV, Loktev AV, Zhang Q, et al. A core complex of BBS proteins cooperates with the GTPase Rab8 to promote ciliary membrane biogenesis. Cell. 2007;129:1201–13.View ArticlePubMedGoogle Scholar
- Su X, Driscoll K, Yao G, Raed A, Wu M, Beales PL, Zhou J. Bardet-Biedl syndrome proteins 1 and 3 regulate the ciliary trafficking of polycystic kidney disease 1 protein. Hum Mol Genet. 2014;23:5441–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Marshall JD, Paisey RB, Carey C, et al. Alström Syndrome. 2003 Feb 7 [Updated 2012 May 31]. In: Pagon RA, Adam MP, Ardinger HH, et al., editors. GeneReviews® [Internet]. Seattle (WA): University of Washington, Seattle; 1993-2016. Available from: http://www.ncbi.nlm.nih.gov/books/NBK1267/.
- Leitch CC, Lodh S, Prieto-Echagüe V, Badano JL, Zaghloul NA. Basal body proteins regulate Notch signaling through endosomal trafficking. J Cell Sci. 2014;127:2407–19.View ArticlePubMedPubMed CentralGoogle Scholar
- Lodh S, Hostelley TL, Leitch CC, O’Hare EA, Zaghloul NA. Differential effects on β-cell mass by disruption of Bardet-Biedl syndrome or Alstrom syndrome genes. Hum Mol Genet. 2015.Google Scholar
- Kamburov A, Stelzl U, Lehrach H, Herwig R. The ConsensusPathDB interaction database: 2013 update. Nucleic Acids Res. 2013;41:D793–800.View ArticlePubMedPubMed CentralGoogle Scholar
- The Gene Ontology Consortium. Gene Ontology Consortium: going forward. Nucleic Acids Res. 2014;43:D1049–56.View ArticlePubMed CentralGoogle Scholar
- Seo S, Guo D-F, Bugge K, Morgan DA, Rahmouni K, Sheffield VC. Requirement of Bardet-Biedl syndrome proteins for leptin receptor signaling. Hum Mol Genet. 2009;18:1323–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Mok CA, Héon E, Zhen M. Ciliary dysfunction and obesity. Clin Genet. 2010;77:18–27.View ArticlePubMedGoogle Scholar
- Guo D, Rahmouni K. Bardet-Biedl Syndrome (BBS) Proteins in POMC neurons are required for energy homeostasis. FASEB J. 2015;29:655.12.Google Scholar
- Eckstein N, Haas B, Hass M, Pfeifer V. Systemic therapy of Cushing’s syndrome. Orphanet J Rare Dis. 2014;9:122.View ArticlePubMedPubMed CentralGoogle Scholar
- Vaisse C, Clement K, Durand E, Hercberg S, Guy-Grand B, Froguel P. Melanocortin-4 receptor mutations are a frequent and heterogeneous cause of morbid obesity. J Clin Invest. 2000;106:253–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Loktev AV, Jackson PK. Neuropeptide Y family receptors traffic via the Bardet-Biedl syndrome pathway to signal in neuronal primary cilia. Cell Rep. 2013;5:1316–29.View ArticlePubMedGoogle Scholar
- Brinckman DD, Keppler-Noreuil KM, Blumhorst C, Biesecker LG, Sapp JC, Johnston JJ, Wiggs EA. Cognitive, sensory, and psychosocial characteristics in patients with Bardet-Biedl syndrome. Am J Med Genet A. 2013;161A:2964–71.View ArticlePubMedGoogle Scholar
- Gerdes JM, Christou-Savina S, Xiong Y, et al. Ciliary dysfunction impairs beta-cell insulin secretion and promotes development of type 2 diabetes in rodents. Nat Commun. 2014;5:5308.View ArticlePubMedGoogle Scholar
- Lee BH, Liu J, Wong D, Srinivasan S, Ashrafi K. Hyperactive neuroendocrine secretion causes size, feeding, and metabolic defects of C. elegans Bardet-Biedl syndrome mutants. PLoS Biol. 2011;9:e1001219.View ArticlePubMedPubMed CentralGoogle Scholar
- Manara R, Citton V, Maffei P, et al. Degeneration and plasticity of the optic pathway in Alström syndrome. AJNR Am J Neuroradiol. 2015;36:160–5.View ArticlePubMedGoogle Scholar
- Liu S, Young SM, Varisco BM. Dynamic expression of chymotrypsin-like elastase 1 over the course of murine lung development. Am J Physiol Lung Cell Mol Physiol. 2014;306:L1104–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Smyth MJ, O’Connor MD, Trapani JA. Granzymes: a variety of serine protease specificities encoded by genetically distinct subfamilies. J Leukoc Biol. 1996;60:555–62.PubMedGoogle Scholar
- Gullo L, Ventrucci M, Tomassetti P, Migliori M, Pezzilli R. Fecal elastase 1 determination in chronic pancreatitis. Dig Dis Sci. 1999;44:210–3.View ArticlePubMedGoogle Scholar
- Kukor Z, Toth M, Sahin-Toth M. Human anionic trypsinogen. Properties of autocatalytic activation and degradation and implications in pancreatic diseases. Eur J Biochem. 2003;270:2047–58.View ArticlePubMedGoogle Scholar
- Datta P, Allamargot C, Hudson JS, Andersen EK, Bhattarai S, Drack AV, Sheffield VC, Seo S. Accumulation of non-outer segment proteins in the outer segment underlies photoreceptor degeneration in Bardet-Biedl syndrome. Proc Natl Acad Sci U S A. 2015;112:E4400–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Gerdes JM, Davis EE, Katsanis N. The vertebrate primary cilium in development, homeostasis, and disease. Cell. 2009;137:32–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Gascue C, Tan PL, Cardenas-Rodriguez M, et al. Direct role of Bardet-Biedl syndrome proteins in transcriptional regulation. J Cell Sci. 2012;125:362–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu YP, Tsai I-C, Morleo M, et al. Ciliopathy proteins regulate paracrine signaling by modulating proteasomal degradation of mediators. J Clin Invest. 2014;124:2059–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Novarino G, Akizu N, Gleeson JG. Modeling human disease in humans: the ciliopathies. Cell. 2011;147:70–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Hearn T, Spalluto C, Phillips VJ, Renforth GL, Copin N, Hanley NA, Wilson DI. Subcellular localization of ALMS1 supports involvement of centrosome and basal body dysfunction in the pathogenesis of obesity, insulin resistance, and type 2 diabetes. Diabetes. 2005;54:1581–7.View ArticlePubMedGoogle Scholar
- Ansley SJ, Badano JL, Blacque OE, et al. Basal body dysfunction is a likely cause of pleiotropic Bardet-Biedl syndrome. Nature. 2003;425:628–33.View ArticlePubMedGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Anders S, Pyl PT, Huber W. HTSeq - a python framework to work with high-throughput sequencing data. Bioinformatics. 2014;31:166–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Flicek P, Amode MR, Barrell D, et al. Ensembl 2014. Nucleic Acids Res. 2014;42:D749–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26:136–8.View ArticlePubMedGoogle Scholar