Next Article in Journal
Morphological Signatures of Neurogenesis and Neuronal Migration in Hypothalamic Vasopressinergic Magnocellular Nuclei of the Adult Rat
Next Article in Special Issue
Rare Germline Variants in the Adenomatous Polyposis Coli Gene Associated with Dental and Osseous Anomalies
Previous Article in Journal
TLR2 and TLR4 Polymorphisms Are Not Associated with Dental Caries in Polish Children
Previous Article in Special Issue
Dysfunction in IGF2R Pathway and Associated Perturbations in Autophagy and WNT Processes in Beckwith–Wiedemann Syndrome Cell Lines
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 13
10.3390/ijms25136986
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association Study with Three Control Cohorts of Japanese Patients with Esotropia and Exotropia of Comitant Strabismus and Idiopathic Superior Oblique Muscle Palsy

by
Toshihiko Matsuo
1,2,*,
Ichiro Hamasaki
2,
Yoichiro Kamatani
3,
Takahisa Kawaguchi
4,
Izumi Yamaguchi
4,
Fumihiko Matsuda
4,
Akira Saito
5,
Kazuyuki Nakazono
5 and
Shigeo Kamitsuji
5
1
Graduate School of Interdisciplinary Science and Engineering in Health Systems, Okayama University, Okayama City 700-8558, Japan
2
Department of Ophthalmology, Okayama University Hospital, Okayama City 700-8558, Japan
3
Department of Computational Biology and Medical Sciences, Graduate School of Frontier Sciences, The University of Tokyo, Tokyo 108-8639, Japan
4
Center for Genomic Medicine, Graduate School of Medicine, Kyoto University, Kyoto 606-8507, Japan
5
StaGen Co., Ltd., Tokyo 111-0051, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(13), 6986; https://doi.org/10.3390/ijms25136986
Submission received: 13 May 2024 / Revised: 10 June 2024 / Accepted: 15 June 2024 / Published: 26 June 2024
(This article belongs to the Collection Feature Papers in Molecular Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Esotropia and exotropia in the entity of comitant strabismus are multifactorial diseases with both genetic and environmental backgrounds. Idiopathic superior oblique muscle palsy, as the predominant entity of non-comitant (paralytic) strabismus, also has a genetic background, as evidenced by varying degrees of muscle hypoplasia. A genome-wide association study (GWAS) was conducted of 711 Japanese patients with esotropia (n= 253), exotropia (n = 356), and idiopathic superior oblique muscle palsy (n = 102). The genotypes of single nucleotide polymorphisms (SNPs) were determined by Infinium Asian Screening Array. Three control cohorts from the Japanese population were used: two cohorts from BioBank Japan (BBJ) and the Nagahama Cohort. BBJ (180K) was genotyped by a different array, Illumina Infinium OmniExpressExome or HumanOmniExpress, while BBJ (ASA) and the Nagahama Cohort were genotyped by the same Asian array. After quality control of SNPs and individuals, common SNPs between the case cohort and the control cohort were chosen in the condition of genotyping by different arrays, while all SNPs genotyped by the same array were used for SNP imputation. The SNPs imputed with R-square values ≥ 0.3 were used to compare the case cohort of each entity or the combined entity with the control cohort. In comparison with BBJ (180K), the esotropia group and the exotropia group showed CDCA7 and HLA-F, respectively, as candidate genes at a significant level of p < 5 × 10−8, while the idiopathic superior oblique muscle palsy group showed DAB1 as a candidate gene which is involved in neuronal migration. DAB1 was also detected as a candidate in comparison with BBJ (ASA) and the Nagahama Cohort at a weak level of significance of p < 1 × 10−6. In comparison with BBJ (180K), RARB (retinoic acid receptor-β) was detected as a candidate at a significant level of p < 5 × 10−8 in the combined group of esotropia, exotropia, and idiopathic superior oblique muscle palsy. In conclusion, a series of GWASs with three different control cohorts would be an effective method with which to search for candidate genes for multifactorial diseases such as strabismus.

1. Introduction

Both eyes are located at the front of the head in humans and some animals, and this results in the overlapped visual fields of two eyes, and thus generates binocular vision, which consists of three levels: simultaneous perception, binocular fusion, and stereopsis. Strabismus is a disease of binocular misalignment in which both eyes do not align with the same visual target, leading to abnormal binocular vision, such as monocular suppression, double vision, and a lack of stereopsis. Two entities of strabismus are comitant strabismus and non-comitant (paralytic) strabismus. Comitant strabismus shows the same degree of eye deviation in all directions of the gaze, while non-comitant strabismus shows different degrees of deviation according to the different directions of the gaze. Major types of comitant strabismus are esotropia and exotropia, which present as inward and outward deviations, respectively, of one eye relative to the other eye along the horizontal axis [1]. Major types of non-comitant strabismus involve the paralysis of ocular motor nerves, including the abducens, oculomotor, and trochlear nerves, which innervate the extraocular muscles [2].
Esotropia and exotropia, as forms of comitant strabismus, are fundamentally early-onset in life [3,4,5,6,7] and are multifactorial diseases which have both genetic and environmental background factors [8,9,10]. The prevalence of comitant strabismus is estimated at roughly 1–2% of the population [4,5]. The genetic background for comitant strabismus is supported by the presence of family history [11,12] and the concordance of phenotypes in monozygotic twins [13,14], while an environmental background is suggested by perinatal hypoxia [15,16]. Idiopathic superior oblique muscle palsy, as a form of non-comitant strabismus, manifests in the earlier years of life as congenital palsy or later in life as the decompensation of binocular vision in the process of aging [2,17,18]. The prevalence of idiopathic superior oblique muscle palsy remains to be elucidated and is empirically considered to be lower by an order of magnitude than comitant strabismus. A genetic background is highly suspected in idiopathic superior oblique muscle palsy since varying degrees of hypoplasia of the superior oblique muscle have been disclosed by orbital magnetic resonance imaging [19,20,21,22,23,24].
Linkage analyses by chromosomal mapping to find susceptibility loci for esotropia and exotropia of comitant strabismus have been conducted in families from different populations [25,26,27,28,29,30,31,32]. Linkage analysis has also been performed in small Japanese families with idiopathic superior oblique muscle palsy [2,23]. In contrast, only two reports have described genome-wide association studies of Caucasian populations with esotropia [33] and self-reported strabismus of an unknown entity [34], respectively. These studies showed GET1 (WRB) as a susceptibility gene for esotropia [33] and TSPAN10 as one for self-reported strabismus [34]. In the present research, we conducted a genome-wide association study (GWAS) of Japanese patients with esotropia and exotropia of comitant strabismus as well as idiopathic superior oblique muscle palsy in comparison with three different Japanese cohorts: BBJ (180K) and BBJ (ASA) as two cohorts of BioBank Japan (BBJ), and also the Nagahama Cohort.

2. Results

2.1. Eligible Case Samples

A call rate equal to or greater than 0.99 was obtained in 711 of 723 case samples, which passed the quality check for genomic DNA. SNPs which did not fit for analysis were then deleted based on the following parameters and cutoff values: the mean normalized intensity (R Mean) (AA, AB, and BB) < 0.25, where AB indicated the heterozygote cluster, and cluster separation (Cluster Sep) < 0.4 for SNPs on the Y chromosome; or the R Mean (AA, AB, and BB) < 0.25 and Cluster Sep < 0.4, and call frequency (Call Freq) < 0.99 on the other chromosomes except for the Y chromosome. The number of SNPs determined was 657,060 on the Infinium Asian Screening Array-24 v.1.0 BeadChip.
Of the 711 patients with quality checks, 253 patients had esotropia, 356 had exotropia, and 102 had idiopathic superior oblique muscle palsy. Among the 253 patients with esotropia, infantile esotropia was diagnosed in 84 patients, accommodative esotropia including partially accommodative esotropia in 24, acute-onset (acquired) esotropia in 6, and esotropia (unclassified) in 139. In 356 patients with exotropia, intermittent exotropia was diagnosed in 131 patients, and constant exotropia in 225 patients.
In identity-by-descent (IBD) estimation in the 711 patients (Figure 1), 20 pairs (8 pairs with esotropia, 9 pairs with exotropia, and 3 pairs with idiopathic superior oblique muscle palsy) showed PI_HAT = 0.5, while 19 pairs (3 pairs with esotropia, 12 pairs with exotropia, 2 pairs with esotropia in one individual and exotropia in the other, and 2 pairs with idiopathic superior oblique muscle palsy) showed PI_HAT = 1. After the exclusion of patients with PI_HAT > 0.1, 666 case samples in a GWAS with BBJ (180K) (Figure 2) or BBJ (ASA) (Figure 3) consisted of 239 patients with esotropia, 331 patients with exotropia, and 96 patients with idiopathic superior oblique muscle palsy (Figure 1). In a GWAS with the Nagahama Cohort, the number of case samples was 662 after the different flow of exclusion (Figure 4).

2.2. Imputed SNP Dosage in Three Different Control Cohorts

The imputed SNP dosages in three control cohorts were compared with the SNP dosage in the Japanese population (1KG JPT) of the 1000 Genomes Project to show no deviation in allele frequencies (Figure 5). The imputed SNP dosages were also compared among the three case–control groups to show no deviation in allele frequencies (Figure 6).

2.3. GWAS with Three Different Control Cohorts

Manhattan plots were drawn with chromosomal location on the horizontal axis and statistical value, −log10(p-value), on the vertical axis for esotropia (Figure 7A), exotropia (Figure 7B), idiopathic superior oblique muscle palsy (Figure 7C), the combination of esotropia and exotropia (Figure 7D), and all entities of esotropia, exotropia, and superior oblique muscle palsy (Figure 7E). In three series of GWASs with three different control cohorts, meaningful peaks with a loose standard of p < 1 × 10−6 formed by the gathering of proximal SNPs in the linkage are indicated by red circles in each plot. Q–Q (log quantile–quantile p-value) plots which were depicted with the expected −log10(p) on the horizontal axis versus the observed −log10(p) on the vertical axis are shown next to each Manhattan plot. Only the group of exotropia was suggested to have a structured bias from the Q–Q plots.
As common peaks in three series of GWASs with three different control cohorts, three loci were suggested for esotropia (SFPQ and HLX-AS1 in chromosome 1, and MIR4471 in chromosome 8 as candidate genes; Figure 8A–C, Table 1) and two loci for idiopathic superior oblique muscle palsy (DAB1 in chromosome 1 and PELO in chromosome 5 as candidate genes; Figure 8D,E, Table 1).
At the standard significance of p < 5 × 10−8 (Table 2, Supplementary Figure S1), CDCA7 in chromosome 2 for esotropia, HLA-F in chromosome 6 for exotropia, and DAB1-AS1 in chromosome 1 for idiopathic superior oblique muscle palsy were suggested as candidates in a GWAS with BBJ (180K). Additionally, in a GWAS with BBJ (180K), ACAP2 in chromosome 3 and HLA-F in chromosome 6 were suggested as candidates in the combined group of esotropia and exotropia. In a GWAS with BBJ (180K), RARB in chromosome 3 and HLA-F in chromosome 6 were suggested as candidates in the combined three-entities group with esotropia, exotropia, and idiopathic superior oblique muscle palsy. In contrast, no peak with p < 5 × 10−8 was found in a GWAS with BBJ (ASA) or the Nagahama Cohort.

3. Discussion

The long-term goal of our studies up until now is to search for strabismus susceptibility genes. In our initial attempt, we used linkage analyses in small-sized Japanese families with comitant strabismus including esotropia and exotropia [28,29]. In our previous study of linkage analyses, we have reached two potential candidate genes for strabismus susceptibility, MGST2 and WNT2 in chromosomal loci 4q28.3 and 7q31.2, respectively [32]. We then proceeded to make MGST2-deficient mice by CRISPR-Cas9 technology and proved that eyeballs in the homozygous mice were in an ovoid shape and significantly larger in size and volume than the wild-type eyeballs, by measurements with magnetic resonance imaging [35]. The changes in the shape and size of the eyeballs might result in a different torque of the extraocular muscles which attach to the surface of the eyeballs, and, thus, might lead to the different alignment of both eyes.
Since the alignment of both eyes would be influenced by the innervation of the extraocular muscles with ocular motor nerves such as the oculomotor, trochlear, and abducens nerves from the brainstem, candidate genes for strabismus susceptibility might be expressed in the brain and have an effect on the function of the brain. Furthermore, binocularity such as stereopsis belongs to a higher function of the central nervous system to integrate visual sensory input from both eyes. In this context, we had also carried out whole-genome sequencing in three small families with idiopathic superior oblique muscle palsy as a major form of non-comitant strabismus [2]. In that previous study, we combined the two methods of linkage analysis and SNP sorting, and we reached several candidate genes for idiopathic superior oblique muscle palsy [2].
Based on the rationale that a different approach in genetic statistics would help reach candidate genes for strabismus susceptibility, we performed a series of the genome-wide association study (GWAS) for esotropia and exotropia of comitant strabismus as well as idiopathic superior oblique muscle palsy in the entity of non-comitant strabismus in the present research. A critical point in the GWAS is what kind of a genetic cohort should be used as a control. We first chose the preexisting BBJ cohort of the Japanese population which was genotyped by a different array from that used in strabismic patients in this study. In the condition that different arrays are used separately to genotype the case and control, the set of loaded SNPs in each array was not exactly the same. Therefore, the common SNPs between the different arrays have to be chosen and the number of SNPs used for a comparison of the case with the control is reduced to a large extent. SNP imputation based on the smaller number of common SNPs would be less accurate. The cohort in the Tohoku Medical Megabank (ToMMo) which was typed by the same Asian array was not used in this study since the population in the northern part of Japan is somewhat different from that in the western part of Japan where the patients in this study resided. Therefore, we chose the cohort of the Nagahama Study with participants living in western Japan, although the number is small. Meanwhile, BBJ released a new series of a cohort of the Japanese population which was genotyped by the Asian array.
In the GWAS, related individuals are excluded from the case population and control population, based on the identity-by-descent (IBD) analysis. In an attempt to test different conditions in statistical analyses, we initially excluded only one of the pairs with IBD = 1.0 in the case cohorts and retained pairs with IBD = 0.5 to perform the GWAS since the case samples were collected not only from independent individuals but also from family members with strabismus for the linkage analysis [28,29,32]. We also performed a GWAS with the case samples with the exclusion of IBD > 0.1 in the usual setting. The GWAS results were basically the same whether the case samples did or did not contain pairs with IBD = 0.5.
In the case groups, patients with esotropia, exotropia, or idiopathic superior oblique muscle palsy were defined as each independent group to perform a GWAS with three different control cohorts. In addition, we defined, as the case group, the combined group of patients with esotropia and exotropia of comitant strabismus, and also defined the all-combined samples with esotropia, exotropia, and idiopathic superior oblique muscle palsy. A rationale for combining the different entities of comitant and non-comitant strabismus is that these different phenotypes are indeed present in a pedigree [28,29,36]. A merit of combining patients with different phenotypes would be that the number of patients in a case group becomes larger, large enough to detect a significant locus which would serve as a common factor for comitant and non-comitant strabismus with a genetic background. In contrast, a demerit would be that the combined case group with heterogeneous phenotypes would dilute the significant locus which would be specific for each group with different phenotypes.
The Nagahama Cohort genotyped by the same Asian array which was used for the genotyping of the case group had too small a number of individuals to detect significant loci. The GWAS using BBJ (ASA) as a control cohort could not detect significant loci even though this control cohort had a large number of individuals and shared the same typing array with the case group in this study. In contrast, the GWAS with the largest number of control individuals in BBJ (180K) in this study did reach strabismus susceptibility loci at the significant level of p = 5 × 10−8. Since the typing arrays were different between the case group and the control group of BBJ (180K), common SNPs were chosen first to proceed to SNP imputation. Therefore, the SNP imputation based on the smaller number of SNPs in the GWAS with BBJ (180K) would be less accurate than the GWAS with BBJ (ASA). Under the circumstances, we tried to find common significant loci which were shared by three series of GWASs with three different sets of control cohorts. In order not to miss any interested loci, the significant level was loosened to p = 1 × 10−6 from the standard level of p = 5 × 10−8 in this research attempt. This approach would be expected to reinforce the findings in a single GWAS with BBJ (180K) which was genotyped by the different array.
The loci for strabismus susceptibility which were detected as significant in this GWAS did not overlap with the loci which were reported in previous linkage analyses [25,26,27,28,29,30,31,32,36] and a preceding GWAS in other populations [33,34]. MGST2 (microsomal glutathione S-transferase 2 at 4q31.1) was detected as a susceptibility gene for comitant strabismus of esotropia and exotropia by our previous linkage study [32] but was not replicated by the present GWAS. WRB (tryptophan-rich basic protein, also known as GET1, guided entry of tail-anchored proteins factor 1, at 21q22.2) and TSPAN10 (tetraspanin 10 at 17q25.3), which were detected as a susceptibility gene for esotropia in the U.S.A. [33] and for self-reported strabismus of unknown entity in the U.K. [34], respectively, were also not replicated by the present GWAS (see regional plots in Supplementary Figure S2). Since strabismus is a multifactorial disease, the loci detected in this GWAS and the other loci detected by the preceding linkage studies and GWAS would play a role in the development of strabismus like polygenetic risk scores which have been described in GWASs for other diseases such as diabetes mellitus [37,38,39,40]. As strabismus susceptibility genes, this GWAS suggests non-coding RNAs, which have been also suggested by the other studies [41,42].
Notably, DAB1 was detected as a susceptibility gene at a significant level of p < 5 × 10−8 for idiopathic superior oblique muscle palsy in the GWAS with BBJ (180K), and was also detected as such at a weak level of significance of p < 1 × 10−6 in the GWAS with BBJ (ASA) and the Nagahama Cohort (Table 1 and Table 2). DAB1 is involved in neuronal migration and is a responsible gene for spinocerebellar ataxia type 37 [43]. Trochlear nerve hypoplasia is considered to be a cause for idiopathic superior oblique muscle palsy which manifests as varying levels of the hypoplasia of the superior oblique muscle [2,17,19,20]. Under the circumstances, aberrant neuronal migration would serve as the background for the trochlear nerve hypoplasia [41]. In comparison with the control cohort, BBJ (180K), RARB (retinoic acid receptor-β) was detected as a candidate at a significant level of p < 5 × 10−8 in the combined group of esotropia, exotropia, and idiopathic superior oblique muscle palsy (Table 2). Since retinoic acids are involved in neural crest cell migration in eye morphogenesis [44,45,46], retinoic acids would have a role in the development of extraocular muscles of neural crest cell origin [47,48].
A major limitation in this study was that the three control cohorts used in this study would contain unidentified individuals with strabismus, although in a small number. This limitation would be cancelled by the large number of individuals in the control cohorts, compared with the case cohort. Another limitation is the fact that clinical manifestations as phenotypes vary largely from individual to individual in each entity of esotropia, exotropia, and idiopathic superior oblique muscle palsy. For instance, the degree of deviation is different and the level of stereopsis is different from individual to individual. The diagnosis of each entity in the case samples was, of course, made by ophthalmologists with a subspecialty in strabismus, based on the diagnostic criteria. Interestingly, a candidate gene, DAB1, was detected reproducibly in idiopathic superior oblique muscle palsy by the series of GWASs with all three control cohorts. The results suggest the entity of idiopathic superior oblique muscle palsy would be more homogeneous, compared with the entities of esotropia and exotropia, from the viewpoint of the diagnostic criteria.
In conclusion, this study is the first to conduct a GWAS for comitant and non-comitant (paralytic) strabismus in the Japanese population. To the best of our knowledge, this study is the first, especially, to conduct a GWAS for idiopathic superior oblique muscle palsy. The series of GWASs with three different control cohorts would be an effective method to search for candidate genes for multifactorial diseases such as strabismus. The susceptibility genes detected in this study remain candidates, and, thus, future studies will be required for the functional validation of these candidate genes in comitant and non-comitant strabismus.

4. Materials and Methods

4.1. Participants

Participants in this study were patients diagnosed with esotropia or exotropia of comitant strabismus, or idiopathic superior oblique muscle palsy at Okayama University Hospital (Figure 1). Genomic DNA was isolated from leukocytes of 10 mL peripheral blood in each patient after written consent to be involved in the study. The study conformed to the tenets of the Declaration of Helsinki and was approved by the Ethics Committee of Okayama University Graduate School of Medicine, Dentistry, and Pharmaceutical Sciences and Okayama University Hospital (Identifier: 1512-027, date of approval on updated version: 14 August 2020).

4.2. Genotyping of Participants

Genome-wide SNP genotyping was carried out with Infinium Asian Screening Array-24 v.1.0 BeadChip Kit (Illumina, San Diego, CA, USA) at Riken Genesis, Co., Ltd. (Tokyo, Japan). The quality of genomic DNA was assessed by agarose gel electrophoresis and the PicoGreen assay (Thermo Fisher Scientific, Waltham, MA, USA) using a plate reader (SpectraMax, Molecular Devices, San Hose, CA, USA). Whole genome of each sample which passed the quality check was alkali-denatured, neutralized, and amplified enzymatically by Multi-Sample Amplification Master Mix (MSM, Illumina) overnight in an incubation oven (Illumina). The amplified DNA was enzymatically fragmented by Fragmentation Solution (FMS, Illumina) for one hour in a Hybex microsample incubator (SciGene, Sunnyvale, CA, USA). The fragmented DNA was precipitated with 2-propanol by centrifugation, dissolved in hybridization buffer, and denatured at 48 °C overnight in a hybridization oven (Illumina). The fragmented and denatured DNA in a microsample incubator was kept at 95 °C for 20 min, and, then, was dispensed with Tecan system (Maennedorf, Switzerland) to the BeadChip, and hybridized at 48 °C overnight in the hybridization oven. Unhybridized free DNA on the BeadChip was washed away, and primers which hybridized with the sample DNA were extended with labelled nucleotides by Extension Mix Long (EML, Illumina). The BeadChip was coated with an anti-decay reagent, dried in a vacuum dessicator, and placed on a carrier of iScan system (Illumina) for scan.
The imaging data were analyzed by genotyping analysis software GenomeStudio 2011.1 (Illumina) to calculate the call rate of each sample. Data were built at 4 steps in GenomeStudio: at step 1 to use standard cluster file before quality check, at step 2 to perform quality check to delete samples with low call rates, at step 3 to perform cluster optimization, and at step 4 to perform SNP quality check by parameter-setting to delete SNPs which did not fit for analysis. Clustering Algorithm was set at GenTrain 2.0 and GenCall Threshold was set at 0.15 in GenomeStudio. Genotyping data output was carried out by the software GenomeStudio Framework version 2013 and Genotyping Module version 2.0.3.

4.3. Control Samples

Three Japanese cohorts were used as control samples in genome-wide association study, and were designated in this study as BBJ (180K), BBJ (ASA), and Nagahama. In the first cohort (Figure 1), BBJ (180K), the SNP data of 217,038 samples in BioBank Japan which were registered as 47 target common diseases [49] were obtained from the National Bioscience Database Center (NBDC) after ethical approval (Identifier: 0016, date of approval: 15 February 2019) by the Institutional Review Board affiliated with Tsukuba International Clinical Pharmacology Clinic (Tsukuba City, Japan) for the use as control, based on the Ethical Guidelines for Medical and Health Research Involving Human Subjects, issued by the Government of Japan. These BioBank Japan data were genotyped by Illumina Infinium OmniExpressExome-8 v1.0 or v1.2 or Illumina HumanOmniExpress-12 v1.0. Samples with call rate less than 95% were excluded to obtain dataset of 182,476 samples for analysis.
In the second cohort (Figure 1), BBJ (ASA), genome-wide SNP data of 54,405 samples were obtained from NBDC after ethical approval (Identifier: 0026, date of approval: 28 March 2022) by the Institutional Review Board affiliated with Tsukuba International Clinical Pharmacology Clinic (Tsukuba City, Japan) for the use as control. The genotype data were determined by Infinium Asian Screening Array-24 v1.0 BeadChip Kit in 11,716 individuals from BBJ first cohort and 42,689 individuals from BBJ second cohort (Dataset ID JGAS000412, Release Date: 30 November 2021). The controlled access of the dataset was approved as Data in Use ID JGAD000529 (Principal Investigator: Shigeo Kamitsuji, Affiliation: Statistical Analysis Division, Stagen Co., Ltd., Country/Region: Japan, Research Title: Genome-Wide Association Study to identify genetic factors for strabismus in Japanese population, Period of Data Use: 13 February 2023–28 February 2027). Samples with call rate less than 95% were excluded to obtain dataset of 53,409 samples for analysis.
In the third cohort, Nagahama, the genotype data were determined by Infinium Asian Screening Array-24 v1.0 BeadChip Kit in 3570 individuals who resided in Nagahama City, Shiga Prefecture in western Japan. The project is a collaboration between Kyoto University and Nagahama City [50] and this study was approved by Ethics Committee in Kyoto University (Identifier: G1289, date of approval: 26 November 2020).

4.4. Whole-Genome Imputation

Quality control of SNPs was carried out by allele frequency and R-square value. In the case group, expected minimum allele frequency was defined as 1/(the number of cases × 2). In the control group, the minimum allele frequency was evaluated by 95% confidence interval of alternative allele frequency in the 1000 Genomes Project JPT (Japanese Population in Tokyo) [51,52]. The fitness with Hardy–Weinberg equilibrium principle was evaluated by R-square values ≥ 0.3 which were used for quality control in SNP imputation. Cryptic relatedness of samples to each other was estimated by identity-by-descent (IBD) PI_HAT values in plink [53].
Based on quality control to exclude SNPs with call rate less than 99%, minor allele frequency less than 1%, and p-value of Hardy–Weinberg equilibrium test less than 1 × 10−6, 459,558 of 657,060 SNPs were chosen in 711 samples of comitant strabismus or idiopathic superior oblique muscle palsy (Figure 2). According to the same standard, 495,375 of 779,476 SNPs were chosen in 182,476 samples of the BBJ (180K). Based on the reference panels for Japanese population [51,52], imputation was carried out by the software Beagle v5.1 and plink2 [53] in 1269 of 1304 chunks with 2 Mbp length for the total 127,315 SNPs which were common between the case and the control [54,55,56,57]. The remaining chunks could not be imputed because of no SNP data on those chunks. The data of 17,792,106 SNPs were used for genome-wide association study after quality control of 47,109,465 imputed SNPs by imputation quality metric R-square (estimated r2, specific to each SNP) equal to or greater than 0.3 and minor allele frequency greater than zero (Figure 2).
In the case of BBJ (ASA) as control samples, 435,851 of 657,060 SNPs were chosen in 53,409 samples of BBJ (ASA), according to the same standard as described above. Based on the reference panels for Japanese population [51,52], imputation was carried out by a method of pre-phasing imputation with Minimac4 for 415,352 SNPs which were common between the case and the control [58,59,60]. The data of 17,883,408 SNPs were used for genome-wide association study after quality control of 47,109,465 imputed SNPs by imputation quality metric R-square (estimated r2, specific to each SNP) equal to or greater than 0.3. In GWAS, the number of SNPs was 12,597,304 after the final exclusion of 5,286,104 SNPs with minor allele frequency = 0 (Figure 3).
In the case of Nagahama Cohort as control samples, SNP imputation was carried out as shown in the flow chart (Figure 4).

4.5. Genome-Wide Association Study

As SNP dosage, the genotype value was defined as an expected value (continuous value from 0 to 2) for having the number of the alternative allele against the reference allele on the UCSC hg19 Genome Browser. Odds ratio for each SNP was estimated by one increment in SNP dosage in the case group, compared with the control group. The association was evaluated by Wald test and the distribution of p-values were depicted as a Manhattan plot. In logistic regression analysis, case group or control group was defined as an explanatory variable and SNP dosage as a response variable. When regression coefficient (BETA) is a positive value, the case group was interpreted to have the alternative allele (target allele in population) at a higher rate compared with the control group. Odds ratio estimated by one increment in SNP dosage in the case group, compared with the control group, could be obtained by calculation of eBETA [61,62]. In addition to the program, plink2 [53], regenie [61] was used for calculation to reduce other biases in the combined analysis of genotypes and covariates. The chromosomal areas in interest were depicted by the software Locus Zoom js v0.12 [63].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25136986/s1.

Author Contributions

T.M. made the diagnosis, collected genomic DNA, and designed and conducted the study. I.H. made the diagnosis and collected genomic DNA. Y.K. advised the study design and conducted the GWAS with the Nagahama Cohort. T.K., I.Y. and F.M. maintained the Nagahama Cohort. A.S., K.N. and S.K. conducted the GWAS with BioBank Japan. Conceptualization, T.M., Y.K. and S.K.; methodology, Y.K., A.S., K.N. and S.K.; software, A.S., K.N. and S.K.; validation, T.M., Y.K. and K.N.; formal analysis, T.M.; investigation, K.N.; resources, T.M. and I.H.; data curation, K.N.; writing—original draft preparation, T.M.; writing—review and editing, Y.K. and S.K.; visualization, T.M. and K.N.; supervision, T.M. and S.K.; project administration, T.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of Okayama University Graduate School of Medicine, Dentistry, and Pharmaceutical Sciences and Okayama University Hospital (Identifier: 1512-027, date of approval on updated version: 14 August 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Additional data are available upon reasonable request from the corresponding author.

Conflicts of Interest

Authors Akira Saito, Kazuyuki Nakazono and Shigeo Kamitsuji were employed by the company StaGen Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Zhang, J.; Matsuo, T.; Hamasaki, I.; Sato, K. Whole exome-sequencing of pooled genomic DNA samples to detect quantitative trait loci in esotropia and exotropia of strabismus in Japanese. Life 2022, 12, 41. [Google Scholar] [CrossRef] [PubMed]
  2. Matsuo, T.; Chaomulige, M.M.; Hosoya, O.; Saito, A.; Nakazono, K. Candidate genes in testing strategies for linkage analysis and bioinformatic sorting of whole genome sequencing data in three small Japanese families with idiopathic superior oblique muscle palsy. Int. J. Mol. Sci. 2022, 23, 8626. [Google Scholar] [CrossRef]
  3. Matsuo, T.; Matsuo, C.; Matsuoka, H.; Kio, K. Detection of strabismus and amblyopia in 1.5- and 3-year-old children by a preschool vision-screening program in Japan. Acta Med. Okayama 2007, 61, 9–16. [Google Scholar] [PubMed]
  4. Matsuo, T.; Matsuo, C. The prevalence of strabismus and amblyopia in Japanese elementary school children. Ophthalmic Epidemiol. 2005, 12, 31–36. [Google Scholar] [CrossRef] [PubMed]
  5. Matsuo, T.; Matsuo, C. Comparison of prevalence rates of strabismus and amblyopia in Japanese elementary school children between the years 2003 and 2005. Acta Med. Okayama 2007, 61, 329–334. [Google Scholar] [PubMed]
  6. Matsuo, T.; Matsuo, C.; Kio, K.; Ichiba, N.; Matsuoka, H. Is refraction with a hand-held autorefractometer useful in addition to visual acuity testing and questionnaires in preschool vision screening at 3.5 years in Japan? Acta Med. Okayama 2009, 63, 195–202. [Google Scholar] [PubMed]
  7. Matsuo, T.; Matsuo, C.; Kayano, M.; Mitsufuji, A.; Satou, C.; Matsuoka, H. Photorefraction with spot vision screener versus visual acuity testing as community-based preschool vision screening at the age of 3.5 years in Japan. Int. J. Environ. Res. Public Health 2022, 19, 8655. [Google Scholar] [CrossRef] [PubMed]
  8. Michaelides, M.; Moore, A.T. The genetics of strabismus. J. Med. Genet. 2004, 41, 641–646. [Google Scholar] [CrossRef] [PubMed]
  9. Maconachie, G.D.E.; Gottlob, I.; McLean, R.J. Risk factors and genetics in common comitant strabismus: A systematic review of the literature. JAMA Ophthalmol. 2013, 131, 1179–1186. [Google Scholar] [CrossRef] [PubMed]
  10. Ye, X.C.; Pegado, V.; Patel, M.S.; Wasserman, W.W. Strabismus genetics across a spectrum of eye misalignment disorders. Clin. Genet. 2014, 86, 103–111. [Google Scholar] [CrossRef] [PubMed]
  11. Matsuo, T.; Yamane, T.; Ohtsuki, H. Heredity versus abnormalities in pregnancy and delivery as risk factors for different types of comitant strabismus. J. Pediatr. Ophthalmol. Strabismus 2001, 38, 78–82. [Google Scholar] [CrossRef] [PubMed]
  12. Taira, Y.; Matsuo, T.; Yamane, T.; Hasebe, S.; Ohtsuki, H. Clinical features of comitant strabismus related to family history of strabismus or abnormalities in pregnancy and delivery. Jpn. J. Ophthalmol. 2003, 47, 208–213. [Google Scholar] [CrossRef] [PubMed]
  13. Podgor, M.J.; Remaley, N.A.; Chew, E. Associations between siblings for esotropia and exotropia. Arch. Ophthalmol. 1996, 114, 739–744. [Google Scholar] [CrossRef] [PubMed]
  14. Matsuo, T.; Hayashi, M.; Fujiwara, H.; Yamane, T.; Ohtsuki, H. Concordance of strabismic phenotypes in monozygotic versus multizygotic twins and other multiple births. Jpn. J. Ophthalmol. 2002, 46, 59–64. [Google Scholar] [CrossRef] [PubMed]
  15. Cotter, S.A.; Varma, R.; Tarczy-Hornoch, K.; McKean-Cowdin, R.; Lin, J.; Wen, G.; Wei, J.; Borchert, M.; Azen, S.P.; Torres, M.; et al. Risk factors associated with childhood strabismus: The Multi-Ethnic Pediatric Eye Disease and Baltimore Pediatric Eye Disease Studies. Ophthalmology 2011, 118, 2251–2261. [Google Scholar] [CrossRef] [PubMed]
  16. Jacobson, L.K.; Dutton, G.N. Periventricular leukomalacia: An important cause of visual and ocular motility dysfunction in children. Surv. Ophthalmol. 2000, 45, 1–13. [Google Scholar] [CrossRef] [PubMed]
  17. Engle, E.C. Genetic basis of congenital strabismus. Arch. Ophthalmol. 2007, 125, 189–195. [Google Scholar] [CrossRef] [PubMed]
  18. Graeber, C.P.; Hunter, D.G.; Engle, E.C. The genetic basis of incomitant strabismus: Consolidation of the current knowledge of the genetic foundations of disease. Semin. Ophthalmol. 2013, 28, 427–437. [Google Scholar] [CrossRef] [PubMed]
  19. Matsuo, T.; Ohtsuki, H.; Sogabe, Y.; Konishi, H.; Takenawa, K.; Watanabe, Y. Vertical abnormal retinal correspondence in three patients with congenital absence of the superior oblique muscle. Am. J. Ophthalmol. 1988, 106, 341–345. [Google Scholar] [CrossRef] [PubMed]
  20. Uchiyama, E.; Matsuo, T.; Imai, S.; Itoshima, E. Paretic side/normal side ratios of cross-sectional areas of the superior oblique muscle vary largely in idiopathic superior oblique palsy. Am. J. Ophthalmol. 2010, 149, 508–512. [Google Scholar] [CrossRef] [PubMed]
  21. Jiang, Y.; Matsuo, T.; Fujiwara, H.; Hasebe, S.; Ohtsuki, H.; Yasuda, T. ARIX gene polymorphisms in patients with congenital superior oblique muscle palsy. Br. J. Ophthalmol. 2004, 88, 263–267. [Google Scholar] [CrossRef] [PubMed]
  22. Jiang, Y.; Matsuo, T.; Fujiwara, H.; Hasebe, S.; Ohtsuki, H.; Yasuda, T. ARIX and PHOX2B polymorphisms in patients with congenital superior oblique muscle palsy. Acta Med. Okayama 2005, 59, 55–62. [Google Scholar] [PubMed]
  23. Imai, S.; Matsuo, T.; Itoshima, E.; Ohtsuki, H. Clinical features, ARIX and PHOX2B nucleotide changes in three families with congenital superior oblique muscle palsy. Acta Med. Okayama 2008, 62, 45–53. [Google Scholar] [PubMed]
  24. Ohkubo, S.I.; Matsuo, T.; Hasebe, K.; Shira, Y.H.; Itoshima, E.; Ohtsuki, H. Phenotype-phenotype and genotype-phenotype correlations in patients with idiopathic superior oblique muscle palsy. J. Hum. Genet. 2012, 57, 122–129. [Google Scholar] [CrossRef] [PubMed]
  25. Parikh, V.; Shugart, Y.Y.; Doheny, K.F.; Zhang, J.; Li, L.; Williams, J.; Hayden, D.; Craig, B.; Capo, H.; Chamblee, D.; et al. A strabismus susceptibility locus on chromosome 7p. Proc. Natl. Acad. Sci. USA 2003, 100, 12283–12288. [Google Scholar] [CrossRef] [PubMed]
  26. Rice, A.; Nsengimana, J.; Simmons, I.G.; Toomes, C.; Hoole, J.; Willoughby, C.E.; Cassidy, F.; Williams, G.A.; George, N.D.; Sheridan, E.; et al. Replication of the recessive STBMS1 locus but with dominant inheritance. Investig. Ophthalmol. Vis. Sci. 2009, 50, 3210–3217. [Google Scholar] [CrossRef] [PubMed]
  27. Fujiwara, H.; Matsuo, T.; Sato, M.; Yamane, T.; Kitada, M.; Hasebe, S.; Ohtsuki, H. Genome-wide search for strabismus susceptibility loci. Acta Med. Okayama 2003, 57, 109–116. [Google Scholar] [PubMed]
  28. Shaaban, S.; Matsuo, T.; Fujiwara, H.; Itoshima, E.; Furuse, T.; Hasebe, S.; Zhang, Q.; Ott, J.; Ohtsuki, H. Chromosomes 4q28.3 and 7q31.2 as new susceptibility loci for comitant strabismus. Investig. Ophthalmol. Vis. Sci. 2009, 50, 654–661. [Google Scholar] [CrossRef] [PubMed]
  29. Shaaban, S.; Matsuo, T.; Strauch, K.; Ohtsuki, H. Investigation of parent-of-origin effect in comitant strabismus using MOD score analysis. Mol. Vis. 2009, 15, 1351–1358. [Google Scholar] [PubMed]
  30. Khan, A.O.; Shinwari, J.; Abu Dhaim, N.; Khalil, D.; Al Sharif, L.; Al Tassan, N. Potential linkage of different phenotypic forms of childhood strabismus to a recessive susceptibility locus (16p13.12-p12.3). Mol. Vis. 2011, 17, 971–976. [Google Scholar] [PubMed]
  31. Bosten, J.M.; Hogg, R.E.; Bargary, G.; Goodbourn, P.T.; Lawrance-Owen, A.J.; Mollon, J.D. Suggestive association with ocular phoria at chromosome 6p22. Investig. Ophthalmol. Vis. Sci. 2014, 55, 345–352. [Google Scholar]
  32. Zhang, J.; Matsuo, T. MGST2 and WNT2 are candidate genes for comitant strabismus susceptibility in Japanese patients. PeerJ 2017, 5, e3935. [Google Scholar] [PubMed]
  33. Shaaban, S.; MacKinnon, S.; Andrews, C.; Staffieri, S.E.; Maconachie, G.D.E.; Chan, W.M.; Whitman, M.C.; Morton, S.U.; Yazar, S.; MacGregor, S.; et al. Genome-wide association study identifies a susceptibility locus for comitant esotropia and suggests a parent-of-origin effect. Investig. Ophthalmol. Vis. Sci. 2018, 59, 4054–4064. [Google Scholar]
  34. Plotnikov, D.; Shah, R.L.; Rodrigues, J.N.; Cumberland, P.M.; Rahi, J.S.; Hysi, P.G.; Atan, D.; Williams, C.; Guggenheim, J.A.; UK Biobank Eye and Vision Consortium. A commonly occurring genetic variant within the NPLOC4-TSPAN10-PDE6G gene cluster is associated with the risk of strabismus. Hum. Genet. 2019, 138, 723–737. [Google Scholar] [PubMed]
  35. Chaomulige; Matsuo, T.; Sugimoto, K.; Miyaji, M.; Hosoya, O.; Ueda, M.; Kobayashi, R.; Horii, T.; Hatada, I. Morphometric analysis of the eye by magnetic resonance imaging in MGST2-gene-deficient mice. Biomedicines 2024, 12, 370. [Google Scholar] [CrossRef] [PubMed]
  36. Ye, X.C.; Roslin, N.M.; Paterson, A.D.; Lyons, C.J.; Pegado, V.; Richmond, P.; Shyr, C.; Fornes, O.; Han, X.H.; Higginson, M.; et al. Linkage analysis identifies an isolated strabismus locus at 14q12 overlapping with FOXG1 syndrome region. J. Med. Genet. 2022, 59, 46–55. [Google Scholar] [PubMed]
  37. Uffelmann, E.; Huang, Q.Q.; Munung, N.S.; de Vries, J.; Okada, Y.; Martin, A.R.; Martin, H.C.; Lappalainen, T.; Posthuma, D. Genome-wide association studies. Nat. Rev. 2021, 1, 59. [Google Scholar]
  38. Kanai, M.; Akiyama, M.; Takahashi, A.; Matoba, N.; Momozawa, Y.; Ideda, M.; Iwata, N.; Ikegawa, S.; Hirata, M.; Matsuda, K.; et al. Genetic analysis of quantitative traits in the Japanese population links cell types to complex human diseases. Nat. Genet. 2018, 50, 390–400. [Google Scholar] [CrossRef] [PubMed]
  39. Morita, Y.; Kamatani, Y.; Ito, H.; Ikegawa, S.; Kawaguchi, T.; Kawaguchi, S.; Takahashi, M.; Terao, C.; Ito, S.; Nishitani, K.; et al. Improved genetic prediction of the risk of knee osteoarthritis using the risk factor-based polygenic score. Arthritis Res. Ther. 2023, 25, 103. [Google Scholar] [PubMed]
  40. Shoaib, M.; Ye, Q.; IglayReger, H.; Tan, M.H.; Boehnke, M.; Burant, C.F.; Soleimanpour, S.A.; Talium, A.S.G. Evaluation of polygenic risk scores to differentiate between type 1 and type 2 diabetes. Genet. Epidemiol. 2023, 47, 303–313. [Google Scholar] [CrossRef] [PubMed]
  41. Ye, X.C.; van der Lee, R.; Wasserman, W.W. Curation and bioinformatic analysis of strabismus genes supports functional heterogeneity and proposes candidate genes with connections to RASopathies. Gene 2019, 697, 213–226. [Google Scholar] [CrossRef] [PubMed]
  42. Ma, W.X.; Huang, X.G.; Yang, T.K.; Yao, J.Y. Involvement of dysregulated coding and long non-coding RNAs in the pathogenesis of strabismus. Mol. Med. Rep. 2018, 17, 7737–7745. [Google Scholar] [CrossRef] [PubMed]
  43. Corral-Juan, M.; Serrano-Munuera, C.; Rabano, A.; Cota-Gonzalez, D.; Segarra-Roca, A.; Ispierto, L.; Cano-Orgaz, A.T.; Adarmes, A.D.; Mendez-del-Barrio, C.; Jesus, S.; et al. Clinical, genetic and neuropathological characterization of spinocerebellar ataxia type 37. Brain 2018, 141, 1981–1997. [Google Scholar] [CrossRef] [PubMed]
  44. Matsuo, T. The genes involved in the morphogenesis of the eye. Jpn. J. Ophthalmol. 1993, 37, 215–251. [Google Scholar] [PubMed]
  45. Cvekl, A.; Wang, W.L. Retinoic acid signaling in mammalian eye development. Exp. Eye Res. 2009, 89, 280–291. [Google Scholar] [CrossRef] [PubMed]
  46. Thompson, B.; Katsanis, N.; Apostolopoulos, N.; Thompson, D.C.; Nebert, D.W.; Vasiliou, V. Genetics and functions of the retinoic acid pathway, with special emphasis on the eye. Hum. Genom. 2019, 13, 61. [Google Scholar] [CrossRef] [PubMed]
  47. Hebert, S.L.; Fitzpatrick, K.P.; McConnell, S.A.; Cucak, A.; Yuan, C.; McLoon, L.K. Effects of retinoic acid signaling on extraocular muscle myogenic precursor cells in vitro. Exp. Cell Res. 2017, 361, 101–111. [Google Scholar] [CrossRef] [PubMed]
  48. Comai, G.E.; Tesařová, M.; Dupé, V.; Rhinn, M.; Vallecillo-García, P.; da Silva, F.; Feret, B.; Exelby, K.; Dollé, P.; Carlsson, L.; et al. Local retinoic acid signaling directs emergence of the extraocular muscle functional unit. PLoS Biol. 2020, 18, e3000902. [Google Scholar] [CrossRef] [PubMed]
  49. Nagai, A.; Hirata, M.; Kamatani, Y.; Muto, K.; Matsuda, K.; Kiyohara, Y.; Ninomiya, T.; Tamakoshi, A.; Yamagata, Z.; Mushiroda, T.; et al. Overview of the BioBank Japan project: Study design and profile. J. Epidemiol. 2017, 27, S2–S8. [Google Scholar] [CrossRef] [PubMed]
  50. Setoh, K.; Matsuda, F. Chapter 7. Cohort profile: The Nagahama prospective genome cohort for comprehensive human bioscience (The Nagahama Study). In Socio-Life Science and the COVID-19 Outbreak: Public Health and Public Policy; Series of Economics, Law, and Institutions on Asia Pacific; Yano, M., Matsuda, F., Sakuntabhai, A., Hirota, S., Eds.; eBook; Springer Nature: Singapore, 2022; pp. 127–143. [Google Scholar]
  51. The 1000 Genomes Project Consortium. A map of human genome variation from population-scale sequencing. Nature 2010, 467, 1061–1073. [Google Scholar] [CrossRef] [PubMed]
  52. Sudmant, P.H.; Rausch, T.; Gardner, E.J.; Handsaker, R.E.; Abyzov, A.; Huddleston, J.; Zhang, Y.; Ye, K.; Jun, G.; Hsi-Yang Fritz, M.; et al. An integrated map of structural variation in 2504 human genomes. Nature 2015, 526, 75–81. [Google Scholar] [CrossRef] [PubMed]
  53. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [PubMed]
  54. Li, Y.; Willer, C.; Sanna, S.; Abecasis, G. Genotype imputation. Annu. Rev. Genom. Hum. Genet. 2009, 10, 387–406. [Google Scholar] [CrossRef] [PubMed]
  55. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef] [PubMed]
  56. Browning, B.L.; Zhou, Y.; Browning, S.R. A one-penny imputed genome from next-generation reference panels. Am. J. Hum. Genet. 2018, 103, 338–348. [Google Scholar] [CrossRef] [PubMed]
  57. Browning, S.R.; Browning, B.L. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am. J. Hum. Genet. 2007, 81, 1084–1097. [Google Scholar] [CrossRef] [PubMed]
  58. Loh, P.R.; Palamara, P.F.; Price, A.L. Fast and accurate long-range phasing in a UK Biobank cohort. Nat. Genet. 2016, 48, 811–816. [Google Scholar] [CrossRef] [PubMed]
  59. Loh, P.R.; Danecek, P.; Palamara, P.F.; Fuchsberger, C.; Reshef, Y.A.; Finucane, H.K.; Schoenherr, S.; Forer, L.; McCarthy, S.; Abecasis, G.R.; et al. Reference-based phasing using the Haplotype Reference Consortium panel. Nat. Genet. 2016, 48, 1443–1448. [Google Scholar] [CrossRef] [PubMed]
  60. Das, S.; Forer, L.; Schonherr, S.; Sidore, C.; Locke, A.E.; Kwong, A.; Vrieze, S.I.; Chew, E.Y.; Levy, S.; McGue, M.; et al. Next-generation genotype imputation service and methods. Nat. Genet. 2016, 48, 1284–1287. [Google Scholar] [CrossRef] [PubMed]
  61. Mbatchou, J.; Barnard, L.; Backman, J.; Marcketta, A.; Kosmicki, J.A.; Ziyatdinov, A.; Benner, C.; O’Dushlaine, C.; Barber, M.; Boutkov, B.; et al. Computationally efficient whole-genome regression for quantitative and binary traits. Nat. Genet. 2021, 53, 1097–1103. [Google Scholar] [CrossRef] [PubMed]
  62. Abraham, G.; Qiu, Y.; Inouye, M. FlashPCA2: Principal component analysis of Biobank-scale genotype datasets. Bioinformatics 2017, 33, 2776–2778. [Google Scholar] [CrossRef] [PubMed]
  63. Pruim, R.J.; Welch, R.P.; Sanna, S.; Teslovich, T.M.; Chines, P.S.; Gliedt, T.P.; Boehnke, M.; Abecasis, G.R.; Willer, C.J. LocusZoom: Regional visualization of genome-wide association scan results. Bioinformatics 2010, 26, 2336–2337. [Google Scholar] [CrossRef] [PubMed]
Ijms 25 06986 g001
Figure 1. Flow charts for selection process of case samples and control samples in BBJ (180K) and BBJ (ASA). BBJ, BioBank Japan; SO palsy, idiopathic superior oblique muscle palsy.
Figure 1. Flow charts for selection process of case samples and control samples in BBJ (180K) and BBJ (ASA). BBJ, BioBank Japan; SO palsy, idiopathic superior oblique muscle palsy.
Ijms 25 06986 g001
Ijms 25 06986 g002
Figure 2. SNP selection and SNP imputation in case samples and control samples in BBJ (180K). BBJ, BioBank Japan.
Figure 2. SNP selection and SNP imputation in case samples and control samples in BBJ (180K). BBJ, BioBank Japan.
Ijms 25 06986 g002
Ijms 25 06986 g003
Figure 3. SNP selection and SNP imputation in case samples and control samples in BBJ (ASA). BBJ, BioBank Japan.
Figure 3. SNP selection and SNP imputation in case samples and control samples in BBJ (ASA). BBJ, BioBank Japan.
Ijms 25 06986 g003
Ijms 25 06986 g004
Figure 4. SNP selection and SNP imputation in case samples and control samples in Nagahama Cohort.
Figure 4. SNP selection and SNP imputation in case samples and control samples in Nagahama Cohort.
Ijms 25 06986 g004
Ijms 25 06986 g005
Figure 5. Comparison of SNP dosage after SNP imputation. Note that there is no difference between 1KG JPT in the vertical axis and imputed SNPs in the horizontal axis. SNPs in chromosome 6 (red box) were used as shown in Manhattan plot where red line indicates p = 5 × 10−8 and blue line indicates p = 1 × 10−6. Upper panels: GWAS with BBJ (180K). Middle panels: GWAS with BBJ (ASA). Lower panels: GWAS with Nagahama Cohort. 1KG JPT, Japanese population of 1000 Genomes Project.
Figure 5. Comparison of SNP dosage after SNP imputation. Note that there is no difference between 1KG JPT in the vertical axis and imputed SNPs in the horizontal axis. SNPs in chromosome 6 (red box) were used as shown in Manhattan plot where red line indicates p = 5 × 10−8 and blue line indicates p = 1 × 10−6. Upper panels: GWAS with BBJ (180K). Middle panels: GWAS with BBJ (ASA). Lower panels: GWAS with Nagahama Cohort. 1KG JPT, Japanese population of 1000 Genomes Project.
Ijms 25 06986 g005
Ijms 25 06986 g006
Figure 6. Comparison of SNP dosage after SNP imputation. Note that there is no difference between BBJ (180K) and BBJ (ASA) in the case (top left panel) and the control (top right panel). Furthermore, note that there is no difference between BBJ (180K) and Nagahama Cohort (bottom left panel) and between BBJ (ASA) and Nagahama Cohort (bottom right panel) in the case.
Figure 6. Comparison of SNP dosage after SNP imputation. Note that there is no difference between BBJ (180K) and BBJ (ASA) in the case (top left panel) and the control (top right panel). Furthermore, note that there is no difference between BBJ (180K) and Nagahama Cohort (bottom left panel) and between BBJ (ASA) and Nagahama Cohort (bottom right panel) in the case.
Ijms 25 06986 g006
Ijms 25 06986 g007aIjms 25 06986 g007bIjms 25 06986 g007cIjms 25 06986 g007dIjms 25 06986 g007e
Figure 7. Meaningful peaks (red circles) formed by proximal SNPs in each Manhattan plot with control cohort of BBJ (180K), BBJ (ASA), and Nagahama. Chromosomal location on the horizontal axis versus −log10(p) on the vertical axis of each Manhattan plot. Red line indicates p = 5 × 10−8 and blue line indicates p = 1 × 10−6. Q–Q (log quantile–quantile p-value) plots with expected −log10(p) on the horizontal axis versus observed −log10(p) on the vertical axis, shown next to each Manhattan plot. Red line in each Q–Q plot indicate a correlation line. Only the group of exotropia is suggested to have a structured bias from the Q–Q plots. The case group is esotropia (A), exotropia (B), idiopathic superior oblique muscle palsy (C), comitant strabismus of esotropia and exotropia (D), and all cases of esotropia, exotropia, and idiopathic superior oblique muscle palsy (E).
Figure 7. Meaningful peaks (red circles) formed by proximal SNPs in each Manhattan plot with control cohort of BBJ (180K), BBJ (ASA), and Nagahama. Chromosomal location on the horizontal axis versus −log10(p) on the vertical axis of each Manhattan plot. Red line indicates p = 5 × 10−8 and blue line indicates p = 1 × 10−6. Q–Q (log quantile–quantile p-value) plots with expected −log10(p) on the horizontal axis versus observed −log10(p) on the vertical axis, shown next to each Manhattan plot. Red line in each Q–Q plot indicate a correlation line. Only the group of exotropia is suggested to have a structured bias from the Q–Q plots. The case group is esotropia (A), exotropia (B), idiopathic superior oblique muscle palsy (C), comitant strabismus of esotropia and exotropia (D), and all cases of esotropia, exotropia, and idiopathic superior oblique muscle palsy (E).
Ijms 25 06986 g007aIjms 25 06986 g007bIjms 25 06986 g007cIjms 25 06986 g007dIjms 25 06986 g007e
Ijms 25 06986 g008aIjms 25 06986 g008bIjms 25 06986 g008cIjms 25 06986 g008dIjms 25 06986 g008e
Figure 8. Regional plots by Locus Zoom of common peaks with a loose standard of p < 1 × 10−6 in 3 series of GWASs with 3 different control cohorts of BBJ (180K), BBJ (ASA), and Nagahama Cohort. The case group is esotropia (A–C) and idiopathic superior oblique muscle palsy (D,E).
Figure 8. Regional plots by Locus Zoom of common peaks with a loose standard of p < 1 × 10−6 in 3 series of GWASs with 3 different control cohorts of BBJ (180K), BBJ (ASA), and Nagahama Cohort. The case group is esotropia (A–C) and idiopathic superior oblique muscle palsy (D,E).
Ijms 25 06986 g008aIjms 25 06986 g008bIjms 25 06986 g008cIjms 25 06986 g008dIjms 25 06986 g008e
Table 1. Common peaks with a loose significant level at p < 1 × 10−6 in GWAS with three different control cohorts in Japanese population.
Table 1. Common peaks with a loose significant level at p < 1 × 10−6 in GWAS with three different control cohorts in Japanese population.
GWASSNPChromosomal LocationGenesFull NameFunction
Esotropia vs. BBJ (180K)
Esotropia vs. BBJ (ASA)
Esotropia vs. Nagahama		35,665,716
35,612,917
35,384,605		1p34.3SFPQSplicing factor proline and glutamine richAlternative mRNA splicing
Esotropia vs. BBJ (180K)
Esotropia vs. BBJ (ASA)
Esotropia vs. Nagahama		22,1040,171
22,1040,171
22,1030,589		1q41HLX-AS1HLX antisense RNA 1
Esotropia vs. BBJ (180K)
Esotropia vs. BBJ (ASA)
Esotropia vs. Nagahama		101,419,925
101,406,820
101,488,901		8q22.2MIR4471microRNA 4471Post-transcriptional regulation of gene expression
SO palsy vs. BBJ (180K)
SO palsy vs. BBJ (ASA)
SO palsy vs. Nagahama		58,240,503
58,357,705
58,243,570		1p32.2-p32.1DAB1DAB adaptor protein 1Neuronal migration
SO palsy vs. BBJ (180K)
SO palsy vs. BBJ (ASA)
SO palsy vs. Nagahama		51,793,524
51,842658
51,778,296		5q11.2PELOPelota mRNA surveillance and ribosome rescue factorCell cycle control
Cited from the database “gene” in National Center for Biotechnology Information (NCBI) in U.S.A. GWAS in condition of the exclusion of identity-by-descent (IBD) > 0.1 and no regenie prediction. BBJ, BioBank Japan; Nagahama, Nagahama Cohort; SO palsy, idiopathic superior oblique (SO) muscle palsy. See Figure 8 for regional plots corresponding to each locus and Supplementary Table S1 for original data.
Table 2. Peaks with a standard significant level at p < 5 × 10−8 in GWAS with the control cohort, BBJ (180K).
Table 2. Peaks with a standard significant level at p < 5 × 10−8 in GWAS with the control cohort, BBJ (180K).
GWASSNPChromosomal LocationGenesFull NameFunction
Esotropia vs. BBJ (180K)174,265,7792q31.1CDCA7Cell division cycle associated 7Cell transformation
Exotropia vs. BBJ (180K)29,677,8776p22.1HLA-FMajor histocompatibility complex, class I, F
SO palsy vs. BBJ (180K)58,240,5031p32.2DAB1-AS1DAB1 antisense RNA 1Neuronal migration
ET+XT vs. BBJ (180K)195,019,9263q29ACAP2ArfGAP with coiled-coil, ankyrin repeat, and PH domains 2Actin filament-based process
ET+XT vs. BBJ (180K)29,677,8776p22.1HLA-FMajor histocompatibility complex, class I, F
ET+XT+SO vs. BBJ (180K)25,228,1133p24.2RARBRetinoic acid receptor betaEmbryonic morphogenesis, cell growth and differentiation
ET+XT+SO vs. BBJ (180K)29,677,8776p22.1HLA-FMajor histocompatibility complex, class I, F
Cited from the database “gene” in National Center for Biotechnology Information (NCBI) in U.S.A. GWAS in condition of the exclusion of identity-by-descent (IBD) > 0.1 and no regenie prediction. BBJ, BioBank Japan; SO palsy, idiopathic superior oblique (SO) muscle palsy; ET, esotropia; XT, exotropia. See Supplementary Figure S1 for regional plots corresponding to each locus.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Matsuo, T.; Hamasaki, I.; Kamatani, Y.; Kawaguchi, T.; Yamaguchi, I.; Matsuda, F.; Saito, A.; Nakazono, K.; Kamitsuji, S. Genome-Wide Association Study with Three Control Cohorts of Japanese Patients with Esotropia and Exotropia of Comitant Strabismus and Idiopathic Superior Oblique Muscle Palsy. Int. J. Mol. Sci. 2024, 25, 6986. https://doi.org/10.3390/ijms25136986

AMA Style

Matsuo T, Hamasaki I, Kamatani Y, Kawaguchi T, Yamaguchi I, Matsuda F, Saito A, Nakazono K, Kamitsuji S. Genome-Wide Association Study with Three Control Cohorts of Japanese Patients with Esotropia and Exotropia of Comitant Strabismus and Idiopathic Superior Oblique Muscle Palsy. International Journal of Molecular Sciences. 2024; 25(13):6986. https://doi.org/10.3390/ijms25136986

Chicago/Turabian Style

Matsuo, Toshihiko, Ichiro Hamasaki, Yoichiro Kamatani, Takahisa Kawaguchi, Izumi Yamaguchi, Fumihiko Matsuda, Akira Saito, Kazuyuki Nakazono, and Shigeo Kamitsuji. 2024. "Genome-Wide Association Study with Three Control Cohorts of Japanese Patients with Esotropia and Exotropia of Comitant Strabismus and Idiopathic Superior Oblique Muscle Palsy" International Journal of Molecular Sciences 25, no. 13: 6986. https://doi.org/10.3390/ijms25136986

APA Style

Matsuo, T., Hamasaki, I., Kamatani, Y., Kawaguchi, T., Yamaguchi, I., Matsuda, F., Saito, A., Nakazono, K., & Kamitsuji, S. (2024). Genome-Wide Association Study with Three Control Cohorts of Japanese Patients with Esotropia and Exotropia of Comitant Strabismus and Idiopathic Superior Oblique Muscle Palsy. International Journal of Molecular Sciences, 25(13), 6986. https://doi.org/10.3390/ijms25136986

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop