Multivariate analysis of a missense variant in CREBRF reveals associations with measures of adiposity in people of Polynesian ancestries.

The minor allele of rs373863828, a missense variant in CREB3 Regulatory Factor, is associated with several cardiometabolic phenotypes in Polynesian peoples. To better understand the variant, we tested the association of rs373863828 with a panel of correlated phenotypes (body mass index [BMI], weight, height, HDL cholesterol, triglycerides, and total cholesterol) using multivariate Bayesian association and network analyses in a Samoa cohort (n = 1632), Aotearoa New Zealand cohort (n = 1419), and combined cohort (n = 2976). An expanded set of phenotypes (adding estimated fat and fat-free mass, abdominal circumference, hip circumference, and abdominal-hip ratio) was tested in the Samoa cohort (n = 1496). In the Samoa cohort, we observed significant associations (log10 Bayes Factor [BF] ≥ 5.0) between rs373863828 and the overall phenotype panel (8.81), weight (8.30), and BMI (6.42). In the Aotearoa New Zealand cohort, we observed suggestive associations (1.5 < log10 BF < 5) between rs373863828 and the overall phenotype panel (4.60), weight (3.27), and BMI (1.80). In the combined cohort, we observed concordant signals with larger log10 BFs. In the Samoa-specific expanded phenotype analyses, we also observed significant associations between rs373863828 and fat mass (5.65), abdominal circumference (5.34), and hip circumference (5.09). Bayesian networks provided evidence for a direct association of rs373863828 with weight and indirect associations with height and BMI.

A Polynesian-specific copy number variant encompassing the MICA gene associates with gout.

Gout is of particularly high prevalence in the Maori and Pacific (Polynesian) populations of Aotearoa New Zealand (NZ). Here, we investigated the contribution of common population-specific copy number variation (CNV) to gout in the Aotearoa NZ Polynesian population. Microarray-generated genome-wide genotype data from Aotearoa NZ Polynesian individuals with (n = 1196) and without (n = 1249) gout were analyzed. Comparator population groups were 552 individuals of European ancestry and 1962 of Han Chinese ancestry. Levels of circulating major histocompatibility complex (MHC) class I polypeptide-related sequence A (MICA) were measured by enzyme-linked immunosorbent assay. Fifty-four CNV regions (CNVRs) appearing in at least 10 individuals were detected, of which seven common (>2%) CNVRs were specific to or amplified in Polynesian people. A burden test of these seven revealed associations of insertion/deletion with gout (odds ratio (OR) 95% confidence interval [CI] = 1.80 [1.01; 3.22], P = 0.046). Individually testing of the seven CNVRs for association with gout revealed nominal association of CNVR1 with gout in Western Polynesian (Chr6: 31.36-31.45 Mb, OR = 1.72 [1.03; 2.92], P = 0.04), CNVR6 in the meta-analyzed Polynesian sample sets (Chr1: 196.75-196.92 Mb, OR = 1.86 [1.16; 3.00], P = 0.01) and CNVR9 in Western Polynesian (Chr1: 189.35-189.54 Mb, OR = 2.75 [1.15; 7.13], P = 0.03). Analysis of European gout genetic association data demonstrated a signal of association at the CNVR1 locus that was an expression quantitative trait locus for MICA. The most common CNVR (CNVR1) includes deletion of the MICA gene, encoding an immunomodulatory protein. Expression of MICA was reduced in the serum of individuals with the deletion. In summary, we provide evidence for the association of CNVR1 containing MICA with gout in Polynesian people, implicating class I MHC-mediated antigen presentation in gout.

Association of rare and common genetic variants in MOCOS with inadequate response to allopurinol.

