Novel SIX6 mutations cause recessively inherited congenital cataract, microcornea, and corneal opacification with or without coloboma and microphthalmia

Purpose

To investigate the molecular basis of recessively inherited congenital cataract, microcornea, and corneal opacification with or without coloboma and microphthalmia in two consanguineous families.

Methods

Conventional autozygosity mapping was performed using single nucleotide polymorphism (SNP) microarrays. Whole-exome sequencing was completed on genomic DNA from one affected member of each family. Exome sequence data were also used for homozygosity mapping and copy number variation analysis. PCR and Sanger sequencing were used to confirm the identification of mutations and to screen further patients. Evolutionary conservation of protein sequences was assessed using CLUSTALW, and protein structures were modeled using PyMol.

Results

In family MEP68, a novel homozygous nucleotide substitution in SIX6 was found, c.547G>C, that converts the evolutionarily conserved aspartic acid residue at the 183^rd amino acid in the protein to a histidine, p.(Asp183His). This residue mapped to the third helix of the DNA-binding homeobox domain in SIX6, which interacts with the major groove of double-stranded DNA. This interaction is likely to be disrupted by the mutation. In family F1332, a novel homozygous 1034 bp deletion that encompasses the first exon of SIX6 was identified, chr14:g.60975890_60976923del. Both mutations segregated with the disease phenotype as expected for a recessive condition and were absent from publicly available variant databases.

Conclusions

Our findings expand the mutation spectrum in this form of inherited eye disease and confirm that homozygous human SIX6 mutations cause a developmental spectrum of ocular phenotypes that includes not only the previously described features of microphthalmia, coloboma, and congenital cataract but also corneal abnormalities.

Biallelic variants in DNA2 cause microcephalic primordial dwarfism

Microcephalic primordial dwarfism (MPD) is a group of rare single-gene disorders characterized by the extreme reduction in brain and body size from early development onwards. Proteins encoded by MPD-associated genes play important roles in fundamental cellular processes, notably genome replication and repair. Here we report the identification of four MPD individuals with biallelic variants in DNA2, which encodes an adenosine triphosphate (ATP)-dependent helicase/nuclease involved in DNA replication and repair. We demonstrate that the two intronic variants (c.1764-38_1764-37ins(53) and c.74+4A>C) found in these individuals substantially impair DNA2 transcript splicing. Additionally, we identify a missense variant (c.1963A>G), affecting a residue of the ATP-dependent helicase domain that is highly conserved between humans and yeast, with the resulting substitution (p.Thr655Ala) predicted to directly impact ATP/ADP (adenosine diphosphate) binding by DNA2. Our findings support the pathogenicity of these variants as biallelic hypomorphic mutations, establishing DNA2 as an MPD disease gene.

Gain-of-function DNMT3A mutations cause microcephalic dwarfism and hypermethylation of Polycomb-regulated regions

DNA methylation and Polycomb are key factors in the establishment of vertebrate cellular identity and fate. Here we report de novo missense mutations in DNMT3A, which encodes the DNA methyltransferase DNMT3A. These mutations cause microcephalic dwarfism, a hypocellular disorder of extreme global growth failure. Substitutions in the PWWP domain abrogate binding to the histone modifications H3K36me2 and H3K36me3, and alter DNA methylation in patient cells. Polycomb-associated DNA methylation valleys, hypomethylated domains encompassing developmental genes, become methylated with concomitant depletion of H3K27me3 and H3K4me3 bivalent marks. Such de novo DNA methylation occurs during differentiation of Dnmt3aW326R pluripotent cells in vitro, and is also evident in Dnmt3aW326R/+ dwarf mice. We therefore propose that the interaction of the DNMT3A PWWP domain with H3K36me2 and H3K36me3 normally limits DNA methylation of Polycomb-marked regions. Our findings implicate the interplay between DNA methylation and Polycomb at key developmental regulators as a determinant of organism size in mammals.

DNA Polymerase Epsilon Deficiency Causes IMAGe Syndrome with Variable Immunodeficiency

During genome replication, polymerase epsilon (Pol ε) acts as the major leading-strand DNA polymerase. Here we report the identification of biallelic mutations in POLE, encoding the Pol ε catalytic subunit POLE1, in 15 individuals from 12 families. Phenotypically, these individuals had clinical features closely resembling IMAGe syndrome (intrauterine growth restriction [IUGR], metaphyseal dysplasia, adrenal hypoplasia congenita, and genitourinary anomalies in males), a disorder previously associated with gain-of-function mutations in CDKN1C. POLE1-deficient individuals also exhibited distinctive facial features and variable immune dysfunction with evidence of lymphocyte deficiency. All subjects shared the same intronic variant (c.1686+32C>G) as part of a common haplotype, in combination with different loss-of-function variants in trans. The intronic variant alters splicing, and together the biallelic mutations lead to cellular deficiency of Pol ε and delayed S-phase progression. In summary, we establish POLE as a second gene in which mutations cause IMAGe syndrome. These findings add to a growing list of disorders due to mutations in DNA replication genes that manifest growth restriction alongside adrenal dysfunction and/or immunodeficiency, consolidating these as replisome phenotypes and highlighting a need for future studies to understand the tissue-specific development roles of the encoded proteins.

A Recurrent De Novo Heterozygous COG4 Substitution Leads to Saul-Wilson Syndrome, Disrupted Vesicular Trafficking, and Altered Proteoglycan Glycosylation

The conserved oligomeric Golgi (COG) complex is involved in intracellular vesicular transport, and is composed of eight subunits distributed in two lobes, lobe A (COG1-4) and lobe B (COG5-8). We describe fourteen individuals with Saul-Wilson syndrome, a rare form of primordial dwarfism with characteristic facial and radiographic features. All affected subjects harbored heterozygous de novo variants in COG4, giving rise to the same recurrent amino acid substitution (p.Gly516Arg). Affected individuals' fibroblasts, whose COG4 mRNA and protein were not decreased, exhibited delayed anterograde vesicular trafficking from the ER to the Golgi and accelerated retrograde vesicular recycling from the Golgi to the ER. This altered steady-state equilibrium led to a decrease in Golgi volume, as well as morphologic abnormalities with collapse of the Golgi stacks. Despite these abnormalities of the Golgi apparatus, protein glycosylation in sera and fibroblasts from affected subjects was not notably altered, but decorin, a proteoglycan secreted into the extracellular matrix, showed altered Golgi-dependent glycosylation. In summary, we define a specific heterozygous COG4 substitution as the molecular basis of Saul-Wilson syndrome, a rare skeletal dysplasia distinct from biallelic COG4-CDG.

Mutations in TOP3A Cause a Bloom Syndrome-like Disorder

Bloom syndrome, caused by biallelic mutations in BLM, is characterized by prenatal-onset growth deficiency, short stature, an erythematous photosensitive malar rash, and increased cancer predisposition. Diagnostically, a hallmark feature is the presence of increased sister chromatid exchanges (SCEs) on cytogenetic testing. Here, we describe biallelic mutations in TOP3A in ten individuals with prenatal-onset growth restriction and microcephaly. TOP3A encodes topoisomerase III alpha (TopIIIα), which binds to BLM as part of the BTRR complex, and promotes dissolution of double Holliday junctions arising during homologous recombination. We also identify a homozygous truncating variant in RMI1, which encodes another component of the BTRR complex, in two individuals with microcephalic dwarfism. The TOP3A mutations substantially reduce cellular levels of TopIIIα, and consequently subjects' cells demonstrate elevated rates of SCE. Unresolved DNA recombination and/or replication intermediates persist into mitosis, leading to chromosome segregation defects and genome instability that most likely explain the growth restriction seen in these subjects and in Bloom syndrome. Clinical features of mitochondrial dysfunction are evident in several individuals with biallelic TOP3A mutations, consistent with the recently reported additional function of TopIIIα in mitochondrial DNA decatenation. In summary, our findings establish TOP3A mutations as an additional cause of prenatal-onset short stature with increased cytogenetic SCEs and implicate the decatenation activity of the BTRR complex in their pathogenesis.

The expanding phenotype of RNU4ATAC pathogenic variants to Lowry Wood syndrome

RNU4ATAC pathogenic variants to date have been associated with microcephalic osteodysplastic primordial dwarfism, type 1 and Roifman syndrome. Both conditions are clinically distinct skeletal dysplasias with microcephalic osteodysplastic primordial dwarfism, type 1 having a more severe phenotype than Roifman syndrome. Some of the overlapping features of the two conditions include developmental delay, microcephaly, and immune deficiency. The features also overlap with Lowry Wood syndrome, another rare but well-defined skeletal dysplasia for which the genetic etiology has not been identified. Characteristic features include multiple epiphyseal dysplasia and microcephaly. Here, we describe three patients with Lowry Wood syndrome with biallelic RNU4ATAC pathogenic variants. This report expands the phenotypic spectrum for biallelic RNU4ATAC disorder causing variants and is the first to establish the genetic cause for Lowry Wood syndrome.

Characterizing the morbid genome of ciliopathies

Background

Ciliopathies are clinically diverse disorders of the primary cilium. Remarkable progress has been made in understanding the molecular basis of these genetically heterogeneous conditions; however, our knowledge of their morbid genome, pleiotropy, and variable expressivity remains incomplete.

Results

We applied genomic approaches on a large patient cohort of 371 affected individuals from 265 families, with phenotypes that span the entire ciliopathy spectrum. Likely causal mutations in previously described ciliopathy genes were identified in 85% (225/265) of the families, adding 32 novel alleles. Consistent with a fully penetrant model for these genes, we found no significant difference in their “mutation load” beyond the causal variants between our ciliopathy cohort and a control non-ciliopathy cohort. Genomic analysis of our cohort further identified mutations in a novel morbid gene TXNDC15, encoding a thiol isomerase, based on independent loss of function mutations in individuals with a consistent ciliopathy phenotype (Meckel-Gruber syndrome) and a functional effect of its deficiency on ciliary signaling. Our study also highlighted seven novel candidate genes (TRAPPC3, EXOC3L2, FAM98C, C17orf61, LRRCC1, NEK4, and CELSR2) some of which have established links to ciliogenesis. Finally, we show that the morbid genome of ciliopathies encompasses many founder mutations, the combined carrier frequency of which accounts for a high disease burden in the study population.

