Whole-genome sequencing in patients with ciliopathies uncovers a novel recurrent tandem duplication in IFT140

Pathogenic variants in the IFT140 gene and an intriguing clinical presentation in two pediatric patients. Cases report and review of literature

BACKGROUND

The IFT140 gene is one of many genes involved in the synthesis of proteins needed for cilium function. Ciliopathies are a group of disorders associated with the dysfunction of ciliary structures and express as an individual organ system disease as well as multisystem disorders. Dysfunctional cilia typically manifest as pleiotropic clinical features, reflecting their widespread distribution and varied functionality.

CASES PRESENTATION

We present two cases: Case 1, a male with two pathological variations in IFT140 gene, a compound heterozygote, with kidney failure, retinal dystrophy, cardiomyopathy, and situs inversus and Case 2, a female with an IFT140 pathogenic homozygous variant, presented with nephrotic range proteinuria, retinitis pigmentosa, and pseudotumor cerebri.

CONCLUSIONS

As cilia dysfunction is known to cause pleiotropic clinical features due to the presence of cilia in different organs in the body, the clinical picture of the IFT140 mutation is also very heterogeneous. Our cases reveal unprecedented manifestations - LVNC, situs inversus, and pseudotumor cerebri - not previously documented in IFT140 mutation. These findings underscore the importance of genetic screening in ciliopathies.

Autonomous and non-cell autonomous role of cilia in structural birth defects in mice

Ciliopathies are associated with wide spectrum of structural birth defects (SBDs), indicating important roles for cilia in development. Here, we provide novel insights into the temporospatial requirement for cilia in SBDs arising from deficiency in Ift140, an intraflagellar transport (IFT) protein regulating ciliogenesis. Ift140-deficient mice exhibit cilia defects accompanied by wide spectrum of SBDs including macrostomia (craniofacial defects), exencephaly, body wall defects, tracheoesophageal fistula (TEF), randomized heart looping, congenital heart defects (CHDs), lung hypoplasia, renal anomalies, and polydactyly. Tamoxifen inducible CAGGCre-ER deletion of a floxed Ift140 allele between E5.5 to 9.5 revealed early requirement for Ift140 in left-right heart looping regulation, mid to late requirement for cardiac outflow septation and alignment, and late requirement for craniofacial development and body wall closure. Surprisingly, CHD were not observed with 4 Cre drivers targeting different lineages essential for heart development, but craniofacial defects and omphalocele were observed with Wnt1-Cre targeting neural crest and Tbx18-Cre targeting epicardial lineage and rostral sclerotome through which trunk neural crest cells migrate. These findings revealed cell autonomous role of cilia in cranial/trunk neural crest-mediated craniofacial and body wall closure defects, while non-cell autonomous multi-lineage interactions underlie CHD pathogenesis, revealing unexpected developmental complexity for CHD associated with ciliopathies.

Rare IFT140-Associated Phenotype of Cranioectodermal Dysplasia and Features of Diagnostic Journey in Patients with Suspected Ciliopathies

Here we present a patient with a cranioectodermal phenotype associated with pathogenic variants in the IFT140 gene. Most frequently, pathogenic variants in IFT140 correspond to the phenotype of Mainzer-Saldino syndrome. Only four patients have previously been described with this cranioectodermal phenotype and variants in IFT140. In comparison to other IFT140-cranioectodermal patients, our proband had similar skeletal features among with early onset end-stage renal failure that required kidney transplantation but did not have common ophthalmological features such as retinopathy, optic nerve atrophy, or nystagmus. Following exome sequencing, a splicing variant and exons 27-30 tandem duplication were suspected and further validated. The two other patients with Mainzer-Saldino syndrome that we described displayed a typical clinical picture but a special diagnostic journey. In both cases, at first only one pathogenic variant was detected following panel or exome NGS sequencing. Further WGS was performed for one of them where tandem duplication was found. Screening the third patient for the same tandem duplication was successful and revealed the presence of this duplication. Thus, we suggest that the description of the clinical feature polymorphism in a rare IFT140-cranioectodermal phenotype is extremely important for providing genetic counseling for families, as well as the formation of the correct diagnostic path for patients with a variant in IFT140.

Autonomous and non-cell autonomous etiology of ciliopathy associated structural birth defects

Ciliopathies are associated with wide spectrum of structural birth defects (SBD), indicating important roles for cilia in development. Here we provide novel insights into the temporospatial requirement for cilia in SBDs arising from deficiency in Ift140 , an intraflagellar transport protein regulating ciliogenesis. Ift140 deficient mice exhibit cilia defects accompanied by wide spectrum of SBDs including macrostomia (craniofacial defects), exencephaly, body wall defects, tracheoesophageal fistula, randomized heart looping, congenital heart defects (CHD), lung hypoplasia, renal anomalies, and polydactyly. Tamoxifen inducible CAG-Cre deletion of a floxed Ift140 allele between E5.5 to 9.5 revealed early requirement for Ift140 in left-right heart looping regulation, mid to late requirement for cardiac outflow septation and alignment, and late requirement for craniofacial development and body wall closure. Surprisingly, CHD was not observed with four Cre drivers targeting different lineages essential for heart development, but craniofacial defects and omphalocele were observed with Wnt1-Cre targeting neural crest and Tbx18-Cre targeting epicardial lineage and rostral sclerotome through which trunk neural crest cells migrate. These findings revealed cell autonomous role of cilia in cranial/trunk neural crest mediated craniofacial and body wall closure defects, while non-cell autonomous multi-lineage interactions underlie CHD pathogenesis, revealing unexpected developmental complexity for CHD associated with ciliopathy.

WGS Revealed Novel BBS5 Pathogenic Variants, Missed by WES, Causing Ciliary Structure and Function Defects

Bardet-Biedl syndrome (BBS) is an autosomal recessive ciliopathy that affects multiple organs, leading to retinitis pigmentosa, polydactyly, obesity, renal anomalies, cognitive impairment, and hypogonadism. Until now, biallelic pathogenic variants have been identified in at least 24 genes delineating the genetic heterogeneity of BBS. Among those, BBS5 is a minor contributor to the mutation load and is one of the eight subunits forming the BBSome, a protein complex implied in protein trafficking within the cilia. This study reports on a European BBS5 patient with a severe BBS phenotype. Genetic analysis was performed using multiple next-generation sequencing (NGS) tests (targeted exome, TES and whole exome, WES), and biallelic pathogenic variants could only be identified using whole-genome sequencing (WGS), including a previously missed large deletion of the first exons. Despite the absence of family samples, the biallelic status of the variants was confirmed. The BBS5 protein's impact was confirmed on the patient's cells (presence/absence and size of the cilium) and ciliary function (Sonic Hedgehog pathway). This study highlights the importance of WGS and the challenge of reliable structural variant detection in patients' genetic explorations as well as functional tests to assess a variant's pathogenicity.

Identical IFT140 Variants Cause Variable Skeletal Ciliopathy Phenotypes—Challenges for the Accurate Diagnosis

Ciliopathies are rare congenital disorders, caused by defects in the cilium, that cover a broad clinical spectrum. A subgroup of ciliopathies showing significant phenotypic overlap are known as skeletal ciliopathies and include Jeune asphyxiating thoracic dysplasia (JATD), Mainzer-Saldino syndrome (MZSDS), cranioectodermal dysplasia (CED), and short-rib polydactyly (SRP). Ciliopathies are heterogeneous disorders with >187 associated genes, of which some genes are described to cause more than one ciliopathy phenotype. Both the clinical and molecular overlap make accurate diagnosing of these disorders challenging. We describe two unrelated Polish patients presenting with a skeletal ciliopathy who share the same compound heterozygous variants in IFT140 (NM_014,714.4) r.2765_2768del; p.(Tyr923Leufs*28) and exon 27–30 duplication; p.(Tyr1152_Thr1394dup). Apart from overlapping clinical symptoms the patients also show phenotypic differences; patient 1 showed more resemblance to a Mainzer-Saldino syndrome (MZSDS) phenotype, while patient 2 was more similar to the phenotype of cranioectodermal dysplasia (CED). In addition, functional testing in patient-derived fibroblasts revealed a distinct cilium phenotyps for each patient, and strikingly, the cilium phenotype of CED-like patient 2 resembled that of known CED patients. Besides two variants in IFT140, in depth exome analysis of ciliopathy associated genes revealed a likely-pathogenic heterozygous variant in INTU for patient 2 that possibly affects the same IFT-A complex to which IFT140 belongs and thereby could add to the phenotype of patient 2. Taken together, by combining genetic data, functional test results, and clinical findings we were able to accurately diagnose patient 1 with “IFT140-related ciliopathy with MZSDS-like features” and patient 2 with “IFT140-related ciliopathy with CED-like features”. This study emphasizes that identical variants in one ciliopathy associated gene can lead to a variable ciliopathy phenotype and that an in depth and integrated analysis of clinical, molecular and functional data is necessary to accurately diagnose ciliopathy patients.