OBJECTIVES: The minor allele of the common rs2231142 ABCG2 variant predicts inadequate response to allopurinol urate lowering therapy. We hypothesize that additional variants in genes encoding urate transporters and allopurinol-to-oxypurinol metabolic enzymes also predict allopurinol response. METHODS: This study included a subset of participants with gout from the Long-term Allopurinol Safety Study Evaluating Outcomes in Gout Patients, whose whole genome was sequenced (n = 563). Good responders had a 4:1 or 5:1 ratio of good (serum urate (SU) <0.36 mmol/l on allopurinol ≤300 mg/day) to poor (SU ≥ 0.36 mmol/l despite allopurinol >300 mg/day) responses over 5-6 timepoints, while inadequate responders had a 1:4 or 1:5 ratio of good to poor responses. Adherence to allopurinol was determined by pill counts, and for a subgroup (n = 303), by plasma oxypurinol >20µmol/l. Using the sequence kernel association test (SKAT) we estimated the combined effect of rare and common variants in urate secretory (ABCC4, ABCC5, ABCG2, SLC17A1, SLC17A3, SLC22A6, SLC22A8) and reuptake genes (SLC2A9, SLC22A11) and in allopurinol-to-oxypurinol metabolic genes (AOX1, MOCOS, XDH) on allopurinol response. RESULTS: There was an association of rare and common variants in the allopurinol-to-oxypurinol gene group (PSKAT-C = 0.019), and in MOCOS, encoding molybdenum cofactor sulphurase, with allopurinol response (PSKAT-C = 0.011). Evidence for genetic association with allopurinol response in the allopurinol-to-oxypurinol gene group (PSKAT-C = 0.002) and MOCOS (PSKAT-C < 0.001) was stronger when adherence to allopurinol therapy was confirmed by plasma oxypurinol. CONCLUSION: We provide evidence for common and rare genetic variation in MOCOS associating with allopurinol response.

Development and Internal Validation of a Clinical and Genetic Risk Score for Rheumatoid Arthritis-Associated Interstitial Lung Disease.

OBJECTIVE: Although clinical and genetic risk factors have been identified for rheumatoid arthritis-associated interstitial lung disease (RA-ILD), there are no current tools allowing for risk stratification. We sought to develop and validate an ILD risk model in a large, multicentre, prospective RA cohort. METHODS: Participants in the Veterans Affairs RA (VARA) registry were genotyped for 12 single nucleotide polymorphisms (SNPs) associated with idiopathic pulmonary fibrosis. ILD was validated through systematic record review. A genetic risk score (GRS) was computed from minor alleles weighted by effect size with ILD, using backward selection. The GRS was combined with clinical risk factors within a logistic regression model. Internal validation was completed using bootstrapping, and model performance was assessed by the area under the receiver operating curve (AUC). RESULTS: Of 2,386 participants (89% male, mean age 69.5 years), 9.4% had ILD. Following backward selection, five SNPs contributed to the GRS. The GRS and clinical factors outperformed clinical factors alone in discriminating ILD (AUC 0.675 vs 0.635, p< 0.001). The shrinkage-corrected performance for combined and clinical-only models was 0.667 (95% CI 0.628, 0.712) and 0.623 (95% CI 0.584, 0.651), respectively. Twenty percent of the cohort had a combined risk score below a cut-point with >90% sensitivity. CONCLUSION: A clinical and genetic risk model discriminated ILD in a large, multicentre RA cohort better than a clinical-only model, excluding 20% of the cohort from low-yield testing. These results demonstrate the potential utility of a GRS in RA-ILD and support further investigation into individualized risk stratification and screening.

The development and evaluation of dose-prediction tools for allopurinol therapy (Easy-Allo tools).

AIMS: Dose escalation at the initiation of allopurinol therapy can be protracted and resource intensive. Tools to predict the allopurinol doses required to achieve target serum urate concentrations would facilitate the implementation of more efficient dose-escalation strategies. The aim of this research was to develop and externally evaluate allopurinol dosing tools, one for use when the pre-urate-lowering therapy serum urate is known (Easy-Allo1) and one for when it is not known (Easy-Allo2). METHODS: A revised population pharmacokinetic-pharmacodynamic model was developed using data from 653 people with gout. Maintenance doses to achieve the serum urate target of <0.36 mmol L-1 in >80% of individuals were simulated and evaluated against external data. The predicted and observed allopurinol doses were compared using the mean prediction error (MPE) and root mean square error (RMSE). The proportion of Easy-Allo predicted doses within 100 mg of the observed was quantified. RESULTS: Allopurinol doses were predicted by total body weight, baseline urate, ethnicity and creatinine clearance. Easy-Allo1 produced unbiased and suitably precise dose predictions (MPE 2 mg day-1 95% confidence interval [CI] -13-17, RMSE 91%, 90% within 100 mg of the observed dose). Easy-Allo2 was positively biased by about 70 mg day-1 and slightly less precise (MPE 70 mg day-1 95% CI 52-88, RMSE 131%, 71% within 100 mg of the observed dose). CONCLUSIONS: The Easy-Allo tools provide a guide to the allopurinol maintenance dose requirement to achieve the serum urate target of <0.36 mmol L-1 and will aid in the development of novel dose-escalation strategies for allopurinol therapy.

Genetic and Physiological Effects of Insulin-Like Growth Factor-1 (IGF-1) on Human Urate Homeostasis.

SIGNIFICANCE STATEMENT: Hyperinsulinemia induces hyperuricemia by activating net renal urate reabsorption in the renal proximal tubule. The basolateral reabsorptive urate transporter GLUT9a appears to be the dominant target for insulin. By contrast, IGF-1 infusion reduces serum urate (SU), through mechanisms unknown. Genetic variants of IGF1R associated with reduced SU have increased IGF-1R expression and interact with genes encoding the GLUT9 and ABCG2 urate transporters, in a sex-specific fashion, which controls the SU level. Activation of IGF-1/IGF-1R signaling in Xenopus oocytes modestly activates GLUT9a and inhibits insulin's stimulatory effect on the transporter, which also activates multiple secretory urate transporters-ABCG2, ABCC4, OAT1, and OAT3. The results collectively suggest that IGF-1 reduces SU by activating secretory urate transporters and inhibiting insulin's action on GLUT9a. BACKGROUND: Metabolic syndrome and hyperinsulinemia are associated with hyperuricemia. Insulin infusion in healthy volunteers elevates serum urate (SU) by activating net urate reabsorption in the renal proximal tubule, whereas IGF-1 infusion reduces SU by mechanisms unknown. Variation within the IGF1R gene also affects SU levels. METHODS: Colocalization analyses of a SU genome-wide association studies signal at IGF1R and expression quantitative trait loci signals in cis using COLOC2, RT-PCR, Western blotting, and urate transport assays in transfected HEK 293T cells and in Xenopus laevis oocytes. RESULTS: Genetic association at IGF1R with SU is stronger in women and is mediated by control of IGF1R expression. Inheritance of the urate-lowering homozygous genotype at the SLC2A9 locus is associated with a differential effect of IGF1R genotype between men and women. IGF-1, through IGF-1R, stimulated urate uptake in human renal proximal tubule epithelial cells and transfected HEK 293T cells, through activation of IRS1, PI3/Akt, MEK/ERK, and p38 MAPK; urate uptake was inhibited in the presence of uricosuric drugs, specific inhibitors of protein tyrosine kinase, PI3 kinase (PI3K), ERK, and p38 MAPK. In X. laevis oocytes expressing ten individual urate transporters, IGF-1 through endogenous IGF-1R stimulated urate transport mediated by GLUT9, OAT1, OAT3, ABCG2, and ABCC4 and inhibited insulin's stimulatory action on GLUT9a and OAT3. IGF-1 significantly activated Akt and ERK. Specific inhibitors of PI3K, ERK, and PKC significantly affected IGF-1 stimulation of urate transport in oocytes. CONCLUSIONS: The combined results of infusion, genetics, and transport experiments suggest that IGF-1 reduces SU by activating urate secretory transporters and inhibiting insulin's action.

Genetic Risk Scores for the Clinical Rheumatologist.

BACKGROUND/HISTORICAL PERSPECTIVE: The advent of genome-wide sequencing and large-scale genetic epidemiological studies has led to numerous opportunities for the application of genetics in clinical medicine. Leveraging this information toward the formation of clinically useful tools has been an ongoing research goal in this area. A genetic risk score (GRS) is a measure that attempts to estimate the cumulative contribution of established genetic risk factors toward an outcome of interest, taking into account the cumulative risk that each of these individual genetic risk factors conveys. The purpose of this perspective is to provide a systematic framework to evaluate a GRS for clinical application. SUMMARY OF CURRENT LITERATURE: Since the initial polygenic risk score methodology in 2007, there has been increasing GRS application across the medical literature. In rheumatology, this has included application to rheumatoid arthritis, gout, spondyloarthritis, lupus, and inflammatory arthritis. MAJOR CONCLUSIONS: GRSs are particularly relevant to rheumatology, where common diseases have many complex genetic factors contributing to risk. Despite this, there is no widely accepted method for the critical application of a GRS, which can be a particular challenge for the clinical rheumatologist seeking to clinically apply GRSs. This review provides a framework by which the clinician may systematically evaluate a GRS. FUTURE RESEARCH DIRECTIONS: As genotyping becomes more accessible and cost-effective, it will become increasingly important to recognize the clinical applicability of GRSs and identify those of the highest utility for patient care. This framework for the evaluation of a GRS will also help ensure reliability among GRS research in rheumatology, thereby helping to advance the field.

Association of rs9939609 in FTO with BMI among Polynesian peoples living in Aotearoa New Zealand and other Pacific nations.

The fat mass and obesity associated (FTO) locus consistently associates with higher body mass index (BMI) across diverse ancestral groups. However, previous small studies of people of Polynesian ancestries have failed to replicate the association. In this study, we used Bayesian meta-analysis to test rs9939609, the most replicated FTO variant, for association with BMI with a large sample (n = 6095) of Aotearoa New Zealanders of Polynesian (Maori and Pacific) ancestry and of Samoan people living in the Independent State of Samoa and in American Samoa. We did not observe statistically significant association within each separate Polynesian subgroup. Bayesian meta-analysis of the Aotearoa New Zealand Polynesian and Samoan samples resulted in a posterior mean effect size estimate of +0.21 kg/m2, with a 95% credible interval [+0.03 kg/m2, +0.39 kg/m2]. While the Bayes Factor (BF) of 0.77 weakly favors the null, the BF = 1.4 Bayesian support interval is [+0.04, +0.20]. These results suggest that rs9939609 in FTO may have a similar effect on mean BMI in people of Polynesian ancestries as previously observed in other ancestral groups.

A functional genomics approach reveals suggestive quantitative trait loci associated with combined TLR4 and BCP crystal-induced inflammation and osteoarthritis.

OBJECTIVE: Basic calcium phosphate (BCP) crystals can activate the NLRP3 inflammasome and are potentially involved in the pathogenesis of osteoarthritis (OA). In order to elucidate relevant inflammatory mechanisms in OA, we used a functional genomics approach to assess genetic variation influencing BCP crystal-induced cytokine production. METHOD: Peripheral blood mononuclear cells (PBMCs) were isolated from healthy volunteers who were previously genotyped and stimulated with BCP crystals and/or lipopolysaccharide (LPS) after which cytokines release was assessed. Cytokine quantitative trait locus (cQTL) mapping was performed. For in vitro validation of the cQTL located in anoctamin 3 (ANO3), PBMCs were incubated with Tamoxifen and Benzbromarone prior to stimulation. Additionally, we performed co-localisation analysis of our top cQTLs with the most recent OA meta-analysis of genome-wide association studies (GWAS). RESULTS: We observed that BCP crystals and LPS synergistically induce IL-1ß in human PBMCs. cQTL analysis revealed several suggestive loci influencing cytokine release upon stimulation, among which are quantitative trait locus annotated to ANO3 and GLIS3. As functional validation, anoctamin inhibitors reduced IL-1ß release in PBMCs after stimulation. Co-localisation analysis showed that the GLIS3 locus was shared between LPS/BCP crystal-induced IL-1ß and genetic association with Knee OA. CONCLUSIONS: We identified and functionally validated a new locus, ANO3, associated with LPS/BCP crystal-induced inflammation in PBMCs. Moreover, the cQTL in the GLIS3 locus co-localises with the previously found locus associated with Knee OA, suggesting that this Knee OA locus might be explained through an inflammatory mechanism. These results form a basis for further exploration of inflammatory mechanisms in OA.

Profiling of serum oxylipins identifies distinct spectrums and potential biomarkers in young people with very early onset gout.

OBJECTIVE: Oxylipins modulate inflammation via complex pathways. The oxylipin profile in gout remains unexplored. In this study, we systemically profiled oxylipins in young men and identified new oxylipin biomarkers for clinical use in differentiating gout from hyperuricaemia. MATERIAL AND METHODS: Oxylipin profiling was performed in 90 men (30 very early onset gout, 30 asymptomatic hyperuricaemia [HU] and 30 normouricaemia [NU], all aged <20 years) divided into discovery and validation sample sets. The dataset was analysed based on orthogonal projection to latent structure-discriminant analysis. Correlation network and pathway enrichment were conducted to reveal potential oxylipin-involved pathways of gout. Candidate oxylipins were further evaluated and optimized in the validation cohort, and differential oxylipin biomarkers combined with or without serum urate were applied to construct diagnostic models. RESULTS: In discovery stage, 21 differential oxylipins in the gout vs HU comparisons and 14 differential oxylipins in the gout vs NU comparisons were discovered. Correlation network analysis was performed and 14(S)-HDHA (14S-hydroxy-4Z,7Z,10Z,12E,16Z,19Z-docosahexaenoic acid) was identified as a hub metabolite in both comparisons. Seven down-regulated oxylipins in the gout vs HU group and five down-regulated oxylipins in the gout vs NU group were validated. Diagnostic models were constructed with the above oxylipins, with 14(S)-HDHA alone having an area under the curve of 1 (95% CI, 1, 1) in both comparisons. CONCLUSIONS: Young men with very early onset gout have distinct oxylipin spectrums, especially those derived from arachidonic acid, eicosapentaenoic acid and docosahexaenoic acid. Differential oxylipins could serve as candidate serum biomarkers in differentiating gout from hyperuricaemia.