Conclusions

Our study increases our understanding of the morbid genome of ciliopathies. We also provide the strongest evidence, to date, in support of the classical Mendelian inheritance of Bardet-Biedl syndrome and other ciliopathies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1099-5) contains supplementary material, which is available to authorized users.

Mutations in the pH-Sensing G-protein-Coupled Receptor GPR68 Cause Amelogenesis Imperfecta

Amelogenesis is the process of dental enamel formation, leading to the deposition of the hardest tissue in the human body. This process requires the intricate regulation of ion transport and controlled changes to the pH of the developing enamel matrix. The means by which the enamel organ regulates pH during amelogenesis is largely unknown. We identified rare homozygous variants in GPR68 in three families with amelogenesis imperfecta, a genetically and phenotypically heterogeneous group of inherited conditions associated with abnormal enamel formation. Each of these homozygous variants (a large in-frame deletion, a frameshift deletion, and a missense variant) were predicted to result in loss of function. GPR68 encodes a proton-sensing G-protein-coupled receptor with sensitivity in the pH range that occurs in the developing enamel matrix during amelogenesis. Immunohistochemistry of rat mandibles confirmed localization of GPR68 in the enamel organ at all stages of amelogenesis. Our data identify a role for GPR68 as a proton sensor that is required for proper enamel formation.

Pathogenic Variants in PIGG Cause Intellectual Disability with Seizures and Hypotonia

Glycosylphosphatidylinositol (GPI) is a glycolipid that anchors >150 various proteins to the cell surface. At least 27 genes are involved in biosynthesis and transport of GPI-anchored proteins (GPI-APs). To date, mutations in 13 of these genes are known to cause inherited GPI deficiencies (IGDs), and all are inherited as recessive traits. IGDs mainly manifest as intellectual disability, epilepsy, coarse facial features, and multiple organ anomalies. These symptoms are caused by the decreased surface expression of GPI-APs or by structural abnormalities of GPI. Here, we present five affected individuals (from two consanguineous families from Egypt and Pakistan and one non-consanguineous family from Japan) who show intellectual disability, hypotonia, and early-onset seizures. We identified pathogenic variants in PIGG, a gene in the GPI pathway. In the consanguineous families, homozygous variants c.928C>T (p.Gln310(∗)) and c.2261+1G>C were found, whereas the Japanese individual was compound heterozygous for c.2005C>T (p.Arg669Cys) and a 2.4 Mb deletion involving PIGG. PIGG is the enzyme that modifies the second mannose with ethanolamine phosphate, which is removed soon after GPI is attached to the protein. Physiological significance of this transient modification has been unclear. Using B lymphoblasts from affected individuals of the Egyptian and Japanese families, we revealed that PIGG activity was almost completely abolished; however, the GPI-APs had normal surface levels and normal structure, indicating that the pathogenesis of PIGG deficiency is not yet fully understood. The discovery of pathogenic variants in PIGG expands the spectrum of IGDs and further enhances our understanding of this etiopathogenic class of intellectual disability.

Congenital Myasthenic Syndrome Type 19 Is Caused by Mutations in COL13A1, Encoding the Atypical Non-fibrillar Collagen Type XIII α1 Chain

The neuromuscular junction (NMJ) consists of a tripartite synapse with a presynaptic nerve terminal, Schwann cells that ensheathe the terminal bouton, and a highly specialized postsynaptic membrane. Synaptic structural integrity is crucial for efficient signal transmission. Congenital myasthenic syndromes (CMSs) are a heterogeneous group of inherited disorders that result from impaired neuromuscular transmission, caused by mutations in genes encoding proteins that are involved in synaptic transmission and in forming and maintaining the structural integrity of NMJs. To identify further causes of CMSs, we performed whole-exome sequencing (WES) in families without an identified mutation in known CMS-associated genes. In two families affected by a previously undefined CMS, we identified homozygous loss-of-function mutations in COL13A1, which encodes the alpha chain of an atypical non-fibrillar collagen with a single transmembrane domain. COL13A1 localized to the human muscle motor endplate. Using CRISPR-Cas9 genome editing, modeling of the COL13A1 c.1171delG (p.Leu392Sfs^∗71) frameshift mutation in the C2C12 cell line reduced acetylcholine receptor (AChR) clustering during myotube differentiation. This highlights the crucial role of collagen XIII in the formation and maintenance of the NMJ. Our results therefore delineate a myasthenic disorder that is caused by loss-of-function mutations in COL13A1, encoding a protein involved in organization of the NMJ, and emphasize the importance of appropriate symptomatic treatment for these individuals.

HACE1 deficiency causes an autosomal recessive neurodevelopmental syndrome

Background

The genetic aetiology of neurodevelopmental defects is extremely diverse, and the lack of distinctive phenotypic features means that genetic criteria are often required for accurate diagnostic classification. We aimed to identify the causative genetic lesions in two families in which eight affected individuals displayed variable learning disability, spasticity and abnormal gait.

Methods

Autosomal recessive inheritance was suggested by consanguinity in one family and by sibling recurrences with normal parents in the second. Autozygosity mapping and exome sequencing, respectively, were used to identify the causative gene.

Results

In both families, biallelic loss-of-function mutations in HACE1 were identified. HACE1 is an E3 ubiquitin ligase that regulates the activity of cellular GTPases, including Rac1 and members of the Rab family. In the consanguineous family, a homozygous mutation p.R219* predicted a truncated protein entirely lacking its catalytic domain. In the other family, compound heterozygosity for nonsense mutation p.R748* and a 20-nt insertion interrupting the catalytic homologous to the E6-AP carboxyl terminus (HECT) domain was present; western blot analysis of patient cells revealed an absence of detectable HACE1 protein.

Conclusion

HACE1 mutations underlie a new autosomal recessive neurodevelopmental disorder. Previous studies have implicated HACE1 as a tumour suppressor gene; however, since cancer predisposition was not observed either in homozygous or heterozygous mutation carriers, this concept may require re-evaluation.

Variability of systemic and oro-dental phenotype in two families with non-lethal Raine syndrome with FAM20C mutations

Background

Raine syndrome (RS) is a rare autosomal recessive bone dysplasia typified by osteosclerosis and dysmorphic facies due to FAM20C mutations. Initially reported as lethal in infancy, survival is possible into adulthood. We describe the molecular analysis and clinical phenotypes of five individuals from two consanguineous Brazilian families with attenuated Raine Syndrome with previously unreported features.

Methods

The medical and dental clinical records were reviewed. Extracted deciduous and permanent teeth as well as oral soft tissues were analysed. Whole exome sequencing was undertaken and FAM20C cDNA sequenced in family 1.

Results

Family 1 included 3 siblings with hypoplastic Amelogenesis Imperfecta (AI) (inherited abnormal dental enamel formation). Mild facial dysmorphism was noted in the absence of other obvious skeletal or growth abnormalities. A mild hypophosphataemia and soft tissue ectopic mineralization were present. A homozygous FAM20C donor splice site mutation (c.784 + 5 g > c) was identified which led to abnormal cDNA sequence. Family 2 included 2 siblings with hypoplastic AI and tooth dentine abnormalities as part of a more obvious syndrome with facial dysmorphism. There was hypophosphataemia, soft tissue ectopic mineralization, but no osteosclerosis. A homozygous missense mutation in FAM20C (c.1487C > T; p.P496L) was identified.

Conclusions

The clinical phenotype of non-lethal Raine Syndrome is more variable, including between affected siblings, than previously described and an adverse impact on bone growth and health may not be a prominent feature. By contrast, a profound failure of dental enamel formation leading to a distinctive hypoplastic AI in all teeth should alert clinicians to the possibility of FAM20C mutations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-015-0154-5) contains supplementary material, which is available to authorized users.

Mutation screening of retinal dystrophy patients by targeted capture from tagged pooled DNAs and next generation sequencing

Purpose

Retinal dystrophies are genetically heterogeneous, resulting from mutations in over 200 genes. Prior to the development of massively parallel sequencing, comprehensive genetic screening was unobtainable for most patients. Identifying the causative genetic mutation facilitates genetic counselling, carrier testing and prenatal/pre-implantation diagnosis, and often leads to a clearer prognosis. In addition, in a proportion of cases, when the mutation is known treatment can be optimised and patients are eligible for enrolment into clinical trials for gene-specific therapies.

Methods

Patient genomic DNA was sheared, tagged and pooled in batches of four samples, prior to targeted capture and next generation sequencing. The enrichment reagent was designed against genes listed on the RetNet database (July 2010). Sequence data were aligned to the human genome and variants were filtered to identify potential pathogenic mutations. These were confirmed by Sanger sequencing.

Results

Molecular analysis of 20 DNAs from retinal dystrophy patients identified likely pathogenic mutations in 12 cases, many of them known and/or confirmed by segregation. These included previously described mutations in ABCA4 (c.6088C>T,p.R2030*; c.5882G>A,p.G1961E), BBS2 (c.1895G>C,p.R632P), GUCY2D (c.2512C>T,p.R838C), PROM1 (c.1117C>T,p.R373C), RDH12 (c.601T>C,p.C201R; c.506G>A,p.R169Q), RPGRIP1 (c.3565C>T,p.R1189*) and SPATA7 (c.253C>T,p.R85*) and new mutations in ABCA4 (c.3328+1G>C), CRB1 (c.2832_2842+23del), RP2 (c.884-1G>T) and USH2A (c.12874A>G,p.N4292D).

Conclusions

Tagging and pooling DNA prior to targeted capture of known retinal dystrophy genes identified mutations in 60% of cases. This relatively high success rate may reflect enrichment for consanguineous cases in the local Yorkshire population, and the use of multiplex families. Nevertheless this is a promising high throughput approach to retinal dystrophy diagnostics.

SAMS, a syndrome of short stature, auditory-canal atresia, mandibular hypoplasia, and skeletal abnormalities is a unique neurocristopathy caused by mutations in Goosecoid

Short stature, auditory canal atresia, mandibular hypoplasia, and skeletal abnormalities (SAMS) has been reported previously to be a rare, autosomal-recessive developmental disorder with other, unique rhizomelic skeletal anomalies. These include bilateral humeral hypoplasia, humeroscapular synostosis, pelvic abnormalities, and proximal defects of the femora. To identify the genetic basis of SAMS, we used molecular karyotyping and whole-exome sequencing (WES) to study small, unrelated families. Filtering of variants from the WES data included segregation analysis followed by comparison of in-house exomes. We identified a homozygous 306 kb microdeletion and homozygous predicted null mutations of GSC, encoding Goosecoid homeobox protein, a paired-like homeodomain transcription factor. This confirms that SAMS is a human malformation syndrome resulting from GSC mutations. Previously, Goosecoid has been shown to be a determinant at the Xenopus gastrula organizer region and a segment-polarity determinant in Drosophila. In the present report, we present data on Goosecoid protein localization in staged mouse embryos. These data and the SAMS clinical phenotype both suggest that Goosecoid is a downstream effector of the regulatory networks that define neural-crest cell-fate specification and subsequent mesoderm cell lineages in mammals, particularly during shoulder and hip formation. Our findings confirm that Goosecoid has an essential role in human craniofacial and joint development and suggest that Goosecoid is an essential regulator of mesodermal patterning in mammals and that it has specific functions in neural crest cell derivatives.

Simple and efficient identification of rare recessive pathologically important sequence variants from next generation exome sequence data

Massively parallel ("next generation") DNA sequencing (NGS) has quickly become the method of choice for seeking pathogenic mutations in rare uncharacterized monogenic diseases. Typically, before DNA sequencing, protein-coding regions are enriched from patient genomic DNA, representing either the entire genome ("exome sequencing") or selected mapped candidate loci. Sequence variants, identified as differences between the patient's and the human genome reference sequences, are then filtered according to various quality parameters. Changes are screened against datasets of known polymorphisms, such as dbSNP and the 1000 Genomes Project, in the effort to narrow the list of candidate causative variants. An increasing number of commercial services now offer to both generate and align NGS data to a reference genome. This potentially allows small groups with limited computing infrastructure and informatics skills to utilize this technology. However, the capability to effectively filter and assess sequence variants is still an important bottleneck in the identification of deleterious sequence variants in both research and diagnostic settings. We have developed an approach to this problem comprising a user-friendly suite of programs that can interactively analyze, filter and screen data from enrichment-capture NGS data. These programs ("Agile Suite") are particularly suitable for small-scale gene discovery or for diagnostic analysis.

Founder mutations and genotype-phenotype correlations in Meckel-Gruber syndrome and associated ciliopathies

Background

Meckel-Gruber syndrome (MKS) is an autosomal recessive lethal condition that is a ciliopathy. MKS has marked phenotypic variability and genetic heterogeneity, with mutations in nine genes identified as causative to date.

Methods

Families diagnosed with Meckel-Gruber syndrome were recruited for research studies following informed consent. DNA samples were analyzed by microsatellite genotyping and direct Sanger sequencing.

Results

We now report the genetic analyses of 87 individuals from 49 consanguineous and 19 non-consanguineous families in an unselected cohort with reported MKS, or an associated severe ciliopathy in a kindred. Linkage and/or direct sequencing were prioritized for seven MKS genes (MKS1, TMEM216, TMEM67/MKS3, RPGRIP1L, CC2D2A, CEP290 and TMEM237) selected on the basis of reported frequency of mutations or ease of analysis. We have identified biallelic mutations in 39 individuals, of which 13 mutations are novel and previously unreported. We also confirm general genotype-phenotype correlations.

Conclusions

TMEM67 was the most frequently mutated gene in this cohort, and we confirm two founder splice-site mutations (c.1546 + 1 G > A and c.870-2A > G) in families of Pakistani ethnic origin. In these families, we have also identified two separate founder mutations for RPGRIP1L (c. 1945 C > T p.R649X) and CC2D2A (c. 3540delA p.R1180SfsX6). Two missense mutations in TMEM67 (c. 755 T > C p.M252T, and c. 1392 C > T p.R441C) are also probable founder mutations. These findings will contribute to improved genetic diagnosis and carrier testing for affected families, and imply the existence of further genetic heterogeneity in this syndrome.

Prostaglandin transporter mutations cause pachydermoperiostosis with myelofibrosis

