Somatic mutations in the ARID1A tumor-suppressor gene have been frequently identified in ovarian clear cell carcinoma (CCC) cases. BAF250a encoded by ARID1A is a member of the SWI/SNF complex, but the expression and mutation status of other SWI/SNF subunits have not been explored. The current study aimed to elucidate the biological and clinical significance of the SWI/SNF complex subunits, by assessing the expression and mutation status of SWI/SNF subunits, and distinct genomic aberrations associated with their expression. Of 82 CCC specimens, 38 samples presented no BAF250a expression, and 50 samples exhibited the loss of at least one subunit of the SWI/SNF complex. Cases which lack at least one SWI/SNF complex component exhibited significantly more advanced stages, faster growth and stronger nuclear atypia compared with SWI/SNF-positive samples (p<0.05). Although BAF250a expression is not related to poor prognosis, the group presenting the loss of at least one SWI/SNF complex subunit exhibited significantly shorter overall and progression-free survivals (p<0.05). A multivariate analysis suggested that the expression status of the SWI/SNF complex serves as an independent prognostic factor (p<0.005). The cases positive for all SWI/SNF subunits demonstrated significantly greater DNA copy number alterations, such as amplification at chromosomes 8q.24.3 and 20q.13.2-20q.13.33 (including ZNF217) and deletion at chromosomes 13q12.11-13q14.3 (including RB1), 17p13.2-17p13.1 (including TP53) and 19p13.2-19p13.12. In conclusion, the CCCs exhibiting the loss of one or multiple SWI/SNF complex subunits demonstrated aggressive behaviors and poor prognosis, whereas the CCCs with positive expression for all SWI/SNF components presented more copy number alterations and a favorable prognosis.

Nonalcoholic fatty liver disease (NAFLD) includes a broad range of liver pathologies from simple steatosis to cirrhosis and fibrosis, in which a subtype accompanying hepatocyte degeneration and fibrosis is classified as nonalcoholic steatohepatitis (NASH). NASH accounts for approximately 10-30% of NAFLD and causes a higher frequency of liver-related death, and its progression of NASH has been considered to be complex involving multiple genetic factors interacting with the environment and lifestyle.

HLA-DRB1, especially the shared epitope (SE), is strongly associated with rheumatoid arthritis (RA). However, recent studies have shown that SE is at most weakly associated with RA without anti-citrullinated peptide/protein antibody (ACPA). We have recently reported that ACPA-negative RA is associated with specific HLA-DRB1 alleles and diplotypes. Here, we attempted to detect genetically different subsets of ACPA-negative RA by classifying ACPA-negative RA patients into two groups based on their positivity for rheumatoid factor (RF). HLA-DRB1 genotyping data for totally 954 ACPA-negative RA patients and 2,008 healthy individuals in two independent sets were used. HLA-DRB1 allele and diplotype frequencies were compared among the ACPA-negative RF-positive RA patients, ACPA-negative RF-negative RA patients, and controls in each set. Combined results were also analyzed. A similar analysis was performed in 685 ACPA-positive RA patients classified according to their RF positivity. As a result, HLA-DRB1*04:05 and *09:01 showed strong associations with ACPA-negative RF-positive RA in the combined analysis (p = 8.8×10(-6) and 0.0011, OR: 1.57 (1.28-1.91) and 1.37 (1.13-1.65), respectively). We also found that HLA-DR14 and the HLA-DR8 homozygote were associated with ACPA-negative RF-negative RA (p = 0.00022 and 0.00013, OR: 1.52 (1.21-1.89) and 3.08 (1.68-5.64), respectively). These association tendencies were found in each set. On the contrary, we could not detect any significant differences between ACPA-positive RA subsets. As a conclusion, ACPA-negative RA includes two genetically distinct subsets according to RF positivity in Japan, which display different associations with HLA-DRB1. ACPA-negative RF-positive RA is strongly associated with HLA-DRB1*04:05 and *09:01. ACPA-negative RF-negative RA is associated with DR14 and the HLA-DR8 homozygote.

Tubular injury and interstitial fibrosis are the hallmarks of chronic kidney disease. While recent studies have verified that proximal tubular injury triggers interstitial fibrosis, the impact of fibrosis on tubular injury and regeneration remains poorly understood. We generated a novel mouse model expressing diphtheria toxin receptor on renal fibroblasts to allow for the selective disruption of renal fibroblast function. Administration of diphtheria toxin induced upregulation of the tubular injury marker Ngal and caused tubular proliferation in healthy kidneys, whereas administration of diphtheria toxin attenuated tubular regeneration in fibrotic kidneys. Microarray analysis revealed down-regulation of the retinol biosynthesis pathway in diphtheria toxin-treated kidneys. Healthy proximal tubules expressed retinaldehyde dehydrogenase 2 (RALDH2), a rate-limiting enzyme in retinoic acid biosynthesis. After injury, proximal tubules lost RALDH2 expression, whereas renal fibroblasts acquired strong expression of RALDH2 during the transition to myofibroblasts in several models of kidney injury. The retinoic acid receptor (RAR) RARγ was expressed in proximal tubules both with and without injury, and αB-crystallin, the product of an RAR target gene, was strongly expressed in proximal tubules after injury. Furthermore, BMS493, an inverse agonist of RARs, significantly attenuated tubular proliferation in vitro. In human biopsy tissue from patients with IgA nephropathy, detection of RALDH2 in the interstitium correlated with older age and lower kidney function. These results suggest a role of retinoic acid signaling and cross-talk between fibroblasts and tubular epithelial cells during tubular injury and regeneration, and may suggest a beneficial effect of fibrosis in the early response to injury.

The recently emerged pachychoroid concept has changed the understanding of age-related macular degeneration (AMD), which is a major cause of blindness; recent studies attributed AMD in part to pachychoroid disease central serous chorioretinopathy (CSC), suggesting the importance of elucidating the CSC pathogenesis. Our large genome-wide association study followed by validation studies in three independent Japanese and European cohorts, consisting of 1546 CSC samples and 13,029 controls, identified two novel CSC susceptibility loci: TNFRSF10A-LOC389641 and near GATA5 (rs13278062, odds ratio = 1.35, P = 1.26 × 10-13; rs6061548, odds ratio = 1.63, P = 5.36 × 10-15). A T allele at TNFRSF10A-LOC389641 rs13278062, a risk allele for CSC, is known to be a risk allele for AMD. This study not only identified new susceptibility genes for CSC, but also improves the understanding of the pathogenesis of AMD.

Although intensive efforts have been undertaken to elucidate the genetic background of immunoglobulin A nephropathy (IgAN), genetic factors associated with the pathogenesis of this disease are still not well understood. We designed a case-control association study that was based on linkage disequilibrium among single-nucleotide polymorphisms (SNPs) in the selectin gene cluster on chromosome 1q24-25, and we found two SNPs in the E-selectin gene (SELE8 and SELE13) and six SNPs in the L-selectin gene (SELL1, SELL4, SELL5, SELL6, SELL10, and SELL11) that were significantly associated with IgAN in Japanese patients. All eight SNPs were in almost complete linkage disequilibrium. SELE8 and SELL10 caused amino acid substitutions from His to Tyr and from Pro to Ser (chi2=9.02, P=.0026, odds ratio = 2.73 [95% confidence interval [CI] 1.38--5.38] for His-to-Tyr substitutions; chi2=17.4, P=.000031, odds ratio = 3.61 [95% CI 1.91--6.83] for Pro-to-Ser substitutions), and SELL1 could affect promoter activity of the L-selectin gene (chi2=19.5, P=.000010, odds ratio = 3.77 [95% CI 2.02--7.05]). The TGT haplotype at these three loci was associated significantly with IgAN (chi2=18.67, P=.000016, odds ratio = 1.88 [95% CI 1.41--2.51]). Our results suggest that these eight SNPs in selectin genes may be useful for screening populations susceptible to the IgAN phenotype that involves interstitial infiltration.

To assess the use of plasma free amino acids (PFAAs) as biomarkers for metabolic disorders, it is essential to identify genetic factors that influence PFAA concentrations. PFAA concentrations were absolutely quantified by liquid chromatography-mass spectrometry using plasma samples from 1338 Japanese individuals, and genome-wide quantitative trait locus (QTL) analysis was performed for the concentrations of 21 PFAAs. We next conducted a conditional QTL analysis using the concentration of each PFAA adjusted by the other 20 PFAAs as covariates to elucidate genetic determinants that influence PFAA concentrations. We identified eight genes that showed a significant association with PFAA concentrations, of which two, SLC7A2 and PKD1L2, were identified. SLC7A2 was associated with the plasma levels of arginine and ornithine, and PKD1L2 with the level of glycine. The significant associations of these two genes were revealed in the conditional QTL analysis, but a significant association between serine and the CPS1 gene disappeared when glycine was used as a covariate. We demonstrated that conditional QTL analysis is useful for determining the metabolic pathways predominantly used for PFAA metabolism. Our findings will help elucidate the physiological roles of genetic components that control the metabolism of amino acids.

