Search

Abstract The genetic architecture of white matter lesions (WMLs) in Asian populations has not been well-characterized. Here, we performed a genome-wide association study (GWAS) to identify loci associated with the WML volume. Brain MRI and DNA samples were collected from 9479 participants in the Japan Prospective Studies Collaboration for Aging and Dementia (JPSC-AD). The GWAS confirmed three known WML-associated loci (SH3PXD2A, GFAP, and TRIM47). The lead variant of GFAP was a common missense variant (p.D295N) in East Asians. Meta-GWAS using the publicly available summary statistics of UK Biobank identified one previously unreported locus 6q23.2 (SLC2A12). Integration with expression quantitative trait locus data implied the newly identified locus affects SLC2A12 expression. The effect sizes of 20 lead variants at the WML-associated loci were moderately correlated between JPSC-AD and UK Biobank. These results indicate that the alteration in GFAP protein caused by the common missense variant in East Asians influences the WML volume.

Abstract Pelvic organ prolapse (POP) affects approximately 40% of elderly women, characterized by the descent of the pelvic organs into the vaginal cavity. Here we present the results of a genome-wide association study (GWAS) for susceptibility to POP comprising 771 cases and 76,625 controls in the Japanese population. We identified a significant association of WT1 locus with POP in the Japanese population; rs10742277; odds ratio (OR) = 1.48, 95% confidence interval (CI), 1.29–1.68, P = 6.72 × 10− 9. Subsequent cross-ancestry GWAS meta-analysis combining the Japanese data and previously reported European data, including 28,857 cases and 622,916 controls, identified FGFR2 locus as a novel susceptibility locus to POP (rs7072877; OR = 1.06, 95% CI, 1.04–1.08, P = 4.11 × 10− 8). We also observed consistent directions of the effects for 21 out of 24 European GWAS derived loci (binomial test P = 2.8 × 10− 4), indicating that most of susceptibility loci for POP are shared across the Japanese and European populations.

Abstract Little is known about the prospective association between autosomal mosaic chromosomal alterations (mCAs), a group of large-scale somatic mutations on autosomes, and bladder cancer. Here we utilized data from 99,877 participants who were free of physician-diagnosed cancer at baseline (2004–2008) of the China Kadoorie Biobank to estimate the associations between autosomal mCAs and bladder cancer (ICD-10: C67). A total of 2874 autosomal mCAs events among 2612 carriers (2.6%) were detected. After a median follow-up of 12.4 years, we discovered that participants with all autosomal mCAs exhibited higher risks of bladder cancer, with a multivariable-adjusted hazard ratio (HR) (95% confidence interval [CI]) of 2.60 (1.44, 4.70). The estimate of such association was even stronger for mosaic loss events (HR [95% CI]: 6.68 [2.92, 15.30]), while it was not significant for CN-LOH events. Both expanded (cell fraction ≥10%) and non-expanded autosomal mCAs, as well as mosaic loss, were associated with increased risks of bladder cancer. Of interest, physical activity (PA) significantly modified the associations of autosomal mCAs and mosaic loss (P interaction = 0.038 and 0.012, respectively) with bladder cancer. The increased risks of bladder cancer were only observed with mCAs and mosaic loss among participants with a lower level of PA (HR [95% CI]: 5.11 [2.36, 11.09] and 16.30 [6.06, 43.81]), but not among participants with a higher level of PA. Our findings suggest that peripheral leukocyte autosomal mCAs may represent a novel risk factor for bladder cancer, and PA may serve as a potential intervention target for mCAs carriers.

Abstract Sarcopenia is a common skeletal muscle disease in older people. Lower limb muscle strength is a good predictive value for sarcopenia; however, little is known about its genetic components. Here, we conducted a genome-wide association study (GWAS) for knee extension strength in a total of 3452 Japanese aged 60 years or older from two independent cohorts. We identified a significant locus, rs10749438 which is an intronic variant in TACC2 (transforming acidic coiled-coil-containing 2) (P = 4.2 × 10−8). TACC2, encoding a cytoskeleton-related protein, is highly expressed in skeletal muscle, and is reported as a target of myotonic dystrophy 1-associated splicing alterations. These suggest that changes in TACC2 expression are associated with variations in muscle strength in older people. The association was consistently observed in young and middle-aged subjects. Our findings would shed light on genetic components of lower limb muscle strength and indicate TACC2 as a potential therapeutic target for sarcopenia.

Abstract Short- and long-read sequencing technologies are routinely used to detect DNA variants, including SNVs, indels, and structural variations (SVs). However, the differences in the quality and quantity of variants detected between short- and long-read data are not fully understood. In this study, we comprehensively evaluated the variant calling performance of short- and long-read-based SNV, indel, and SV detection algorithms (6 for SNVs, 12 for indels, and 13 for SVs) using a novel evaluation framework incorporating manual visual inspection. The results showed that indel-insertion calls greater than 10 bp were poorly detected by short-read-based detection algorithms compared to long-read-based algorithms; however, the recall and precision of SNV and indel-deletion detection were similar between short- and long-read data. The recall of SV detection with short-read-based algorithms was significantly lower in repetitive regions, especially for small- to intermediate-sized SVs, than that detected with long-read-based algorithms. In contrast, the recall and precision of SV detection in nonrepetitive regions were similar between short- and long-read data. These findings suggest the need for refined strategies, such as incorporating multiple variant detection algorithms, to generate a more complete set of variants using short-read data.

Abstract Here we report the largest Asian genome-wide association study (GWAS) for systemic sclerosis performed to date, based on data from Japanese subjects and comprising of 1428 cases and 112,599 controls. The lead SNP is in the FCGR/FCRL region, which shows a penetrating association in the Asian population, while a complete linkage disequilibrium SNP, rs10917688, is found in a cis-regulatory element for IRF8. IRF8 is also a significant locus in European GWAS for systemic sclerosis, but rs10917688 only shows an association in the presence of the risk allele of IRF8 in the Japanese population. Further analysis shows that rs10917688 is marked with H3K4me1 in primary B cells. A meta-analysis with a European GWAS detects 30 additional significant loci. Polygenic risk scores constructed with the effect sizes of the meta-analysis suggest the potential portability of genetic associations beyond populations. Prioritizing the top 5% of SNPs of IRF8 binding sites in B cells improves the fitting of the polygenic risk scores, underscoring the roles of B cells and IRF8 in the development of systemic sclerosis. The results also suggest that systemic sclerosis shares a common genetic architecture across populations.

BackgroundIt has been well documented that Takayasu arteritis (TAK) and ulcerative colitis (UC) coexist in the same patients. HLA-B*52 characterizes the co-occurrence, which is one of the common genetic features between these two diseases, indicating shared underlying pathologic mechanisms. Anti-integrin αvβ6 antibody (Ab) is present in sera of UC patients in a highly specific manner. We investigated if there were any associations between anti-integrin αvβ6 Ab and TAK, considering the risk HLA alleles.MethodsA total of 227 Japanese TAK patients were recruited in the current study and their serum samples were subjected to measurement of anti-integrin αvβ6 Ab by ELISA. The clinical information, including the co-occurrence of UC, was collected. The HLA allele carrier status was determined by Luminex or genotype imputation.ResultsThe information about the presence of UC was available for 165 patients, among which eight (4.84%) patients had UC. Anti-integrin αvβ6 antibody was identified in 7 out of 8 TAK subjects with UC (87.5%) while only 5 out of 157 (3.18%) TAK subjects without UC had the antibody (OR 121, p=7.46×10-8). A total of 99 out of 218 (45.4%) patients were HLA-B*52 carriers. There was no significant association between the presence of anti-integrin αvβ6 Ab and HLA-B*52 carrier status in those without UC (OR 2.01, 95% CI 0.33-12.4, p = 0.189).ConclusionsThe prevalence of anti-integrin αvβ6 Ab was high in TAK patients with UC, but not in the absence of concomitant UC. The effect of HLA-B*52 on anti-integrin αvβ6 Ab production would be minimal.

Abstract Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Abstract Prostate cancer (PrCa) is the second most common cancer worldwide in males. While strongly warranted, the prediction of mortality risk due to PrCa, especially before its development, is challenging. Here, we address this issue by maximizing the statistical power of genetic data with multi-ancestry meta-analysis and focusing on binding sites of the androgen receptor (AR), which has a critical role in PrCa. Taking advantage of large Japanese samples ever, a multi-ancestry meta-analysis comprising more than 300,000 subjects in total identifies 9 unreported loci including ZFHX3, a tumor suppressor gene, and successfully narrows down the statistically finemapped variants compared to European-only studies, and these variants strongly enrich in AR binding sites. A polygenic risk scores (PRS) analysis restricting to statistically finemapped variants in AR binding sites shows among cancer-free subjects, individuals with a PRS in the top 10% have a strongly higher risk of the future death of PrCa (HR: 5.57, P = 4.2 × 10−10). Our findings demonstrate the potential utility of leveraging large-scale genetic data and advanced analytical methods in predicting the mortality of PrCa.

Abstract Background Polygenic risk score (PRS) analysis is used to predict disease risk. Although PRS has been shown to have great potential in improving clinical care, PRS accuracy assessment has been mainly focused on European ancestry. This study aimed to develop an accurate genetic risk score for knee osteoarthritis (OA) using a multi-population PRS and leveraging a multi-trait PRS in the Japanese population. Methods We calculated PRS using PRS-CS-auto, derived from genome-wide association study (GWAS) summary statistics for knee OA in the Japanese population (same ancestry) and multi-population. We further identified risk factor traits for which PRS could predict knee OA and subsequently developed an integrated PRS based on multi-trait analysis of GWAS (MTAG), including genetically correlated risk traits. PRS performance was evaluated in participants of the Nagahama cohort study who underwent radiographic evaluation of the knees (n = 3,279). PRSs were incorporated into knee OA integrated risk models along with clinical risk factors. Results A total of 2,852 genotyped individuals were included in the PRS analysis. The PRS based on Japanese knee OA GWAS was not associated with knee OA (p = 0.228). In contrast, PRS based on multi-population knee OA GWAS showed a significant association with knee OA (p = 6.7 × 10−5, odds ratio (OR) per standard deviation = 1.19), whereas PRS based on MTAG of multi-population knee OA, along with risk factor traits such as body mass index GWAS, displayed an even stronger association with knee OA (p = 5.4 × 10−7, OR = 1.24). Incorporating this PRS into traditional risk factors improved the predictive ability of knee OA (area under the curve, 74.4% to 74.7%; p = 0.029). Conclusions This study showed that multi-trait PRS based on MTAG, combined with traditional risk factors, and using large sample size multi-population GWAS, significantly improved predictive accuracy for knee OA in the Japanese population, even when the sample size of GWAS of the same ancestry was small. To the best of our knowledge, this is the first study to show a statistically significant association between the PRS and knee OA in a non-European population. Trial registration No. C278.

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E–11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

ABSTRACT Adolescent idiopathic scoliosis (AIS) is the most common form of pediatric musculoskeletal disorder. Observational studies have pointed to several risk factors for AIS, but almost no evidence exists to support their causal association with AIS. Here, we applied Mendelian randomization (MR), known to limit bias from confounding and reverse causation, to investigate causal associations between body composition and puberty‐related exposures and AIS risk in Europeans and Asians. For our two‐sample MR studies, we used single nucleotide polymorphisms (SNPs) associated with body mass index (BMI), waist‐hip ratio, lean mass, childhood obesity, bone mineral density (BMD), 25‐hydroxyvitamin D (25OHD), age at menarche, and pubertal growth in large European genome‐wide association studies (GWAS), and with adult osteoporosis risk and age of menarche in Biobank Japan. We extracted estimates of the aforementioned SNPs on AIS risk from the European or Asian subsets of the largest multiancestry AIS GWAS (N = 7956 cases/88,459 controls). The results of our inverse variance‐weighted (IVW) MR estimates suggest no causal association between the aforementioned risk factors and risk of AIS. Pleiotropy‐sensitive MR methods yielded similar results. However, restricting our analysis to European females with AIS, we observed a causal association between estimated BMD and the risk of AIS (IVW odds ratio for AIS = 0.1, 95% confidence interval 0.01 to 0.7, p = 0.02 per SD increase in estimated BMD), but this association was no longer significant after adjusting for BMI, body fat mass, and 25OHD and remained significant after adjusting for age at menarche in multivariable MR. In conclusion, we demonstrated a protective causal effect of BMD on AIS risk in females of European ancestry, but this effect was modified by BMI, body fat mass, and 25OHD levels. Future MR studies using larger AIS GWAS are needed to investigate small effects of the aforementioned exposures on AIS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Abstract Background Rheumatoid arthritis (RA) occurs across the globe in different ethnic populations. Most RA patients harbor anti-modified protein antibodies (AMPA); however, it is unclear whether differences exist in autoantibody responses at different geographic locations and between different ethnic groups, which could provide new clues regarding factors underlying autoantibody development. We therefore investigated AMPA prevalence and association with HLA DRB1 alleles and smoking in four ethnically diverse populations on four different continents. Methods Anti-carbamylated (anti-CarP), anti-malondialdehyde acetaldehyde (anti-MAA), and anti-acetylated protein antibodies (anti-AcVim) IgG were determined in anti-citrullinated protein antibody-positive Dutch (NL, n = 103), Japanese (JP, n = 174), First Nations Peoples in Canada (FN, n = 100), and black South African (SA, n = 67) RA patients. Ethnicity-matched local healthy controls were used to calculate cut-offs. Risk factors associated with AMPA seropositivity in each cohort were identified using logistic regression. Results Median AMPA levels were higher in First Nations Peoples in Canada and especially South African patients, as reflected by percentage seropositivity: NL, JP, FN, and SA: anti-CarP: 47%, 43%, 58%, and 76% (p

A genome-wide association study in the BioBank Japan cohort reveals new locus, BET1L, which is associated with survival time in a Japanese population.

Abstract Background Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery. Results To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3–5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism. Conclusions Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Asthma is a heterogeneous, complex syndrome that arises in individuals with various genetic and exposure variations. Here, the authors show that disease comorbidity patterns can serve as a surrogate for these variations, and identify asthma endotypes distinguished by comorbidity patterns, asthma risk loci, gene expression, and health-related phenotypes.

Abstract Copy number (CN) polymorphisms of complement C4 play distinct roles in many conditions, including immune-mediated diseases. We investigated the association of C4 CN with systemic sclerosis (SSc) risk. Imputed total C4, C4A, C4B, and HERV-K CN were analyzed in 26,633 individuals and validated in an independent cohort. Our results showed that higher C4 CN confers protection to SSc, and deviations from CN parity of C4A and C4B augmented risk. The protection contributed per copy of C4A and C4B differed by sex. Stronger protection was afforded by C4A in men and by C4B in women. C4 CN correlated well with its gene expression and serum protein levels, and less C4 was detected for both in SSc patients. Conditioned analysis suggests that C4 genetics strongly contributes to the SSc association within the major histocompatibility complex locus and highlights classical alleles and amino acid variants of HLA-DRB1 and HLA-DPB1 as C4-independent signals.

Ossification of the posterior longitudinal ligament of the spine (OPLL) is an intractable disease leading to severe neurological deficits. Its etiology and pathogenesis are primarily unknown. The relationship between OPLL and comorbidities, especially type 2 diabetes (T2D) and high body mass index (BMI), has been the focus of attention; however, no trait has been proven to have a causal relationship. We conducted a meta-analysis of genome-wide association studies (GWASs) using 22,016 Japanese individuals and identified 14 significant loci, 8 of which were previously unreported. We then conducted a gene-based association analysis and a transcriptome-wide Mendelian randomization approach and identified three candidate genes for each. Partitioning heritability enrichment analyses observed significant enrichment of the polygenic signals in the active enhancers of the connective/bone cell group, especially H3K27ac in chondrogenic differentiation cells, as well as the immune/hematopoietic cell group. Single-cell RNA sequencing of Achilles tendon cells from a mouse Achilles tendon ossification model confirmed the expression of genes in GWAS and post-GWAS analyses in mesenchymal and immune cells. Genetic correlations with 96 complex traits showed positive correlations with T2D and BMI and a negative correlation with cerebral aneurysm. Mendelian randomization analysis demonstrated a significant causal effect of increased BMI and high bone mineral density on OPLL. We evaluated the clinical images in detail and classified OPLL into cervical, thoracic, and the other types. GWAS subanalyses identified subtype-specific signals. A polygenic risk score for BMI demonstrated that the effect of BMI was particularly strong in thoracic OPLL. Our study provides genetic insight into the etiology and pathogenesis of OPLL and is expected to serve as a basis for future treatment development.

IntroductionAdolescent idiopathic scoliosis (AIS) is a disorder with a three-dimensional spinal deformity and is a common disease affecting 1-5% of adolescents. AIS is also known as a complex disease involved in environmental and genetic factors. A relation between AIS and body mass index (BMI) has been epidemiologically and genetically suggested. However, the causal relationship between AIS and BMI remains to be elucidated.Material and methodsMendelian randomization (MR) analysis was performed using summary statistics from genome-wide association studies (GWASs) of AIS (Japanese cohort, 5,327 cases, 73,884 controls; US cohort: 1,468 cases, 20,158 controls) and BMI (Biobank Japan: 173430 individual; meta-analysis of genetic investigation of anthropometric traits and UK Biobank: 806334 individuals; European Children cohort: 39620 individuals; Population Architecture using Genomics and Epidemiology: 49335 individuals). In MR analyses evaluating the effect of BMI on AIS, the association between BMI and AIS summary statistics was evaluated using the inverse-variance weighted (IVW) method, weighted median method, and Egger regression (MR-Egger) methods in Japanese.ResultsSignificant causality of genetically decreased BMI on risk of AIS was estimated: IVW method (Estimate (beta) [SE] = -0.56 [0.16], p = 1.8 × 10-3), weighted median method (beta = -0.56 [0.18], p = 8.5 × 10-3) and MR-Egger method (beta = -1.50 [0.43], p = 4.7 × 10-3), respectively. Consistent results were also observed when using the US AIS summary statistic in three MR methods; however, no significant causality was observed when evaluating the effect of AIS on BMI.ConclusionsOur Mendelian randomization analysis using large studies of AIS and GWAS for BMI summary statistics revealed that genetic variants contributing to low BMI have a causal effect on the onset of AIS. This result was consistent with those of epidemiological studies and would contribute to the early detection of AIS.

Rationale & Objective: Chronic kidney disease, defined by albuminuria and/or reduced estimated glomerular filtration rate (eGFR), has been reported to be associated with brain atrophy and/or higher white matter lesion volume (WMLV), but there are few large-scale population-based studies assessing this issue. This study aimed to examine the associations between the urinary albumin-creatinine ratio (UACR) and eGFR levels and brain atrophy and WMLV in a large-scale community-dwelling older population of Japanese. Study Design: Population-based cross-sectional study. Setting & Participants: A total of 8,630 dementia-free community-dwelling Japanese aged greater than or equal to 65 years underwent brain magnetic resonance imaging scanning and screening examination of health status in 2016-2018. Exposures: UACR and eGFR levels. Outcomes: The total brain volume (TBV)-to-intracranial volume (ICV) ratio (TBV/ICV), the regional brain volume-to-TBV ratio, and the WMLV-to-ICV ratio (WMLV/ICV). Analytical Approach: The associations of UACR and eGFR levels with the TBV/ICV, the regional brain volume-to-TBV ratio, and the WMLV/ICV were assessed by using an analysis of covariance. Results: Higher UACR levels were significantly associated with lower TBV/ICV and higher geometric mean values of the WMLV/ICV (P for trend = 0.009 and

Abstract Background Relapsing polychondritis (RP) is a rare inflammatory disease characterized by recurrent inflammation and destruction of cartilaginous tissues. RP has characteristics of autoimmune disease and some reports have noted co-occurrence with autoimmune thyroid disease (AITD), consisting of Graves’ disease (GD) and Hashimoto thyroiditis (HT). However, there have been no detailed studies on the co-occurrence of RP and AITD. In this study, we aimed to determine whether patients with RP tend to be complicated with AITD. We also analyzed the clinical and genetic profiles of patients in whom these diseases co-occur. Methods We recruited 117 patients with RP and reviewed their medical records. Furthermore, we genotyped Human Leucocyte Antigen (HLA)-A, B Cw, DRB1, DQB1, and DPB1 alleles for 93 of the 117 patients. The prevalence of AITD among the patients with RP was compared with that among the general Japanese population. We also analyzed the clinical and genetic features of the patients with both RP and AITD. Results The prevalence of GD among the patients with RP was 4.3% (5 among 117 patients), significantly higher than that among Japanese (0.11%) (p = 2.44 × 10–7, binomial test). RP patients with GD tended to have nasal involvement (p = 0.023) (odds ratio (OR) 2.58) and HLA-DPB1*02:02 (p = 0.035, OR 10.41). We did not find significant enrichment of HT in patients with RP. Conclusions Patients with RP appear to be at elevated risk of GD. Nasal involvement and HLA-DPB1*02:02 characterize the subset of RP patients with GD, which may guide attempts to characterize a distinct subtype of RP for precision medicine.

Abstract Aims/Introduction We aimed to replicate a new diabetes subclassification based on objective clinical information at admission in a diabetes educational inpatient program. We also assessed the educational outcomes for each cluster. Methods We included diabetes patients who participated in the educational inpatient program during 2009–2020 and had sufficient clinical information for the cluster analysis. We applied a data‐driven clustering method proposed in a previous study and further evaluated the clinical characteristics of each cluster. We investigated the association between the clusters and changes in hemoglobin A1c level from the start of the education program. We also assessed the risk of re‐admission for the educational program. Results We divided a total of 651 patients into five clusters. Their clinical characteristics followed the same pattern as in previous studies. The intercluster ranking of the cluster center coordinates showed strong correlation coefficients with those of the previous studies (mean ρ = 0.88). Patients classified as severe insulin‐resistant diabetes (cluster 3) showed a more pronounced progression of renal dysfunction than patients classified as the other clusters. The patients classified as severe insulin‐deficient diabetes (cluster 2) had the highest rate of reduction in hemoglobin A1c level from the start of the program (P

Abstract Background The management of anxiety and depression symptoms in rheumatoid arthritis (RA) patients is vital. Previous study findings on this topic are conflicting, and the topic remains to be thoroughly investigated. This study aimed to clarify the association of RA disease activity with anxiety and depression symptoms after controlling for physical disability, pain, and medication. Methods We conducted a cross-sectional study of RA patients from the XXX Rheumatoid Arthritis Management Alliance cohort. We assessed patients using the Disease Activity Score (DAS28), Health Assessment Questionnaire Disability Index (HAQ-DI), and Hospital Anxiety and Depression Scale (HADS). Anxiety and depression symptoms were defined by a HADS score ≥ 8. We analyzed the data using multivariable logistic regression analyses. Results Of 517 participants, 17.6% had anxiety symptoms and 27.7% had depression symptoms. The multivariable logistic regression analysis demonstrated that DAS28 was not independently associated with anxiety symptoms (odds ratio [OR] [95% confidence interval; CI] 0.93 [0.48–1.78]; p = 0.82) and depression symptoms (OR [95% CI] 1.45 [0.81–2.61]; p = 0.22). However, DAS28 patient global assessment (PtGA) severity was associated with anxiety symptoms (OR [95% CI] 1.15 [1.02–1.29]; p = 0.03) and depression symptoms (OR [95% CI] 1.21 [1.09–1.35]; p

Abstract HLA-B*52 is an established genetic factor in Takayasu arteritis (TAK). Recently, single nucleotide polymorphisms (SNPs) in IL12B (rs6871626) and LILRA3 (rs103294) were newly identified as non-HLA susceptibility loci in TAK. Here, we examined how these SNPs contribute to clinical characteristics and vascular damage in TAK. We retrospectively reviewed the medical records of 99 TAK patients enrolled in our previous genome-wide association study, and whose genotypes for IL12B rs6871626, LILRA3 rs103294, and HLA-B*52 were available. Incidence of aortic regurgitation (AR) was significantly associated with the A allele (risk allele) of IL12B rs6871626 (CC 42%, AC 61%, AA 81%; p = 0.0052; odds ratio [OR] 2.45), as well as with the incidence of hypertension (p = 0.049; OR 1.82) and the proportion of patients who underwent aortic valve replacement (p = 0.023; OR 3.64). Regarding vascular damage, there was positive correlation between the Takayasu Arteritis Damage Score and the A allele of IL12B rs6871626 (CC 3.42 ± 2.71, AC 4.06 ± 3.25, AA 6.00 ± 2.81; p = 0.0035; β = 1.35) and between the Vasculitis Damage Index and the A allele (CC 3.47 ± 1.98, AC 4.33 ± 2.40, AA 5.37 ± 2.22; p = 0.0054; β = 0.96). Contrarily, no correlation was found between LILRA3 rs103294 and vascular damage. In the present study, IL12B rs6871626 was associated with vascular damage in TAK, whereas LILRA3 rs103294 was not. Genotyping of IL12B rs6871626 may help to identify patients at risk of disease progression.

Abstract Retinitis pigmentosa (RP) is a heterogenous hereditary disorder leading to blindness. Despite using next-generation sequencing technologies, causal variants in about 60% of RP cases remain unknown. The heterogeneous genetic inheritance pattern makes it difficult to pinpoint causal variants. Besides, rare penetrating variants are hardly observed in general case–control studies. Thus, a family-based analysis, specifically in a consanguineous family, is a clinically and genetically valuable approach for RP. We analyzed a Japanese consanguineous family with a member suffering from RP with a typical autosomal recessive pattern. We sequenced five direct descendants and spouse using Whole-exome sequencing (WES) and Whole-genome sequencing (WGS). We identified a homozygous pathogenic missense variant in CNGA1 (NM_000087.3, c.839G > A, p.Arg280His) in the proband, while we found no homozygous genotypes in the other family members. CNGA1 was previously reported to be associated with RP. We confirmed the genotypes by the Sanger sequencing. Additionally, we assessed the homozygous genotype in the proband for the possibility of a founder mutation using homozygosity analysis. Our results suggested the two copies of the variant derived from a founder mutation. In conclusion, we found the homozygotes for c.839G > A in CNGA1 as causal for RP.

Abstract Background The anti-cyclic citrullinated peptide (CCP) antibody is a diagnostic biomarker of rheumatoid arthritis (RA). However, some non-RA connective tissue disease (CTD) patients also test positive for the anti-CCP antibody and, thus, may ultimately develop RA. We retrospectively investigated whether anti-CCP-positive non-RA CTD patients developed RA and attempted to identify factors that may differentiate RA-overlapping CTD from pure CTD. Methods In total, 842 CTD patients with a primary diagnosis that was not RA were selected from our CTD database as of December 2012. Anti-CCP antibody titers were obtained from a retrospective chart review or measured using stored sera. RA was diagnosed according to the 1987 revised American College of Rheumatology classification criteria. Thirty-three anti-CCP-positive non-RA CTD patients were retrospectively followed up for the development of RA. Bone erosions on the hands and feet were assessed by X-ray. Citrullination dependency was evaluated by an in-house ELISA, the HLA-DRB1 allele was typed, and the results obtained were then compared between RA-overlapping and non-RA anti-CCP-positive CTD patients. Results Two out of 33 anti-CCP-positive CTD patients (6.1%) developed RA during a mean follow-up period of 8.9 years. X-rays were examined in 27 out of the 33 patients, and only one (3.7%) showed bone erosions. The frequency of the HLA-DRB1 shared epitope (SE) and anti-CCP antibody titers were both significantly higher in anti-CCP-positive RA-overlapping CTD patients than in anti-CCP-positive non-RA CTD patients, while no significant differences were observed in citrullination dependency. Conclusions Anti-CCP-positive non-RA CTD patients rarely developed RA. HLA-DRB1 SE and anti-CCP antibody titers may facilitate the differentiation of RA-overlapping CTD from anti-CCP-positive non-RA CTD.

Abstract Takayasu arteritis (TAK) affects the aorta and its primary branches, mainly in young women. In its advanced stages, it can cause severe complications, such as cerebral infarction, impaired vision, and valvular heart diseases. In the aortic tissue of TAK, there is increased infiltration of cytotoxic lymphocytes, such as natural killer (NK) cells and CD8+T cells, and enhanced expression of accessory molecules, such as major histocompatibility complex (MHC) and MHC class I chain-related gene (MIC) family. Genome-wide association studies on TAK have identified susceptibility genes, such as IL-12p40, MICA, MICB, leukocyte immunoglobulin-like receptor A3 (LILRA3), and LILRB3. Other studies have also shown their involvement in the pathophysiology of TAK. In addition, we reported the importance of NK cells by enhancer enrichment analysis. These results suggest that the gene polymorphisms that potentially upregulate the expression of cytokines and accessory molecules, which contribute to the activation of cytotoxic lymphocytes, are associated with the development of TAK. Based on these results, new molecular targeted therapies look promising.

Background: The underlying pathology of inguinal hernia is still not fully known; thus, further investigations of genetic backgrounds is needed. Here, we aimed to identify genetic factors attributing to inguinal hernias and explore the polygenic architecture of which some components are population-specific, while others are more common among populations. Methods: We performed a genome-wide association study (GWAS) on subjects with inguinal hernias using BioBank Japan (BBJ) data with 1,983 cases and 172,507 controls, followed by a trans-ethnic meta-analysis with UK Biobank (UKBB) data. We performed downstream analyses in order to identify the mechanisms underlying inguinal hernias supported by genetic findings. Findings: We identified a locus closest to ELN, which encodes elastin, at the GWAS significant level. The trans-ethnic meta-analysis revealed 23 additional significant loci, including five loci newly identified not significant in BBJ or UKBB GWAS: TGFB2, RNA5SP214/VGLL2, LOC646588, HMCN2, and ATP5F1CP1/CDKN3. Downstream analyses revealed the overlap of GWAS significant signals in extracellular components, including elastin fiber formation. We also found a highly shared polygenic architecture across different populations (trans-ethnic genetic-effect correlation = 0•77, standard error = 0•26) and population-specific lead variants in ELN, indicating the critical role of elastin in inguinal hernias. Interpretation: We identified a significant locus of the ELN gene in the Japanese population and five additional loci across different populations. Downstream analyses revealed highly shared genetic architectures across populations and highlighted the important roles of extracellular components in the development of inguinal hernias. These findings deepen our understanding of the mechanisms underlying inguinal hernia. Funding: The Japan Agency for Medical Research and Development (AMED) (Grant Number: JP19km0605001)

Kidney stones form in the presence of overabundance of crystal-forming substances such as Ca2+ and oxalate. Here, the authors report genome-wide association analyses for kidney stone disease, report seven previously unknown loci and find that some of these loci also associate with Ca2+ concentration and excretion.

Mosaic loss of chromosome Y (mLOY) is associated with age and smoking but also genetic factors play a role. Here, Terao et al. perform GWAS for mLOY in 95,380 Japanese men and identify 46 loci that overlap with hematopoietic stem cell enhancers and transcription factor binding sites critical for hematopoiesis.

Adolescent idiopathic scoliosis (AIS) is a common pediatric disease leading to spinal deformities. Here, the authors report GWAS followed by genome-wide meta-analysis in up to 79,211 Japanese individuals, identifying 20 genetic loci for AIS, 14 of which were previously unreported, and perform in vitro validation for rs1978060.

Genetic variants in RP1 can cause hereditary retinal degeneration (HRD). Here, in a genomic screen of 331 Japanese HRD patients, the authors identify a near-polymorphic RP1 variant that causes Mendelian HRD in trans with an Alu insertion and otherwise is associated with HRD according to a complex model of inheritance.

Abstract Background Structural variations (SVs) or copy number variations (CNVs) greatly impact the functions of the genes encoded in the genome and are responsible for diverse human diseases. Although a number of existing SV detection algorithms can detect many types of SVs using whole genome sequencing (WGS) data, no single algorithm can call every type of SVs with high precision and high recall. Results We comprehensively evaluate the performance of 69 existing SV detection algorithms using multiple simulated and real WGS datasets. The results highlight a subset of algorithms that accurately call SVs depending on specific types and size ranges of the SVs and that accurately determine breakpoints, sizes, and genotypes of the SVs. We enumerate potential good algorithms for each SV category, among which GRIDSS, Lumpy, SVseq2, SoftSV, Manta, and Wham are better algorithms in deletion or duplication categories. To improve the accuracy of SV calling, we systematically evaluate the accuracy of overlapping calls between possible combinations of algorithms for every type and size range of SVs. The results demonstrate that both the precision and recall for overlapping calls vary depending on the combinations of specific algorithms rather than the combinations of methods used in the algorithms. Conclusion These results suggest that careful selection of the algorithms for each type and size range of SVs is required for accurate calling of SVs. The selection of specific pairs of algorithms for overlapping calls promises to effectively improve the SV detection accuracy.

Sequences homologous to human herpesvirus 6 (HHV-6) are integrated within the nuclear genome of about 1% of humans, but it is not clear how this came about. It is also uncertain whether integrated HHV-6 can reactivate into an infectious virus. HHV-6 integrates into telomeres, and this has recently been associated with polymorphisms affecting MOV10L1. MOV10L1 is located on the subtelomere of chromosome 22q (chr22q) and is required to make PIWI-interacting RNAs (piRNAs). As piRNAs block germline integration of transposons, piRNA-mediated repression of HHV-6 integration has been proposed to explain this association. In vitro, recombination of the HHV-6 genome along its terminal direct repeats (DRs) leads to excision from the telomere and viral reactivation, but the expected "solo-DR scar" has not been described in vivo. Here we screened for integrated HHV-6 in 7,485 Japanese subjects using whole-genome sequencing (WGS). Integrated HHV-6 was associated with polymorphisms on chr22q. However, in contrast to prior work, we find that the reported MOV10L1 polymorphism is physically linked to an ancient endogenous HHV-6A variant integrated into the telomere of chr22q in East Asians. Unexpectedly, an HHV-6B variant has also endogenized in chr22q; two endogenous HHV-6 variants at this locus thus account for 72% of all integrated HHV-6 in Japan. We also report human genomes carrying only one portion of the HHV-6B genome, a solo-DR, supporting in vivo excision and possible viral reactivation. Together these results explain the recently-reported association between integrated HHV-6 and MOV10L1/piRNAs, suggest potential exaptation of HHV-6 in its coevolution with human chr22q, and clarify the evolution and risk of reactivation of the only intact (non-retro)viral genome known to be present in human germlines.

Abstract Background A previous study revealed the association between susceptibility to Takayasu arteritis (TAK) and a single nucleotide polymorphism (SNP) rs6871626 located in IL12B, which encodes interleukin (IL)-12p40, a common component of IL-12p70 and IL-23. We investigated the expression of these cytokines in patients with TAK, stratifying them into those with or without the risk allele at the rs6871626 SNP. Methods Plasma levels of IL-12p40, IL-12p70, and IL-23 were quantified in 44 patients with TAK and 19 healthy controls (HCs) by enzyme-linked immunosorbent assays. Monocytes were obtained from 20 patients with TAK and 14 HCs, treated with interferon-γ (IFN-γ) and lipopolysaccharide, and then supernatant cytokines were quantified. In addition, the ratio of IFN-γ+ or IL-17A+ cells to CD4+ T cells was measured by flow cytometric analysis of peripheral blood mononuclear cells. Results The levels of plasma IL-12p40, plasma IL-12p70, and supernatant IL-12p70 were significantly higher in patients with TAK than in HCs, whereas there were no significant differences in the levels of plasma IL-23, supernatant IL-23, or supernatant IL-12p40. The levels of plasma IL-12p70, supernatant IL-12p40, and supernatant IL-12p70 were significantly higher in patients with the risk allele than in those without. The ratio of CD4+IFN-γ+ cells was significantly higher in patients with the risk allele, whereas CD4+IL-17A+ cells showed no differences. Conclusions The rs6871626 SNP in IL12B may influence the increased expression of IL-12p40 and IL-12p70. These enhanced cytokines might play roles in the pathophysiology of TAK.

Abstract 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

Abstract Disease activity of rheumatoid arthritis (RA), evaluated as Disease Activity Score (DAS), is associated with joint destruction. Since joint destruction reflects the history of disease activities, we hypothesized that time-averaged disease activity would better correlate with joint destruction than one-time disease activity. We recruited RA patients in IORRA (n = 557) and KURAMA (n = 204) cohorts, and calculated time-averaged DAS28 to model a modified Sharp/van der Heijde score (SHS). We evaluated the fitting of the model using time-averaged DAS28 among 1000 models in which we randomly picked up one-time DAS28. We also used clinical disease activity index (CDAI) or data in the BeSt study (European population). After conditioning on autoantibody and disease duration, time-averaged DAS28 showed significant improvement of model fitting compared with one-time DAS28 in both cohorts (p = 0.001 and 0.034, respectively). Time-averaged CDAI also showed a better fit. Integration of multiple DAS fit SHS better in the BeSt study. A good fit of time-averaged DAS could be observed using five to six time points of DAS. In conclusion, time-averaged disease activity fits the joint destruction model better than one-time disease activity. Usage of time-averaged disease activity as a covariate would increase the power of studies to identify novel correlates of joint destruction.

Abstract Background Pulmonary arterial hypertension (PAH) is a severe lung disease with only few effective treatments available. Familial cases of PAH are usually recognized as an autosomal dominant disease, but incomplete penetrance of the disease makes it difficult to identify pathogenic variants in accordance with a Mendelian pattern of inheritance. Methods To elucidate the complex genetic basis of PAH, we obtained whole exome- or genome-sequencing data of 17 subjects from 9 families with heritable PAH and applied gene-based association analysis with 9 index patients and 300 PAH-free controls. Results A burden of rare variants in BMPR2 significantly contributed to the risk of the disease (p = 6.0 × 10−8). Eight of nine families carried four previously reported single nucleotide variants and four novel insertion/deletion variants in the gene. One of the novel variants was a large 6.5 kilobase-deletion. In the remaining one family, the patient carried a pathogenic variant in a member of potassium channels, KCNK3, which was the first replicative finding of channelopathy in an Asian population. Conclusions The variety of rare pathogenic variants suggests that gene-based association analysis using genome-wide sequencing data from increased number of samples is essential to tracing the genetic heterogeneity and developing an appropriate panel for genetic testing.

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component. We recently identified a novel SLE susceptibility locus near RASGRP1, which governs the ERK/MAPK kinase cascade and B-/T-cell differentiation and development. However, precise causal RASGRP1 functional variant(s) and their mechanisms of action in SLE pathogenesis remain undefined. Our goal was to fine-map this locus, prioritize genetic variants likely to be functional, experimentally validate their biochemical mechanisms, and determine the contribution of these SNPs to SLE risk. We performed a meta-analysis across six Asian and European cohorts (9,529 cases; 22,462 controls), followed by in silico bioinformatic and epigenetic analyses to prioritize potentially functional SNPs. We experimentally validated the functional significance and mechanism of action of three SNPs in cultured T-cells. Meta-analysis identified 18 genome-wide significant (p < 5 × 10−8) SNPs, mostly concentrated in two haplotype blocks, one intronic and the other intergenic. Epigenetic fine-mapping, allelic, eQTL, and imbalance analyses predicted three transcriptional regulatory regions with four SNPs (rs7170151, rs11631591-rs7173565, and rs9920715) prioritized for functional validation. Luciferase reporter assays indicated significant allele-specific enhancer activity for intronic rs7170151 and rs11631591-rs7173565 in T-lymphoid (Jurkat) cells, but not in HEK293 cells. Following up with EMSA, mass spectrometry, and ChIP-qPCR, we detected allele-dependent interactions between heterogeneous nuclear ribonucleoprotein K (hnRNP-K) and rs11631591. Furthermore, inhibition of hnRNP-K in Jurkat and primary T-cells downregulated RASGRP1 and ERK/MAPK signaling. Comprehensive association, bioinformatics, and epigenetic analyses yielded putative functional variants of RASGRP1, which were experimentally validated. Notably, intronic variant (rs11631591) is located in a cell type-specific enhancer sequence, where its risk allele binds to the hnRNP-K protein and modulates RASGRP1 expression in Jurkat and primary T-cells. As risk allele dosage of rs11631591 correlates with increased RASGRP1 expression and ERK activity, we suggest that this SNP may underlie SLE risk at this locus.

Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10-75) and its proxy position 11 (P = 1.1x10-67), followed by HLA-DRB1-37 (P = 4.5x10-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10-8), HLA-DPB1-35 (P = 9.0x10-16), HLA-DQB1-37 (P = 2.7x10-14), and HLA-B-9 (P = 6.5x10-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10-14; DRB1-13, P = 2.9x10-13; DRB1-30, P = 3.9x10-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.

Rheumatoid arthritis (RA) is an autoimmune disease characterized by chronic inflammation and subsequent proliferation of synovial tissues, which eventually leads to cartilage and bone destruction without effective treatments. Anti-citrullinated cyclic peptide/protein antibody (ACPA) and rheumatoid factor (RF) are two main characteristic autoantibodies found in RA patients and are associated with unfavorable disease outcomes. Although etiologies and causes of the disease have not been fully clarified yet, it is likely that interactive contributions of genetic and environmental factors play a main role in RA pathology. Previous works have demonstrated several genetic and environmental factors as risks of RA development and/or autoantibody productions. Among these, cigarette smoking and HLA-DRB1 are the well-established environmental and genetic risks, respectively. In this narrative review, we provide a recent update on genetic contributions to RA and the environmental risks of RA with a special focus on cigarette smoking and its impacts on RA pathology. We also describe gene−environmental interaction in RA pathogenesis with an emphasis on cigarette smoking and HLA-DRB1.

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.

To determine whether the presence of periodontitis (PD) and Porphyromonas gingivalis (Pg) in the subgingival biofilm associates with the development of rheumatoid arthritis (RA) in treatment naïve preclinical stage of arthritis patients.We conducted a prospective cohort study of 72 consecutive patients with arthralgia who had never been treated with any anti-rheumatic drugs or glucocorticoids. Periodontal status at baseline was assessed by dentists. PD was defined stringently by the maximal probing depth≧4 mm, or by the classification by the 5th European Workshop in Periodontology (EWP) in 2005 using attachment loss. Up to eight plaque samples were obtained from each patient and the presence of Pg was determined by Taqman PCR. The patients were followed up for 2 years and introduction rate of methotrexate (MTX) treatment on the diagnosis of RA was compared in patients with or without PD or Pg.Patients with PD (probing depth≧4mm) had higher arthritis activity (p = 0.02) and higher risk for future introduction of MTX treatment on the diagnosis of RA during the follow up than patients without PD (Hazard ratio 2.68, p = 0.03). Arthritis activity and risk for MTX introduction increased with the severity of PD assessed by EWP, although not statistically significant. On the other hand, presence of Pg was not associated with arthritis activity (p = 0.72) or the risk for MTX introduction (p = 0.45).In treatment naïve arthralgia patients, PD, but not the presence of Pg, associates with arthritis activity and future requirement of MTX treatment on the diagnosis of RA.

To investigate clinical and radiological differences between joint destruction in the wrist and the feet in patients with RA.A cross-sectional clinical study was conducted in an RA cohort at a single institution. Clinical data included age, sex and duration of disease. Laboratory data included sero-positivity for anti-cyclic citrullinated peptide (CCP) antibody and RF. Radiological measurements included Larsen grades and the modified Sharp/van der Heijde method (SHS) for the hands/wrists and the feet. Statistical analyses were performed using the Kruskal-Wallis H-test, a dummy variable linear regression model and multivariate logistic regression analysis with 95% confidence interval and odds ratios.A total of 405 patients were enrolled, and 314 patients were analysed in this study. The duration of disease in the foot-dominant group was significantly less than that in the wrist-dominant group. When patients were subdivided by duration of disease, the Larsen grade of the feet was significantly higher than that of the wrist in the first quadrant subgroup, but this was reversed with increasing duration of disease. Anti-CCP status was a significant predictive factor for joint destruction in the wrist but not in the feet, while RF status was not predictive in either the wrist or the feet.Joint destruction in the feet started earlier than in the wrist, but the latter progresses faster with increasing duration of disease. Anti-CCP status predicts joint destruction in the wrist better than in the feet.

To compare the effect of tocilizumab (TCZ) with other biologic therapies in improving anemia of rheumatoid arthritis (RA) patients.We compared the change of hemoglobin (Hb) levels in a cohort of 147 consecutive RA patients who were treated with biologics for more than 12 weeks. Twenty eight patients were treated with TCZ, and 119 patients were treated with biologics other than TCZ (87 with TNF inhibitors and 32 with abatacept). The change of Hb levels from baseline to week 12 was compared between the TCZ and the non-TCZ groups. We performed univariate and multivariate analyses with adjustment of potential confounders such as baseline characteristics, concomitant treatment, and the clinical response to treatment.Hb levels generally increased after biologic therapies both in the TCZ and the non-TCZ groups. The increase of Hb levels was greater in the TCZ group than in the non-TCZ groups (1.1 g/dL in the TCZ group vs 0.3 g/dL in the non-TCZ group, p = 0.009). Univariate analysis revealed that increase of Hb levels was also significantly associated with lower Hb, higher Low Hemoglobin Density, and higher CRP levels at baseline and greater reduction in the clinical disease activity index. TCZ therapy was significantly associated with the increase of Hb levels even after adjustment for these factors by multivariate analysis (p

Rheumatoid arthritis (RA) is a bone destructive autoimmune disease. Many patients with RA recognize fluctuations of their joint synovitis according to changes of air pressure, but the correlations between them have never been addressed in large-scale association studies. To address this point we recruited large-scale assessments of RA activity in a Japanese population, and performed an association analysis. Here, a total of 23,064 assessments of RA activity from 2,131 patients were obtained from the KURAMA (Kyoto University Rheumatoid Arthritis Management Alliance) database. Detailed correlations between air pressure and joint swelling or tenderness were analyzed separately for each of the 326 patients with more than 20 assessments to regulate intra-patient correlations. Association studies were also performed for seven consecutive days to identify the strongest correlations. Standardized multiple linear regression analysis was performed to evaluate independent influences from other meteorological factors. As a result, components of composite measures for RA disease activity revealed suggestive negative associations with air pressure. The 326 patients displayed significant negative mean correlations between air pressure and swellings or the sum of swellings and tenderness (p = 0.00068 and 0.00011, respectively). Among the seven consecutive days, the most significant mean negative correlations were observed for air pressure three days before evaluations of RA synovitis (p = 1.7 × 10(-7), 0.00027, and 8.3 × 10(-8), for swellings, tenderness and the sum of them, respectively). Standardized multiple linear regression analysis revealed these associations were independent from humidity and temperature. Our findings suggest that air pressure is inversely associated with synovitis in patients with RA.

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8 × 10(-8)), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3' UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1 × 10(-11) in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry.