HLA-B27, axial spondyloarthritis and survival.

INTRODUCTION: Ankylosing spondylitis (AS), and carriage of HLA-B27 gene in otherwise healthy individuals, are reportedly associated with increased mortality. We evaluated this hypothesis, using data from both a 35-year AS follow-up study and UK Biobank data. METHODS: In 1985, 363 members of the Swiss AS Patient Society and 806 relatives were screened clinically and then radiographically for AS/axial spondyloarthritis (axSpA). Life expectancy was analysed in 377 axSpA patients having available pelvic radiographs and HLA-B27 status, comparing with matched Swiss population data. Survival in relation to HLA-B27 status in the general population was studied in UK Biobank European-ancestry participants (n=407 480, n=30 419 deaths). RESULTS: AS patients have increased standardised mortality rate (SMR) compared with the general population (1.37, 95% CI 1.11 to 1.62). This increase was significant for HLA-B27-positive AS (SMR 1.38, 95% CI 1.11 to 1.65). Shortened life expectancy was observed among both HLA-B27-positive AS women (SMR 1.77, 95% CI 1.09 to 2.70) and men (SMR 1.31, 95% CI 1.02 to 1.59). Patients with non-radiographic axSpA (nr-axSpA) had significantly lower SMR: 0.44 (95% CI 0.23 to 0.77), compared with the general population. In the UK Biobank European-ancestry population cohort, HLA-B27 carriage was not significantly associated with any change in mortality (HR 1, 95% CI 0.97 to 1.1, p=0.349, adjusted by sex), in either males (HR 1, 95% CI 0.98 to 1.1, p=0.281) or females (HR 0.96, 95% CI 0.9 to 1, p=0.232), and no increase in vascular disease mortality was observed. DISCUSSION: AS patients, but not nr-axSpA patients, have a significantly shortened life expectancy. Increased mortality is particularly significant among women with HLA-B27-positive AS. HLA-B27 carriage in the European-ancestry general population does not influence survival, or the risk of death due to vascular disease.

Correction: Stoll et al. Impact of HLA-B27 and Disease Status on the Gut Microbiome of the Offspring of Ankylosing Spondylitis Patients. Children 2022, 9, 569.

The authors wish to make the following correction to this paper [...].

Recurrence of axial spondyloarthritis among first-degree relatives in a prospective 35-year-follow-up family study.

OBJECTIVE: The lifetime recurrence rate (RR) of axial spondyloarthritis (axSpA) among first-degree relatives (FDR) and the effect of proband's gender, HLA-B27 and radiographic status is unclear. Our 35-year-follow-up family study has enabled these issues to be addressed. METHODS: In 1985, 363 ankylosing spondylitis (AS) probands (members of the Swiss AS Patient Society) and 806 FDR recruited into the study, completed questionnaires regarding axSpA manifestations, underwent a physical examination and most also underwent pelvic radiography and HLA-B27 typing. At follow-up in 2018-2019, of the former participants whose current addresses could be retrieved, 162 had died and 485 (125 patients with AS plus 360 FDR) completed a postal questionnaire. RESULTS: At follow-up, 48 of 177 HLA-B27(+) FDR had developed axSpA, an RR of 27.1% (95% CI 20.6% to 33.7%). 27/148 (18.2%) children of AS probands (modified New York (mNY) criteria) were affected versus 2/50 (4.0%) children of non-radiographic axSpA probands (p=0.0138, OR=5.36; 95% CI 1.23 to 23.40). Children of female probands were more often affected (12/22; 54.5%) than of male probands (15/78; 19.2%) (p=0.0003; OR=4.89; 95% CI 1.96 to 12.23). This increased risk applies equally to sons and daughters. CONCLUSION: The lifetime RR of axSpA for HLA-B27(+) FDR is substantial (27.1%), and disease severity (as defined by radiographic sacroiliitis by the mNY criteria) is an additional risk factor. Affected mothers pass on the disease significantly more often to their offspring than do affected fathers. These findings may lead to better assessment of lifetime risk for axSpA in the offspring. Moreover, investigation of this gender effect may uncover additional putative disease susceptibility factors.

Heterogeneity of axial spondyloarthritis: genetics, sex and structural damage matter.

OBJECTIVE: Axial spondyloarthritis (axSpA) comprises both radiographic and non-radiographic disease. However, the paucity of specific objective measures for the disease and current classification criteria showing suboptimal specificity contribute to disease heterogeneity observed in clinical practice and research. We used a historical cohort of patients with axSpA to assess sources of heterogeneity. METHODS: The study involved 363 axSpA probands recruited from membership of the Swiss Ankylosing Spondylitis Patient Society. Participants underwent examination by a rheumatologist, completed questionnaires and provided blood samples for HLA typing. Patients underwent radiography of sacroiliac joints and were categorised according to the New York (NY) criteria (ankylosing spondylitis (AS) or non-radiographic axSpA (nr-axSpA)) and HLA-B27 status. Genetic characterisation by single nucleotide polymorphism microarray was performed and AS polygenic risk scores (PRS) were calculated. RESULTS: Considerable heterogeneity was observed. The male to female ratio for AS (NY+) was 3:1, but 1:1 for nr-axSpA. For HLA-27(+) AS, the ratio was 2.5:1, but nearly 1:1 for HLA-B27(-) disease. Women with nr-axSpA had strikingly lower mean PRS and lower HLA-B27 prevalence than men with nr-axSpA or NY(+) male and female patients with AS. PRS was able to distinguish male but not female patients with nr-axSpA from related healthy first-degree relatives. Radiographic sacroiliitis was strongly associated with HLA-B27, especially in men. CONCLUSION: Women clinically diagnosed with axSpA but without radiographic sacroiliitis as a group have a disease that is distinct from AS by the modified New York criteria overall and from nr-axSpA in men. Given the high degree of heterogeneity, stratified or adjusted analysis of effectiveness studies is indicated, taking genetics, sex and radiographic damage (sacroiliitis) into account.

Impact of HLA-B27 and Disease Status on the Gut Microbiome of the Offspring of Ankylosing Spondylitis Patients.

Multiple studies have shown the microbiota to be abnormal in patients with spondyloarthritis (SpA). The purpose of this study was to explore the genetic contributions of these microbiota abnormalities. We analyzed the impact of HLA-B27 on the microbiota of children at risk for SpA and compared the microbiota of HLA-B27+ pediatric offspring of ankylosing spondylitis (AS) patients with that of HLA-B27+ children with SpA. Human DNA was obtained from the offspring for determination of HLA-B27 status and polygenic risk score (PRS). Fecal specimens were collected from both groups for sequencing of the V4 region of the 16S ribosomal RNA gene. Among the offspring of AS patients, there was slight clustering by HLA-B27 status. After adjusting for multiple comparisons, five operational taxonomic units (OTUs) representing three unique taxa distinguished the HLA-B27+ from negative children: Blautia and Coprococcus were lower in the HLA-B27+ offspring, while Faecalibacterium prausnitzii was higher. HLA-B27+ offspring without arthritis were compared to children with treatment-naïve HLA-B27+ SpA. After adjustments, clustering by diagnosis was present. A total of 21 OTUs were significantly associated with diagnosis state, including Bacteroides (higher in SpA patients) and F. prausnitzii (higher in controls). Thus, our data confirmed associations with B. fragilis and F. prausnitzii with juvenile SpA, and also suggest that the mechanism by which HLA-B27 is associated with SpA may not involve alterations of the microbiota.

Novel approach by natural language processing for COVID-19 knowledge discovery.

BACKGROUND: The impact of COVID-19 on public health has mandated an 'all hands on deck' scientific response. The current clinical study and basic research on COVID-19 are mainly based on existing publications or our knowledge of coronavirus. However, efficiently retrieval of accurate, relevant knowledge on COVID-19 can pose significant challenges for researchers. METHODS: To improve quality in accessing important literature findings, we developed a novel natural language processing (NLP) method to automatically recognize the associations among potential targeted host organ systems, associated clinical manifestations, and pathways. We further validated these associations through clinician experts' evaluations and prioritize candidate drug targets through bioinformatics network analysis. RESULTS: We found that the angiotensin-converting enzyme 2 (ACE2), a receptor that SARS-CoV-2 required for cell entry, is associated with cardiovascular and endocrine organ system and diseases. Furthermore, we found SARS-CoV-2 is associated with some important pathways such as IL-6, TNF-alpha, and IL-1 beta-induced dyslipidemia, which are related to inflammation, lipogenesis, and oxidative stress mechanisms, suggesting potential drug candidates. CONCLUSION: We prioritized the list of therapeutic targets involved in antiviral and immune modulating drugs for experimental validation, rendering it valuable during public health crises marked by stresses on clinical and research capacity. Our automatic intelligence pipeline also contributes to other novel and emerging disease management and treatments in the future.

Factors predicting axial spondyloarthritis among first-degree relatives of probands with ankylosing spondylitis: a family study spanning 35 years.

OBJECTIVE: Factors predicting axial spondyloarthritis (axSpA) among first-degree relatives (FDRs) of ankylosing spondylitis (AS) patients need to be defined. We investigated the predictive value of the probands' HLA-B27 and radiographic sacroiliitis status on disease occurrence among their FDR. We also assessed the predictive value of features of the clinical history, including chronic inflammatory back pain (CIBP) and acute anterior uveitis (AAU), among the FDR and how they can be used to improve classification and diagnosis of axSpA. METHODS: In 1985, we studied 363 AS probands and 806 FDR who underwent rheumatologic examination, completed questionnaires, provided blood samples for HLA-typing and underwent radiography of sacroiliac joints. At follow-up in 2018-2019, 125 patients and 360 FDR were available for study, and completed a postal questionnaire about axSpA features. FDRs were asked to report whether after 1985 they had been diagnosed by Swiss rheumatologists as having axSpA. RESULTS: Among HLA-B27(+) FDR, axSpA occurred in 25.4%-26.3%, independent of the radiographic sacroiliitis status of the proband. AAU occurred in 13/34 (38.2%) FDR with axSpA vs 29/251 (11.6%) FDR without axSpA (p=0.00004, OR=4.74 95% CI 2.15 to 10.47). The presence of CIBP at baseline did not predict later occurrence of axSpA but combining CIBP and pain/discomfort at the thoracic spine and at anterior (ventral) chest wall ever, assessed at follow-up in 2018-2019, provided 83.1% sensitivity and 87.2% specificity for current axSpA. CONCLUSION: Occurrence of AAU among FDR of axSpA probands should prompt screening for axSpA. Moreover, co-occurrence of CIBP and pain/discomfort in the thoracic spine and at anterior chest wall as a three-question tool may further enhance clinical suspicion of axSpA among these FDR.

Oligofructose-Enriched Inulin Intake, Gut Microbiome Characteristics, and the V̇O2 Peak Response to High-Intensity Interval Training in Healthy Inactive Adults.

BACKGROUND: The gut microbiome has been associated with cardiorespiratory fitness. OBJECTIVES: To assess the effects of oligofructose (FOS)-enriched inulin supplementation on the gut microbiome and the peak oxygen uptake (V̇O2peak) response to high-intensity interval training (HIIT). METHODS: The study was a randomized controlled trial. Forty sedentary and apparently healthy adults [n = 31 women; aged 31.8 ± 9.8 y, BMI (in kg⋅m-2) 25.9 ± 4.3] were randomly allocated to 1) 6 wk of supervised HIIT (4 × 4-min bouts at 85-95% peak heart rate, interspersed with 3 min of active recovery, 3·wk-1) + 12 g·d-1 of FOS-enriched inulin (HIIT-I) or 2) 6 wk of supervised HIIT (3·wk-1, 4 × 4-min bouts) + 12 g·d-1 of maltodextrin/placebo (HIIT-P). Each participant completed an incremental treadmill test to assess V̇O2peak and ventilatory thresholds (VTs), provided a stool and blood sample, and completed a 24-h diet recall questionnaire and FFQ before and after the intervention. Gut microbiome analyses were performed using metagenomic sequencing. Fecal short-chain fatty acids were measured by mass spectrometry. RESULTS: There were no differences in the mean change in V̇O2peak response between groups (P = 0.58). HIIT-I had a greater improvement in VTs than HIIT-P [VT1 (lactate accumulation): mean difference + 4.3% and VT2 (lactate threshold): +4.2%, P < 0.05]. HIIT-I had a greater increase in the abundance of Bifidobacterium taxa [false discovery rate (FDR) < 0.05] and several metabolic processes related to exercise capacity (FDR < 0.05). Exploratory analysis of merged data found participants with a greater response to HIIT (V̇O2peak ≥3.5 mL⋅kg-1⋅min-1) had a 2.2-fold greater mean abundance of gellan degradation pathways (FDR < 0.05) and a greater, but not significant, abundance of Bifidobacterium uniformis species (P < 0.00023, FDR = 0.08). CONCLUSIONS: FOS-enriched inulin supplementation did not potentiate HIIT-induced improvements in V̇O2peak but led to gut microbiome changes possibly associated with greater ventilatory threshold improvements in healthy inactive adults. Gellan degradation pathways and B. uniformis spp. were associated with greater V̇O2peak responses to HIIT.

Intratumoral T-cell receptor repertoire is predictive of interim PET scan results in patients with diffuse large B-cell lymphoma treated with rituximab/cyclophosphamide/doxorubicin/prednisolone/vincristine (R-CHOP) chemoimmunotherapy.

OBJECTIVES: A diverse intratumoral T-cell receptor (TCR) repertoire is associated with improved survival in diffuse large B-cell lymphoma (DLBCL) treated with rituximab/cyclophosphamide/doxorubicin/prednisolone/vincristine (R-CHOP) chemoimmunotherapy. We explored the impact of intratumoral TCR repertoire on interim PET (iPET) done after four cycles of R-CHOP, the relationships between intratumoral and circulating repertoire, and the phenotypes of expanded clonotypes. METHODS: We sequenced the third complementarity-determining region of TCRß in tumor samples, blood at pre-therapy and after four cycles of R-CHOP in 35 patients enrolled in ALLGNHL21 trial in high-risk DLBCL. We correlated the TCR diversity metrics with iPET status, gene expression profiles and HLA-class I genotypes. We then sequenced the FACS-sorted peripheral blood T cells in six patients, and pentamer-sorted EBV-specific CD8+ T cells in one patient from this cohort. RESULTS: Compared with iPET- patients, the intratumoral TCR repertoire in iPET+ patients was characterised by higher cumulative frequency of abundant clonotypes and higher productive clonality. There was a variable overlap between circulating and intratumoral repertoires, with the dominant intratumoral clonotypes more likely to be detected in the blood. The majority of shared clonotypes were CD8+ PD-1HI T cells, and CD8+ T cells had the largest clonal expansions in tumor and blood. In a patient with EBV+ DLBCL, EBV-specific intratumoral clonotypes were trackable in the blood. CONCLUSION: This study demonstrates that clonally expanded intratumoral TCR repertoires are associated with iPET+ and that the blood can be used to track tumor-associated antigen-specific clonotypes. These findings assist the rationale design and therapeutic monitoring of immunotherapeutic strategies in DLBCL.

De novo protein design by an energy function based on series expansion in distance and orientation dependence.

