Mitochondrial disorders have the highest incidence among congenital metabolic disorders characterized by biochemical respiratory chain complex deficiencies. It occurs at a rate of 1 in 5,000 births, and has phenotypic and genetic heterogeneity. Mutations in about 1,500 nuclear encoded mitochondrial proteins may cause mitochondrial dysfunction of energy production and mitochondrial disorders. More than 250 genes that cause mitochondrial disorders have been reported to date. However exact genetic diagnosis for patients still remained largely unknown. To reveal this heterogeneity, we performed comprehensive genomic analyses for 142 patients with childhood-onset mitochondrial respiratory chain complex deficiencies. The approach includes whole mtDNA and exome analyses using high-throughput sequencing, and chromosomal aberration analyses using high-density oligonucleotide arrays. We identified 37 novel mutations in known mitochondrial disease genes and 3 mitochondria-related genes (MRPS23, QRSL1, and PNPLA4) as novel causative genes. We also identified 2 genes known to cause monogenic diseases (MECP2 and TNNI3) and 3 chromosomal aberrations (6q24.3-q25.1, 17p12, and 22q11.21) as causes in this cohort. Our approaches enhance the ability to identify pathogenic gene mutations in patients with biochemically defined mitochondrial respiratory chain complex deficiencies in clinical settings. They also underscore clinical and genetic heterogeneity and will improve patient care of this complex disorder.

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.

The genetic architecture of common traits, including the number, frequency, and effect sizes of inherited variants that contribute to individual risk, has been long debated. Genome-wide association studies have identified scores of common variants associated with type 2 diabetes, but in aggregate, these explain only a fraction of the heritability of this disease. Here, to test the hypothesis that lower-frequency variants explain much of the remainder, the GoT2D and T2D-GENES consortia performed whole-genome sequencing in 2,657 European individuals with and without diabetes, and exome sequencing in 12,940 individuals from five ancestry groups. To increase statistical power, we expanded the sample size via genotyping and imputation in a further 111,548 subjects. Variants associated with type 2 diabetes after sequencing were overwhelmingly common and most fell within regions previously identified by genome-wide association studies. Comprehensive enumeration of sequence variation is necessary to identify functional alleles that provide important clues to disease pathophysiology, but large-scale sequencing does not support the idea that lower-frequency variants have a major role in predisposition to type 2 diabetes.

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.

Despite many years of study into inversions, very little is known about their functional consequences, especially in humans. A common hypothesis is that the selective value of inversions stems in part from their effects on nearby genes, although evidence of this in natural populations is almost nonexistent. Here we present a global analysis of a new 415-kb polymorphic inversion that is among the longest ones found in humans and is the first with clear position effects. This inversion is located in chromosome 19 and has been generated by non-homologous end joining between blocks of transposable elements with low identity. PCR genotyping in 541 individuals from eight different human populations allowed the detection of tag SNPs and inversion genotyping in multiple populations worldwide, showing that the inverted allele is mainly found in East Asia with an average frequency of 4.7%. Interestingly, one of the breakpoints disrupts the transcription factor gene ZNF257, causing a significant reduction in the total expression level of this gene in lymphoblastoid cell lines. RNA-Seq analysis of the effects of this expression change in standard homozygotes and inversion heterozygotes revealed distinct expression patterns that were validated by quantitative RT-PCR. Moreover, we have found a new fusion transcript that is generated exclusively from inverted chromosomes around one of the breakpoints. Finally, by the analysis of the associated nucleotide variation, we have estimated that the inversion was generated ~40,000-50,000 years ago and, while a neutral evolution cannot be ruled out, its current frequencies are more consistent with those expected for a deleterious variant, although no significant association with phenotypic traits has been found so far.

Thicker or thinner central corneas may lead to either overestimation or underestimation of intraocular pressure, which is the most important causal and treatable risk factor for glaucoma. However, the findings on the associations between diabetes, random glucose, and glycated hemoglobin A1c (HbA1c) with central corneal thickness (CCT) are conflicting.

Lipoproteins are major players in the development and progression of atherosclerotic plaques leading to coronary stenosis and myocardial infarction. Epidemiological, genetic and experimental observations have implicated the association of sphingolipids and intermediates of sphingolipid synthesis in atherosclerosis. We aimed to investigate relationships between quantitative changes in serum sphingolipids, the regulation of the metabolism of lipoproteins (LDL, HDL), and endophenotypes of coronary artery disease (CAD).