An allopurinol adherence tool using plasma oxypurinol concentrations.

AIMS: This study aimed to develop and evaluate an allopurinol adherence tool based on steady-state oxypurinol plasma concentrations, allopurinol's active metabolite. METHODS: Plasma oxypurinol concentrations were simulated stochastically from an oxypurinol pharmacokinetic model for allopurinol doses of 100-800 mg daily, accounting for differences in renal function, diuretic use and ethnicity. For each scenario, the 20th percentile for peak and trough concentrations defined the adherence threshold, below which imperfect adherence was assumed. Predictive performance was evaluated using both simulated low adherence and against data from 146 individuals with paired oxypurinol plasma concentrations and adherence measures. Sensitivity and specificity (S&S), negative and positive predictive values (NPV, PPV) and receiver operating characteristic (ROC) area under the curve (AUC) were determined. The predictive performance of the tool was evaluated using adherence data from an external study (CKD-FIX). RESULTS: The allopurinol adherence tool produced S&S values for trough thresholds of 89-98% and 76-84%, respectively, and 90%-98% and 76-83% for peak thresholds. PPV and NPV were 79-84% and 88-94%, respectively, for trough and 80-85% and 89-98%, respectively, for peak concentrations. The ROC AUC values ranged from 0.84 to 0.88 and from 0.86 to 0.89 for trough and peak concentrations, respectively. S&S values for the external evaluation were found to be 75.8% and 86.5%, respectively, producing an ROC AUC of 0.8113. CONCLUSION: A tool to identify people with gout who require additional support to maintain adherence using plasma oxypurinol concentrations was developed and evaluated. The predictive performance of the tool is suitable for adherence screening in clinical trials and may have utility in some clinical practice settings.

Cadmium exposure dysregulates purine metabolism and homeostasis across the gut-liver axis in a mouse model.

Cadmium (Cd) exposure has been associated with the development of enterohepatic circulation disorders and hyperuricemia, but the possible contribution of chronic low-dose Cd exposure to disease progression is still need to be explored. A mouse model of wild-type mice (WT) and Uox-knockout mice (Uox-KO) to find out the toxic effects of chronic low-dose Cd exposure on liver purine metabolism by liquid chromatography-mass spectrometry (LC-MS) platform and associated intestinal flora. High throughput omics analysis including metabolomics and transcriptomics showed that Cd exposure can cause disruption of purine metabolism and energy metabolism. Cd changes several metabolites associated with purine metabolism (xanthine, hypoxanthine, adenosine, uridine, inosine) and related genes, which are associated with elevated urate levels. Microbiome analysis showed that Cd exposure altered the disturbance of homeostasis in the gut. Uox-KO mice were more susceptible to Cd than WT mice. Our findings extend the understanding of potential toxicological interactions between liver and gut microbiota and shed light on the progression of metabolic diseases caused by Cd exposure.

Genomic dissection of 43 serum urate-associated loci provides multiple insights into molecular mechanisms of urate control.

