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.

Exome-wide benchmark of difficult-to-sequence regions using short-read next-generation DNA sequencing.

Next-generation DNA sequencing (NGS) in short-read mode has recently been used for genetic testing in various clinical settings. NGS data accuracy is crucial in clinical settings, and several reports regarding quality control of NGS data, primarily focusing on establishing NGS sequence read accuracy, have been published thus far. Variant calling is another critical source of NGS errors that remains unexplored at the single-nucleotide level despite its established significance. In this study, we used a machine-learning-based method to establish an exome-wide benchmark of difficult-to-sequence regions at the nucleotide-residue resolution using 10 genome sequence features based on real-world NGS data accumulated in The Genome Aggregation Database (gnomAD) of the human reference genome sequence (GRCh38/hg38). The newly acquired metric, designated the 'UNMET score,' along with additional lines of structural information from the human genome, allowed us to assess the sequencing challenges within the exonic region of interest using conventional short-read NGS. Thus, the UNMET score could provide a basis for addressing potential sequential errors in protein-coding exons of the human reference genome sequence GRCh38/hg38 in clinical sequencing.

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.

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.

Biallelic missense CEP55 variants cause prenatal MARCH syndrome.

CEP55 encodes centrosomal protein 55 kDa, which plays a crucial role in mitosis, particularly cytokinesis. Biallelic CEP55 variants cause MARCH syndrome (multinucleated neurons, anhydramnios, renal dysplasia, cerebellar hypoplasia and hydranencephaly). Here, we describe a Japanese family with two affected siblings harboring novel compound heterozygous CEP55 variants, NM_001127182: c.[1357 C > T];[1358 G > A] p.[(Arg453Cys)];[(Arg453His)]. Both presented clinically with typical lethal MARCH syndrome. Although a combination of missense and nonsense variants has been reported previously, this is the first report of biallelic missense CEP55 variants. These variants biallelically affected the same amino acid, Arg453, in the last 40 amino acids of CEP55. These residues are functionally important for CEP55 localization to the midbody during cell division, and may be associated with severe clinical outcomes. More cases of pathogenic CEP55 variants are needed to establish the genotype-phenotype correlation.

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.

Pregnancy outcomes in patients with familial Mediterranean fever: systematic review and meta-analysis.

OBJECTIVE: The relationship between FMF and pregnancy outcomes remains unclear. This systematic review and meta-analysis aimed to clarify this association. METHODS: Electronic databases-PubMed, Web of Science, Cochrane, and EMBASE-were searched on 20 December 2022, using specific search terms. Case-control, cohort, and randomized clinical trial studies comparing patients with FMF and healthy controls were considered eligible. We excluded systematic reviews, meta-analyses, case series with fewer than five cases, republished articles without new findings on pregnancy outcomes, studies targeting paternal FMF, and those not published in English. The results were summarized in the form of odds ratios (ORs) and 95% CIs, using a random-effects model. This study was registered in the University hospital Medical Information Network Clinical Trials Registry (Japan) as UMIN000049827. RESULTS: The initial electronic search identified 611 records, of which 9 were included in this meta-analysis (177 735 pregnancies, 1242 with FMF, and 176 493 healthy controls). FMF was significantly associated with increased odds of preterm deliveries (OR, 1.67; 95% CI, 1.05-2.67; I2 = 22%) and insignificantly associated with increased odds of fetal growth restriction (OR, 1.45; 95% CI, 0.90-2.34; I2 = 0%) and hypertensive disorders during pregnancy (OR, 1.28; 95% CI, 0.87-1.87; I2 = 0%). CONCLUSION: FMF was significantly associated with preterm delivery and insignificantly associated with fetal growth restriction and hypertensive disorders. All of the included studies were observational studies. Treatment characteristics were not fully collected from the articles, and further analysis of treatments for FMF in pregnancy is still warranted.

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.

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).

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.

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.

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.

Efficient detection of somatic UBA1 variants and clinical scoring system predicting patients with variants in VEXAS syndrome.

OBJECTIVES: To efficiently detect somatic UBA1 variants and establish a clinical scoring system predicting patients with pathogenic variants in VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome. METHODS: Eighty-nine Japanese patients with clinically suspected VEXAS syndrome were recruited [81 males and 8 females; median onset age (IQR) 69.3 years (62.1-77.6)]. Peptide nucleic acid-clamping PCR (PNA-PCR), regular PCR targeting exon 3 clustering UBA1 variants, and subsequent Sanger sequencing were conducted for variant screening. Partitioning digital PCR (pdPCR) or targeted amplicon deep sequencing (TAS) was also performed to evaluate the variant allele frequency (VAF). We developed our clinical scoring system to predict UBA1 variant-positive and ­negative patients and assessed the diagnostic value of our system using receiver operating characteristic (ROC) curve analysis. RESULTS: Forty patients with reported pathogenic UBA1 variants (40/89, 44.9%) were identified, including a case having a variant with VAF of 1.7%, using a highly sensitive method. Our clinical scoring system considering >50 years of age, cutaneous lesions, lung involvement, chondritis, and macrocytic anaemia efficiently predicted patients with UBA1 variants (the area under the curve for the scoring total was 0.908). CONCLUSIONS: Genetic screening with the combination of regular PCR and PNA-PCR detected somatic UBA1 variants with high sensitivity and specificity. Our scoring system could efficiently predict patients with UBA1 variants.

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.