Mismatch repair (MMR)-deficient or microsatellite instability (MSI) colorectal cancer includes two subtypes; Lynch syndrome and sporadic MSI cancer, both of which generate multiple neoantigens due to unrepaired mutations. Although such patients respond very well to immune checkpoint therapy, their diagnosis can be confused by low quality DNA samples owing to formalin fixation and/or low cancer cell content. Here we prepared high-quality DNA samples from in vitro-cultured cancer spheroids that consisted of the pure cell population. We evaluated their diagnostic power by on-chip electrophoresis, mutational burden assessment, and direct sequencing. Because formalin-fixed paraffin-embedded (FFPE) tissues are widely used as the DNA source, we compared such samples with spheroid DNA. Additionally, we performed immunohistochemistry (IHC) for MMR proteins on spheroids as well as primary tumor sections. Of 111 cases of colorectal cancer patients, we found seven MSI-high cases in which all diagnostic results agreed on spheroid-based assays, whereas the results with the FFPE DNA were less reliable though analyzable. Importantly, there was an MSS case that appeared as MSI by IHC on primary tumor sections. Based on these results, we propose to employ cultured cancer spheroids as the source of both DNA and IHC specimens for more reliable clinical diagnosis.

The role of leptin (a hormone related to fat mass) in cognition remains equivocal. Our objective was to investigate the relationship between circulating leptin concentration and cognition in older adults, accounting for potential confounders. We categorized 1061 community-dwelling older participants ≥60 years (mean ± SD, 70.6 ± 6.4 years; 41.6% female) from the Singapore Kidney Eye Study according to quintiles of leptin concentration (≤2.64; 2.64⁻5.1; 5.2⁻8.6; 8.7⁻17.96; ≥18 ng/mL). Cognition was assessed using the total and domain scores of the Abbreviated Mental Test (AMT). Age, gender, body mass index, mean arterial pressure, smoking, alcohol, education, memory complaint, anxiodepressive disorders, circulating concentrations of 25-hydroxyvitamin D, glycosylated hemoglobin, low-density lipoprotein cholesterol, and estimated glomerular filtration rate were used as potential confounders. Participants within the lowest (Q1) and highest (Q5) leptin quintiles exhibited lower (i.e., worse) mean total AMT scores compared to those within the intermediate quintiles (Q2, Q3, and Q4). Compared to Q3 as the reference, Q1 and Q5 were associated with decreased total AMT score (respectively, β = -0.53 p = 0.018; β = -0.60 p = 0.036). Compared to Q3, Q5 was also associated with decreased subscores on anterograde (β = -0.19 p = 0.020) and retrograde episodic memories (β = -0.18 p = 0.039). We found a non-linear U-shaped relationship between circulating leptin and cognition, with both lower and higher concentrations of leptin being associated with more severe cognitive impairment in community-dwelling older Asians.

Glaucoma is the leading cause of irreversible blindness globally 1 . Despite its gravity, the disease is frequently undiagnosed in the community 2 . Raised intraocular pressure (IOP) is the most important risk factor for primary open-angle glaucoma (POAG)3,4. Here we present a meta-analysis of 139,555 European participants, which identified 112 genomic loci associated with IOP, 68 of which are novel. These loci suggest a strong role for angiopoietin-receptor tyrosine kinase signaling, lipid metabolism, mitochondrial function and developmental processes underlying risk for elevated IOP. In addition, 48 of these loci were nominally associated with glaucoma in an independent cohort, 14 of which were significant at a Bonferroni-corrected threshold. Regression-based glaucoma-prediction models had an area under the receiver operating characteristic curve (AUROC) of 0.76 in US NEIGHBORHOOD study participants and 0.74 in independent glaucoma cases from the UK Biobank. Genetic-prediction models for POAG offer an opportunity to target screening and timely therapy to individuals most at risk.

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

Glycated hemoglobin A1C (HbA1C) level is used as a diagnostic marker for diabetes mellitus and a predictor of diabetes associated complications. Genome-wide association studies have identified genetic variants associated with HbA1C level. Most of these studies have been conducted in populations of European ancestry. Here we report the findings from a meta-analysis of genome-wide association studies of HbA1C levels in 6,682 non-diabetic subjects of Chinese, Malay and South Asian ancestries. We also sought to examine the associations between HbA1C associated SNPs and microvascular complications associated with diabetes mellitus, namely chronic kidney disease and retinopathy. A cluster of 6 SNPs on chromosome 17 showed an association with HbA1C which achieved genome-wide significance in the Malays but not in Chinese and Asian Indians. No other variants achieved genome-wide significance in the individual studies or in the meta-analysis. When we investigated the reproducibility of the findings that emerged from the European studies, six loci out of fifteen were found to be associated with HbA1C with effect sizes similar to those reported in the populations of European ancestry and P-value ≤ 0.05. No convincing associations with chronic kidney disease and retinopathy were identified in this study.

