Unravelling undiagnosed rare disease cases by HiFi long-read genome sequencing.

Solve-RD is a pan-European rare disease (RD) research program that aims to identify disease-causing genetic variants in previously undiagnosed RD families. We utilised 10-fold coverage HiFi long-read sequencing (LRS) for detecting causative structural variants (SVs), single nucleotide variants (SNVs), insertion-deletions (InDels), and short tandem repeat (STR) expansions in extensively studied RD families without clear molecular diagnoses. Our cohort includes 293 individuals from 114 genetically undiagnosed RD families selected by European Rare Disease Network (ERN) experts. Of these, 21 families were affected by so-called 'unsolvable' syndromes for which genetic causes remain unknown, and 93 families with at least one individual affected by a rare neurological, neuromuscular, or epilepsy disorder without genetic diagnosis despite extensive prior testing. Clinical interpretation and orthogonal validation of variants in known disease genes yielded thirteen novel genetic diagnoses due to de novo and rare inherited SNVs, InDels, SVs, and STR expansions. In an additional four families, we identified a candidate disease-causing SV affecting several genes including an MCF2 / FGF13 fusion and PSMA3 deletion. However, no common genetic cause was identified in any of the 'unsolvable' syndromes. Taken together, we found (likely) disease-causing genetic variants in 13.0% of previously unsolved families and additional candidate disease-causing SVs in another 4.3% of these families. In conclusion, our results demonstrate the added value of HiFi long-read genome sequencing in undiagnosed rare diseases.

A recurrent ACTA1 amino acid change in mosaic form causes milder asymmetric myopathy.

We describe three patients with asymmetric congenital myopathy without definite nemaline bodies and one patient with severe nemaline myopathy. In all four patients, the phenotype had been caused by pathogenic missense variants in ACTA1 leading to the same amino acid change, p.(Gly247Arg). The three patients with milder myopathy were mosaic for their variants. In contrast, in the severely affected patient, the missense variant was present in a de novo, constitutional form. The grade of mosaicism in the three mosaic patients ranged between 20 % and 40 %. We speculate that the milder clinical and histological manifestations of the same ACTA1 variant in the patients with mosaicism reflect the lower abundance of mutant actin in their muscle tissue. Similarly, the asymmetry of body growth and muscle weakness may be a consequence of the affected cells being unevenly distributed. The partial improvement in muscle strength with age in patients with mosaicism might be due to an increased proportion over time of nuclei carrying and expressing two normal alleles.

Novel Compound Heterozygous Splice-Site Variants in TPM3 Revealed by RNA Sequencing in a Patient with an Unusual Form of Nemaline Myopathy: A Case Report.

BACKGROUND: Pathogenic variants in the TPM3 gene, encoding slow skeletal muscle α-tropomyosin account for less than 5% of nemaline myopathy cases. Dominantly inherited or de novo missense variants in TPM3 are more common than recessive loss-of-function variants. The recessive variants reported to date seem to affect either the 5' or the 3' end of the skeletal muscle-specific TPM3 transcript. OBJECTIVES: The aim of the study was to identify the disease-causing gene and variants in a Finnish patient with an unusual form of nemaline myopathy. METHODS: The genetic analyses included Sanger sequencing, whole-exome sequencing, targeted array-CGH, and linked-read whole genome sequencing. RNA sequencing was done on total RNA extracted from cultured myoblasts and myotubes of the patient and controls. TPM3 protein expression was assessed by Western blot analysis. The diagnostic muscle biopsy was analyzed by routine histopathological methods. RESULTS: The patient had poor head control and failure to thrive, but no hypomimia, and his upper limbs were clearly weaker than his lower limbs, features which in combination with the histopathology suggested TPM3-caused nemaline myopathy. Muscle histopathology showed increased fiber size variation and numerous nemaline bodies predominantly in small type 1 fibers. The patient was found to be compound heterozygous for two splice-site variants in intron 1a of TPM3: NM_152263.4:c.117+2_5delTAGG, deleting the donor splice site of intron 1a, and NM_152263.4:c.117 + 164 C>T, which activates an acceptor splice site preceding a non-coding exon in intron 1a. RNA sequencing revealed inclusion of intron 1a and the non-coding exon in the transcripts, resulting in early premature stop codons. Western blot using patient myoblasts revealed markedly reduced levels of the TPM3 protein. CONCLUSIONS: Novel biallelic splice-site variants were shown to markedly reduce TPM3 protein expression. The effects of the variants on splicing were readily revealed by RNA sequencing, demonstrating the power of the method.

Array Comparative Genomic Hybridisation and Droplet Digital PCR Uncover Recurrent Copy Number Variation of the TTN Segmental Duplication Region.

