Genome-wide identification of tandem repeats associated with splicing variation across 49 tissues in humans.

Tandem repeats (TRs) are one of the largest sources of polymorphism, and their length is associated with gene regulation. Although previous studies reported several tandem repeats regulating gene splicing in cis (spl-TRs), no large-scale study has been conducted. In this study, we established a genome-wide catalog of 9537 spl-TRs with a total of 58,290 significant TR-splicing associations across 49 tissues (false discovery rate 5%) by using Genotype-Tissue expression (GTex) Project data. Regression models explaining splicing variation by using spl-TRs and other flanking variants suggest that at least some of the spl-TRs directly modulate splicing. In our catalog, two spl-TRs are known loci for repeat expansion diseases, spinocerebellar ataxia 6 (SCA6) and 12 (SCA12). Splicing alterations by these spl-TRs were compatible with those observed in SCA6 and SCA12. Thus, our comprehensive spl-TR catalog may help elucidate the pathomechanism of genetic diseases.

Whole-exome sequencing reveals causative genetic variants for several overgrowth syndromes in molecularly negative Beckwith-Wiedemann spectrum.

Background Beckwith-Wiedemann syndrome (BWS) is an imprinting disorder caused by (epi)genetic alterations at 11p15. Because approximately 20% of patients test negative via molecular testing of peripheral blood leukocytes, the concept of Beckwith-Wiedemann spectrum (BWSp) was established to encompass a broader cohort with diverse and overlapping phenotypes. The prevalence of other overgrowth syndromes concealed within molecularly negative BWSp remains unexplored. Methods We conducted whole-exome sequencing (WES) on 69 singleton patients exhibiting molecularly negative BWSp. Variants were confirmed by Sanger sequencing or quantitative genomic PCR. We compared BWSp scores and clinical features between groups with classical BWS (cBWS), atypical BWS or isolated lateralised overgrowth (aBWS+ILO) and overgrowth syndromes identified via WES. Results Ten patients, one classified as aBWS and nine as cBWS, showed causative gene variants for Simpson-Golabi-Behmel syndrome (five patients), Sotos syndrome (two), Imagawa-Matsumoto syndrome (one), glycosylphosphatidylinositol biosynthesis defect 11 (one) or 8q duplication/9p deletion (one). BWSp scores did not distinguish between cBWS and other overgrowth syndromes. Birth weight and height in other overgrowth syndromes were significantly larger than in aBWS+ILO and cBWS, with varying intergroup frequencies of clinical features. Conclusion Molecularly negative BWSp encapsulates other syndromes, and considering both WES and clinical features may facilitate accurate diagnosis.

Complex chromosomal 6q rearrangements revealed by combined long-molecule genomics technologies.

Technologies for detecting structural variation (SV) have advanced with the advent of long-read sequencing, which enables the validation of SV at a nucleotide level. Optical genome mapping (OGM), a technology based on physical mapping, can also provide comprehensive SVs analysis. We applied long-read whole genome sequencing (LRWGS) to accurately reconstruct breakpoint (BP) segments in a patient with complex chromosome 6q rearrangements that remained elusive by conventional karyotyping. Although all BPs were precisely identified by LRWGS, there were two possible ways to construct the BP segments in terms of their orders and orientations. Thus, we also used OGM analysis. Notably, OGM recognized entire inversions exceeding 500 kb in size, which LRWGS could not characterize. Consequently, here we successfully unveil the full genomic structure of this complex chromosomal 6q rearrangement and cryptic SVs through combined long-molecule genomic analyses, showcasing how LRWGS and OGM can complement each other in SV analysis.

Long-term course of a case with a novel homozygous kyphoscoliosis peptidase variant.

We herein report a case with a novel homozygous variant in the kyphoscoliosis peptidase (KY) gene. A 58-year-old Japanese female was referred to our hospital with a gait disturbance that gradually worsened after the age of 50. She had bilateral equinus foot deformity since early childhood. Neurological examination revealed moderate weakness of the neck, trunk, femoral, and brachial muscles, mild respiratory failure, and areflexia. Whole-exome sequencing revealed a novel homozygous frameshift variant of the KY gene, NM_178554.6:c.824del p.(Glu275Glyfs*53). Our case demonstrated that KY-associated neuromuscular disease can present with extremely slow progressive muscle weakness and respiratory failure over a long natural course.

Detection of hidden intronic DDC variant in aromatic L-amino acid decarboxylase deficiency by adaptive sampling.

Aromatic l-amino acid decarboxylase (AADC) deficiency is an autosomal recessive neurotransmitter disorder caused by pathogenic DOPA decarboxylase (DDC) variants. We previously reported Japanese siblings with AADC deficiency, which was confirmed by the lack of enzyme activity; however, only a heterozygous missense variant was detected. We therefore performed targeted long-read sequencing by adaptive sampling to identify any missing variants. Haplotype phasing and variant calling identified a novel deep intronic variant (c.714+255 C > A), which was predicted to potentially activate the noncanonical splicing acceptor site. Minigene assay revealed that wild-type and c.714+255 C > A alleles had different impacts on splicing. Three transcripts, including the canonical transcript, were detected from the wild-type allele, but only the noncanonical cryptic exon was produced from the variant allele, indicating that c.714+255 C > A was pathogenic. Target long-read sequencing may be used to detect hidden pathogenic variants in unresolved autosomal recessive cases with only one disclosed hit variant.

A heterozygous germline deletion within USP8 causes severe neurodevelopmental delay with multiorgan abnormalities.

Ubiquitin-specific protease 8 (USP8) is a deubiquitinating enzyme involved in deubiquitinating the enhanced epidermal growth factor receptor for escape from degradation. Somatic variants at a hotspot in USP8 are a cause of Cushing's disease, and a de novo germline USP8 variant at this hotspot has been described only once previously, in a girl with Cushing's disease and developmental delay. In this study, we investigated an exome-negative patient with severe developmental delay, dysmorphic features, and multiorgan dysfunction by long-read sequencing, and identified a 22-kb de novo germline deletion within USP8 (chr15:50469966-50491995 [GRCh38]). The deletion involved the variant hotspot, one rhodanese domain, and two SH3 binding motifs, and was presumed to be generated through nonallelic homologous recombination through Alu elements. Thus, the patient may have perturbation of the endosomal sorting system and mitochondrial autophagy through the USP8 defect. This is the second reported case of a germline variant in USP8.

Novel missense variants cause intermediate phenotypes in the phenotypic spectrum of SLC5A6-related disorders.

SLC5A6 encodes the sodium-dependent multivitamin transporter, a transmembrane protein that uptakes biotin, pantothenic acid, and lipoic acid. Biallelic SLC5A6 variants cause sodium-dependent multivitamin transporter deficiency (SMVTD) and childhood-onset biotin-responsive peripheral motor neuropathy (COMNB), which both respond well to replacement therapy with the above three nutrients. SMVTD usually presents with various symptoms in multiple organs, such as gastrointestinal hemorrhage, brain atrophy, and global developmental delay, at birth or in infancy. Without nutrient replacement therapy, SMVTD can be lethal in early childhood. COMNB is clinically milder and has a later onset than SMVTD, at approximately 10 years of age. COMNB symptoms are mostly limited to peripheral motor neuropathy. Here we report three patients from one Japanese family harboring novel compound heterozygous missense variants in SLC5A6, namely NM_021095.4:c.[221C>T];[642G>C] p.[(Ser74Phe)];[(Gln214His)]. Both variants were predicted to be deleterious through multiple lines of evidence, including amino acid conservation, in silico predictions of pathogenicity, and protein structure considerations. Drosophila analysis also showed c.221C>T to be pathogenic. All three patients had congenital brain cysts on neonatal cranial imaging, but no other morphological abnormalities. They also had a mild motor developmental delay that almost completely resolved despite no treatment. In terms of severity, their phenotypes were intermediate between SMVTD and COMNB. From these findings we propose a new SLC5A6-related disorder, spontaneously remitting developmental delay with brain cysts (SRDDBC) whose phenotypic severity is between that of SMVTD and COMNB. Further clinical and genetic evidence is needed to support our suggestion.

Novel compound heterozygous ABCA2 variants cause IDPOGSA, a variable phenotypic syndrome with intellectual disability.

The gene for ATP binding cassette subfamily A member 2 (ABCA2) is located at chromosome 9q34.3. Biallelic ABCA2 variants lead to intellectual developmental disorder with poor growth and with or without seizures or ataxia (IDPOGSA). In this study, we identified novel compound heterozygous ABCA2 variants (NM_001606.5:c.[5300-17C>A];[6379C>T]) by whole exome sequencing in a 28-year-old Korean female patient with intellectual disability. These variants included intronic and nonsense variants of paternal and maternal origin, respectively, and are absent from gnomAD. SpliceAI predicted that the intron variant creates a cryptic acceptor site. Reverse transcription-PCR using RNA extracted from a lymphoblastoid cell line of the patient confirmed two aberrant transcripts. Her clinical features are compatible with those of IDPOGSA.

Complete nanopore repeat sequencing of SCA27B (GAA-FGF14 ataxia) in Japanese.

BACKGROUND: Although pure GAA expansion is considered pathogenic in SCA27B, non-GAA repeat motif is mostly mixed into longer repeat sequences. This study aimed to unravel the complete sequencing of FGF14 repeat expansion to elucidate its repeat motifs and pathogenicity. METHODS: We screened FGF14 repeat expansion in a Japanese cohort of 460 molecularly undiagnosed adult-onset cerebellar ataxia patients and 1022 controls, together with 92 non-Japanese controls, and performed nanopore sequencing of FGF14 repeat expansion. RESULTS: In the Japanese population, the GCA motif was predominantly observed as the non-GAA motif, whereas the GGA motif was frequently detected in non-Japanese controls. The 5'-common flanking variant was observed in all Japanese GAA repeat alleles within normal length, demonstrating its meiotic stability against repeat expansion. In both patients and controls, pure GAA repeat was up to 400 units in length, whereas non-pathogenic GAA-GCA repeat was larger, up to 900 units, but they evolved from different haplotypes, as rs534066520, located just upstream of the repeat sequence, completely discriminated them. Both (GAA)≥250 and (GAA)≥200 were enriched in patients, whereas (GAA-GCA)≥200 was similarly observed in patients and controls, suggesting the pathogenic threshold of (GAA)≥200 for cerebellar ataxia. We identified 14 patients with SCA27B (3.0%), but their single-nucleotide polymorphism genotype indicated different founder alleles between Japanese and Caucasians. The low prevalence of SCA27B in Japanese may be due to the lower allele frequency of (GAA)≥250 in the Japanese population than in Caucasians (0.15% vs 0.32%-1.26%). CONCLUSIONS: FGF14 repeat expansion has unique features of pathogenicity and allelic origin, as revealed by a single ethnic study.

Brain mosaicism of hedgehog signalling and other cilia genes in hypothalamic hamartoma.

Hypothalamic hamartoma (HH) is a rare benign developmental brain lesion commonly associated with a well characterized epilepsy phenotype. Most individuals with HH are non-syndromic without additional developmental anomalies nor a family history of disease. Nonetheless, HH is a feature of Pallister-Hall (PHS) and Oro-Facial-Digital Type VI (OFD VI) syndromes, both characterized by additional developmental anomalies. Initial genetic of analysis HH began with syndromic HH, where germline inherited or de novo variants in GLI3, encoding a central transcription factor in the sonic hedgehog (Shh) signalling pathway, were identified in most individuals with PHS. Following these discoveries in syndromic HH, the hypothesis that post-zygotic mosaicism in related genes may underly non-syndromic HH was tested. We discuss the identified mosaic variants within individuals with non-syndromic HH, review the analytical methodologies and diagnostic yields, and explore understanding of the functional role of the implicated genes with respect to Shh signalling, and cilia development and function. We also outline future challenges in studying non-syndromic HH and suggest potential novel strategies to interrogate brain mosaicism in HH.

De novo ARF3 variants cause neurodevelopmental disorder with brain abnormality.

An optimal Golgi transport system is important for mammalian cells. The adenosine diphosphate (ADP) ribosylation factors (ARF) are key proteins for regulating cargo sorting at the Golgi network. In this family, ARF3 mainly works at the trans-Golgi network (TGN), and no ARF3-related phenotypes have yet been described in humans. We here report the clinical and genetic evaluations of two unrelated children with de novo pathogenic variants in the ARF3 gene: c.200A > T (p.Asp67Val) and c.296G > T (p.Arg99Leu). Although the affected individuals presented commonly with developmental delay, epilepsy and brain abnormalities, there were differences in severity, clinical course and brain lesions. In vitro subcellular localization assays revealed that the p.Arg99Leu mutant localized to Golgi apparatus, similar to the wild-type, whereas the p.Asp67Val mutant tended to show a disperse cytosolic pattern together with abnormally dispersed Golgi localization, similar to that observed in a known dominant negative variant (p.Thr31Asn). Pull-down assays revealed that the p.Asp67Val had a loss-of-function effect and the p.Arg99Leu variant had increased binding of the adaptor protein, Golgi-localized, γ-adaptin ear-containing, ARF-binding protein 1 (GGA1), supporting the gain of function. Furthermore, in vivo studies revealed that p.Asp67Val transfection led to lethality in flies. In contrast, flies expressing p.Arg99Leu had abnormal rough eye, as observed in the gain-of-function variant p.Gln71Leu. These data indicate that two ARF3 variants, the possibly loss-of-function p.Asp67Val and the gain-of-function p.Arg99Leu, both impair the Golgi transport system. Therefore, it may not be unreasonable that they showed different clinical features like diffuse brain atrophy (p.Asp67Val) and cerebellar hypoplasia (p.Arg99Leu).

De Novo Truncating Variants in the Last Exon of SEMA6B Cause Progressive Myoclonic Epilepsy.

De novo variants (DNVs) cause many genetic diseases. When DNVs are examined in the whole coding regions of genes in next-generation sequencing analyses, pathogenic DNVs often cluster in a specific region. One such region is the last exon and the last 50 bp of the penultimate exon, where truncating DNVs cause escape from nonsense-mediated mRNA decay [NMD(-) region]. Such variants can have dominant-negative or gain-of-function effects. Here, we first developed a resource of rates of truncating DNVs in NMD(-) regions under the null model of DNVs. Utilizing this resource, we performed enrichment analysis of truncating DNVs in NMD(-) regions in 346 developmental and epileptic encephalopathy (DEE) trios. We observed statistically significant enrichment of truncating DNVs in semaphorin 6B (SEMA6B) (p value: 2.8 × 10-8; exome-wide threshold: 2.5 × 10-6). The initial analysis of the 346 individuals and additional screening of 1,406 and 4,293 independent individuals affected by DEE and developmental disorders collectively identified four truncating DNVs in the SEMA6B NMD(-) region in five individuals who came from unrelated families (p value: 1.9 × 10-13) and consistently showed progressive myoclonic epilepsy. RNA analysis of lymphoblastoid cells established from an affected individual showed that the mutant allele escaped NMD, indicating stable production of the truncated protein. Importantly, heterozygous truncating variants in the NMD(+) region of SEMA6B are observed in general populations, and SEMA6B is most likely loss-of-function tolerant. Zebrafish expressing truncating variants in the NMD(-) region of SEMA6B orthologs displayed defective development of brain neurons and enhanced pentylenetetrazole-induced seizure behavior. In summary, we show that truncating DNVs in the final exon of SEMA6B cause progressive myoclonic epilepsy.

Complete SAMD12 repeat expansion sequencing in a four-generation BAFME1 family with anticipation.

Benign adult familial myoclonic epilepsy type 1 (BAFME1) is an autosomal dominant, adult-onset neurological disease caused by SAMD12 repeat expansion. In BAFME1, anticipation, such as the earlier onset of tremor and/or seizures in the next generation, was reported. This could be explained by intergenerational repeat instability, leading to larger expansions in successive generations. We report a four-generation BAFME1-affected family with anticipation. Using Nanopore long-read sequencing, detailed information regarding the sizes, configurations, and compositions of the expanded SAMD12 repeats across generations was obtained. Unexpectedly, a grandmother-mother-daughter triad showed similar repeat structures but with slight repeat expansions, despite quite variable age of onset of seizures (range: 52-14 years old), implying a complex relationship between the SAMD12 repeat expansion sequence and anticipation. This study suggests that different factor(s) from repeat expansion could modify the anticipation in BAFME1.

Distal arthrogryposis in a girl arising from a novel TNNI2 variant inherited from paternal somatic mosaicism.

TNNI2 at 11p15.5 encodes troponin I2, fast skeletal type, which is a member of the troponin I gene family and a component of the troponin complex. Distal arthrogryposis (DA) is characterized by congenital limb contractures without primary neurological or muscular effects. DA is inherited in an autosomal dominant fashion and is clinically and genetically heterogeneous. Exome sequencing identified a causative variant in TNNI2 [NM_003282.4:c.532T>C p.(Phe178Leu)] in a Japanese girl with typical DA2b. Interestingly, the familial study using Sanger sequencing suggested a mosaic variant in her healthy father. Subsequent targeted amplicon-based deep sequencing detected the TNNI2 variant with variant allele frequencies of 9.4-17.7% in genomic DNA derived from peripheral blood leukocytes, saliva, hair, and nails in the father. We confirmed a disease-causing variant in TNNI2 in the proband inherited from her asymptomatic father with its somatic variant. Our case demonstrates that careful clinical and genetic evaluation is required in DA.

A novel homozygous CHMP1A variant arising from segmental uniparental disomy causes pontocerebellar hypoplasia type 8.

Pontocerebellar hypoplasia (PCH) is currently classified into 16 subgroups. Using mostly next-generation sequencing, pathogenic variants have been identified in as many as 24 PCH-associated genes. PCH type 8 (PCH8) is a rare heterogeneous disorder. Its clinical presentation includes severe development delay, increased muscle tone, microcephaly, and magnetic resonance imaging (MRI) abnormalities such as reduced cerebral white matter, a thin corpus callosum, and brainstem and cerebellar hypoplasia. To date, only two variants in the CHMP1A gene (MIM: 164010), NM_002768.5: c.88 C > T (p.Glu30*) and c.28-13 G > A, have been identified homozygously in seven patients with PCH8 from four families (MIM: 614961). CHMP1A is a subunit of the endosomal sorting complex required for transport III (ESCRT-III), which regulates the formation and release of extracellular vesicles. Biallelic CHMP1A loss of function impairs the ESCRT-III-mediated release of extracellular vesicles, which causes impaired progenitor proliferation in the developing brain. Herein, we report a patient with PCH8 who had a homozygous CHMP1A variant, c.122delA (p.Asn41Metfs*2), which arose from segmental uniparental disomy. Although our patient had similar MRI findings to those of previously reported patients, with no progression, we report some novel neurological and developmental findings that expand our knowledge of the clinical consequences associated with CHMP1A variants.

Long-read sequencing revealing intragenic deletions in exome-negative spastic paraplegias.

Hereditary spastic paraplegias (HSPs) are a heterogeneous group of neurodegenerative disorders characterized by progressive spasticity and weakness in the lower extremities. To date, a total of 88 types of SPG are known. To diagnose HSP, multiple technologies, including microarray, direct sequencing, multiplex ligation-dependent probe amplification, and short-read next-generation sequencing, are often chosen based on the frequency of HSP subtypes. Exome sequencing (ES) is commonly used. We used ES to analyze ten cases of HSP from eight families. We identified pathogenic variants in three cases (from three different families); however, we were unable to determine the cause of the other seven cases using ES. We therefore applied long-read sequencing to the seven undetermined HSP cases (from five families). We detected intragenic deletions within the SPAST gene in four families, and a deletion within PSEN1 in the remaining family. The size of the deletion ranged from 4.7 to 12.5 kb and involved 1-7 exons. All deletions were entirely included in one long read. We retrospectively performed an ES-based copy number variation analysis focusing on pathogenic deletions, but were not able to accurately detect these deletions. This study demonstrated the efficiency of long-read sequencing in detecting intragenic pathogenic deletions in ES-negative HSP patients.

Three KINSSHIP syndrome patients with mosaic and germline AFF3 variants.

AFF3 at 2q11.2 encodes the nuclear transcriptional activator AF4/FMR2 Family Member 3. AFF3 constitutes super elongation complex like 3, which plays a role in promoting the expression of genes involved in neurogenesis and development. The degron motif in AFF3 with nine highly conserved amino acids is recognized by E3 ubiquitin ligase to induce protein degradation. Recently, AFF3 missense variants in this region and variants featuring deletion including this region were identified and shown to cause KINSSHIP syndrome. In this study, we identified two novel and one previously reported missense variants in the degron of AFF3 in three unrelated Japanese patients. Notably, two of these three variants exhibited mosaicism in the examined tissues. This study suggests that mosaic variants also cause KINSSHIP syndrome, showing various phenotypes.

Repeat conformation heterogeneity in cerebellar ataxia, neuropathy, vestibular areflexia syndrome.

Cerebellar ataxia, neuropathy, vestibular areflexia syndrome (CANVAS) is a late-onset, slow-progressing multisystem neurodegenerative disorder. Biallelic AAGGG repeat expansion in RFC1 has been identified as causative of this disease, and repeat conformation heterogeneity (ACAGG repeat) was also recently implied. To molecularly characterize this disease in Japanese patients with adult-onset ataxia, we accumulated and screened 212 candidate families by an integrated approach consisting of flanking PCR, repeat-primed PCR, Southern blotting and long-read sequencing using Sequel II, GridION or PromethION. We identified 16 patients from 11 families, of whom seven had ACAGG expansions [(ACAGG)exp/(ACAGG)exp] (ACAGG homozygotes), two had ACAGG and AAGGG expansions [(ACAGG)exp/(AAGGG)exp] (ACAGG/AAGGG compound heterozygotes) and seven had AAGGG expansions [(AAGGG)exp/(AAGGG)exp] (AAGGG homozygotes). The overall detection rate was 5.2% (11/212 families including one family having two expansion genotypes). Long-read sequencers revealed the entire sequence of both AAGGG and ACAGG repeat expansions at the nucleotide level of resolution. Clinical assessment and neuropathology results suggested that patients with ACAGG expansions have similar clinical features to previously reported patients with homozygous AAGGG expansions, although motor neuron involvement was more notable in patients with ACAGG expansions (even if one allele was involved). Furthermore, a later age of onset and slower clinical progression were implied in patients with ACAGG/AAGGG compound heterozygous expansions compared with either ACAGG or AAGGG homozygotes in our very limited cohort. Our study clearly shows the occurrence of repeat conformation heterogeneity, with possible different impacts on the affected nervous systems. The difference in disease onset and progression between compound heterozygotes and homozygotes might also be suspected but with very limited certainty due to the small sample number of cases in our study. Studies of additional patients are needed to confirm this.

Patients with biallelic GGC repeat expansions in NOTCH2NLC exhibiting a typical neuronal intranuclear inclusion disease phenotype.

We report two patients with autosomal dominant neuronal intranuclear inclusion disease (NIID) harboring the biallelic GGC repeat expansion in NOTCH2NLC to uncover the impact of repeat expansion zygosity on the clinical phenotype. The zygosity of the entire NOTCH2NLC GGC repeat expansion and DNA methylation were comprehensively evaluated using fluorescent amplicon length PCR (AL-PCR), Southern blotting and targeted long-read sequencing, and detailed genetic/epigenetic and clinical features were described. In AL-PCR, we could not recognize the wild-type allele in both patients. Targeted long-read sequencing revealed that one patient harbored a homozygous repeat expansion. The other patient harbored compound heterozygous repeat expansions. The GGC repeats and the nearest CpG island were hypomethylated in all expanded alleles in both patients. Both patients harboring the biallelic GGC repeat expansion showed a typical dementia-dominant NIID phenotype. In conclusion, the biallelic GGC repeat expansion in two typical NIID patients indicated that NOTCH2NLC-related diseases could be completely dominant.

Pathogenic variants detected by RNA sequencing in Cornelia de Lange syndrome.

Recent studies suggest that transcript isoforms significantly overlap (approximately 60%) between brain tissue and Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs). Interestingly, 14 cohesion-related genes with variants that cause Cornelia de Lange Syndrome (CdLS) are highly expressed in the brain and LCLs. In this context, we first performed RNA sequencing of LCLs from 22 solved (with pathogenic variants) and 19 unsolved (with no confirmed variants) CdLS cases. Next, an RNA sequencing pipeline was developed using solved cases with two different methods: short variant analysis (for single-nucleotide and indel variants) and aberrant splicing detection analysis. Then, 19 unsolved cases were subsequently applied to our pipeline, and four pathogenic variants in NIPBL (one inframe deletion and three intronic variants) were newly identified. Two of three intronic variants were located at Alu elements in deep-intronic regions, creating cryptic exons. RNA sequencing with LCLs was useful for identifying hidden variants in exome-negative cases.