Chronic obstructive pulmonary disease (COPD) is a possible risk factor for cardiovascular disease. The association of COPD with the pathogenicity of infection with Chlamydia pneumoniae and Mycoplasma pneumoniae is controversial. We conducted a cross-sectional study to clarify the association between atypical pneumoniae seropositivity and COPD in a general population. We also investigated genetic polymorphisms conferring susceptibility to a pneumonia titer. The study included 9040 Japanese subjects (54 ± 13 years). COPD was defined as a ratio of forced expiratory volume in 1 second to forced vital capacity of less than 70%. Serum levels of IgA and IgG antibodies to C pneumoniae were determined using an enzyme-linked immunoassay, and M pneumoniae seropositivity was assessed by a particle agglutination test. Subjects seropositive for C pneumoniae (26.1%) had a higher prevalence of COPD (seropositive, 5.8%; seronegative, 3.1%; P < 0.001) after adjustment for age, sex, height, weight, and smoking status. The association between M pneumoniae seropositivity (20.4%) and COPD was also significant in covariate-adjusted analysis (P < 0.001). A genome-wide association analysis of the C pneumoniae IgA index identified a susceptible genotype (rs17634369) near the IKZF1 gene, and the seropositive rate of C pneumoniae significantly differed among genotypes (AA, 22.5; AG, 25.3; GG, 29.7%, P < 0.001). On multiple regression analysis, seropositivity for both C pneumoniae (odds ratio = 1.41, P = 0.004) and M pneumoniae (odds ratio = 1.60, P = 0.002) was an independent determinant for COPD, while no direct association was found with the rs17634369 genotype. Seropositivity for both C pneumoniae and M pneumoniae is an independent risk factor for COPD in the general population.

Whole-genome and -exome resequencing using next-generation sequencers is a powerful approach for identifying genomic variations that are associated with diseases. However, systematic strategies for prioritizing causative variants from many candidates to explain the disease phenotype are still far from being established, because the population-specific frequency spectrum of genetic variation has not been characterized. Here, we have collected exomic genetic variation from 1208 Japanese individuals through a collaborative effort, and aggregated the data into a prevailing catalog. In total, we identified 156 622 previously unreported variants. The allele frequencies for the majority (88.8%) were lower than 0.5% in allele frequency and predicted to be functionally deleterious. In addition, we have constructed a Japanese-specific major allele reference genome by which the number of unique mapping of the short reads in our data has increased 0.045% on average. Our results illustrate the importance of constructing an ethnicity-specific reference genome for identifying rare variants. All the collected data were centralized to a newly developed database to serve as useful resources for exploring pathogenic variations. Public access to the database is available at http://www.genome.med.kyoto-u.ac.jp/SnpDB/.

The incidence of high myopia is increasing worldwide with myopic maculopathy, a complication of myopia, often progressing to blindness. Our two-stage genome-wide association study of myopic maculopathy identifies a susceptibility locus at rs11873439 in an intron of CCDC102B (P = 1.77 × 10-12 and Pcorr = 1.61 × 10-10). In contrast, this SNP is not significantly associated with myopia itself. The association between rs11873439 and myopic maculopathy is further confirmed in 2317 highly myopic patients (P = 2.40 × 10-6 and Pcorr = 1.72 × 10-4). CCDC102B is strongly expressed in the retinal pigment epithelium and choroids, where atrophic changes initially occur in myopic maculopathy. The development of myopic maculopathy thus likely exhibits a unique background apart from the development of myopia itself; elucidation of the roles of CCDC102B in myopic maculopathy development may thus provide insights into preventive methods for blindness in patients with high myopia.