Intragenic segmental duplication regions are potential hotspots for recurrent copy number variation and possible pathogenic aberrations. Two large sarcomeric genes, nebulin and titin, both contain such segmental duplication regions. Using our custom Comparative Genomic Hybridisation array, we have previously shown that a gain or loss of more than one copy of the repeated block of the nebulin triplicate region constitutes a recessive pathogenic mutation. Using targeted array-CGH, similar copy number variants can be detected in the segmental duplication region of titin. Due to the limitations of the array-CGH methodology and the repetitiveness of the region, the exact copy numbers of the blocks could not be determined. Therefore, we developed complementary custom Droplet Digital PCR assays for the titin segmental duplication region to confirm true variation. Our combined methods show that the titin segmental duplication region is subject to recurrent copy number variation. Gains and losses were detected in samples from healthy individuals as well as in samples from patients with different muscle disorders. The copy number variation observed in our cohort is likely benign, but pathogenic copy number variants in the segmental duplication region of titin cannot be excluded. Further investigations are needed, however, this region should no longer be neglected in genetic analyses.

A custom ddPCR method for the detection of copy number variations in the nebulin triplicate region.

The human genome contains repetitive regions, such as segmental duplications, known to be prone to copy number variation. Segmental duplications are highly identical and homologous sequences, posing a specific challenge for most mutation detection methods. The giant nebulin gene is expressed in skeletal muscle. It harbors a large segmental duplication region composed of eight exons repeated three times, the so-called triplicate region. Mutations in nebulin are known to cause nemaline myopathy and other congenital myopathies. Using our custom targeted Comparative Genomic Hybridization arrays, we have previously shown that copy number variations in the nebulin triplicate region are pathogenic when the copy number of the segmental duplication block deviates two or more copies from the normal number, which is three per allele. To complement our Comparative Genomic Hybridization arrays, we have established a custom Droplet Digital PCR method for the detection of copy number variations within the nebulin triplicate region. The custom Droplet Digital PCR assays allow sensitive, rapid, high-throughput, and cost-effective detection of copy number variations within this region and is ready for implementation a screening method for disease-causing copy number variations of the nebulin triplicate region. We suggest that Droplet Digital PCR may also be used in the study and diagnostics of other segmental duplication regions of the genome.

Congenital asymmetric distal myopathy with hemifacial weakness caused by a heterozygous large de novo mosaic deletion in nebulin.

We report the first mosaic mutation, a deletion of exons 11-107, identified in the nebulin gene in a Finnish patient presenting with a predominantly distal congenital myopathy and asymmetric muscle weakness. The female patient is ambulant and currently 26 years old. Muscle biopsies showed myopathic features with type 1 fibre predominance, strikingly hypotrophic type 2 fibres and central nuclei, but no nemaline bodies. The deletion was detected in a copy number variation analysis based on next-generation sequencing data. The parents of the patient did not carry the deletion. Mosaicism was detected using a custom, targeted comparative genomic hybridisation array. Expression of the truncated allele, less than half the size of full-length nebulin, was confirmed by Western blotting. The clinical and histological picture resembled that of a family with a slightly smaller deletion, and that in patients with recessively inherited distal forms of nebulin-caused myopathy. Asymmetry, however, was a novel feature.

Improving Copy Number Variant Detection from Sequencing Data with a Combination of Programs and a Predictive Model.

Bioinformatics tools for analyzing copy number variants (CNVs) from massively parallel sequencing (MPS) data are less well developed compared with other variant types. We present an efficient bioinformatics pipeline for CNV detection from gene panel MPS data in neuromuscular disorders. CNVs were generated in silico into samples sequenced with a previously published MPS gene panel. The in silico CNVs from these samples were analyzed with four programs having complementary CNV detection ranges: CoNIFER, XHMM, ExomeDepth, and CODEX. A logistic regression model was trained with the obtained set of in silico CNV detections to predict true-positive CNV detections among all CNV detections from samples. This model was validated using 66 control samples with a verified true-positive (n = 58) or false-positive (n = 8) CNV detection. Applying all four programs together provided more sensitive detection results with in silico CNVs than other program combinations or any program alone. Furthermore, a model with CNV detection-specific scores from all four programs as variables performed overall best in the validation. No single program could detect all CNV sizes and types equally or with enough accuracy. Therefore, a combination of carefully selected programs should be used to maximize detection accuracy. In addition, the detected CNVs should be reviewed with a statistical model to streamline and standardize the filtering of the detections for annotation.

Dominantly inherited distal nemaline/cap myopathy caused by a large deletion in the nebulin gene.

We report the first family with a dominantly inherited mutation of the nebulin gene (NEB). This ∼100 kb in-frame deletion encompasses NEB exons 14-89, causing distal nemaline/cap myopathy in a three-generation family. It is the largest deletion characterized in NEB hitherto. The mutated allele was shown to be expressed at the mRNA level and furthermore, for the first time, a deletion was shown to cause the production of a smaller mutant nebulin protein. Thus, we suggest that this novel mutant nebulin protein has a dominant-negative effect, explaining the first documented dominant inheritance of nebulin-caused myopathy. The index patient, a young man, was more severely affected than his mother and grandmother. His first symptom was foot drop at the age of three, followed by distal muscle atrophy, slight hypomimia, high-arched palate, and weakness of the neck and elbow flexors, hands, tibialis anterior and toe extensors. Muscle biopsies showed myopathic features with type 1 fibre predominance in the index patient and nemaline bodies and cap-like structures in biopsies from his mother and grandmother. The muscle biopsy findings constitute a further example of nemaline bodies and cap-like structures being part of the same spectrum of pathological changes.