High serum urate is a prerequisite for gout and associated with metabolic disease. Genome-wide association studies (GWAS) have reported dozens of loci associated with serum urate control; however, there has been little progress in understanding the molecular basis of the associated loci. Here, we employed trans-ancestral meta-analysis using data from European and East Asian populations to identify 10 new loci for serum urate levels. Genome-wide colocalization with cis-expression quantitative trait loci (eQTL) identified a further five new candidate loci. By cis- and trans-eQTL colocalization analysis, we identified 34 and 20 genes, respectively, where the causal eQTL variant has a high likelihood that it is shared with the serum urate-associated locus. One new locus identified was SLC22A9 that encodes organic anion transporter 7 (OAT7). We demonstrate that OAT7 is a very weak urate-butyrate exchanger. Newly implicated genes identified in the eQTL analysis include those encoding proteins that make up the dystrophin complex, a scaffold for signaling proteins and transporters at the cell membrane; MLXIP that, with the previously identified MLXIPL, is a transcription factor that may regulate serum urate via the pentose-phosphate pathway and MRPS7 and IDH2 that encode proteins necessary for mitochondrial function. Functional fine mapping identified six loci (RREB1, INHBC, HLF, UBE2Q2, SFMBT1 and HNF4G) with colocalized eQTL containing putative causal SNPs. This systematic analysis of serum urate GWAS loci identified candidate causal genes at 24 loci and a network of previously unidentified genes likely involved in control of serum urate levels, further illuminating the molecular mechanisms of urate control.

The CREBRF diabetes-protective rs373863828-A allele is associated with enhanced early insulin release in men of Maori and Pacific ancestry.

AIMS/HYPOTHESIS: The minor A allele of rs373863828 (CREBRF p.Arg457Gln) is associated with increased BMI, but reduced risk of type 2 and gestational diabetes in Polynesian (Pacific peoples and Aotearoa New Zealand Maori) populations. This study investigates the effect of the A allele on insulin release and sensitivity in overweight/obese men without diabetes. METHODS: A mixed meal tolerance test was completed by 172 men (56 with the A allele) of Maori or Pacific ancestry, and 44 (24 with the A allele) had a frequently sampled IVGTT and hyperinsulinaemic-euglycaemic clamp. Mixed linear models with covariates age, ancestry and BMI were used to analyse the association between the A allele of rs373863828 and markers of insulin release and blood glucose regulation. RESULTS: The A allele of rs373863828 is associated with a greater increase in plasma insulin 30 min following a meal challenge without affecting the elevation in plasma glucose or incretins glucagon-like polypeptide-1 or gastric inhibitory polypeptide. Consistent with this point, following an i.v. infusion of a glucose bolus, participants with an A allele had higher early (p < 0.05 at 2 and 4 min) plasma insulin and C-peptide concentrations for a similar elevation in blood glucose as those homozygous for the major (G) allele. Despite increased plasma insulin, rs373863828 genotype was not associated with a significant difference (p > 0.05) in insulin sensitivity index or glucose disposal during hyperinsulinaemic-euglycaemic clamp. CONCLUSIONS/INTERPRETATION: rs373863828-A allele associates with increased glucose-stimulated insulin release without affecting insulin sensitivity, suggesting that CREBRF p.Arg457Gln may increase insulin release to reduce the risk of type 2 diabetes.

Mid-pass whole genome sequencing enables biomedical genetic studies of diverse populations.

BACKGROUND: Historically, geneticists have relied on genotyping arrays and imputation to study human genetic variation. However, an underrepresentation of diverse populations has resulted in arrays that poorly capture global genetic variation, and a lack of reference panels. This has contributed to deepening global health disparities. Whole genome sequencing (WGS) better captures genetic variation but remains prohibitively expensive. Thus, we explored WGS at "mid-pass" 1-7x coverage. RESULTS: Here, we developed and benchmarked methods for mid-pass sequencing. When applied to a population without an existing genomic reference panel, 4x mid-pass performed consistently well across ethnicities, with high recall (98%) and precision (97.5%). CONCLUSION: Compared to array data imputed into 1000 Genomes, mid-pass performed better across all metrics and identified novel population-specific variants with potential disease relevance. We hope our work will reduce financial barriers for geneticists from underrepresented populations to characterize their genomes prior to biomedical genetic applications.

Assessing the Relationship Between Serum Urate and Urolithiasis Using Mendelian Randomization: An Analysis of the UK Biobank.