Pachydermoperiostosis, or primary hypertrophic osteoarthropathy (PHO), is an inherited multisystem disorder, whose features closely mimic the reactive osteoarthropathy that commonly accompanies neoplastic and inflammatory pathologies. We previously described deficiency of the prostaglandin-degrading enzyme 15-hydroxyprostaglandin dehydrogenase (HPGD) as a cause of this condition, implicating elevated circulating prostaglandin E(2) (PGE(2)) as causative of PHO, and perhaps also as the principal mediator of secondary HO. However, PHO is genetically heterogeneous. Here, we use whole-exome sequencing to identify recessive mutations of the prostaglandin transporter SLCO2A1, in individuals lacking HPGD mutations. We performed exome sequencing of four probands with severe PHO, followed by conventional mutation analysis of SLCO2A1 in nine others. Biallelic SLCO2A1 mutations were identified in 12 of the 13 families. Affected individuals had elevated urinary PGE(2), but unlike HPGD-deficient patients, also excreted considerable quantities of the PGE(2) metabolite, PGE-M. Clinical differences between the two groups were also identified, notably that SLCO2A1-deficient individuals have a high frequency of severe anemia due to myelofibrosis. These findings reinforce the key role of systemic or local prostaglandin excess as the stimulus to HO. They also suggest that the induction or maintenance of hematopoietic stem cells by prostaglandin may depend upon transporter activity.

Evolutionarily assembled cis-regulatory module at a human ciliopathy locus

Neighboring genes are often coordinately expressed within cis-regulatory modules, but evidence that nonparalogous genes share functions in mammals is lacking. Here, we report that mutation of either TMEM138 or TMEM216 causes a phenotypically indistinguishable human ciliopathy, Joubert syndrome. Despite a lack of sequence homology, the genes are aligned in a head-to-tail configuration and joined by chromosomal rearrangement at the amphibian-to-reptile evolutionary transition. Expression of the two genes is mediated by a conserved regulatory element in the noncoding intergenic region. Coordinated expression is important for their interdependent cellular role in vesicular transport to primary cilia. Hence, during vertebrate evolution of genes involved in ciliogenesis, nonparalogous genes were arranged to a functional gene cluster with shared regulatory elements.

CEP41 is mutated in Joubert syndrome and is required for tubulin glutamylation at the cilium

Tubulin glutamylation is a post-translational modification that occurs predominantly in the ciliary axoneme and has been suggested to be important for ciliary function. However, its relationship to disorders of the primary cilium, termed ciliopathies, has not been explored. Here we mapped a new locus for Joubert syndrome (JBTS), which we have designated as JBTS15, and identified causative mutations in CEP41, which encodes a 41-kDa centrosomal protein. We show that CEP41 is localized to the basal body and primary cilia, and regulates ciliary entry of TTLL6, an evolutionarily conserved polyglutamylase enzyme. Depletion of CEP41 causes ciliopathy-related phenotypes in zebrafish and mice and results in glutamylation defects in the ciliary axoneme. Our data identify CEP41 mutations as a cause of JBTS and implicate tubulin post-translational modification in the pathogenesis of human ciliary dysfunction.

TMEM237 is mutated in individuals with a Joubert syndrome related disorder and expands the role of the TMEM family at the ciliary transition zone

Joubert syndrome related disorders (JSRDs) have broad but variable phenotypic overlap with other ciliopathies. The molecular etiology of this overlap is unclear but probably arises from disrupting common functional module components within primary cilia. To identify additional module elements associated with JSRDs, we performed homozygosity mapping followed by next-generation sequencing (NGS) and uncovered mutations in TMEM237 (previously known as ALS2CR4). We show that loss of the mammalian TMEM237, which localizes to the ciliary transition zone (TZ), results in defective ciliogenesis and deregulation of Wnt signaling. Furthermore, disruption of Danio rerio (zebrafish) tmem237 expression produces gastrulation defects consistent with ciliary dysfunction, and Caenorhabditis elegans jbts-14 genetically interacts with nphp-4, encoding another TZ protein, to control basal body-TZ anchoring to the membrane and ciliogenesis. Both mammalian and C. elegans TMEM237/JBTS-14 require RPGRIP1L/MKS5 for proper TZ localization, and we demonstrate additional functional interactions between C. elegans JBTS-14 and MKS-2/TMEM216, MKSR-1/B9D1, and MKSR-2/B9D2. Collectively, our findings integrate TMEM237/JBTS-14 in a complex interaction network of TZ-associated proteins and reveal a growing contribution of a TZ functional module to the spectrum of ciliopathy phenotypes.

A meckelin-filamin A interaction mediates ciliogenesis

MKS3, encoding the transmembrane receptor meckelin, is mutated in Meckel-Gruber syndrome (MKS), an autosomal-recessive ciliopathy. Meckelin localizes to the primary cilium, basal body and elsewhere within the cell. Here, we found that the cytoplasmic domain of meckelin directly interacts with the actin-binding protein filamin A, potentially at the apical cell surface associated with the basal body. Mutations in FLNA, the gene for filamin A, cause periventricular heterotopias. We identified a single consanguineous patient with an MKS-like ciliopathy that presented with both MKS and cerebellar heterotopia, caused by an unusual in-frame deletion mutation in the meckelin C-terminus at the region of interaction with filamin A. We modelled this mutation and found it to abrogate the meckelin-filamin A interaction. Furthermore, we found that loss of filamin A by siRNA knockdown, in patient cells, and in tissues from Flna(Dilp2) null mouse embryos results in cellular phenotypes identical to those caused by meckelin loss, namely basal body positioning and ciliogenesis defects. In addition, morpholino knockdown of flna in zebrafish embryos significantly increases the frequency of dysmorphology and severity of ciliopathy developmental defects caused by mks3 knockdown. Our results suggest that meckelin forms a functional complex with filamin A that is disrupted in MKS and causes defects in neuronal migration and Wnt signalling. Furthermore, filamin A has a crucial role in the normal processes of ciliogenesis and basal body positioning. Concurrent with these processes, the meckelin-filamin A signalling axis may be a key regulator in maintaining correct, normal levels of Wnt signalling.

Mutations in MEGF10, a regulator of satellite cell myogenesis, cause early onset myopathy, areflexia, respiratory distress and dysphagia (EMARDD)

Infantile myopathies with diaphragmatic paralysis are genetically heterogeneous, and clinical symptoms do not assist in differentiating between them. We used phased haplotype analysis with subsequent targeted exome sequencing to identify MEGF10 mutations in a previously unidentified type of infantile myopathy with diaphragmatic weakness, areflexia, respiratory distress and dysphagia. MEGF10 is highly expressed in activated satellite cells and regulates their proliferation as well as their differentiation and fusion into multinucleated myofibers, which are greatly reduced in muscle from individuals with early onset myopathy, areflexia, respiratory distress and dysphagia.

Next generation sequencing identifies mutations in Atonal homolog 7 (ATOH7) in families with global eye developmental defects

The atonal homolog 7 (ATOH7) gene encodes a transcription factor involved in determining the fate of retinal progenitor cells and is particularly required for optic nerve and ganglion cell development. Using a combination of autozygosity mapping and next generation sequencing, we have identified homozygous mutations in this gene, p.E49V and p.P18RfsX69, in two consanguineous families diagnosed with multiple ocular developmental defects, including severe vitreoretinal dysplasia, optic nerve hypoplasia, persistent fetal vasculature, microphthalmia, congenital cataracts, microcornea, corneal opacity and nystagmus. Most of these clinical features overlap with defects in the Norrin/β-catenin signalling pathway that is characterized by dysgenesis of the retinal and hyaloid vasculature. Our findings document Mendelian mutations within ATOH7 and imply a role for this molecule in the development of structures at the front as well as the back of the eye. This work also provides further insights into the function of ATOH7, especially its importance in retinal vascular development and hyaloid regression.

Illuminator, a desktop program for mutation detection using short-read clonal sequencing

Current methods for sequencing clonal populations of DNA molecules yield several gigabases of data per day, typically comprising reads of < 100 nt. Such datasets permit widespread genome resequencing and transcriptome analysis or other quantitative tasks. However, this huge capacity can also be harnessed for the resequencing of smaller (gene-sized) target regions, through the simultaneous parallel analysis of multiple subjects, using sample "tagging" or "indexing". These methods promise to have a huge impact on diagnostic mutation analysis and candidate gene testing. Here we describe a software package developed for such studies, offering the ability to resolve pooled samples carrying barcode tags and to align reads to a reference sequence using a mutation-tolerant process. The program, Illuminator, can identify rare sequence variants, including insertions and deletions, and permits interactive data analysis on standard desktop computers. It facilitates the effective analysis of targeted clonal sequencer data without dedicated computational infrastructure or specialized training.

TTC21B contributes both causal and modifying alleles across the ciliopathy spectrum

Ciliary dysfunction leads to a broad range of overlapping phenotypes, termed collectively as ciliopathies. This grouping is underscored by genetic overlap, where causal genes can also contribute modifying alleles to clinically distinct disorders. Here we show that mutations in TTC21B/IFT139, encoding a retrograde intraflagellar transport (IFT) protein, cause both isolated nephronophthisis (NPHP) and syndromic Jeune Asphyxiating Thoracic Dystrophy (JATD). Moreover, although systematic medical resequencing of a large, clinically diverse ciliopathy cohort and matched controls showed a similar frequency of rare changes, in vivo and in vitro evaluations unmasked a significant enrichment of pathogenic alleles in cases, suggesting that TTC21B contributes pathogenic alleles to ∼5% of ciliopathy patients. Our data illustrate how genetic lesions can be both causally associated with diverse ciliopathies, as well as interact in trans with other disease-causing genes, and highlight how saturated resequencing followed by functional analysis of all variants informs the genetic architecture of disorders.