An Alu-mediated duplication in NMNAT1, involved in NAD biosynthesis, causes a novel syndrome, SHILCA, affecting multiple tissues and organs

We investigated the genetic origin of the phenotype displayed by three children from two unrelated Italian families, presenting with a previously unrecognized autosomal recessive disorder that included a severe form of spondylo-epiphyseal dysplasia, sensorineural hearing loss, intellectual disability and Leber congenital amaurosis (SHILCA), as well as some brain anomalies that were visible at the MRI. Autozygome-based analysis showed that these children shared a 4.76 Mb region of homozygosity on chromosome 1, with an identical haplotype. Nonetheless, whole-exome sequencing failed to identify any shared rare coding variants, in this region or elsewhere. We then determined the transcriptome of patients' fibroblasts by RNA sequencing, followed by additional whole-genome sequencing experiments. Gene expression analysis revealed a 4-fold downregulation of the gene NMNAT1, residing indeed in the shared autozygous interval. Short- and long-read whole-genome sequencing highlighted a duplication involving 2 out of the 5 exons of NMNAT1 main isoform (NM_022787.3), leading to the production of aberrant mRNAs. Pathogenic variants in NMNAT1 have been previously shown to cause non-syndromic Leber congenital amaurosis (LCA). However, no patient with null biallelic mutations has ever been described, and murine Nmnat1 knockouts show embryonic lethality, indicating that complete absence of NMNAT1 activity is probably not compatible with life. The rearrangement found in our cases, presumably causing a strong but not complete reduction of enzymatic activity, may therefore result in an intermediate syndromic phenotype with respect to LCA and lethality.

AnnotSV and knotAnnotSV: a web server for human structural variations annotations, ranking and analysis

With the dramatic increase of pangenomic analysis, Human geneticists have generated large amount of genomic data including millions of small variants (SNV/indel) but also thousands of structural variations (SV) mainly from next-generation sequencing and array-based techniques. While the identification of the complete SV repertoire of a patient is getting possible, the interpretation of each SV remains challenging. To help identifying human pathogenic SV, we have developed a web server dedicated to their annotation and ranking (AnnotSV) as well as their visualization and interpretation (knotAnnotSV) freely available at the following address: https://www.lbgi.fr/AnnotSV/. A large amount of annotations from >20 sources is integrated in our web server including among others genes, haploinsufficiency, triplosensitivity, regulatory elements, known pathogenic or benign genomic regions, phenotypic data. An ACMG/ClinGen compliant prioritization module allows the scoring and the ranking of SV into 5 SV classes from pathogenic to benign. Finally, the visualization interface displays the annotated SV in an interactive way including popups, search fields, filtering options, advanced colouring to highlight pathogenic SV and hyperlinks to the UCSC genome browser or other public databases. This web server is designed for diagnostic and research analysis by providing important resources to the user.

Novel TTLL5 Variants Associated with Cone-Rod Dystrophy and Early-Onset Severe Retinal Dystrophy

Variants of the TTLL5 gene, which encodes tubulin tyrosine ligase-like family member five, are a rare cause of cone dystrophy (COD) or cone-rod dystrophy (CORD). To date, only a few TTLL5 patients have been clinically and genetically described. In this study, we report five patients harbouring biallelic variants of TTLL5. Four adult patients presented either COD or CORD with onset in the late teenage years. The youngest patient had a phenotype of early onset severe retinal dystrophy (EOSRD). Genetic analysis was performed by targeted next generation sequencing of gene panels and assessment of copy number variants (CNV). We identified eight variants, of which six were novel, including two large multiexon deletions in patients with COD or CORD, while the EOSRD patient harboured the novel homozygous p.(Trp640*) variant and three distinct USH2A variants, which might explain the observed rod involvement. Our study highlights the role of TTLL5 in COD/CORD and the importance of large deletions. These findings suggest that COD or CORD patients lacking variants in known genes may harbour CNVs to be discovered in TTLL5, previously undetected by classical sequencing methods. In addition, variable phenotypes in TTLL5-associated patients might be due to the presence of additional gene defects.

Molecular genetics of renal ciliopathies

