Exocrine pancreas proteases regulate ß-cell proliferation in zebrafish ciliopathy models and in murine systems.

Pancreatic ß-cells are a critical cell type in the pathology of diabetes. Models of genetic syndromes featuring diabetes can provide novel mechanistic insights into regulation of ß-cells in the context of disease. We previously examined ß-cell mass in models of two ciliopathies, Alström Syndrome (AS) and Bardet-Biedl Syndrome (BBS), which are similar in the presence of metabolic phenotypes, including obesity, but exhibit strikingly different rates of diabetes. Zebrafish models of these disorders show deficient ß-cells with diabetes in AS models and an increased ß-cells absent diabetes in BBS models, indicating ß-cell generation or maintenance that correlates with disease prevalence. Using transcriptome analyses, differential expression of several exocrine pancreas proteases with directionality that was consistent with ß-cell numbers were identified. Based on these lines of evidence, we hypothesized that pancreatic proteases directly impact ß-cells. In the present study, we examined this possibility and found that pancreatic protease genes contribute to proper maintenance of normal ß-cell numbers, proliferation in larval zebrafish, and regulation of AS and BBS ß-cell phenotypes. Our data suggest that these proteins can be taken up directly by cultured ß-cells and ex vivo murine islets, inducing proliferation in both. Endogenous uptake of pancreatic proteases by ß-cells was confirmed in vivo using transgenic zebrafish and in intact murine pancreata. Taken together, these findings support a novel proliferative signaling role for exocrine pancreas proteases through interaction with endocrine ß-cells.

A BBS1 SVA F retrotransposon insertion is a frequent cause of Bardet-Biedl syndrome.

Bardet-Biedl syndrome (BBS) is a ciliopathy characterized by retinitis pigmentosa, obesity, polydactyly, cognitive impairment and renal failure. Pathogenic variants in 24 genes account for the molecular basis of >80% of cases. Toward saturated discovery of the mutational basis of the disorder, we carefully explored our cohorts and identified a hominid-specific SINE-R/VNTR/Alu type F (SVA-F) insertion in exon 13 of BBS1 in eight families. In six families, the repeat insertion was found in trans with c.1169 T > G, p.Met390Arg and in two families the insertion was found in addition to other recessive BBS loci. Whole genome sequencing, de novo assembly and SNP array analysis were performed to characterize the genomic event. This insertion is extremely rare in the general population (found in 8 alleles of 8 BBS cases but not in >10 800 control individuals from gnomAD-SV) and due to a founder effect. Its 2435 bp sequence contains hallmarks of LINE1 mediated retrotransposition. Functional studies with patient-derived cell lines confirmed that the BBS1 SVA-F is deleterious as evidenced by a significant depletion of both mRNA and protein levels. Such findings highlight the importance of dedicated bioinformatics pipelines to identify all types of variation.

ZRANB3 is an African-specific type 2 diabetes locus associated with beta-cell mass and insulin response.

Genome analysis of diverse human populations has contributed to the identification of novel genomic loci for diseases of major clinical and public health impact. Here, we report a genome-wide analysis of type 2 diabetes (T2D) in sub-Saharan Africans, an understudied ancestral group. We analyze ~18 million autosomal SNPs in 5,231 individuals from Nigeria, Ghana and Kenya. We identify a previously-unreported genome-wide significant locus: ZRANB3 (Zinc Finger RANBP2-Type Containing 3, lead SNP p = 2.831 × 10-9). Knockdown or genomic knockout of the zebrafish ortholog results in reduction in pancreatic ß-cell number which we demonstrate to be due to increased apoptosis in islets. siRNA transfection of murine Zranb3 in MIN6 ß-cells results in impaired insulin secretion in response to high glucose, implicating Zranb3 in ß-cell functional response to high glucose conditions. We also show transferability in our study of 32 established T2D loci. Our findings advance understanding of the genetics of T2D in non-European ancestry populations.

Genomic knockout of alms1 in zebrafish recapitulates Alström syndrome and provides insight into metabolic phenotypes.

