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

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.

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.

Idiopathic superior oblique muscle palsy is a major type of paralytic, non-comitant strabismus and presents vertical and cyclo-torsional deviation of one eye against the other eye, with a large vertical fusion range and abnormal head posture such as head tilt. Genetic background is considered to play a role in its development, as patients with idiopathic superior oblique muscle palsy have varying degrees of muscle hypoplasia and, rarely, the complete absence of the muscle, that is, aplasia. In this study, whole genome sequencing was performed, and single nucleotide variations and short insertions/deletions (SNVs/InDels) were annotated in two patients each in three small families (six patients in total) with idiopathic superior oblique muscle palsy, in addition to three normal individuals in one family. At first, linkage analysis was carried out in the three families and SNVs/InDels in chromosomal loci with negative LOD scores were excluded. Next, SNVs/InDels shared by the six patients, but not by the three normal individuals, were chosen. SNVs/InDels were further narrowed down by choosing low-frequency (<1%) or non-registered SNVs/InDels in four databases for the Japanese population, and then by choosing SNVs/InDels with functional influence, leading to one candidate gene, SSTR5-AS1 in chromosome 16. The six patients were heterozygous for 13-nucleotide deletion in SSTR5-AS1, except for one homozygous patient, while the three normal individuals were wild type. Targeted polymerase chain reaction (PCR) and direct sequencing of PCR products confirmed the 13-nucleotide deletion in SSTR5-AS1. In the face of newly-registered SSTR5-AS1 13-nucleotide deletion at a higher frequency in a latest released database for the Japanese population, the skipping of low-frequency and non-registration sorting still resulted in only 13 candidate genes including SSTR5-AS1 as common variants. The skipping of linkage analysis also led to the same set of 13 candidate genes. Different testing strategies that consisted of linkage analysis and simple unintentional bioinformatics could reach candidate genes in three small families with idiopathic superior oblique muscle palsy.

Japanese GWAS identifies variants for bust-size, dysmenorrhea, and menstrual fever that are eQTLs for relevant protein-coding or long non-coding RNAs.

Traits related to primary and secondary sexual characteristics greatly impact females during puberty and day-to-day adult life. Therefore, we performed a GWAS analysis of 11,348 Japanese female volunteers and 22 gynecology-related phenotypic variables, and identified significant associations for bust-size, menstrual pain (dysmenorrhea) severity, and menstrual fever. Bust-size analysis identified significant association signals in CCDC170-ESR1 (rs6557160; P = 1.7 × 10-16) and KCNU1-ZNF703 (rs146992477; P = 6.2 × 10-9) and found that one-third of known European-ancestry associations were also present in Japanese. eQTL data points to CCDC170 and ZNF703 as those signals' functional targets. For menstrual fever, we identified a novel association in OPRM1 (rs17181171; P = 2.0 × 10-8), for which top variants were eQTLs in multiple tissues. A known dysmenorrhea signal near NGF replicated in our data (rs12030576; P = 1.1 × 10-19) and was associated with RP4-663N10.1 expression, a putative lncRNA enhancer of NGF, while a novel dysmenorrhea signal in the IL1 locus (rs80111889; P = 1.9 × 10-16) contained SNPs previously associated with endometriosis, and GWAS SNPs were most significantly associated with IL1A expression. By combining regional imputation with colocalization analysis of GWAS/eQTL signals along with integrated annotation with epigenomic data, this study further refines the sets of candidate causal variants and target genes for these known and novel gynecology-related trait loci.

Genome-wide association study in Japanese females identifies fifteen novel skin-related trait associations.

Skin trait variation impacts quality-of-life, especially for females from the viewpoint of beauty. To investigate genetic variation related to these traits, we conducted a GWAS of various skin phenotypes in 11,311 Japanese women and identified associations for age-spots, freckles, double eyelids, straight/curly hair, eyebrow thickness, hairiness, and sweating. In silico annotation with RoadMap Epigenomics epigenetic state maps and colocalization analysis of GWAS and GTEx Project eQTL signals provided information about tissue specificity, candidate causal variants, and functional target genes. Novel signals for skin-spot traits neighboured AKAP1/MSI2 (rs17833789; P = 2.2 × 10-9), BNC2 (rs10810635; P = 2.1 × 10-22), HSPA12A (rs12259842; P = 7.1 × 10-11), PPARGC1B (rs251468; P = 1.3 × 10-21), and RAB11FIP2 (rs10444039; P = 5.6 × 10-21). HSPA12A SNPs were the only protein-coding gene eQTLs identified across skin-spot loci. Double edged eyelid analysis identified that a signal around EMX2 (rs12570134; P = 8.2 × 10-15) was also associated with expression of EMX2 and the antisense-RNA gene EMX2OS in brain putamen basal ganglia tissue. A known hair morphology signal in EDAR was associated with both eyebrow thickness (rs3827760; P = 1.7 × 10-9) and straight/curly hair (rs260643; P = 1.6 × 10-103). Excessive hairiness signals' top SNPs were also eQTLs for TBX15 (rs984225; P = 1.6 × 10-8), BCL2 (rs7226979; P = 7.3 × 10-11), and GCC2 and LIMS1 (rs6542772; P = 2.2 × 10-9). For excessive sweating, top variants in two signals in chr2:28.82-29.05 Mb (rs56089836; P = 1.7 × 10-11) were eQTLs for either PPP1CB or PLB1, while a top chr16:48.26-48.45 Mb locus SNP was a known ABCC11 missense variant (rs6500380; P = 6.8 × 10-10). In total, we identified twelve loci containing sixteen association signals, of which fifteen were novel. These findings will help dermatologic researchers better understand the genetic underpinnings of skin-related phenotypic variation in human populations.

Genome-wide association study of self-reported food reactions in Japanese identifies shrimp and peach specific loci in the HLA-DR/DQ gene region.

Food allergy is an increasingly important health problem in the world. Several genome-wide association studies (GWAS) focused on European ancestry samples have identified food allergy-specific loci in the HLA class II region. We conducted GWAS of self-reported reactivity with common foods using the data from 11011 Japanese women and identified shrimp and peach allergy-specific loci in the HLA-DR/DQ gene region tagged by rs74995702 (P = 6.30 × 10-17, OR = 1.91) and rs28359884 (P = 2.3 × 10-12, OR = 1.80), respectively. After HLA imputation using a Japanese population-specific reference, the most strongly associated haplotype was HLA-DRB1*04:05-HLA-DQB1*04:01 for shrimp allergy (P = 3.92 × 10-19, OR = 1.99) and HLA-DRB1*09:01-HLA-DQB1*03:03 for peach allergy (P = 1.15 × 10-7, OR = 1.68). Additionally, both allergies' associated variants were eQTLs for several HLA genes, with HLA-DQA2 the single eQTL gene shared between the two traits. Our study suggests that allergy to certain foods may be related to genetic differences that tag both HLA alleles having particular epitope binding specificities as well as variants modulating expression of particular HLA genes. Investigating this further could increase our understanding of food allergy aetiology and potentially lead to better therapeutic strategies for allergen immunotherapies.

Japan PGx Data Science Consortium Database: SNPs and HLA genotype data from 2994 Japanese healthy individuals for pharmacogenomics studies.

Japan Pharmacogenomics Data Science Consortium (JPDSC) has assembled a database for conducting pharmacogenomics (PGx) studies in Japanese subjects. The database contains the genotypes of 2.5 million single-nucleotide polymorphisms (SNPs) and 5 human leukocyte antigen loci from 2994 Japanese healthy volunteers, as well as 121 kinds of clinical information, including self-reports, physiological data, hematological data and biochemical data. In this article, the reliability of our data was evaluated by principal component analysis (PCA) and association analysis for hematological and biochemical traits by using genome-wide SNP data. PCA of the SNPs showed that all the samples were collected from the Japanese population and that the samples were separated into two major clusters by birthplace, Okinawa and other than Okinawa, as had been previously reported. Among 87 SNPs that have been reported to be associated with 18 hematological and biochemical traits in genome-wide association studies (GWAS), the associations of 56 SNPs were replicated using our data base. Statistical power simulations showed that the sample size of the JPDSC control database is large enough to detect genetic markers having a relatively strong association even when the case sample size is small. The JPDSC database will be useful as control data for conducting PGx studies to explore genetic markers to improve the safety and efficacy of drugs either during clinical development or in post-marketing.

Japanese population structure, based on SNP genotypes from 7003 individuals compared to other ethnic groups: effects on population-based association studies.

Because population stratification can cause spurious associations in case-control studies, understanding the population structure is important. Here, we examined Japanese population structure by "Eigenanalysis," using the genotypes for 140,387 SNPs in 7003 Japanese individuals, along with 60 European, 60 African, and 90 East-Asian individuals, in the HapMap project. Most Japanese individuals fell into two main clusters, Hondo and Ryukyu; the Hondo cluster includes most of the individuals from the main islands in Japan, and the Ryukyu cluster includes most of the individuals from Okinawa. The SNPs with the greatest frequency differences between the Hondo and Ryukyu clusters were found in the HLA region in chromosome 6. The nonsynonymous SNPs with the greatest frequency differences between the Hondo and Ryukyu clusters were the Val/Ala polymorphism (rs3827760) in the EDAR gene, associated with hair thickness, and the Gly/Ala polymorphism (rs17822931) in the ABCC11 gene, associated with ear-wax type. Genetic differentiation was observed, even among different regions in Honshu Island, the largest island of Japan. Simulation studies showed that the inclusion of different proportions of individuals from different regions of Japan in case and control groups can lead to an inflated rate of false-positive results when the sample sizes are large.