Evolution of genome diagnostics in epidermolysis bullosa: Unveiling the power of next-generation sequencing.

BACKGROUND: Genome diagnostics is considered gold standard diagnostics for epidermolysis bullosa (EB), a phenotypically and genetically heterogeneous group of rare disorders characterized by blistering and wounding of mucocutaneous tissues. EB is caused by pathogenic variants in genes encoding proteins of the dermo-epidermal junction. Accurate genetic diagnosis of EB is crucial for prognostication, counselling and precision-medicine. Genome diagnostics for EB started in 1991 with the introduction of Sanger sequencing (SS), analysing one gene at a time. In 2013, SS was superseded by next-generation sequencing (NGS), that allow for high-throughput sequencing of multiple genes in parallel. Several studies have shown a beneficial role for NGS in EB diagnostics, but its true benefit has not been quantified. OBJECTIVES: To determine the benefit of NGS in EB by systematically evaluating the performance of different genome diagnostics used over time based on robust data from the Dutch EB Registry. METHODS: The diagnostic performances of SS and NGS were systematically evaluated in a retrospective observational study including all index cases with a clinical diagnosis of EB in whom genome diagnostics was performed between 01 January 1994 and 01 January 2022 (n = 308), registered at the Dutch EB Expertise Centre. RESULTS: Over time, a genetic diagnosis was made in 289/308 (94%) EB cases. The diagnostic yield increased from 89% (SS) to 95% (NGS). Most importantly, NGS significantly reduced diagnostic turnaround time (39 days vs. 211 days, p < 0.001). The likelihood of detecting variants of uncertain significance and additional findings increased from 5% and 1% (SS) to 22% and 13% (NGS) respectively. CONCLUSIONS: Our study quantifies the benefit of NGS-based methods and demonstrate they have had a major impact on EB diagnostics through an increased diagnostic yield and a dramatically decreased turnaround time (39 days). Although our diagnostic yield is high (95%), further improvement of genome diagnostics is urgently needed to provide a genetic diagnosis in all EB patients.

Junctional epidermolysis bullosa of late onset explained by mutations in COL17A1.

BACKGROUND: Junctional epidermolysis bullosa of late onset (JEB-lo) is a rare disease characterized by blistering of primarily the hands and feet starting in childhood. The pathogenesis remains unclear. OBJECTIVES: To clarify the pathogenesis of JEB-lo. METHODS: Two patients with JEB-lo, a brother and a sister, were examined using electron microscopy (EM), immunofluorescence (IF) antigen mapping and molecular analysis. RESULTS: We found subtle changes in IF antigen mapping and EM. The most remarkable changes were loss of the apical-lateral staining of monoclonal antibodies (mAbs) against type XVII collagen (Col17), and a broadened distribution of mAb staining against the ectodomain of Col17, laminin-332 and type VII collagen. Mutation analysis of COL17A1, encoding Col17, showed a compound heterozygosity for a novel mutation c.1992_1995delGGGT and the known mutation c.3908G>A in both patients. The deletion c.1992_1995delGGGT results in a premature termination codon and mRNA decay, leaving the patients functionally hemizygous for the missense mutation c.3908G>A (p.R1303Q) in the noncollagenous 4 domain of Col17. CONCLUSIONS: JEB-lo is an autosomal recessive disorder caused by mutations in COL17A1, and subtle aberrations in EM and IF antigen mapping are clues to diagnosis.

Herlitz junctional epidermolysis bullosa: diagnostic features, mutational profile, incidence and population carrier frequency in the Netherlands.

BACKGROUND: Junctional epidermolysis bullosa, type Herlitz (JEB-H) is a lethal, autosomal recessive blistering disease caused by null mutations in the genes coding for the lamina lucida/densa adhesion protein laminin-332 (LAMB3, LAMA3 and LAMC2). OBJECTIVES: To present the diagnostic features and molecular analyses of all 22 patients with JEB-H in the Dutch Epidermolysis Bullosa Registry between 1988 and 2011, and to calculate the disease incidence and carrier frequency in the Netherlands. METHODS: All patients were analysed with immunofluorescence antigen mapping (IF), electron microscopy (EM) and molecular analysis. RESULTS: The mean lifespan of our patients with JEB-H was 5·8 months (range 0·5-32·6). IF showed absent (91%) or strongly reduced (9%) staining for laminin-332 with monoclonal antibody GB3. In EM the hemidesmosomes and sub-basal dense plates were hypoplastic or absent. We identified mutations in all 22 patients: in 19 we found LAMB3 mutations, in two LAMA3 mutations, and in one LAMC2 mutations. We found three novel splice site mutations in LAMB3: (i) c.29-2A>G resulting in an out-of-frame skip of exon 3 and a premature termination codon (PTC); (ii) c.1289-2_1296del10 leading to an out-of-frame skip of exon 12 and a PTC; and (iii) c.3228+1G>T leading to an exon 21 skip. CONCLUSIONS: All diagnostic tools should be evaluated to clarify the diagnosis of JEB-H. We have identified 11 different mutations in 22 patients with JEB-H, three of them novel. In the Netherlands the incidence rate of JEB-H is 4·0 per one million live births. The carrier frequency of a JEB-H mutation in the Dutch population is 1 in 249.

The human E48 antigen, highly homologous to the murine Ly-6 antigen ThB, is a GPI-anchored molecule apparently involved in keratinocyte cell-cell adhesion.

The E48 antigen, a putative human homologue of the 20-kD protein present in desmosomal preparations of bovine muzzle, and formerly called desmoglein III (dg4), is a promising target antigen for antibody-based therapy of squamous cell carcinoma in man. To anticipate the effect of high antibody dose treatment, and to evaluate the possible biological involvement of the antigen in carcinogenesis, we set out to molecularly characterize the antigen. A cDNA clone encoding the E48 antigen was isolated by expression cloning in COS cells. Sequence analysis revealed that the clone contained an open reading frame of 128 amino acids, encoding a core protein of 13,286 kD. Database searching showed that the E48 antigen has a high level of sequence similarity with the mouse ThB antigen, a member of the Ly-6 antigen family. Phosphatidylinositol-specific (PI-specific) phospholipase-C treatment indicated that the E48 antigen is glycosylphosphatidylinositol-anchored (GPI-anchored) to the plasma membrane. The gene encoding the E48 antigen is a single copy gene, located on human chromosome 8 in the 8q24-qter region. The expression of the gene is confined to keratinocytes and squamous tumor cells. The putative mouse homologue, the ThB antigen, originally identified as an antigen on cells of the lymphocyte lineage, was shown to be highly expressed in squamous mouse epithelia. Moreover, the ThB expression level is in keratinocytes, in contrast to that in lymphocytes, not mouse strain related. Transfection of mouse SV40-polyoma transformed mouse NIH/3T3 cells with the E48 cDNA confirmed that the antigen is likely to be involved in cell-cell adhesion.

An association screen of myelin-related genes implicates the chromosome 22q11 PIK4CA gene in schizophrenia.

Several lines of evidence, including expression analyses, brain imaging and genetic studies suggest that the integrity of myelin is disturbed in schizophrenia patients. In this study, we first reconstructed a pathway of 138 myelin-related genes, all involved in myelin structure, composition, development or maintenance. Then we performed a two-stage association analysis on these 138 genes using 771 single nucleotide polymorphisms (SNPs). Analysis of our data from 310 cases vs 880 controls demonstrated association of 10 SNPs from six genes. Specifically, we observed highly significant P-values for association in PIK4CA (observed P=6.1 x 10(-6)). These findings remained significant after Bonferroni correction for 771 tests. The PIK4CA gene is located in the chromosome 22q11 deletion syndrome region, which is of particular interest because it has been implicated in schizophrenia. We also report weak association of SNPs in PIK3C2G, FGF1, FGFR1, ARHGEF10 and PSAP (observed P

An amino-terminal DAX1 (NROB1) missense mutation associated with isolated mineralocorticoid deficiency.

CONTEXT: Mutations in DAX1 (dosage-sensitive sex reversal-adrenal hypoplasia congenita critical region on the X chromosome gene 1; NR0B1) cause X-linked adrenal hypoplasia congenita, a disease characterized by primary adrenal failure, testicular dysgenesis, and gonadotropin deficiency. Most DAX1 mutations are deletions, nonsense, or frameshift mutations that markedly impair its transcriptional activity. Missense mutations have been restricted to the carboxy-terminal domain and are associated with more variable clinical phenotypes. OBJECTIVE: The objective was to identify novel clinical phenotypes associated with DAX1 missense mutations. PATIENTS AND DESIGN: We investigated the genetic basis of isolated mineralocorticoid deficiency in a patient who carries a unique missense mutation (W105C) in the amino-terminal region of DAX1. RESULTS: The W105C DAX1 mutation in the proband was present in three asymptomatic hemizygous males, but it was not detected in the general population. Using in vitro studies of DAX1 expression and function in transfected cells, we demonstrate that the mutant DAX1 protein exhibits mild loss of function, whether studied for genes it represses or for genes it activates. Structure-function studies suggest that the W105C and other mutations in the aminoterminus are compensated by the presence of repeated LXXLL motifs that mediate DAX1 interactions with other proteins. CONCLUSIONS: We describe the first missense mutation in the aminoterminus of DAX1 and conclude that mutations in this region may be partially compensated by redundant functional domains. Mild DAX1 mutations may be a cause of isolated mineralocorticoid deficiency.

The PIP5K2A and RGS4 genes are differentially associated with deficit and non-deficit schizophrenia.

Several putative schizophrenia susceptibility genes have recently been reported, but it is not clear whether these genes are associated with schizophrenia in general or with specific disease subtypes. In a previous study, we found an association of the neuregulin 1 (NRG1) gene with non-deficit schizophrenia only. We now report an association study of four schizophrenia candidate genes in patients with and without deficit schizophrenia, which is characterized by severe and enduring negative symptoms. Single-nucleotide polymorphisms (SNPs) were genotyped in the DTNBP1 (dysbindin), G72/G30 and RGS4 genes, and the relatively unknown PIP5K2A gene, which is located in a region of linkage with both schizophrenia and bipolar disorder. The sample consisted of 273 Dutch schizophrenia patients, 146 of whom were diagnosed with deficit schizophrenia and 580 controls. The strongest evidence for association was found for the A-allele of SNP rs10828317 in the PIP5K2A gene, which was associated with both clinical subtypes (P = 0.0004 in the entire group; non-deficit P = 0.016, deficit P = 0.002). Interestingly, this SNP leads to a change in protein composition. In RGS4, the G-allele of the previously reported SNP RGS4-1 (single and as part of haplotypes with SNP RGS4-18) was associated with non-deficit schizophrenia (P = 0.03) but not with deficit schizophrenia (P = 0.79). SNPs in the DTNBP1 and G72/G30 genes were not significantly associated in any group. In conclusion, our data provide further evidence that specific genes may be involved in different schizophrenia subtypes and suggest that the PIP5K2A gene deserves further study as a general susceptibility gene for schizophrenia.

Novel Algorithms for Improved Sensitivity in Non-Invasive Prenatal Testing.

Non-invasive prenatal testing (NIPT) of cell-free DNA in maternal plasma, which is a mixture of maternal DNA and a low percentage of fetal DNA, can detect fetal aneuploidies using massively parallel sequencing. Because of the low percentage of fetal DNA, methods with high sensitivity and precision are required. However, sequencing variation lowers sensitivity and hampers detection of trisomy samples. Therefore, we have developed three algorithms to improve sensitivity and specificity: the chi-squared-based variation reduction (χ2VR), the regression-based Z-score (RBZ) and the Match QC score. The χ2VR reduces variability in sequence read counts per chromosome between samples, the RBZ allows for more precise trisomy prediction, and the Match QC score shows if the control group used is representative for a specific sample. We compared the performance of χ2VR to that of existing variation reduction algorithms (peak and GC correction) and that of RBZ to trisomy prediction algorithms (standard Z-score, normalized chromosome value and median-absolute-deviation-based Z-score). χ2VR and the RBZ both reduce variability more than existing methods, and thereby increase the sensitivity of the NIPT analysis. We found the optimal combination of algorithms was to use both GC correction and χ2VR for pre-processing and to use RBZ as the trisomy prediction method.

Sublocalization of the synovial sarcoma-associated t(X;18) chromosomal breakpoint in Xp11.2 using cosmid cloning and fluorescence in situ hybridization.

In a previous study we localized the synovial sarcoma-associated t(X;18)(p11;q11) breakpoint within the ornithine aminotransferase-like 1 (OATL1) cluster on the X chromosome. This localization was delineated from both somatic cell hybrid and fluorescence in situ hybridization (FISH) analysis of patient material, using OAT-specific cDNA and YAC probes. Simultaneously, Knight et al. (1992, Mol. Hum. Genet, in press) mapped this same breakpoint in their patient material adjacent to the more proximal OATL2 region on the X chromosome. Here we report the analysis of two additional tumors and demonstrate that again in these cases the chromosomal break occurs within the OATL1 cluster. In order to further specify the breakpoint, we subcloned the OATL1 YAC (no. 2) into cosmids. At least one of these cosmids (0.38) hybridizes to sequences that bracket the translocation breakpoint, as demonstrated by both Southern blot and FISH analysis. These observations confirm and substantiate our previous findings. In addition, cosmid 0.38 should be a valuable instrument for the ultimate isolation and identification of the gene(s) involved in the development of synovial sarcoma.

Mapping of the SCA23 locus involved in autosomal dominant cerebellar ataxia to chromosome region 20p13-12.3.

We report upon a Dutch autosomal dominant cerebellar ataxia (ADCA) family, clinically characterized by a late-onset (>40 years), slowly progressive, isolated spinocerebellar ataxia (SCA). Neuropathological examination in one affected subject showed neuronal loss in the Purkinje cell layer, dentate nuclei and inferior olives, thinning of cerebellopontine tracts, demyelination of posterior and lateral columns in the spinal cord, as well as ubiquitin-positive intranuclear inclusions in nigral neurons that were considered to be Marinesco bodies. Data obtained from the genome-wide linkage analysis revealed a maximal lod score of 3.46 at = 0.00 for marker D20S199. This new SCA locus, on chromosome region 20p13-p12.3, was designated SCA23 after approval by the HUGO Nomenclature Committee. Currently, candidate genes are being screened for mutations within the SCA23 interval. In addition to the recently identified SCA14, SCA19 and FGF14 families, SCA23 is yet another novel SCA locus in the Dutch ADCA population, which further defines the genetic heterogeneity of ADCA families in the Netherlands.

Clinical and molecular correlations in spinocerebellar ataxia type 6: a study of 24 Dutch families.

BACKGROUND: Autosomal dominant cerebellar ataxias (ADCAs), or spinocerebellar ataxias (SCAs), are a heterogeneous group of neurodegenerative disorders. Mild CAG repeat expansions in the alpha(1A) voltage-dependent calcium channel gene are associated with SCA type 6 (SCA6). OBJECTIVE: To obtain further insight into the contribution of SCA6 mutations to the phenotypic variability in Dutch patients with ataxia. DESIGN: Survey and case series. SETTING: Hospitalized care, referral center. PATIENTS AND METHODS: The SCA6 locus was analyzed for CAG repeat expansions in a referred sample of 220 Dutch families with progressive cerebellar ataxia. Clinical characteristics of patients with SCA6 were investigated and correlated with molecular findings. RESULTS: The diagnosis SCA6 was confirmed in 24 families comprising 30 familial and 4 sporadic cases. Mean +/- SD age at onset was 50.1 +/- 11.1 years. Expanded CAG repeats with sizes 22, 23, and 25 were found. These sizes correlated inversely with age at onset. No intergenerational changes in CAG repeat size were detected. Despite this, 2 families showed clinical anticipation. CONCLUSIONS: This study provides the first detailed description of Dutch patients with SCA6. Clinical analysis identifies SCA6 as a late-onset ataxia in which eye movement abnormalities are prominent and consistent early manifestations. No single clinical sign can be considered specific for SCA6. Some patients have ataxia combined with episodic headaches or nausea, suggesting an overlap among SCA6, eposidic ataxia type 2, and familial hemiplegic migraine. Spinocerebellar ataxia type 6 accounts for approximately 11% of all Dutch families with ADCA. Analysis of SCA6 contributes further to the genetic classification of patients with ADCA, including patients without a clear family history of the disease.

Spinocerebellar ataxias in the Netherlands: prevalence and age at onset variance analysis.

BACKGROUND: International prevalence estimates of autosomal dominant cerebellar ataxias (ADCA) vary from 0.3 to 2.0 per 100,000. The prevalence of ADCA in the Netherlands is unknown. Fifteen genetic loci have been identified (SCA-1-8, SCA-10-14, SCA-16, and SCA-17) and nine of the corresponding genes have been cloned. In SCA-1, SCA2, SCA3, SCA6, SCA7, SCA-12 and SCA-17 the mutation has been shown to be an expanded CAG repeat. Previously, the length of the CAG repeat was found to account for 50 to 80% of variance in age at onset. Because of heterogeneity in encoded proteins, different pathophysiologic mechanisms leading to neurodegeneration could be involved. The relationship between CAG repeat length and age at onset would then differ accordingly. METHOD: Based on the results of SCA mutation analysis in the three DNA diagnostic laboratories that serve the entire Dutch population, the authors surveyed the number of families and affected individuals per SCA gene, as well as individual repeat length and age at onset. Regression analysis was applied to study the relationship between CAG repeat length and age at onset per SCA gene. The slopes of the different regression curves were compared. RESULTS: On November 1, 2000, mutations were found in 145 ADCA families and 391 affected individuals were identified. The authors extrapolated a minimal prevalence of 3.0 per 100,000 (range 2.8 to 3.8/100,000). SCA3 was the most frequent mutation. CAG repeat length contributed to 52 to 76% of age at onset variance. Regression curve slopes for SCA-1, SCA2, SCA3, and SCA7 did not differ significantly. CONCLUSIONS: The estimated minimal prevalence of ADCA in the Netherlands is 3.0 per 100,000 inhabitants. Except for SCA6, the relationship between age at onset and CAG repeat expansion does not differ significantly between SCA-1, SCA2, SCA3, and SCA7 patient groups in our population, indicating that these SCA subtypes share similar mechanisms of polyglutamine-induced neurotoxicity, despite heterogeneity in gene products.

Towards the isolation of a human malignant extragonadal germ cell tumour-associated breakpoint in chromosome 11q13.

In a previous study we have defined a subgroup of human malignant extragonadal germ cell tumours that is characterized by complex translocations involving chromosomes 6 and 11 (Echten et al. 1995). Here we report (i) the use of fluorescent in situ hybridization, pulsed field gel electrophoresis and direct visual hybridization techniques to localize the tumour-associated breakpoint within band 11q13, and (ii) the construction of a phage library enriched for this region to facilitate genomic walks towards the breakpoint. Extensive breakpoint-flanking contigs were generated and within these contigs six candidate genes could be identified.

Amplification of chromosome subregion 12p11.2-p12.1 in a metastasis of an i(12p)-negative seminoma: relationship to tumor progression?

Cytogenetic analysis of a metastasis of a human testicular germ cell tumor (seminoma) revealed multiple numerical and structural anomalies, including an abnormally banding region (ABR) present on the short arm of one of the chromosome 12 homologs. Fluorescence in situ- and comparative genomic hybridization experiments revealed that the ABR results from the amplification of 12p11.2-p12.1 derived sequences. We speculate that this particular region may harbor gene(s) relevant for testicular germ cell tumor progression.

Overrepresentation of chromosome 12p sequences and karyotypic evolution in i(12p)-negative testicular germ-cell tumors revealed by fluorescence in situ hybridization.

Human testicular germ-cell tumors (TGCTs) comprise a heterogeneous group of solid neoplasms. These tumors are characterized by the presence of a highly specific chromosomal abnormality, i.e., an isochromosome of the short arm of chromosome 12. At present, this i(12p) chromosome is found in more than 80% of TGCTs. Isochromosome 12p has also been observed in some ovarian and extragonadal germ cell tumors. In the remaining so-called i(12p)-negative TGCTs other abnormalities involving chromosome 12, mainly 12p, can be found. In order to establish whether 12p abnormalities other than i(12p) are a common phenomenon in TGCTs, a panel of 11 i(12p)-negative tumors was investigated using multicolor fluorescence in situ hybridization. All TGCTs examined appeared to contain chromosomal abnormalities involving 12p, resulting in a distinct overrepresentation of short arm sequences. In addition, indications were obtained for a clonal evolution in one of the tumors. Our data suggest that the occurrence of 12p abnormalities is a common phenomenon in i(12p)-negative TGCTs and that these abnormalities, analogous to i(12p), may contribute to the process of tumor development.

Identification of a yeast artificial chromosome that spans the human papillary renal cell carcinoma-associated t(X;1) breakpoint in Xp11.2.

Recently, a specific chromosome abnormality, t(X;1)(p11;q21), was described for a subgroup of human papillary renal cell carcinomas. The translocation breakpoint in Xp11 is located in the same region as that in t(X;18)(p11;q11)-positive synovial sarcoma. We used fluorescence in situ hybridization (FISH) and somatic cell hybridization techniques to demonstrate 1) that the Xp11 translocation breakpoint in papillary renal cell carcinoma differs from that observed in synovial sarcoma and has a more proximal location, and 2) that an ornithine aminotransferase (OAT)L2 containing yeast artificial chromosome (YAC) spans the X;1 translocation. This YAC provides an ideal starting point from which the breakpoint itself and the gene(s) involved can be isolated and characterized.

Fine mapping of the human renal oncocytoma-associated translocation (5;11)(q35;q13) breakpoint.

Recent cytogenetic analysis of a series of human renal oncocytomas revealed the presence of a recurring chromosomal translocation (5;11)(q35;q13) as sole anomaly in a subset of the tumors. The molecular characterization of this translocation was initiated using two primary t(5;11)-positive renal oncocytomas and a panel of somatic cell hybrids derived from one of these tumors, in conjunction with fluorescence in situ hybridization (FISH) and Southern blot analysis. The breakpoint in chromosome band 11q13 could be located within a genomic interval of at maximum 400 Kb immediately centromeric to the BCL1 locus.

Molecular characterization of a recurring complex chromosomal translocation in two human extragonadal germ cell tumors.

The molecular characterization of a recurring complex chromosomal translocation involving 6p21, 6p22, 6q23, and 11q13 in two independent but similar extragonadal human germ cell tumors was initiated using fluorescence in situ hybridization (FISH) and pulse field gel electrophoresis (PFGE) techniques. By using a series of specific probes from the 11q13 region, the translocation breakpoint in this chromosomal band could be located within a long-range restriction enzyme map in between the markers D11S457 and D11S546. In addition, aberrantly hybridizing restriction fragments were revealed by PFGE in both tumors, indicating that the breakpoint region must be located within a distance of at maximum 200 kilobase pairs (kbp) from the nearest DNA marker (D11S546).