We encountered a family of Japanese descent in which multiple members developed lung cancer. Using whole-exome sequencing, we identified a novel germline mutation in the transmembrane domain of the human epidermal growth factor receptor 2 (HER2) gene (G660D). A novel somatic mutation (V659E) was also detected in the transmembrane domain of HER2 in one of 253 sporadic lung adenocarcinomas. Because the transmembrane domain of HER2 is considered to be responsible for the dimerization and subsequent activation of the HER family and downstream signaling pathways, we performed functional analyses of these HER2 mutants. Mutant HER2 G660D and V659E proteins were more stable than wild-type protein. Both the G660D and V659E mutants activated Akt. In addition, they activated p38, which is thought to promote cell proliferation in lung adenocarcinoma. Our findings strongly suggest that mutations in the transmembrane domain of HER2 may be oncogenic, causing hereditary and sporadic lung adenocarcinomas.

The contribution of smoking to rheumatoid arthritis (RA) is hypothesized to be mediated through formation of anti-citrullinated protein antibodies (ACPA). In RA, however, autoantibodies such as ACPA, rheumatoid factor (RF), and anti-carbamylated protein antibodies (anti-CarP) often occur together, and it is thus unclear whether smoking is specifically associated with some autoantibodies rather than others. We therefore investigated whether smoking is only associated with ACPA or with the presence of multiple RA-related autoantibodies.

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

Previous studies showed an independent association of low ankle-brachial index (ABI) with cognitive impairment. However, the association between low ABI and cognition in the presence of both cerebrovascular disease (CeVD) and neurodegeneration is lacking. We aimed at investigating a) the association of low ABI with markers of CeVD and cortical thickness, and b) whether the association of low ABI with cognition is influenced by these markers.

Hypoxia augments inflammatory responses and osteoclastogenesis by incompletely understood mechanisms. We identified COMMD1 as a cell-intrinsic negative regulator of osteoclastogenesis that is suppressed by hypoxia. In human macrophages, COMMD1 restrained induction of NF-κB signaling and a transcription factor E2F1-dependent metabolic pathway by the cytokine RANKL. Downregulation of COMMD1 protein expression by hypoxia augmented RANKL-induced expression of inflammatory and E2F1 target genes and downstream osteoclastogenesis. E2F1 targets included glycolysis and metabolic genes including CKB that enabled cells to meet metabolic demands in challenging environments, as well as inflammatory cytokine-driven target genes. Expression quantitative trait locus analysis linked increased COMMD1 expression with decreased bone erosion in rheumatoid arthritis. Myeloid deletion of Commd1 resulted in increased osteoclastogenesis in arthritis and inflammatory osteolysis models. These results identify COMMD1 and an E2F-metabolic pathway as key regulators of osteoclastogenic responses under pathological inflammatory conditions and provide a mechanism by which hypoxia augments inflammation and bone destruction.

Anti-centromere antibody (ACA) is one of the classical anti-nuclear antibody (ANA) staining patterns. However, characteristics of ACA in comparison with the other ANA patterns and clinical features of ACA-positive subjects have not been elucidated. Here, we examined all ANA patterns by indirect immunofluorescence for 859 rheumatoid arthritis (RA) patients. Together with the ANA data of 9,575 healthy volunteers, we compared distributions of the ANA levels. ACA was the only ANA that demonstrated a definite bimodal distribution of levels. ACA showed significantly higher levels than the other ANA staining patterns in both RA and healthy population (p < 0.0001). ACA-positivity was associated with old age and was observed more in females. We further recruited another cohort of 3,353 RA patients and confirmed the findings. ACA was also associated with Raynaud's phenomenon (p = 6.8 × 10-11) in RA. As a conclusion, ACA displays a specific ANA staining pattern with a bimodal distribution, and ACA-positive RA may constitute a distinct subset with specific clinical features.

SIX1 and SIX6 are glaucoma susceptibility genes. Previous reports indicate that the single nucleotide polymorphism (SNP) rs33912345 in SIX6 is associated with inferior circumpapillary retinal nerve fibre layer (cpRNFL) thickness (cpRNFLT). Although the region of visual field defect in glaucoma patients is directly related to cpRNFL thinning, a detailed sector analysis has not been performed in genetic association studies. In the present study, we evaluated 26 tagging SNPs in the SIX1/SIX6 locus ±50 kb region in a population of 2,306 Japanese subjects with 4- and 32-sector cpRNFLT analysis. While no SNPs showed a significant association with cpRNFLT in the 4-sectored analysis, the finer 32-sector assessment clearly showed a significant association between rs33912345 in the SIX1/SIX6 locus with inferior cpRNFL thinning at 292.5-303.8° (β = -4.55, P = 3.0 × 10-5). Furthermore, the fine-sectored cpRNFLT analysis indicated that SIX1/SIX6 polymorphisms would affect cpRNFL thinning at 281.3-303.8°, which corresponds to parafoveal scotoma in a visual field test of glaucoma patients.

Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes.