Deciphering the impact of coding and non-coding SCN1A gene variants on RNA splicing.

Variants that disrupt normal pre-mRNA splicing are increasingly being recognized as a major cause of monogenic disorders. The SCN1A gene, a key epilepsy gene that is linked to various epilepsy phenotypes, is no exception. Approximately 10% of all reported variants in the SCN1A gene are designated as splicing variants, with many located outside of the canonical donor and acceptor splice sites, and most have not been functionally investigated. However, given its restricted expression pattern, functional analysis of splicing variants in the SCN1A gene could not be routinely performed. In this study, we conducted a comprehensive analysis of all reported SCN1A variants and their potential to impact SCN1A splicing and conclude that splicing variants are substantially misannotated and under-represented. We created a splicing reporter system consisting of 18 splicing vectors covering all 26 protein-coding exons with different genomic contexts and several promoters of varying strengths in order to reproduce the wild-type splicing pattern of the SCN1A gene, revealing cis-regulatory elements essential for proper recognition of SCN1A exons. Functional analysis of 95 SCN1A variants was carried out, including all 68 intronic variants reported in the literature, located outside of the splice sites canonical dinucleotides; 21 exonic variants of different classes (synonymous, missense, nonsense and in-frame deletion) and six variants observed in patients with epilepsy. Interestingly, almost 20% of tested intronic variants had no influence on SCN1A splicing, despite being reported as causative in the literature. Moreover, we confirmed that the majority of predicted exonic variants affect splicing unravelling their true molecular mechanism. We used functional data to perform genotype-phenotype correlation, revealing distinct distribution patterns for missense and splice-affecting 'missense' variants and observed no difference in the phenotype severity of variants leading to in-frame and out-of-frame isoforms, indicating that the Nav1.1 protein is highly intolerant to structural variations. Our work demonstrates the importance of functional analysis in proper variant annotation and provides a tool for high-throughput delineation of splice-affecting variants in SCN1A in a whole-gene manner.

Annotation of uORFs in the OMIM genes allows to reveal pathogenic variants in 5'UTRs.

An increasing number of studies emphasize the role of non-coding variants in the development of hereditary diseases. However, the interpretation of such variants in clinical genetic testing still remains a critical challenge due to poor knowledge of their pathogenicity mechanisms. It was previously shown that variants in 5'-untranslated regions (5'UTRs) can lead to hereditary diseases due to disruption of upstream open reading frames (uORFs). Here, we performed a manual annotation of upstream translation initiation sites (TISs) in human disease-associated genes from the OMIM database and revealed ∼4.7 thousand of TISs related to uORFs. We compared our TISs with the previous studies and provided a list of 'high confidence' uORFs. Using a luciferase assay, we experimentally validated the translation of uORFs in the ETFDH, PAX9, MAST1, HTT, TTN,GLI2 and COL2A1 genes, as well as existence of N-terminal CDS extension in the ZIC2 gene. Besides, we created a tool to annotate the effects of genetic variants located in uORFs. We revealed the variants from the HGMD and ClinVar databases that disrupt uORFs and thereby could lead to Mendelian disorders. We also showed that the distribution of uORFs-affecting variants differs between pathogenic and population variants. Finally, drawing on manually curated data, we developed a machine-learning algorithm that allows us to predict the TISs in other human genes.

Functional investigation of SCN1A deep-intronic variants activating poison exons inclusion.

Dravet syndrome is a devastating epileptic syndrome characterized by intractable epilepsy with an early age of onset, regression of developmental milestones, ataxia, and motor deficits. Loss-of-function pathogenic variants in the SCN1A gene are found in the majority of patients with Dravet syndrome; however, a significant number of patients remain undiagnosed even after comprehensive genetic testing. Previously, it was shown that intronic elements in the SCN1A gene called poison exons can incorporate into SCN1A mRNA, leading to haploinsufficiency and potentially causing Dravet syndrome. Here, we developed a splicing reporter assay for all described poison exons of the SCN1A gene and validated it using previously reported and artificially introduced variants. Overall, we tested 18 deep-intronic single nucleotide variants and one complex allele in the SCN1A gene. Our approach is capable of evaluating the effect of both variants affecting cis-regulatory sequences and splice-site variants, with the potential to functionally annotate every possible variant within these elements. Moreover, using antisense-modified uridine-rich U7 small nuclear RNAs, we were able to block poison exon incorporation in mutant constructs, an approach that could be used as a promising therapeutic intervention in Dravet syndrome patients with deep-intronic variants.

Splice-site variant in the RPS7 5'-UTR leads to a decrease in the mRNA level and development of Diamond-Blackfan anemia.

Diamond-Blackfan anemia (DBA) is an inherited bone marrow failure syndrome characterized by erythroid aplasia. Pathogenic variants in ribosomal protein (RP) genes, GATA1, TSR2, and EPO, are considered to be the etiology of DBA. Variants in 5'-untranslated regions (UTRs) of these genes are poorly studied and can complicate the variant interpretation. We investigated the functional consequences NM_001011.4:c.-19 + 1G > T variant in the donor splice-site of the RPS7 5'-UTR. This variant was found in a family where two sons with DBA were carriers. Father, who also had this variant, developed myelodysplastic syndrome, which caused his death. Search for candidate causal variants and copy number variations in DBA-associated genes left RPS7 variant as the best candidate. Trio whole exome sequencing analysis revealed no pathogenic variants in other genes. Functional analysis using luciferase expression system revealed that this variant leads to disruption of splicing. Also, a decrease in the levels of mRNA and protein expression was detected. In conclusion, the established consequences of 5'-UTR splice-site variant c.-19 + 1G > T in the RPS7 gene provide evidence that it is likely pathogenic.

Evaluation of Pathogenicity and Causativity of Variants in the MPZ and SH3TC2 Genes in a Family Case of Hereditary Peripheral Neuropathy.

The implementation of NGS methods into clinical practice allowed researchers effectively to establish the molecular cause of a disorder in cases of a genetically heterogeneous pathology. In cases of several potentially causative variants, we need additional analysis that can help in choosing a proper causative variant. In the current study, we described a family case of hereditary motor and sensory neuropathy (HMSN) type 1 (Charcot-Marie-Tooth disease). DNA analysis revealed two variants in the SH3TC2 gene (c.279G>A and c.1177+5G>A), as well as a previously described variant c.449-9C>T in the MPZ gene, in a heterozygous state. This family segregation study was incomplete because of the proband's father's unavailability. To evaluate the variants' pathogenicity, minigene splicing assay was carried out. This study showed no effect of the MPZ variant on splicing, but the c.1177+5G>A variant in the SH3TC2 gene leads to the retention of 122 nucleotides from intron 10 in the RNA sequence, causing a frameshift and an occurrence of a premature stop codon (NP_078853.2:p.Ala393GlyfsTer2).

Clinical and Genetic Characteristics of Calvarial Doughnut Lesions with Bone Fragility in Three Families with a Reccurent SGMS2 Gene Variant.

