A case of Bloom syndrome manifesting with therapy-related myelodysplastic syndromes harboring a novel BLM gene variant.

Bloom syndrome (BS) is an autosomal recessive genetic disorder caused by variants in the BLM gene. BS is characterized by distinct facial features, elongated limbs, and various dermatological complications including photosensitivity, poikiloderma, and telangiectatic erythema. The BLM gene encodes a RecQ helicase critical for genome maintenance, stability, and repair, and a deficiency in functional BLM protein leads to genomic instability and high predisposition to various types of cancers, particularly hematological and gastrointestinal malignancies. Here, we report a case of BS with a previously unreported variant in the BLM gene. The patient was a 34-year-old woman who presented with short stature, prominent facial features, and a history of malignancies, including lymphoma, breast cancer, and myelodysplastic syndromes (MDS). She was initially treated with azacitidine for MDS and showed transient improvement, but eventually died at age of 35 due to progression of MDS. Genetic screening revealed compound heterozygous variants in the BLM gene, with a recurrent variant previously reported in BS in one allele and a previously unreported variant in the other allele. Based on her characteristic clinical features and the presence of heterozygous variants in the BLM gene, she was diagnosed with BS harboring compound heterozygous BLM variants.

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.

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.

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.

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.

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.

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.

Protective effect of hydroxychloroquine on infections in patients with systemic lupus erythematosus: an observational study using the LUNA registry.

Objectives: Infection is a leading cause of death in patients with systemic lupus erythematosus (SLE). Alt hough hydroxychloroquine (HCQ) has been reported to inhibit infection, evidence from Asian populations remains insufficient. We investigated this effect in Japanese SLE patients. Methods: Data from the Lupus Registry of Nationwide Institutions were used in this study. The patients were ≥20 years old and met the American College of Rheumatology (ACR) classification criteria revised in 1997. We defined "severe infections" as those requiring hospitalization. We analyzed the HCQ's effect on infection suppression using a generalized estimating equation (GEE) logistic regression model as the primary endpoint and performed a survival analysis for the duration until the first severe infection. Results: Data from 925 patients were used (median age, 45 [interquartile range 35-57] years; female, 88.1%). GEE analysis revealed that severe infections were significantly associated with glucocorticoid dose (odds ratio [OR] 1.968 [95% confidence interval, 1.379-2.810], p<0.001), immunosuppressants (OR 1.561 [1.025-2.380], p=0.038), and baseline age (OR 1.043 [1.027-1.060], p<0.001). HCQ tended to suppress severe infections, although not significantly (OR 0.590 [0.329-1.058], p=0.077). Survival time analysis revealed a lower incidence of severe infections in the HCQ group than in the non-HCQ group (p<0.001). In a Cox proportional hazards model, baseline age (hazard ratio [HR] 1.029 [1.009-1.050], p=0.005) and HCQ (HR 0.322 [0.142-0.728], p=0.006) were significantly related to incidence. Conclusion: HCQ may help extend the time until the occurrence of infection complications and tends to decrease infection rates.

Long-term remission of VEXAS syndrome achieved by a single course of CHOP therapy: A case report.

We herein describe the case of a 52-year-old male patient who presented with fever, arthritis, and neutrophilic dermatosis in 2013 and subsequently experienced macrophage activation syndrome treated with high-dose glucocorticoid therapy. Due to the persistent symptoms refractory to several immunomodulatory and immunosuppressive (IS) drug therapies with dapsone, methotrexate, tacrolimus, infliximab (IFX), and tocilizumab (TCZ), he received prednisolone (PSL) ≥20 mg/day to suppress disease activity. In 2017, Epstein-Barr virus (EBV)-associated haemophagocytic lymphohistiocytosis (HLH) was diagnosed and initially treated with immunochemotherapy consisting of dexamethasone, cyclosporine (CyA), and etoposide (ET). Because of the suboptimal response to the initial therapy, cytoreduction therapy consisting of CHOP (combination chemotherapy consisting of cyclophosphamide, doxorubicin, vincristine, and PSL) was administered. This regimen improved the EBV-associated HLH. Later, the patient's condition stabilised with methylprednisolone 1 mg/day and CyA 100 mg/day. In 2022, ubiquitylation-initiating E1 enzyme (UBA1) variant analysis using Sanger sequencing of peripheral blood leukocytes detected a previously reported somatic variant (NM_003334.3: c.118-1G>C), confirming the diagnosis of vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic (VEXAS) syndrome. The clinical course in the present case suggested the possibility that CHOP could be a potential treatment option for VEXAS syndrome, in the pathophysiology of which the expansion of clones with UBA1 variant seems to play a pivotal role.

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.

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.

X-linked intellectual disability related to a novel variant of KLHL15.

Kelch-like (KLHL) 15, localized on chromosome Xp22.11, was recently identified as an X-linked intellectual disability gene. Herein, we report a case of a male patient with a novel nonsense variant, c.736 C > T p.(Arg246*), in KLHL15, who presented with impaired intelligence, short stature, frequent hypoglycemia, and periodic fever. Patients with nonsense variants in KLHL15 may develop intellectual disabilities, minor skeletal anomalies, and facial dysmorphisms.

Biallelic structural variations within FGF12 detected by long-read sequencing in epilepsy.

We discovered biallelic intragenic structural variations (SVs) in FGF12 by applying long-read whole genome sequencing to an exome-negative patient with developmental and epileptic encephalopathy (DEE). We also found another DEE patient carrying a biallelic (homozygous) single-nucleotide variant (SNV) in FGF12 that was detected by exome sequencing. FGF12 heterozygous recurrent missense variants with gain-of-function or heterozygous entire duplication of FGF12 are known causes of epilepsy, but biallelic SNVs/SVs have never been described. FGF12 encodes intracellular proteins interacting with the C-terminal domain of the alpha subunit of voltage-gated sodium channels 1.2, 1.5, and 1.6, promoting excitability by delaying fast inactivation of the channels. To validate the molecular pathomechanisms of these biallelic FGF12 SVs/SNV, highly sensitive gene expression analyses using lymphoblastoid cells from the patient with biallelic SVs, structural considerations, and Drosophila in vivo functional analysis of the SNV were performed, confirming loss-of-function. Our study highlights the importance of small SVs in Mendelian disorders, which may be overlooked by exome sequencing but can be detected efficiently by long-read whole genome sequencing, providing new insights into the pathomechanisms of human diseases.

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.

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.

A missense variant at the RAC1-PAK1 binding site of RAC1 inactivates downstream signaling in VACTERL association.

RAC1 at 7p22.1 encodes a RAC family small GTPase that regulates actin cytoskeleton organization and intracellular signaling pathways. Pathogenic RAC1 variants result in developmental delay and multiple anomalies. Here, exome sequencing identified a rare de novo RAC1 variant [NM_018890.4:c.118T > C p.(Tyr40His)] in a male patient. Fetal ultrasonography indicated the patient to have multiple anomalies, including persistent left superior vena cava, total anomalous pulmonary venous return, esophageal atresia, scoliosis, and right-hand polydactyly. After birth, craniofacial dysmorphism and esophagobronchial fistula were confirmed and VACTERL association was suspected. One day after birth, the patient died of respiratory failure caused by tracheal aplasia type III. The molecular mechanisms of pathogenic RAC1 variants remain largely unclear; therefore, we biochemically examined the pathophysiological significance of RAC1-p.Tyr40His by focusing on the best characterized downstream effector of RAC1, PAK1, which activates Hedgehog signaling. RAC1-p.Tyr40His interacted minimally with PAK1, and did not enable PAK1 activation. Variants in the RAC1 Switch II region consistently activate downstream signals, whereas the p.Tyr40His variant at the RAC1-PAK1 binding site and adjacent to the Switch I region may deactivate the signals. It is important to accumulate data from individuals with different RAC1 variants to gain a full understanding of their varied clinical presentations.

Clinical and genetic features of Japanese cases of MDS associated with VEXAS syndrome.

VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome is a new disease entity with autoinflammatory disorders (AID) driven by somatic variants in UBA1 that frequently co-exists with myelodysplastic syndromes (MDS). Clinicopathological and molecular features of Japanese cases with VEXAS-associated MDS remain elusive. We previously reported high prevalence of UBA1 variants in Japanese patients with relapsing polychondritis, in which 5 cases co-occurred with MDS. Here, we report clinicopathological and variant profiles of these 5 cases and 2 additional cases of MDS associated with VEXAS syndrome. Clinical characteristics of these cases included high prevalence of macrocytic anemia with marked cytoplasmic vacuoles in myeloid/erythroid precursors and low bone marrow (BM) blast percentages. All cases were classified as low or very low risk by the revised international prognostic scoring system (IPSS-R). Notably, 4 out of 7 cases showed significant improvement of anemia by treatment with prednisolone (PSL) or cyclosporin A (CsA), suggesting that an underlying inflammatory milieu induced by VEXAS syndrome may aggravate macrocytic anemia in VEXAS-associated MDS. Targeted deep sequencing of blood samples suggested that MDS associated with VEXAS syndrome tends to involve a smaller number of genes and lower risk genetic lesions than classical MDS.

An integrated genetic analysis of epileptogenic brain malformed lesions.

Focal cortical dysplasia is the most common malformation during cortical development, sometimes excised by epilepsy surgery and often caused by somatic variants of the mTOR pathway genes. In this study, we performed a genetic analysis of epileptogenic brain malformed lesions from 64 patients with focal cortical dysplasia, hemimegalencephy, brain tumors, or hippocampal sclerosis. Targeted sequencing, whole-exome sequencing, and single nucleotide polymorphism microarray detected four germline and 35 somatic variants, comprising three copy number variants and 36 single nucleotide variants and indels in 37 patients. One of the somatic variants in focal cortical dysplasia type IIB was an in-frame deletion in MTOR, in which only gain-of-function missense variants have been reported. In focal cortical dysplasia type I, somatic variants of MAP2K1 and PTPN11 involved in the RAS/MAPK pathway were detected. The in-frame deletions of MTOR and MAP2K1 in this study resulted in the activation of the mTOR pathway in transiently transfected cells. In addition, the PTPN11 missense variant tended to elongate activation of the mTOR or RAS/MAPK pathway, depending on culture conditions. We demonstrate that epileptogenic brain malformed lesions except for focal cortical dysplasia type II arose from somatic variants of diverse genes but were eventually linked to the mTOR pathway.