RATIONALE & OBJECTIVE: The association between hyperuricemia and urolithiasis has been previously reported. However, this association is based on observational data, which are prone to residual confounding. The aim of this work was to use Mendelian randomization (MR) to evaluate if this relationship represents a causal effect of hyperuricemia. STUDY DESIGN: MR analysis using 2 approaches: 2-stage MR and 2-sample MR. SETTING & PARTICIPANTS: Participants aged 40-69 years from the UK Biobank Resource. EXPOSURE: Serum urate. OUTCOME: Urolithiasis. ANALYTICAL APPROACH: An observational analysis testing for an association between serum urate level and urolithiasis was performed using logistic regression. For MR analyses, serum urate-associated single-nucleotide polymorphisms, identified from genome-wide association data, were used as instrumental variables for serum urate. In the 2-stage MR analysis, a weighted genetic urate score was calculated from the instrumental variables, and a control function estimation model was fit. In the 2-sample MR analysis, multiple-instrument MR via the inverse-variance weighted method was performed. RESULTS: Individual-level data were available for 359,827 participants, of whom 6,398 (1.8%) reported urolithiasis. In the observational analysis, serum urate was positively associated with urolithiasis in an unadjusted analysis (odds ratio [OR], 1.47 [95% CI, 1.42-1.51]); however, after adjustment for relevant confounders, no association was observed (OR, 1.03 [95% CI, 0.99-1.08]). In the 2-stage MR analysis, no significant causal effect of serum urate level on urolithiasis was observed in the unadjusted (OR, 0.93 [95% CI, 0.81-1.08]) or adjusted (OR, 0.94 [95% CI, 0.80-1.09]) models. In the 2-sample MR analysis, multiple-instrument MR did not indicate a causal effect of serum urate on urolithiasis. LIMITATIONS: Stone composition and urinalysis data, including urine pH, were not available for this study. CONCLUSIONS: Our analyses do not support a causal effect of serum urate level on urolithiasis. The association between serum urate level and urolithiasis reported in observational studies is likely due to residual confounding.

Trans-ancestral dissection of urate- and gout-associated major loci SLC2A9 and ABCG2 reveals primate-specific regulatory effects.

Gout is a complex inflammatory arthritis affecting ~20% of people with an elevated serum urate level (hyperuricemia). Gout and hyperuricemia are essentially specific to humans and other higher primates, with varied prevalence across ancestral groups. SLC2A9 and ABCG2 are major loci associated with both urate and gout in multiple ancestral groups. However, fine mapping has been challenging due to extensive linkage disequilibrium underlying the associated regions. We used trans-ancestral fine mapping integrated with primate-specific genomic information to address this challenge. Trans-ancestral meta-analyses of GWAS cohorts of either European (EUR) or East Asian (EAS) ancestry resulted in single-variant resolution mappings for SLC2A9 (rs3775948 for urate and rs4697701 for gout) and ABCG2 (rs2622621 for gout). Tests of colocalization of variants in both urate and gout suggested existence of a shared candidate causal variant for SLC2A9 only in EUR and for ABCG2 only in EAS. The fine-mapped gout variant rs4697701 was within an ancient enhancer, whereas rs2622621 was within a primate-specific transposable element, both supported by functional evidence from the Roadmap Epigenomics project in human primary tissues relevant to urate and gout. Additional primate-specific elements were found near both loci and those adjacent to SLC2A9 overlapped with known statistical epistatic interactions associated with urate as well as multiple super-enhancers identified in urate-relevant tissues. We conclude that by leveraging ancestral differences trans-ancestral fine mapping has identified ancestral and functional variants for SLC2A9 or ABCG2 with primate-specific regulatory effects on urate and gout.

Effects of fenofibrate therapy on renal function in primary gout patients.

OBJECTIVE: To investigate the incidence and potential risk factors for development of fenofibrate-associated nephrotoxicity in gout patients. METHODS: A total of 983 gout patients on fenofibrate treatment who visited the dedicated Gout Clinic at the Affiliated Hospital of Qingdao University between September 2016 and June 2020 were retrospectively enrolled from the electronic records system. Fenofibrate-associated nephrotoxicity was defined as an increase in serum creatinine (SCr) ≥0.3 mg/dl within 6 months of fenofibrate initiation. The change trend of SCr and uric acid levels during the treatment period were assessed by a generalised additive mixed model (GAMM). Multivariate analysis was performed for risk factors affecting elevated SCr. RESULTS: A total of 100 (10.2%) patients experienced an increase in SCr ≥0.3 mg/dl within 6 months after fenofibrate initiation. The median change of SCr in the whole cohort was 0.11 mg/dl [interquartile range (IQR) 0.03-0.20], whereas it was 0.36 (0.33-0.45) in the fenofibrate-associated nephrotoxicity group. In a multivariable regression model, chronic kidney disease (CKD) [odds ratio (OR) 2.39 (95% CI 1.48, 3.86)] and tophus [OR 2.29 (95% CI 1.39, 3.78)] were identified to be risk predictors, independent of measured covariates, of fenofibrate-associated nephrotoxicity. During the treatment period, although SCr temporarily increased, serum urate and triglyceride concentrations decreased using the interaction analysis of GAMM. Of those with fenofibrate withdrawal records, the SCr increase in 65% of patients was reversed after an average of 49 days off the drug. CONCLUSIONS: This observational study implied that fenofibrate-associated nephrotoxicity occurs frequently in gout patients, especially in patients with tophi or CKD. The potential renal risks of fenofibrate usage in gout needs additional research.

