Published papers

Transient receptor potential canonical 6 (TRPC6) variants, which were initially detected in adult-onset familial focal segmental glomerulosclerosis (FSGS), were also identified in pediatric-onset one. Here, we present a patient with adult-onset steroid-resistant nephrotic syndrome (SRNS) who harbored a likely pathogenic TRPC6 variant and partially responded to calcineurin inhibitors (CNIs). A 44-year-old woman with stable rheumatoid arthritis, systemic lupus erythematosus, and Sjögren’s syndrome was presented with nephrotic syndrome. Her renal biopsy results showed minor glomerular abnormalities. Upon admission, she was treated with steroids for around 4 weeks, but it was ineffective. After 1–2 weeks of cyclosporine A (CyA) administration, urine output increased, renal function improved without a decrease in proteinuria, and she was discharged. Her renal function was maintained for 2 months, but after a CyA dose reduction, she was again admitted to the hospital due to relapsing edema, decreased urine output, and worsening renal function. CyA was replaced by tacrolimus (TAC). A second renal biopsy showed nearly the same findings as the first except for tubulointerstitial lesions. After 1–2 weeks of TAC administration, urine output increased, and renal function improved. However, urinary protein levels did not decrease as before. After discharge, a whole exome analysis revealed a heterozygous splice donor site variant NM_004621.6;c.2644 + 1G > A in TRPC6. Genetic testing identified a novel splice donor site variant of TRPC6 in a patient with adult-onset SRNS, which prevented unnecessary steroid continuation. The safety and efficacy of CNI in TRPC6 glomerulopathy must be evaluated in future larger studies with longer follow-up.

Aberrations in the COL4A1 gene lead to the faulty secretion of type IV collagen, resulting in damage to vascular basement membrane and various symptoms affecting multiple organs. Epilepsy associated with COL4A1-related disorder is often challenging to treat. We report a 5-year-old boy with developmental delay, congenital ocular abnormalities diagnosed with Peters anomaly, and abnormal neuroimaging findings. The patient had epileptic spasms starting at the age of 10 months and was diagnosed with West syndrome. Initial treatment, including adrenocorticotropic hormone(ACTH), was only transiently effective. At 1 year 9 months of age, vigabatrin(VGB)administration was introduced cautiously considering its potential adverse effects that might worsen ocular lesions. His seizures were suppressed with a modest dose(100 mg/kg/day)of VGB with no signs of electroretinogram abnormalities on repeated ophthalmological examinations. VGB was discontinued successfully at 3 years 6 months of age with no recurrence of seizures. A de novo pathogenic variant NM_001845.6:c.3620G> A(p.G1207E)was found in COL4A1. These findings suggest that VGB may be used effectively and safely in treating refractory West syndrome, even in patients with congenital ocular abnormalities associated with COL4A1-related disorder.

Type I interferon (IFN-I) plays an important role in the innate immune response through inducing IFN-I-stimulated genes (ISGs). However, how alternative splicing (AS) events, especially over time, affect their function remains poorly understood. We generated an annotation (113,843 transcripts) for IFN-I-stimulated human B cells called isoISG using high-accuracy long-read sequencing data from PacBio Sequel II/IIe. Transcript isoform profiling using isoISG revealed that isoform switching occurred in the early response to IFN-I so that ISGs would gain functional domains (e.g., C4B) or higher protein production (e.g., IRF3). Conversely, isoforms lacking functional domains increased during the late phase of IFN-I response, mainly due to intron retention events. This suggests that isoform switching both triggers and terminates IFN-I responses at the translation and protein levels. Furthermore, genetic variants influencing the isoform ratio of ISGs were associated with immunological and infectious diseases. AS has essential roles in regulating innate immune response and associated diseases.

Background: An alternative approach to investigate associations between genetic variants and disease is to examine deviations from the Hardy–Weinberg equilibrium (HWE) in genotype frequencies within a case population, instead of case–control association analysis. The HWE analysis requires disease cases and demonstrates a notable ability in mapping recessive variants. Allelic heterogeneity is a common phenomenon in diseases. While gene-based case–control association analysis successfully incorporates this heterogeneity, there are no such approaches for HWE analysis. Therefore, we proposed a gene-based HWE test (gene-HWT) by aggregating single-nucleotide polymorphism (SNP)-level HWE test statistics in a gene to address allelic heterogeneity. Results: This method used only genotype count data and publicly available linkage disequilibrium information and has a very low computational cost. Extensive simulations demonstrated that gene-HWT effectively controls the type I error at a low significance level and outperforms SNP-level HWE test in power when there are multiple causal variants within a gene. Using gene-HWT, we analyzed genotype count data from a genome-wide association study of six cancer types in Japanese individuals and suggest DGKE and ANO3 as potential germline factors in colorectal cancer. Furthermore, FSTL4 was suggested through a combined analysis across the six cancer types, with particularly notable associations observed in colorectal and prostate cancers. Conclusions: These findings indicate the potential of gene-HWT to elucidate the genetic basis of complex diseases, including cancer.

There are hundreds of rare syndromic diseases involving hearing loss, many of which are not targeted for clinical genetic testing. We systematically explored the genetic causes of undiagnosed syndromic hearing loss using a combination of whole exome sequencing (WES) and a phenotype similarity search system called PubCaseFinder. Fifty-five families with syndromic hearing loss of unknown cause were analyzed using WES after prescreening of several deafness genes depending on patient clinical features. Causative genes were identified in 22 families, including both established genes associated with syndromic hearing loss (PTPN11, CHD7, KARS1, OPA1, DLX5, MITF, SOX10, MYO7A, and USH2A) and those associated with nonsyndromic hearing loss (STRC, EYA4, and KCNQ4). Association of a DLX5 variant with incomplete partition type I (IP-I) anomaly of the inner ear was identified in a patient with cleft lip and palate and acetabular dysplasia. The study identified COL1A1, CFAP52, and NSD1 as causative genes through phenotype similarity search or by analogy. ZBTB10 was proposed as a novel candidate gene for syndromic hearing loss with IP-I. A mouse model with homozygous Zbtb10 frameshift variant resulted in embryonic lethality, suggesting the importance of this gene for early embryonic development. Our data highlight a wide spectrum of rare causative genes in patients with syndromic hearing loss, and demonstrate that WES analysis combined with phenotype similarity search is a valuable approach for clinical genetic testing of undiagnosed disease.

Since the publication of this article, it has been brought to our attention, that the identified mutation (NM_015277: c.2617 G > A; p.Glu873Lys) is identical with the mutation (NM_001144967: c.2677 G > A; p.Glu893Lys) reported by Broix et al (Nature Genetics 48, 1349–1358, 2016 https://doi.org/10.1038/ng.3676). Therefore the mutation is not novel but recurrent. Accordingly, the word “novel” should be deleted throughout the article including the title. Thus, the title should read “A missense mutation in the HECT domain of NEDD4L identified in a girl with periventricular nodular heterotopia, polymicrogyria, and cleft palate”.

Collapsin response mediator protein 2 (Crmp2) is an evolutionarily well-conserved tubulin-binding cytosolic protein that plays critical roles in the formation of neural circuitry in model organisms including zebrafish and rodents. No clinical evidence that CRMP2 variants are responsible for monogenic neurogenic disorders in humans presently exists. Here, we describe two patients with de novo non-synonymous variants (S14R and R565C) of CRMP2 and intellectual disability associated with hypoplasia of the corpus callosum. We further performed various functional assays of CRMP2 variants using zebrafish and zebrafish Crmp2 (abbreviated as z-CRMP2 hereafter) and an antisense morpholino oligonucleotide [AMO]-based experimental system in which crmp2-morphant zebrafish exhibit the ectopic positioning of caudal primary (CaP) motor neurons. Whereas the co-injection of wild-type z-CRMP2 mRNA suppressed the ectopic positioning of CaP motor neurons in Crmp2-morphant zebrafish, the co-injection of R566C or S15R, z-CRMP2, which corresponds to R565C and S14R of human CRMP2, failed to rescue the ectopic positioning. Transfection experiments of zebrafish or rat Crmp2 using plasmid vectors in HeLa cells, with or without a proteasome inhibitor, demonstrated that the expression levels of mutant Crmp2 protein encoded by R565C and S14R CRMP2 variants were decreased, presumably because of increased degradation by proteasomes. When we compared CRMP2-tubulin interactions using co-immunoprecipitation and cellular localization studies, the R565C and S14R mutations weakened the interactions. These results collectively suggest that the CRMP2 variants detected in the present study consistently led to the loss-of-function of CRMP2 protein and support the notion that pathogenic variants in CRMP2 can cause intellectual disabilities in humans.

When a de novo balanced reciprocal translocation is identified in patients with multiple congenital abnormalities, attempts are often made to infer the relationship between the phenotype of the patient and genes in the proximity of the breakpoint. Here, we report a patient with intellectual disability, atrial septal defect, syndactyly, and cleft lip and palate who had an “apparently balanced” de novo reciprocal translocation t(4:18)(q31;q11.2) as well as a 7-Mb cryptic deletion spanning the HOXD cluster on chromosome 2q31 that was unrelated to the reciprocal translocation. Further analysis using a nanopore long-read sequencer showed complex rearrangements on both derivative chromosomes 4 and 18 and the deleted chromosome 2. First, the TLL1 locus, which is associated with atrial septal defect, was disrupted by the rearrangement involving chromosome 4. Second, the deleted interval at 2q31 included the entire HOXD cluster, the deletion of which is known to cause toe syndactyly, and the DLX1 and DLX2 loci, which are responsible for cleft lip and palate. Among the haplo-sensitive genes within the deleted interval on 2q31, only the RAPGEF4 gene is known to be associated with an autistic phenotype. Hence, most of the clinical features of the patient could be ascribed to specific genomic rearrangements. We have shown the effectiveness of long-read sequencing in defining, in detail, the likely effects of an apparently balanced translocation and cryptic deletion. The results of the present analysis suggest the possibility of phenotypic prediction through a detailed analysis of structural abnormalities, including balanced translocations and deletions.

Objective: Adaptive sampling, a nanopore sequencing method that enriches regions of interest (ROI), is cost-effective and useful. However, the process of defining targeted regions and creating the corresponding definition file (.bed file) are time-consuming and laborious. To simplify this process, we have developed a tool to easily create a.bed file for adaptive sampling directly from gene names. Results: The tool is freely available on GitHub at https://github.com/medicalbioinfo/BED-Craft. The input is a text file containing one or more gene names (symbols), and even with a large number of input genes (e.g., thousands), the resulting.bed file is generated in less than a second. The length of the buffer region added upstream and downstream of the ROI is designed to account for genome strand orientation, ensuring efficient adaptive sampling. The buffer length can also be modified by the user. The tool supports the genomes of human hg19, hg38, T2T-CHM13, and other species. For researchers unfamiliar with command-line input, a GUI version of the tool is also available at https://keio-cmg.jp/BED-Craft/. This easy-to-use.bed file generation tool enables adaptive sampling by easily changing the target genes of interest in nanopore sequencing, and provide great benefits to researchers and diagnostic laboratories.

CCP110 (centriolar coiled coil protein 110, also known as CP110) is one of the essential proteins localized in the centrosome that plays critical roles in the regulation of the cell cycle and also in the initiation of ciliogenesis. So far, no human congenital disorders have been identified to be associated with pathogenic variants of CCP110. Mice with biallelic loss-of-function variants of Ccp110 (Ccp110−/−) are known to manifest multiple organ defects, including a small body size, polydactyly, omphalocele, congenital heart defects, cleft palate, short ribs, and a small thoracic cage, a pattern of abnormalities closely resembling that in “ciliopathies” in humans. Herein, we report a 7-month-old male infant who presented with growth failure and skeletal abnormalities, including a narrow thorax and severe brachydactyly. Trio exome analysis of the genomic DNA of the patient and his parents showed that the patient was a compound heterozygote for truncating variants of CCP110, including a frameshift variant NM_001323572.2:c.856_857del, p.(Val286Leufs*5) inherited from the father, and a nonsense variant NM_001323572.2:c.1129C>T, p.(Arg377*) inherited from the mother. The strikingly similar pattern of malformations between Ccp110−/− mice and the 7-month-old male infant reported herein carrying unequivocal truncating CCP110 variants strongly supports the contention that CCP110 is a novel disease-causative gene.

Patatin-like phospholipase domain-containing lipase 8 (PNPLA8), one of the calcium-independent phospholipase A2 enzymes, is involved in various physiological processes through the maintenance of membrane phospholipids. Biallelic variants in PNPLA8 have been associated with a range of paediatric neurodegenerative disorders. However, the phenotypic spectrum, genotype–phenotype correlations and the underlying mechanisms are poorly understood. Here, we newly identified 14 individuals from 12 unrelated families with biallelic ultra-rare variants in PNPLA8 presenting with a wide phenotypic spectrum of clinical features. Analysis of the clinical features of current and previously reported individuals (25 affected individuals across 20 families) showed that PNPLA8-related neurological diseases manifest as a continuum ranging from variable developmental and/or degenerative epileptic–dyskinetic encephalopathy to childhood-onset neurodegeneration. We found that complete loss of PNPLA8 was associated with the more profound end of the spectrum, with congenital microcephaly. Using cerebral organoids generated from human induced pluripotent stem cells, we found that loss of PNPLA8 led to developmental defects by reducing the number of basal radial glial cells and upper-layer neurons. Spatial transcriptomics revealed that loss of PNPLA8 altered the fate specification of apical radial glial cells, as reflected by the enrichment of gene sets related to the cell cycle, basal radial glial cells and neural differentiation. Neural progenitor cells lacking PNPLA8 showed a reduced amount of lysophosphatidic acid, lysophosphatidylethanolamine and phosphatidic acid. The reduced number of basal radial glial cells in patient-derived cerebral organoids was rescued, in part, by the addition of lysophosphatidic acid. Our data suggest that PNPLA8 is crucial to meet phospholipid synthetic needs and to produce abundant basal radial glial cells in human brain development.

Cockayne syndrome (CS) is a rare autosomal recessive disorder caused by mutations in ERCC8 or ERCC6. Most pathogenic variants in ERCC8 are single nucleotide substitutions. Structural variants (SVs) have been reported in patients with ERCC8-related CS. However, comprehensive molecular detection, including SVs of ERCC8, in CS patients remains problematic. Herein, we present three Japanese patients with ERCC8-related CS in whom causative SVs were identified using whole-exome-based copy number variation (CNV) detection tools. One patient showed compound heterozygosity for a 259-kb deletion and a deletion of exon 4 which has previously been reported as an Asia-specific variant. The other two patients were homozygous for the same exon 4 deletion. The exon 4 deletion was detected only by the ExomeDepth software. Intrigued by the discrepancy in the detection capability of various tools for the SVs, we evaluated the analytic performance of four whole-exome-based CNV detection tools using an exome data set from 337 healthy individuals. A total of 1,278,141 exons were predicted as being affected by the 4 CNV tools. Interestingly 95.1% of these affected exons were detected by one tool alone. Thus, we expect that the use of multiple tools may improve the detection rate of SVs from aligned exome data.

Congenital Heart Defect and Ectodermal Dysplasia (CHDED) is an autosomal dominant disorder caused by the PRKD1 gene. CHDED is characterized by heart defects and ectodermal dysplasia. To date, eight patients with CHDED have been described. Calcifications were present in three patients with CHDED. (two patients; renal calcifications, one patient; brain calcifications). The organ distribution of calcifications in CHDED has been unclear. We report here another patient with CHDED and brain calcifications. The patient was a 9-month-old Japanese girl. She presented with heart defects and ectodermal dysplasia. At 6 months of age, she had generalized seizures, and a CT scan revealed calcifications in the bilateral deep cerebral white matter. The seizures resolved with the administration of levetiracetam. The patient had a de novo, heterozygous pathogenic variant, c.1808G > A, p.(Arg603His), in the PRKD1 gene. Together with the previously reported patients mentioned above, we demonstrated the role of the PRKD1 variant in brain calcification. We propose that PRKD1 and two genes, ITGB2 and JAM2, which are known to be associated with brain calcification, act through a common signaling pathway abnormality. In support of our hypothesis, there are some experimental results that link PRKD1 and JAM2. PRKD1 functions with the integrin ITGB2 as a partner. JAM2, which is associated with brain calcification and is critical for maintaining of the tight junction of the endothelial cells, interacts with integrins including ITGB2. Therefore, PRKD1 could lead to the pathological phenotype of brain calcification.

The recognition of syndromic forms of cleft palate is important for condition-specific management. Here, we report a patient with cleft palate, congenital heart disease, intellectual disability, and café-au-lait spots who had a deletion of chromosome 15q14. The identification of the precise breakpoints using a Nanopore-based long-read sequencer showed that the deletion spanned MEIS2 and SPRED1 loci. Cleft palate and café-au-lait spots can be ascribed to MEIS2 and SPRED1, respectively. Patients with cleft palate and café-au-lait spots should be encouraged to undergo a detailed genomic evaluation, including screening for a 15q14 deletion, to enable appropriate anticipatory medico-surgical management and genetic counseling.

BACKGROUND: The cause of chronic thromboembolic pulmonary hypertension (CTEPH) remains largely unknown. Recently, clonal hematopoiesis (CH) has been reported to be associated with cardiovascular and thromboembolic diseases. Here, we investigated the prevalence and clinical impact of CH in patients with CTEPH. METHODS AND RESULTS: Whole-exome sequencing and deep-panel sequencing were performed in 214 patients with CTEPH. Clinical data before and after treatment were compared between patients with and without CH. RNA sequencing and serum analysis were performed to explore the pathogenesis that CH contributes to CTEPH. Among the enrolled patients, 20.1%, notably 44.4% who were 80 to 89 years old, had variants in CH-associated genes. In regard to clinical impact, B-type natriuretic peptide levels and home oxygen therapy rate were significantly higher, and 6-minute walk distance was significantly shorter after treatment in patients with CH than in those without CH. Moreover, novel clot reformation in the pulmonary artery despite the use of anticoagulants and additional angioplasty events after treatment completion were more frequent in patients with CH. RNA sequencing analysis revealed that blood coagulation and neutrophil extracellular trap formation pathways were enriched in patients with CH. Additionally, serum citrullinated histone H3 levels were higher in patients with CH than those without CH. These results were consistent in the subgroup of patients who did not have the history of hematological disorders. CONCLUSIONS: The findings in this study raise the possibility that CH will induce a more prothrombotic state through neutrophil activation and neutrophil extracellular trap formation, contributing to pathogenesis and poor treatment response in patients with CTEPH.