MOTIVATION: Despite many successes, de novo protein design is not yet a solved problem as its success rate remains low. The low success rate is largely because we do not yet have an accurate energy function for describing the solvent-mediated interaction between amino acid residues in a protein chain. Previous studies showed that an energy function based on series expansions with its parameters optimized for side-chain and loop conformations can lead to one of the most accurate methods for side chain (OSCAR) and loop prediction (LEAP). Following the same strategy, we developed an energy function based on series expansions with the parameters optimized in four separate stages (recovering single-residue types without and with orientation dependence, selecting loop decoys and maintaining the composition of amino acids). We tested the energy function for de novo design by using Monte Carlo simulated annealing. RESULTS: The method for protein design (OSCAR-Design) is found to be as accurate as OSCAR and LEAP for side-chain and loop prediction, respectively. In de novo design, it can recover native residue types ranging from 38% to 43% depending on test sets, conserve hydrophobic/hydrophilic residues at ∼75%, and yield the overall similarity in amino acid compositions at more than 90%. These performance measures are all statistically significantly better than several protein design programs compared. Moreover, the largest hydrophobic patch areas in designed proteins are near or smaller than those in native proteins. Thus, an energy function based on series expansion can be made useful for protein design. AVAILABILITY AND IMPLEMENTATION: The Linux executable version is freely available for academic users at http://zhouyq-lab.szbl.ac.cn/resources/.

Polygenic Risk Scores have high diagnostic capacity in ankylosing spondylitis.

OBJECTIVE: We sought to test the hypothesis that Polygenic Risk Scores (PRSs) have strong capacity to discriminate cases of ankylosing spondylitis (AS) from healthy controls and individuals in the community with chronic back pain. METHODS: PRSs were developed and validated in individuals of European and East Asian ethnicity, using data from genome-wide association studies in 15 585 AS cases and 20 452 controls. The discriminatory values of PRSs in these populations were compared with other widely used diagnostic tests, including C-reactive protein (CRP), HLA-B27 and sacroiliac MRI. RESULTS: In people of European descent, PRS had high discriminatory capacity with area under the curve (AUC) in receiver operator characteristic analysis of 0.924. This was significantly better than for HLA-B27 testing alone (AUC=0.869), MRI (AUC=0.885) or C-reactive protein (AUC=0.700). PRS developed and validated in individuals of East Asian descent performed similarly (AUC=0.948). Assuming a prior probability of AS of 10% such as in patients with chronic back pain under 45 years of age, compared with HLA-B27 testing alone, PRS provides higher positive values for 35% of patients and negative predictive values for 67.5% of patients. For PRS, in people of European descent, the maximum positive predictive value was 78.2% and negative predictive value was 100%, whereas for HLA-B27, these values were 51.9% and 97.9%, respectively. CONCLUSIONS: PRS have higher discriminatory capacity for AS than CRP, sacroiliac MRI or HLA-B27 status alone. For optimal performance, PRS should be developed for use in the specific ethnic groups to which they are to be applied.

Genome wide association study of response to interval and continuous exercise training: the Predict-HIIT study.