Calvarial doughnut lesions (CDL) with bone fragility with or without spondylometaphyseal dysplasia (MIM: #126550) is a rare autosomal dominant skeletal disorder characterized by low bone mineral density, spinal and peripheral fractures, and specific sclerotic lesions of the cranial bones. In the current classification of skeletal disorders, the disease is included in the group of bone fragility disorders along with osteogenesis imperfecta. The disease is caused by pathogenic variants in the SGMS2 gene, the protein product of which is sphingomyelin synthase 2, which primarily contributes to sphingomyelin (SM) synthesis-the main lipid component of the plasma membrane essential for bone mineralization. To date, 15 patients from eight families with CDL with bone fragility have been described in the literature, and a recurrent variant c.148C>T (p.Arg50Ter) in the SGMS2 gene has been identified, which was found in patients from six families. We diagnosed the disease in 11 more patients from three unrelated families, caused by the same heterozygous nonsense variant c.148C>T (p.Arg50Ter) in the SGMS2 gene. Our results show wide interfamilial and intrafamilial phenotypic variability in patients with a detected recurrent variant in the SGMS2 gene, the presence of which must be taken into consideration in the diagnosis of the disease. The primary analysis of this variant will contribute to optimal molecular genetic diagnostics, which can reduce diagnostic costs and time.

Clinical and Functional Characteristics of the E92K CFTR Gene Variant in the Russian and Turkish Population of People with Cystic Fibrosis.

The pathogenic variant E92K (c.274G > A) of the CFTR gene is rare in America and Europe, but it is common for people with cystic fibrosis from Russia and Turkey. We studied the effect of the E92K genetic variant on the CFTR function. The function of the CFTR channel was studied using the intestinal current measurements (ICM) method. The effects of CFTR modulators on the restoration of the CFTR function were studied in the model of intestinal organoids. To assess the effect of E92K on pre-mRNA splicing, the RT-PCR products obtained from patients' intestinal organoid cultures were analyzed. Patients with the genetic variant E92K are characterized by an older age of diagnosis compared to homozygotes F508del and a high frequency of pancreatic sufficiency. The results of the sweat test and the ICM method showed partial preservation of the function of the CFTR channel. Functional analysis of CFTR gene expression revealed a weak effect of the E92K variant on mRNA-CFTR splicing. Lumacaftor (VX-809) has been shown to restore CFTR function in an intestinal organoid model, which allows us to consider the E92K variant as a promising target for therapy with CFTR correctors.

Congenital myopathy as a new phenotype caused by two undescribed variants in ASCC1 gene.

We present a patient with congenital myopathy and an inborn epiphysiolysis of the ulna. Whole-exome sequencing analysis revealed two novel mutations in Activation Signal Cointegrator Complex 1 (ASCC1) gene in a compound heterozygous state-a splicing variant c.395-2A>G and a deletion of the first two coding exons. Homozygous and compound heterozygous LoF variants in ASCC1 gene lead to a severe phenotype of spinal muscular atrophy with congenital bone fractures 2 (SMABF2). All patients described to date presented with a severe muscular hypotony, inborn fractures, and passed away shortly after birth while our proband had moderate hypotony, no fractures, but epiphysiolysis and he was 3.5 years old at the time of examination. To explain the phenotype of our patient, we performed an RNA analysis of all family members. We discovered that the c.395-2A>G variant results in two aberrant mRNA isoforms. We also validated the deletion of two exons in ASCC1 gene that lead to the increased expression of this truncated transcript by 1.8 times. To investigate the possible impact of this deletion on the phenotype we predicted a new Kozak sequence in exon 4 that could lead to the formation of a truncated protein with shortened KH domain and a full RNA ligase-like domain. We suggest that this unexpectedly different phenotype of the proband with ASCC1-related disorder could be explained by the presence of the truncated protein with an increased expression.

Meier-Gorlin Syndrome: Clinical Misdiagnosis, Genetic Testing and Functional Analysis of ORC6 Mutations and the Development of a Prenatal Test.

Meier−Gorlin syndrome (MGS) is a rare genetic developmental disorder that causes primordial proportional dwarfism, microtia, the absence of or hypoplastic patellae and other skeletal anomalies. Skeletal symptoms overlapping with other syndromes make MGS difficult to diagnose clinically. We describe a 3-year-old boy with short stature, recurrent respiratory infections, short-rib dysplasia, tower head and facial dysmorphisms who was admitted to the Tomsk Genetic Clinic to verify a clinical diagnosis of Jeune syndrome. Clinical exome sequencing revealed two variants (compound heterozygosity) in the ORC6 gene: c.2T>C(p.Met1Thr) and c.449+5G>A. In silico analysis showed the pathogenicity of these two mutations and predicted a decrease in donor splicing site strength for c.449+5G>A. An in vitro minigene assay indicated that variant c.449+5G>A causes complete skipping of exon 4 in the ORC6 gene. The parents requested urgent prenatal testing for MGS for the next pregnancy, but it ended in a miscarriage. Our results may help prevent MGS misdiagnosis in the future. We also performed in silico and functional analyses of ORC6 mutations and developed a restriction fragment length polymorphism and haplotype-based short-tandem-repeat assay for prenatal genetic testing for MGS. These findings should elucidate MGS etiology and improve the quality of genetic counselling for affected families.

Hyperammonemia in Russia Due to Carbonic Anhydrase VA Deficiency Caused by Homozygous Mutation p.Lys185Lys (c.555G>A) of the CA5A Gene

Hyperammonemia due to carbonic anhydrase VA deficiency (OMIM# 615751) is a rare, life-threatening hereditary disease caused by biallelic mutations in the CA5A gene, presenting as encephalopathic hyperammonemia of unexplained origin during the neonatal period and infancy. Here, we present a detailed description of a 5-year-old patient with the homozygous mutation p.Lys185Lys (c.555G>A) in the CA5A gene. This variant was previously described by van Karnebeek et al. in 2014 in a boy of Russian origin. We found a high frequency of carriers of this mutation in Russia; 1:213, which is 7 times higher than the expected frequency calculated based on data on Western European populations. Thus, targeted testing for the mutation p.Lys185Lys (c.555G>A) in the CA5A gene should be useful for early detection by selective screening in neonatal intensive care units.

Upstream ORF frameshift variants in the PAX6 5'UTR cause congenital aniridia.

Congenital aniridia (AN) is a severe autosomal dominant panocular disorder associated with pathogenic variants in the PAX6 gene. Previously, we performed a molecular genetic study of a large cohort of Russian patients with AN and revealed four noncoding nucleotide variants in the PAX6 5'UTR. 14 additional PAX6-5'UTR variants were also reported in the literature, but the mechanism of their pathogenicity remained unclear. In the present study, we experimentally analyze five patient-derived PAX6 5'UTR-variants: four variants that we identified in Russian patients (c.-128-2delA, c.-125dupG, c.-122dupG, c.-118_-117del) and one previously reported (c.-52+5G>C). We show that the variants lead to a decrease in the protein translation efficiency, while mRNA expression level is not significantly reduced. Two of these variants also affect splicing. Furthermore, we predict and experimentally validate the presence of an evolutionarily conserved small uORF in the PAX6 5'UTR. All studied variants lead to the frameshift of the uORF, resulting in its extension. This extended out-of-frame uORF overlaps with the downstream CDS and thereby reduces its translation efficiency. We conclude that the uORF frameshift may be the main mechanism of pathogenicity for at least 15 out of 18 known PAX6 5'UTR variants. Moreover, we predict additional uORFs in the PAX6 5'UTR.

Mutation in PHACTR1 associated with multifocal epilepsy with infantile spasms and hypsarrhythmia.

A young boy with multifocal epilepsy with infantile spasms and hypsarrhythmia with minimal organic lesions of brain structures underwent DNA diagnosis using whole-exome sequencing. A heterozygous amino-acid substitution p.L519R in a PHACTR1 gene was identified. PHACTR1 belongs to a protein family of G-actin binding protein phosphatase 1 (PP1) cofactors and was not previously associated with a human disease. The missense single nucleotide variant in the proband was shown to occur de novo in the paternal allele. The mutation was shown in vitro to reduce the affinity of PHACTR1 for G-actin, and to increase its propensity to form complexes with the catalytic subunit of PP1. These properties are associated with altered subcellular localization of PHACTR1 and increased ability to induce cytoskeletal rearrangements. Although the molecular role of the PHACTR1 in neuronal excitability and differentiation remains to be defined, PHACTR1 has been previously shown to be involved in Slack channelopathy pathogenesis, consistent with our findings. We conclude that this activating mutation in PHACTR1 causes a severe type of sporadic multifocal epilepsy in the patient.

Investigation of LINC00493/SMIM26 Gene Suggests Its Dual Functioning at mRNA and Protein Level.

The amount of human long noncoding RNA (lncRNA) genes is comparable to protein-coding; however, only a small number of lncRNAs are functionally annotated. Previously, it was shown that lncRNAs can participate in many key cellular processes, including regulation of gene expression at transcriptional and post-transcriptional levels. The lncRNA genes can contain small open reading frames (sORFs), and recent studies demonstrated that some of the resulting short proteins could play an important biological role. In the present study, we investigate the widely expressed lncRNA LINC00493. We determine the structure of the LINC00493 transcript, its cell localization and influence on cell physiology. Our data demonstrate that LINC00493 has an influence on cell viability in a cell-type-specific manner. Furthermore, it was recently shown that LINC00493 has a sORF that is translated into small protein SMIM26. The results of our knockdown and overexpression experiments suggest that both LINC00493/SMIM26 transcript and protein affect cell viability, but in the opposite manner.

Functional Analysis of the PCCA and PCCB Gene Variants Predicted to Affect Splicing.

It is estimated that up to one-third of all variants causing inherited diseases affect splicing; however, their deleterious effects and roles in disease pathogenesis are often not fully characterized. Given their prevalence and the development of various antisense-based splice-modulating approaches, pathogenic splicing variants have become an important object of genomic medicine. To improve the accuracy of variant interpretation in public mutation repositories, we applied the minigene splicing assay to study the effects of 24 variants that were predicted to affect normal splicing in the genes associated with propionic acidemia (PA)-PCCA and PCCB. As a result, 13 variants (including one missense and two synonymous variants) demonstrated a significant alteration of splicing with the predicted deleterious effect at the protein level and were characterized as spliceogenic loss-of-function variants. The analysis of the available data for the studied variants and application of the American College of Medical Genetics and the Association for Molecular Pathology (ACMG/AMP) guidelines allowed us to precisely classify five of the variants and change the pathogenic status of nine. Using the example of the PA genes, we demonstrated the utility of the minigene splicing assay in the fast and effective assessment of the spliceogenic effect for identified variants and highlight the necessity of their standardized classification.

Complex Transposon Insertion as a Novel Cause of Pompe Disease.

Pompe disease (OMIM#232300) is an autosomal recessive lysosomal storage disorder caused by mutations in the GAA gene. According to public mutation databases, more than 679 pathogenic variants have been described in GAA, none of which are associated with mobile genetic elements. In this article, we report a novel molecular genetic cause of Pompe disease, which could be hardly detected using routine molecular genetic analysis. Whole genome sequencing followed by comprehensive functional analysis allowed us to discover and characterize a complex mobile genetic element insertion deep in the intron 15 of the GAA gene in a patient with infantile onset Pompe disease.

Recessive myotonia congenita caused by a homozygous splice site variant in CLCN1 gene: a case report.

BACKGROUND: Myotonia congenita is a rare neuromuscular disease, which is characterized by a delay in muscle relaxation after evoked or voluntary contraction. Myotonia congenita can be inherited in a dominant (Thomsen disease) and recessive form (Becker disease) and both are caused by pathogenic variants in the CLCN1 gene. Noncanonical splice site variants are often classified as variants of uncertain significance, due to insufficient accuracy of splice-predicting tools. Functional analysis using minigene plasmids is widely used in such cases. Moreover, functional analysis is very useful in investigation of the disease pathogenesis, which is necessary for development of future therapeutic approaches. To our knowledge only one noncanonical splice site variant in the CLCN1 gene was functionally characterized to date. We further contribute to this field by evaluation the molecular mechanism of splicing alteration caused by the c.1582 + 5G > A in a homozygous state. CASE PRESENTATION: We report a clinical case of an affected 6-y.o boy with athletic appearance due to muscle hypertrophy, calf muscle stiffness, cramping and various myotonic signs in a consanguineous family with no history of neuromuscular disorders. The neurological examination showed percussion-activated myotonia in the hands and legs. Plasma creatine kinase enzyme and transaminases levels were normal. Electromyography at the time of examination shows myotonic runs in the upper and lower extremities. CONCLUSIONS: Functional analysis of the variant in a minigene system showed alteration of splicing leading to loss of function, thereby confirming that the variant is pathogenic.

Analysis of candidate genes expected to be essential for melanoma surviving.

INTRODUCTION: Cancers may be treated by selective targeting of the genes vital for their survival. A number of attempts have led to discovery of several genes essential for surviving of tumor cells of different types. In this work, we tried to analyze genes that were previously predicted to be essential for melanoma surviving. Here we present the results of transient siRNA-mediated knockdown of the four of such genes, namely, UNC45A, STK11IP, RHPN2 and ZNFX1, in melanoma cell line A375, then assayed the cells for their viability, proliferation and ability to migrate in vitro. In our study, the knockdown of the genes predicted as essential for melanoma survival does not lead to statistically significant changes in cell viability. On the other hand, for each of the studied genes, mobility assays showed that the knockdown of each of the target genes accelerates the speed of cells migrating. Possible explanation for such counterintuitive results may include insufficiency of the predicting computational models or the necessity of a multiplex knockdown of the genes. AIMS: To examine the hypothesis of essentiality of hypomutated genes for melanoma surviving we have performed knockdown of several genes in melanoma cell line and analyzed cell viability and their ability to migrate. METHODS: Knockdown was performed by siRNAs transfected by Metafectene PRO. The levels of mRNAs before and after knockdown were evaluated by RT-qPCR analysis. Cell viability and proliferation were assessed by MTT assay. Cell migration was assessed by wound healing assay. RESULTS: The knockdown of the genes predicted as essential for melanoma survival does not lead to statistically significant changes in cell viability. On the other hand, for each of the studied genes, mobility assays showed that the knockdown of each of the target genes accelerates the speed of cells migrating. CONCLUSION: Our results do not confirm initial hypothesis that the genes predicted essential for melanoma survival as a matter of fact support the survival of melanoma cells.

Noncompaction cardiomyopathy is caused by a novel in-frame desmin (DES) deletion mutation within the 1A coiled-coil rod segment leading to a severe filament assembly defect.

Mutations in DES, encoding desmin protein, are associated with different kinds of skeletal and/or cardiac myopathies. However, it is unknown, whether DES mutations are associated with left ventricular hypertrabeculation (LVHT). Here, we performed a clinical examination and subsequent genetic analysis in a family, with two individuals presenting LVHT with conduction disease and skeletal myopathy. The genetic analysis revealed a novel small in-frame deletion within the DES gene, p.Q113_L115del, affecting the α-helical rod domain. Immunohistochemistry analysis of explanted myocardial tissue from the index patient revealed an abnormal cytoplasmic accumulation of desmin and a degraded sarcomeric structure. Cell transfection experiments with wild-type and mutant desmin verified the cytoplasmic aggregation and accumulation of mutant desmin. Cotransfection experiments were performed to model the heterozygous state of the patients and revealed a dominant negative effect of the mutant desmin on filament assembly. DES:p.Q113_L115del is classified as a pathogenic mutation associated with dilated cardiomyopathy with prominent LVHT.

LINC01420 RNA structure and influence on cell physiology.

BACKGROUND: It was shown that the major part of human genome is transcribed and produces a large number of long noncoding RNAs (lncRNAs). Today there are many evidences that lncRNAs play important role in the regulation of gene expression during different cellular processes. Moreover, lncRNAs are involved in the development of various human diseases. However, the function of the major part of annotated transcripts is currently unknown, whereas different lncRNAs annotations tend to have low overlap. Recent studies revealed that some lncRNAs have small open reading frames (smORFs), that produce the functional microproteins. However, the question whether the function of such genes is determined by microprotein or RNA itself or both remains open. Thus, the study of new lncRNA genes is important to understanding the functional role of such a heterogeneous class of genes. RESULTS: In the present study, we used reverse transcription PCR and rapid amplification of cDNA ends (RACE) analysis to determine the structure of the LINC01420 transcript. We revealed that LINC01420 has two isoforms that differ in length of the last exon and are localized predominantly in the cytoplasm. We showed that expression of the short isoform is much higher than the long. Besides, MTT and wound-healing assays revealed that LINC01420 inhibited cell migration in human melanoma cell line A375, but does not influence on cell viability. CONCLUSION: During our work, D'Lima et al. found smORF in the first exon of the LINC01420 gene. This smORF produces functional microprotein named non-annotated P-body dissociating polypeptide (NoBody). However, our results provide new facts about LINC01420 transcript and its function.

CTNS mRNA molecular analysis revealed a novel mutation in a child with infantile nephropathic cystinosis: a case report.

BACKGROUND: Cystinosis is an autosomal recessive lysosomal storage disorder characterized by accumulation of cystine in lysosomes throughout the body. Cystinosis is caused by mutations in the CTNS gene that encodes the lysosomal cystine carrier protein cystinosin. CTNS mutations result in either complete absence or reduced cystine transporting function of the protein. The diagnosis of nephropathic cystinosis is generally based on measuring leukocyte cystine level, demonstration of corneal cystine crystals by the slit lamp examination and confirmed by genetic analysis of the CTNS gene. CASE PRESENTATION: A boy born to consanguineous Caucasian parents had the characteristic clinical features of the infantile nephropathic cystinosis including renal Fanconi syndrome (polydipsia/polyuria, metabolic acidosis, hypokalemia, hypophosphatemia, low molecular weight proteinuria, glycosuria, cystine crystals in the cornea) and elevated WBC cystine levels. Initially we performed RFLP analysis of the common in the Northern European population 57-kb deletion of proband's DNA, then a direct Sanger sequencing which revealed no mutations in the coding part of the CTNS gene. To confirm the diagnosis we performed RT-PCR analysis of total RNA obtained from patient-derived fibroblasts in combination with cDNA sequencing. This revealed the skipping of exon 4 and exon 5 in the CTNS in our patient. Therefore, we detected a novel 9-kb homozygous deletion in the CTNS gene at genomic DNA level, spanning region from intron 3 to intron 5. In order to identify the inheritance pattern of the deletion we analyzed DNA of proband's mother and father. Both parents were found to be heterozygous carriers of the CTNS mutation. CONCLUSIONS: Analysis of CTNS gene transcript allowed to identify a large homozygous deletion in the patient with infantile nephropathic cystinosis. Mutational detection at RNA level may be an efficient tool to establish the genetic defect in some cystinosis patients.