Molecular genetics and pathogenic mechanisms for the severe ciliopathies: insights into neurodevelopment and pathogenesis of neural tube defects

Meckel-Gruber syndrome (MKS) is a severe autosomal recessively inherited disorder characterized by developmental defects of the central nervous system that comprise neural tube defects that most commonly present as occipital encephalocele. MKS is considered to be the most common syndromic form of neural tube defect. MKS is genetically heterogeneous with six known disease genes: MKS1, MKS2/TMEM216, MKS3/TMEM67, RPGRIP1L, CEP290, and CC2D2A with the encoded proteins all implicated in the correct function of primary cilia. Primary cilia are microtubule-based organelles that project from the apical surface of most epithelial cell types. Recent progress has implicated the involvement of cilia in the Wnt and Shh signaling pathways and has led to an understanding of their role in normal mammalian neurodevelopment. The aim of this review is to provide an overview of the molecular genetics of the human disorder, and to assess recent insights into the etiology and molecular cell biology of severe ciliopathies from mammalian animal models of MKS.

Nesprin-2 interacts with meckelin and mediates ciliogenesis via remodelling of the actin cytoskeleton

Meckel-Gruber syndrome (MKS) is a severe autosomal recessively inherited disorder caused by mutations in genes that encode components of the primary cilium and basal body. Here we show that two MKS proteins, MKS1 and meckelin, that are required for centrosome migration and ciliogenesis interact with actin-binding isoforms of nesprin-2 (nuclear envelope spectrin repeat protein 2, also known as Syne-2 and NUANCE). Nesprins are important scaffold proteins for maintenance of the actin cytoskeleton, nuclear positioning and nuclear-envelope architecture. However, in ciliated-cell models, meckelin and nesprin-2 isoforms colocalized at filopodia prior to the establishment of cell polarity and ciliogenesis. Loss of nesprin-2 and nesprin-1 shows that both mediate centrosome migration and are then essential for ciliogenesis, but do not otherwise affect apical-basal polarity. Loss of meckelin (by siRNA and in a patient cell-line) caused a dramatic remodelling of the actin cytoskeleton, aberrant localization of nesprin-2 isoforms to actin stress-fibres and activation of RhoA signalling. These findings further highlight the important roles of the nesprins during cellular and developmental processes, particularly in general organelle positioning, and suggest that a mechanistic link between centrosome positioning, cell polarity and the actin cytoskeleton is required for centrosomal migration and is essential for early ciliogenesis.

Loss of the metalloprotease ADAM9 leads to cone-rod dystrophy in humans and retinal degeneration in mice

Cone-rod dystrophy (CRD) is an inherited progressive retinal dystrophy affecting the function of cone and rod photoreceptors. By autozygosity mapping, we identified null mutations in the ADAM metallopeptidase domain 9 (ADAM9) gene in four consanguineous families with recessively inherited early-onset CRD. We also found reduced photoreceptor responses in Adam9 knockout mice, previously reported to be asymptomatic. In 12-month-old knockout mice, photoreceptors appear normal, but the apical processes of the retinal pigment epithelium (RPE) cells are disorganized and contact between photoreceptor outer segments (POSs) and the RPE apical surface is compromised. In 20-month-old mice, there is clear evidence of progressive retinal degeneration with disorganized POS and thinning of the outer nuclear layer (ONL) in addition to the anomaly at the POS-RPE junction. RPE basal deposits and macrophages were also apparent in older mice. These findings therefore not only identify ADAM9 as a CRD gene but also identify a form of pathology wherein retinal disease first manifests at the POS-RPE junction.

Mutations in the cilia gene ARL13B lead to the classical form of Joubert syndrome

Joubert syndrome (JS) and related disorders are a group of autosomal-recessive conditions sharing the "molar tooth sign" on axial brain MRI, together with cerebellar vermis hypoplasia, ataxia, and psychomotor delay. JS is suggested to be a disorder of cilia function and is part of a spectrum of disorders involving retinal, renal, digital, oral, hepatic, and cerebral organs. We identified mutations in ARL13B in two families with the classical form of JS. ARL13B belongs to the Ras GTPase family, and in other species is required for ciliogenesis, body axis formation, and renal function. The encoded Arl13b protein was expressed in developing murine cerebellum and localized to the cilia in primary neurons. Overexpression of human wild-type but not patient mutant ARL13B rescued the Arl13b scorpion zebrafish mutant. Thus, ARL13B has an evolutionarily conserved role mediating cilia function in multiple organs.

Recent advances in the molecular pathology, cell biology and genetics of ciliopathies

Primary cilia have a broad tissue distribution and are present on most cell types in the human body. Until recently, they were considered to be redundant organelles, but progress over the past 5 years has led to an understanding of their role in normal mammalian development. The class of inherited disorders that involve aberrant ciliary function are known as ciliopathies, and although their range of severity can vary, they share some common and unexpected clinical phenotypes. The aim of this review is to assess recent insights into the structure, function and formation of primary cilia, and relate this to the pathology, molecular genetics and cell biology of the ciliopathies.