BACKGROUND: Low cardiorespiratory fitness (V̇O2peak) is highly associated with chronic disease and mortality from all causes. Whilst exercise training is recommended in health guidelines to improve V̇O2peak, there is considerable inter-individual variability in the V̇O2peak response to the same dose of exercise. Understanding how genetic factors contribute to V̇O2peak training response may improve personalisation of exercise programs. The aim of this study was to identify genetic variants that are associated with the magnitude of V̇O2peak response following exercise training. METHODS: Participant change in objectively measured V̇O2peak from 18 different interventions was obtained from a multi-centre study (Predict-HIIT). A genome-wide association study was completed (n = 507), and a polygenic predictor score (PPS) was developed using alleles from single nucleotide polymorphisms (SNPs) significantly associated (P < 1 × 10-5) with the magnitude of V̇O2peak response. Findings were tested in an independent validation study (n = 39) and compared to previous research. RESULTS: No variants at the genome-wide significance level were found after adjusting for key covariates (baseline V̇O2peak, individual study, principal components which were significantly associated with the trait). A Quantile-Quantile plot indicates there was minor inflation in the study. Twelve novel loci showed a trend of association with V̇O2peak response that reached suggestive significance (P < 1 × 10-5). The strongest association was found near the membrane associated guanylate kinase, WW and PDZ domain containing 2 (MAGI2) gene (rs6959961, P = 2.61 × 10-7). A PPS created from the 12 lead SNPs was unable to predict V̇O2peak response in a tenfold cross validation, or in an independent (n = 39) validation study (P > 0.1). Significant correlations were found for beta coefficients of variants in the Predict-HIIT (P < 1 × 10-4) and the validation study (P < × 10-6), indicating that general effects of the loci exist, and that with a higher statistical power, more significant genetic associations may become apparent. CONCLUSIONS: Ongoing research and validation of current and previous findings is needed to determine if genetics does play a large role in V̇O2peak response variance, and whether genomic predictors for V̇O2peak response trainability can inform evidence-based clinical practice. Trial registration Australian New Zealand Clinical Trials Registry (ANZCTR), Trial Id: ACTRN12618000501246, Date Registered: 06/04/2018, http://www.anzctr.org.au/Trial/Registration/TrialReview.aspx?id=374601&isReview=true .

Genetics and the axial spondyloarthritis spectrum.

The axial SpAs (axSpAs) are clearly clinically a heterogeneous set of diseases with markedly varying extra-articular features. These diseases are all highly heritable and have overlapping but differing genetic origins. Shared features include association with HLA class I alleles and genes of the IL-23 pathway, among other things. Significant differences do exist however, both in the genetic loci involved and at specific loci in the individual genetic variants associated with each disease. These similarities and differences are of great interest in regards to disease pathogenesis and treatment development, although individually they are too small in effect to be of prognostic or diagnostic value. Polygenic risk scores, which capture a high proportion of the genetic variation between disorders, have been shown to have clinically useful discriminatory capacity in axSpA. This suggests they have the potential to enable improved disease classification, incorporating basic pathogenic features such as genomics, and ultimately benefitting clinical care. The aim of this article is to review the genetic characteristics of the spectrum of axSpAs and to discuss how this influences our understanding of the disease pathogenesis and the clinical implications of this understanding.

Biomarker development for axial spondyloarthritis.

The term axial spondyloarthritis (axSpA) encompasses a heterogeneous group of diseases that have variable presentations, extra-articular manifestations and clinical outcomes, and that will respond differently to treatments. The prototypical type of axSpA, ankylosing spondylitis, is thought to be caused by interaction between the genetically primed host immune system and gut microbiota. Currently used biomarkers such as HLA-B27 status, C-reactive protein and erythrocyte sedimentation rate have, at best, moderate diagnostic and predictive value. Improved biomarkers are needed for axSpA to assist with early diagnosis and to better predict treatment responses and long-term outcomes. Advances in a range of 'omics' technologies and statistical approaches, including genomics approaches (such as polygenic risk scores), microbiome profiling and, potentially, transcriptomic, proteomic and metabolomic profiling, are making it possible for more informative biomarker sets to be developed for use in such clinical applications. Future developments in this field will probably involve combinations of biomarkers that require novel statistical approaches to analyse and to produce easy to interpret metrics for clinical application. Large publicly available datasets from well-characterized case-cohort studies that use extensive biological sampling, particularly focusing on early disease and responses to medications, are required to establish successful biomarker discovery and validation programmes.

Identification of susceptibility variants to benign childhood epilepsy with centro-temporal spikes (BECTS) in Chinese Han population.

BACKGROUND: Benign Childhood Epilepsy with Centro-temporal Spikes (BECTS) is the most common form of idiopathic epilepsy in children, accounting for up to 23% of pediatric epilepsy. The pathogenesis of BECTS is unknown, but it is thought that genetic factors play a role in susceptibility to the disease. METHODS: To investigate the role of common genetic variants in BECTS pathogenesis, a 2-stage genome-wide association study (GWAS) was performed in 1,800 Chinese Han BECTS patients, and 7,090 healthy controls. Genetic findings were used in a Mendelian Randomization study in the UK Biobank dataset to investigate the potential role of smoking in BECTS. FINDINGS: Definitive evidence of a role for common-variant heritability was demonstrated, with heritability of BECTS of >10% observed even with conservative disease prevalence assumptions. Although no individual locus achieved genome-wide significance, twelve loci achieved suggestive evidence of association (5 × 10-8

A missense mutation in the MLKL brace region promotes lethal neonatal inflammation and hematopoietic dysfunction.

MLKL is the essential effector of necroptosis, a form of programmed lytic cell death. We have isolated a mouse strain with a single missense mutation, MlklD139V, that alters the two-helix 'brace' that connects the killer four-helix bundle and regulatory pseudokinase domains. This confers constitutive, RIPK3 independent killing activity to MLKL. Homozygous mutant mice develop lethal postnatal inflammation of the salivary glands and mediastinum. The normal embryonic development of MlklD139V homozygotes until birth, and the absence of any overt phenotype in heterozygotes provides important in vivo precedent for the capacity of cells to clear activated MLKL. These observations offer an important insight into the potential disease-modulating roles of three common human MLKL polymorphisms that encode amino acid substitutions within or adjacent to the brace region. Compound heterozygosity of these variants is found at up to 12-fold the expected frequency in patients that suffer from a pediatric autoinflammatory disease, chronic recurrent multifocal osteomyelitis (CRMO).

Genomewide Association Study of Acute Anterior Uveitis Identifies New Susceptibility Loci.

Purpose: Acute anterior uveitis (AAU) is a common intraocular inflammatory disease. AAU occurs in 30% to 50% of patients with ankylosing spondylitis (AS), and both conditions are strongly associated with human leukocyte antigen (HLA)-B27, implying a shared etiology. This study aims to apply genomewide association study (GWAS) to characterize the genetic associations of AAU and their relationship to the genetics of AS. Methods: We undertook the GWAS analyses in 2752 patients with AS with AAU (cases) and 3836 patients with AS without AAU (controls). There were 7,436,415 single-nucleotide polymorphisms (SNPs) available after SNP microarray genotyping, imputation, and quality-control filtering. Results: We identified one locus associated with AAU at genomewide significance: rs9378248 (P = 2.69 × 10-8, odds ratio [OR] = 0.78), lying close to HLA-B. Suggestive association was observed at 11 additional loci, including previously reported AS loci ERAP1 (rs27529, P = 2.19 × 10-7, OR = 1.22) and NOS2 (rs2274894, P = 8.22 × 10-7, OR = 0.83). Multiple novel suggestive associations were also identified, including MERTK (rs10171979, P = 2.56 × 10-6, OR = 1.20), KIFAP3 (rs508063, P = 5.64 × 10-7, OR = 1.20), CLCN7 (rs67412457, P = 1.33 × 10-6, OR = 1.25), ACAA2 (rs9947182, P = 9.70 × 10-7, OR = 1.37), and 5 intergenic loci. The SNP-based heritability is approximately 0.5 for AS alone, and is much higher (approximately 0.7) for AS with AAU. Consistent with the high heritability, a genomewide polygenic risk score shows strong power in identifying individuals at high risk of either AS with AAU or AS alone. Conclusions: We report here the first GWAS for AAU and identify new susceptibility loci. Our findings confirm the strong overlap in etiopathogenesis of AAU with AS, and also provide new insights into the genetic basis of AAU.