Alström syndrome (OMIM #203800) is an autosomal recessive obesity ciliopathy caused by loss-of-function mutations in the ALMS1 gene. In addition to multi-organ dysfunction, such as cardiomyopathy, retinal degeneration and renal dysfunction, the disorder is characterized by high rates of obesity, insulin resistance and early-onset type 2 diabetes mellitus (T2DM). To investigate the underlying mechanisms of T2DM phenotypes, we generated a loss-of-function deletion of alms1 in the zebrafish. We demonstrate conservation of hallmark clinical characteristics alongside metabolic syndrome phenotypes, including a propensity for obesity and fatty livers, hyperinsulinemia and glucose response defects. Gene expression changes in ß-cells isolated from alms1-/- mutants revealed changes consistent with insulin hypersecretion and glucose sensing failure, which were corroborated in cultured murine ß-cells lacking Alms1. We also found evidence of defects in peripheral glucose uptake and concomitant hyperinsulinemia in the alms1-/- animals. We propose a model in which hyperinsulinemia is the primary and causative defect underlying generation of T2DM associated with alms1 deficiency. These observations support the alms1 loss-of-function zebrafish mutant as a monogenic model for mechanistic interrogation of T2DM phenotypes.

Survey of Human Chromosome 21 Gene Expression Effects on Early Development in Danio rerio.

Trisomy for human chromosome 21 (Hsa21) results in Down syndrome (DS), one of the most genetically complex conditions compatible with human survival. Assessment of the physiological consequences of dosage-driven overexpression of individual Hsa21 genes during early embryogenesis and the resulting contributions to DS pathology in mammals are not tractable in a systematic way. A recent study looked at loss-of-function of a subset of Caenorhabditis elegans orthologs of Hsa21 genes and identified ten candidates with behavioral phenotypes, but the equivalent over-expression experiment has not been done. We turned to zebrafish as a developmental model and, using a number of surrogate phenotypes, we screened Hsa21 genes for effects on early embyrogenesis. We prepared a library of 164 cDNAs of conserved protein coding genes, injected mRNA into early embryos and evaluated up to 5 days post-fertilization (dpf). Twenty-four genes produced a gross morphological phenotype, 11 of which could be reproduced reliably. Seven of these gave a phenotype consistent with down regulation of the sonic hedgehog (Shh) pathway; two showed defects indicative of defective neural crest migration; one resulted consistently in pericardial edema; and one was embryonic lethal. Combinatorial injections of multiple Hsa21 genes revealed both additive and compensatory effects, supporting the notion that complex genetic relationships underlie end phenotypes of trisomy that produce DS. Together, our data suggest that this system is useful in the genetic dissection of dosage-sensitive gene effects on early development and can inform the contribution of both individual loci and their combinatorial effects to phenotypes relevant to the etiopathology of DS.

TM6SF2 rs58542926 impacts lipid processing in liver and small intestine.

The transmembrane 6 superfamily member 2 (TM6SF2) loss-of-function variant rs58542926 is a genetic risk factor for nonalcoholic fatty liver disease and progression to fibrosis but is paradoxically associated with lower levels of hepatically derived triglyceride-rich lipoproteins. TM6SF2 is expressed predominantly in liver and small intestine, sites for triglyceride-rich lipoprotein biogenesis and export. In light of this, we hypothesized that TM6SF2 may exhibit analogous effects on both liver and intestine lipid homeostasis. To test this, we genotyped rs58542926 in 983 bariatric surgery patients from the Geisinger Medical Center for Nutrition and Weight Management, Geisinger Health System, in Pennsylvania and from 3,556 study participants enrolled in the Amish Complex Disease Research Program. Although these two cohorts have different metabolic profiles, carriers in both cohorts had improved fasting lipid profiles. Importantly, following a high-fat challenge, carriers in the Amish Complex Disease Research Program cohort exhibited significantly lower postprandial serum triglycerides, suggestive of a role for TM6SF2 in the small intestine. To gain further insight into this putative role, effects of TM6SF2 deficiency were studied in a zebrafish model and in cultured human Caco-2 enterocytes. In both systems TM6SF2 deficiency resulted in defects in small intestine metabolism in response to dietary lipids, including significantly increased lipid accumulation, decreased lipid clearance, and increased endoplasmic reticulum stress. CONCLUSIONS: These data strongly support a role of TM6SF2 in the regulation of postprandial lipemia, potentially through a similar function for TM6SF2 in the lipidation and/or export of both hepatically and intestinally derived triglyceride-rich lipoproteins. (Hepatology 2017;65:1526-1542).

FOXE3 contributes to Peters anomaly through transcriptional regulation of an autophagy-associated protein termed DNAJB1.

FOXE3 is a lens-specific transcription factor that has been associated with anterior segment ocular dysgenesis. To determine the transcriptional target(s) of FOXE3 that are indispensable for the anterior segment development, we examined the transcriptome and the proteome of cells expressing truncated FOXE3 responsible for Peters anomaly identified through linkage-coupled next-generation whole-exome sequencing. We found that DNAJB1, an autophagy-associated protein, was the only candidate exhibiting differential expression in both screens. We confirmed the candidacy of DNAJB1 through chromatin immunoprecipitation and luciferase assays while knockdown of DNAJB1 in human lens epithelial cells resulted in a mitotic arrest. Subsequently, we targeted dnajb1a in zebrafish through injection of a splice-blocking morpholino. The dnajb1a morphants exhibited underdeveloped cataractous lenses with persistent apoptotic nuclei. In conclusion, here we report DNAJB1 is a transcriptional target of FOXE3 in a novel pathway that is crucial for the development of the anterior segment of the eye.

Differential effects on ß-cell mass by disruption of Bardet-Biedl syndrome or Alstrom syndrome genes.

Rare genetic syndromes characterized by early-onset type 2 diabetes have revealed the importance of pancreatic ß-cells in genetic susceptibility to diabetes. However, the role of genetic regulation of ß-cells in disorders that are also characterized by highly penetrant obesity, a major additional risk factor, is unclear. In this study, we investigated the contribution of genes associated with two obesity ciliopathies, Bardet-Biedl Syndrome and Alstrom Syndrome, to the production and maintenance of pancreatic ß-cells. Using zebrafish models of these syndromes, we identified opposing effects on production of ß-cells. Loss of the Alstrom gene, alms1, resulted in a significant decrease in ß-cell production whereas loss of BBS genes, bbs1 or bbs4, resulted in a significant increase. Examination of the regulatory program underlying ß-cell production suggested that these effects were specific to ß-cells. In addition to the initial production of ß-cells, we observed significant differences in their continued maintenance. Under prolonged exposure to high glucose conditions, alms1-deficient ß-cells were unable to continually expand as a result of decreased proliferation and increased cell death. Although bbs1-deficient ß-cells were similarly susceptible to apoptosis, the overall maintenance of ß-cell number in those animals was sustained likely due to increased proliferation. Taken together, these findings implicate discrepant production and maintenance of ß-cells in the differential susceptibility to diabetes found between these two genetic syndromes.

RFX transcription factors are essential for hearing in mice.

Sensorineural hearing loss is a common and currently irreversible disorder, because mammalian hair cells (HCs) do not regenerate and current stem cell and gene delivery protocols result only in immature HC-like cells. Importantly, although the transcriptional regulators of embryonic HC development have been described, little is known about the postnatal regulators of maturating HCs. Here we apply a cell type-specific functional genomic analysis to the transcriptomes of auditory and vestibular sensory epithelia from early postnatal mice. We identify RFX transcription factors as essential and evolutionarily conserved regulators of the HC-specific transcriptomes, and detect Rfx1,2,3,5 and 7 in the developing HCs. To understand the role of RFX in hearing, we generate Rfx1/3 conditional knockout mice. We show that these mice are deaf secondary to rapid loss of initially well-formed outer HCs. These data identify an essential role for RFX in hearing and survival of the terminally differentiating outer HCs.

BBS4 is necessary for ciliary localization of TrkB receptor and activation by BDNF.

Primary cilia regulate an expanding list of signaling pathways in many different cell types. It is likely that identification of the full catalog of pathways associated with cilia will be necessary to fully understand their role in regulation of signaling and the implications for diseases associated with their dysfunction, ciliopathies. Bardet-Biedl Syndrome (BBS) is one such ciliopathy which is characterized by a spectrum of phenotypes. These include neural defects such as impaired cognitive development, centrally mediated hyperphagia and peripheral sensory defects. Here we investigate potential defects in a signaling pathway associated with neuronal function, brain derived neurotrophic factor (BDNF) signaling. Upon loss of BBS4 expression in cultured cells, we observed decreased phosphorylation and activation by BDNF of its target receptor, TrkB. Assessment of ciliary localization revealed that, TrkB localized to the axonemes or basal bodies of cilia only in the presence of BDNF. Axonemal localization, specifically, was abrogated with loss of BBS4. Finally, we present evidence that loss of the ciliary axoneme through depletion of KIF3A impedes activation of TrkB. Taken together, these data suggest the possibility of a previously uninvestigated pathway associated with perturbation of ciliary proteins.

Basal body proteins regulate Notch signaling through endosomal trafficking.

Proteins associated with primary cilia and basal bodies mediate numerous signaling pathways, but little is known about their role in Notch signaling. Here, we report that loss of the Bardet-Biedl syndrome proteins BBS1 or BBS4 produces increased Notch-directed transcription in a zebrafish reporter line and in human cell lines. Pathway overactivation is accompanied by reduced localization of Notch receptor at both the plasma membrane and the cilium. In Drosophila mutants, overactivation of Notch can result from receptor accumulation in endosomes, and recent studies implicate ciliary proteins in endosomal trafficking, suggesting a possible mechanism by which overactivation occurs in BBS mutants. Consistent with this, we observe genetic interaction of BBS1 and BBS4 with the endosomal sorting complexes required for transport (ESCRT) gene TSG101 and accumulation of receptor in late endosomes, reduced endosomal recycling and reduced receptor degradation in lysosomes. We observe similar defects with disruption of BBS3. Loss of another basal body protein, ALMS1, also enhances Notch activation and the accumulation of receptor in late endosomes, but does not disrupt recycling. These findings suggest a role for these proteins in the regulation of Notch through endosomal trafficking of the receptor.

Ciliopathy proteins regulate paracrine signaling by modulating proteasomal degradation of mediators.

Cilia are critical mediators of paracrine signaling; however, it is unknown whether proteins that contribute to ciliopathies converge on multiple paracrine pathways through a common mechanism. Here, we show that loss of cilopathy-associated proteins Bardet-Biedl syndrome 4 (BBS4) or oral-facial-digital syndrome 1 (OFD1) results in the accumulation of signaling mediators normally targeted for proteasomal degradation. In WT cells, several BBS proteins and OFD1 interacted with proteasomal subunits, and loss of either BBS4 or OFD1 led to depletion of multiple subunits from the centrosomal proteasome. Furthermore, overexpression of proteasomal regulatory components or treatment with proteasomal activators sulforaphane (SFN) and mevalonolactone (MVA) ameliorated signaling defects in cells lacking BBS1, BBS4, and OFD1, in morphant zebrafish embryos, and in induced neurons from Ofd1-deficient mice. Finally, we tested the hypothesis that other proteasome-dependent pathways not known to be associated with ciliopathies are defective in the absence of ciliopathy proteins. We found that loss of BBS1, BBS4, or OFD1 led to decreased NF-κB activity and concomitant IκBß accumulation and that these defects were ameliorated with SFN treatment. Taken together, our data indicate that basal body proteasomal regulation governs paracrine signaling pathways and suggest that augmenting proteasomal function might benefit ciliopathy patients.

Genetic deletion of trkB.T1 increases neuromuscular function.

Neurotrophin-dependent activation of the tyrosine kinase receptor trkB.FL modulates neuromuscular synapse maintenance and function; however, it is unclear what role the alternative splice variant, truncated trkB (trkB.T1), may have in the peripheral neuromuscular axis. We examined this question in trkB.T1 null mice and demonstrate that in vivo neuromuscular performance and nerve-evoked muscle tension are significantly increased. In vitro assays indicated that the gain-in-function in trkB.T1(-/-) animals resulted specifically from an increased muscle contractility, and increased electrically evoked calcium release. In the trkB.T1 null muscle, we identified an increase in Akt activation in resting muscle as well as a significant increase in trkB.FL and Akt activation in response to contractile activity. On the basis of these findings, we conclude that the trkB signaling pathway might represent a novel target for intervention across diseases characterized by deficits in neuromuscular function.

Brain-derived neurotrophic factor modulates antiretroviral-induced mechanical allodynia in the mouse.

Nucleoside reverse transcriptase inhibitors (NRTIs) are key components of HIV/AIDS treatment to reduce viral load. However, these drugs can induce chronic neuropathic pain, leading to increased morbidity in HIV patients. This study examines the role of brain-derived neurotrophic factor (BDNF) in the spinal dorsal horn (SDH) in development of mechanical allodynia in male C57BL/6J mice treated with the NRTI stavudine (d4T). After d4T administration, mice developed increased neuronal activity and BDNF expression in the SDH and hind paw mechanical allodynia that was exacerbated by intrathecal BDNF administration. Intrathecal BDNF alone also increased neuronal activity and caused mechanical allodynia. Because excess BDNF amplified d4T-induced mechanical allodynia and neuronal activity, the impact of decreasing BDNF in the SDH was investigated. After d4T, BDNF heterozygous mice were less allodynic than wild-type littermates, which was negated by intrathecal BDNF administration. Finally, pretreatment with intrathecal trkB-Fc chimera prior to d4T or administration of the tyrosine kinase inhibitor K252a 3 days after d4T blocked BDNF-mediated signaling, significantly attenuated the development of mechanical allodynia (trkB-Fc), and decreased neuronal activity (trkB-Fc and K252a). Taken together, these findings provide evidence that BDNF in the SDH contributes to the development of NRTI-induced painful peripheral neuropathy and may represent a new therapeutic opportunity.

Functional analyses of variants reveal a significant role for dominant negative and common alleles in oligogenic Bardet-Biedl syndrome.

Technological advances hold the promise of rapidly catalyzing the discovery of pathogenic variants for genetic disease. However, this possibility is tempered by limitations in interpreting the functional consequences of genetic variation at candidate loci. Here, we present a systematic approach, grounded on physiologically relevant assays, to evaluate the mutational content (125 alleles) of the 14 genes associated with Bardet-Biedl syndrome (BBS). A combination of in vivo assays with subsequent in vitro validation suggests that a significant fraction of BBS-associated mutations have a dominant-negative mode of action. Moreover, we find that a subset of common alleles, previously considered to be benign, are, in fact, detrimental to protein function and can interact with strong rare alleles to modulate disease presentation. These data represent a comprehensive evaluation of genetic load in a multilocus disease. Importantly, superimposition of these results to human genetics data suggests a previously underappreciated complexity in disease architecture that might be shared among diverse clinical phenotypes.

In vivo evidence that truncated trkB.T1 participates in nociception.

Brain-Derived Neurotrophic Factor (BDNF) is a central nervous system modulator of nociception. In animal models of chronic pain, BDNF exerts its effects on nociceptive processing by binding to the full-length receptor tropomyosin-related kinase B (trkB.FL) and transducing intracellular signaling to produce nocifensive behaviors. In addition to trkB.FL, the trkB locus also produces a widely-expressed alternatively-spliced truncated isoform, trkB.T1. TrkB.T1 binds BDNF with high affinity; however the unique 11 amino acid intracellular cytoplasmic tail lacks the kinase domain of trkB.FL. Recently, trkB.T1 was shown to be specifically up-regulated in a model of HIV-associated neuropathic pain, potentially implicating trkB.T1 as a modulator of nociception. Here, we report that trkB.T1 mRNA and protein is up-regulated in the spinal dorsal horn at times following antiretroviral drug treatment and hind paw inflammation in which nocifensive behaviors develop. While genetic depletion of trkB.T1 did not affect baseline mechanical and thermal thresholds, the absence of trkB.T1 resulted in significant attenuation of inflammation- and antiretroviral-induced nocifensive behaviors. Our results suggest that trkB.T1 up-regulation following antiretroviral treatment and tissue inflammation participates in the development and maintenance of nocifensive behavior and may represent a novel therapeutic target for pain treatment.

Genome-wide screen identifies drug-induced regulation of the gene giant axonal neuropathy (Gan) in a mouse model of antiretroviral-induced painful peripheral neuropathy.

Painful peripheral neuropathy is a debilitating complication of the treatment of HIV with nucleoside reverse transcriptase inhibitors (NRTIs). Patients are living longer with these drugs; however many develop excruciating, unremitting, and often treatment-limiting neuropathy that is resistant to conventional pain management therapies. Improving patient comfort and quality of life is paramount and depends on a clearer understanding of this devastating side effect. The mechanisms underlying the development of NRTI-induced neuropathy, however, remain unclear. Using a mouse model of NRTI-induced neuropathy, the authors conducted an unbiased whole-genome microarray screen to identify molecular targets in the spinal dorsal horn, which is the location where integration of ascending sensory transmission and descending modulatory effects occur. Analysis of the microarray data identified a change in the gene giant axonal neuropathy 1 (Gan1). Mutation of this gene has been linked to the development of giant axonal neuropathy (GAN), a rare autosomal recessive condition characterized by a progressive sensorimotor neuropathy. Gan1 has not been previously linked to nerve pathologies in other populations. In this study, downregulation of the Gan1 gene and the gene protein product, gigaxonin, was validated via quantitative polymerase chain reaction ([qPCR] gene expression) and Western blot analyses (protein level). Our report is the first to suggest that Gan1 might be a novel molecular target in the development of NRTI-induced peripheral neuropathy with implications for new therapeutic approaches to preventing or reducing a significant side effect of HIV treatment.

An essential role for DYF-11/MIP-T3 in assembling functional intraflagellar transport complexes.

MIP-T3 is a human protein found previously to associate with microtubules and the kinesin-interacting neuronal protein DISC1 (Disrupted-in-Schizophrenia 1), but whose cellular function(s) remains unknown. Here we demonstrate that the C. elegans MIP-T3 ortholog DYF-11 is an intraflagellar transport (IFT) protein that plays a critical role in assembling functional kinesin motor-IFT particle complexes. We have cloned a loss of function dyf-11 mutant in which several key components of the IFT machinery, including Kinesin-II, as well as IFT subcomplex A and B proteins, fail to enter ciliary axonemes and/or mislocalize, resulting in compromised ciliary structures and sensory functions, and abnormal lipid accumulation. Analyses in different mutant backgrounds further suggest that DYF-11 functions as a novel component of IFT subcomplex B. Consistent with an evolutionarily conserved cilia-associated role, mammalian MIP-T3 localizes to basal bodies and cilia, and zebrafish mipt3 functions synergistically with the Bardet-Biedl syndrome protein Bbs4 to ensure proper gastrulation, a key cilium- and basal body-dependent developmental process. Our findings therefore implicate MIP-T3 in a previously unknown but critical role in cilium biogenesis and further highlight the emerging role of this organelle in vertebrate development.

Hypomorphic mutations in syndromic encephalocele genes are associated with Bardet-Biedl syndrome.

Meckel-Gruber syndrome (MKS) is a genetically heterogeneous, neonatally lethal malformation and the most common form of syndromic neural tube defect (NTD). To date, several MKS-associated genes have been identified whose protein products affect ciliary function. Here we show that mutations in MKS1, MKS3 and CEP290 (also known as NPHP6) either can cause Bardet-Biedl syndrome (BBS) or may have a potential epistatic effect on mutations in known BBS-associated loci. Five of six families with both MKS1 and BBS mutations manifested seizures, a feature that is not a typical component of either syndrome. Functional studies in zebrafish showed that mks1 is necessary for gastrulation movements and that it interacts genetically with known bbs genes. Similarly, we found two families with missense or splice mutations in MKS3, in one of which the affected individual also bears a homozygous nonsense mutation in CEP290 that is likely to truncate the C terminus of the protein. These data extend the genetic stratification of ciliopathies and suggest that BBS and MKS, although distinct clinically, are allelic forms of the same molecular spectrum.

Disruption of the basal body compromises proteasomal function and perturbs intracellular Wnt response.

Primary cilia and basal bodies are evolutionarily conserved organelles that mediate communication between the intracellular and extracellular environments. Here we show that bbs1, bbs4 and mkks (also known as bbs6), which encode basal body proteins, are required for convergence and extension in zebrafish and interact with wnt11 and wnt5b. Suppression of bbs1, bbs4 and mkks transcripts results in stabilization of beta-catenin with concomitant upregulation of T-cell factor (TCF)-dependent transcription in both zebrafish embryos and mammalian ciliated cells, a defect phenocopied by the silencing of the axonemal kinesin subunit KIF3A but not by chemical disruption of the cytoplasmic microtubule network. These observations are attributable partly to defective degradation by the proteasome; suppression of BBS4 leads to perturbed proteasomal targeting and concomitant accumulation of cytoplasmic beta-catenin. Cumulatively, our data indicate that the basal body is an important regulator of Wnt signal interpretation through selective proteolysis and suggest that defects in this system may contribute to phenotypes pathognomonic of human ciliopathies.