The efficacy and safety of citrate mixture vs sodium bicarbonate on urine alkalization in Chinese primary gout patients with benzbromarone: a prospective, randomized controlled study.

OBJECTIVES: To compare the efficacy and safety of citrate mixture and sodium bicarbonate on urine alkalization in gout patients under benzbromarone treatment. METHODS: A prospective, randomized, parallel controlled trial was conducted among 200 gout patients in the dedicated gout clinic of the Affiliated Hospital of Qingdao University. The participants were randomly divided into two groups (1:1), sodium bicarbonate group (3 g/day) and citrate mixture group (7 g/day). All patients were prescribed with 25 mg/day benzbromarone at initiation and maintained at a dose of 50 mg/day. Clinical and biochemical data were collected at each follow-up time point (baseline, weeks 2, 4, 8 and 12). RESULTS: A total of 182 patients completed the 12-week urine alkalization study. The urine pH value of both groups increased significantly from the baseline to the final follow-up time point (sodium bicarbonate group, 5.50-6.00, P < 0.05; citrate mixture group, 5.53-5.93, P < 0.05). While the comparisons regarding urine pH between treatment groups showed no significant differences for each time point. The estimated glomerular filtration rate (eGFR) dropped significantly after 12 weeks' trial in the sodium bicarbonate group (P < 0.01), while it was comparable between baseline and the last follow-up (P > 0.05) in the citrate mixture group. Results of urine analysis showed that the incident rate of occult blood in the sodium bicarbonate group was higher than that in the citrate mixture group (38 vs 24%, P < 0.05), accompanied by a similar occurrence of kidney stones. After 12-week follow-up, the frequency of twice gout flare in the citrate mixture group was significantly lower than that in sodium bicarbonate group (4 vs 12%, P = 0.037). No treatment-emergent adverse events occurred. CONCLUSION: The efficacy of citrate mixture on urine alkalization is comparable to sodium bicarbonate under benzbromarone treatment without significant adverse events. Citrate mixture is superior to sodium bicarbonate in lowering the incidence of urine occult blood and the frequency of gout attacks. TRIAL REGISTRATION: Registered with ChiCTR (http://www.chictr.org.cn), No. ChiCTR1800018518.

Common variants of EDA are associated with non-syndromic hypodontia.

OBJECTIVE: The aim of this case-control study was to investigate the association between non-syndromic hypodontia and nineteen common variants of candidate genes ectodysplasin A (EDA), paired box 9 (PAX9), msh homeobox 1 (MSX1) and axis inhibition protein 2 (AXIN2). SETTINGS AND SAMPLE POPULATION: Sixty-one hypodontia cases were frequency-matched to 253 controls with no missing teeth (excluding the third molars). MATERIAL AND METHODS: Self-report data and DNA samples were collected from each participant. RESULTS: The sample had a mean age of 16.6 years (SD = 7.3), with most participants being female (59.6%), and of New Zealand European origin (75.4%). Using multiple logistic regression analysis, it was found that the T-allele of rs12853659 (EDA) and the G-allele of rs2428151 (EDA) were both associated with a higher risk of hypodontia (odds ratio, OR = 2.79, 95% CI = 1.11-7.01; and OR = 2.87, 95% CI = 1.04-7.94, respectively). The G-allele of rs2520378 (EDA) showed a protective effect with an OR of 0.61 (95% CI = 0.38-0.99). The EDA SNP findings were consistent with previous reports included in a meta-analysis. No associations were found with the PAX9, AXIN2 and MSX1 genes, after adjusting for sex and ethnicity. CONCLUSIONS: Common variants of the EDA genes are associated with specific phenotypes of non-syndromic hypodontia, thus confirming their role in the regulatory pathways of normal tooth development. However, larger samples are needed to investigate the association further.