Myopia is one of the most common ocular disorders worldwide. Pathological myopia, also called high myopia, comprises 1% to 5% of the general population and is one of the leading causes of legal blindness in developed countries. To identify genetic determinants associated with pathological myopia in Japanese, we conducted a genome-wide association study, analyzing 411,777 SNPs with 830 cases and 1,911 general population controls in a two-stage design (297 cases and 934 controls in the first stage and 533 cases and 977 controls in the second stage). We selected 22 SNPs that showed P-values smaller than 10(-4) in the first stage and tested them for association in the second stage. The meta-analysis combining the first and second stages identified an SNP, rs577948, at chromosome 11q24.1, which was associated with the disease (P = 2.22x10(-7) and OR of 1.37 with 95% confidence interval: 1.21-1.54). Two genes, BLID and LOC399959, were identified within a 200-kb DNA encompassing rs577948. RT-PCR analysis demonstrated that both genes were expressed in human retinal tissue. Our results strongly suggest that the region at 11q24.1 is a novel susceptibility locus for pathological myopia in Japanese.

The cell's movement and morphological change are two interrelated cellular processes. An integrated analysis is needed to explore the relationship between them. However, it has been challenging to investigate them as a whole. The cell's trajectory can be described by its speed, curvature, and torsion. On the other hand, the three-dimensional (3D) cell shape can be studied by using a shape descriptor such as spherical harmonic (SH) descriptor, which is an extension of a Fourier transform in 3D space. We propose a novel method using parallel-transport (PT) to integrate these shape-movement data by using moving frames as the 3D-shape coordinate system. This moving frame is purely determined by the velocity vector. On this moving frame, the movement change will influence the coordinate system for shape analysis. By analyzing the change of the SH coefficients over time in the moving frame, we can observe the relationship between shape and movement. We illustrate the application of our approach using simulated and real datasets in this paper.

Various conjugative plasmids were obtained by exogenous plasmid capture, biparental mating, and/or triparental mating methods from different environmental samples in Japan. Based on phylogenetic analyses of their whole-nucleotide sequences, new IncP/P-1 plasmids that could be classified into novel subgroups were obtained. Mini-replicons of the plasmids were constructed, and each of them was incompatible with at least one of the IncP/P-1 plasmids, although they showed diverse iteron sequences in their oriV regions. There were two large clades of IncP/P-1 plasmids, clade I and II. Plasmids in clade I and II included antibiotic resistance genes. Notably, nucleotide compositions of newly found plasmids exhibited different tendencies compared with those of the previously well-studied IncP/P-1 plasmids. Indeed, the host range of plasmids of clade II was different from that of clade I. Although few PromA plasmids have been reported, the number of plasmids belonging to PromAβ, and -γ subgroups detected in this study was close to that of IncP/P-1 plasmids. The host ranges of PromAγ and PromAδ plasmids were broad and transferred to different and distinct classes of Proteobacteria. Interestingly, PromA plasmids and many IncP/P-1 plasmids do not carry any accessory genes. These findings indicate the presence of "hitherto-unnoticed" conjugative plasmids, including IncP/P-1 or PromA derivative ones in nature. These plasmids would have important roles in the exchange of various genes, including antibiotic resistance genes, among different bacteria in nature. IMPORTANCE Plasmids are known to spread among different bacteria. However, which plasmids spread among environmental samples and in which environments they are present is still poorly understood. This study showed that unidentified conjugative plasmids were present in various environments. Different novel IncP/P-1 plasmids were found, whose host ranges were different from those of known plasmids, showing wide diversity of IncP/P-1 plasmids. PromA plasmids, exhibiting a broad host range, were diversified into several subgroups and widely distributed in varied environments. These findings are important for understanding how bacteria naturally exchange their genes, including antibiotic resistance genes, a growing threat to human health worldwide.

Keratoconus is a common ocular disorder that causes progressive corneal thinning and is the leading indication for corneal transplantation. Central corneal thickness (CCT) is a highly heritable characteristic that is associated with keratoconus. In this two-stage genome-wide association study (GWAS) of CCT, we identified a locus for CCT, namely STON2 rs2371597 (P = 2.32 × 10-13), and confirmed a significant association between STON2 rs2371597 and keratoconus development (P = 0.041). Additionally, strong STON2 expression was observed in mouse corneal epithelial basal cells. We also identified SMAD3 rs12913547 as a susceptibility locus for keratoconus development using predictive analysis with IBM's Watson question answering computer system (P = 0.001). Further GWAS analyses combined with Watson could effectively reveal detailed pathways underlying keratoconus development.

A model that can predict treatment response for a patient with specific baseline characteristics would help decision-making in personalized medicine. The aim of the study is to develop such a model in the treatment of rheumatoid arthritis (RA) patients who receive certolizumab (CTZ) plus methotrexate (MTX) therapy, using individual participant data meta-analysis (IPD-MA).