Renal ciliopathies are a heterogenous group of inherited disorders leading to an array of phenotypes that include cystic kidney disease and renal interstitial fibrosis leading to progressive chronic kidney disease and end-stage kidney disease. The renal tubules are lined with epithelial cells that possess primary cilia that project into the lumen and act as sensory and signalling organelles. Mutations in genes encoding ciliary proteins involved in the structure and function of primary cilia cause ciliopathy syndromes and affect many organ systems including the kidney. Recognised disease phenotypes associated with primary ciliopathies that have a strong renal component include autosomal dominant and recessive polycystic kidney disease and their various mimics, including atypical polycystic kidney disease and nephronophthisis. The molecular investigation of inherited renal ciliopathies often allows a precise diagnosis to be reached where renal histology and other investigations have been unhelpful and can help in determining kidney prognosis. With increasing molecular insights, it is now apparent that renal ciliopathies form a continuum of clinical phenotypes with disease entities that have been classically described as dominant or recessive at both extremes of the spectrum. Gene-dosage effects, hypomorphic alleles, modifier genes and digenic inheritance further contribute to the genetic complexity of these disorders. This review will focus on recent molecular genetic advances in the renal ciliopathy field with a focus on cystic kidney disease phenotypes and the genotypes that lead to them. We discuss recent novel insights into underlying disease mechanisms of renal ciliopathies that might be amenable to therapeutic intervention.

IFT144 and mild retinitis pigmentosa in Mainzer-Saldino syndrome: A new association

Ciliopathies are a wide and heterogeneous group of diseases affecting intraflagellar transport. Among them, Mainzer-Saldino syndrome (MSS) shows phalangeal cone-shaped epiphysis, renal disease and retinal involvement. Short stature, cerebellar ataxia and hepatic fibrosis might also be found. IFT140 is the most commonly reported mutation in MSS. We will report on the case of a patient with a clinical diagnosis of Mainzer-Saldino syndrome due to IFT144 dysfunction. This mutation has not been previously related to MSS but it has been found in other ciliopathies and both syndromic and non-syndromic retinitis pigmentosa. At birth our patient showed trigonocephaly, early progressive renal failure requiring transplant, intrahepatic biliary duct dilation, cone-shaped epiphyses, growth retardation and retinitis pigmentosa with mild ophthalmic impairment. The best corrected visual acuity reached 0.15/0.22 LogMAR. The posterior pole showed abnormal macular reflex, mild vascular attenuation in the periphery and diffuse pigmentary changes. Autofluorescence showed bull's eye signal increase. Computerized optic tomography assessed the absence of external retinal layers in the extrafoveal macula. In conclusion, IFT144 genetic study may be involved in MSS and thus must be considered for diagnosis. Mild ophthalmic symptomatology despite early onset retinitis pigmentosa in the context of MSS has been found in this case caused by IFT144 mutation.

DINTD: Detection and Inference of Tandem Duplications From Short Sequencing Reads

Tandem duplication (TD) is an important type of structural variation (SV) in the human genome and has biological significance for human cancer evolution and tumor genesis. Accurate and reliable detection of TDs plays an important role in advancing early detection, diagnosis, and treatment of disease. The advent of next-generation sequencing technologies has made it possible for the study of TDs. However, detection is still challenging due to the uneven distribution of reads and the uncertain amplitude of TD regions. In this paper, we present a new method, DINTD (Detection and INference of Tandem Duplications), to detect and infer TDs using short sequencing reads. The major principle of the proposed method is that it first extracts read depth and mapping quality signals, then uses the DBSCAN (Density-Based Spatial Clustering of Applications with Noise) algorithm to find the possible TD regions. The total variation penalized least squares model is fitted with read depth and mapping quality signals to denoise signals. A 2D binary search tree is used to search the neighbor points effectively. To further identify the exact breakpoints of the TD regions, split-read signals are integrated into DINTD. The experimental results of DINTD on simulated data sets showed that DINTD can outperform other methods for sensitivity, precision, F1-score, and boundary bias. DINTD is further validated on real samples, and the experiment results indicate that it is consistent with other methods. This study indicates that DINTD can be used as an effective tool for detecting TDs.

Proteasome subunit PSMC3 variants cause neurosensory syndrome combining deafness and cataract due to proteotoxic stress

