DNA methylation changes in the genome of patients with hypogonadotropic hypogonadism.

Although some Mendelian neurodevelopmental disorders have been shown to entail specific DNA methylation changes designated as epi-signatures, it remains unknown whether epi-signatures are consistent features of other genetic disorders. Here, we analyzed DNA methylation profiles of patients with hypogonadotropic hypogonadism (HH), a rare neuroendocrine disorder typically caused by monogenic or oligogenic mutations. First, we performed microarray-based genome-wide methylation analyses of nine patients with HH due to ANOS1, SOX2, or SOX10 variants and 12 control individuals. The results showed that 1118 probes were differentially methylated in one or more patients. The differentially methylated probes were highly variable among patients. No significant methylation changes were observed in genes functionally associated with ANOS1, SOX2, or SOX10. Then, we performed pyrosequencing of six selected CpG sites in the nine patients and 35 additional HH patients. The results of the patients were compared with those of 48 fertile men. There were no common methylation changes among these patients, with the exception of hypermethylation of two CpG sites in the ZNF245 promoter of three patients. Hypermethylation of the promoter has previously been reported as a very rare epigenetic polymorphism in the general population. These results indicate that genomes of HH patients have considerable DNA methylation changes; however, these changes are more likely to be physiological epigenetic variations than disease-specific epi-signatures. Our data suggest a possible association between hypermethylation of the ZNF254 promoter and HH, which needs to be examined in future studies.

Comprehensive study on central precocious puberty: molecular and clinical analyses in 90 patients.

CONTEXT: Defects in MKRN3, DLK1, KISS1, and KISS1R and some disorders, such as Temple syndrome (TS14), cause central precocious puberty (CPP). Recently, pathogenic variants (PVs) in MECP2 have been reported to be associated with CPP. OBJECTIVE: We aimed to clarify the contribution of (epi)genetic abnormalities to CPP and clinical and hormonal features in each etiology. SUBJECTS AND METHODS: We conducted targeted sequencing for MKRN3, DLK1, MECP2, KISS1, and KISS1R and methylation analysis for screening of imprinting disorders such as TS14 associated with CPP in 90 patients with CPP (no history of brain injuries and negative brain MRI) and collected their clinical and laboratory data. We measured serum DLK1 levels in three patients with TS14 and serum MKRN3 levels in two patients with MKRN3 genetic defects, together with some etiology-unknown patients with CPP and controls. RESULTS: We detected eight patients with TS14 (six, epimutation; one, mosaic maternal uniparental disomy chromosome 14; one, microdeletion) and three patients with MKRN3 genetic defects (one, PV; one, 13-bp deletion in the 5'-untranslated region (5'-UTR); one, microdeletion) with family histories of paternal early puberty. There were no patients with PVs identified in MECP2, KISS1, or KISS1R. We confirmed low serum MKRN3 level in the patient with a deletion in 5'-UTR. The median height at initial evaluation of TS14 patients was lower than that of all patients. Six patients with TS14 were born small for gestational age (SGA). CONCLUSION: (Epi)genetic causes were identified in 12.2% of patients with CPP at our center. For patients with CPP born SGA or together with family histories of paternal early puberty, (epi)genetic testing for TS14 and MKRN3 genetic defects should be considered. (271/250).

Compound heterozygous KCTD19 variants in a man with isolated nonobstructive azoospermia.

Case: A 40-year-old Japanese man with nonobstructive azoospermia (NOA) was found to carry rare variants in KCTD19, a newly identified causative gene for spermatogenic failure. This patient was identified through mutation screening of KCTD19 in 97 men with etiology-unknown isolated NOA. Outcome: The patient had two heterozygous variants in KCTD19 that affect consensus sequences of splice-donor sites [c.300+2T>A and c.2667C>T (p.E889E)]. Both variants were predicted to cause exon skipping. Long-read sequencing confirmed the compound heterozygosity of the variants. The patient exhibited small testes and a mildly elevated level of follicle-stimulating hormone but no other phenotypic abnormalities. Testicular histology showed borderline findings between spermatocyte maturation arrest and severe hypospermatogenesis. Conclusion: These results provide evidence that biallelic loss-of-function variants of KCTD19 represent rare causes of isolated NOA.

PTPN11 and FLNA variants in a boy with ambiguous genitalia, short stature, and non-specific dysmorphic features.

Noonan syndrome is a congenital disorder characterized by distinctive facial appearance, congenital heart defects, short stature, and skeletal dysplasia. Although boys with Noonan syndrome frequently exhibit cryptorchidism, a mild form of 46,XY disorders of sex development (DSD), they barely manifest more severe genital abnormalities. Here, we report a boy with ambiguous genitalia, short stature, and non-specific dysmorphic features. He had no cardiac abnormalities or skeletal dysplasia. His score in the Noonan syndrome diagnostic criteria (36 of 157 points, 23%) was lower than the cutoff for diagnosis (50%). Whole-exome sequencing identified a de novo heterozygous variant (c.922A>G: p.Asn308Asp) in PTPN11 and a maternally inherited hemizygous variant (c.1439C>T: p.Pro480Leu) in FLNA. The PTPN11 variant was a known causative mutation for Noonan syndrome. FLNA is a causative gene for neurodevelopmental and skeletal abnormalities and has also been implicated in 46,XY DSD. The p.Pro480Leu variant of FLNA was assessed as deleterious by in silico analyses. These results provide evidence that whole-exome sequencing is a powerful tool for diagnosing patients with atypical disease manifestations. Furthermore, our data suggest a possible role of digenic mutations as phenotypic modifiers of Noonan syndrome.

Homozygous 6-bp deletion of IGFALS in a prepubertal boy with short stature.

Biallelic IGFALS variants lead to acid‒labile subunit (ALS) deficiency characterized by growth hormone resistance with or without delayed puberty. Here, we report a prepubertal boy with a homozygous 2-amino acid deletion within the fourth N-glycosylation motif (c.1103_1108del, p.N368_S370delinsT) associated with parental consanguinity. He showed short stature consistent with ALS deficiency. This case expands the mutation spectrum of IGFALS to include the elimination of only one N-glycosylation motif of ALS.

Long-read next-generation sequencing for molecular diagnosis of pediatric endocrine disorders.

Recent advances in long-read next-generation sequencing (NGS) have enabled researchers to identify several pathogenic variants overlooked by short-read NGS, array-based comparative genomic hybridization, and other conventional methods. Long-read NGS is particularly useful in the detection of structural variants and repeat expansions. Furthermore, it can be used for mutation screening in difficultto- sequence regions, as well as for DNA-methylation analyses and haplotype phasing. This mini-review introduces the usefulness of long-read NGS in the molecular diagnosis of pediatric endocrine disorders.

SHOX and sex difference in height: a hypothesis.

The mean height is taller in males than in females, except for early teens. In this regard, previous studies have revealed that (1) distribution of the mean adult heights in subjects with disorders accompanied by discordance between sex chromosome complement and bioactive sex steroids and in control subjects (the British height standards) indicates that, of the ~12.5 cm of sex difference in the mean adult height, ~9 cm is accounted for by the difference in the sex chromosome complement and the remaining ~3.5 cm is explained by the dimorphism in sex steroids (primarily due to the growth-promoting effect of gonadal androgens); (2) according to the infancy-childhood-puberty growth model, the sex difference in the childhood growth function produces height differences of ~1 cm in childhood and 8-10 cm at 18-20 years of age, whereas the sex difference in the pubertal growth function yields height difference of ~4.5 cm at 18-20 years of age; and (3) SHOX expression and methylation analyses using knee cartilage tissues and cultured chondrocytes have shown lower SHOX expression levels in female samples than in male samples and methylation patterns consistent with partial spreading of X-inactivation affecting SHOX in female samples. These findings suggest that small but persistent sex difference in SHOX expression dosage leads to the variation in the sex steroid independent childhood growth function, thereby yielding the sex difference in height which remains small in childhood but becomes obvious in adulthood.

Isodicentric Y Chromosome with Multiple Breakpoints in the Pseudoautosomal Region 1.

INTRODUCTION: Isodicentric Y chromosomes are relatively common structural variants of the human genome. The underlying mechanism of isodicentric Y chromosomes with short arm breakpoints [idic(Yq)] remains to be clarified. CASE PRESENTATION: We encountered a Japanese man with azoospermia and mild short stature. G-banding and array-based comparative genomic hybridization indicated that his karyotype was 45,X/46,X,idic(Y)(qter→p11.32::p11.32→qter) with a ∼1.8 Mb terminal deletion. Whole-genome sequencing suggested that the Y chromosome had four breakpoints in a ∼7 kb region of the pseudoautosomal region 1 (PAR1). CONCLUSION: This case was assumed to have an idic(Yq) resulting from multiple DNA double-strand breaks in PAR1. This rearrangement may have been facilitated by the PAR1-specific chromatin architecture. The clinical features of the patient can be ascribed to SHOX haploinsufficiency and the presence of a 45,X cell line, although copy-number gains of some Yq genes and the size reduction of PAR1 may also contribute to his spermatogenic failure.

Machine learning trial to detect sex differences in simple sticker arts of 1606 preschool children.

BACKGROUND: Previous studies suggested that drawings made by preschool boys and girls show distinguishable differences. However, children's drawings on their own are too complexly determined and inherently ambiguous to be a reliable indicator. In the present study, we attempted to develop a machine learning algorithm for classification of sex of the subjects using children's artworks. METHODS: We studied three types of simple sticker artworks from 1606 Japanese preschool children aged 51-83 months (803 boys and 803 girls). Those artworks were processed into digitalized data. Simulated data based on the original data were also generated. Logistic regression approach was applied to each dataset to make a classifier, and run on each dataset in a stratified ten-fold cross-validation with hyperparameter tuning. A probability score was calculated in each sample and utilized for sex classification. Prediction performance was evaluated using accuracy, recall, and precision scores, as well as learning curves. RESULTS: Two models created from the original and simulated data showed comparably low metrics. The distributions of probability scores in the samples from boys and girls mostly overlapped and were indistinguishable. Learning curves of the models showed an extremely under-fitted pattern. CONCLUSIONS: Our machine learning algorithm was unable to distinguish simple sticker arts created by boys and girls. More complex tasks will enable to develop an accurate classifier.

Molecular and Clinical Features of Congenital Hypothyroidism Due to Multiple DUOX2 Variants.

Background: DUOX2 is one of the major causative genes of congenital hypothyroidism (CH). Still, the mutation spectrum and clinical outcomes of biallelic DUOX2 variants are not fully understood. This study aimed to elucidate the molecular features and long-term clinical manifestations of CH caused by multiple pathogenic DUOX2 variants. Methods: A total of 255 patients with CH were screened for rare variants of 11 known causative genes. DUOX2 variants were classified according to their protein structure and residual activity. In vitro assays were performed for several variants of unknown functions. Clinical analyses were conducted for patients with multiple pathogenic variants of DUOX2 but not of other genes. Results: We identified 24 pathogenic variants of DUOX2, together with two benign variants and seven variants of uncertain significance, in 63 patients. The pathogenic variants included three missense substitutions and one frameshift variant that have not yet been linked to CH. Twenty-one patients carried multiple pathogenic DUOX2 variants without any other pathogenic gene variants. Three of the 21 patients harbored homozygous variants. Family analysis, long-read amplicon sequencing, and haplotype phasing confirmed compound heterozygosity of the DUOX2 variants in 14 patients, whereas the allelic positions of the variants in the remaining four patients could not be determined. Of the 21 patients, 19 were treated with levothyroxine; their ages at drug withdrawal ranged from 9 months to 21.4 years. Three patients required retreatment after drug-free intervals of 6 months, 8 months, and 10 years. There were no differences in clinical severity among patients with DUOX2 amorphic/amorphic, amorphic/hypomorphic, and hypomorphic/hypomorphic variants. Conclusions: These results broaden the mutational spectrum of DUOX2. Furthermore, our data imply that patients with multiple pathogenic DUOX2 variants typically exhibit transient CH without significant genotype-phenotype correlations. Most importantly, this study demonstrated for the first time that these patients are at risk of developing recurrent hypothyroidism after a long drug-free interval.

Functional variants in a TTTG microsatellite on 15q26.1 cause familial nonautoimmune thyroid abnormalities.

Insufficient thyroid hormone production in newborns is referred to as congenital hypothyroidism. Multinodular goiter (MNG), characterized by an enlarged thyroid gland with multiple nodules, is usually seen in adults and is recognized as a separate disorder from congenital hypothyroidism. Here we performed a linkage analysis of a family with both nongoitrous congenital hypothyroidism and MNG and identified a signal at 15q26.1. Follow-up analyses with whole-genome sequencing and genetic screening in congenital hypothyroidism and MNG cohorts showed that changes in a noncoding TTTG microsatellite on 15q26.1 were frequently observed in congenital hypothyroidism (137 in 989) and MNG (3 in 33) compared with controls (3 in 38,722). Characterization of the noncoding variants with epigenomic data and in vitro experiments suggested that the microsatellite is located in a thyroid-specific transcriptional repressor, and its activity is disrupted by the variants. Collectively, we presented genetic evidence linking nongoitrous congenital hypothyroidism and MNG, providing unique insights into thyroid abnormalities.

Efficacy and safety of GH treatment in Japanese children with short stature due to SHOX deficiency: a randomized phase 3 study.

We conducted a randomized phase 3 study to investigate the efficacy and safety of GH treatment in prepubertal Japanese patients with short stature due to SHOX deficiency. The patients were randomly allocated to the GH-GH group (n = 10), in which the patients were treated with GH (0.35 mg/kg/wk) subcutaneously once daily for 24 mo, or the no-treatment (NT)-GH group (n = 9), in which the patients were untreated for the first 12 mo and then administered the same dosage of GH for the next 12 mo. At month 12, the ∆height standard deviation score (SDS) for chronological age (CA) and serum IGF-1 level were significantly higher in the GH-GH group than those in the NT-GH group. In contrast, bone age (BA) and ΔBA/ΔCA were numerically higher in the GH-GH group but were not statistically significant. At month 24, these parameters were comparable between the two groups. The height velocity was significantly larger in the GH-GH group during the first year and in the NT-GH group during the second year. No serious adverse drug reactions were observed; however, one patient in the GH-GH group exhibited increased insulin resistance at month 24. These results indicated that GH is a promising treatment option for short stature in patients with SHOX deficiency.

Expression levels and DNA methylation profiles of the growth gene SHOX in cartilage tissues and chondrocytes.

All attempts to identify male-specific growth genes in humans have failed. This study aimed to clarify why men are taller than women. Microarray-based transcriptome analysis of the cartilage tissues of four adults and chondrocytes of 12 children showed that the median expression levels of SHOX, a growth gene in the pseudoautosomal region (PAR), were higher in male samples than in female samples. Male-dominant SHOX expression was confirmed by quantitative RT-PCR for 36 cartilage samples. Reduced representation bisulfite sequencing of four cartilage samples revealed sex-biased DNA methylation in the SHOX-flanking regions, and pyrosequencing of 22 cartilage samples confirmed male-dominant DNA methylation at the CpG sites in the SHOX upstream region and exon 6a. DNA methylation indexes of these regions were positively correlated with SHOX expression levels. These results, together with prior findings that PAR genes often exhibit male-dominant expression, imply that the relatively low SHOX expression in female cartilage tissues reflects the partial spread of X chromosome inactivation into PAR. Altogether, this study provides the first indication that sex differences in height are ascribed, at least in part, to the sex-dependent epigenetic regulation of SHOX. Our findings deserve further validation.

Systematic molecular analyses for 115 karyotypically normal men with isolated non-obstructive azoospermia.

STUDY QUESTION: Do copy-number variations (CNVs) in the azoospermia factor (AZF) regions and monogenic mutations play a major role in the development of isolated (non-syndromic) non-obstructive azoospermia (NOA) in Japanese men with a normal 46, XY karyotype? SUMMARY ANSWER: Deleterious CNVs in the AZF regions and damaging sequence variants in eight genes likely constitute at least 8% and approximately 8% of the genetic causes, respectively, while variants in other genes play only a minor role. WHAT IS KNOWN ALREADY: Sex chromosomal abnormalities, AZF-linked microdeletions, and monogenic mutations have been implicated in isolated NOA. More than 160 genes have been reported as causative/susceptibility/candidate genes for NOA. STUDY DESIGN, SIZE, DURATION: Systematic molecular analyses were conducted for 115 patients with isolated NOA and a normal 46, XY karyotype, who visited our hospital between 2017 and 2021. PARTICIPANTS/MATERIALS, SETTING, METHODS: We studied 115 unrelated Japanese patients. AZF-linked CNVs were examined using sequence-tagged PCR and multiplex ligation-dependent probe amplification, and nucleotide variants were screened using whole exome sequencing (WES). An optimized sequence kernel association test (SKAT-O), a gene-based association study using WES data, was performed to identify novel disease-associated genes in the genome. The results were compared to those of previous studies and our in-house control data. MAIN RESULTS AND THE ROLE OF CHANCE: Thirteen types of AZF-linked CNVs, including the hitherto unreported gr/gr triplication and partial AZFb deletion, were identified in 63 (54.8%) cases. When the gr/gr deletion, a common polymorphism in Japan, was excluded from data analyses, the total frequency of CNVs was 23/75 (30.7%). This frequency is higher than that of the reference data in Japan and China (11.1% and 14.7%, respectively). Known NOA-causative AZF-linked CNVs were found in nine (7.8%) cases. Rare damaging variants in known causative genes (DMRT1, PLK4, SYCP2, TEX11, and USP26) and hemizygous/multiple-heterozygous damaging variants in known spermatogenesis-associated genes (TAF7L, DNAH2, and DNAH17) were identified in nine cases (7.8% in total). Some patients carried rare damaging variants in multiple genes. SKAT-O detected no genes whose rare damaging variants were significantly accumulated in the patient group. LIMITATIONS, REASONS FOR CAUTION: The number of participants was relatively small, and the clinical information of each patient was fragmentary. Moreover, the pathogenicity of identified variants was assessed only by in silico analyses. WIDER IMPLICATIONS OF THE FINDINGS: This study showed that various AZF-linked CNVs are present in more than half of Japanese NOA patients. These results broadened the structural variations of AZF-linked CNVs, which should be considered for the molecular diagnosis of spermatogenic failure. Furthermore, the results of this study highlight the etiological heterogeneity and possible oligogenicity of isolated NOA. STUDY FUNDING/COMPETING INTEREST(S): This study was supported by Grants from the Japan Society for the Promotion of Science (21K19283 and 21H0246), the Japan Agency for Medical Research and Development (22ek0109464h0003), the National Center for Child Health and Development, the Canon Foundation, the Japan Endocrine Society, and the Takeda Science Foundation. The results of this study were based on samples and patient data obtained from the International Center for Reproductive Medicine, Dokkyo Medical University Saitama Medical Center, Koshigaya, Japan. The authors have no conflicts of interest to disclose. TRIAL REGISTRATION NUMBER: N/A.

A novel variant of IGSF1 in siblings with congenital central hypothyroidism whose diagnosis was prompted by school health checkups.

Following the partial revision of the enforcement regulations of the School Health and Safety Act, school health checkups incorporated growth evaluation of schoolchildren in April 2016 using growth charts. We report cases of congenital central hypothyroidism (C-CH) in siblings with a novel nonsense variant in the immunoglobulin superfamily member 1 gene (IGSF1); their diagnoses were prompted by school health checkups. School checkups revealed that the older brother was overweight and had a reduced growth rate at the age of 11 yr, whereas the younger brother was overweight and had short stature at the age of 8 yr. They were diagnosed with C-CH because of normal thyroid-stimulating hormone (TSH) levels despite a low free thyroxine level and low TSH response in the thyrotropin-releasing hormone stress test. Only the older brother had prolactin deficiency and testicular growth without elevated testosterone levels. The siblings harbored a novel nonsense variant in exon 16 of IGSF1 (NM_001555.5: c.3056G>A: p.Trp1019Ter) and were diagnosed with IGSF1 deficiency. In Japan, C-CH may be overlooked because TSH-based newborn screening alone is usually performed for patients with congenital hypothyroidism. The implementation of growth monitoring using growth charts in school health checkups may prompt new C-CH diagnoses.

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.

(Epi)genetic and clinical characteristics in 84 patients with pseudohypoparathyroidism type 1B.

OBJECTIVE: Pseudohypoparathyroidism type 1B (PHP1B) caused by methylation defects of differentially methylated regions (DMRs) on the GNAS locus can be categorized into groups according to etiologies and methylation defect patterns of the DMRs. The aim of this study was to clarify the clinical characteristics of each group. DESIGN: Comprehensive molecular analyses consisting of methylation, copy number, and microsatellite analyses. METHODS: Eighty-four patients with PHP1B were included in this study. We classified them into 5 groups, namely, autosomal dominant inheritance-PHP1B (Group 1, G1), sporadic-PHP1B (G2), and atypical-PHP1B (G3-G5), based on the methylation defect patterns in 4 DMRs on the GNAS locus and etiologies and evaluated the clinical findings in each group and compared them among the groups. RESULTS: G2 had the youngest age and the highest serum intact parathyroid hormone levels among the 5 groups at the time of diagnosis. The most common symptoms at the time of diagnosis were tetany in G1, and seizures or loss of consciousness in G2. Albright's hereditary osteodystrophy and PHP-suggestive features were most frequently observed in the G2 proband. Nine patients had neurodevelopmental disorders (NDs) consisting of mild to borderline intellectual disability and/or developmental delay. There were no significant correlations between the average methylation ratios of 7 CpG sites in the GNAS-A/B:TSS-DMR and hormonal and biochemical findings. CONCLUSION: This study revealed the differences in some clinical characteristics, particularly clinical features, and ages at the time of diagnosis between G2 and other groups and detailed NDs observed in some patients with PHP1B.

Variations in gender identity and sexual orientation of university students.

Background: Previous studies have shown that a small percentage of people in the general population have atypical gender identity and/or sexual orientation. Aim: This study aimed to explore variations in gender identity and sexual orientation in university students and determine genetic factors associated with these variations. Methods: Deviations from complete gender congruence and exclusive heterosexual orientation in 736 Japanese university students were quantitatively assessed with self-assessment questionnaires. Next, we conducted genetic tests for 80 participants who showed relatively low gender identity scores and/or atypical sexual orientation. These genetic tests consisted of repeat number analysis of the androgen receptor gene (AR) and a SKAT-O: an optimal unified sequence kernel association test, which is an exome-based rare variant association study. The results of the genetic tests were compared with the Japanese reference data and the results of our 637 control samples. Outcomes: We calculated the gender identity and sexual orientation scores of all participants and analyzed the molecular data of 80 selected participants. Results: The gender identity scores of 736 participants were broadly distributed: only ~15% of natal males and ~5% of natal females had the maximum score that corresponds to complete gender congruence. The sexual orientation scores also varied: ~80% of natal males and ~60% of natal females showed exclusive heterosexual orientation. We found no association between gender characteristics and AR repeat numbers. The SKAT-O showed that rare damaging variants of TDRP and 3 other genes were more common in the 80 participants than in the control group. Clinical Implications: Our data support the view that gender is a phenotypic continuum rather than a binary trait. Strength and Limitations: This study quantitatively assessed the gender characteristics of a large cohort of university students. Moreover, we conducted systematic screening for genetic factors associated with gender variations. The weaknesses of the study were the limited analytic power of the questionnaires, the relatively small sample for molecular analyses, and incomplete clinical information and relatively advanced ages of the control group. Conclusion: This study revealed significant variations in gender identity and sexual orientation in university students, which may be partly associated with variants in TDRP or other genes.

The Structural Abnormalities Are Deeply Involved in the Cause of RPGRIP1-Related Retinal Dystrophy in Japanese Patients.

Leber congenital amaurosis (LCA) is the most severe form of inherited retinal dystrophy. RPGRIP1-related LCA accounts for 5-6% of LCA. We performed whole-exome sequencing and whole-genome sequencing (WGS) on 29 patients with clinically suspected LCA and examined ophthalmic findings in patients with biallelic pathogenic variants of RPGRIP1. In addition to five previously reported cases, we identified five cases from four families with compound heterozygous RPGRIP1 variants using WGS. Five patients had null variants comprising frameshift variants, an Alu insertion, and microdeletions. A previously reported 1339 bp deletion involving exon 18 was found in four cases, and the deletion was relatively prevalent in the Japanese population (allele frequency: 0.002). Microdeletions involving exon 1 were detected in four cases. In patients with RPGRIP1 variants, visual acuity remained low, ranging from light perception to 0.2, and showed no correlation with age. In optical coherence tomography images, the ellipsoid zone (EZ) length decreased with age in all but one case of unimpaired EZ. The retinal structure was relatively preserved in all cases; however, there were cases with great differences in visual function compared to their siblings and a 56-year-old patient who still had a faint EZ line. Structural abnormalities may be important genetic causes of RPGRIP1-related retinal dystrophy in Japanese patients, and WGS was useful for detecting them.