Group 2 innate lymphoid cells (ILC2s) are critical for the immune response against parasite infection and tissue homeostasis and involved in the pathogenesis of allergy and inflammatory diseases. Although multiple molecules positively regulating ILC2 development and activation have been extensively investigated, the factors limiting their population size and response remain poorly studied. Here, we found that CD45, a membrane-bound tyrosine phosphatase essential for T cell development, negatively regulated ILC2s in a cell-intrinsic manner. ILC2s in CD45-deficient mice exhibited enhanced proliferation and maturation in the bone marrow and hyperactivated phenotypes in the lung with high glycolytic capacity. Furthermore, CD45 signaling suppressed the type 2 inflammatory response by lung ILC2s and alleviated airway inflammation and pulmonary fibrosis. Finally, the interaction with galectin-9 influenced CD45 signaling in ILC2s. These results demonstrate that CD45 is a cell-intrinsic negative regulator of ILC2s and prevents lung inflammation and fibrosis via ILC2s.

Rheumatoid arthritis (RA) is a major cause of adult chronic inflammatory arthritis and a typical complex trait. Although several genetic determinants have been identified, they account for only a part of the genetic susceptibility. We conducted a genome-wide association study of RA in Japanese using 225,079 SNPs genotyped in 990 cases and 1,236 controls from two independent collections (658 cases and 934 controls in collection1; 332 cases and 302 controls in collection2), followed by replication studies in two additional collections (874 cases and 855 controls in collection3; 1,264 cases and 948 controls in collection4). SNPs showing p<0.005 in the first two collections and p<10(-4) by meta-analysis were further genotyped in the latter two collections. A novel risk variant, rs2000811, in intron2 of the myelin basic protein (MBP) at chromosome 18q23 showed strong association with RA (p = 2.7×10(-8), OR 1.23, 95% CI: 1.14-1.32). The transcription of MBP was significantly elevated with the risk allele compared to the alternative allele (p<0.001). We also established by immunohistochemistry that MBP was expressed in the synovial lining layer of RA patients, the main target of inflammation in the disease. Circulating autoantibody against MBP derived from human brain was quantified by ELISA between patients with RA, other connective tissue diseases and healthy controls. As a result, the titer of anti-MBP antibody was markedly higher in plasma of RA patients compared to healthy controls (p<0.001) and patients with other connective tissue disorders (p<0.001). ELISA experiment using citrullinated recombinant MBP revealed that a large fraction of anti-MBP antibody in RA patients recognized citrullinated MBP. This is the first report of a genetic study in RA implicating MBP as a potential autoantigen and its involvement in pathogenesis of the disease.

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes--to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

One of the major issues in genome-wide association studies is to solve the missing heritability problem. While considering epistatic interactions among multiple SNPs may contribute to solving this problem, existing software cannot detect statistically significant high-order interactions. We propose software named LAMPLINK, which employs a cutting-edge method to enumerate statistically significant SNP combinations from genome-wide case-control data. LAMPLINK is implemented as a set of additional functions to PLINK, and hence existing procedures with PLINK can be applicable. Applied to the 1000 Genomes Project data, LAMPLINK detected a combination of five SNPs that are statistically significantly accumulated in the Japanese population.

Quantitative high-throughput screenings (qHTSs) for genotoxicity are conducted as part of comprehensive toxicology screening projects. The most widely used method is to compare the dose-response data of a wild-type and DNA repair gene knockout mutants, using model-fitting to the Hill equation (HE). However, this method performs poorly when the observed viability does not fit the equation well, as frequently happens in qHTS. More capable methods must be developed for qHTS where large data variations are unavoidable. In this study, we applied an isotonic regression (IR) method and compared its performance with HE under multiple data conditions. When dose-response data were suitable to draw HE curves with upper and lower asymptotes and experimental random errors were small, HE was better than IR, but when random errors were big, there was no difference between HE and IR. However, when the drawn curves did not have two asymptotes, IR showed better performance (p < 0.05, exact paired Wilcoxon test) with higher specificity (65% in HE vs. 96% in IR). In summary, IR performed similarly to HE when dose-response data were optimal, whereas IR clearly performed better in suboptimal conditions. These findings indicate that IR would be useful in qHTS for comparing dose-response data.