Copy number variants from 4800 exomes contribute to ~7% of genetic diagnoses in movement disorders, muscle disorders and neuropathies.

Various groups of neurological disorders, including movement disorders and neuromuscular diseases, are clinically and genetically heterogeneous. Diagnostic panel-based exome sequencing is a routine test for these disorders. Despite the success rates of exome sequencing, it results in the detection of causative sequence variants in 'only' 25-30% of cases. Copy number variants (CNVs), i.e. deletion or duplications, explain 10-20% of individuals with multisystemic phenotypes, such as co-existing intellectual disability, but may also have a role in disorders affecting a single system (organ), like neurological disorders with normal intelligence. In this study, CNVs were extracted from clinical exome sequencing reports of 4800 probands primarily with a movement disorder, myopathy or neuropathy. In 88 (~2%) probands, phenotype-matching CNVs were detected, representing ~7% of genetically confirmed cases. CNVs varied from involvement of over 100 genes to single exons and explained X-linked, autosomal dominant, or - recessive disorders, the latter due to either a homozygous CNV or a compound heterozygous CNV with a sequence variant on the other allele. CNVs were detected affecting genes where deletions or duplications are established as a common mechanism, like PRKN (in Parkinson's disease), DMD (in Duchenne muscular dystrophy) and PMP22 (in neuropathies), but also genes in which no intragenic CNVs have been reported to date. Analysis of CNVs as part of panel-based exome sequencing for genetically heterogeneous neurological diseases provides an additional diagnostic yield of ~2% without extra laboratory costs. Therefore it is recommended to perform CNV analysis for movement disorders, muscle disease, neuropathies, or any other single-system disorder.

Systematic analysis of short tandem repeats in 38,095 exomes provides an additional diagnostic yield.

PURPOSE: Expansions of a subset of short tandem repeats (STRs) have been implicated in approximately 30 different human genetic disorders. Despite extensive application of exome sequencing (ES) in routine diagnostic genetic testing, STRs are not routinely identified from these data. METHODS: We assessed diagnostic utility of STR analysis in exome sequencing by applying ExpansionHunter to 2,867 exomes from movement disorder patients and 35,228 other clinical exomes. RESULTS: We identified 38 movement disorder patients with a possible aberrant STR length. Validation by polymerase chain reaction (PCR) and/or repeat-primed PCR technologies confirmed the presence of aberrant expansion alleles for 13 (34%). For seven of these patients the genotype was compatible with the phenotypic description, resulting in a molecular diagnosis. We subsequently tested the remainder of our diagnostic ES cohort, including over 30 clinically and genetically heterogeneous disorders. Optimized manual curation yielded 167 samples with a likely aberrant STR length. Validations confirmed 93/167 (56%) aberrant expansion alleles, of which 48 were in the pathogenic range and 45 in the premutation range. CONCLUSION: Our work provides guidance for the implementation of STR analysis in clinical ES. Our results show that systematic STR evaluation may increase diagnostic ES yield by 0.2%, and recommend making STR evaluation a routine part of ES interpretation in genetic testing laboratories.

KIF1A variants are a frequent cause of autosomal dominant hereditary spastic paraplegia.

Variants in the KIF1A gene can cause autosomal recessive spastic paraplegia 30, autosomal recessive hereditary sensory neuropathy, or autosomal (de novo) dominant mental retardation type 9. More recently, variants in KIF1A have also been described in a few cases with autosomal dominant spastic paraplegia. Here, we describe 20 KIF1A variants in 24 patients from a clinical exome sequencing cohort of 347 individuals with a mostly 'pure' spastic paraplegia. In these patients, spastic paraplegia was slowly progressive and mostly pure, but with a highly variable disease onset (0-57 years). Segregation analyses showed a de novo occurrence in seven cases, and a dominant inheritance pattern in 11 families. The motor domain of KIF1A is a hotspot for disease causing variants in autosomal dominant spastic paraplegia, similar to mental retardation type 9 and recessive spastic paraplegia type 30. However, unlike these allelic disorders, dominant spastic paraplegia was also caused by loss-of-function variants outside this domain in six families. Finally, three missense variants were outside the motor domain and need further characterization. In conclusion, KIF1A variants are a frequent cause of autosomal dominant spastic paraplegia in our cohort (6-7%). The identification of KIF1A loss-of-function variants suggests haploinsufficiency as a possible mechanism in autosomal dominant spastic paraplegia.

Natural Exon Skipping Sets the Stage for Exon Skipping as Therapy for Dystrophic Epidermolysis Bullosa.

Dystrophic epidermolysis bullosa (DEB) is a devastating blistering disease affecting skin and mucous membranes. It is caused by pathogenic variants in the COL7A1 gene encoding type VII collagen, and can be inherited dominantly or recessively. Recently, promising proof-of-principle has been shown for antisense oligonucleotide (AON)-mediated exon skipping as a therapeutic approach for DEB. However, the precise phenotypic effect to be anticipated from exon skipping, and which patient groups could benefit, is not yet clear. To answer these questions, we studied new clinical and molecular data on seven patients from the Dutch EB registry and reviewed the literature on COL7A1 exon skipping variants. We found that phenotypes associated with dominant exon skipping cannot be distinguished from phenotypes caused by other dominant DEB variants. Recessive exon skipping phenotypes are generally relatively mild in the spectrum of recessive DEB. Therefore, for dominant DEB, AON-mediated exon skipping is unlikely to ameliorate the phenotype. In contrast, the overall severity of phenotypes associated with recessive natural exon skipping pivots toward the milder end of the spectrum. Consequently, we anticipate AON-mediated exon skipping for recessive DEB caused by bi-allelic null variants should lead to a clinically relevant improvement of this devastating phenotype.

Panel-Based Exome Sequencing for Neuromuscular Disorders as a Diagnostic Service.

BACKGROUND: Neuromuscular disorders (NMDs) are clinically and genetically heterogeneous. Accurate molecular genetic diagnosis can improve clinical management, provides appropriate genetic counseling and testing of relatives, and allows potential therapeutic trials. OBJECTIVE: To establish the clinical utility of panel-based whole exome sequencing (WES) in NMDs in a population with children and adults with various neuromuscular symptoms. METHODS: Clinical exome sequencing, followed by diagnostic interpretation of variants in genes associated with NMDs, was performed in a cohort of 396 patients suspected of having a genetic cause with a variable age of onset, neuromuscular phenotype, and inheritance pattern. Many had previously undergone targeted gene testing without results. RESULTS: Disease-causing variants were identified in 75/396 patients (19%), with variants in the three COL6-genes (COL6A1, COL6A2 and COL6A3) as the most common cause of the identified muscle disorder, followed by variants in the RYR1 gene. Together, these four genes account for almost 25% of cases in whom a definite genetic cause was identified. Furthermore, likely pathogenic variants and/or variants of uncertain significance were identified in 95 of the patients (24%), in whom functional and/or segregation analysis should be used to confirm or reject the pathogenicity. In 18% of the cases with a disease-causing variant of which we received additional clinical information, we identified a genetic cause in genes of which the associated phenotypes did not match that of the patients. Hence, the advantage of panel-based WES is its unbiased approach. CONCLUSION: Whole exome sequencing, followed by filtering for NMD genes, offers an unbiased approach for the genetic diagnostics of NMD patients. This approach could be used as a first-tier test in neuromuscular disorders with a high suspicion of a genetic cause. With uncertain results, functional testing and segregation analysis are needed to complete the evidence.

Clinical exome sequencing for cerebellar ataxia and spastic paraplegia uncovers novel gene-disease associations and unanticipated rare disorders.

Cerebellar ataxia (CA) and hereditary spastic paraplegia (HSP) are two of the most prevalent motor disorders with extensive locus and allelic heterogeneity. We implemented clinical exome sequencing, followed by filtering data for a 'movement disorders' gene panel, as a generic test to increase variant detection in 76 patients with these disorders. Segregation analysis or phenotypic re-evaluation was utilized to substantiate findings. Disease-causing variants were identified in 9 of 28 CA patients, and 8 of 48 HSP patients. In addition, possibly disease-causing variants were identified in 1 and 8 of the remaining CA and HSP patients, respectively. In 10 patients with CA, the total disease-causing or possibly disease-causing variants were detected in 8 different genes, whereas 16 HSP patients had such variants in 12 different genes. In the majority of cases, the identified variants were compatible with the patient phenotype. Interestingly, in some patients variants were identified in genes hitherto related to other movement disorders, such as TH variants in two siblings with HSP. In addition, rare disorders were uncovered, for example, a second case of HSP caused by a VCP variant. For some patients, exome sequencing results had implications for treatment, exemplified by the favorable L-DOPA treatment in a patient with HSP due to ATP13A2 variants (Parkinson type 9). Thus, clinical exome sequencing in this cohort of CA and HSP patients suggests broadening of disease spectra, revealed novel gene-disease associations, and uncovered unanticipated rare disorders. In addition, clinical exome sequencing results have shown their value in guiding practical patient management.

Recurrent FXYD2 p.Gly41Arg mutation in patients with isolated dominant hypomagnesaemia.

BACKGROUND: Magnesium (Mg(2+)) is an essential ion for cell growth, neuroplasticity and muscle contraction. Blood Mg(2+) levels <0.7 mmol/L may cause a heterogeneous clinical phenotype, including muscle cramps and epilepsy and disturbances in K(+) and Ca(2+) homeostasis. Over the last decade, the genetic origin of several familial forms of hypomagnesaemia has been found. In 2000, mutations in FXYD2, encoding the γ-subunit of the Na(+)-K(+)-ATPase, were identified to cause isolated dominant hypomagnesaemia (IDH) in a large Dutch family suffering from hypomagnesaemia, hypocalciuria and chondrocalcinosis. However, no additional patients have been identified since then. METHODS: Here, two families with hypomagnesaemia and hypocalciuria were screened for mutations in the FXYD2 gene. Moreover, the patients were clinically and genetically characterized. RESULTS: We report a p.Gly41Arg FXYD2 mutation in two families with hypomagnesaemia and hypocalciuria. Interestingly, this is the same mutation as was described in the original study. As in the initial family, several patients suffered from muscle cramps, chondrocalcinosis and epilepsy. Haplotype analysis revealed an overlapping haplotype in all families, suggesting a founder effect. CONCLUSIONS: The recurrent p.Gly41Arg FXYD2 mutation in two new families with IDH confirms that FXYD2 mutation causes hypomagnesaemia. Until now, no other FXYD2 mutations have been reported which could indicate that other FXYD2 mutations will not cause hypomagnesaemia or are embryonically lethal.

Mechanisms of natural gene therapy in dystrophic epidermolysis bullosa.

Revertant mosaicism has been reported in several inherited diseases, including the genetic skin fragility disorder epidermolysis bullosa (EB). Here, we describe the largest cohort of seven patients with revertant mosaicism and dystrophic EB (DEB), associated with mutations in the COL7A1 gene, and determine the underlying molecular mechanisms. We show that revertant mosaicism occurs both in autosomal dominantly and recessively inherited DEB. We found that null mutations resulting in complete loss of collagen VII and severe disease, as well as missense or splice-site mutations associated with some preserved collagen VII function and a milder phenotype, were corrected by revertant mosaicism. The mutation, subtype, and severity of the disease are thus not decisive for the presence of revertant mosaicism. Although collagen VII is synthesized and secreted by both keratinocytes and fibroblasts, evidence for reversion was only found in keratinocytes. The reversion mechanisms included back mutations/mitotic recombinations in 70% of the cases and second-site mutations affecting splicing in 30%. We conclude that revertant mosaicism is more common than previously assumed in patients with DEB, and our findings will have implications for future therapeutic strategies using the patient's naturally corrected cells as a source for cell-based therapies.

Mutations in BICD2, which encodes a golgin and important motor adaptor, cause congenital autosomal-dominant spinal muscular atrophy.

Spinal muscular atrophy (SMA) is a heterogeneous group of neuromuscular disorders caused by degeneration of lower motor neurons. Although functional loss of SMN1 is associated with autosomal-recessive childhood SMA, the genetic cause for most families affected by dominantly inherited SMA is unknown. Here, we identified pathogenic variants in bicaudal D homolog 2 (Drosophila) (BICD2) in three families afflicted with autosomal-dominant SMA. Affected individuals displayed congenital slowly progressive muscle weakness mainly of the lower limbs and congenital contractures. In a large Dutch family, linkage analysis identified a 9q22.3 locus in which exome sequencing uncovered c.320C>T (p.Ser107Leu) in BICD2. Sequencing of 23 additional families affected by dominant SMA led to the identification of pathogenic variants in one family from Canada (c.2108C>T [p.Thr703Met]) and one from the Netherlands (c.563A>C [p.Asn188Thr]). BICD2 is a golgin and motor-adaptor protein involved in Golgi dynamics and vesicular and mRNA transport. Transient transfection of HeLa cells with all three mutant BICD2 cDNAs caused massive Golgi fragmentation. This observation was even more prominent in primary fibroblasts from an individual harboring c.2108C>T (p.Thr703Met) (affecting the C-terminal coiled-coil domain) and slightly less evident in individuals with c.563A>C (p.Asn188Thr) (affecting the N-terminal coiled-coil domain). Furthermore, BICD2 levels were reduced in affected individuals and trapped within the fragmented Golgi. Previous studies have shown that Drosophila mutant BicD causes reduced larvae locomotion by impaired clathrin-mediated synaptic endocytosis in neuromuscular junctions. These data emphasize the relevance of BICD2 in synaptic-vesicle recycling and support the conclusion that BICD2 mutations cause congenital slowly progressive dominant SMA.

The inversa type of recessive dystrophic epidermolysis bullosa is caused by specific arginine and glycine substitutions in type VII collagen.

BACKGROUND: The inversa type of recessive dystrophic epidermolysis bullosa (RDEB-I) is a rare variant of dystrophic epidermolysis bullosa, characterised by blistering in the body flexures, trunk, and mucosa. The cause of this specific distribution is unknown. So far, 20 COL7A1 genotypes have been described in RDEB-I and genotype-phenotype correlations have not been studied extensively. The aim of the study was to gain more insight into the pathophysiology of this intriguing RDEB-I phenotype. METHODS: Twenty Dutch and British RDEB-I patients, and full genotypes in 18 of them, were identified. The literature on RDEB-I genotypes was reviewed and an extensive genotype-phenotype correlation study for RDEB-I was conducted. RESULTS: All 20 patients had generalised blistering at birth and during early infancy. In most patients, the age of transition from generalised to inversa distribution was before the age of 4 years. A spectrum of disease severity, ranging from the mildest 'mucosal only' phenotype to the severest phenotype with limited acral involvement, was noted. The 29 genotypes of these RDEB-I patients and those reported in the literature revealed that RDEB-I is associated with specific recessive arginine and glycine substitutions in the triple helix domain of type VII collagen. DISCUSSION AND CONCLUSION: Why these specific arginine and glycine substitutions cause the inversa distribution remains unknown. It was not possible to identify clear differences in location and nature of substituting amino acids between these mutations and missense mutations causing other RDEB phenotypes. It is hypothesised that the higher skin temperature in the affected areas plays an important role in the pathophysiology of RDEB-I.

Targeted next-generation sequencing of a 12.5 Mb homozygous region reveals ANO10 mutations in patients with autosomal-recessive cerebellar ataxia.

Autosomal-recessive cerebellar ataxias comprise a clinically and genetically heterogeneous group of neurodegenerative disorders. In contrast to their dominant counterparts, unraveling the molecular background of these ataxias has proven to be more complicated and the currently known mutations provide incomplete coverage for genotyping of patients. By combining SNP array-based linkage analysis and targeted resequencing of relevant sequences in the linkage interval with the use of next-generation sequencing technology, we identified a mutation in a gene and have shown its association with autosomal-recessive cerebellar ataxia. In a Dutch consanguineous family with three affected siblings a homozygous 12.5 Mb region on chromosome 3 was targeted by array-based sequence capture. Prioritization of all detected sequence variants led to four candidate genes, one of which contained a variant with a high base pair conservation score (phyloP score: 5.26). This variant was a leucine-to-arginine substitution in the DUF 590 domain of a 16K transmembrane protein, a putative calcium-activated chloride channel encoded by anoctamin 10 (ANO10). The analysis of ANO10 by Sanger sequencing revealed three additional mutations: a homozygous mutation (c.1150_1151del [p.Leu384fs]) in a Serbian family and a compound-heterozygous splice-site mutation (c.1476+1G>T) and a frameshift mutation (c.1604del [p.Leu535X]) in a French family. This illustrates the power of using initial homozygosity mapping with next-generation sequencing technology to identify genes involved in autosomal-recessive diseases. Moreover, identifying a putative calcium-dependent chloride channel involved in cerebellar ataxia adds another pathway to the list of pathophysiological mechanisms that may cause cerebellar ataxia.

Massively parallel sequencing of ataxia genes after array-based enrichment.

Massively parallel sequencing has tremendous diagnostic potential but requires enriched templates for sequencing. Here we report the validation of an array-based sequence capture method in genetically heterogeneous disorders. The model disorder selected was AR ataxia, using five subjects with known mutations and two unaffected controls. The genomic sequences of seven disease genes, together with two control loci were targeted on a 2-Mb sequence-capture array. After enrichment, the patients' DNA samples were analyzed using one-quarter Roche GS FLX Titanium sequencing run, resulting in an average of 65 Mb of sequence data per patient. This was sufficient for an average 25-fold coverage/base in all targeted regions. Enrichment showed high specificity; on average, 80% of uniquely mapped reads were on target. Importantly, this approach enabled automated detection of deletions and hetero- and homozygous point mutations for 6/7 mutant alleles, and greater than 99% accuracy for known SNP variants. Our results also clearly show reduced coverage for sequences in repeat-rich regions, which significantly impacts the reliable detection of genomic variants. Based on these findings we recommend a minimal coverage of 15-fold for diagnostic implementation of this technology. We conclude that massive parallel sequencing of enriched samples enables personalized diagnosis of heterogeneous genetic disorders and qualifies for rapid diagnostic implementation.

The EEC syndrome and SHFM: report of two cases and mutation analysis of p63 gene.

The p63 gene is a transcription factor and a member of the p53 family. Heterozygote mutation of the p63 gene is suggested in a number of human syndromes including limb development and/or ectodermal dysplasia. The EEC syndrome, consisting of ectrodactyly (E), ectodermal dysplasia (E) and cleft lip (C) with or without cleft palate, is the prototype of these syndromes with the presence of heterozygote mutation in the p63 gene in most of the patients. Nonsyndromic split hand/foot malformation (SHFM) is one of the EEC-like syndromes, and the p63 gene mutation was reported in only a few patients. Five different loci have been mapped to date, but the etiology is yet to be explained in the rest of the patients. Here, we report two cases. Case 1, diagnosed with EEC syndrome, had type 2 urogenital sinus and a new heterozygous mutation of 934G>A (D312N) in exon 8 of the p63 gene. Case 2 was diagnosed as SHFM with no mutation in the p63 gene. Genotype and phenotype correlation of these two cases among the reported patients is discussed in this report.

Design and validation of a conformation-sensitive capillary electrophoresis system for mutation identification of the COL7A1 gene with automated peak comparison.

Dystrophic epidermolysis bullosa is a heritable skin disease in which blisters occur because of a defect in type VII collagen resulting from mutations in the COL7A1 gene that is composed of 118 exons. Although a few mutations are specific to certain populations owing to founder effects, and although a few mutational hotspots exist, most mutations are unique to families and can be found scattered throughout the entire COL7A1 gene. This emphasizes the need for a sensitive, reliable, and efficient mutation scanning technique. Therefore, we developed a conformation-sensitive capillary electrophoresis (CSCE) system for COL7A1 mutation scanning. Here we report on the design and validation of this system. The CSCE technique is based on the principle of heteroduplex formation when polymerase chain reaction-amplified DNA fragments containing heterozygous sequence changes are slowly reannealed. These fluorescently labeled fragments have different migration characteristics and can be detected on a multi-capillary automated sequencer. Validation was performed by analysis of 29 known COL7A1 sequence changes, covering 33% of amplicons. After optimization of the conditions, all 29 sequence changes were detected by the CSCE system, irrespective of length or CG-content of amplicons and position of sequence changes, reflecting an analytical sensitivity of 90.2-100% (95% confidence interval). We conclude that this CSCE system is a rapid, reliable, cost-effective, and highly sensitive way of mutation scanning for COL7A1 in a molecular genetics service laboratory.

Design and validation of a conformation sensitive capillary electrophoresis-based mutation scanning system and automated data analysis of the more than 15 kbp-spanning coding sequence of the SACS gene.

In this study, we developed and analytically validated a fully automated, robust confirmation sensitive capillary electrophoresis (CSCE) method to perform mutation scanning of the large SACS gene. This method facilitates a rapid and cost-effective molecular diagnosis of autosomal recessive spastic ataxia of Charlevoix-Saguenay. Critical issues addressed during the development of the CSCE system included the position of a DNA variant relative to the primers and the CG-content of the amplicons. The validation was performed in two phases; a retrospective analysis of 32 samples containing 41 different known DNA variants and a prospective analysis of 20 samples of patients clinically suspected of having autosomal recessive spastic ataxia of Charlevoix-Saguenay. These 20 samples appeared to contain 73 DNA variants. In total, in 32 out of the 45 amplicons, a DNA variant was present, which allowed verification of the detection capacity during the validation process. After optimization of the original design, the overall analytical sensitivity of CSCE for the SACS gene was 100%, and the analytical specificity of CSCE was 99.8%. In conclusion, CSCE is a robust technique with a high analytical sensitivity and specificity, and it can readily be used for mutation scanning of the large SACS gene. Furthermore this technique is less time-consuming and less expensive, as compared with standard automated sequencing.

Long-term follow-up of patients with recessive dystrophic epidermolysis bullosa in the Netherlands: expansion of the mutation database and unusual phenotype-genotype correlations.

BACKGROUND: The current classification of recessive dystrophic epidermolysis bullosa (RDEB) comprises two major subtypes: 'severe generalized RDEB' (RDEB-sev gen) with early-onset, extensive, generalized blistering and scarring, complete absence of type VII collagen, and bi-allelic COL7A1 null mutations; milder 'generalized other RDEB' (RDEB-O) with reduced-to-normal type VII collagen expression, and non-null genotypes. OBJECTIVE: To search for previously unrecognized phenotype-genotype correlations in 33 Dutch RDEB families. METHODS: We analyzed extensive clinical follow-up data, available for all patients up to 19 years, detailed type VII collagen immunostaining and genotypes, and correlated clinical phenotype to molecular phenotype and genotype. RESULTS: We identified 20 novel COL7A1 mutations. In 14 of 15 RDEB-sev gen patients type VII collagen was completely absent, one had strongly reduced type VII collagen, and all carried bi-allelic null mutations. Five of 11 RDEB-O patients developed pseudosyndactyly of the fingers preceded by skin atrophy and flexion contractures later in childhood and adolescence. All five had esophageal involvement and growth retardation. Type VII collagen immunostaining ranged from strongly reduced to slightly reduced in RDEB-O patients with pseudosyndactyly, whereas RDEB-O patients without pseudosyndactyly had slightly reduced to normal type VII collagen staining. There was no difference in genotypes between both groups, although we unexpectedly found bi-allelic null mutations in two of five RDEB-O patients with pseudosyndactyly. CONCLUSION: Pseudosyndactyly occurs in approximately half of RDEB-O patients when type VII collagen is strongly reduced. The prognosis in RDEB cannot always be simply predicted from the COL7A1 genotype.

Spectrum of p63 mutations in a selected patient cohort affected with ankyloblepharon-ectodermal defects-cleft lip/palate syndrome (AEC).

Heterozygous mutations in the p63 gene underlie a group of at least seven allelic syndromes, including ankyloblepharon-ectodermal defects-cleft lip/palate syndrome (AEC) and Rapp Hodgkin syndrome (RHS), which involves varying degrees of ectodermal dysplasia, orofacial clefting and limb malformations. Mutations in the AEC and Rapp Hodgkin syndromes cluster in the 3' end of the p63 gene. Previously reported mutations are mainly missense and frameshift mutations in exons 13 and 14, affecting the p63alpha-specific SAM (sterile alpha motif) and TI (transactivation inhibitory) domains. A patient cohort affected by AEC syndrome was evaluated during International Research Symposium supported by the National Foundation for Ectodermal Dysplasias. Nineteen patients underwent full clinical evaluations and 18 had findings consistent with a diagnosis of AEC syndrome. These 19 patients, along with 5 additional relatives had genomic DNA analysis. Twenty-one of the 24 participants from 12 families were found to have mutations in the p63 gene. Eleven different mutations were identified; 10 were novel mutations. Eight were missense mutations within the coding region of the SAM domain. Three other mutations were located in exon 14 sequences, which encode the TI domain. The effects of the mutations in the SAM and TI domains are poorly understood and functional studies are required to understand the pathological mechanisms. However, AEC and RHS mutations in the 5' and 3' ends of the p63 gene point towards a critical role of the DeltaNp63alpha isoform for the AEC/RHS phenotype.