The ubiquitin–proteasome system degrades ubiquitin‐modified proteins to maintain protein homeostasis and to control signalling. Whole‐genome sequencing of patients with severe deafness and early‐onset cataracts as part of a neurological, sensorial and cutaneous novel syndrome identified a unique deep intronic homozygous variant in the PSMC3 gene, encoding the proteasome ATPase subunit Rpt5, which lead to the transcription of a cryptic exon. The proteasome content and activity in patient's fibroblasts was however unaffected. Nevertheless, patient's cells exhibited impaired protein homeostasis characterized by accumulation of ubiquitinated proteins suggesting severe proteotoxic stress. Indeed, the TCF11/Nrf1 transcriptional pathway allowing proteasome recovery after proteasome inhibition is permanently activated in the patient's fibroblasts. Upon chemical proteasome inhibition, this pathway was however impaired in patient's cells, which were unable to compensate for proteotoxic stress although a higher proteasome content and activity. Zebrafish modelling for knockout in PSMC3 remarkably reproduced the human phenotype with inner ear development anomalies as well as cataracts, suggesting that Rpt5 plays a major role in inner ear, lens and central nervous system development.

Compound heterozygous IFT140 variants in two Polish families with Sensenbrenner syndrome and early onset end-stage renal disease

Background

Sensenbrenner syndrome, which is also known as cranioectodermal dysplasia (CED), is a rare, autosomal recessive ciliary chondrodysplasia characterized by a variety of clinical features including a distinctive craniofacial appearance as well as skeletal, ectodermal, liver and renal anomalies. Progressive renal disease can be life-threatening in this condition. CED is a genetically heterogeneous disorder. Currently, variants in any of six genes (IFT122, WDR35, IFT140, IFT43, IFT52 and WDR19) have been associated with this syndrome. All of these genes encode proteins essential for intraflagellar transport (IFT) a process that is required for cilium assembly, maintenance and function. Intra- and interfamilial clinical variability has been reported in CED, which is consistent with CED’s genetic heterogeneity and is indicative of genetic background effects.

Results

Two male CED patients from two unrelated Polish families were included in this study. Clinical assessment revealed distinctive clinical features of Sensenbrenner syndrome, such as dolichocephaly, shortening of long bones and early onset renal failure. Ectodermal anomalies also included thin hair, short and thin nails, and small teeth in both patients. Next generation sequencing (NGS) techniques were performed in order to determine the underlying genetic cause of the disorder using whole exome sequencing (WES) for patient 1 and a custom NGS-based panel for patient 2. Subsequent qPCR and duplex PCR analysis were conducted for both patients. Genetic analyses identified compound heterozygous variants in the IFT140 gene in both affected individuals. Both patients harbored a tandem duplication variant p.Tyr1152_Thr1394dup on one allele. In addition, a novel missense variant, p.(Leu109Pro), and a previously described p.(Gly522Glu) variant were identified in the second allele in patients 1 and 2, respectively. Segregation analysis of the variants was consistent with the expected autosomal recessive disease inheritance pattern. Both patients had severe renal failure requiring kidney transplantation in early childhood.

Conclusion

The finding of compound heterozygous IFT140 mutations in two unrelated CED patients provide further evidence that IFT140 gene mutations are associated with this syndrome. Our studies confirm that IFT140 changes in patients with CED are associated with early onset end-stage renal disease. Moreover, this report expands our knowledge of the clinical- and molecular genetics of Sensenbrenner syndrome and it highlights the importance of multidisciplinary approaches in the care of CED patients.

Novel IQCE variations confirm its role in postaxial polydactyly and cause ciliary defect phenotype in zebrafish

Polydactyly is one of the most frequent inherited defects of the limbs characterized by supernumerary digits and high-genetic heterogeneity. Among the many genes involved, either in isolated or syndromic forms, eight have been implicated in postaxial polydactyly (PAP). Among those, IQCE has been recently identified in a single consanguineous family. Using whole-exome sequencing in patients with uncharacterized ciliopathies, including PAP, we identified three families with biallelic pathogenic variations in IQCE. Interestingly, the c.895_904del (p.Val301Serfs*8) was found in all families without sharing a common haplotype, suggesting a recurrent mechanism. Moreover, in two families, the systemic phenotype could be explained by additional pathogenic variants in known genes (TULP1, ATP6V1B1). RNA expression analysis on patients' fibroblasts confirms that the dysfunction of IQCE leads to the dysregulation of genes associated with the hedgehog-signaling pathway, and zebrafish experiments demonstrate a full spectrum of phenotypes linked to defective cilia: Body curvature, kidney cysts, left-right asymmetry, misdirected cilia in the pronephric duct, and retinal defects. In conclusion, we identified three additional families confirming IQCE as a nonsyndromic PAP gene. Our data emphasize the importance of taking into account the complete set of variations of each individual, as each clinical presentation could finally be explained by multiple genes.