The gout epidemic in French Polynesia: a modelling study of data from the Ma'i u'u epidemiological survey

BACKGROUND

Gout is the most common cause of inflammatory arthritis worldwide, particularly in Pacific regions. We aimed to establish the prevalence of gout and hyperuricaemia in French Polynesia, their associations with dietary habits, their comorbidities, the prevalence of the HLA-B*58:01 allele, and current management of the disease.

METHODS

The Ma'i u'u survey was epidemiological, prospective, cross-sectional, and gout-focused and included a random sample of adults from the general adult population of French Polynesia. It was conducted and data were collected between April 13 and Aug 16, 2021. Participants were randomly selected to represent the general adult population of French Polynesia on the basis of housing data collected during the 2017 territorial census. Each selected household was visited by a research nurse from the Ma'i u'u survey who collected data via guided, 1-h interviews with participants. In each household, the participant was the individual older than 18 years with the closest upcoming birthday. To estimate the frequency of HLA-B*58:01, we estimated HLA-B haplotypes on individuals who had whole-genome sequencing to approximately 5× average coverage (mid-pass sequencing). A subset of individuals who self-reported Polynesian ancestry and not European, Chinese, or other ancestry were used to estimate Polynesian-ancestry specific allele frequencies. Bivariate associations were reported for weighted participants; effect sizes were estimated through the odds ratio (OR) of the association calculated on the basis of a logistic model fitted with weighted observations.

FINDINGS

Among the random sample of 2000 households, 896 participants were included, 140 individuals declined, and 964 households could not be contacted. 22 participants could not be weighted due to missing data, so the final weighted analysis included 874 participants (449 [51·4%] were female and 425 [48·6%] were male) representing the 196 630 adults living in French Polynesia. The estimated prevalence of gout was 14·5% (95% CI 9·9-19·2), representing 28 561 French Polynesian adults, that is 25·5% (18·2-32·8) of male individuals and 3·5% (1·0-6·0) of female individuals. The prevalence of hyperuricaemia was estimated at 71·6% (66·7-76·6), representing 128 687 French Polynesian adults. In multivariable analysis, age (OR 1·5, 95% CI 1·2-1·8 per year), male sex (10·3, 1·8-60·7), serum urate (1·6, 1·3-2·0 per 1 mg/dL), uraturia (0·8, 0·8-0·8 per 100 mg/L), type 2 diabetes (2·1, 1·4-3·1), BMI more than 30 kg/m2 (1·1, 1·0-1·2 per unit), and percentage of visceral fat (1·7, 1·1-2·7 per 1% increase) were associated with gout. There were seven heterozygous HLA-B*58:01 carriers in the full cohort of 833 individuals (seven [0·4%] of 1666 total alleles) and two heterozygous carriers in a subset of 696 individuals of Polynesian ancestry (two [0·1%]).

INTERPRETATION

French Polynesia has an estimated high prevalence of gout and hyperuricaemia, with gout affecting almost 15% of adults. Territorial measures that focus on increasing access to effective urate-lowering therapies are warranted to control this major public health problem.

FUNDING

Variant Bio, the French Polynesian Health Administration, Lille Catholic University Hospitals, French Society of Rheumatology, and Novartis.

A Polynesian-specific missense CETP variant alters the lipid profile

Identifying population-specific genetic variants associated with disease and disease-predisposing traits is important to provide insights into the genetic determinants of health and disease between populations, as well as furthering genomic justice. Various common pan-population polymorphisms at CETP associate with serum lipid profiles and cardiovascular disease. Here, sequencing of CETP identified a missense variant rs1597000001 (p.Pro177Leu) specific to Māori and Pacific people that associates with higher HDL-C and lower LDL-C levels. Each copy of the minor allele associated with higher HDL-C by 0.236 mmol/L and lower LDL-C by 0.133 mmol/L. The rs1597000001 effect on HDL-C is comparable with CETP Mendelian loss-of-function mutations that result in CETP deficiency, consistent with our data, which shows that rs1597000001 lowers CETP activity by 27.9%. This study highlights the potential of population-specific genetic analyses for improving equity in genomics and health outcomes for population groups underrepresented in genomic studies.

Genetic testing for misclassified monogenic diabetes in Māori and Pacific peoples in Aōtearoa New Zealand with early-onset type 2 diabetes

Aims

Monogenic diabetes accounts for 1-2% of diabetes cases yet is often misdiagnosed as type 2 diabetes. The aim of this study was to examine in Māori and Pacific adults clinically diagnosed with type 2 diabetes within 40 years of age, (a) the prevalence of monogenic diabetes in this population (b) the prevalence of beta-cell autoantibodies and (c) the pre-test probability of monogenic diabetes.

Methods

Targeted sequencing data of 38 known monogenic diabetes genes was analyzed in 199 Māori and Pacific peoples with BMI of 37.9 ± 8.6 kg/m² who had been diagnosed with type 2 diabetes between 3 and 40 years of age. A triple-screen combined autoantibody assay was used to test for GAD, IA-2, and ZnT8. MODY probability calculator score was generated in those with sufficient clinical information (55/199).

Results

No genetic variants curated as likely pathogenic or pathogenic were found. One individual (1/199) tested positive for GAD/IA-2/ZnT8 antibodies. The pre-test probability of monogenic diabetes was calculated in 55 individuals with 17/55 (31%) scoring above the 20% threshold considered for diagnostic testing referral.

Discussion

Our findings suggest that monogenic diabetes is rare in Māori and Pacific people with clinical age, and the MODY probability calculator likely overestimates the likelihood of a monogenic cause for diabetes in this population.

Riboceine Rescues Auranofin-Induced Craniofacial Defects in Zebrafish

Craniofacial abnormalities are a common group of congenital developmental disorders that can require intensive oral surgery as part of their treatment. Neural crest cells (NCCs) contribute to the facial structures; however, they are extremely sensitive to high levels of oxidative stress, which result in craniofacial abnormalities under perturbed developmental environments. The oxidative stress-inducing compound auranofin (AFN) disrupts craniofacial development in wildtype zebrafish embryos. Here, we tested whether the antioxidant Riboceine (RBC) rescues craniofacial defects arising from exposure to AFN. RBC rescued AFN-induced cellular apoptosis and distinct defects of the cranial cartilage in zebrafish larvae. Zebrafish embryos exposed to AFN have higher expression of antioxidant genes gstp1 and prxd1, with RBC treatment partially rescuing these gene expression profiles. Our data suggest that antioxidants may have utility in preventing defects in the craniofacial cartilage owing to environmental or genetic risk, perhaps by enhancing cell survival.

Evolution and genomic organization of the insect sHSP gene cluster and coordinate regulation in phenotypic plasticity

Background

Conserved syntenic gene complexes are rare in Arthropods and likely only retained due to functional constraint. Numerous sHSPs have been identified in the genomes of insects, some of which are located clustered in close proximity. Previous phylogenetic analyses of these clustered sHSP have been limited to a small number of holometabolous insect species and have not determined the pattern of evolution of the clustered sHSP genes (sHSP-C) in insect or Arthropod lineages.

Results

Using eight genomes from representative insect orders and three non-insect arthropod genomes we have identified that a syntenic cluster of sHSPs (sHSP-C) is a hallmark of most Arthropod genomes. Using 11 genomes from Hymenopteran species our phylogenetic analyses have refined the evolution of the sHSP-C in Hymenoptera and found that the sHSP-C is order-specific with evidence of birth-and-death evolution in the hymenopteran lineage. Finally we have shown that the honeybee sHSP-C is co-ordinately expressed and is marked by genomic features, including H3K27me3 histone marks consistent with coordinate regulation, during honeybee ovary activation.

Conclusions

The syntenic sHSP-C is present in most insect genomes, and its conserved coordinate expression and regulation implies that it is an integral genomic component of environmental response in arthropods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01885-8.

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.

Gout, Rheumatoid Arthritis, and the Risk of Death Related to Coronavirus Disease 2019: An Analysis of the UK Biobank

Objectives

The objectives for this study were to assess whether gout and/or rheumatoid arthritis (RA) are risk factors for coronavirus disease 2019 (COVID‐19) diagnosis and to assess whether gout and/or RA are risk factors for death from COVID‐19.

Methods

We used data from the UK Biobank. Multivariable‐adjusted logistic regression was employed in the following analyses: analysis A, to test for association between gout and/or RA and COVID‐19 diagnosis (n = 473,139); analysis B, to test for association between gout and/or RA and death from COVID‐19 in a case‐control cohort of people who died of or survived COVID‐19 (n = 2059); analysis C, to test for association between gout and/or RA and death from COVID‐19 in the entire UK Biobank cohort (n = 473,139).

Results

RA, but not gout, was associated with COVID‐19 diagnosis in analysis A. Neither RA nor gout was associated with risk of death in the group diagnosed with COVID‐19 in analysis B. However, RA was associated with risk of death related to COVID‐19 by using the UK Biobank cohort in analysis C, independent of comorbidities and other measured risk factors (odds ratio [OR] 1.9; 95% confidence interval CI 1.2–3.0). Gout was not associated with death related to COVID‐19 in the same UK Biobank analysis (OR 1.2; 95% CI 0.8–1.7).

Conclusion

RA is a risk factor for death from COVID‐19 by using the UK Biobank cohort. These findings require replication in larger data sets that also allow for inclusion of a wider range of factors.

Corrigendum: Functional Urate-Associated Genetic Variants Influence Expression of lincRNAs LINC01229 and MAFTRR

[This corrects the article DOI: 10.3389/fgene.2018.00733.].

A non-coding genetic variant maximally associated with serum urate levels is functionally linked to HNF4A-dependent PDZK1 expression

The precise molecular mechanisms by which urate-associated genetic variants affect urate levels are unknown. Here, we tested for functional linkage of the maximally associated genetic variant rs1967017 at the PDZK1 locus to elevated PDZK1 expression. We performed expression quantitative trait loci (eQTL) and likelihood analyses and gene expression assays. Zebrafish were used to evaluate tissue-specific gene expression. Luciferase assays in HEK293 and HepG2 cells measured the effect of rs1967017 on transcription amplitude. Probabilistic Annotation Integrator analysis revealed rs1967017 as most likely to be causal and rs1967017 was an eQTL for PDZK1 in the intestine. The region harboring rs1967017 was capable of directly driving green fluorescent protein expression in the kidney, liver and intestine of zebrafish embryos, consistent with a conserved ability to confer tissue-specific expression. Small interfering RNA depletion of HNF4A reduced endogenous PDZK1 expression in HepG2 cells. Luciferase assays showed that the T allele of rs1967017 gains enhancer activity relative to the urate-decreasing C allele, with T allele enhancer activity abrogated by HNF4A depletion. HNF4A physically binds the rs1967017 region, suggesting direct transcriptional regulation of PDZK1 by HNF4A. Computational prediction of increased motif strength, together with our functional assays, suggests that the urate-increasing T allele of rs1967017 strengthens a binding site for the transcription factor HNF4A. Our and other data predict that the urate-raising T allele of rs1967017 enhances HNF4A binding to the PDZK1 promoter, thereby increasing PDZK1 expression. As PDZK1 is a scaffold protein for many ion channel transporters, increased expression can be predicted to increase activity of urate transporters and alter excretion of urate.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Antioxidant treatment ameliorates phenotypic features of SMC1A-mutated Cornelia de Lange syndrome in vitro and in vivo

Cornelia de Lange syndrome (CdLS) is a rare disease characterized by cognitive impairment, multisystemic alterations and premature aging. Furthermore, CdLS cells display gene expression dysregulation and genomic instability. Here, we demonstrated that treatment with antioxidant drugs, such as ascorbic acid and riboceine, reduced the level of genomic instability and extended the in vitro lifespan of CdLS cell lines. We also found that antioxidant treatment partially rescued the phenotype of a zebrafish model of CdLS. Gene expression profiling showed that antioxidant drugs caused dysregulation of gene transcription; notably, a number of genes coding for the zinc finger (ZNF)-containing Krueppel-associated box (KRAB) protein domain (KRAB-ZNF) were found to be downregulated. Taken together, these data suggest that antioxidant drugs have the potential to ameliorate the developmental phenotype of CdLS.

Functional Urate-Associated Genetic Variants Influence Expression of lincRNAs LINC01229 and MAFTRR

Genetic variation in the genomic regulatory landscape likely plays a crucial role in the pathology of disease. Non-coding variants associated with disease can influence the expression of long intergenic non-coding RNAs (lincRNAs), which in turn function in the control of protein-coding gene expression. Here, we investigate the function of two independent serum urate-associated signals (SUA1 and SUA2) in close proximity to lincRNAs and an enhancer that reside ∼60 kb and ∼300 kb upstream of MAF, respectively. Variants within SUA1 are expression quantitative trait loci (eQTL) for LINC01229 and MAFTRR, both co-expressed with MAF. We have also identified that variants within SUA1 are trans-eQTL for genes that are active in kidney- and serum urate-relevant pathways. Serum urate-associated variants rs4077450 and rs4077451 within SUA2 lie within an enhancer that recruits the transcription factor HNF4α and forms long range interactions with LINC01229 and MAFTRR. The urate-raising alleles of rs4077450 and rs4077451 increase enhancer activity and associate with increased expression of LINC01229. We show that the SUA2 enhancer region drives expression in the zebrafish pronephros, recapitulating endogenous MAF expression. Depletion of MAFTRR and LINC01229 in HEK293 cells in turn lead to increased MAF expression. Collectively, our results are consistent with serum urate variants mediating long-range transcriptional regulation of the lincRNAs LINC01229 and MAFTRR and urate relevant genes (e.g., SLC5A8 and EHHADH) in trans.

An ovine hepatorenal fibrocystic model of a Meckel-like syndrome associated with dysmorphic primary cilia and TMEM67 mutations

Meckel syndrome (MKS) is an inherited autosomal recessive hepatorenal fibrocystic syndrome, caused by mutations in TMEM67, characterized by occipital encephalocoele, renal cysts, hepatic fibrosis, and polydactyly. Here we describe an ovine model of MKS, with kidney and liver abnormalities, without polydactyly or occipital encephalocoele. Homozygous missense p.(Ile681Asn; Ile687Ser) mutations identified in ovine TMEM67 were pathogenic in zebrafish phenotype rescue assays. Meckelin protein was expressed in affected and unaffected kidney epithelial cells by immunoblotting, and in primary cilia of lamb kidney cyst epithelial cells by immunofluorescence. In contrast to primary cilia of relatively consistent length and morphology in unaffected kidney cells, those of affected cyst-lining cells displayed a range of short and extremely long cilia, as well as abnormal morphologies, such as bulbous regions along the axoneme. Putative cilia fragments were also consistently located within the cyst luminal contents. The abnormal ciliary phenotype was further confirmed in cultured interstitial fibroblasts from affected kidneys. These primary cilia dysmorphologies and length control defects were significantly greater in affected cells compared to unaffected controls. In conclusion, we describe abnormalities involving primary cilia length and morphology in the first reported example of a large animal model of MKS, in which we have identified TMEM67 mutations.

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0623-3) contains supplementary material, which is available to authorized users.