A nebulin super-repeat panel reveals stronger actin binding toward the ends of the super-repeat region.

INTRODUCTION: Nebulin is a giant actin-binding protein in the thin filament of the skeletal muscle sarcomere. Studies of nebulin interactions are limited by the size, complexity, and poor solubility of the protein. We divided the nebulin super-repeat region into a super-repeat panel, and studied nebulin/actin interactions. METHODS: Actin binding was studied using a co-sedimentation assay with filamentous actin and 26 different nebulin super-repeats. RESULTS: The panel revealed notable differences in actin binding between the super-repeats. Both ends of the super-repeat region bound actin significantly more strongly, whereas the central part of the protein bound actin weakly. Thus, the binding between nebulin and actin formed a location-dependent pattern of strong vs. weak binding. DISCUSSION: The nebulin super-repeat panel allowed us to study the actin binding of each super-repeat individually. The panel will be a powerful tool in elucidating nebulin function in health and disease. Muscle Nerve 59:116-121, 2019.

New mutations found by Next-Generation Sequencing screening of Spanish patients with Nemaline Myopathy.

Nemaline Myopathy (NM) is a rare genetic disorder that encompasses a large spectrum of myopathies characterized by hypotonia and generalized muscle weakness. To date, mutations in thirteen different genes have been associated with NM. The most frequently responsible genes are NEB (50% of cases) and ACTA1 (15-25% of cases). In this report all known NM related genes were screened by Next Generation Sequencing in five Spanish patients in order to genetically confirm the clinical and histological diagnosis of NM. Four mutations in NEB (c.17779_17780delTA, c.11086A>C, c.21076C>T and c.2310+5G>A) and one mutation in ACTA1 (c.871A>T) were found in four patients. Three of the four mutations in NEB were novel. A cDNA sequencing assay of the novel variants c.17779_17780delTA, c.11086A>C and c.2310+5G>A revealed that the intronic variant c.2310+5G>A affected the splicing process. Mutations reported here could help clinicians and geneticists in NM diagnosis.

An Extended Targeted Copy Number Variation Detection Array Including 187 Genes for the Diagnostics of Neuromuscular Disorders.

BACKGROUND: Our previous array, the Comparative Genomic Hybridisation design (CGH-array) for nemaline myopathy (NM), named the NM-CGH array, revealed pathogenic copy number variation (CNV) in the genes for nebulin (NEB) and tropomyosin 3 (TPM3), as well as recurrent CNVs in the segmental duplication (SD), i.e. triplicate, region of NEB (TRI, exons 82-89, 90-97, 98-105). In the light of this knowledge, we have designed and validated an extended CGH array, which includes a selection of 187 genes known to cause neuromuscular disorders (NMDs). OBJECTIVE: Our aim was to develop a reliable method for CNV detection in genes related to neuromuscular disorders for routine mutation detection and analysis, as a much-needed complement to sequencing methods. METHODS: We have developed a novel custom-made 4×180 k CGH array for the diagnostics of NMDs. It includes the same tiled ultra-high density coverage of the 12 known or putative NM genes as our 8×60 k NM-CGH-array but also comprises a selection of 175 additional genes associated with NMDs, including titin (TTN), at a high to very high coverage. The genes were divided into three coverage groups according to known and potential pathogenicity in neuromuscular disorders. RESULTS: The array detected known and putative CNVs in all three gene coverage groups, including the repetitive regions of NEB and TTN. CONCLUSIONS: The targeted neuromuscular disorder 4×180 k array-CGH (NMD-CGH-array v1.0) design allows CNV detection for a broader spectrum of neuromuscular disorders at a high resolution.

Copy number variation analysis increases the diagnostic yield in muscle diseases.

OBJECTIVE: Copy number variants (CNVs) were analyzed from next-generation sequencing data, with the aim of improving diagnostic yield in skeletal muscle disorder cases. METHODS: Four publicly available bioinformatic analytic tools were used to analyze CNVs from sequencing data from patients with muscle diseases. The patients were previously analyzed with a targeted gene panel for single nucleotide variants and small insertions and deletions, without achieving final diagnosis. Variants detected by multiple CNV analysis tools were verified with either array comparative genomic hybridization or PCR. The clinical significance of the verified CNVs was interpreted, considering previously identified variants, segregation studies, and clinical information of the patient cases. RESULTS: Combining analysis of all different mutation types enabled integration of results and identified the final cause of the disease in 9 myopathy cases. Complex effects like compound heterozygosity of different mutation types and compound disease arising from variants of different genes were unraveled. We identified the first large intragenic deletion of the titin (TTN) gene implicated in the pathogenesis of a severe form of myopathy. Our work also revealed a "double-trouble" effect in a patient carrying a single heterozygous insertion/deletion mutation in the TTN gene and a Becker muscular dystrophy causing deletion in the dystrophin gene. CONCLUSIONS: Causative CNVs were identified proving that analysis of CNVs is essential for increasing the diagnostic yield in muscle diseases. Complex severe muscular dystrophy phenotypes can be the result of different mutation types but also of the compound effect of 2 different genetic diseases.